ArcGIS小图斑根据相邻地类属性融合。

原则：

优先级从上到下：（以地理国情普查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中进行