本文主要是介绍arcpy 解析土地整治报备坐标点为是矢量面，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！
数据格式：
代码片段：
# coding: utf-8import os
import arcpy
import json
import sysdef set_encoding():a, b, c = sys.stdin, sys.stdout, sys.stderrreload(sys)sys.setdefaultencoding('utf-8')sys.stdin, sys.stdout, sys.stderr = a, b, cdef get_files(root):files = os.listdir(root)return filesdef create_feature(root, workspace, sp_ref):temp_dic_root, database = get_path(workspace)temp_dic = os.path.join(temp_dic_root, "temp_data")if not os.path.exists(temp_dic):os.makedirs(temp_dic)for sourcefile in get_files(root):if not sourcefile.endswith(".txt"):continuefilename = os.path.splitext(sourcefile)[0]out_file = workspace + "\\v_" + filenamejsonfile = os.path.join(temp_dic, filename + ".json")sourcefile = root + "\\" + sourcefilecoordinates = analysis(sourcefile)if create_geojson(jsonfile, coordinates, sp_ref):arcpy.JSONToFeatures_conversion(jsonfile, out_file)arcpy.AddMessage("info: " + filename)else:arcpy.AddWarning("waring: " + filename)def create_geojson(targetfile, coordinates, sp_ref):features_template = {"displayFieldName": "","fieldAliases": {"FID": "FID"},"geometryType": "esriGeometryPolygon","spatialReference": {"wkt": ""},"fields": [{"name": "FID","type": "esriFieldTypeOID","alias": "FID"}],"features": []}features_template["spatialReference"]["wkt"] = sp_reffor fid, item in enumerate(coordinates):features_template["features"].append(get_geometry(fid, item))with open(targetfile, "w") as out_file:out_file.write(json.dumps(features_template))return Truedef get_geometry(fid, rings):feature_template = {"attributes": {"FID": 0}, "geometry": {"rings": []}}feature_template["attributes"]["FID"] = fidtemp = [ring for index, ring in rings.items()]feature_template["geometry"]["rings"].extend(temp)return feature_templatedef analysis(sourcefile):geometry = []rings = {}with open(sourcefile, "r") as in_file:line = in_file.readline()while line:flag1 = "@" in line.replace("\n", "").split(",")flag2 = line.startswith("J")if not (flag1 or flag2):line = in_file.readline()continueif flag1:geometry.append(rings)rings = {}elif flag2:temp = line.split(",")ring_index = temp[1]ring_part = list(map(float, reversed(temp[2:])))if not rings.has_key(ring_index):rings[ring_index] = []rings[ring_index].append(ring_part)line = in_file.readline()geometry.append(rings)geometry.pop(0)return geometrydef get_path(path):temps = path.split('\\')result = []if '.' in temps[-1]:root = "\\".join(temps[0:-1])result.append(root)result.append(temps[-1])return resultif __name__ == "__main__":set_encoding()root = arcpy.GetParameterAsText(0)workspace = arcpy.GetParameterAsText(1)sp_ref = arcpy.GetParameterAsText(2)sp_ref = sp_ref.replace("\'", "\"")create_feature(root, workspace, sp_ref)
工具参数：
工具界面：
这篇关于arcpy 解析土地整治报备坐标点为是矢量面的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！
