原则:
优先级从上到下:(以地理国情普查LCA为例,500平以下作为小面积)
1、如果最小地类一致,则根据小地类进行融合(如0360跟0360融合)
2、如果没有最小地类,则根据相同大的地类进行融合(0360跟0311融合)
3、如果大的地类也不一致,则与最大面积相邻面进行融合。
思路:
1、暂时不会做工具箱,所以需要一步一步的做;
2、第一、二步的核心是根据地类代码来融合。假设小图斑的数据为原始数据,而相邻大图斑为目标数据,就是你把原始地类代码替换成目标地类代码,同时应保留一个目标数据的唯一码,防止目标数据在融合过程中与其他相邻地类代码进行融合了;
3、根据第一二步的数据筛选出第三步的数据进行eliminate,这里需要注意 你是要与相邻地类中线最长的还是面积最大的融合。
4、过程中需要select load数据 所以需要耐心替换代码根据spatialjoin的工具+python代码来实现
电脑环境:
win7,Arcmap10.3,ArcGIS内置python2.7.8
步骤如下:
分离500平面积
注意:这里要在面积大于500的图层中新建一个字段将OID挂接过去当做唯一码
通过spatialjoin获取相邻图斑属性
这里要注意是用小面积挂接大面积图层,另外注意新增两个TEXT字段来获取大面积图层的信息以便判读地类编码CC字段是否复合融合原则。工具运行成功后,在输出的图层中新建两个文本字段,融合CC码及融合唯一码
工具运行成功后,在输出的图层中新建两个文本字段,融合CC码及融合唯一码
到这里,属性已经获取成功了,现在就是来判断这些属性跟面积小于500属性的关系了,其中面积大于500图斑的数量为9795-2127=7688条,最终数量需要与之进行验证。
利用python字段计算器进行判断
代码如下:
def aaa(aa,bb,cc):
i= 0
for x in aa.split(","):
if x == cc:
return "".join(x)+"-"+bb.split(",")[i]
elif x[:3] == cc[:3]:
return "".join(x)+"-"+bb.split(",")[i]
elif x[:2] == cc[:2]:
return "".join(x)+"-"+bb.split(",")[i]
i += 1
计算完之后,可以看出其实复合第一第二原则的很少,而null值即为不符合原则的这里先分离唯一码和CC码,
代码如下:!融合CC码!.split("-")[1]
再计算CC码
代码如下:!融合CC码!.split("-")[0]
接下里分离复合第一二步原则的的图形,为空值的即不符合按第三原则进行融合,然后将符合第一二原则的图形LOAD进入面积大于500的图层,其中要注意融合CC码对应面积大于500图层的CC,融合唯一值面积对应大于500的唯一码。
导入大于500面积的图层中进行融合处理
进行dissolve融合,按照唯一码和CC,融合后的数量应该与导入前一致(7668)
效果如下图:
然后将符合原则三的图斑load进入融合后的图层进行面积Eliminate消除
到这里图形的数据已经处理完了,可以看出 还是7668条数据,接下来要进行属性挂接
最后属性挂接
通过唯一码字段将原来的属性信息挂接到图形数据中。你可以在access或者arcmap中进行