Estou serializando minhas geometrias do arco-íris como geojson, para poder 'hidratá-las novamente como geometrias mais tarde e estou tendo 2 problemas no ciclo .:
PROBLEMA 1: Precisão
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
Se eu verificar a representação da string, as coordenadas mudaram ligeiramente:
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)]]}
PROBLEMA 2: Furos Se o polígono possui furos, geo_interface gera um erro:
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'
Alguma idéia de como resolver esses problemas?
arcpy
geojson
attributeerror
Víctor Velarde
fonte
fonte
Respostas:
OK - bem, eu pensei que tinha resolvido.
substitua a linha ~ 80 deste arquivo C: \ Python26 \ ArcGIS10.0 \ Lib \ arcpy \ arcobjects \ geometries.py a partir disso:
para isso (ou algo que seja mais conciso e elegante e faça a mesma coisa):
Basicamente, eles esqueceram de considerar rosquinhas na forma que são marcadas por valores de pontos nulos. Isso gera um bom geoJson (partes separadas), mas o método arcpy.AsShape destrói o GeoJSON.
este código:
gera isso:
Desisto. ;)
Atualização O problema dos buracos foi resolvido na versão 10.1 com este pedaço de python:
fonte
C:\Program Files\ArcGIS\Server\arcpy\arcpy\arcobjects\geometries.py
em 10.1, mas se você estiver em 10.0, poderá corrigi-lo..py
arquivo. Eu pensei que ele fez um service pack para 10, mas acho que não.