一种利用公共网络地图辅助土地调查的方法

摘要

本发明公开了一种利用公共网络地图辅助土地调查的方法，当在屏幕上显示较大区域时，系统自动调用公共网络地图(如谷歌、百度、天地图)的网络地图或遥感图作为调查底图；当在屏幕显示被调查区域较小范围或在土地调查测量地物点坐标时，选用预装的当地土地利用现状数据利用R-TREE空间索引快速调入相应区域的底图。本方法的优点在于采用移动测量设备(如GPS/PDA一体机)进行较大区域土地调查时，不需要对土地利用现状数据作复杂、费时、技术门槛高的自动综合，具有快速部署、使用方便的优点。

说明书

技术领域

本发明属于计算机信息技术领域，涉及地理信息科学。

背景技术

随着第二次全国土地调查的顺利完成，年度的土地变更调查逐渐提上日程。年度的 土地变更调查是对土地调查数据库维护的关键，土地变更调查是县级国土资源管理部门 根据下发的上时点土地利用现状数据和影像数据，对现状土地利用情况和最新影像进行 外业实地调查，获得变化地类图斑、权属(宗地)、行政区划数据，从而生成增量数据 包，实现对土地调查数据库进行更新的过程。

外业调查是土地更新调查的重要环节，传统的方式往往是携带纸质土地利用现状 图、影像、表格到现场手工填写，利用皮尺和测量仪器进行野外测量、相机拍照，操作 复杂，内业整理量大，效率低，难以满足管理要求。近年来，随着3S技术的发展，以 手持GPS/PDA为代表的移动测量设备逐渐应用于土地调查外业实地调查工作中并成为外 业土地利用现状调查的主流设备。但是，受限于移动GIS平台硬件的处理能力和存储空 间，必须使用配套的数据预处理软件对数据进行自动综合、建立空间索引，过程繁琐复 杂，不利于推广应用。采用移动测量设备进行土地调查时，目前的做法一般用土地利用 现状或者遥感图做调查底图，而土地利用现状一般是预装在移动测量设备的。

土地利用现状数据是一种矢量结构的数据，目前在移动测量设备上显示大规模土地 利用现状数据，为了同时支持大区域全貌和小区域局部底图的快速显示，必须对土地利 用现状数据进行预处理。由于土地利用现状数据中最主要的数据是面状地物，而这些面 状地物具有致密覆盖，不存在等级属性的特点，要实现快速LOD显示，通常需要对底图 数据进行如下预处理，主要包括：

(1)对土地利用现状进行多个级别的自动综合

(2)对每个级别综合处的土地利用现状数据建立空间索引

而土地利用现状数据的致密覆盖和无等级特性使得其综合算法缺乏可靠的依据，计 算复杂，处理起来费时费力，土地利用现状数据综合一直缺乏好的解决方案，技术门槛 高，从而限制了土地利用现状矢量数据在土地调查中的应用，尤其是大范围土地调查。

由此，避免土地利用现状自动综合又能在不同比例下均能给调查人员快速、忠实展 示调查底图至关重要。

发明内容

发明目的：针对上述现有技术，提出一种利用公共网络地图辅助土地调查的方法， 能够快速支持大区域全貌和小区域局部底图的快速显示。

技术方案：一种利用公共网络地图辅助土地调查的方法，包括如下步骤：

步骤1)，移动测量设备中土地利用现状数据的预处理：把土地利用现状数据用嵌入 式数据库sqlite存储，并对土地利用现状数据添加R-tree索引；

步骤2)，获取移动测量设备的土地利用现状数据显示最小比例S_min；

步骤3)，在移动测量设备每次更新显示调查底图时，当移动测量设备屏幕对应显示 范围大于或等于由所述土地利用现状数据显示最小比例S_min对应的地理区域时，通过所 述R-tree索引调入移动测量设备本地的土地利用现状数据作为调查底图；否则，直接调 用远程公共底图服务器上的公共地图或者遥感地图作为移动测量设备中土地调查的底 图。

进一步的，所述步骤2)中，移动测量设备的土地利用现状数据显示最小比例S_min由用户预先设定，并存储到移动测量设备的配置文件中。

进一步的，当所述移动测量设备不能读取配置文件中的土地利用现状数据显示最小 比例S_min时，移动测量设备通过自适应方式自动获取所述S_min值，包括如下步骤：

假设第i次调整移动测量设备显示范围时，显示的土地利用现状数据大小为东西长 为L_i，南北长为L_i*H/W；初始时，i取0；

步骤a)，在显示的土地利用现状数据上随机选取P个点，以所述P个点中的每个点为 中心，对东西长L_i，南北长L_i*H/W的土地利用现状数据，逐个进行空间地物检索并将空 间地物调入内存后显示在移动测量设备的屏幕上，同时统计每个点对应区域的空间地物 调入内存至屏幕显示的时间，得到的最长时间记为t_max；

步骤b)，假设每次屏幕进行放大、缩小、移动、全局显示操作时，地图更新时间的 最大允许显示延迟为T_max，比较所述步骤a)得到的t_max和T_max：如果t_max＜T_max，按比例系数 k扩大以所述每个点为中心的显示区域，即L_i+1＝L_i*k，H_i+1＝H_i*k，其中k取值范 围为1.4-2.0；

步骤c)，重复步骤a)至步骤b)，直至第k次调整时t_max＞T_max，则取前一次对应的显示 比例作为S_min，，S_min＝L_k-1/W。

进一步的，所述步骤3)中调用远程公共底图服务器上的公共地图或者遥感地图作 为土地调查的底图时采用缓冲机制，包括如下具体步骤：

步骤a)，设置缓冲区域中固定缓冲地图块数目为M，并且对每块缓冲地图块设置淘 汰指数N_t＝0；

步骤b)，检索缓冲区域是否存在当前要显示的缓冲地图块：如果存在，则将缓冲地 图块调入内存并显示在屏幕，同时将该缓冲地图块对应的淘汰指数N_t增1；如果不存在， 则从远程公共底图服务器上下载要显示的公共地图块或者遥感地图块到本地，如果此时 缓冲区未满，则将公共地图块或者遥感地图块作为缓冲地图块添加到缓冲区，并对加添 的缓冲地图块设置淘汰指数N_t＝0；如果缓冲区已满，则将下载的公共地图块或者遥感地 图块替换缓冲区中淘汰指数N_t最小值对应的缓冲地图块，并设置替换后的缓冲地图块的 淘汰指数N_t＝0。

进一步的，还包括通过用户界面选择公共地图种类的步骤，从远程公共底图服务器 上下载公共地图块时，所述公共地图包括谷歌地图、百度地图或天地图。

进一步的，所述步骤3)中，通过所述R-tree索引调入本地的土地利用现状数据作为 调查底图并显示的时间t大于T_max，则依据S_min＝S_min*sqrt(T_max/t)更新所述S_min。

有益效果：本发明提出的一种利用公共网络地图辅助土地调查的方法，在移动测量 设备每次更新显示调查底图时，当移动测量设备屏幕对应显示范围大于或等于由土地利 用现状数据显示最小比例S_min对应的地理区域时，即移动设备需要侧重于微观层面的具 体土地利用信息时，调用本地的土地利用现状数据实时显示，局部区域的土地利用现状 细节的展示对土地调查是必要的，并且对土地利用现状不需要做自动综合。当移动测量 设备屏幕对应显示范围小于由土地利用现状数据显示最小比例S_min对应的地理区域时， 即移动测量设备侧重显示的是宏观的信息时，由于公共地图比土地利用现状数据的引导 功能更强，此时屏幕显示内容用于引导工作人员快速到达现场的效果好于土地利用现状 数据。因此，采用移动测量设备进行较大区域土地调查时，不需要对土地利用现状数据 作复杂、费时、技术门槛高的自动综合，具有快速部署、使用方便的优点。

土地利用现状数据显示最小比例S_min的获取可采用自适应方式自动获取方式。获取 方式创造性的通过随机统计土地利用现状数据上一定区域的空间地物调入内存至屏幕 显示的时间，将该时间与地图更新时间的最大允许显示延迟为T_max作比较后，通过迭代 方法，逐次增加统计区域大小直至获取合适的S_min值，此时的S_min值是在空间地物调入 内存至屏幕显示的时间不超过T_max前提下能够最大显示土地利用现状数据底图。

调用远程公共底图服务器上的公共地图或者遥感地图作为土地调查的底图时采用 缓冲机制，通过对缓冲区中的缓冲地图块设置淘汰指数，将新下载的公共地图块或者遥 感地图块替换缓冲区中淘汰指数最小值对应的缓冲地图块，在实现下载地图块存储功能 的同时又保证了缓冲区中的地图块是调查最常用到的，从而保证之后使用时能够最快显 示所需底图。

附图说明

图一：一种利用公共网络地图辅助土地调查的方法流程图；

图二：自适应获取显示土地利用现状数据显示最小比例的方法流程图；

图三：调用远程底图的缓冲流程图。

具体实施方式

下面结合附图对本发明做更进一步的解释。

如图1所示，一种利用公共网络地图辅助土地调查的方法，包括如下步骤：

步骤1)，移动测量设备中土地利用现状数据的预处理：

把土地利用现状数据用嵌入式数据库sqlite存储，并对土地利用现状数据中的空 间地物添加R-tree索引。如果原始土地利用现状数据为shape或VCT等其他格式，需要 预先将之转化为sqlite格式；在转换的过程中，统计整个土地利用现状数的面积S和 该范围内空间地物的个数N，计算得到单位面积内平均空间地物为N/S。

步骤2)，获取移动测量设备的土地利用现状数据显示最小比例S_min：

土地利用现状数据显示最小比例S_min由用户预先设定，并存储到移动测量设备的配 置文件中，一般初始显示地理区域的宽度2000m～3000m，如果屏幕宽度W为480，那么 S_min取值设定在4.17～6.25之间。当移动测量设备不能读取配置文件中的土地利用现状数 据显示最小比例S_min，即配置文件中不存在S_min信息或被损坏时，如图2所示，移动测 量设备在程序初次运行时通过自适应方式自动获取S_min值，包括如下步骤：

假设移动测量设备初始显示的土地利用现状数据大小为东西长L₀，南北长L₀*H/W， L₀可以依据经验设定(如1000m)，其中W为移动测量设备屏幕的宽度，H为移动测量设备 屏幕的高度，后续第一次调整显示的土地利用现状数据大小后，东西长表示为L₁，南北长 表示为L₁*H/W，第二次调整后东西长表示为L₂，南北长表示为L₂*H/W，以此类推第i次调 整时东西长表示为L_i，南北长表示为L_i*H/W。

步骤a)，在显示的土地利用现状数据上随机选取P个点，P一般取5。以P个点中的每 个点为中心。假设第i次时，对东西长L_i，南北长L_i*H/W的土地利用现状数据，逐个进行 空间地物检索并将空间地物调入内存后显示在移动测量设备的屏幕上，同时统计每个点 对应区域的空间地物调入内存至屏幕显示的时间，得到的最长时间记为t_max；

步骤b)，假设每次屏幕进行放大、缩小、移动、全局显示操作时，地图更新时间的 最大允许显示延迟为T_max，比较步骤a)得到的t_max和T_max：如果t_max＜T_max，按比例系数k扩 大以上述每个点为中心的显示区域，即L_i+1＝L_i*k，H_i+1＝H_i*k，其中k取值范围为 1.4-2.0；

步骤c)，重复步骤a)至步骤b)，直至第k次调整时t_max＞T_max，则取前一次对应的显示 比例作为S_min，S_min＝L_k-1/W。

需要说明的是，通过用户预设的设置L_i初始取值以及最大允许显示延迟T_max值，初次 计算得到的t_max应小于T_max。

步骤3)，在移动测量设备每次更新显示调查底图时，计算当前移动测量设备的地图 显示比例S，S＝Lg/W，其中W为移动测量设备屏幕宽度，Lg是屏幕宽度为W时的地理长 度。当移动测量设备屏幕对应显示范围大于或等于由土地利用现状数据显示最小比例 S_min对应的地理区域，即S＞S_min时，通过R-tree索引调入移动测量设备本地的土地利用 现状数据作为调查底图；否则，直接调用远程公共底图服务器上的公共地图或者遥感地 图作为移动测量设备中土地调查的底图，并在移动测量设备屏幕上显示。如图3所示， 调用远程公共底图服务器上的公共地图或者遥感地图作为土地调查的底图时采用缓冲 机制，包括如下具体步骤：

步骤a)，设置缓冲区域中固定缓冲地图块数目为M，并且对每块缓冲地图块设置淘 汰指数N_t＝0；

步骤b)，检索缓冲区域是否存在当前要显示的缓冲地图块：如果存在，则将缓冲地 图块调入内存并显示在屏幕，同时将该缓冲地图块对应的淘汰指数N_t增1；如果不存在， 则采用HTTP等协议从远程公共底图服务器上下载要显示的公共地图块或者遥感地图块 到本地，如果此时缓冲区未满，则将公共地图块或者遥感地图块作为缓冲地图块添加到 缓冲区，并对加添的缓冲地图块设置淘汰指数N_t＝0；如果缓冲区已满，则将下载的公共 地图块或者遥感地图块替换缓冲区中淘汰指数N_t最小值对应的缓冲地图块，并设置替换 后的缓冲地图块的淘汰指数N_t＝0。其中，从远程公共底图服务器上下载公共地图块时， 可通过移动测量设备的用户界面选择公共地图种类的，公共地图包括谷歌地图、百度地 图或天地图。

步骤3)中，当通过R-tree索引调入本地的土地利用现状数据作为调查底图并显示的 时间t大于T_max时，则依据S_min＝S_min*sqrt(T_max/t)更新上述S_min值。此时显示最小比例 S_min值变小，移动测量设备能够较快的显示所需底图，即通过更新配置文件中的S_min值， 确保底图更新速度。当移动测量设备屏幕根据实际需要显示公共地图或土地利用现状数 据后，用户可以进行与移动测量设备的人机交互，在下一次进行屏幕进行放大、缩小、 移动、全局显示操作时等更新显示操作时，重新进行步骤3)。

以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员 来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也 应视为本发明的保护范围。