阅读说明：本技术 一种遥感数据角度异常信息提取方法及系统 (Method and system for extracting angle abnormal information of remote sensing data ) 是由 姚佛军 耿新霞 杨建民 吴胜华 于 2020-04-24 设计创作，主要内容包括：本发明公开了一种遥感数据角度异常信息提取方法及系统,涉及矿产勘查领域,主要包括首先将获取的工作区的遥感数据按波段形式存储,形成遥感数据集；然后对遥感数据集中的数据进行角度变换处理生成光谱角度数据集；最后采用正交分解技术从光谱角度数据集的数据中提取遥感数据角度异常信息。本发明集合光谱特征和光谱差异两个方面,从角度上进行异常提取,从而获得比遥感数据普通异常提取效果更好的结果。(The invention discloses a method and a system for extracting angle anomaly information of remote sensing data, which relate to the field of mineral exploration and mainly comprise the steps of firstly storing the acquired remote sensing data of a working area according to a wave band form to form a remote sensing data set; then carrying out angle transformation processing on the data in the remote sensing data set to generate a spectrum angle data set; and finally, extracting the angle abnormal information of the remote sensing data from the data of the spectrum angle data set by adopting an orthogonal decomposition technology. The method integrates two aspects of spectral characteristics and spectral differences, and performs abnormal extraction from an angle, so that a result with a better effect than the common abnormal extraction effect of the remote sensing data is obtained.)

1. A method for extracting angle abnormal information of remote sensing data is characterized by comprising the following steps:

obtaining multiband remote sensing data of a working area;

storing the remote sensing data according to a wave band form to form a remote sensing data set; each waveband remote sensing data of the working area consists of a two-dimensional matrix with coordinates;

carrying out angle transformation processing on the data in the remote sensing data set to generate a spectrum angle data set;

and extracting remote sensing data angle abnormal information from the data of the spectrum angle data set by adopting an orthogonal decomposition technology.

2. The method for extracting angular anomaly information of remote sensing data according to claim 1, wherein after extracting angular anomaly information of remote sensing data from data of the spectral angular dataset by using an orthogonal decomposition technique, the method further comprises:

and superposing the remote sensing angle abnormal information on the multiband remote sensing data, and outputting and displaying a remote sensing target position information image.

3. The method for extracting angular anomaly information of remote sensing data according to claim 1, wherein the performing angular transformation processing on the data in the remote sensing data set to generate a spectrum angular data set specifically comprises:

calculating the average value of the remote sensing data in the remote sensing data set;

and carrying out angle transformation processing on the remote sensing data set by taking the mean value as an origin to obtain angle data, and further generating a spectrum angle data set.

4. The method for extracting angular anomaly information of remote sensing data according to claim 1, wherein the step of processing the data in the spectral angular data set by using an orthogonal decomposition technology to determine the angular anomaly information specifically comprises the steps of:

screening an angle abnormal information extraction area in the spectrum angle data set by using a skewness coefficient and a kurtosis coefficient;

processing the data in the angle abnormal information extraction area by adopting an orthogonal decomposition algorithm, and extracting the angle abnormal information of the remote sensing data;

and optimizing and segmenting the remote sensing data angle abnormal information by utilizing a normal distribution principle to obtain final remote sensing data angle abnormal information.

5. The system for extracting the abnormal information of the remote sensing data angle is characterized by comprising the following steps:

the remote sensing data acquisition module of the working area is used for acquiring multiband remote sensing data of the working area;

the remote sensing data set generation module is used for storing the remote sensing data according to a wave band form to form a remote sensing data set; each waveband remote sensing data of the working area consists of a two-dimensional matrix with coordinates;

the spectrum angle data set generating module is used for carrying out angle transformation processing on the data in the remote sensing data set to generate a spectrum angle data set;

and the angle abnormal information extraction module is used for extracting the angle abnormal information of the remote sensing data from the data of the spectrum angle data set by adopting an orthogonal decomposition technology.

6. The remote sensing data angle anomaly information extraction system according to claim 5, further comprising:

and the image output module is used for superposing the remote sensing data angle abnormal information on the multiband remote sensing data and outputting and displaying a remote sensing target position information image.

7. The system for extracting angular anomaly information of remote sensing data according to claim 5, wherein the spectral angular dataset generation module specifically comprises:

the mean value calculating unit is used for calculating the mean value of the remote sensing data in the remote sensing data set;

and the spectrum angle data set generating unit is used for carrying out angle transformation processing on the remote sensing data set by taking the mean value as an origin to obtain angle data so as to generate a spectrum angle data set.

8. The remote sensing data angle anomaly information extraction system according to claim 5, wherein the angle anomaly information extraction module specifically comprises:

an extraction area determination unit, configured to screen an angle abnormal information extraction area in the spectrum angle data set by using a skewness coefficient and a kurtosis coefficient;

the angle abnormal information extraction unit is used for processing the data in the angle abnormal information extraction area by adopting an orthogonal decomposition algorithm and extracting the angle abnormal information of the remote sensing data;

and the angle abnormal information optimizing unit is used for optimizing and segmenting the remote sensing data angle abnormal information to obtain the final remote sensing data angle abnormal information.

Technical Field

The invention relates to the field of mineral exploration, in particular to a method and a system for extracting angle anomaly information of remote sensing data.

Background

A large number of practices prove that the extraction of the common abnormal information of the remote sensing data in the gobi desert region has a certain effect, but the extracted common abnormal information of the remote sensing data has the problems of much interference, difficult identification of weak information and the like. The reason is that the common abnormal information extraction of the remote sensing data only focuses on the data spectrum characteristics, and the angle information of the spectrum difference cannot be utilized.

Disclosure of Invention

The invention aims to provide a method and a system for extracting angle abnormal information of remote sensing data, which integrate two aspects of spectral characteristics and spectral differences to extract the abnormal information from the angle, thereby obtaining a result with better effect than the common abnormal extraction effect of the remote sensing data.

In order to achieve the purpose, the invention provides the following scheme:

a method for extracting angle abnormal information of remote sensing data comprises the following steps:

obtaining multiband remote sensing data of a working area;

storing the remote sensing data according to a wave band form to form a remote sensing data set; each waveband remote sensing data of the working area consists of a two-dimensional matrix with coordinates;

carrying out angle transformation processing on the data in the remote sensing data set to generate a spectrum angle data set;

and extracting remote sensing data angle abnormal information from the data of the spectrum angle data set by adopting an orthogonal decomposition technology.

Optionally, after extracting the remote sensing data angle anomaly information from the data of the spectrum angle data set by using an orthogonal decomposition technique, the method further includes:

and superposing the remote sensing angle abnormal information on the multiband remote sensing data, and outputting and displaying a remote sensing target position information image.

A remote sensing data angle anomaly information extraction system comprises:

the remote sensing data acquisition module of the working area is used for acquiring multiband remote sensing data of the working area;

the remote sensing data set generation module is used for storing the remote sensing data according to a wave band form to form a remote sensing data set; each waveband remote sensing data of the working area consists of a two-dimensional matrix with coordinates;

the spectrum angle data set generating module is used for carrying out angle transformation processing on the data in the remote sensing data set to generate a spectrum angle data set;

and the angle abnormal information extraction module is used for extracting the angle abnormal information of the remote sensing data from the data of the spectrum angle data set by adopting an orthogonal decomposition technology.

Optionally, the method further includes:

and the image output module is used for superposing the remote sensing data angle abnormal information on the multiband remote sensing data and outputting and displaying a remote sensing target position information image.

According to the specific embodiment provided by the invention, the invention discloses the following technical effects:

the invention provides a method and a system for extracting angle abnormal information of remote sensing data, which mainly comprise the steps of firstly storing the acquired remote sensing data of a working area according to a wave band form to form a remote sensing data set; then carrying out angle transformation processing on the data in the remote sensing data set to generate a spectrum angle data set; and finally, extracting the angle abnormal information of the remote sensing data from the data of the spectrum angle data set by adopting an orthogonal decomposition technology. The method integrates two aspects of spectral characteristics and spectral differences, and performs abnormal extraction from an angle, so that a result with a better effect than the common abnormal extraction effect of the remote sensing data is obtained.

Drawings

In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.

FIG. 1 is a flow chart of the method for extracting the abnormal angle information of the remote sensing data according to the present invention;

FIG. 2 is a histogram of angle frequency domain corresponding to the abnormal angle information extraction area according to the present invention;

FIG. 3 is a schematic diagram of the location of the angular anomaly information according to the present invention;

FIG. 4 is a diagram of remote sensing target information according to the present invention;

FIG. 5 is a structural diagram of the remote sensing data angle anomaly information extraction system of the present invention.

Detailed Description

The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.

The invention aims to provide a method and a system for extracting angle abnormal information of remote sensing data, which integrate two aspects of spectral characteristics and spectral differences to extract the abnormal information from the angle, thereby obtaining a result with better effect than the common abnormal extraction effect of the remote sensing data.

In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.

The angle abnormity is extracted based on the spectrum angle data set abnormity extraction technology, errors caused by common abnormity extraction of data values can be overcome, and the purpose of accurate abnormity identification is achieved by combining two advantages of geometric abnormity extraction advantage and statistical classification advantage. The method mainly comprises the following three steps: constructing angles, extracting angle abnormity, cutting abnormity and optimizing.

As shown in fig. 1, the method for extracting angle anomaly information of remote sensing data provided by the invention comprises the following steps.

Step 101: and acquiring multiband remote sensing data of the working area.

The obtained multiband remote sensing data of the working area (mining area) can be various optical remote sensing data such as hyperspectral data, multispectral data and the like, such as ASTER, LANDSAT and the like.

Step 102: storing the remote sensing data according to a wave band form to form a remote sensing data set; wherein, each wave band remote sensing data of the working area is composed of a two-dimensional matrix with coordinates.

The obtained remote sensing data is stored (output) according to wave BANDs, for example, BAND1 represents a first wave BAND, BAND2 represents a second wave BAND, and the like, and the total number of M wave BANDs is N. Each wave band is composed of a two-dimensional matrix with coordinates, and a remote sensing data set is formed, and is shown as the following formula.

Step 103: carrying out angle transformation processing on the data in the remote sensing data set to generate a spectrum angle data set; the method specifically comprises the following steps:

calculating the average value of the remote sensing data in the remote sensing data set, namely:

BANDI represents two-dimensional matrix data (namely, remote sensing data of the I wave band) of the I wave band, and angle transformation processing is carried out on the remote sensing data set by taking the mean value as an origin to obtain angle data, namely, the angle data

Bi is a specific value in the two-dimensional matrix data of the ith waveband,the spectral angle data after the j wave band is converted is represented, and the newly produced spectral angle data set is

Φ₁The spectrum angle data after 1 st wave band conversion can also be subjected to angle transformation in one wave band, such as angle transformation in BANDI wave bandAs follows

Bi is a specific numerical value in the two-dimensional matrix data of the I wave band, Bi_jkIs a specific numerical value of the jth row and kth column in the two-dimensional matrix data of the I waveband, theta I is spectrum angle data converted in the waveband, and a newly produced spectrum angle data set in a waveband is

Θ₁The spectrum angle data converted from the remote sensing data in the 1 st wave band is obtained.

Step 104: and extracting remote sensing data angle abnormal information from the data of the spectrum angle data set by adopting an orthogonal decomposition technology.

(1) And screening out an angle abnormal information extraction area shown in the figure 3 in the spectrum angle data set by using the skewness coefficient and the kurtosis coefficient.

On the spectral angle data sets ANGSET and ACRSET, an area p × q is assumed to exist, and a certain wave band pixel phi of the remote sensing image in the area_stAn interval (s 1.. p; t 1.. q) is [ x₀，x_p]Mean value of picture elements ofAnd the standard deviation is sigma, and the histogram is judged by using the skewness coefficient and the kurtosis coefficient. The angle frequency domain histogram corresponding to the angle anomaly information extraction area is shown in fig. 2.

The skewness coefficient satisfies the formula:

wherein the content of the first and second substances,₁given a very small positive number.

The kurtosis coefficient satisfies the formula:

wherein the content of the first and second substances,₂given a very small positive number.

Then this p × q region is the selected data X.

(2) And processing the data in the angle abnormal information extraction area by adopting an orthogonal decomposition algorithm, and extracting the angle abnormal information of the remote sensing data.

And for the selected data X, extracting the angle abnormal information of the remote sensing data by adopting an orthogonal decomposition algorithm.

The principle is as follows: the first step is to move the origin of coordinates to make the average zero. After this step, the coordinates can be rotated so that one axis coincides with the direction in which the data has the largest distribution, this new axis after rotation being the first orthogonal basis, which occupies the first large share of the total information quantity, and the other axis perpendicular to it represents the direction of the remaining information quantity, which is the second orthogonal basis. In a multidimensional space of more than two dimensions, such processing continues to determine a set of cartesian axes which gradually allocate (consume) all the information, which is not entirely contained in a one-level-one orthogonal basis, but rather which of the original parameters has several orthogonal bases. The sum of the information amounts of the orthogonal bases is equal to the sum of the information amounts before conversion, that is, the information amount is conserved.

The original data with several bands are mapped onto several new orthogonal bases, each of which is formed by linear additive combination of eigenvectors, and in mathematics, some new variables ξ are found₁，ξ₂，……，ξ_uSo that they are linear functions of X and are uncorrelated with each other, i.e.

In practice, u is determined²A constant L_fg(f, g ═ 1, …, u) is expressed in a matrix:

CL＝λL；

in the formula: l is an eigenmatrix, each L_fgIs a component of this eigenvector; λ is the eigenvalue of the C matrix. λ and L have the following characteristics:

referred to as traces, or total information content;

l (i.e., the principal components) corresponding to different λ are linearly uncorrelated.

From linear algebra, the eigen-polynomial of the covariance matrix C is known as det (λ I-C), and the root λ of the eigen-polynomial is the eigenvalue of the covariance matrix C.

The calculation process is as follows:

covariance matrix C:

eigenvalue λ: i λ I-C ═ 0; i is a unit matrix ═

Eigenvector L: (λ I-C) L ═ 0;

when the coordinate axes of the N-band data are transposed, the covariance matrix is also transformed, and the covariance between the bands after transformation becomes zero.

The sum of the squares of the distances of each point from its center of gravity is the sum of the eigenvalues, and this sum can be expressed as S. In a sense, it can be said that the ratio of the amount of information "made up" of the first orthogonal basis to the total amount of information is λ₁S, the ratio of the information quantity "formed" by the first two orthogonal bases to the total information quantity is (lambda)₁+λ₂) S, and so on. For example, "the first 4 components constitute u% of the information amount" may be said for convenience.

The eigenvalues of an orthogonal basis are the mean square error values introduced into the corresponding eigenvectors if the orthogonal basis is eliminated.

The eigenvectors obtained are considered corresponding to the respective bands participating in the orthogonal transformation, and those eigenvectors corresponding to the abnormal characteristics of alteration, generally the 4 th vector, are considered. The correspondence is shown in table 1.

Table 1 table of correspondence

If an anomaly is characterized by V_a4>V_b4<V_c4>V_d4Then V_a4、V_c4A constant and V_b4、V_d4Are of opposite sign, and V_a4And V_c4、V_b4And V_d4Are the same as the symbols in (a). Eigenvector 4 requirement V for abnormal cleavage_c4Is positive, if negative, it needs to be converted into positive, and the formula is as follows:

wherein, V_c4 ^TIs a V_c4Results after negative conversion.

(3) And optimizing and segmenting the remote sensing data angle abnormal information by utilizing a normal distribution principle to obtain the remote sensing data angle abnormal information shown in the figure 3 finally.

Before orthogonal decomposition transformation, the histogram of each wave band is made to be normal distribution after processing, the transformed abnormal principal component (namely, a certain eigenvector) histogram is also normal distribution, and abnormal cutting is carried out by utilizing the related theory of normal distribution. The normal distribution formula is as follows:

where X is a random variable and σ is referred to as the standard error. For an orthogonal transformation, σ is called the standard deviation, and is defined as follows:

n is the number of samples and n is the number of samples,is a mean value, x_iIs the value of each sample. The scale of the sigma representing the normal distribution curve can be borrowed when abnormal cutting or data cutting is carried out. For example, the orthogonal transformation result may be the mean valueAs understood as representing the regional background, the lower abnormality limit is determined and the abnormality intensity level is classified using (X + k σ). Typically, a minimum and maximum limit of ± 4 σ is taken.

When cutting abnormity, the scale can reduce subjective arbitrariness, and the abnormity grading is calculated according to the following formula:

l127.5 + k SF; or L127.5 + k 127.5/4; h ═ L + 1;

h, L are cut high and low threshold values, respectively; k is a multiple; σ is the standard deviation; SF is a scale factor; σ and SF are given by the orthogonal transform report.

Step 105: and superposing the remote sensing data angle abnormal information on the multiband remote sensing data, and outputting and displaying a remote sensing target position information image.

And superposing the angle abnormal information on the original image, and outputting a superposed image suitable for human eyes to observe, namely the remote sensing target position information image shown in figure 4. The final image is output by software in a JPG or TIF format.

As shown in fig. 5, the present invention further provides a system for extracting angle anomaly information of remote sensing data, including:

and the working area remote sensing data acquisition module 201 is used for acquiring multiband remote sensing data of the working area.

The remote sensing data set generation module 202 is used for storing the remote sensing data according to a wave band form to form a remote sensing data set; wherein, each wave band remote sensing data of the working area is composed of a two-dimensional matrix with coordinates.

And the spectrum angle data set generating module 203 is used for performing angle transformation processing on the data in the remote sensing data set to generate a spectrum angle data set.

And the angle abnormal information extraction module 204 is used for extracting the remote sensing data angle abnormal information from the data of the spectrum angle data set by adopting an orthogonal decomposition technology.

And the image output module 205 is used for superposing the remote sensing data angle abnormal information on the multiband remote sensing data and outputting and displaying a remote sensing target position information image.

The spectrum angle data set generating module 203 specifically includes:

and the mean value calculating unit is used for calculating the mean value of the remote sensing data in the remote sensing data set.

And the spectrum angle data set generating unit is used for carrying out angle transformation processing on the remote sensing data set by taking the mean value as an origin to obtain angle data so as to generate a spectrum angle data set.

The angle anomaly information extraction module 204 specifically includes:

and the extraction area determining unit is used for screening an angle abnormal information extraction area in the spectrum angle data set by using the skewness coefficient and the kurtosis coefficient.

And the angle abnormal information extraction unit is used for processing the data in the angle abnormal information extraction area by adopting an orthogonal decomposition algorithm and extracting the angle abnormal information of the remote sensing data.

And the angle abnormal information optimizing unit is used for optimizing and segmenting the remote sensing data angle abnormal information to obtain the final remote sensing data angle abnormal information.

The invention aims to provide a method and a system for extracting angle abnormal information of remote sensing data, and provides a new technical method for extracting potential, hidden and weak abnormal information. According to the method, the remote sensing target information is determined by constructing the spectrum angle data set on the basis of the remote sensing original data and extracting the angle abnormal information by adopting the orthogonal decomposition technology on the basis of the spectrum angle data set, the problems of high difficulty in extracting the hidden, potential and weak abnormal target information, more false abnormality, inaccurate extraction effect, deviation and the like can be solved, the abnormality is improved in a targeted manner, and the technical support is provided for the purpose of accurately searching. The method has the advantages of saving time, manpower and material resources and achieving the effect of double results with half the effort on mineral exploration work, and is a new technology for promoting production development by scientific and technical progress. The invention finds a plurality of mine points in the shallow coverage area, the arid area and the vegetation coverage area of China, and makes contribution to the exploration of mines in China.

The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.

The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.