原則: 優(yōu)先級(jí)從上到下:(以地理國(guó)情普查L(zhǎng)CA為例,500平以下作為小面積) 1、如果最小地類(lèi)一致,則根據(jù)小地類(lèi)進(jìn)行融合(如0360跟0360融合) 2、如果沒(méi)有最小地類(lèi),則根據(jù)相同大的地類(lèi)進(jìn)行融合(0360跟0311融合) 3、如果大的地類(lèi)也不一致,則與最大面積相鄰面進(jìn)行融合。 思路: 1、暫時(shí)不會(huì)做工具箱,所以需要一步一步的做; 2、第一、二步的核心是根據(jù)地類(lèi)代碼來(lái)融合。假設(shè)小圖斑的數(shù)據(jù)為原始數(shù)據(jù),而相鄰大圖斑為目標(biāo)數(shù)據(jù),就是你把原始地類(lèi)代碼替換成目標(biāo)地類(lèi)代碼,同時(shí)應(yīng)保留一個(gè)目標(biāo)數(shù)據(jù)的唯一碼,防止目標(biāo)數(shù)據(jù)在融合過(guò)程中與其他相鄰地類(lèi)代碼進(jìn)行融合了; 3、根據(jù)第一二步的數(shù)據(jù)篩選出第三步的數(shù)據(jù)進(jìn)行eliminate,這里需要注意 你是要與相鄰地類(lèi)中線最長(zhǎng)的還是面積最大的融合。 4、過(guò)程中需要select load數(shù)據(jù) 所以需要耐心替換代碼根據(jù)spatialjoin的工具+python代碼來(lái)實(shí)現(xiàn) 電腦環(huán)境: win7,Arcmap10.3,ArcGIS內(nèi)置python2.7.8 步驟如下: 分離500平面積 注意:這里要在面積大于500的圖層中新建一個(gè)字段將OID掛接過(guò)去當(dāng)做唯一碼 通過(guò)spatialjoin獲取相鄰圖斑屬性這里要注意是用小面積掛接大面積圖層,另外注意新增兩個(gè)TEXT字段來(lái)獲取大面積圖層的信息以便判讀地類(lèi)編碼CC字段是否復(fù)合融合原則。工具運(yùn)行成功后,在輸出的圖層中新建兩個(gè)文本字段,融合CC碼及融合唯一碼 工具運(yùn)行成功后,在輸出的圖層中新建兩個(gè)文本字段,融合CC碼及融合唯一碼 到這里,屬性已經(jīng)獲取成功了,現(xiàn)在就是來(lái)判斷這些屬性跟面積小于500屬性的關(guān)系了,其中面積大于500圖斑的數(shù)量為9795-2127=7688條,最終數(shù)量需要與之進(jìn)行驗(yàn)證。 接下里分離復(fù)合第一二步原則的的圖形,為空值的即不符合按第三原則進(jìn)行融合,然后將符合第一二原則的圖形LOAD進(jìn)入面積大于500的圖層,其中要注意融合CC碼對(duì)應(yīng)面積大于500圖層的CC,融合唯一值面積對(duì)應(yīng)大于500的唯一碼。 進(jìn)行dissolve融合,按照唯一碼和CC,融合后的數(shù)量應(yīng)該與導(dǎo)入前一致(7668) 最后屬性掛接通過(guò)唯一碼字段將原來(lái)的屬性信息掛接到圖形數(shù)據(jù)中。你可以在access或者arcmap中進(jìn)行
版權(quán)聲明:本文為CSDN博主「GIS小菜鳥(niǎo)」 |
|
來(lái)自: huajiang0954 > 《林業(yè)》