麻辣GIS微信平台

更多 GIS 干货

微信关注不错过

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

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

A: 具体步骤如下

假设 GeoJSON 数据如下:

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

这里的坐标为经纬度,可以使用 pyprojShapely 来定位这个区域,代码如下:

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问答

没有下文

麻辣GIS-Sailor
作者:
GIS爱好者,学GIS,更爱玩GIS。

声明

1.本文所分享的所有需要用户下载使用的内容(包括但不限于软件、数据、图片)来自于网络或者麻辣GIS粉丝自行分享,版权归该下载资源的合法拥有者所有,如有侵权请第一时间联系本站删除。

2.下载内容仅限个人学习使用,请切勿用作商用等其他用途,否则后果自负。

3.如果本文对您有所帮助,请您对我们团队进行 打赏捐助让我们在传播3S的路上可以走得更远,不胜感激。PS:如果你是学生党,请优先把经费用于购买学习资料 以及 与小哥哥/小姐姐约会上:-)

如您有疑问,可在文末留言,或到QQ群中提问。

本站QQ群(一):291616564 加入QQ群

本站QQ群(二):166408035 加入QQ群

本站QQ群(三):627853279 加入QQ群

本站QQ群(四):436386604 加入QQ群

本站QQ群(五):606176554 加入QQ群

微信公众号:malagis,扫描右边二维码直接关注。

发表评论