python计算椭球面积代码_使用ArcPy如何计算面图层的椭球面积？

59 阅读 0 评论 39 点赞

我是靠谱客的博主活泼鱼，最近开发中收集的这篇文章主要介绍python计算椭球面积代码_使用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", &#

如果觉得靠谱客网站的内容还不错，欢迎将靠谱客网站推荐给程序员好友。

本图文内容来源于网友提供，作为学习参考使用，或来自网络收集整理，版权属于原作者所有。