图斑地类面积理论值计算方法说明 下载本文

图斑地类面积理论值

计算方法说明

1 计算原则

1.1 面积平差原则

任何平差结果必须满足以下原则:

(一)涉及图形变化的图斑,分割或合并操作后,图斑地类面积总值应保持不变。 (二)图斑地类面积没有发生变化的,不允许参与平差。

1.2 面积平差方法

(一)将误差优先平差到图斑地类面积有变化的图斑。 (二)将面积残差以有效位数的最小值(面积为0.01),按参与变化的图斑地类面积数值,由大及小依次分配。

(三)当平差方法产生的结果和原则相抵触时,需调整平差方法以遵循原则要求。

1.3 面积计算方法

将地类图斑图层的图斑地类面积平差给临时数据,然后将临时数据的面积之和赋给调查成果数据。如下图(A1、A2为地类图斑图层的图斑,B1、B2、B3、B4为用于辅助计算的临时数据,C1、C2为调查成果中的图斑):

计算过程为:

1) A1、A2的边界将C1和C2打散形成临时数据B1、B2、B3、B4,再将地类图斑图层A1

的图斑地类面积按照B1、B2图形面积的比例分配,如果存在未分配的面积值,平差方法按照以0.01平均分配,余数按照面积由大及小排序分摊,得到平差后的B1、B2的图斑地类面积。

例如:B1所占比例为1/3、B2所占比例为2/3,将A1的面积的0.33分配给B1,0.66分配给B2,未分配的残差0.01分配给所占比例大的B2。

2) C1、C2的图斑地类面积分别为B1、 B3之和与B2、B4之和。

2 计算过程(含平差)

以HBTC为例,该图层图斑地类面积理论值计算的依据是调查底图(DLTB),下面介绍该图层理论值计算的过程。

说明:

耕地后备资源图斑:参考国家下发的调查底图,计算理论面积 不稳定耕地、新增耕地:参考二次调查地图图斑,计算理论面积

1)使用DLTB分割HBTC。采用的工具为arcgis中的Intersect,得到分割结果HBTC_Intersect。(DLTB为图2中的“地类图斑”、HBTC为图2中的“调查成果”、HBTC_Intersect为图3中的“临时数据”);

图1

图2 2)根据面积比例计算图斑地类面积(TBDLMJ)。利用HBTC_Intersect计算碎斑(B1、B2、B3、B4)的图斑地类面积(TBDLMJ)。A1、A2的边界将C1和C2打散形成HBTC_Intersect中的B1、B2、B3、B4,再将地类图斑图层A1的图斑地类面积(TBDLMJ)按照B1、B2图形面积(shapearea)的比例分配,得到碎斑(B1、B2)的图斑地类面积理论值(未平差)该

值暂时命名为TBDLMJ_BL:

B1.TBDLMJ_BL = round ( ( B1.shapearea / A1. Shapearea ) * A1.TBDLMJ , 2 )

3)平差。将HBTC_Intersect中相同DLTB_BSM的图斑对TBDLMJ_BL和shapearea进行求和,如:B1、B2,若abs ( ( sum(TBDLMJ_BL)-A1. TBDLMJ ) ) >= 0.01 and abs ( ( sum(shapearea)- A1. shapearea) )<0.01,则该DLTB图斑(A1)相关的碎面图斑(B1、B2)的TBDLMJ需要进行平差,其平差的残差CANCHA= sum(TBDLMJ_BL)- A1. TBDLMJ,将残差以0.01平均分配给相关碎斑(B1、B2)的TBDLMJ_BL,余数按照碎斑(B1、B2)面积由大及小排序分摊,得到平差后B1、B2的图斑地类面积(TBDLMJ_BL),同理计算B3、B4的图斑地类面积(TBDLMJ_BL)。

4)HBTC图斑的图斑地类面积理论值计算。HBTC_Intersect中按照HBTC_BSM对平差后的TBDLMJ_BL面积进行分组求和(B1、B3的TBDLMJ_BL求和即为C1的图斑地类面积理论值),即为HBTC图层中各图斑的图斑地类面积理论值。

上述过程描述质检软件平差的思路,从第3步开始,需要通过开发实现,若想初步验证图斑地类面积的准确性,可执行“3计算过程(不含平差)”的方法。

3 计算过程(不含平差)

说明:

下文详细描述了“2计算过程(含平差)”中第1步和第2步的过程,可得到未平差的图斑地类面积的理论值。通常情况下,一个图斑需要平差的数值不会很大,所以通过下面的方法得到的图斑地类面积与理论值的差异一般不会超出10平方米。

1) 在DLTB中创建DLTB_AREA字段(double类型),并将shapearea字段值赋值给该字段

2)使用DLTB分割HBTC。采用的工具为arcgis中的Intersect,得到分割结果HBTC_Intersect;

图1

3)得到叠加结果集HBTC_Intersect。

4)打开HBTC_Intersect属性表,创建两个字段BL(比例,double),BL_MJ(比例面积,double), 5)计算BL的值,即专题数据图斑所占地类图斑的几何面积比例,公式如下:[Shape_Area] / [DLTB_AREA]

6)计算BL_MJ,根据几何面积比例乘以地类图斑的TBDLMJ再保留两位小数得到的数值即为分割后的专题图斑的TBDLMJ未平差的值,公式如下:ROUND( [BL] * [TBDLMJ_1] ,2)

7)将HBTC_Intersect根据BSM,统计BL_MJ的值(即:将相同标识码的图斑进行合并,然后计算合并后的BL_MJ),得到结果集HBTC_Intersect_Dissolve1

8)SUM_BL_MJ字段即为HBTC中各图斑尚未平差的理论面积值。