Question

感谢这个网站的帮助，我现在设法编写了我的代码，将气象站降雨数据插入1公里的网格，然后使用 R 进行绘图。最后一段代码消除了不需要的数据，然后使用quilt.plot绘制。要绘制的数据（mydf.final）是具有以下标题的值列表：索引，东向，北向，降雨和平均降雨量。如果我执行dput：

，则要绘制的数据样本如下所示

dput(mydf.final(head, 10))

structure(list(easting = c(101000, 101000, 101000, 101000, 101000, 
101000, 101000, 101000, 101000, 101000), northing = c(740000, 
741000, 742000, 743000, 744000, 745000, 746000, 747000, 748000, 
749000), rainfall = c(40.0997151135538, 40.3344163486536, 40.5681562532368, 
40.8009001543464, 41.0326137107226, 41.2632629265481, 41.4928141651644, 
41.721234162219, 41.9484900390941, 42.1745493154662), saar = c(1081.62, 
1081.62, 1081.62, 1081.62, 1081.62, 1078.68, 1078.68, 1078.68, 
1078.68, 1078.68)), .Names = c("easting", "northing", "rainfall", 
"saar"), row.names = 4717:4726, class = "data.frame")

我可以使用以下方式绘制此数据：

quilt.plot(cbind(mydf.final$easting,mydf.final$northing),mydf.final$rainfall, add.legend=TRUE,nx=599,ny=1209,xlim=c(0,700000),ylim=c(0,1250000))

这给了我一张我无法发布的图片，因为我没有足够的声望点。 Amyway，该地块是英国的地图，显示不同颜色的降雨量，最大降雨量显示红色和最小蓝色。

虽然这没关系，但我想做两件事：

确保要绘制的所有降雨数据均为0或更大。一些降雨量数据是负值（例如-1），那么在绘图之前如何将所有负值更改为0？
我不想使用最小值和最大值（蓝色到红色）之间的颜色范围，而是想要对这些值进行分级，例如，蓝色为0mm至5mm，浅蓝色为5mm至10mm，黄色为10mm至20mm，红色大于30mm等等。我是如何实现这一目标的？

Answer 1

从罗马的想法中窃取，但改为使用findInterval：

y<- mydf.final$rainfall  # simplify typing
ycuts <- findInterval(y,c(0,5,10,20,30,max(y)))
quilt.plot(yourxdata, y , col=c('blue','lightblue','yellow','purple','red')[ycuts], {the rest of your arguments})

如何使用分箱数据生成R中的被子图？

1 个答案:

如何使用分箱数据生成R中​​的被子图？

1 个答案:

如何使用分箱数据生成R中的被子图？