问题描述
我正在尝试查找人口普查区和县的区域,其中有一些不连续。我正在使用 python area
库。这是我适用于所有连续形状的代码。
obj = {'type':'polygon','coordinates': feature['geometry']['coordinates']}
current_area = area(obj) / square_meter_to_square_mile
这是一个坐标的例子:
{'type': 'polygon','coordinates': [[[-79.891153,40.497034,0],[-79.890443,40.494777,[-79.889832,40.491527,[-79.884604,40.490922,[-79.884484,40.488949,[-79.884592,40.487329,[-79.888175,40.48644,[-79.8882,40.483948,[-79.883772,40.483406,[-79.879559,40.483247,[-79.876632,40.483316,[-79.874238,40.48355,[-79.865498,40.485069,[-79.865786,40.491525,[-79.865968,40.495515,[-79.862811,40.49765,[-79.859216,40.497693,[-79.859175,40.496435,[-79.856864,40.498169,[-79.851113,40.493601,[-79.850339,40.49594,[-79.84968,40.497589,[-79.848853,40.503658,[-79.849326,40.506821,[-79.849451,40.507425,[-79.850546,40.514188,[-79.850665,40.51721,[-79.850494,40.520508,[-79.851755,40.520607,[-79.85488,40.520812,[-79.85527,40.528789,[-79.867579,40.528512,[-79.867259,40.522247,[-79.871734,40.526388,[-79.872237,40.527164,[-79.871867,40.517058,[-79.871676,40.509746,[-79.87278,40.50971,[-79.872672,40.503888,[-79.878623,40.503694,[-79.879079,40.499454,[-79.883633,40.499245,[-79.887353,40.498669,[-79.891153,0]]]}
这是它的样子:
解决方法
错误不是我认为的地方。那组坐标都是一个形状,没有问题。对于有两个形状的那些,我必须遍历子形状列表并分别处理每个区域。