三调与二调图斑叠加分析,筛选不同地类面积占比,筛选举证图斑

主要步骤:

  1. 标识数据,叠加分析,用标识,生成所有相交图斑,属性有原图斑的地类和国家的地类,以及原图斑的面积。,
  2. 生成的面域,增加4个字段(图斑的三调一级类,图斑的国家NYYPDL,是否相同,标识后的图斑面积),转换三调地类为二调的一级类,转换国家地类为二调一级类。
  3. 计算两个地类是否相同,相同就把是否相同这个字段填写相同,不同就填不同。
  4. 导出所有不同的细碎图斑,计算图斑面积等于SHAPE面积
  5. 根把原图斑号,融合细碎图斑,汇总细碎图斑的面积
  6. 添加一个字段,计算细碎图斑的面积占比。或者面积周长比等,用于进一步筛选,是否需要举证。

三调与二调图斑叠加分析,筛选不同地类面积占比,筛选举证图斑

处理三调图斑属性,删掉不要的字段,把三调图斑的图斑编号,用计算字段,计算一下,赋ID,这样不重复,用于最后的合并,增中一个三调图斑面积字段,计算等于SHAPE.AREA

三调与二调图斑叠加分析,筛选不同地类面积占比,筛选举证图斑

叠加分析,标识工具

三调与二调图斑叠加分析,筛选不同地类面积占比,筛选举证图斑

添加转换地类的字段,用字段计算器

三调与二调图斑叠加分析,筛选不同地类面积占比,筛选举证图斑

{"0101":"01","0102":"01","0103":"01","05H1":"20","0508":"20","06H1":"20","0601":"20","0602":"20","0701":"20","0702":"20","08H1":"20","08H2":"20","0808":"20","0809":"20","09":"20","1004":"20","1005":"20","1201":"20","0201":"02","0202":"02","0203":"02","0204":"02","0301":"03","0302":"03","0303":"03","0304":"03","0305":"03","0306":"03","0307":"03","0401":"04","0402":"04","0403":"04","0404":"04","1001":"10","1002":"10","1003":"10","1006":"10","1007":"10","1008":"10","1009":"10","1101":"SM","1102":"SM","1103":"SM","1104":"SM","1105":"SM","1106":"SM","1107":"SM","1109":"SM","1110":"SM","1201":"12","1202":"12","1203":"12","1204":"12","1108":"12","1205":"12","1206":"12","1207":"12"}[ !DLBM! ]

 

代码中括号内为三调转二调一级地类对应字典表,如果有缺少的,可以自己添加,

三调与二调图斑叠加分析,筛选不同地类面积占比,筛选举证图斑

{"0101":"01","0102":"01","0103":"01","05H1":"20","0508":"20","06H1":"20","0601":"20","0602":"20","0701":"20","0702":"20","08H1":"20","08H2":"20","0808":"20","0809":"20","09":"20","1004":"20","1005":"20","1201":"20","0201":"02","0202":"02","0203":"02","0204":"02","0301":"03","0302":"03","0303":"03","0304":"03","0305":"03","0306":"03","0307":"03","0401":"04","0402":"04","0403":"04","0404":"04","1001":"10","1002":"10","1003":"10","1006":"10","1007":"10","1008":"10","1009":"10","1101":"SM","1102":"SM","1103":"SM","1104":"SM","1105":"SM","1106":"SM","1107":"SM","1109":"SM","1110":"SM","1201":"12","1202":"12","1203":"12","1204":"12","1108":"12","1205":"12","1206":"12","1207":"12","011":"01","012":"01","013":"01","20":"20","021":"02","022":"02","023":"02","031":"03","032":"03","033":"03","041":"04","042":"04","043":"04","101":"10","102":"10","104":"10","105":"10","106":"10","107":"10","111":"SM","112":"SM","113":"SM","114":"SM","115":"SM","116":"SM","117":"SM","118":"SM","119":"SM","121":"12","122":"12","123":"12","124":"12","125":"12","126":"12","127":"12","CK":"CK","DT":"DT","JZ":"JZ","SM":"SM","YH":"YH","03/04":"03/04","01":"01","02":"02","03":"03","04":"04","203":"20","204":"20","1003":"10","1109":"SM",}[ !NYYPDL! ]

计算国家地类转换,我加了三调转二调一级的,和本身NYYPD出现的一些地类。如果还有缺少,可以手动添加字典表,国家NYYPD的填写,好像没有什么规律,如果报错了,可以看一下是缺少哪一个,手动添加就可以。

 

三调与二调图斑叠加分析,筛选不同地类面积占比,筛选举证图斑

计算是否相同字段,勾上代码块,定义函数,下面调用函数

def demo(dldm, gjdm):

   answer="不同"

   if dldm ==gjdm:

      answer="相同"

   return answer

demo( !SDDL_ZH! , !GJDL_ZH! )

 

粘贴时注意代码的缩进

 

HZMJ字段,赋值图形的面积,用于后面融合时汇总

 

导出所有不同的细碎图斑,如果沟渠和路不要的话,也可以去掉。用定义查询,然后导出SHP。

三调与二调图斑叠加分析,筛选不同地类面积占比,筛选举证图斑

导出后,对地类不同的细碎图斑进行融合,融合字段选tbbh,这样原来相同的图斑号的细碎图斑就可以重新融合,统计字段,把TBBH,统计FIRST,HZMJ,统计SUM。其它看情况需要。

 

三调与二调图斑叠加分析,筛选不同地类面积占比,筛选举证图斑

融合后,添加数据,添加字段,计算面积占比或者其它的。如果不需要计算的话,可以根据SUM_HZMJ字段,排序,根据不同地类的面积要求,来筛选需要举证的范围。