0

Rを使用してGoogleマップにポイント（緯度と経度を指定）をプロットする簡単な方法を見つけようとしています。 -maps-in-r ... より具体的には、「meuse」データセットを使用して、そこで説明されている 2 番目の方法です。

ただし、コードを理解できません。正確には何ですか -

coordinates(meuse)<-~x+y # convert to SPDF
proj4string(meuse) <- CRS('+init=epsg:28992')

行う？

http://rwiki.sciviews.org/doku.php?id=tips:spatial-data:change_crsをチェックしましたが、あまり光を当てませんでした。

4

1 に答える 1

2

これはそれがどのようになるかです：

library(sp)
data(meuse)
str(meuse)

> str(meuse)
'data.frame':   155 obs. of  14 variables:
 $ x      : num  181072 181025 181165 181298 181307 ...
 $ y      : num  333611 333558 333537 333484 333330 ...
 $ cadmium: num  11.7 8.6 6.5 2.6 2.8 3 3.2 2.8 2.4 1.6 ...
 $ copper : num  85 81 68 81 48 61 31 29 37 24 ...

meusedata.frame です。

この行を実行すると

coordinates(meuse)<-~x+y

x 変数と y 変数 (列) を使用して SpatialPolygonsDataFrame オブジェクトを作成します。このオブジェクトには、(x と y からの) 座標と、さまざまな金属や土地の特徴に関する属性があります。

曲面 (私たちは大きな球体に住んでいます) を平面に投影するのは簡単ではないため、これを実現するにはいくつかの計算が必要です。これは、データがどの投影であるかをRに伝える場所です。したがって、行

proj4string(meuse) <- CRS('+init=epsg:28992')
于 2013-11-10T08:29:27.953 に答える