一种基于反距离权重插值反函数的等值线图像数据提取方法技术

本发明专利技术涉及图像数据信息的提取方法，具体说是一种基于反距离权重插值反函数的等值线图像数据提取方法。具体在等值线图矢量化基础上，由等值线节点值构成已知点矩阵，由待提取数据点值构成未知点矩阵，及已知点和未知点空间关系构成的反距离权重矩阵，共同构建反距离权重插值的反函数，再根据反距离权重矩阵的秩，判断并求解未知点矩阵的实数解。距离权重插值原理基础，可计算等值线图像范围内任意点位的数值信息，且准确性高、可批量计算，提取的数值信息可为进一步的多源数据融合提供数据基础。

【技术实现步骤摘要】
一种基于反距离权重插值反函数的等值线图像数据提取方法
本专利技术涉及图像数据信息的提取方法，具体说是一种基于反距离权重插值反函数的等值线图像数据提取方法。
技术介绍
随着大数据时代的来临，2015年国务院印发《促进大数据发展行动纲要》，将发展大数据上升为国家战略，落实“互联网+”，建设数据强国。在环境领域，大数据除了依靠传感器高频接收与互联网实时产生的数据外，不同尺度、定点监测、长周期调查也是数据的重要来源，而图像是环境领域历史信息最为丰富的数据源。现有数据提取方法多为根据图像像素颜色等特征进行解译，但提取结果仅为分类数据。具体的说，对于等值线图而言，提取的分类数据对应着两条等值线之间的数值范围，而非可参与运算的数值。该类方法提取的数据只能用于分类信息的归类比较，无法用于数值数据的多源数据融合，更无法满足大数据精确分析与深度挖掘的需求。等值线图是环境领域中除遥感影像数据外最为直观的图像。等值线图是利用已知的点位数值，采用插值方法推算区域内各点数值，然后将数值相等点联成曲线。插值方法主要是由反距离权重插值、克里格插值、泰森多边形插值等，由于方法适用范围及环境领域原始信息不足等问题，反距离权重插值是环境领域最常用的插值方法。因此，如何根据反距离权重插值方法原理，以及等值线上已知点数值信息，提取出等值线图中任意位置的数据，对环境领域图像的数据提取及多源数据融合具有重要的理论和技术指导意义。
技术实现思路
本专利技术目的在于提供一种基于反距离权重插值反函数的等值线图像数据提取方法。为实现上述目的本专利技术采用的技术方案为：一种基于反距离权重插值反函数的等值线图像数据提取方法，包括以下步骤：1)在等值线图F矢量化基础上，提取等值线节点t(xt,yt,s)信息，即大地经度、大地纬度与等值线值，由等值线节点s值构成已知点矩阵T；2)在等值线图F范围内，设定待提取数据点p(xp,yp,v)，其大地经度和大地纬度为已知，待提取数值v为未知，由待提取数据点v值构成未知点矩阵P；3)利用已知点t和未知点p的大地经纬度，构成反距离权重矩阵W，与已知点矩阵T和未知点矩阵P，构建反距离权重插值反函数P＝W-1T；4)计算反距离权重矩阵W的秩，若W为满秩，则根据反距离权重插值反函数，求未知点矩阵P的唯一实数解；若W不为满秩，则调整已知点数量或位置，直至W的秩为满秩为止。所述已知点t与未知点p数量分别为m和n，需满足m≥n。所述反距离权重矩阵W，是根据未知点p与已知点t的距离关系，计算得到的权重值，具体表达为：其中，wij是第i个未知点与第j个已知点之间的权重值，dij是第i个未知点与第j个已知点之间的距离，dik是第i个未知点与第k个已知点之间的距离，h为幂数。所述反距离权重插值反函数，具体是反距离权重插值所构建多元一次方程组W·P＝T的逆运算，即P＝W-1T，具体表达为：其中，W＝[wij]m×n是权重矩阵，T＝[sij]m×1是已知点矩阵，P＝[vij]n×1是未知点矩阵。本专利技术具有以下优点及有益效果：本专利技术在等值线图矢量化基础上，由等值线节点值构成已知点矩阵，由待提取数据点值构成未知点矩阵，及已知点和未知点空间关系构成的反距离权重矩阵，共同构建反距离权重插值的反函数，再根据反距离权重矩阵的秩，判断并求解未知点矩阵的实数解。距离权重插值原理基础，可计算等值线图像范围内任意点位的数值信息，且准确性高、可批量计算，提取的数值信息可为进一步的多源数据融合提供数据基础。本专利技术以反距离权重插值原理基础，利用等值线上点位数值已知的特点，采用反距离权重插值的反函数，计算等值线图像范围内任意点位的数值信息，且准确性高、可批量计算，提取的数值信息可为进一步的多源数据融合提供数据基础。附图说明图1是实施例1生成等值线图的测试点位分布。图2是实施例1反距离权重插值生成的等值线图。图3是实施例1等值线图中选择参与运算的节点。图4是实施例1提取数据值与测试值相关性。图5是实施例1已知点数量对提取数据值准确性的影响。图6是实施例2反距离权重插值生成的等值线图。图7是实施例2提取数据值与测试值相关性。具体实施方式下面结合附图及实施例对本专利技术做进一步的详细说明。一种基于反距离权重插值反函数的等值线图像数据提取方法。步骤是：1)在等值线图F矢量化基础上，提取m个等值线节点t信息，由等值线节点值构成已知点矩阵T；2)在等值线图F范围内，设定n个待提取数据点p，则由待提取数据点值构成未知点矩阵P；3)根据已知点和未知点的大地经纬度构成的反距离权重矩阵W，与已知点矩阵T和未知点矩阵P，构建反距离权重插值的反函数，即P＝W-1T；4)根据反距离权重矩阵W的秩，求解未知点矩阵P的实数解，或者调整已知点数量或位置，重新建立方程求解。步骤1：构建已知点矩阵T1)利用ArcGIS软件，根据等值线图F复杂程度，通过手动或ArcScan工具对其矢量化，设置参与计算等值线节点t共m个，m＞0；2)利用ArcGIS软件的AddXYCoordinates功能，提取等值线节点t的大地经度xt和大地纬度yt，并根据对应等值线赋予节点t数值s，三者组成完成的节点信息t(xt,yt,s)；3)由等值线节点s值构成了已知点矩阵T，即T＝[sij]m×1。步骤2：构建未知点矩阵P1)利用ArcGIS软件，在等值线图F范围内，手动或自动设置待提取数据点，设置参与计算的待提取数据点等值线节点p共n个，n＞0且n≤m；2)与已知点操作相同，利用ArcGIS软件的AddXYCoordinates功能，提取等值线节点p的大地经度xp和大地纬度yp，而待提取的未知数值设为v，三者组成完成的待提取点位信息p(xp,yp,v)；3)由待提取数据点v值构成了未知点矩阵P，即P＝[vij]n×1。步骤3：构建反距离权重插值反函数1)反距离权重矩阵W，是根据未知点p与已知点t的距离关系，计算得到的权重值，具体表达为：其中，wij是第i个未知点与第j个已知点之间的权重值，dij是第i个未知点与第j个已知点之间的距离，dik是第i个未知点与第k个已知点之间的距离，h为幂数。2)构成的反距离权重矩阵W为[wij]m×n，其中i表示第i个已知点，i＝1,2,3……m，j表示第j个未知点，j＝1,2,3……n。3)根据反距离权重插值原理，将反距离权重矩阵W、已知点矩阵T和未知点矩阵P，构建多元一次方程组W·P＝T，具体表达式为：其中，W＝[wij]m×n是权重矩阵，T＝[sij]m×1是已知点矩阵，P＝[vij]n×1是未知点矩阵。4)利用矩阵运算特点，构建反距离权重插值的反函数，具体是反距离权重插值方法构建多元一次方程组W·P＝T的逆运算，即P＝W-1T，具体表达为：步骤4：求解未知点矩阵P1)计算反距离权重矩阵W的秩，记为R(W)；2)若矩阵W为满秩，即R(W)＝n，则根据反距离权重插值的反函数(公式2)，应用MatLAB自主编写程序，求未知点矩阵P的唯一实数解；3)若矩阵W不为满秩，即R(W)<n，则调整已知点数量或位置，重复步骤一到步骤三，直至R(W)＝n为止。实施例1：本实施例中等值线插值图是在25km2范围内，按照5km×5km网格布设离散点(图1)，并给每个离散点赋予0～100的随机数。利用反距离权重插值法，大本文档来自技高网...

【技术保护点】
一种基于反距离权重插值反函数的等值线图像数据提取方法，其特征在于包括以下步骤：1)在等值线图F矢量化基础上，提取等值线节点t(x

【技术特征摘要】
1.一种基于反距离权重插值反函数的等值线图像数据提取方法，其特征在于包括以下步骤：1)在等值线图F矢量化基础上，提取等值线节点t(xt,yt,s)信息，即大地经度、大地纬度与等值线值，由等值线节点s值构成已知点矩阵T；2)在等值线图F范围内，设定待提取数据点p(xp,yp,v)，其大地经度和大地纬度为已知，待提取数值v为未知，由待提取数据点v值构成未知点矩阵P；3)利用已知点t和未知点p的大地经纬度，构成反距离权重矩阵W，与已知点矩阵T和未知点矩阵P，构建反距离权重插值反函数P＝W-1T；4)计算反距离权重矩阵W的秩，若W为满秩，则根据反距离权重插值反函数，求未知点矩阵P的唯一实数解；若W不为满秩，则调整已知点数量或位置，直至W的秩为满秩为止。2.按权利要求1所述的基于反距离权重插值反函数的等值线图像...

【专利技术属性】
技术研发人员：吴波，张玲妍，郭书海，
申请(专利权)人：中国科学院沈阳应用生态研究所，
类型：发明
国别省市：辽宁,21