日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當(dāng)前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

Arcpy计算字段

發(fā)布時(shí)間:2023/12/10 编程问答 35 豆豆
生活随笔 收集整理的這篇文章主要介紹了 Arcpy计算字段 小編覺得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.
# coding=utf-8 # 字段分類 import arcpy import os# 獲取當(dāng)前shp的所有字段 def getFieldNames(fc):fieldList = []for f in arcpy.ListFields(fc):fieldList.append(str(f.name))return fieldListRecalss_area = """ def Reclass(area1):if (area1 <= 250):return 0elif (area1 > 250 and area1 <= 500):return 1elif (area1 > 500 and area1 <= 1000):return 2elif (area1 > 1000 and area1 <= 2000):return 3elif (area1 > 2000 and area1 <= 4000):return 4elif (area1 > 4000 and area1 <= 7000):return 5elif (area1 > 7000 and area1 <= 10000):return 6elif (area1 > 10000):return 7 """Reclass_shape = """ def Reclass(shp1):if (shp1 <= 18):return 0elif (shp1 <= 25):return 1elif (shp1 > 25 and shp1 <= 30):return 2elif (shp1 > 30 and shp1 <= 40):return 3elif (shp1 > 40 and shp1 <= 55):return 4elif (shp1 > 55 and shp1 <= 70):return 5elif (shp1 > 70 and shp1 <=90):return 6elif (shp1 > 90):return 7 """Reclass_angle = """ def Reclass(angle):if (angle <= -77.5):return 0elif (angle > -77.5 and angle <= -45):return 1elif (angle > -45 and angle <= -22.5):return 2elif (angle > -22.5 and angle <=0):return 3elif (angle > 0 and angle <= 22.5):return 4elif (angle > 22.5 and angle <= 45):return 5elif (angle > 45 and angle <=77.5):return 6elif (angle > 77.5):return 7 """Reclass_extn = """ def Reclass(extn):if (extn <= 1.5):return 0elif (extn > 1.5 and extn <= 2):return 1elif (extn > 2 and extn <= 3):return 2elif (extn > 3 and extn <=5):return 3elif (extn > 5 and extn <= 7):return 4elif (extn > 7 and extn <= 9):return 5elif (extn > 9 and extn <=10):return 6elif (extn > 10):return 7 """list = ['Area_8', 'Shape_8', 'Angle_8', 'extn_8'] # list = ["nghr", "nghr_1"] arcpy.env.workspace = r"C:\Users\user\Downloads\data2\b1" fcl = arcpy.ListFeatureClasses() for shp in fcl:print shpfieldList = getFieldNames(shp)print fieldList# 如果shp中已有分類字段,則刪除for p in list:if p in fieldList:arcpy.DeleteField_management(shp, p)# 添加用于分類的字段arcpy.AddField_management(shp, "Area_8", "SHORT", "", "", "", "", "NULLABLE", "NON_REQUIRED", "")arcpy.AddField_management(shp, "Shape_8", "SHORT", "", "", "", "", "NULLABLE", "NON_REQUIRED", "")arcpy.AddField_management(shp, "Angle_8", "SHORT", "", "", "", "", "NULLABLE", "NON_REQUIRED", "")arcpy.AddField_management(shp, "Extn_8", "SHORT", "", "", "", "", "NULLABLE", "NON_REQUIRED", "")# 開始分類try:arcpy.CalculateField_management(shp, "Area_8", "Reclass(!area1!)", "PYTHON_9.3", Recalss_area)except arcpy.ExecuteError:print arcpy.GetMessages()try:arcpy.CalculateField_management(shp, "Shape_8", "Reclass(!shp1!)", "PYTHON_9.3", Reclass_shape)except arcpy.ExecuteError:print arcpy.GetMessages()try:arcpy.CalculateField_management(shp, "Angle_8", "Reclass(!angle!)", "PYTHON_9.3", Reclass_angle)except arcpy.ExecuteError:print arcpy.GetMessages()try:arcpy.CalculateField_management(shp, "Extn_8", "Reclass(!extn!)", "PYTHON_9.3", Reclass_extn)except arcpy.ExecuteError:print arcpy.GetMessages()fieldList2 = getFieldNames(shp)print fieldList2

總結(jié)

以上是生活随笔為你收集整理的Arcpy计算字段的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò),歡迎將生活随笔推薦給好友。