3

オブジェクトの正確な地図の縮尺を m または km 単位で調べ、SpatialPolygonsそれをプロットしてspplot、読者が距離を把握できるようにしたいと考えています。???つまり、プロットの に入力する数値を知りたいということです。

library("sp")
library("CARBayes")

data("spatialhousedata")
proj4string(spatialhousedata) <- CRS("+proj=cass +lat_0=52.41864827777778 +lon_0=13.62720366666667 
                                      +x_0=40000 +y_0=10000 +datum=potsdam +units=m +no_defs 
                                      +ellps=bessel +towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7 ")

scalebar=list("SpatialPolygonsRescale", layout.scale.bar(),
              offset = c(225000,647000), scale = 10000, fill=c("transparent","black"))
text1=list("sp.text",  c(225000,649000), "0")
text2=list("sp.text",  c(230000,649000), "??? m")

spplot(spatialhousedata, "price",
       sp.layout=list(scalebar, text1, text2))

多分何かsp::coordinatesまたはImap::gdist？SpatialPolygonsこれを見つけるために必要な事前の再投影はありますか?

4

1 に答える 1

3

データが投影されます (試してくださいis.projected(spatialhousedata))。私はのために行きます

text2=list("sp.text",  c(235000,649000), "10 km")

スケール バーは、軸スケールを描画することで確認できます。

spplot(spatialhousedata, "price", scales = list(draw=TRUE), 
    sp.layout=list(scalebar, text1, text2))
于 2015-03-27T14:54:18.353 に答える