3

R / ggplot を使用して、世界の地震帯を含む世界地図を作成しています。ここから変更されたシェープファイルを使用しています: http://gmo.gfz-potsdam.de/。

プロットには 9 つのレイヤーと 4400 を超えるグループがあります。ただし、一部のレイヤーには穴があるはずですが、穴がないように見えます。したがって、これを回避するには、最初にこれらのポリゴンをレンダリングして、背面に表示されるようにする必要があります。geom_polygon を使用する場合、ポリゴンのレンダリング順序を変更するにはどうすればよいですか?

require("rgdal") # requires sp, will use proj.4 if installed
require("maptools")
require("ggplot2")
gpclibPermit() # required for fortify method
require("plyr")

fill_colours <- c("#dddddd", "#8FC463FF", "#46A742FF", "#FFFC00FF", "#F2B212FF", "#E9726AFF", "#D92C50FF", "#FF1B23FF", "#B1681FFF")

eq.shp = readOGR(dsn=".", "PGA_PolygonsPY v2_region")

# Create a discrete scale, there are probably more efficient ways of doing this
eq.shp$MAXVALUE[eq.shp$MAXVALUE==10.0] = 9
eq.shp$MAXVALUE[eq.shp$MAXVALUE== 4.8] = 8
eq.shp$MAXVALUE[eq.shp$MAXVALUE== 4.0] = 7
eq.shp$MAXVALUE[eq.shp$MAXVALUE== 3.2] = 6
eq.shp$MAXVALUE[eq.shp$MAXVALUE== 2.4] = 5
eq.shp$MAXVALUE[eq.shp$MAXVALUE== 1.6] = 4
eq.shp$MAXVALUE[eq.shp$MAXVALUE== 0.8] = 3
eq.shp$MAXVALUE[eq.shp$MAXVALUE== 0.4] = 2
eq.shp$MAXVALUE[eq.shp$MAXVALUE== 0.2] = 1
eq.shp$MAXVALUE  <- sprintf("%d", eq.shp$MAXVALUE)

eq.shp@data$id = rownames(eq.shp@data)
eq.points = fortify(eq.shp, region='id')

# Join the data frames back together
eq.df = join(eq.points, eq.shp@data, by='id')

p <- ggplot()
p <- p + geom_polygon(data = eq.df, aes(long, lat, group=group, fill=MAXVALUE))
p <- p + scale_fill_manual(values=fill_colours) 
p <- p + coord_equal(ratio = 1, xlim=c(65, 85), ylim=c(35, 45))
p

このコードは、次のプロットを生成します。
ggplot出力

上記のプロットを元のデータ データ (GIS ソフトウェアを使用してプロット) と比較すると、次のアーティファクトが見られます。これは、茶色のポリゴンが最後にレンダリングされるためです。- 同様に、周囲の大きな緑色の領域 (経度 = 80、緯度 = 38) は、すべて黄色です。これは、黄色のポリゴンが緑色のポリゴンの後にレンダリングされるためです。- 欠けている穴、例えば (経度 = 74、緯度 = 39) の近くの領域についてはあまり心配していません。

GIS ソフトウェアからの元のデータ

4

1 に答える 1

3

さて、私はこれを解決することができました。これはちょっとした回避策であり、レンダリング時間が長くなりますが、私の要件は満たしています。

ggplot を使用してポリゴンのレンダリングをソート/並べ替えようとする代わりに、for ループを使用してこれを強制し、必要に応じてレイヤーをレンダリングしました。

ソリューション：

次のコードを置き換えます。

p <- p + geom_polygon(data = eq.df, aes(long, lat, group=group, fill=MAXVALUE))

次のように:

for (i in seq(9, 1, -1)) {
    p <- p + geom_polygon(data = eq.df[eq.df$MAXVALUE==sprintf("%d", i),], aes(long, lat, group=group, fill=MAXVALUE))
}

これにより、希望する順序でポリゴンがレンダリングされます。次に、元のものと非常によく似たものを取得します。

ここに画像の説明を入力

ただし、新しい問題が発生したことがわかります。穴のないポリゴンがいくつかあるため、それらをレンダリングすると、他のレイヤーの一部が覆われてしまいます。したがって、ポリゴンを正しくレンダリングするには、ポリゴンに穴を作成する必要があります。

コメントを残してくれたすべての人に感謝します。

更新: 上記の解決策は、ポリゴンのレンダリングを並べ替えるために機能しますが、根本的な問題を完全に解決するわけではありません。そのためには、より洗練されたアルゴリズムが必要になるようです。また、ポリゴンの面積とポリゴンの加重面積で並べ替えてみましたが、どちらも理想的な解決策を達成することはできませんでした。最終的な解決策には、手動での介入が必要です。

于 2012-08-22T06:40:03.110 に答える