一种农地权属公示图自动生成方法

摘要

本发明公开了一种农地权属公示图自动生成方法，实现了在农村土地确权工作中承包地块权属公示图的自动化制作。目前多数专业软件对承包地块公示图的制作大都需要相应行政区划图层作为范围控制图层进行出图操作。但在实际工作中，承包地块公示范围往往是以小组为单位，多数地区行政区划数据大多精确至村边界，现有出图软件还未很好地解决范围控制图层缺少情况下难以自动出图这一问题。另外，由于出图打印设备限制，公示图的实际出图范围同样存在限制，在很多时候进行大比例尺公示图制作时需要对图层进行分幅处理，现有专业软件需要大量的人工操作来实现图层分幅，自动化程度不高。针对以上两点不足，本发明将计算机技术与地理信息技术结合，实现了农地权属公示图制作流程的自动化和智能化并提高了制作的效率和准确性，节约了大量的时间和人力成本。

说明书

技术领域

本发明涉及基础地理信息领域，具体是指一种高效的农地权属公示图自动生成方法。

背景技术

随着农村土地确权工作的不断推进，各类历史遗留的土地权属问题不断凸显。农村土地 所有权、使用权等权利混乱，没有得到明确规定和有效保障。农村土地确权工作旨在摸清土 地状况，明确土地权属，推动农村土地管理的规范化。在这个过程中，解决农村土地存在的 地块不实、四至不清、面积不准等问题是首要任务。在此基础上，面向农户制作农地权属公 示图是进行下一步工作的重要实施环节。农地权属公示图主要包括承包地块权属信息、图饰 元素、制图信息和农地公示信息统计表，其中统计表信息由承包方姓名、实测总面积、总数 量等信息组成。只有经过农户在农地权属公示图上的确认签字，农地外业调查成果才能作为 土地确权今后工作开展的基础。因此，农地权属公示图的制作作为土地确权工作中的一个实 施环节显得尤为重要。

现有专业软件在对农地权属公示图制作流程中，由于对公示范围控制和图层分幅处理等 方面支持的不足，导致无法实现制图的自动化和智能化。另外，针对公示信息统计表中记录 往往也是通过人工结合软件的半自动化方式进行图层权属信息的统计汇总。对于土地确权工 作中面对的海量承包地块数据，若还利用现有的制图方式制作农地权属公示图，将需投入大 量的人力和物力且还不能保障工作效率的提高和制图结果的准确。因此，在农地权属公示图 制作方面，特别是组单位行政区划数据获取、图层分幅和公示信息统计汇总方面，存在诸多 的不便且易产生错误。农地权属公示图制作流程无法实现自动化和智能化，这将给今后的土 地承包经营权管理和分析带来极大的不便。

发明内容

本发明的目的在于提供一种自动生成农地权属信息制作公示图的方法。

本发明通过下述技术方案实现：一种农地权属公示图自动生成方法，包括以下步骤：

(1)导入包括承包地块的空间位置、几何形状和地块标注信息的地图数据，选择承包地 块图层，通过对图层中各个要素的遍历操作获取该图层的最小外包矩形信息；

(2)利用图层最小外包矩形信息并结合公示比例尺标准，动态计算该图层的视图范围；

(3)将上一步计算出的视图范围的宽度与实际出图范围中宽度限值进行比较，判断图层 是否需要进行分幅处理；如果不需要分幅，直接输出作为公示数据源；否则将分幅处理后的 承包地块数据作为公示数据源；

(4)权属信息关联：基于上述操作所获公示数据源，借助于地理信息系统(GIS)空间 数据模型的数据组织功能，利用地块标注的关联作用，实现权属信息与承包地块的自动关联。

(5)当前图幅权属信息统计：以农户为统计单元，对当前图幅中农户所有相同权属信息 的地块进行统计汇总，将统计结果作为公示图层中权属统计表的数据源，并最终添加至公示 图中。

(6)智能化图层渲染与图饰添加操作：根据公示数据源中承包地块的几何形状，实现权 属信息的动态标注与可视化表达；实现对公示数据源的批量渲染；动态、自适应添加图饰信 息，包括图饰元素和制图信息；

(7)综合上述操作结果并根据用户选择的输出分辨率和输出格式，输出最终的承包地块 制图成果。

其中，最小外包矩形是指包围图元，且平行于坐标轴的最小外接矩形，在GIS和计算机 图形学中应用十分广泛。

为了更好的实现本发明所述的方法，进一步地，所述步骤(2)中，根据所获承包地块图 层最小外包矩形信息，可计算承包地块图层在X、Y轴方向范围分别为：XMax-XMin、 YMax-YMin。结合公示比例尺标准(如1:1000)动态计算该图层视图范围的宽度W1，W1＝ (XMax-XMin)/1000。

为了更好的实现本发明所述的方法，进一步地，所述步骤(3)中，将上一步计算出的视 图范围的宽度(记为W1)与实际出图范围中宽度限值(W)进行比较，计算Math.ceil(W1/W)， 结果记为n。若n<1则将该图层直接作为公示数据源。反之，则对图层进行分幅处理并保存 成新图层作为公示数据源。例如，实际出图范围宽度限值为91.4cm，承包地块图层在公示比 例尺下宽度计算为120.34cm。Math.ceil(W1/W)＝2，n>1，则其结果即为图层分幅数。接下来 需对承包地块图层进行分幅处理，将处理结果另存为新的承包地块图层，并将其作为公示数 据源。

为了更好的实现本发明所述的方法，进一步地，所述步骤(4)中，权属信息关联以所获 公示数据源为基础，借助于地理信息系统(GIS)空间数据模型在数据组织方面的优越性，以 承包地块中地块标注字段作为连接的主键，将在外部独立存储的权属信息文件导入承包地块 图层中，实现了承包地块图形数据和属性数据的一体化管理。

为了更好的实现本发明所述的方法，进一步地，所述步骤(5)中，对当前图幅权属信息 进行统计，动态汇总图层中相同权属地块的属性信息，包括承包地块的承包方姓名、实测总 面积、总数量。如在一张公示图上农户张三拥有三块地，面积分别为1.21、0.5、0.7，则统计 汇总后将形成一条数据记录，即农户张三拥有地块总面积为2.41、总数量为3。

为了更好的实现本发明所述的方法，进一步地，所述步骤(5)中，公示信息统计表的生 成算法如下：以图层中承包方编码字段为主键，对相同承包方编码(即同一农户)地块的权 属信息进行统计汇总；通过对图层中要素的遍历，将相同承包方编码的权属信息统计为一条 记录，对于存在承包方编码但承包方姓名异常的记录无公示意义，因此本方法中对于这种情 况记录采取不进行统计原则；图层中所有要素遍历完成后，将统计汇总结果存入临时数据集 中作为公示信息统计表的数据源。计算统计表中公示字段宽度可得单个统计表中表头宽度 W3，公示字段包括承包方姓名、实测总面积、地块数量、确认签字，结合公示比例尺下图层 视图宽度W4，计算表格布局列数c1＝Math.floor(W3/W4)，Math.floor方法是将计算结果向 下取整，如Math.floor(1.12)＝1；通过c1和统计汇总结果总数n1，计算表格布局行数 r1＝Math.ceil(n1/c1)，Math.ceil方法是将计算结果向上取整，如Math.ceil(1.12)＝2；公示 信息统计表设计完成后，将统计汇总结果写入表格中，并最终添加至公示图层。

为了更好的实现本发明所述的方法，进一步地，所述步骤(6)中，利用预定义的图层渲 染样式，包括地块填充色、边界线颜色和线型，实现了对公示数据源的批量渲染。操作过程 中可动态设置渲染参数，以期满足了不同公示要求。

为了更好的实现本发明所述的方法，进一步地，所述步骤(6)中，利用VBScript编写 的权属信息标注功能，实现了包括承包方姓名、地块编码、实测面积等权属信息的动态标注 和可视化表达，代码编写如下：

FunctionFindLabel([地块编码],[承包方姓名],[实测面积])

dima

a＝len([承包方姓名])

FindLabel＝[承包方姓名]&vbCrLf&string(int(a),"—")&space(1)&round([实测面积],2)&

vbCrLf&[地块编码]

EndFunction

其中，权属信息标注内容能够准确的添加至对应承包地块空间位置，公示过程中农户可 通过标注信息中的地块编码查看相应地块的详细属性信息。利用VBScript编写的权属标注代 码简单易懂，实际出图过程中可针对特殊情况，也可快速进行调整。

为了更好的实现本发明所述的方法，进一步地，所述步骤(6)中，基于上述操作结果， 可动态、自适应的添加图饰信息，自适应是指可根据出图要求调整图饰信息的在公示图中的 位置。其中图饰元素包括图例、指北针、比例尺等，图例信息由承包方姓名、地块编码和实 测面积信息组成，是通过读取图层中随机一条完整的属性记录获得；比例尺即为公示比例尺； 指北针用于提供方位信息。制图信息包括制图标题、制图人、审核员及制图单位等，制图标 题通过获取公示图层名动态生成并添加至公示图中；制图人、审核员和制图单位可通过动态 输入并自适应添加至公示图中。

为了更好的实现本发明所述的方法，进一步地，所述步骤(7)中，综合上述操作结果， 并根据用户选择选择的输出分辨率和输出格式，其中输出格式包括JPG、PDF、GIF、EPS等， 输出最终的承包地块制图成果。

本发明与现有技术相比，具有以下优点和有益效果：

本发明实现了在土地确权工作中，农地权属公示图的自动生成，克服了现有专业软件在 农地权属公示图制作中存在的不足。将计算机技术和地理信息技术结合，充分利用空间数据 模型在数据组织方面的优越性。通过遍历图层中所有要素获取的图层外包矩形信息，计算公 示比例尺下图层视图范围的宽度并与实际出图范围的宽度限值进行比较，实现自动分幅操作， 避免了打印设备对图幅宽度的限制；实现了公示数据源中权属信息的关联，利用主键关联字 段将勾绘的承包地块空间数据与独立存储的权属信息文件进行自动关联，满足了承包地块图 形数据和属性数据的一体化管理，同时图层中丰富的权属信息可为后期地块标注和权属信息 的统计汇总提供数据基础；以当前图幅为基础，对图层中权属信息进行动态统计并以表格形 式添加至公示图中，农户通过表格信息可高效快捷的查看自家所有地块总面积和数量信息并 进行确认签字；针对公示数据源数据量大的特点，采用外部预定义的渲染样式，对公示数据 源进行批量渲染利用简单易懂的轻量级脚本语言VBScript，实现了权属信息的动态标注。根 据出图要求，实现了图饰信息的动态、自适应添加。从承包地块数据处理到公示图制图成果 的输出，本方法借鉴了现有专业软件的优点，同时克服其存在的不足，经应用实践证明，显 著提高了公示图的制图效率，节约了大量的时间和人力成本。

附图说明

图1为发明方法流程示意图；

图2为本发明方法中图层分幅实现流程；

图3为本发明方法中公示信息统计表生成流程；

图4为本发明方法中地块标注与权属信息关联方法示意图；

图5.1为本发明未进行分幅处理的技术效果图；

图5.2.1和图5.2.2是分幅处理后的技术效果图。

具体实施方式

以下将结合附图和具体实施方式对本发明的一种农地权属公示图自动生成方法中的关键 模块进行详细描述。下面通过参考附件描述的实施例是示例性的，仅用于解释本发明，而不 能理解为对本发明的限制。

实施例1：

本实施例的具体实现一种农地权属公示图自动生成方法，具体流程图如图1所示。

步骤(1)导入包括承包地块的空间位置、几何形状和地块标注信息的地图数据，选择承 包地块图层，通过对图层中各个要素的遍历操作获取其最小外包矩形信息，包括XMax、XMin、 YMax、YMin的值，其中XMax为所述最小外包矩形在X轴的最大值，XMin为所述最小外 包矩形在X轴范围的最小值，YMax为所述最小外包矩形在Y轴范围的最大值，YMin为所 述最小外包矩形在Y轴范围的最小值。计算所获取的最小外包矩形在X、Y轴方向范围，分 别为：XMax-XMin、YMax-YMin。

步骤(2)结合图层公示比例尺标准，动态计算该图层视图范围(W1,H1)，W1为视图宽 度，H1为视图高度。

步骤(3)将上一步计算出的视图范围的宽度(记为W1)与实际出图范围中宽度限值W 进行比较，计算Math.ceil(W1/W)，结果记为n，此方法是将计算结果向上取整，如计算值 为Math.ceil(0.12)，则n＝1。若n<1则将该图层直接作为公示数据源。反之，则需进行图层 分幅处理，分幅过程中根据图层视图宽度构建要素选择框实现图层的分幅流程，并保存成新 图层发送至公示数据源。

步骤(4)权属信息关联：如图4所示，基于上步骤所获公示数据源，借助于地理信息系 统(GIS)的空间数据模型的数据组织功能，利用地块标注的关联作用，实现包括承包方姓名、 地块编码、实测面积等权属信息与承包地块空间信息的自动关联，以承包地块中地块标注作 为连接的主键将独立存储的权属信息文件导入承包地块图层中；对当前图幅权属信息统计， 以图层中所设主键字段为标准实现关联。

步骤(5)当前图幅权属信息统计：动态汇总图层中相同权属地块的属性信息，包括承包 地块的承包方姓名、实测总面积、总数量，作为公示数据源中公示信息统计表的数据来源。

步骤(6)智能化图层渲染与图饰添加操作：根据公示数据源中承包地块的几何形状，利 用VBJavaScript编写的权属信息标注功能，实现了权属信息的动态标注与可视化表达，其中 权属信息能够准确的添加至相应承包地块空间位置，并可通过标注中的地块编码查看该地块 的详细属性信息；实现了对公示数据源的批量渲染，包括统一设置图层中地块填充色、边界 线颜色和线型等。

利用VBScript编写的权属信息标注功能，实现包括承包方姓名、地块编码、实测面积等 权属信息的动态标注和可视化表达，如下。

FunctionFindLabel([地块编码],[承包方姓名],[实测面积])

dima

a＝len([承包方姓名])

FindLabel＝[承包方姓名]&vbCrLf&string(int(a),"—")&space(1)&round([实测面积],2)&

vbCrLf&[地块编码]

EndFunction

其中，权属信息标注内容能够准确的添加至对应承包地块空间位置，公示过程中农户可 通过标注信息中的地块编码查看相应地块的详细属性信息。利用VBScript编写的权属标注代 码简单易懂，实际出图过程中可针对特殊情况，也可快速进行调整。

根据出图要求和权属信息统计结果，可以自动添加图饰元素和制图信息，其中图饰元素 包括图例、指北针、比例尺等，图例信息由承包方姓名、地块编码和实测面积组成，是通过 读取图层中随机一条完整的属性记录获得；比例尺即为公示比例尺；指北针用于提供方位信 息。制图信息包括制图标题、制图人、审核员及制图单位等，制图标题通过获取公示图层名 动态生成并添加至公示图中；制图标题根据承包地块图层名自动添加，制图人、审核员和制 图单位可可通过动态输入并自适应添加至公示图中。步骤(7)将输出的图层根据用户选择的 输出分辨率和输出格式，输出最终的承包地块制图成果，输出格式包括JPG、PDF、GIF、EPS。

图2为图层分幅实现流程，其中首先需要确定承包地块图层；通过对图层中各个要素的 遍历操作获取其最小外包矩形信息，包括XMax、XMin、YMax、YMin的值，并计算其在X、 Y轴方向范围分别为：XMax-XMin(记为layerWidth)、YMax-YMin。结合图层公示比例尺 标准(假设为1：1000)，动态计算该图层视图范围的宽度(记为W1)：1:1000＝W1/layerWidth； 将上一步计算出的视图范围的宽度W1与实际出图范围中宽度限值(记为W)进行比较，计 算Math.ceil(W1/W)，计算结果记为n，此方法是将计算结果向上取整，如计算值为Math.ceil (0.12)，则n＝1；上述计算结果若不大于1，则将其直接作为公示数据。反之，则以n作为 图层分幅数，制定图层分幅策略；本文中图层分幅原理为以承包地块图层最小外包矩形中 XMin值为X轴方向的起点，YMax-YMin的值为固定高度，构建数据选择框，X轴方向起始 终点为layerWidth/n，记为Wx。循环分幅过程中，以上一次X轴方向终点作为下一次分幅X 轴方向起点，依次构建数据选择框。其中，将每次选择的承包地块要素，保存成新图层并输 出为公示数据源，直到分幅完成。

图3为公示信息统计表生成流程，以图层中承包方编码字段为主键，对相同承包方编码， 即同一农户，对地块的权属信息进行统计汇总，图层中所有要素遍历完成后，将统计汇总结 果作为公示信息统计表的数据源，汇总后的数据记录包括承包地块的承包方姓名、实测总面 积、总数量。如在一张公示图上农户张三拥有三块地，面积分别为1.21、0.5、0.7，则统计汇 总后将形成一条数据记录，即农户张三拥有地块总面积为2.41、总数量为3，如图3中公示 信息统计表所示。将统计记录读取至数据集中进行保存，基于统计记录进行公示信息统计表 的制作。首先通过公示字段(承包方姓名、实测总面积、地块数量、确认签字)，计算表头宽 度，记为W3，结合公示比例尺下图层宽度(记为W4)，计算表格布局列数 c1＝Math.floor(w3/w4)，此方法是将计算结果向下取整，如Math.floor(1.12)＝1；通过统计记 录总数(记为n1)并结合上一步中所计算的布局列数，计算表格布局行数r1＝Math.ceil(n1/c1)， 此方法是将计算结果向上取整，如Math.ceil(1.12)＝2；公示信息统计表设计完成后，则可 将之前数据集中统计汇总数据写入设计表格中，并最终将表格添加至公示图层。

最终公示效果如图5.1，图5.2.1和图5.2.2所示。