2017-10-10 20 views
0

私はbroom :: tidy関数に問題があります。これは、シェイプファイル(アラスカのすべての川の地図です)に​​固有のものです。fortify、broom :: tidyを使用して間違った緯度/経度

私は

readOGR機能を使用してファイルをロードし、その後、ほうき::整頓機能を使用して、データフレームに空間データを強化。

shp_points <- tidy(shp, region="group") 

いずれもうまくいくようです。 ggplotを使って点をプロットすると、空間的に意味のある地図が得られます。

ggplot(data=shape_points, aes(x=long, y=lat, group=group)) + 
    geom_path() + 
    coord_fixed() 

map of alaska rivers

あなたは緯度と長軸を見ればしかし、私は認識していない形式です。私はそれらのための明白な変換要因も見ません。誰もがこれを経験する?

セッション情報:
Rバージョン3.4.2(2017年9月28日)
プラットフォーム:x86_64版-W64-MINGW32/x64(64ビット)
の下で実行します。Windows> = 8のx64(9200を構築します)
他の取り付けパッケージ:
1 broom_0.4.2 dplyr_0.7.4 rgeos_0.3-25 ggplot2_2.2.1 rgdal_1.2-13 sp_1.2-5

+1

データは投影されている可能性があります。メタデータは何を言いましたか? –

+0

ソリューションは機能しましたか? –

答えて

0

第一歩

後4326

library(sp) 

shp <- spTransform(shp,"+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0") 

そして、あなたの残りのステップに従ってください:

はEPSGにあなたのCRSを変換してみます。あなたの軸は十進数の緯度と経度でなければなりません

+0

あなたの助けてくれてありがとう! spTransform(xSP、CRSobj、...)のエラー: NA参照システムからの変換ができません –