如何使用sf :: st_centroid来计算多边形的质心？

我想使用新的“sf”包在R中操作一些巴西人口普查数据。我能够导入数据，但我得到一个错误，当我尝试创建原始多边形如何使用sf :: st_centroid来计算多边形的质心？

library(sf) 

#Donwload data 
filepath <- 'ftp://geoftp.ibge.gov.br/organizacao_do_territorio/malhas_territoriais/malhas_de_setores_censitarios__divisoes_intramunicipais/censo_2010/setores_censitarios_shp/ac/ac_setores_censitarios.zip' 
download.file(filepath,'ac_setores_censitarios.zip') 
unzip('ac_setores_censitarios.zip') 
d <- st_read('12SEE250GC_SIR.shp',stringsAsFactors = F)

现在我试图创建一个包含列“几何”的心一个新的几何列的重心，但出现错误：

d$centroid <- st_centroid(d$geometry) 
Warning message: 
In st_centroid.sfc(d$geometry) : 
    st_centroid does not give correct centroids for longitude/latitude data

我该如何解决这个问题？

来源

2017-09-12 LucasMation

这不是一个错误，它是一个警告。值被创建。 –

底层为sf的所有GEOS功能需要投影坐标才能正常工作，因此您应该在适当投影的数据上运行st_centroid。我对巴西可用的CRS不太了解，但EPSG：29101似乎工作正常：

library(tidyverse) 

d$centroids <- st_transform(d, 29101) %>% 
    st_centroid() %>% 
    # this is the crs from d, which has no EPSG code: 
    st_transform(., '+proj=longlat +ellps=GRS80 +no_defs') %>% 
    # since you want the centroids in a second geometry col: 
    st_geometry() 

# check with 
plot(st_geometry(d)) 
plot(d[, 'centroids'], add = T, col = 'red', pch = 19)

来源

2017-10-04 07:15:35

如何使用sf :: st_centroid来计算多边形的质心？

回答

相关问题