[英]How to find out the map scale of a SpatialPolygons object in R?
我想找出我的SpatialPolygons
对象的m或km的精确地图比例,并将其绘制在我的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
需要找到它吗?
您的数据是预测的(尝试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))
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.