1

このウェブサイトのおかげで、気象観測所の降雨データを 1 km のグリッドに補間し、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))

評判ポイントが足りないので載せられない画像になってしまいました。エイミーウェイ、このプロットは、さまざまな色の降雨量を示す英国の地図で、最大降雨量が赤で最小降水量が青で示されています。

これでもかまいませんが、次の 2 つのことを行いたいと考えています。

  1. プロットするすべての降雨データが 0 以上であることを確認してください。一部の降雨データは負の値 (例: -1) ですが、プロットする前にすべての負の値を 0 に変更するにはどうすればよいですか?

  2. 最小値と最大値 (青から赤) の間の色範囲を使用するのではなく、値をビニングしたいので、たとえば、青は 0mm から 5mm、水色は 5mm から 10mm、黄色は 10mm から 20mm、赤はそれより大きい30mmなど。これを達成する方法はありますか?

4

1 に答える 1

0

ローマンのアイデアから盗むが、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})
于 2013-12-11T13:05:28.293 に答える