我已经将一个 python 脚本与 ArcGIS arcpy 放在一起，用于创建 polygonzo json 多边形 ( http://code.google.com/p/polygonzo/ )。这是我的python脚本...
import os, string, arcpy
arcpy.env.overwriteOutput = True
layer = "C:\\Other\\Shapefiles\\Geo500K_JSON\\GEOLOGY_500K_Project.shp"
output = "C:\\Other\\Shapefiles\\Geo500K_JSON\\"
outfile = output + "Geo500K.json"
jsonFile = open(outfile,'w')
jsonFile.write('var geo = {\n')
jsonFile.write('\t"type": "FeatureCollection",\n')
jsonFile.write('\t"features": [\n')
idfield = "ORIG_LABEL"
shape_field = arcpy.Describe(layer).shapeFieldName
rows = arcpy.SearchCursor(layer,"","","",idfield + " A")
row = rows.next()
while row:
geostring = '' #for each lat/lng pt
geolist = [] # array for storing individual geostrings
ringList = [] #array for storing geolist array with geostrings separated by commas
partList = [] #array for storing partlist, final array used
shapeString = ''
jsonFile.write('\t\t{"type": "Feature", ')
extent = row.Shape.extent
ne = str(extent.XMax) + ',' + str(extent.YMax)
sw = str(extent.XMin) + ',' + str(extent.YMin)
jsonFile.write('"bbox": [' + sw + ', ' + ne + '],')
jsonFile.write('"properties":{')
geoLabel = str(row.getValue(idfield))
jsonFile.write('"label": "' + geoLabel + '", ')
geoName = str(row.getValue("FM_NAME"))
jsonFile.write('"name": "' + geoName + '", ')
lithType = str(row.getValue("LithType"))
jsonFile.write('"lithType": "' + lithType + '", ')
rank = str(row.getValue("Rank"))
jsonFile.write('"rank": "' + rank + '", ')
lithName = str(row.getValue("LithName"))
jsonFile.write('"lithName": "' + lithName + '", ')
ageType = str(row.getValue("AgeType"))
jsonFile.write('"ageType": "' + ageType + '", ')
minAge = str(row.getValue("MinAge"))
jsonFile.write('"minAge": "' + minAge + '", ')
maxAge = str(row.getValue("MaxAge"))
jsonFile.write('"maxAge": "' + maxAge + '", ')
part = row.getValue(shape_field).centroid
jsonFile.write('"center":[' + str(part.X) + ',' + str(part.Y) + '],')
jsonFile.write('"centroid":[' + str(part.X) + ',' + str(part.Y) + ']},')
jsonFile.write('"geometry":{"type":"MultiPolygon","coordinates":[[[')
feat = row.shape
for p in range(feat.partCount):
pInt = p
part = feat.getPart(p)
pt = part.next()
while pt:
lat = str(round(pt.Y,6))
lon = str(round(pt.X,6))
geostring = '[' + lon + ',' + lat + ']'
geolist.append(geostring)
pt = part.next()
#if now following point go to the next part which should be an interior ring.
if not pt:
ringList.append(',' .join(geolist))
geostring = ''
geolist = []
pt = part.next()
if pt:
print 'Interior Ring: ' + geoLabel
partList.append(',' .join(ringList))
ringList = []
shapeString = ']], [[' .join(partList)
jsonFile.write(shapeString)
jsonFile.write(']]]}},\n')
row = rows.next()
#jsonFile.seek(-1, os.SEEK_END)
#jsonFile.truncate()
jsonFile.write('\t]\n')
jsonFile.write('}')
jsonFile.close()
del row, rows
当脚本遇到内部环时，它只打印一个警告。我不知道如何处理它们。不幸的是，我使用的许多多边形都有内环。我使用一个具有内环的多边形组合了一张测试地图。这是它的样子...... http://www.geology.ar.gov/test/test-polygonzo.html
polyzo 可以处理内环吗？
更新：我非常感谢您的回复 Michael Geary 先生！但是，我无法使用 json 模块让您的 python 脚本工作。其中有一些错误，我在上面对其进行了编辑，但它吐出了一个空白文档。可能是我不够努力。在查看了带有内环的多面体在 json 格式中应该是什么样子的示例之后，我回到了我的 python 脚本的工作上（是的，如果不使用 json 模块，让 json 有效有点困难）。我添加了更多评论，所以如果你有时间，你可以使用 json 模块让你的脚本工作——我想看看一个工作示例。这是我最终的python脚本....
import os, string, arcpy
arcpy.env.overwriteOutput = True
layer = "C:\\Other\\Shapefiles\\Geo500K_JSON\\GEOLOGY_500K_kn.shp"
output = "C:\\Other\\Shapefiles\\Geo500K_JSON\\"
outfile = output + "Geo500K_knTest.json"
jsonFile = open(outfile,'w')
jsonFile.write('var geo = {\n')
jsonFile.write('\t"type": "FeatureCollection",\n')
jsonFile.write('\t"features": [\n')
idfield = "ORIG_LABEL"
shape_field = arcpy.Describe(layer).shapeFieldName
rows = arcpy.SearchCursor(layer,"","","",idfield + " A")
row = rows.next()
#loop through the attribute table
while row:
jsonFile.write('\t\t{"type": "Feature", \n')
extent = row.Shape.extent
ne = str(extent.XMax) + ',' + str(extent.YMax)
sw = str(extent.XMin) + ',' + str(extent.YMin)
jsonFile.write('\t\t"bbox": [' + sw + ', ' + ne + '],\n')
jsonFile.write('\t\t"properties":{\n')
geoLabel = str(row.getValue(idfield))
jsonFile.write('\t\t\t"label": "' + geoLabel + '", \n')
geoName = str(row.getValue("FM_NAME"))
jsonFile.write('\t\t\t"name": "' + geoName + '", \n')
lithType = str(row.getValue("LithType"))
jsonFile.write('\t\t\t"lithType": "' + lithType + '", \n')
rank = str(row.getValue("Rank"))
jsonFile.write('\t\t\t"rank": "' + rank + '", \n')
lithName = str(row.getValue("LithName"))
jsonFile.write('\t\t\t"lithName": "' + lithName + '", \n')
ageType = str(row.getValue("AgeType"))
jsonFile.write('\t\t\t"ageType": "' + ageType + '", \n')
minAge = str(row.getValue("MinAge"))
jsonFile.write('\t\t\t"minAge": "' + minAge + '", \n')
maxAge = str(row.getValue("MaxAge"))
jsonFile.write('\t\t\t"maxAge": "' + maxAge + '", \n')
centroid = row.getValue(shape_field).centroid
jsonFile.write('\t\t\t"center":[' + str(centroid.X) + ',' + str(centroid.Y) + '], \n')
jsonFile.write('\t\t\t"centroid":[' + str(centroid.X) + ',' + str(centroid.Y) + '] \n')
jsonFile.write('\t\t\t}, \n') #end of properties
jsonFile.write('\t\t"geometry":{\n\t\t\t"type":"MultiPolygon",\n\t\t\t"coordinates":[\n')
feat = row.shape #get the shape/geography of the row in the attribute table
partnum = 1
#loop through the parts of the polygon (some may have more that one part)
for p in range(feat.partCount):
jsonFile.write('\t\t\t\t[\n\t\t\t\t\t[\n')
jsonFile.write('\t\t\t\t\t\t//Part ' + str(partnum) + '\n')
jsonFile.write('\t\t\t\t\t\t//Outer ring of Part ' + str(partnum) + '\n')
part = feat.getPart(p) #return an array of point objects for particular part
pt = part.next() #return specific pt object of array
innerRingNum = 1
#loop through each pt object/vertex of part
while pt:
lat = round(pt.Y,7) #get latitude of pt object and round to 7 decimal places
lon = round(pt.X,7) #get longitude of pt object and round to 7 decimal places
jsonFile.write('\t\t\t\t\t\t[' + str(lon) + ',' + str(lat) + '],\n') #assemble [lon,lat]
pt = part.next() #go to next pt object to continue loop
#if no following point go to the next part which should be an interior ring.
if not pt:
#we've got an interior ring so let's loop through the vertices of the ring
pt = part.next()
if pt:
jsonFile.seek(-3, os.SEEK_END)
jsonFile.truncate() #remove trailing comma
jsonFile.write('\n\t\t\t\t\t],\n')
jsonFile.write('\t\t\t\t\t[\n')
jsonFile.write('\t\t\t\t\t\t//Inner ring ' + str(innerRingNum) + ' of Part ' + str(partnum) + '\n')
print 'Interior Ring: ' + geoLabel
innerRingNum += 1
partnum += 1
jsonFile.seek(-3, os.SEEK_END)
jsonFile.truncate() #remove trailing comma
jsonFile.write('\n\t\t\t\t\t]\n\t\t\t\t],\n')
jsonFile.seek(-3, os.SEEK_END)
jsonFile.truncate() #remove trailing comma
jsonFile.write('\n\t\t\t]\n\t\t\t}\n\t\t},\n')
row = rows.next()
jsonFile.seek(-3, os.SEEK_END)
jsonFile.truncate() #remove trailing comma
jsonFile.write('\n\t]\n')
jsonFile.write('}')
jsonFile.close()
del row, rows
我还要补充一点，polygonzo 以及您与他人分享的意愿给我留下了深刻的印象。但是，您提供的 javascript 和 python 确实可以使用更多注释来更快地理解这一切。