Rでロッジポールパインシードプランニングユニットのシェープファイルを読み込みましたreadShapePoly。
spu <- readShapePoly("spus.shp")
棘だけをマッピングすると、それらが表示されます。
plot(spu, border=TRUE)
シェイプファイルをカナダのマップに追加しようとすると、spus が表示されるはずですが、表示されません。
map("worldHires","Canada", xlim=c(-141,-110), ylim=c(45,65), col="gray90", fill=TRUE)
plot(spu, add=TRUE, border=TRUE)
これは、シェープファイルに .prj ファイルに含まれる投影情報が欠落しているためだと思います。
PROJCS["NAD_1983_Albers",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],
PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Albers"],PARAMETER["False_Easting",1000000.0],
PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-126.0],PARAMETER["Standard_Parallel_1",50.0],
PARAMETER["Standard_Parallel_2",58.5],PARAMETER["Latitude_Of_Origin",45.0],UNIT["Meter",1.0]]
の部分的な出力
summary(spu)
Object of class SpatialPolygonsDataFrame
Coordinates:
min max
x 642162.8 1870556
y 457057.2 1421478
Is projected: NA
proj4string : [NA]
Data attributes:
proj4stringプロジェクション情報が R の spu オブジェクトに関連付けられていないことを示しています。正しいwithを使用してこれを修正できると思いますreadShapePoly,が、.prj ファイルの情報を proj4 文字列に変換する方法が見つかりませんでした。どうすればこれを行うことができますか/これは正しいアプローチですか?