CN101504315A - 扩展形态学与正交子空间投影结合的端元自动提取方法 - Google Patents
扩展形态学与正交子空间投影结合的端元自动提取方法 Download PDFInfo
- Publication number
- CN101504315A CN101504315A CNA2009100784427A CN200910078442A CN101504315A CN 101504315 A CN101504315 A CN 101504315A CN A2009100784427 A CNA2009100784427 A CN A2009100784427A CN 200910078442 A CN200910078442 A CN 200910078442A CN 101504315 A CN101504315 A CN 101504315A
- Authority
- CN
- China
- Prior art keywords
- end member
- projection
- expansion
- structural element
- pixel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
一种扩展形态学与正交子空间投影结合的端元自动提取方法,步骤如下:(1)高光谱数据读入;(2)确定结构元素尺寸、初始与最大迭代次数;(3)由扩展膨胀、腐蚀操作计算结构元素邻域内混合度最低和最高的像元;(4)由步骤(3)得到的结果计算得到形态离心率指数值;(5)步骤(3)、(4)随迭代次数的增加不断重复,并利用步骤(3)扩展膨胀操作结果更新图像数据,直到达到最大迭代次数;(6)二值化MEI图像,得到端元数据集;(7)采用光谱角匹配方法计算得到第一个端元,并向得到端元的正交子空间投影更新端元数据集;(8)重复步骤(7),直到满足误差要求。该方法是一种稳定性强、可靠性高、精确度高的高光谱端元自动提取方法。
Description
技术领域
本发明涉及扩展形态学与正交子空间投影结合的端元自动提取方法,属于高光谱数据处理方法与应用技术领域,适用于高光谱数据端元提取的理论方法以及在农作物精细分类、矿物蚀变信息提取等领域的应用技术研究。
背景技术
高光谱数据处理与传统的遥感数据处理的一个本质的区别在于后者在对地物分类时通常是将单一像元根据某种规则分到某一不同的种类中,从而更加准确地反映遥感对象的状况。混合像元分解不仅是提高地物识别精度的重要方法,更是定量探测地物成分的有效途径,因此,混合像元分解是高光谱遥感数据应用的主要方向之一,该类技术的发展能够大大推动高光谱数据的应用范围与水平,并不断扩展高光谱数据的应用深度和广度。
目前端元提取的方法主要包括纯像元指数(Pure Pixel Index,PPI)、最小体积变换、迭代误差分析、顶点分量分析等,上述方法均只利用了高光谱数据提供的光谱/特征空间信息,针对这一问题,国外提出了基于数学形态学的方法,该方法同时利用了高光谱数据提供的光谱和空间信息,提高了方法的可靠性。但是该类方法存在着一定程度的问题:采用光谱角制图的方法度量光谱之间的相似度,受噪声影响较大;不同类型的端元采用了大津法,对于光谱相似的不同类型端元很难区分;并且算法的性能直接受结构元素形状与尺寸的影响。
发明内容
本发明的目的是针对现有端元提取方法从数据光谱或空间或特征信息的单一方面出发进行混合像元分解、不同类型端元难以区分等不足,提出一种扩展形态学与正交子空间投影结合的端元自动提取方法。
本发明的技术解决方案是:利用扩展形态学操作与正交子空间投影理论实现高光谱数据端元自动提取的方法,该方法主要是基于扩展数学形态学的理论,利用扩展膨胀和腐蚀操作,并通过计算形态离心率指数进行高光谱数据的端元数据集计算,然后利用光谱角匹配方法提取不同类型的端元,并通过向端元正交子空间投影消除已经提取端元的影响,实现高光谱数据所有端元的自动提取。形态偏心率指数是利用扩展膨胀和腐蚀的结果计算数据的正交投影散度,从而将图像中端元提取出来,引入正交子空间投影方法实现不同类型端元的有效区分,并且在本发明中为了避免由于形态学中结构元素的大小对方法性能的限制和降低计算过程中引入的人为误差的影响,采用迭代的方法,利用扩展膨胀操作不断更新参与计算的图像数据来实现。
本发明扩展形态学与正交子空间投影结合的端元自动提取方法,其步骤如下:
(1)高光谱遥感数据读入;
(2)根据高光谱数据特点,确定结构元素的尺寸、初始迭代次数以及最大迭代次数;
(3)由扩展数学形态学膨胀和腐蚀操作计算每一个结构元素邻域内混合度最低和最高的像元;
(4)由步骤(3)得到的混合度最低和最高的像元结果计算得到形态离心率指数值;
(5)步骤(3)、(4)随迭代次数的增加不断重复,并利用步骤(3)得到的混合度最低的结果不断更新图像数据,直到达到最大的迭代次数;
(6)二值化形态离心率指数MEI图像,得到端元数据集;
(7)利用由步骤(6)得到的端元数据集中任一元素作为参考光谱,采用光谱角匹配方法计算得到参考光谱与端元数据集中所有元素光谱角最大的元素,作为第一个端元,并向第一个端元的正交子空间投影,更新端元数据集;
(8)不断重复步骤(7),直到满足误差要求,从而得到高光谱遥感图像中所有的端元。
其中,步骤(2)中所述的结构元素尺寸为5×5,初始迭代次数t=1,最大迭代次数tmax=5;
其中,步骤(3)所述的“由扩展数学形态学膨胀和腐蚀操作计算每一个结构元素邻域内混合度最低和最高的像元”,其扩展数学形态学膨胀和腐蚀操作定义如下:
其中,d(x,y)、e(x,y)分别表示像元(x,y)在结构元素b(x,y)内的扩展膨胀、腐蚀操作结果;f(x,y)是图像像元(x,y)的值,b(x,y)是结构元素;arg_Max,arg_Mn分别表示使得D达到最大和最小的像元向量;D是为了确定像元混合度的多维向量的排序关系,引入的一个多维向量的度量算子。该度量算子由结构元素内各个像元累加距离计算得到,定义如下:
其中,(s,t)∈Db表示(s,t)是属于结构元素b(x,y)内的点,dist是测量N维向量的逐点线性距离,采用正交投影散度(Orthogonal Projection Divergence,OPD)计算该距离,考虑两个N维光谱信号si=[si1,si2,...,siN]T,sj=[sj1,sj2,...,sjN]T,则N维光谱信号si和sj之间的正交投影散度OPD表示为:
因此,累加距离D能够根据像元的混合程度大小排序结构元素中的向量。通过以上的分析表明,扩展到高光谱图像的膨胀结果得到的是在结构元素内纯度最高的像元,腐蚀结果得到的是在结构元素内混合度最高的像元。
其中,步骤(4)中所述的“由步骤(3)得到的混合度最低和最高的像元结果计算得到形态离心率指数值”,其含义说明如下:通过正交投影散度融合结构元素邻域内膨胀和腐蚀结果,计算得到形态离心率指数MEI,计算公式如下:
MEI(x,y)=OPD[d(x,y),e(x,y)]
其中,步骤(5)中所述的“步骤(3)、(4)随迭代次数的增加不断重复,并利用步骤(3)得到的混合度最低的结果不断更新图像数据,直到达到最大的迭代次数”,其含义说明如下:设ft(x,y)为第t次迭代计算得到的像元点(x,y)的图像数据,ft-1(x,y)为第i-1次迭代计算得到的(x,y)的图像数据,更新图像数据的目的是避免人为引入误差或结构元素尺寸对提取端元纯度的影响,因此采用扩展膨胀操作实现图像数据更新,即
其中,步骤(6)中所述的“计算迭代后的形态离心率指数MEI图像,并进行二值化,得到端元数据集”,其含义说明如下:利用[m-3σ,m+3σ]进行MEI图像的二值化,得到所有端元的空间位置,其中,m为MEI图像均值,σ为MEI图像的标准差。
其中,步骤(7)中所述的“利用由步骤(6)得到的端元数据集中任一元素作为参考光谱,采用光谱角匹配方法计算得到参考光谱与端元数据集中所有元素光谱角最大的元素,作为第一个端元,并向第一个端元的正交子空间投影,更新端元数据集”,其含义说明如下:基于线性混合模型,利用正交子空间投影(Orthogonal Subspace Projection,OSP)概念进行混合像元分解时,仍可以继续将端元矩阵分解,则混合模型可以表示为:
r=dαd+UαU+n
其中r为端元数据集,d为得到的端元光谱信号,αd为端元光谱信号的混合百分比,U为剩余端元光谱信号,αU为剩余端元光谱信号的混合百分比,n为噪声信号;通过MEI图像中提取的端元数据,利用光谱角匹配方法得到其中一个端元光谱信号d=e1,然后利用向其正交子空间投影得到来消除端元数据中的信号d=e1,其中 将代入基于OSP的混合模型中得到 可以看出,不仅消除了d=e1的影响而且有效的压制了噪声然后利用更新端元数据集。
其中,步骤(8)中所述的“不断重复步骤(7),直到满足误差要求,从而得到图像中所有端元”,其含义说明如下:经过m-1次投影之后得到d={e1,e2,...,em},投影循环结束的条件由分解误差决定,从而得到了m个端元。
本发明与现有技术相比的优点在于:本方法同时利用了扩展形态学与正交子空间投影的理论,从数学机理上有效的克服了现有端元提取方法无法将不同类型端元有效区分的问题,实现了无先验信息情况下高光谱数据端元自动提取。它具有以下的优点:(1)通过扩展形态学膨胀和腐蚀操作实现图像数据中纯度最高和混合度最高的像元提取,并利用正交投影散度融合扩展膨胀和腐蚀结果得到形态离心率指数,实现了高光谱数据端元提取,有效地抑制了背景等因素对端元提取结果的影响,提高了方法的稳定性;(2)采用迭代的方法,利用扩展膨胀操作不断更新图像数据来克服由于形态学中结构元素的尺寸等对方法性能的限制;(3)引入正交子空间投影和光谱角匹配,实现不同类型端元的有效区分,提高了方法的精确性。
附图说明
图1为本发明扩展形态学与正交子空间投影结合的端元自动提取方法的流程图。
具体实施方式
为了更好的说明本发明涉及的扩展形态学与正交子空间投影结合的端元自动提取方法,利用星载Hyperion高光谱成像仪数据,采用本发明方法进行西藏驱龙地区蚀变矿物端元提取。本发明一种扩展形态学与正交子空间投影结合的端元自动提取方法,具体实现步骤如下:
(1)高光谱遥感数据读入:读入Hyperion高光谱数据,数据维数为255×550×196;
(2)根据高光谱数据特点以及应用目的,确定结构元素的尺寸和初始迭代次数:结构元素尺寸为5×5,初始迭代次数t=1,最大迭代次数tmax=5;
(3)由扩展数学形态学膨胀和腐蚀操作计算每一个结构元素邻域内混合度最低和最高的像元;
为了实现更加可靠、稳定、精确的高光谱数据端元自动提取,本发明中的方法同时利用了光谱和空间信息,由此引入了扩展数学形态学膨胀和腐蚀操作,定义如下:
其中,d(x,y)、e(x,y)分别表示像元(x,y)在结构元素b(x,y)内的扩展膨胀、腐蚀操作结果;f(x,y)是图像像元(x,y)的反射率值,b(x,y)是结构元素;arg_Max,arg_Min分别表示使得D达到最大和最小的像元向量;D是为了确定像元混合度的多维向量的排序关系,引入的一个高维向量的度量算子。该度量算子由结构元素内各个像元累加距离计算得到,定义如下:
其中,(s,t)∈Db表示(s,t)是属于结构元素b(x,y)内的点,dist是测量N维向量的逐点线性距离。为了有效地利用高光谱数据提供的光谱和空间信息,本发明采用正交投影散度OPD计算该距离,考虑两个196维光谱信号si=[si1,si2,...,si196]T,sj=[sj1,sj2,...,sj196]T,则196维光谱信号si和sj之间的正交投影散度OPD表示为:
因此,累加距离D能够根据像元的混合程度大小排序结构元素中的向量。通过以上的分析表明,扩展到高光谱图像的膨胀结果得到的是在结构元素内纯度最高的像元,腐蚀结果得到的是在结构元素内混合度最高的像元。
(4)由步骤(3)得到的混合度最低和最高的像元结果计算得到形态离心率指数值;
为了克服利用单一膨胀或腐蚀操作人为引入的噪声和背景影响,以及为了定量化图像数据像元的纯度与混合度,采用了形态离心率指数。通过正交投影散度融合结构元素邻域内膨胀和腐蚀结果,计算得到形态离心率指数MEI,实现高光谱数据的端元提取。
MEI(x,y)=OPD[d(x,y),e(x,y)]
(5)步骤(3)、(4)随迭代次数的增加不断重复,并利用步骤(3)得到的混合度最低的结果不断更新图像数据,直到达到最大的迭代次数;
为了提取图像数据中纯度最高的像元,利用能够得到结构元素内纯度最高像元的扩展膨胀操作不断更新图像数据,从而有效地提高端元提取精确性。
设ft(x,y)为第t次迭代计算得到的像元点(x,y)的参与计算图像数据,ft-1(x,y)为第i-1次迭代计算得到的(x,y)的参与计算图像数据,更新图像数据的准则是避免人为引入误差或结构元素尺寸对提取端元纯度的影响,因此采用扩展膨胀操作实现图像数据更新,即
(6)计算迭代后的形态离心率指数MEI图像,并进行二值化,得到端元数据集;
为了自动将形态离心率指数MEI图像二值化,实现所有端元数据的提取,利用[m-3σ,m+3σ]进行MEI图像的二值化,得到所有端元的空间位置,其中,m为MEI图像均值,σ为MEI图像的标准差。
(7)利用由步骤(6)得到的端元数据集中任一元素作为参考光谱,采用光谱角匹配方法计算得到参考光谱与端元数据集中所有元素光谱角最大的元素,作为第一个端元,并向第一个端元的正交子空间投影,更新端元数据集;
为了克服现有基于扩展数学形态学方法无法准确区分不同类型端元的局限,本发明方法提出了基于正交子空间投影的方法,实现了不同类型端元的有效提取。
基于线性混合模型,利用正交子空间投影OSP概念进行混合像元分解时,仍可以继续将端元矩阵分解,则混合模型可以表示为:
r=dαd+UαU+n
其中,r为端元数据集,d为得到的端元光谱信号,αd为端元光谱信号的混合百分比,U为剩余端元光谱信号,αU为剩余端元光谱信号的混合百分比,n为噪声信号;通过MEI图像中提取的端元数据,利用光谱角匹配方法得到其中一个端元光谱信号d=e1,然后利用向其正交子空间投影得到来消除端元数据中的信号d=e1,其中 将代入基于OSP的混合模型中得到 可以看出,不仅消除了d=e1的影响而且有效的压制了噪声然后利用更新端元数据集。
(8)不断重复步骤(7),直到满足误差要求,从而得到图像中所有端元。
经过m-1次投影之后得到d={e1,e2,...,em},投影循环结束的条件由分解误差决定,从而得到了m个端元。分解误差为当前提取的端元与前一次提取的端元光谱变化值,本发明中的循环终止条件为:‖ei-ei-1‖2≤0.01。
以上整个实现过程不需要人为干预或任何先验信息,因此,大大提高了端元提取方法的自动性。通过本发明方法,得到了西藏驱龙地区的五种主要的蚀变矿物端元光谱:贫铝绢云母,富铝绢云母,高岭石化,绿泥石化,褐铁矿化。
Claims (6)
1、一种扩展形态学与正交子空间投影结合的端元自动提取方法,其特征在于:它包含以下步骤:
(1)高光谱遥感数据读入;
(2)根据高光谱数据特点,确定结构元素的尺寸、初始迭代次数以及最大迭代次数;
(3)由扩展数学形态学膨胀和腐蚀操作计算每一个结构元素邻域内混合度最低和最高的像元;
(4)由步骤(3)得到的混合度最低和最高的像元结果计算得到形态离心率指数值;
(5)步骤(3)、(4)随迭代次数的增加不断重复,并利用步骤(3)得到的混合度最低的结果不断更新图像数据,直到达到最大的迭代次数;
(6)二值化形态离心率指数MEI图像,得到端元数据集;
(7)利用由步骤(6)得到的端元数据集中任一元素作为参考光谱,采用光谱角匹配方法计算得到参考光谱与端元数据集中所有元素光谱角最大的元素,作为第一个端元,并向第一个端元的正交子空间投影,更新端元数据集;
(8)不断重复步骤(7),直到满足误差要求,从而得到高光谱遥感图像中所有的端元。
2、根据权利要求1所述的扩展形态学与正交子空间投影结合的端元自动提取方法,其特征在于:步骤(3)中所述的由扩展数学形态学膨胀和腐蚀操作计算每一个结构元素邻域内混合度最低和最高的像元,计算公式如下:
其中,d(x,y)、e(x,y)分别表示像元(x,y)在结构元素b(x,y)内的扩展膨胀、腐蚀操作结果;f(x,y)是图像像元(x,y)的值,b(x,y)是结构元素;arg_Max,arg_Min分别表示使得D达到最大和最小的像元向量;D是一个高维向量的度量算子,计算公式如下:
其中,(s,t)∈Db表示(s,t)是属于结构元素b(x,y)内的点,采用正交投影散度计算两个N维光谱信号Si=[Si1,Si2,...,SiN]T和Sj,=[Sj1,Sj2,...,SjN]T之间的dist距离,正交投影散度OPD表示为:
其中, k=i,j,并且IN×N是N*N维的单位矩阵。
3、根据权利要求1所述的扩展形态学与正交子空间投影结合的端元自动提取方法,其特征在于:其步骤(4)中所述的计算形态离心率指数MEI值计算公式如下:
MEI(x,y)=OPD[d(x,y),e(x,y)]
其中,OPD表示正交投影散度,d(x,y)表示像元(x,y)在结构元素b(x,y)内的扩展膨胀操作结果,e(x,y)表示像元(x,y)在结构元素b(x,y)内的扩展腐蚀操作结果。
4、根据权利要求1所述的扩展形态学与正交子空间投影结合的端元自动提取方法,其特征在于:步骤(5)中所述的更新图像数据采用扩展膨胀操作实现:
其中,ft(x,y)为第t次迭代计算得到的像元点(x,y)的图像数据,ft-1(x,y)为第i-1次迭代计算得到的(x,y)的图像数据,b(x,y)为结构元素。
5、根据权利要求1所述的扩展形态学与正交子空间投影结合的端元自动提取方法,其特征在于:步骤(6)中所述的二值化形态离心率指数MEI图像方法如下:利用[m-3σ,m+3σ]进行MEI图像的二值化,得到所有端元的空间位置,其中,m为MEI图像均值,σ为MEI图像的标准差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100784427A CN101504315B (zh) | 2009-02-23 | 2009-02-23 | 扩展形态学与正交子空间投影结合的端元自动提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100784427A CN101504315B (zh) | 2009-02-23 | 2009-02-23 | 扩展形态学与正交子空间投影结合的端元自动提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101504315A true CN101504315A (zh) | 2009-08-12 |
CN101504315B CN101504315B (zh) | 2011-01-05 |
Family
ID=40976622
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100784427A Expired - Fee Related CN101504315B (zh) | 2009-02-23 | 2009-02-23 | 扩展形态学与正交子空间投影结合的端元自动提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101504315B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101806898A (zh) * | 2010-03-19 | 2010-08-18 | 武汉大学 | 基于端元可变的高光谱遥感影像目标探测方法 |
CN101950361A (zh) * | 2010-09-06 | 2011-01-19 | 中国科学院遥感应用研究所 | 一种基于光谱匹配度的遥感影像专题信息自适应提取方法 |
CN102663402A (zh) * | 2012-04-21 | 2012-09-12 | 西北工业大学 | 基于修正扩展形态学算子的高光谱遥感图像端元提取方法 |
CN103926203A (zh) * | 2014-04-29 | 2014-07-16 | 中国科学院遥感与数字地球研究所 | 一种针对地物光谱不确定性的光谱角度制图方法 |
CN102117483B (zh) * | 2009-12-31 | 2014-11-05 | 核工业北京地质研究院 | 不同空间分辨率的多光谱遥感图像融合方法 |
CN105068136A (zh) * | 2015-07-27 | 2015-11-18 | 中国地质调查局武汉地质调查中心 | 基于多源信息的中南半岛示范区铜金矿潜力定位评价方法 |
CN105160642A (zh) * | 2015-08-29 | 2015-12-16 | 西安科技大学 | 高光谱成像的潜指纹图像增强方法 |
CN105427319A (zh) * | 2015-11-27 | 2016-03-23 | 杭州电子科技大学 | 一种基于自动形态学端元提取的高光谱图像数据解混方法 |
CN107133954A (zh) * | 2017-04-10 | 2017-09-05 | 北京环境特性研究所 | 基于数学形态学的光谱角度匹配方法 |
CN108021874A (zh) * | 2017-11-23 | 2018-05-11 | 北京环境特性研究所 | 一种基于空-谱结合的高光谱端元提取预处理方法 |
CN108073895A (zh) * | 2017-11-22 | 2018-05-25 | 杭州电子科技大学 | 一种基于解混预处理的高光谱目标检测方法 |
CN109840544A (zh) * | 2018-12-19 | 2019-06-04 | 嘉兴市恒创电力设备有限公司 | 一种高光谱图像多端元光谱混合分析方法及装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6724940B1 (en) * | 2000-11-24 | 2004-04-20 | Canadian Space Agency | System and method for encoding multidimensional data using hierarchical self-organizing cluster vector quantization |
CN100483147C (zh) * | 2007-11-02 | 2009-04-29 | 北京航空航天大学 | 高光谱亚像元目标探测方法及装置 |
CN100590402C (zh) * | 2008-01-11 | 2010-02-17 | 北京航空航天大学 | 一种光谱和空间信息结合的高光谱数据分类方法 |
-
2009
- 2009-02-23 CN CN2009100784427A patent/CN101504315B/zh not_active Expired - Fee Related
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102117483B (zh) * | 2009-12-31 | 2014-11-05 | 核工业北京地质研究院 | 不同空间分辨率的多光谱遥感图像融合方法 |
CN101806898A (zh) * | 2010-03-19 | 2010-08-18 | 武汉大学 | 基于端元可变的高光谱遥感影像目标探测方法 |
CN101950361A (zh) * | 2010-09-06 | 2011-01-19 | 中国科学院遥感应用研究所 | 一种基于光谱匹配度的遥感影像专题信息自适应提取方法 |
CN102663402A (zh) * | 2012-04-21 | 2012-09-12 | 西北工业大学 | 基于修正扩展形态学算子的高光谱遥感图像端元提取方法 |
CN102663402B (zh) * | 2012-04-21 | 2013-06-26 | 西北工业大学 | 基于修正扩展形态学算子的高光谱遥感图像端元提取方法 |
CN103926203A (zh) * | 2014-04-29 | 2014-07-16 | 中国科学院遥感与数字地球研究所 | 一种针对地物光谱不确定性的光谱角度制图方法 |
CN103926203B (zh) * | 2014-04-29 | 2016-04-20 | 中国科学院遥感与数字地球研究所 | 一种针对地物光谱不确定性的光谱角度制图方法 |
CN105068136A (zh) * | 2015-07-27 | 2015-11-18 | 中国地质调查局武汉地质调查中心 | 基于多源信息的中南半岛示范区铜金矿潜力定位评价方法 |
CN105160642B (zh) * | 2015-08-29 | 2017-10-03 | 西安科技大学 | 高光谱成像的潜指纹图像增强方法 |
CN105160642A (zh) * | 2015-08-29 | 2015-12-16 | 西安科技大学 | 高光谱成像的潜指纹图像增强方法 |
CN105427319A (zh) * | 2015-11-27 | 2016-03-23 | 杭州电子科技大学 | 一种基于自动形态学端元提取的高光谱图像数据解混方法 |
CN105427319B (zh) * | 2015-11-27 | 2018-10-02 | 杭州电子科技大学 | 一种基于自动形态学端元提取的高光谱图像数据解混方法 |
CN107133954A (zh) * | 2017-04-10 | 2017-09-05 | 北京环境特性研究所 | 基于数学形态学的光谱角度匹配方法 |
CN108073895A (zh) * | 2017-11-22 | 2018-05-25 | 杭州电子科技大学 | 一种基于解混预处理的高光谱目标检测方法 |
CN108073895B (zh) * | 2017-11-22 | 2021-12-14 | 杭州电子科技大学 | 一种基于解混预处理的高光谱目标检测方法 |
CN108021874A (zh) * | 2017-11-23 | 2018-05-11 | 北京环境特性研究所 | 一种基于空-谱结合的高光谱端元提取预处理方法 |
CN109840544A (zh) * | 2018-12-19 | 2019-06-04 | 嘉兴市恒创电力设备有限公司 | 一种高光谱图像多端元光谱混合分析方法及装置 |
CN109840544B (zh) * | 2018-12-19 | 2023-05-16 | 国网浙江省电力有限公司嘉兴供电公司 | 一种高光谱图像多端元光谱混合分析方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN101504315B (zh) | 2011-01-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101504315B (zh) | 扩展形态学与正交子空间投影结合的端元自动提取方法 | |
CN110849374B (zh) | 地下环境定位方法、装置、设备及存储介质 | |
Yang et al. | An automated method to register airborne and terrestrial laser scanning point clouds | |
CN102436652B (zh) | 一种多源遥感图像自动配准方法 | |
CN104952107A (zh) | 基于车载LiDAR点云数据的桥梁三维重建方法 | |
CN103093459B (zh) | 利用机载LiDAR点云数据辅助影像匹配的方法 | |
CN103208118B (zh) | 一种高光谱遥感影像端元提取方法 | |
Wang et al. | Window detection from mobile LiDAR data | |
CN104732553A (zh) | 一种基于多台激光辅助靶标的特征点提取方法 | |
CN107909018B (zh) | 一种稳健的多模态遥感影像匹配方法和系统 | |
CN1963390A (zh) | 一种便捷高效的三维测量方法 | |
CN101556694B (zh) | 一种旋转图像的匹配方法 | |
CN102163333B (zh) | 谱聚类的sar图像变化检测方法 | |
Jia et al. | Bidirectional interaction between BIM and construction processes using a multisource geospatial data enabled point cloud model | |
Zhu et al. | Robust registration of aerial images and LiDAR data using spatial constraints and Gabor structural features | |
CN105354841A (zh) | 一种快速遥感影像匹配方法及系统 | |
CN103578105B (zh) | 一种基于区域特征的多源异类图像配准方法 | |
Jia et al. | A cross-correction LiDAR SLAM method for high-accuracy 2D mapping of problematic scenario | |
Xu et al. | Feature curve-based registration for airborne LiDAR bathymetry point clouds | |
Parmehr et al. | Automatic registration of optical imagery with 3d lidar data using local combined mutual information | |
CN102938148A (zh) | 一种基于v-glcm的高光谱影像纹理分析方法 | |
Parmehr et al. | Automatic registration of multi-source data using mutual information | |
Cord et al. | Three-dimensional building detection and modeling using a statistical approach | |
CN104035131A (zh) | 一种地壳偏移量估计方法及装置 | |
Gan et al. | A random forest based method for urban object classification using lidar data and aerial imagery |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110105 Termination date: 20140223 |