一种农地权属公示图自动生成方法
【技术领域】
[0001]本发明涉及基础地理信息领域，具体是指一种高效的农地权属公示图自动生成方法。
【背景技术】
[0002]随着农村土地确权工作的不断推进，各类历史遗留的土地权属问题不断凸显。农村土地所有权、使用权等权利混乱，没有得到明确规定和有效保障。农村土地确权工作旨在摸清土地状况，明确土地权属，推动农村土地管理的规范化。在这个过程中，解决农村土地存在的地块不实、四至不清、面积不准等问题是首要任务。在此基础上，面向农户制作农地权属公示图是进行下一步工作的重要实施环节。农地权属公示图主要包括承包地块权属信息、图饰元素、制图信息和农地公示信息统计表，其中统计表信息由承包方姓名、实测总面积、总数量等信息组成。只有经过农户在农地权属公示图上的确认签字，农地外业调查成果才能作为土地确权今后工作开展的基础。因此，农地权属公示图的制作作为土地确权工作中的一个实施环节显得尤为重要。
[0003]现有专业软件在对农地权属公示图制作流程中，由于对公示范围控制和图层分幅处理等方面支持的不足，导致无法实现制图的自动化和智能化。另外，针对公示信息统计表中记录往往也是通过人工结合软件的半自动化方式进行图层权属信息的统计汇总。对于土地确权工作中面对的海量承包地块数据，若还利用现有的制图方式制作农地权属公示图，将需投入大量的人力和物力且还不能保障工作效率的提高和制图结果的准确。因此，在农地权属公示图制作方面，特别是组单位行政区划数据获取、图层分幅和公示信息统计汇总方面，存在诸多的不便且易产生错误。农地权属公示图制作流程无法实现自动化和智能化，这将给今后的土地承包经营权管理和分析带来极大的不便。
【发明内容】
[0004]本发明的目的在于提供一种自动生成农地权属信息制作公示图的方法。
[0005]本发明通过下述技术方案实现:一种农地权属公示图自动生成方法，包括以下步骤:
[0006](I)导入包括承包地块的空间位置、几何形状和地块标注信息的地图数据，选择承包地块图层，通过对图层中各个要素的遍历操作获取该图层的最小外包矩形信息；
[0007](2)利用图层最小外包矩形信息并结合公示比例尺标准，动态计算该图层的视图范围；
[0008](3)将上一步计算出的视图范围的宽度与实际出图范围中宽度限值进行比较，判断图层是否需要进行分幅处理;如果不需要分幅，直接输出作为公示数据源;否则将分幅处理后的承包地块数据作为公示数据源；
[0009](4)权属信息关联:基于上述操作所获公示数据源，借助于地理信息系统(GIS)空间数据模型的数据组织功能，利用地块标注的关联作用，实现权属信息与承包地块的自动关联。
[0010](5)当前图幅权属信息统计:以农户为统计单元，对当前图幅中农户所有相同权属信息的地块进行统计汇总，将统计结果作为公示图层中权属统计表的数据源，并最终添加至公示图中。
[0011](6)智能化图层渲染与图饰添加操作:根据公示数据源中承包地块的几何形状，实现权属信息的动态标注与可视化表达;实现对公示数据源的批量渲染;动态、自适应添加图饰信息，包括图饰元素和制图信息；
[0012](7)综合上述操作结果并根据用户选择的输出分辨率和输出格式，输出最终的承包地块制图成果。
[0013]其中，最小外包矩形是指包围图元，且平行于坐标轴的最小外接矩形，在GIS和计算机图形学中应用十分广泛。
[0014]为了更好的实现本发明所述的方法，进一步地，所述步骤(2)中，根据所获承包地块图层最小外包矩形信息，可计算承包地块图层在X、Y轴方向范围分别为:XMax-XMin、YMax-YMin。结合公示比例尺标准(如1:1000)动态计算该图层视图范围的宽度11，￥1 =(XMax-XMin)/1000o
[0015]为了更好的实现本发明所述的方法，进一步地，所述步骤(3)中，将上一步计算出的视图范围的宽度(记为Wl)与实际出图范围中宽度限值(W)进行比较，计算Math.ceiKWl/W)，结果记为η。若n〈l则将该图层直接作为公示数据源。反之，则对图层进行分幅处理并保存成新图层作为公示数据源。例如，实际出图范围宽度限值为91.4cm，承包地块图层在公示比例尺下宽度计算为I^O^cmHceiKWlVW) = ],]^,则其结果即为图层分幅数。接下来需对承包地块图层进行分幅处理，将处理结果另存为新的承包地块图层，并将其作为公示数据源。
[0016]为了更好的实现本发明所述的方法，进一步地，所述步骤(4)中，权属信息关联以所获公示数据源为基础，借助于地理信息系统(GIS)空间数据模型在数据组织方面的优越性，以承包地块中地块标注字段作为连接的主键，将在外部独立存储的权属信息文件导入承包地块图层中，实现了承包地块图形数据和属性数据的一体化管理。
[0017]为了更好的实现本发明所述的方法，进一步地，所述步骤(5)中，对当前图幅权属信息进行统计，动态汇总图层中相同权属地块的属性信息，包括承包地块的承包方姓名、实测总面积、总数量。如在一张公示图上农户张三拥有三块地，面积分别为1.21、0.5、0.7，则统计汇总后将形成一条数据记录，即农户张三拥有地块总面积为2.41、总数量为3。
[0018]为了更好的实现本发明所述的方法，进一步地，所述步骤(5)中，公示信息统计表的生成算法如下:以图层中承包方编码字段为主键，对相同承包方编码(即同一农户)地块的权属信息进行统计汇总;通过对图层中要素的遍历，将相同承包方编码的权属信息统计为一条记录，对于存在承包方编码但承包方姓名异常的记录无公示意义，因此本方法中对于这种情况记录采取不进行统计原则；图层中所有要素遍历完成后，将统计汇总结果存入临时数据集中作为公示信息统计表的数据源。计算统计表中公示字段宽度可得单个统计表中表头宽度W3，公示字段包括承包方姓名、实测总面积、地块数量、确认签字，结合公示比例尺下图层视图宽度W4，计算表格布局列数cl =Math.floor(W3/W4) ,Math.floor方法是将计算结果向下取整，如Math.floor(l.12) = 1;通过cl和统计汇总结果总数nl，计算表格布局行数rl =Math.ceil(nl/cI)，Math.ceil方法是将计算结果向上取整，如Math.ceil(1.12)=2;公示信息统计表设计完成后，将统计汇总结果写入表格中，并最终添加至公示图层。
[0019]为了更好的实现本发明所述的方法，进一步地，所述步骤(6)中，利用预定义的图层渲染样式，包括地块填充色、边界线颜色和线型，实现了对公示数据源的批量渲染。操作过程中可动态设置渲染参数，以期满足了不同公示要求。
[0020]为了更好的实现本发明所述的方法，进一步地，所述步骤(6)中，利用VBScript编写的权属信息标注功能，实现了包括承包方姓名、地块编码、实测面积等权属信息的动态标注和可视化表达，代码编写如下:
[0021]Funct1n FindLabel([地块编码]，[承包方姓名]，[实测面积])
[0022]dim a
[0023]a = len([承包方姓名])
[0024]FindLabel =[承包方姓名]&vbCrLf&string(int(a),〃一")&space(I)&round([实测面积]，2)&
[0025]vbCrLf&[地块编码]
[0026]End Funct1n
[0027]其中，权属信息标注内容能够准确的添加至对应承包地块空间位置，公示过程中农户可通过标注信息中的地块编码查