geom_sf（）でsfオブジェクトをプロットし、latlong以外の投影で描画する

Question

geom_sf（）でポリゴンをlatlong以外のプロジェクションにプロットしようとしています。

私は、データセットのインポート manual pages for geom_sf()

で見つかった例使用しています

：ようやくのCRSを定義ggplotでプロット3857

nc_3857 <- sf::st_transform(nc, "+init=epsg:3857") 

：

nc <- sf::st_read(system.file("shape/nc.shp", package = "sf"), quiet = TRUE) 

EPSGにLatLongブログからの転換をプロット：

ggplot() + geom_sf(data = nc_3857, colour = "red", fill = NA) + coord_sf(crs=st_crs(3857)) 

私はwags84（すなわちepsg：4326）の地図を緯度経度軸で取得しています。私はメートル単位で軸を持っているので、投影されたポリゴンをプロットするにはggplotが必要です。私は間違って何をしていますか？事前 ゼ

Answer 1

で

おかげでそれはちょうどLAT-長い目盛りを重ね、要求された投影でそれをプロットしています。

ノルウェーと同様のことを試してみると、たとえば、北極に近いと、ディスプレイのXY座標は変形の座標ですが、オーバーレイする目盛りは緯度経度です。これはEPSG 3035（円錐形）であるノルウェーのマップでの位置座標：

だからある

が投影ポリゴンをプロットします。ここでの緯線の線がグリッドの場合、緯度経度投影で座標をプロットしていたはずです。

ドキュメント内目盛の唯一の言及はcoord_sfの引数である：

本当に多くを言うことはありません

datum: CRS that provides datum to use when generating graticules 

。

デカルト座標系が必要ですか？ああ試すことができます：

> ggplot() + geom_sf(data = rp, colour = "red", fill = NA) + coord_cartesian() 
Error: geom_sf() must be used with coord_sf() 

はgeom_sfと交互目盛用ggplot2の問題を確認し、そこに何もない場合は、問題を追加します。

Answer 2

もhttps://github.com/tidyverse/ggplot2/issues/2200を参照してみてください

を与える

ggplot() + geom_sf(data = nc_3857, colour = "red", fill = NA) + 
    coord_sf(datum=st_crs(3857))