使用ArcPy如何计算面图层的椭球面积？

Q: 请问谁有计算面图层的椭球面积 的代码吗？

A: 具体步骤如下

假设 GeoJSON 数据如下：

{"type": "Polygon", 
"coordinates": [[
  [-102.05, 41.0], 
  [-102.05, 37.0], 
  [-109.05, 37.0], 
  [-109.05, 41.0]
]]}

这里的坐标为经纬度，可以使用 pyproj 和 Shapely 来定位这个区域，代码如下：

co = {"type": "Polygon", "coordinates": [
    [(-102.05, 41.0),
    (-102.05, 37.0),
    (-109.05, 37.0),
    (-109.05, 41.0)]]}
lon, lat = zip(*co['coordinates'][0])
from pyproj import Proj
pa = Proj("+proj=aea +lat_1=37.0 +lat_2=41.0 +lat_0=39.0 +lon_0=-106.55")

这是一个相等的区域投影，以感兴趣的区域为中心并将其包围。 现在制作新的投影GeoJSON表示，转换为Shapely几何对象，并取以下区域：

x, y = pa(lon, lat)
cop = {"type": "Polygon", "coordinates": [zip(x, y)]}
from shapely.geometry import shape
shape(cop).area  # 268952044107.43506

怎样提问？

目前主要的方式是通过QQ群来提问，大家也可以互相交流互相解决问题。已经开通的QQ群请见文末最下方的声明。感兴趣的可以加入，有的如果已经满员，请换一个再加。你所提问的问题，也可能会被整理到麻辣GIS中，也期待您的参与。

所有问答

本站整理的所有问题请参考：所有GIS问答