计算地理质心的经度/纬度 [英] Calculate longitude/latitude for geographic centroid

问题描述

我想对内布拉斯加州的县级作物产量数据进行一些空间统计分析，用于STAT类。为此，我需要每个县的地理质心的经度和纬度。有人知道在R里怎么做吗？我知道它可以在ArcGIS中完成，但我现在无法访问它。

推荐答案

您没有提供任何详细信息，您的shapefile是从here获得的，但我从here获得了一个，您可以使用gCentroidrgeos：

library(rgdal)
library(sp)
library(rgeos)

nebraska <- readOGR("CountyBoundsUTM/", "CountyUTM")

gCentroid(nebraska, byid=TRUE)

## SpatialPoints:
##           x       y
## 0  721768.5 4636738
## 1  430938.8 4524651
## 2  698036.4 4566570
## 3  370970.6 4641340
## ...
## 89 623301.6 4603228
## 90 618883.0 4486931
## 91 439295.3 4582756
## 92 493680.8 4522680
## Coordinate Reference System (CRS) arguments: +proj=utm +zone=14 +datum=NAD83
## +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0 

这篇关于计算地理质心的经度/纬度的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持IT屋！