CN111552004A - 一种遥感数据角度异常信息提取方法及系统 - Google Patents

一种遥感数据角度异常信息提取方法及系统 Download PDF

Info

Publication number
CN111552004A
CN111552004A CN202010331000.5A CN202010331000A CN111552004A CN 111552004 A CN111552004 A CN 111552004A CN 202010331000 A CN202010331000 A CN 202010331000A CN 111552004 A CN111552004 A CN 111552004A
Authority
CN
China
Prior art keywords
remote sensing
angle
sensing data
data set
data
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
Application number
CN202010331000.5A
Other languages
English (en)
Other versions
CN111552004B (zh
Inventor
姚佛军
耿新霞
杨建民
吴胜华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Mineral Resources of Chinese Academy of Geological Sciences
Original Assignee
Institute of Mineral Resources of Chinese Academy of Geological Sciences
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Mineral Resources of Chinese Academy of Geological Sciences filed Critical Institute of Mineral Resources of Chinese Academy of Geological Sciences
Priority to CN202010331000.5A priority Critical patent/CN111552004B/zh
Priority to CH00819/20A priority patent/CH717361A8/de
Publication of CN111552004A publication Critical patent/CN111552004A/zh
Application granted granted Critical
Publication of CN111552004B publication Critical patent/CN111552004B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28—Investigating the spectrum
    • G01J3/2823—Imaging spectrometer
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V8/00—Prospecting or detecting by optical means
    • G01V8/02—Prospecting
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00—Scenes; Scene-specific elements
    • G06V20/10—Terrestrial scenes
    • G06V20/13—Satellite images
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28—Investigating the spectrum
    • G01J3/2823—Imaging spectrometer
    • G01J2003/2826—Multispectral imaging, e.g. filter imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N2021/1793—Remote sensing
    • G01N2021/1797—Remote sensing in landscape, e.g. crops
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N2021/3129—Determining multicomponents by multiwavelength light
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2201/00—Features of devices classified in G01N21/00
    • G01N2201/12—Circuits of general importance; Signal processing

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Remote Sensing (AREA)
  • Astronomy & Astrophysics (AREA)
  • Chemical & Material Sciences (AREA)
  • Pathology (AREA)
  • Multimedia (AREA)
  • Theoretical Computer Science (AREA)
  • Immunology (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Geology (AREA)
  • Electromagnetism (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Testing And Monitoring For Control Systems (AREA)

Abstract

本发明公开了一种遥感数据角度异常信息提取方法及系统，涉及矿产勘查领域，主要包括首先将获取的工作区的遥感数据按波段形式存储，形成遥感数据集；然后对遥感数据集中的数据进行角度变换处理生成光谱角度数据集；最后采用正交分解技术从光谱角度数据集的数据中提取遥感数据角度异常信息。本发明集合光谱特征和光谱差异两个方面，从角度上进行异常提取，从而获得比遥感数据普通异常提取效果更好的结果。

Description

一种遥感数据角度异常信息提取方法及系统
技术领域
本发明涉及矿产勘查领域，特别是涉及一种遥感数据角度异常信息提取方法及系统。
背景技术
大量实践证明，遥感数据普通异常信息提取在戈壁荒漠区找矿已经具有一定的效果，但提取的遥感数据普通异常信息存在干扰异常多、弱信息难识别等问题。究其原因是因为遥感数据普通异常信息提取仅仅注重了数据光谱特征，而光谱差异的角度信息未能利用。
发明内容
本发明的目的是提供一种遥感数据角度异常信息提取方法及系统，集合光谱特征和光谱差异两个方面，从角度上进行异常提取，从而获得比遥感数据普通异常提取效果更好的结果。
为实现上述目的，本发明提供了如下方案：
一种遥感数据角度异常信息提取方法，包括：
获取工作区的多波段遥感数据；
对所述遥感数据按波段形式存储，形成遥感数据集；其中，工作区的每个波段遥感数据均是由一个带坐标的二维矩阵组成；
对所述遥感数据集中的数据进行角度变换处理，生成光谱角度数据集；
采用正交分解技术，从所述光谱角度数据集的数据中提取遥感数据角度异常信息。
可选的，在采用正交分解技术，从所述光谱角度数据集的数据中提取遥感数据角度异常信息之后，还包括：
在所述多波段的遥感数据上叠加所述遥感角度异常信息，输出显示遥感靶位信息图像。
一种遥感数据角度异常信息提取系统，包括：
工作区遥感数据获取模块，用于获取工作区的多波段遥感数据；
遥感数据集生成模块，用于对所述遥感数据按波段形式存储，形成遥感数据集；其中，工作区的每个波段遥感数据均是由一个带坐标的二维矩阵组成；
光谱角度数据集生成模块，用于对所述遥感数据集中的数据进行角度变换处理，生成光谱角度数据集；
角度异常信息提取模块，用于采用正交分解技术，从所述光谱角度数据集的数据中提取遥感数据角度异常信息。
可选的，还包括：
图像输出模块，用于在所述多波段的遥感数据上叠加所述遥感数据角度异常信息，输出显示遥感靶位信息图像。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明提供了一种遥感数据角度异常信息提取方法及系统，主要包括首先将获取的工作区的遥感数据按波段形式存储，形成遥感数据集；然后对遥感数据集中的数据进行角度变换处理生成光谱角度数据集；最后采用正交分解技术从光谱角度数据集的数据中提取遥感数据角度异常信息。本发明集合光谱特征和光谱差异两个方面，从角度上进行异常提取，从而获得比遥感数据普通异常提取效果更好的结果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明遥感数据角度异常信息提取方法的流程图；
图2为本发明角度异常信息提取区域对应的角度频域直方图；
图3为本发明角度异常信息位置示意图；
图4为本发明遥感靶位信息图；
图5为本发明遥感数据角度异常信息提取系统的结构图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种遥感数据角度异常信息提取方法及系统，集合光谱特征和光谱差异两个方面，从角度上进行异常提取，从而获得比遥感数据普通异常提取效果更好的结果。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
“角度异常”是基于光谱角度数据集异常提取技术提取的异常，能克服数据值普通异常提取造成的错误，结合几何异常提取优势和统计分类优势两种优势，从而达到异常精准识别的目的。主要分为三个步骤：构建角度、角度异常提取、异常切割和优化。
如图1所示，本发明提供的一种遥感数据角度异常信息提取方法包括如下步骤。
步骤101：获取工作区的多波段遥感数据。
获取的工作区(矿区)的多波段遥感数据，可以是高光谱数据，多光谱数据等各类光学遥感数据，如ASTER，LANDSAT等。
步骤102：对所述遥感数据按波段形式存储，形成遥感数据集；其中，工作区的每个波段遥感数据均是由一个带坐标的二维矩阵组成。
对获取的遥感数据按波段存储(输出)，如BAND1表示第一波段，BAND2表示第二波段等等，共计M*N个波段。每个波段都是由一个带坐标的二维矩阵组成，形成的遥感数据集，如下式所示。
Figure BDA0002464957910000031
Figure BDA0002464957910000041
步骤103：对所述遥感数据集中的数据进行角度变换处理，生成光谱角度数据集；具体为：
计算所述遥感数据集中的遥感数据的均值，即：
Figure BDA0002464957910000042
BANDI表示第I波段的二维矩阵数据(即第I波段的遥感数据)，以此均值为原点对所述遥感数据集进行角度变换处理，获得角度数据，即
Figure BDA0002464957910000043
Bi为第i波段的二维矩阵数据中的一个具体数值，
Figure BDA0002464957910000047
表示第j波段转后的光谱角度数据，新生产的光谱角度数据集为
Figure BDA0002464957910000044
Φ1表示第1波段转后的光谱角度数据，也可在一个波段内进行角度变换，如在BANDI波段内进行角度变换如下
Figure BDA0002464957910000045
Bi为第I波段的二维矩阵数据中的一个具体数值，Bijk为第I波段的二维矩阵数据中第j行第k列的一个具体数值，θi为波段内转换的光谱角度数据，新生产的一个波段内的光谱角度数据集为
Figure BDA0002464957910000046
Θ1为第1个波段内遥感数据转后的光谱角度数据。
步骤104：采用正交分解技术，从所述光谱角度数据集的数据中提取遥感数据角度异常信息。
(1)利用偏度系数和峰度系数，在所述光谱角度数据集中筛选出如图3所示的角度异常信息提取区域。
在光谱角度数据集ANGSET和ACRSET上，假设存在一块区域p×q，这个区域内遥感图像的某一个波段像元Φst(s＝1...p；t＝1...q)区间为[x0，xp]，像元均值为
Figure BDA0002464957910000051
标准差为σ，利用偏度系数和峰度系数进行直方图进行判断。其角度异常信息提取区域对应的角度频域直方图如图2所示。
偏度系数满足公式：
Figure BDA0002464957910000052
其中，ε1为给定一个很小的正数。
峰度系数满足公式：
Figure BDA0002464957910000053
其中，ε2为给定一个很小的正数。
那么这个p×q区域就是所选择的数据X。
(2)采用正交分解算法对所述角度异常信息提取区域中的数据进行处理，提取遥感数据角度异常信息。
对选择的数据X，采用正交分解算法提取遥感数据角度异常信息。
原理为：第一步是移动坐标原点，使平均值为零。在这一步骤后，可将坐标旋转，使一个坐标轴与数据具有最大分布的方向相符合，这个旋转后的新轴即第一正交基，它占有总信息量的第一大份额，垂直于它的另一个坐标轴则代表其余信息量的方向，这就是第二正交基。在两维以上的多维空间里，这样的处理将继续进行，以确定一组直角坐标轴，这些直角坐标轴逐渐将全部信息量分配(消耗)掉，它并不能全部包含在一个次一级正交基中，而是有多少个原始参数就会有几个正交基。各正交基的信息量总和与变换前的信息量总和相等，即信息量守恒。
原始有几个波段的数据就影射到几个新的正交基上。各正交基是本征向量线性相加组合而成的。在数学中，就是找出一些新变量ξ1，ξ2，……，ξu，使它们是X的线性函数，而且彼此不相关，即
Figure BDA0002464957910000061
实际上，就是求出u2个常数Lfg(f，g＝1,…，u)按矩阵表示：
Figure BDA0002464957910000062
CL＝λL；
式中：L为本征矩阵，各Lfg是此本征向量的分量；λ为C矩阵的本征值。λ和L有以下特点：
Figure BDA0002464957910000063
称为迹，或总信息量；
对应于不同λ的L(即各主分量)线性不相关。
从线性代数中知道协方差矩阵C的本征多项式为det(λI-C)，此本征多项式的根λ都是协方差矩阵C的本征值。
计算过程为：
协方差矩C：
Figure BDA0002464957910000064
本征值λ：|λI|-C＝0；I为单位矩阵＝
Figure BDA0002464957910000065
本征向量L：(λI-C)L＝0；
在N波段数据坐标轴转置的时候，协方差矩阵也将被变换，变换之后各波段之间的协方差变为零。
各点与其重心距离的平方和即为本征值之和，这个和可表示为S。在一定意义上，可以说第一个正交基“构成”的信息量与总信息量的比例是λ1/S，前两个正交基“构成”的信息量与总信息量之比是(λ1+λ2)/S，如此等等。有时为了方便譬如可以说“前4个分量构成信息量的u％”。
某正交基的本征值就是若消去该正交基后，对相应的本征向量所引入的均方误差值。
对于求出的本征向量，按照与参与正交变换的各个波段进行对应，考虑符合蚀变异常特征的那个本征向量，一般为第4个向量。对应关系如表1所示。
表1 对应关系表
Figure BDA0002464957910000071
如果某一异常的特征为Va4>Vb4<Vc4>Vd4，那么Va4、Vc4一定与Vb4、Vd4的符号相反，而Va4与Vc4、Vb4与Vd4的符号相同。用于异常切割的本征向量4要求Vc4为正号，如果为负号，需要经过转换变成正号，公式如下：
Figure BDA0002464957910000072
其中，Vc4 T为Vc4经过负号转化后的结果。
(3)利用正态分布原理，对所述遥感数据角度异常信息进行优化和分割，得到最终如图3所示的遥感数据角度异常信息。
在正交分解变换前，经过处理后使得每个波段的直方图为正态分布，转化后的异常主分量(即某本征向量)直方图也为正态分布，利用正态分布的有关理论进行异常的切割。正态分布公式如下：
Figure BDA0002464957910000073
其中，X为随机变量，σ称为标准误差。对于正交变换来说，σ称为标准离差，定义如下：
Figure BDA0002464957910000081
n为样本数，
Figure BDA0002464957910000082
为均值，xi为每个样本的值。在做异常切割或数据切割时便可借用σ这个表征正态分布曲线的尺度。例如，正交变换结果可以把均值
Figure BDA0002464957910000083
理解为代表区域背景，利用(X+kσ)确定异常下限和划分异常强度等级。一般取±4σ作为最小与最大限值。
切割异常时有了这一尺度可以减少主观任意性，异常分级是按下式计算的：
L＝127.5+k*SF；或L＝127.5+k*127.5/4；H＝L+1；
式中H、L分别为切割高、低门限值；k为倍数；σ为标准离差；SF为比例因子；σ和SF由正交变换报告给出。
步骤105：在所述多波段的遥感数据上叠加所述遥感数据角度异常信息，输出显示遥感靶位信息图像。
在原始图像上，叠加角度异常信息，输出一幅适合人眼观察的叠加图像，即如图4所示的遥感靶位信息图像。其中，是通过软件输出成为JPG或者TIF格式的最终图像。
如图5所示，本发明还提供了一种遥感数据角度异常信息提取系统，包括：
工作区遥感数据获取模块201，用于获取工作区的多波段遥感数据。
遥感数据集生成模块202，用于对所述遥感数据按波段形式存储，形成遥感数据集；其中，工作区的每个波段遥感数据均是由一个带坐标的二维矩阵组成。
光谱角度数据集生成模块203，用于对所述遥感数据集中的数据进行角度变换处理，生成光谱角度数据集。
角度异常信息提取模块204，用于采用正交分解技术，从所述光谱角度数据集的数据中提取遥感数据角度异常信息。
图像输出模块205，用于在所述多波段的遥感数据上叠加所述遥感数据角度异常信息，输出显示遥感靶位信息图像。
所述光谱角度数据集生成模块203，具体包括：
均值计算单元，用于计算所述遥感数据集中的遥感数据的均值。
光谱角度数据集生成单元，用于以所述均值为原点对所述遥感数据集进行角度变换处理，获得角度数据，进而生成光谱角度数据集。
所述角度异常信息提取模块204，具体包括：
提取区域确定单元，用于利用偏度系数和峰度系数，在所述光谱角度数据集中筛选角度异常信息提取区域。
角度异常信息提取单元，用于采用正交分解算法对所述角度异常信息提取区域中的数据进行处理，提取遥感数据角度异常信息。
角度异常信息优化单元，用于对所述遥感数据角度异常信息进行优化和分割，得到最终的遥感数据角度异常信息。
本发明的目的是一种遥感数据角度异常信息提取方法及系统，为潜在的、隐伏的、微弱的异常信息提取提供了一种新的技术方法。本发明在遥感原始数据，通过构建光谱角度数据集，在光谱角度数据集基础上，采用正交分解技术提取角度异常信息，从而确定遥感靶位信息，能够解决隐伏的、潜在的、微弱的异常靶位信息异常提取难度大、假异常多、提取效果不准确、有偏差等问题，有针对性的改善异常，从而为寻找精确目的提供技术支撑。对于矿产勘查工作可以起到节省时间、节约人力物力，事半功倍的作用，是科学技术进步促进生产发展的新技术。本发明在中国浅覆盖区、干旱区和植被覆盖区发现多个矿点，为中国找矿做出了贡献。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (8)

1.一种遥感数据角度异常信息提取方法，其特征在于，包括：
获取工作区的多波段遥感数据；
对所述遥感数据按波段形式存储，形成遥感数据集；其中，工作区的每个波段遥感数据均是由一个带坐标的二维矩阵组成；
对所述遥感数据集中的数据进行角度变换处理，生成光谱角度数据集；
采用正交分解技术，从所述光谱角度数据集的数据中提取遥感数据角度异常信息。
2.根据权利要求1所述的一种遥感数据角度异常信息提取方法，其特征在于，在采用正交分解技术，从所述光谱角度数据集的数据中提取遥感数据角度异常信息之后，还包括：
在所述多波段的遥感数据上叠加所述遥感角度异常信息，输出显示遥感靶位信息图像。
3.根据权利要求1所述的一种遥感数据角度异常信息提取方法，其特征在于，所述对所述遥感数据集中的数据进行角度变换处理，生成光谱角度数据集，具体包括：
计算所述遥感数据集中的遥感数据的均值；
以所述均值为原点对所述遥感数据集进行角度变换处理，获得角度数据，进而生成光谱角度数据集。
4.根据权利要求1所述的一种遥感数据角度异常信息提取方法，其特征在于，所述采用正交分解技术，对所述光谱角度数据集中的数据进行处理，确定角度异常信息，具体包括：
利用偏度系数和峰度系数，在所述光谱角度数据集中筛选角度异常信息提取区域；
采用正交分解算法对所述角度异常信息提取区域中的数据进行处理，提取遥感数据角度异常信息；
利用正态分布原理，对所述遥感数据角度异常信息进行优化和分割，得到最终的遥感数据角度异常信息。
5.一种遥感数据角度异常信息提取系统，其特征在于，包括：
工作区遥感数据获取模块，用于获取工作区的多波段遥感数据；
遥感数据集生成模块，用于对所述遥感数据按波段形式存储，形成遥感数据集；其中，工作区的每个波段遥感数据均是由一个带坐标的二维矩阵组成；
光谱角度数据集生成模块，用于对所述遥感数据集中的数据进行角度变换处理，生成光谱角度数据集；
角度异常信息提取模块，用于采用正交分解技术，从所述光谱角度数据集的数据中提取遥感数据角度异常信息。
6.根据权利要求5所述的一种遥感数据角度异常信息提取系统，其特征在于，还包括：
图像输出模块，用于在所述多波段的遥感数据上叠加所述遥感数据角度异常信息，输出显示遥感靶位信息图像。
7.根据权利要求5所述的一种遥感数据角度异常信息提取系统，其特征在于，所述光谱角度数据集生成模块，具体包括：
均值计算单元，用于计算所述遥感数据集中的遥感数据的均值；
光谱角度数据集生成单元，用于以所述均值为原点对所述遥感数据集进行角度变换处理，获得角度数据，进而生成光谱角度数据集。
8.根据权利要求5所述的一种遥感数据角度异常信息提取系统，其特征在于，所述角度异常信息提取模块，具体包括：
提取区域确定单元，用于利用偏度系数和峰度系数，在所述光谱角度数据集中筛选角度异常信息提取区域；
角度异常信息提取单元，用于采用正交分解算法对所述角度异常信息提取区域中的数据进行处理，提取遥感数据角度异常信息；
角度异常信息优化单元，用于对所述遥感数据角度异常信息进行优化和分割，得到最终的遥感数据角度异常信息。
CN202010331000.5A 2020-04-24 2020-04-24 一种遥感数据角度异常信息提取方法及系统 Active CN111552004B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202010331000.5A CN111552004B (zh) 2020-04-24 2020-04-24 一种遥感数据角度异常信息提取方法及系统
CH00819/20A CH717361A8 (de) 2020-04-24 2020-07-02 Verfahren und System zum Extrahieren von Winkelanomalie-Informationen aus Fernerkundungsdaten.

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010331000.5A CN111552004B (zh) 2020-04-24 2020-04-24 一种遥感数据角度异常信息提取方法及系统

Publications (2)

Publication Number Publication Date
CN111552004A true CN111552004A (zh) 2020-08-18
CN111552004B CN111552004B (zh) 2023-04-18

Family

ID=72003966

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010331000.5A Active CN111552004B (zh) 2020-04-24 2020-04-24 一种遥感数据角度异常信息提取方法及系统

Country Status (2)

Country Link
CN (1) CN111552004B (zh)
CH (1) CH717361A8 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112288016A (zh) * 2020-10-30 2021-01-29 上海淇玥信息技术有限公司 基于主成分分析算法的渠道反作弊方法、装置和电子设备

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102313526A (zh) * 2010-07-07 2012-01-11 中国科学院地理科学与资源研究所 一种多角度与多光谱遥感数据定量融合反演叶面积指数的方法
CN103926203A (zh) * 2014-04-29 2014-07-16 中国科学院遥感与数字地球研究所 一种针对地物光谱不确定性的光谱角度制图方法
CN104851091A (zh) * 2015-04-28 2015-08-19 中山大学 一种基于卷积增强和hcs变换的遥感影像融合方法
CN107169946A (zh) * 2017-04-26 2017-09-15 西北工业大学 基于非负稀疏矩阵与超球面彩色变换的图像融合方法
CN109145881A (zh) * 2018-10-09 2019-01-04 中国地质科学院矿产资源研究所 一种遥感图像膏盐信息提取方法及装置
CN109961087A (zh) * 2019-02-01 2019-07-02 中国地质科学院矿产资源研究所 基于空间数据集分析的异常遥感信息提取方法及装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102313526A (zh) * 2010-07-07 2012-01-11 中国科学院地理科学与资源研究所 一种多角度与多光谱遥感数据定量融合反演叶面积指数的方法
CN103926203A (zh) * 2014-04-29 2014-07-16 中国科学院遥感与数字地球研究所 一种针对地物光谱不确定性的光谱角度制图方法
CN104851091A (zh) * 2015-04-28 2015-08-19 中山大学 一种基于卷积增强和hcs变换的遥感影像融合方法
CN107169946A (zh) * 2017-04-26 2017-09-15 西北工业大学 基于非负稀疏矩阵与超球面彩色变换的图像融合方法
CN109145881A (zh) * 2018-10-09 2019-01-04 中国地质科学院矿产资源研究所 一种遥感图像膏盐信息提取方法及装置
CN109961087A (zh) * 2019-02-01 2019-07-02 中国地质科学院矿产资源研究所 基于空间数据集分析的异常遥感信息提取方法及装置

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
刁海;张达;狄永军;王振;王浩然;熊光强;: "基于主成分分析和分形模型的ASTER蚀变异常信息提取" *
刘素红,马建文,蔺启忠: "利用掩膜和多因子逐步正交变换区分遥感数据中的岩性信息" *
张玉君，曾朝铭，陈薇: "ETM+ (TM)蚀变遥感异常提取方法研究与应用-方法选择和技术流程", 《国土资源遥感》 *
白彬;: "激光雷达与遥感数据的山区地理信息处理技术" *
马建文，马超飞: "基于空间角度理论的卫星光学遥感数据认知与挖掘", 《中国图象图形学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112288016A (zh) * 2020-10-30 2021-01-29 上海淇玥信息技术有限公司 基于主成分分析算法的渠道反作弊方法、装置和电子设备
CN112288016B (zh) * 2020-10-30 2023-10-31 上海淇玥信息技术有限公司 基于主成分分析算法的渠道反作弊方法、装置和电子设备

Also Published As

Publication number Publication date
CH717361A8 (de) 2021-12-30
CH717361A2 (de) 2021-10-29
CN111552004B (zh) 2023-04-18

Similar Documents

Publication Publication Date Title
CN112101381B (zh) 张量协作图判别分析遥感图像特征提取方法
Nielsen The regularized iteratively reweighted MAD method for change detection in multi-and hyperspectral data
Hou et al. Three-order tucker decomposition and reconstruction detector for unsupervised hyperspectral change detection
CN107992891B (zh) 基于光谱矢量分析多光谱遥感图像变化检测方法
Broadwater et al. Hybrid detectors for subpixel targets
Mundt et al. Partial unmixing of hyperspectral imagery: theory and methods
CN103208011B (zh) 基于均值漂移和组稀疏编码的高光谱图像空谱域分类方法
CN103456020B (zh) 基于treelet特征融合的遥感图像变化检测方法
Lee et al. Exploring Land Cover Classification Accuracy of Landsat 8 Image Using Spectral Index Layer Stacking in Hilly Region of South Korea.
CN104408705A (zh) 一种高光谱图像的异常检测方法
CN103226826A (zh) 基于局部熵视觉注意模型的遥感图像变化检测方法
CN104268561B (zh) 基于结构先验低秩表示的高光谱图像解混方法
CN108898070A (zh) 一种基于无人机平台的高光谱遥感提取薇甘菊装置及方法
CN111552004B (zh) 一种遥感数据角度异常信息提取方法及系统
CN109145881B (zh) 一种遥感图像膏盐信息提取方法及装置
CN113421198B (zh) 一种基于子空间的非局部低秩张量分解的高光谱图像去噪方法
CN112949657B (zh) 一种基于遥感影像纹理特征的林地分布提取方法及装置
An et al. Object-oriented urban dynamic monitoring—A case study of Haidian district of Beijing
Liu et al. Point cloud based iterative segmentation technique for 3d plant phenotyping
Khalid et al. Evaluation the accuracy of oil palm tree detection using deep learning and support vector machine classifiers
Mansourmoghaddam et al. Evaluating the capability of spatial and spectral fusion in land-cover mapping enhancement
CN106033545B (zh) 行列式点过程的波段选择方法
Shen et al. Fractal dimension of irregular region of interest application to corn phenology characterization
CN111985501B (zh) 一种基于自适应高阶张量分解的高光谱图像特征提取方法
CN111091087B (zh) 基于多视图协同典型相关森林遥感影像土地覆盖提取算法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant