10

所以我在 R 中有一个 spatialpolygons 对象；但我不确定为什么我无法从中检索“区域”插槽。

这是我的 R 会话：

> spatialpolygons
An object of class "SpatialPolygons"
Slot "polygons":
[[1]]
An object of class "Polygons"
Slot "Polygons":
[[1]]
An object of class "Polygon"
Slot "labpt":
[1] 20.50516 57.72918

Slot "area":
[1] 36.85484

Slot "hole":
[1] FALSE

Slot "ringDir":
[1] 1

Slot "coords":
         [,1]     [,2]
[1,] 16.48438 59.73633
[2,] 22.59277 61.14258
[3,] 24.74609 55.03418
[4,] 17.49512 55.12207
[5,] 16.48438 59.73633



Slot "plotOrder":
[1] 1

Slot "labpt":
[1] 20.50516 57.72918

Slot "ID":
[1] "myMultiPolygons"

Slot "area":
[1] 36.85484



Slot "plotOrder":
[1] 1

Slot "bbox":
       min      max
x 16.48438 24.74609
y 55.03418 61.14258

Slot "proj4string":
CRS arguments:
 +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 

> spatialpolygons@bbox
       min      max
x 16.48438 24.74609
y 55.03418 61.14258
> spatialpolygons@area
Error: no slot of name "area" for this object of class "SpatialPolygons"
> slotNames(spatialpolygons)
[1] "polygons"    "plotOrder"   "bbox"        "proj4string"
> names(spatialpolygons)
[1] "myMultiPolygons"
4

2 回答 2

19

首先，您应该知道@area插槽不是关于SpatialPolygons*对象实际区域的可靠信息来源。如 中所述?"Polygons-class"，@area插槽仅用作绘图的辅助工具（防止较小的多边形被较大的多边形覆盖）并且不尊重投影或正确解释多边形中的孔。

要获得准确的区域，您应该改为使用rgeos::gArea()具有投影坐标参考系的图层或geosphere::areaPolygon()用于经纬度坐标参考系（即CRS(+proj=longlat)）的图层。

排除这些警告后，下面将展示@area如果您确实需要插槽的内容，如何获取它们。

主要的复杂性是区域槽属于Polygon对象，而不是SpatialPolygons对象（其中Polygon对象是一个元素）。因此，您需要首先深入研究SpatialPolygons对象以提取到单个Polygon对象。

您已经这样做了，您可以只使用@运算符来提取区域插槽的内容。

以下示例使用在包 vignette (warning, pdf)的第 7 节中创建的SpatialPolygons对象：sp

require(sp)
# Example pasted in from Section 7 of the sp vignette
Sr1 = Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2)))
Sr2 = Polygon(cbind(c(5,4,2,5),c(2,3,2,2)))
Sr3 = Polygon(cbind(c(4,4,5,10,4),c(5,3,2,5,5)))
Sr4 = Polygon(cbind(c(5,6,6,5,5),c(4,4,3,3,4)), hole = TRUE)
Srs1 = Polygons(list(Sr1), "s1")
Srs2 = Polygons(list(Sr2), "s2")
Srs3 = Polygons(list(Sr3, Sr4), "s3/4")
SpP = SpatialPolygons(list(Srs1,Srs2,Srs3), 1:3)

# To extract the area of the first (or in your case only) Polygon
SpP@polygons[[1]]@area
# [1] 5.5

# Extract the areas of all three component Polygons
sapply(SpP@polygons, function(x) x@area)
# [1]  5.5  1.5 10.0

## For areas, rgeos::gArea() or geosphere::areaPolygons() are generally more appropriate
## (Note, for instance, that it properly accounts for the hole in the 3rd polygon.)
rgeos::gArea(SpP, byid=TRUE)
#  s1   s2 s3/4 
# 5.5  1.5  9.0
于 2012-01-03T06:45:06.040 回答
8

您可以使用包中的函数计算面积rgeos，下面的示例使用 Josh 的示例数据。这可能更合适，因为area插槽仅用于绘图。

library(rgeos)
gArea(SpP[1,])
## [1] 5.5
gArea(SpP[2,])
##[1] 1.5
gArea(SpP[3,])
## [1] 10

一次全部：

gArea(SpP)
[1] 17

应考虑使用的坐标系，这只是原始几何区域。

帮助页面讨论了该area插槽。

?gArea
....

请注意，该值可能与“多边形”类的“面积”槽不同，因为该值不会减去几何体中任何孔的面积。

?"Polygons-class"
...

'area'：类'"numeric"'的对象；多边形列表的总平面面积，但不是重复计算孔（从 0.9-58 更改 - 对岛求和，孔被忽略而不是减去）；这些值用于确保在较大区域的多边形之后绘制较小区域的多边形，不尊重投影，因为此类的对象没有定义投影

于 2012-01-03T08:09:31.273 回答