arcpy.geometry __geo_interface__ 및 AsShape () 함수 : 정밀도 및 구멍 손실

나는 arcpy 기하학을 geojson으로 직렬화하여 나중에 기하학으로 다시 수화 할 수 있고주기에 2 가지 문제가 있습니다. :

문제 1 : 정밀도

    R0 = arcpy.SearchCursor(self.shpTest, "FID=0").next().getValue("Shape")          
    geojson = R0.__geo_interface__                        
    R1 = arcpy.AsShape(geojson)
    self.assertTrue(R0.equals(R1)) <<< THIS FAILS

문자열 표현을 확인하면 좌표가 약간 변경되었습니다.

    geojson2 = R1.__geo_interface__
    print geojson
    print geojson2  

    {'type': 'Polygon', 'coordinates': [[(442343.5516410945, 4814166.6184399202), (442772.17749834526, 4811610.7383281607), (441565.67508534156, 4811499.6131059099), (440772.50052100699, 4814184.7808806188), (442343.5516410945, 4814166.6184399202)]]}
    {'type': 'Polygon', 'coordinates': [[(442343.55169677734, 4814166.6185302734), (442772.17749023438, 4811610.73828125), (441565.67510986328, 4811499.6130981445), (440772.50048828125, 4814184.7808837891), (442343.55169677734, 4814166.6185302734)]]}

문제 2 : 구멍 다각형에 구멍이 있으면 geo_interface 에서 오류가 발생합니다.

    R0_WithHoles = arcpy.SearchCursor(self.shpTest, "FID=0").next().getValue("Shape")          
    geojson = R0.__geo_interface__  <<< generates this ERROR:

    File "C:\Program Files\ArcGIS\Desktop10.0\arcpy\arcpy\arcobjects\geometries.py", line 68, in __geo_interface__
        return {'type': 'Polygon', 'coordinates': [[(pt.X, pt.Y) for pt in part] for part in self]}
    AttributeError: 'NoneType' object has no attribute 'X'

이러한 문제를 해결하는 방법에 대한 아이디어가 있습니까?

알았어. 잘 풀 었다고 생각 했어.

이 파일 C : \ Python26 \ ArcGIS10.0 \ Lib \ arcpy \ arcobjects \ geometries.py의 ~ 80 행을 다음과 같이 바꾸십시오.

return {'type': 'Polygon', 'coordinates': [[(pt.X, pt.Y) for pt in part] for part in self]}

이것 (또는 더 간결하고 우아하고 같은 일을하는 것)에 :

  obj = {"type": "Polygon"}
    coordinates = []
    for part in self:
        _part = []
        for pt in part:
            if pt is not None:
                print pt
                _part.append([pt.X,pt.Y])
            else:
                print "none"
                coordinates.append(_part)
                _part=[]
        coordinates.append(_part)
    obj["coordinates"]=coordinates
    return obj

기본적으로 그들은 널 포인트 값으로 표시된 모양의 도넛을 고려하는 것을 잊었습니다. 이것은 좋은 geoJson (별도의 부품)을 뱉어 내지 만 arcpy.AsShape 메소드는 GeoJSON을 폐기합니다.

이 코드 :

import arcpy
gj = {
  'type': 'Polygon', 'coordinates': [
   [[-122.803764, 45.509158], [-122.796246, 45.500050], [-122.808193, 45.500109],
      [-122.803764, 45.509158]],
   [[-122.804206, 45.504509], [-122.802882, 45.502522], [-122.801866, 45.504479], 
      [-122.804206, 45.504509]]
   ]
 }

 p = arcpy.AsShape(gj)
 print p.__geo_interface__

이것을 출력합니다 :

    {'type': 'Polygon', 'coordinates': [[[-122.8037109375, 45.50927734375],  
    [-122.79620361328125, 45.5001220703125], [-122.80810546875, 45.5001220703125],
    [-122.8037109375, 45.50927734375]]]}

나는 포기한다. ;)

업데이트 구멍 문제는 10.1 에서이 파이썬 덩어리로 해결되었습니다.

return {'type': 'Polygon', 'coordinates': [[((pt.X, pt.Y) if pt else None)
                                                    for pt in part]
                                                        for part in self]}