亚洲精品久久久中文字幕-亚洲精品久久片久久-亚洲精品久久青草-亚洲精品久久婷婷爱久久婷婷-亚洲精品久久午夜香蕉

您的位置:首頁(yè)技術(shù)文章
文章詳情頁(yè)

GDAL 矢量屬性數(shù)據(jù)修改方式(python)

瀏覽:23日期:2022-08-02 17:48:46

Case:需要給一個(gè)現(xiàn)有的shp數(shù)據(jù)創(chuàng)建一個(gè)字段,并將屬性表中原有的一個(gè)文本類(lèi)型的屬性轉(zhuǎn)換為整型后填入新創(chuàng)建的字段。

Problem:新字段創(chuàng)建成功,但是賦值操作無(wú)效,即無(wú)法成功給字段寫(xiě)入值。

solution:對(duì)字段進(jìn)行賦值后需要,重新寫(xiě)入Feature,否則賦值無(wú)效,即layer0.SetFeature(feature)。

特別注意:在對(duì)數(shù)據(jù)進(jìn)行讀寫(xiě)操作,一定要以讀寫(xiě)的方式打開(kāi),即Open(filePath,1),該方法的原型為Open(pszName,int bUpdate = false),并且具有返回值,參數(shù)說(shuō)明如下:

名稱(chēng) 說(shuō)明 pszName 需要打開(kāi)文件或數(shù)據(jù)源的路徑 bUpdate 是否需要更新數(shù)據(jù)集,默認(rèn)為只讀,如果需要對(duì)數(shù)據(jù)進(jìn)行讀寫(xiě)操作,需要給此參數(shù)賦值1 返回值 返回?cái)?shù)據(jù)集的指針,如果為NULL,則表明打開(kāi)數(shù)據(jù)集失敗

以下將給出創(chuàng)建字段和獲取字段值、子段賦值的實(shí)例。

1.為矢量數(shù)據(jù)創(chuàng)建字段

# 添加字段 defn = layer.GetLayerDefn() fieldIndex=defn.GetFieldIndex(’SSSS’) if fieldIndex<0: # 添加字段 fieldDefn = ogr.FieldDefn(’SSSS’, ogr.OFTInteger) fieldDefn.SetPrecision(9) layer0.CreateField(fieldDefn,1); fieldIndex2 = defn.GetFieldIndex(’SSSS’) if fieldIndex2>0: print('字段創(chuàng)建成功:',fieldIndex)

2.獲取字段值及字段賦值

feature = layer.GetNextFeature() indexA = defn.GetFieldIndex(’code’) indexB = defn.GetFieldIndex(’SSSS’) oField = defn.GetFieldDefn(indexB) fieldName = oField.GetNameRef() while feature is not None: valueA= feature.GetFieldAsInteger(indexA) if valueA is None: feature.SetFieldNull(indexB) continue feature.SetField2(fieldName, valueA) layer0.SetFeature(feature) feature = layer0.GetNextFeature() #feature.Destroy() ds.Destroy()

補(bǔ)充知識(shí):Python批量修改shapefile屬性表字段名(arcpy增刪字段)

嘗試了3種方法,時(shí)間緊迫屢敗屢戰(zhàn),最后終于成功。

方法1和2是失敗記錄,希望有類(lèi)似經(jīng)歷成功的同學(xué)分享下經(jīng)驗(yàn)。

方法3是成功記錄,修改 = 新增 + 計(jì)算 +刪除相關(guān)字段,因?yàn)闆](méi)有刪除改名前字段需求,故刪除部分沒(méi)有另做。

方法1 —— 【將shapefile的dbf文件按csv文件讀寫(xiě)】

——失敗,dbf中有空值及編碼問(wèn)題

不同文件間通過(guò)改后綴簡(jiǎn)單粗暴改寫(xiě),操作方法存在風(fēng)險(xiǎn)

過(guò)于依賴熟悉的領(lǐng)域,由于時(shí)間緊迫沒(méi)有試用dbf第三方庫(kù) GDAL 矢量屬性數(shù)據(jù)修改方式(python)

#-*- coding: utf-8 -*- import osimport shutilimport csv #批量修改shp中dbf文件中的字段名#LANE_WIDTH改為L(zhǎng)ANEWIDTHL#LANE_WID_1改為L(zhǎng)ANEWIDTHR#ORIGIN_LIN改為ORG_LINKID roadDir = r'D:20Q100DATAModifyTitleROAD_TEST'modifyRoadDir = r'D:20Q100DATAModifyTitleModifiedLink'csvDir = r'D:20Q100DATAModifyTitleROAD_CSV'fileList = os.listdir(roadDir) num = 0for file in fileList: #新建修改后的城市文件夾 city = file.split('_')[0] os.mkdir(csvDir + city) os.mkdir(modifyRoadDir + city) shutil.copy(roadDir + file + '/RD_LINK.dbf', csvDir + city + '/RD_LINK.csv') print('正在生成{}新dbf文件......'.format(city)) dbfData = [] #讀取dbf數(shù)據(jù)為csv文件,讀存內(nèi)容部分 csvFile = open(csvDir + city + '/RD_LINK.csv',encoding=’gbk’,errors=’ignore’) csvReader = csv.reader(csvFile) for row in csvReader: titleLine = [] if csvReader.line_num == 1: titleLine.append(row) continue #跳過(guò)第1行——列名 dbfData.append(row) #生成正確的字段名行 newTitleLine = [] modifyDic = {'LANE_WIDTH':'LANEWIDTHL','LANE_WID_1':'LANEWIDTHR','ORIGIN_LIN':'ORG_LINKID'} for fieldName in titleLine: if fieldName in modifyDic: newTitleLine.append(modifyDic[fieldName]) else: newTitleLine.append(fieldName) #寫(xiě)入新dbf文件 newDBF = open(modifyRoadDir + city + '/RD_LINK.dbf',’w’) csvWriter = csv.writer(newDBF) csvWriter.writerow(newTitleLine) for row in dbfData: csvWriter.writerow(row) newDBF.close() print('{}新dbf文件已生成!'.format(city)) print('{}城市全部完成'.format(num))

方法2——直接使用修改字段名函數(shù)——失敗,arcpy模塊沒(méi)有AlterField_management方法

——失敗,但發(fā)現(xiàn)直接探尋官方方法還是比網(wǎng)搜野路子要節(jié)約時(shí)間

GDAL 矢量屬性數(shù)據(jù)修改方式(python)

使用Arcgis10.2 - Advanced浮動(dòng)版,符合許可信息但Arcpy調(diào)用函數(shù)失敗,存疑。附官網(wǎng)AlterField函數(shù)用法介紹:

https://pro.arcgis.com/zh-cn/pro-app/tool-reference/data-management/alter-field-properties.htm

#-*- coding: utf-8 -*- import osimport arcpy #批量修改shp中dbf文件中的字段名#LANE_WIDTH改為L(zhǎng)ANEWIDTHL#LANE_WID_1改為L(zhǎng)ANEWIDTHR#ORIGIN_LIN改為ORG_LINKID roadDir = r'D:20Q100DATAModifyTitleROAD_TEST'fileList = os.listdir(roadDir) num = 0for file in fileList: #截取城市名 city = file.split('_')[0] print('正在修改{}的shp字段......'.format(city)) #修改前后字段名對(duì)照字典 modifyDic = {'LANE_WIDTH':'LANEWIDTHL','LANE_WID_1':'LANEWIDTHR','ORIGIN_LIN':'ORG_LINKID'} #讀取shp數(shù)據(jù) dbfFile = roadDir + file + '/RD_LINK.shp' fieldList = arcpy.ListFields(dbfFile) for field in fieldList: #遍歷字段名 if field.name.upper() in modifyDic: #找到待修改字段名 arcpy.AlterField_management(dbfFile, field = field.name, new_field_name = modifyDic[field.name]) print('{}字段修改成功!'.format(city)) print('{}城市全部完成'.format(num))

方法3——添加改名后新字段(字段屬性與前保持一致)、計(jì)算改名后字段值 = 改名前字段值,成功!

#-*- coding: utf-8 -*- import osimport arcpy #批量修改shp屬性表中的字段名#LANE_WIDTH改為L(zhǎng)ANEWIDTHL#LANE_WID_1改為L(zhǎng)ANEWIDTHR#ORIGIN_LIN改為ORG_LINKID roadDir = r'D:20Q100DATAModifyTitleROAD_TEST'fileList = os.listdir(roadDir) num = 0for file in fileList: city = file.split('_')[0] #截取城市名 print('正在添加和計(jì)算{}的shp正確字段......'.format(city)) #修改前后字段名對(duì)照屬性字典 modifyDic = {'LANE_WIDTH':['LANEWIDTHL','1'],#錯(cuò)誤字段名:[正確字段名,長(zhǎng)度] 'LANE_WID_1':['LANEWIDTHR','1'], 'ORIGIN_LIN':['ORG_LINKID','10']} #讀取shp文件 shpFile = roadDir + file + '/RD_LINK.shp' for wrongfieldName in modifyDic: correctfieldName = modifyDic[wrongfieldName][0] field_length = modifyDic[wrongfieldName][1] # Process: 添加字段 arcpy.AddField_management(shpFile, correctfieldName, 'TEXT', '', '', field_length) # Process: 計(jì)算字段字段 arcpy.CalculateField_management(shpFile, correctfieldName, '['+wrongfieldName+']', 'VB', '') ## # Process: 添加字段——LANEWIDTHL## arcpy.AddField_management(shpFile, 'LANEWIDTHL', 'TEXT', '', '', '1')## # Process: 計(jì)算字段字段——LANEWIDTHL## arcpy.CalculateField_management(shpFile, 'LANEWIDTHL', '['++]', 'VB', '')#### # Process: 添加字段——LANEWIDTHR## arcpy.AddField_management(shpFile, 'LANEWIDTHR', 'TEXT', '', '', '1')## # Process: 計(jì)算字段字段——LANEWIDTHR## arcpy.CalculateField_management(shpFile, 'LANEWIDTHR', '[LANE_WID_1]', 'VB', '')#### # Process: 添加字段——ORG_LINKID## arcpy.AddField_management(shpFile, 'ORG_LINKID', 'TEXT', '', '', '10')## # Process: 計(jì)算字段字段——ORG_LINKID## arcpy.CalculateField_management(shpFile, 'ORG_LINKID', '[ORIGIN_LIN]', 'VB', '')#### # Process: 刪除字段## #arcpy.DeleteField_management(in_table, 'LANE_WIDTH') print('{}修改完成!請(qǐng)確認(rèn)!'.format(city)) num += 1 print('{}城市全部完成!'.format(num))

以上這篇GDAL 矢量屬性數(shù)據(jù)修改方式(python)就是小編分享給大家的全部?jī)?nèi)容了,希望能給大家一個(gè)參考,也希望大家多多支持好吧啦網(wǎng)。

標(biāo)簽: Python 編程
相關(guān)文章:
主站蜘蛛池模板: 日韩一级免费视频 | 国产日韩一区二区三区在线观看 | 青娱乐成人 | 久久久999国产精品 久久久99视频 | 精品日产一区二区三区手机 | 一级欧美一级日韩 | 在线观看亚洲精品国产 | 中文字幕日韩哦哦哦 | 久青草免费视频手机在线观看 | 在线播放免费人成毛片乱码 | 国产乱码精品一区二区三 | 999热这里只有精品 999热精品这里在线观看 | 欧美人成一本免费观看视频 | 国产伦精品一区二区 | 黄色大片黄色大片 | 国产日产精品_国产精品毛片 | 一级国产20岁美女毛片 | 97视频久久| 国产精品a区 | 欧美精品在线视频 | 日韩精品一区二区三区免费视频 | 中文xxx视频| 国产成人精品午夜二三区 | 久久国产热视频 | 黄色视屏在线 | 国产福利精品在线 | 欧美亚洲另类视频 | 免费在线a| 91久久色| 97玖玖| 看全大色黄大色黄大片一级爽 | 国产人妖自拍 | a级高清免费 | 欧美在线观看免费一区视频 | 久久网国产 | 99久久免费精品视频 | 国产疯狂露脸对白 | 国产乱理论片在线观看理论 | 丁香婷婷综合五月六月 | 国外成人免费高清激情视频 | chinese麻豆自制国产 |