建设用地占用永久基本农田比对算法

建设用地占用永久基本农田比对算法

一、建设用地占用永久基本农田比对算法
采用逐基本农田保护图斑计算的方法进行计算，数据精度和计算容差均为数据的默认精度和容差，具体算法及计算顺序如下：
1.建设项目图斑占用永久基本农田面积=（相交部分图形面积/基本农田保护图斑图形面积）* 基本农田保护图斑属性面积（即JBNTMJ）。
2.建设项目用地图斑占用永久基本农田保护图斑面积等于建设项目占用每块永久基本农田保护图斑面积之和。注：若占用单个永久基本农田保护图斑面积小于1平米，不参与汇总；
具体计算流程示意图如下：

如图所示：假设永久基本农田保护图斑A、B的图形面积（Shape_Are）均为100㎡，A、B两图斑的基本农田面积（JBNTMJ）均为80㎡，建设用地项目图斑M与图斑A的重叠图形面积为1㎡，与图斑B的重叠图形面积为30㎡；

计算建设用地项目图斑M占用永久基本农田保护图斑面积具体计算方法如下：
1.建设项目图斑M占用永久基本农田保护图斑B的面积=（相交部分图形面积/基本农田保护图斑图形面积）* 基本农田保护图斑属性面积（即JBNTMJ）；即：建设项目图斑占用永久基本农田保护图斑B的面积=（30/100）*80=24㎡;
2.同理计算得到建设项目图斑M占用永久基本农田保护图斑A的面积=(1/100)*80=0.8㎡；（由于占用面积小于1㎡，故舍去。）
3.建设项目图斑M占用永久基本农田保护图斑面积=24㎡；

二、补划基本农田占用现状比对算法
采用逐地类图斑计算的方法进行计算，具体算法及计算顺序如下：
1.单图斑内零星地物面积=（范围图斑与地类图斑重叠范围所占LXDW要素）LXDWMJ
2.单图斑内线状地物面积=（范围图斑与地类图斑重叠范围所占XZDW要素）XZDWMJ*（所占长度shape_len/原长度shape_len）(计算范围所占XZDW要素)KCBL
3.单图斑内扣除地类面积 =【（范围图斑与地类图斑重叠范围所占DLTB要素）TBMJ(所占面积shape_Area/原图斑shape_Area)- 零星地物面积（步骤1结果）- 线状地物面积（步骤2结果）】（范围图斑与地类图斑重叠范围所占DLTB要素）TKXS
4.单图斑内地类图斑面积=（范围图斑与地类图斑重叠范围所占DLTB要素）TBMJ(所占面积shape_Area/原图斑shape_Area)-零星地物面积（步骤1结果）- 线状地物面积（步骤2结果）-扣除地类面积（步骤3结果）
5.计算范围所占总面积=零星地物面积+线状地物面积+扣除地类面积+地类图斑面积
6.计算相应的地类面积可按地类编码分组统计。
具体计算流程示意如下：

如图所示，已知地类图斑A、B的图形面积(shape_area)均为100，B图斑图斑面积(TBMJ)为100，范围图层与B图斑重叠图形面积为30，B图斑田坎系数（TKXS）为0.05，零星地物面积为5，线状地物面积(XZDWMJ)为30，线状地物长度(shape_len)为20，在范围图斑内的长度为10，由于线状地物位于A、B两图斑交界处，故扣除比例为0.5;

计算范围图层M占地类图斑B的地类图斑面积计算流程如下：
1.B图斑内零星地物面积=（范围图斑与地类图斑重叠范围所占LXDW要素）LXDWMJ ；即：零星地物面积=5;
2.B图斑内线状地物面积=（范围图斑与地类图斑重叠范围所占XZDW要素）XZDWMJ*（所占长度shape_len/原长度shape_len）(范围图斑与地类图斑重叠范围所占XZDW要素)KCBL；即：线状地物面积=30（10/20）0.5=7.5;
3.B图斑内扣除地类面积 =【（范围图斑与地类图斑重叠范围所占DLTB要素）TBMJ(所占面积shape_Area/原图斑shape_Area)- 零星地物面积（步骤1结果）- 线状地物面积（步骤2结果）】（范围图斑与地类图斑重叠范围所占DLTB要素）TKXS;即：扣除地类面积=[100(30/100）-5-7.5]0.05=0.875；
4.B图斑内地类图斑面积=（范围图斑与地类图斑重叠范围所占DLTB要素）TBMJ(所占面积shape_Area/原图斑shape_Area)-零星地物面积（步骤1结果）- 线状地物面积（步骤2结果）-扣除地类面积（步骤3结果）;即：B图斑内地类图斑面积=（100*(30/100））-5-7.5-0.875=16.625；
同理计算得到A图斑内地类图斑面积
求范围图层内占用各地类面积则将范围内所占各图斑内地类图斑面积相加得到。