本发明涉及空间数据处理技术领域,尤其涉及一种基于遥感的矢量电子地图质量提升自动化方法。
背景技术:
GIS地理空间数据是空间信息服务平台的基石,其数据质量直接影响用户对空间信息服务平台的使用感受。每年更新数据与原始数据间存在一定误差偏移,导致已有底图数据与更新电子地图数据存在不匹配,直接影响行业生产和业务应用。因此提升数据质量是GIS地理空间数据的基础业务。
现有的矢量电子地图质量提升技术路线,基于空间精度提升的影像进行矢量数据分区域的空间精度提升。当矢量数据与遥感影像数据不完全套和时,纠偏流程:首先在已有影像和精度提升后的影像上手工选取同名点,生成已有矢量数据与精度提升后影像数据同名点的误差连线;然后根据误差连线计算误差距离和误差方向,基于误差距离和误差方向进行空间聚类形成误差均质区域;对不同区域,根据误差连线进行几何纠正。
传统的矢量影像套和方法在ARCGIS软件中处理,人工操作重复导入影像,凭视觉经验手工选取同名点,然后利用空间纠正功能对矢量数据进行纠偏。显然现有方法存在着自动化程度低,人工干预多,人为干扰因素大等缺陷,很难提高生产效率。
技术实现要素:
本发明的目的在于提供一种基于遥感的矢量电子地图质量提升自动化方法,从而解决现有技术中存在的前述问题。
为了实现上述目的,本发明所述基于遥感的矢量电子地图质量提升自动化方法,所述方法包括:
S1,获取原始影像和精度提升后影像自动匹配同名点;
S2,在匹配后的同名点的基础上建立模型,设置阈值进行错误点剔除;
S3,变换同名点坐标系,并将任意一个区域的多个同名点文件合并为一个文件A;
S4,用文件A对矢量电子地图进行纠正,得到质量提升后的矢量电子地图,并与精度提升后的影像套和,完成基于遥感的矢量电子地图质量提升。
优选地,步骤S1,采用基于区域灰度匹配方法自动产生匹配点,具体为:比较原始影像和精度提升后影像相同局部区域的灰度值,判断灰度值相似程度是否符合预先设置的阈值,如果是,则得到匹配点,继续比较下一个局部区域;如果否,则为找到匹配点,继续比较下一个局部区域,最终得到多对同名点对。
优选地,步骤S1,具体按照下述步骤实现:
S11,将原始影像和精度提升后影像一个作为左影像,另一个作为右影像;
S12,在所述左影像中选定任意一个待定点W作为目标点,以目标点为终点,选取灰度阵列作为目标区域;
S13,在右影像上建立一个大于所述目标区域的灰度阵列作为搜索区;然后依次在搜索区中选取与所述目标区域等大的灰度阵列作为搜索窗口,计算所述搜索窗口与所述目标区域的相似性测度,当相似性测度A为所有相似性测度中的最大值时,与所述相似性测度A对应的搜索窗口A的中心像素作为待定点W的同名点。
更优选地,相似性测度的度量方法采为相关系数法。
优选地,步骤S1中,进行自动匹配同名点涉及到的参数包括:匹配方法、同名点个数的阈值、最小相关系数、移动窗口大小和搜索窗口大小。
优选地,步骤S2,具体按照下述步骤实现:
S21,建立原始影像和精度提升后影像之间的变换模型,在同名点的基础上,采用最小二乘多项式,使用矩阵反演计算模型系数;
具体流程:将原始影像和精度提升后影像一个作为左影像,另一个作为右影像,左、右两幅影像的同名点对分别表示为(x,y)、(x’,y’),根据最小二乘多项式为式(1),使用矩阵反演计算模型系数:
a1,b1,c1,a2,b2,c2均为变换模型的变形系数;
S22,由变换模型计算预测点位置,即根据计算出的变换模型的变形系数a1,b1,c1,a2,b2,c2,并结合公式(2),将左方影像同名点坐标(x,y)变换至右方影像,即为预测点位置(x0,y0);
S23,采用欧氏距离公式(3)计算预测点(x0,y0)和实际匹配同名点(x’,y’)的误差距离dij;
S24,设定误差距离的阈值,当某对同名点对计算得到的误差距离大于预先阈值时,去除该对同名点;剩下的同名点集合即为所求的同名点集。
优选地,步骤S3,具体按照下述步骤实现:
S31,同名点坐标系转换,将同名点的像素平面坐标转化为同名点的经纬度的大地坐标系;
S32,多个同名点文件拼接为一个全区域的同名点文件,具体为:编程实现打开多个同名点文件,逐一读取文件内容,并将读取内容写入到一个新建的文件中。
优选地,步骤S4,用文件A用空间校正中的橡皮拉伸方法配准矢量数据坐标和调整形状采对矢量电子地图进行纠正,使纠正后的矢量电子地图与精度提升的影像套和一致。
优选地,所述矢量电子地图比例尺不设限。
本发明的有益效果是:
本发明提出一种针对矢量电子地图和影像套和的电子地图质量提升的自动化方法,整套生产工艺流程采用编程实现,克服现有方法的人工操作,效率低下的问题,无需任何附加信息和人工干预,自动化程度高,且套和精度满足效果,其结果直接面向应用。
附图说明
图1是基于遥感的矢量电子地图质量提升自动化方法的流程示意图一;
图2是基于遥感的矢量电子地图质量提升自动化方法的流程示意图二。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施方式仅仅用以解释本发明,并不用于限定本发明。
本实施例所述基于遥感的矢量电子地图质量提升自动化方法,所述方法包括:
S1,原始影像和精度提升后影像自动匹配同名点;
采用基于区域灰度匹配方法自动产生匹配点,具体为:比较原始影像和精度提升后影像相同局部区域的灰度值,判断灰度值相似程度是否符合预先设置的阈值,如果是,则得到匹配点,继续比较下一个局部区域;如果否,则为找到匹配点,继续比较下一个局部区域,最终得到多对同名点对。
S2,在匹配后的同名点的基础上建立模型,设置阈值进行错误点剔除;
S3,变换同名点坐标系,并将任意一个区域的多个同名点文件合并为一个文件A;
S4,用文件A对矢量电子地图进行纠正,得到质量提升后的矢量电子地图,并与精度提升后的影像套和,完成基于遥感的矢量电子地图质量提升。
更详细的解释说明为:
(一)步骤S1,具体按照下述步骤实现:
S11,将原始影像和精度提升后影像一个作为左影像,另一个作为右影像;
S12,在所述左影像中选定任意一个待定点W作为目标点,以目标点为终点,选取灰度阵列作为目标区域;
S13,在右影像上建立一个大于所述目标区域的灰度阵列作为搜索区;然后依次在搜索区中选取与所述目标区域等大的灰度阵列作为搜索窗口,计算所述搜索窗口与所述目标区域的相似性测度,当相似性测度A为所有相似性测度中的最大值时,与所述相似性测度A对应的搜索窗口A的中心像素作为待定点W的同名点。相似性测度的度量方法采为相关系数法,即:
其中
在公式(1)的基础上导出了归一化标准相关系数公式(2):
式(2)中的和分别表示矢量x和y的均值。
步骤S1中,进行自动匹配同名点涉及到的参数包括:匹配方法、同名点个数的阈值、最小相关系数、移动窗口大小和搜素窗口大小。
(二)步骤S2,具体按照下述步骤实现:
S21,建立原始影像和精度提升后影像之间的变换模型,在同名点的基础上,采用最小二乘多项式,使用矩阵反演计算模型系数;
具体流程:将原始影像和精度提升后影像一个作为左影像,另一个作为右影像,左、右两幅影像的同名点对分别表示为(x,y)、(x’,y’),根据最小二乘多项式为式(3),使用矩阵反演计算模型系数:
a1,b1,c1,a2,b2,c2均为变换模型的变形系数;
由于同名点对数据量较大且存在一定误差,引入拟合的思想,使得在所有数据点处的残差都较小,即所有误差组合最小,即利用上述原则确定拟合多项式的方法为最小二乘法多项式拟合,确定上述多项式的过程也就是确定多项式中的系数的过程。
S22,由变换模型计算预测点位置,即根据计算出的变换模型的变形系数a1,b1,c1,a2,b2,c2,并结合公式(4),将左方影像同名点坐标(x,y)变换至右方影像,即为预测点位置(x0,y0);
S23,采用欧氏距离公式(5)计算预测点(x0,y0)和实际匹配同名点(x’,y’)的误差距离dij;
S24,设定误差距离的阈值,当某对同名点对计算得到的误差距离大于预先阈值时,去除该对同名点;剩下的同名点集合即为所求的同名点集。
(三)步骤S3,具体按照下述步骤实现:
S31,同名点坐标系转换,将同名点的像素平面坐标转化为同名点的经纬度的大地坐标系;
S32,多个同名点文件拼接为一个全区域的同名点文件,具体为:编程实现打开多个同名点文件,逐一读取文件内容,并将读取内容写入到一个新建的文件中。
步骤S3中保存每对影像同名点,所有影像匹配点生产一个文件;由于整个省区或更大区域的影像分幅分别对应匹配生成同名点,所以需要保证整个区域为一个同名点文件。
(四)步骤S4,用文件A用空间校正中的橡皮拉伸方法配准矢量数据坐标和调整形状采对矢量电子地图进行纠正,使纠正后的矢量电子地图与精度提升的影像套和一致。所述矢量电子地图比例尺不设限。
(五)纠正后矢量与精度提升后影像套和一致性检查;
矢量影像套和一致性检查,以遥感影像作为对照底图,人工目视遍历图层检查。由于整个省级区域偶素矢量数据覆盖范围广,矢量要素复杂繁多,为了准确评估省整体矢量数据质量,检查采取随机抽样的方法。
矢量图层空间套和精度检查主要内容包括:1.图层要素相对遥感影像的整体偏移;2.图层要素相对遥感影像的不规则偏移。检查结果是矢量电子地图与遥感影像基本套和。
通过采用本发明公开的上述技术方案,得到了如下有益的效果:
本发明提出一种针对矢量电子地图和影像套和的电子地图质量提升的自动化方法,整套生产工艺流程采用编程实现,克服现有方法的人工操作,效率低下的问题,无需任何附加信息和人工干预,自动化程度高,且套和精度满足效果,其结果直接面向应用。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视本发明的保护范围。