A method and system for processing remote sensing mapping data of a building house

Abstract

The present invention relates to the technical field of remote sensing mapping, and particularly relates to a method and system for processing remote sensing mapping data of building houses. By acquiring hyperspectral data at different sampling times in the area to be measured, combining the fluctuation of the wavelength absorption rate at different sampling points with the difference in vegetation density, a density difference parameter is obtained. By quantifying the contribution degree of the spectral data at different sampling points to the building house density, and based on the density difference parameter, analyzing the difference in the wavelength absorption rate to obtain the retention degree value of each wavelength, which is used to evaluate the importance of each wavelength. Using this information to correct the initial covariance matrix, increase the weight of the key wavelengths, and obtain an updated covariance matrix. Finally, based on principal component analysis, dimensionality reduction analysis is performed on the updated covariance matrix to ensure that the dimensionality reduction result can more accurately reflect the building house density situation, providing more valuable information for urban management.

Description

Technical Field

[0001] The present invention relates to the technical field of remote sensing mapping, and particularly relates to a method and system for processing remote sensing mapping data of building houses. Background Art

[0002] With the rapid development of remote sensing technology, hyperspectral imaging technology has become an important means for obtaining information on building houses in the field of mapping. Hyperspectral data provides rich spectral information with its high-dimensional characteristics. However, the dimension of hyperspectral data is usually very high. For example, there may be hundreds or even thousands of wavelength data for each location point. Therefore, the dimension of hyperspectral data can be reduced to extract the core information of the original data, thereby improving the efficiency and accuracy of remote sensing monitoring of building houses.

[0003] When the prior art reduces the dimension of hyperspectral data, it usually performs undifferentiated dimension reduction on the hyperspectral data by using the principal component analysis method. However, due to various environmental factors at different times, the absorption rate of wavelengths in the hyperspectral data will be affected. If undifferentiated dimension reduction is used, the dimension reduction result will lose too much key information, resulting in the dimension reduction result being unable to accurately characterize the density status of the building house itself. Summary of the Invention

[0004] In order to solve the technical problem that undifferentiated dimension reduction will cause the dimension reduction result of the hyperspectral data of building houses to lose too much key information, resulting in the dimension reduction result being unable to accurately characterize the density status of the building house itself, the purpose of the present invention is to provide a method and system for processing remote sensing mapping data of building houses, and the specific technical solutions adopted are as follows:

[0005] Obtain the hyperspectral data of the area to be measured at different sampling times;

[0006] For the hyperspectral data at any one sampling time, according to the difference in the vegetation density between different sampling points, the difference in the absorption rate of wavelengths between different sampling points, and the fluctuation of the absorption rate of all sampling points at the same wavelength, obtain the density difference parameter between different sampling points in the hyperspectral data at this sampling time;

[0007] For any one sampling point in the area to be measured, fuse and analyze the positional relationship, the difference in spectral data, and the density difference parameter between this sampling point and the remaining other sampling points at all sampling times to obtain the abnormal degree value of the building house density at this sampling point; for any one wavelength, fuse and analyze the difference in the absorption rate of this wavelength at all sampling points at different times and the abnormal degree value of the building house density of each sampling point, and then analyze the fluctuation situation to obtain the retention degree value of this wavelength;

[0008] Based on the hyperspectral data of all sampling points at all sampling times, an initial covariance matrix corresponding to the spectral data is calculated; the element values in the initial covariance matrix are corrected according to the retention degree values of each wavelength to obtain an updated covariance matrix; principal component analysis is performed on the updated covariance matrix to obtain a dimensionality reduction result.

[0009] Further, the method for obtaining the density difference parameter includes:

[0010] Obtain the normalized difference vegetation index at each sampling point in the area to be measured;

[0011] At any one sampling time, the difference in the normalized difference vegetation index between two different sampling points is used as the vegetation condition difference value;

[0012] Analyze the difference between the absorption rates of two different sampling points at the same wavelength and the fluctuation of the absorption rates of all sampling points at the same wavelength to obtain the density difference factor between different sampling points at this sampling time;

[0013] According to the vegetation condition difference value and the density difference factor between two different sampling points, the density difference parameter between two different sampling points at this sampling time is obtained. The density difference parameter is positively correlated with the density difference factor, the density difference parameter is negatively correlated with the vegetation condition difference value, and the value of the density difference parameter is a normalized value.

[0014] Further, the method for obtaining the density difference factor includes:

[0015] Optionally select a wavelength as the wavelength to be measured. At this sampling time, the difference in the absorption rate of two different sampling points at the wavelength to be measured is used as the absorption rate difference value; the standard deviation of the absorption rates of all sampling points at the wavelength to be measured is used as the absorption rate fluctuation value;

[0016] According to the absorption rate fluctuation value and the absorption rate difference value, an absorption rate change index of two sampling points at the wavelength to be measured is obtained. The absorption rate fluctuation value and the absorption rate difference value both change in the same direction as the absorption rate change index; the average value of the absorption rate change indices of two different sampling points at all wavelengths is used as the density difference factor.

[0017] Further, the method for obtaining the abnormal degree value of the building housing density includes:

[0018] Optionally select a sampling point as the target point, and according to the distance between the target point and each remaining sampling point, screen out the comparison points of the target point among the remaining sampling points;

[0019] For any sampling moment, the difference in the absorbance at the same wavelength in the spectral data of the target point and each comparison point is used as the absorbance distance value; the sum value of the absorbance distance values at all wavelengths is used as the sum distance value. Based on the density difference parameter between the target point and each comparison point and the sum distance value, the deviation value of the spectral data of the target point and each comparison point at this sampling moment is obtained. The deviation value is negatively correlated with the density difference parameter, and the sum distance value and the deviation value are positively correlated;

[0020] The mean value of the deviation values of the spectral data of the target point and each comparison point at all sampling moments is used as the abnormal factor of the building housing density of the target point; the mean value of the abnormal factors of the building housing density of the target point and all comparison points is used as the abnormal degree value of the building housing density at the target point.

[0021] Further, the method for obtaining the comparison points includes:

[0022] Calculate the Euclidean distance between the target point and the position of each remaining sampling point as the distance factor, and use the sampling points with the distance factor less than the preset distance threshold as the comparison points of the target point.

[0023] Further, the method for obtaining the retention degree value includes:

[0024] Combine all sampling moments in pairs to obtain all non-repeating moment combinations;

[0025] For any sampling point, in each moment combination, fuse the abnormal degree value of the building housing density at this sampling point with the difference in the absorbance of the wavelength to be measured at the two sampling moments of this sampling point to obtain the absorbance change characteristic value of this sampling point;

[0026] The standard deviation of the absorbance change characteristic values of the wavelength to be measured at the same moment combination for all sampling points is used as the retention factor of the wavelength to be measured at the same moment combination; the mean value of the retention factors of the wavelength to be measured at all moment combinations is used as the retention degree value of the wavelength to be measured.

[0027] Further, the method for obtaining the absorbance change characteristic value includes:

[0028] The value obtained by performing negative correlation mapping and normalization on the abnormal degree value of the building housing density at each sampling point is used as the retention weight; for any sampling point, in each moment combination, the difference in the absorbance of the wavelength to be measured at the two sampling moments of this sampling point is used as the absorbance change factor; the product of the absorbance change factor of this sampling point in each moment combination and the corresponding retention weight is used as the absorbance change characteristic value of this sampling point at the wavelength to be measured.

[0029] Further, the method for obtaining the initial covariance matrix includes:

[0030] The hyperspectral data includes two position data and one spectral data;

[0031] Taking each spectral data as a column vector, the length of the column vector is equal to the number of wavelengths, and the element value is the wavelength absorption rate;

[0032] Forming an initial matrix with the column vectors corresponding to the spectral data of all sampling points at all sampling times, and calculating the initial matrix to obtain an initial covariance matrix.

[0033] Further, the method for obtaining the updated covariance matrix includes:

[0034] Constructing a correction matrix based on the retention degree values of each wavelength, the correction matrix is a diagonal matrix, and each element value on the diagonal is the retention degree value of each wavelength, and the element values at the remaining positions are 0;

[0035] Multiplying the matrix obtained by multiplying the correction matrix by the initial covariance matrix by the correction matrix again to obtain the updated covariance matrix.

[0036] The present invention also provides a system for processing remote sensing mapping data of a building house, and the system includes:

[0037] A memory, a processor, and a computer program stored in the memory and executable on the processor, and when the processor executes the computer program, the steps of any one of the methods can be implemented.

[0038] The present invention has the following beneficial effects:

[0039] The object of the present invention is to improve the importance of wavelengths reflecting building houses, so as to improve the accuracy of the dimensionality reduction result in characterizing the density status of building houses in the subsequent dimensionality reduction process. First, hyperspectral data of the area to be measured at different sampling times is obtained. Since even at the same sampling time, if only the difference in wavelength absorption rates between sampling points is analyzed, it is impossible to distinguish well whether the difference is affected by the density of vegetation. Therefore, when analyzing the difference in wavelength absorption rates between different sampling points and the fluctuation of absorption rates of all sampling points at the same wavelength, the difference in the density of vegetation between different sampling points is combined therewith to obtain the difference parameter of the building house density between different sampling points at each sampling time. This parameter more accurately reflects the difference in the density of building houses between different sampling points considering the influence of the density of vegetation. Further, in order to determine whether each wavelength determines the density of building houses, it is necessary to analyze the contribution degree of the spectral data of different sampling points to the density of building houses. Since the difference in the spectral data of sampling points with close positions may be caused by the different density of roads between building houses, on the basis of the above-mentioned difference parameter of building house density, combining the positional relationship between each sampling point and the difference between spectral data, the abnormal degree value of the building house density at each sampling point is obtained, which is used to quantify the contribution degree of the spectral data at this sampling point to the density of building houses. Further, by performing a fusion analysis on the difference in absorption rates of each wavelength at different sampling points at different times and the abnormal degree value of building house density and evaluating its fluctuation, the retention degree value of each wavelength can be obtained, which is used to distinguish the wavelengths with high research value and low research value for the density status of building houses. Then, based on all the hyperspectral data, the initial covariance matrix is calculated, and the element values in the initial covariance matrix are corrected based on the retention degree values of each wavelength, so that the wavelengths with high retention degree values have higher weights, and the updated covariance matrix is obtained. Finally, principal component analysis is continued on the updated covariance matrix, and the dimensionality reduction result obtained will contain more key information and can more accurately characterize the density status of building houses. Description of the Drawings

[0040] In order to more clearly illustrate the technical solutions and advantages in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required for the description of the embodiments or the prior art. Obviously, the following drawings are only some embodiments of the present invention. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0041] Figure 1 Flowchart of a method for processing remote sensing mapping data of building houses provided by an embodiment of the present invention;

[0042] Figure 2Flowchart of a method for obtaining density difference parameters provided by an embodiment of the present invention;

[0043] Figure 3 Flowchart of a method for obtaining the abnormal degree value of the density of a building provided by an embodiment of the present invention;

[0044] Figure 4 System block diagram of a building remote sensing mapping data processing system provided by an embodiment of the present invention;

[0045] Figure 5 Schematic diagram of a system provided by an embodiment of the present invention;

[0046] Reference numerals: 401, data acquisition module; 402, density difference analysis module; 403, wavelength retention degree analysis module; 404, dimensionality reduction module; 500, processor; 501, memory; 502, bus; 503 communication interface. Detailed implementation manners

[0047] In order to further elaborate on the technical means and effects adopted by the present invention to achieve the predetermined invention purpose, the following specifically describes, with reference to the accompanying drawings and preferred embodiments, a method and system for processing building remote sensing mapping data according to the present invention, including its specific implementation manners, structures, features, and effects. In the following description, different "one embodiment" or "another embodiment" do not necessarily refer to the same embodiment. In addition, the specific features, structures, or characteristics in one or more embodiments can be combined in any suitable form.

[0048] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by those skilled in the technical field to which the present invention belongs.

[0049] The following specifically describes the specific solutions of a method and system for processing building remote sensing mapping data provided by the present invention with reference to the accompanying drawings.

[0050] Please refer to Figure 1 , which shows a flowchart of a method for processing building remote sensing mapping data provided by an embodiment of the present invention. The method includes the following steps:

[0051] Step S1: Obtain hyperspectral data of the area to be measured at different sampling times.

[0052] With the rapid development of remote sensing technology, hyperspectral imaging technology has evolved into a core tool for obtaining building information in the field of surveying and mapping. Although hyperspectral data provides rich spectral information due to its high-dimensional characteristics, it also brings challenges such as large data volume, high information redundancy, and increased computational complexity. Machine learning technology has achieved remarkable results in the field of data mining, especially in feature extraction and dimensionality reduction. Therefore, the dimensionality reduction method based on machine learning can automatically extract the most representative feature information by learning the inherent laws and features in the original data, thereby reducing the data dimension while retaining the core information of the original data. This method can be applied to the dimensionality reduction process of building hyperspectral data.

[0053] First, we can use aerial photography to obtain hyperspectral data of the area to be tested at different sampling times using a hyperspectral imager, where the sampling interval is set to 2 hours, the wavelength range of the spectrum is 400-1000 nanometers, and the wavelength resolution is set to 1 nanometer.

[0054] It should be noted that hyperspectral data contains data in three dimensions, of which two spatial dimensions can be used to determine the location of the sampling point, and the spectral dimension contains the spectral data of the sampling point, the horizontal axis of the spectral data is the wavelength, and the vertical axis is the absorbance. The specific hyperspectral data collection method, wavelength range, resolution, sampling time interval, etc. can be adjusted according to the implementation scenario and are not limited here.

[0055] Step S2: For the hyperspectral data at any sampling time, the building density difference parameters between different sampling points in the hyperspectral data at the sampling time are obtained according to the differences in vegetation density between different sampling points, the differences in the absorbance of different wavelengths at different sampling points, and the fluctuations in the absorbance of all sampling points at the same wavelength.

[0056] In the embodiment of the present invention, the principal component analysis method is used to perform data dimensionality reduction. However, since not every wavelength in the spectral data has a good reflection effect on the density of buildings, if all wavelengths in the spectral data are subjected to indiscriminate dimensionality reduction, the final dimensionality reduction result will not be able to well characterize the density of buildings, so different wavelengths need to be distinguished.

[0057] Since the hyperspectral data of the same area to be measured is collected at different times, the sampling points in the hyperspectral data at different times will correspond one by one. That is, the same sampling point will have spectral data at multiple different times. It should be noted that the SIFT feature point matching method can also be used to match the sampling points in the hyperspectral data at different times to obtain the spectral data at the same sampling point at different times. And the SIFT feature point matching algorithm is a well-known technical means in the art, so it will not be elaborated here.

[0058] Because in the hyperspectral data of the area to be measured, there are some vegetation parts without buildings at some sampling points, so just based on the difference in spectral data between two points, it is impossible to better distinguish whether the difference in spectral data is caused by the different density of buildings or the density of vegetation. Therefore, considering the difference in the density of vegetation between sampling points, then analyzing the difference in wavelength absorption rate between sampling points, and the fluctuation of the absorption rate of all sampling points at the same wavelength, so as to obtain the density difference parameter between different sampling points, which is used to measure the difference in the density of building houses between sampling points.

[0059] Preferably, in an embodiment of the present invention, the method for obtaining the density difference parameter includes:

[0060] Please refer to Figure 2 , which shows the flowchart of a method for obtaining a density difference parameter provided by an embodiment of the present invention. The method includes the following steps:

[0061] Step S201: Based on the difference in the density of vegetation at different sampling points in the area to be measured, determine the difference value of the vegetation status between different sampling points.

[0062] According to the spectral characteristics of vegetation, the satellite visible light and near-infrared bands are combined to form various vegetation indices. The vegetation index is a simple, effective and empirical measure of the surface vegetation status. Therefore, the normalized vegetation index at each sampling point in the area to be measured can be obtained , which is used to quantify the density of vegetation at each sampling point.

[0063] Then at any sampling time, the difference between the normalized vegetation indices between two different sampling points is used as the difference value of the vegetation status. The formula model of the difference value of the vegetation status is:

[0064]

[0065] Wherein, represents at the th sampling time, sampling point and sampling point The difference value of the vegetation condition between them; It represents the th sampling moment, and the normalized difference vegetation index at the sampling point ; It represents the th sampling moment, and the normalized difference vegetation index at the sampling point ;

[0066] In the formula model of the vegetation condition difference value, if there is a significant difference in the vegetation density at two sampling points. For example, if the vegetation is dense at one sampling point while the other sampling point belongs to the construction area, the vegetation condition difference value between these two sampling points will be larger; conversely, if the vegetation density at the two sampling points is relatively consistent, the vegetation condition difference value will be correspondingly smaller.

[0067] It should be noted that the method for obtaining the normalized difference vegetation index is a well-known technical means to those skilled in the art and will not be elaborated here. And in other embodiments of the present invention, when calculating the vegetation condition difference value, the normalized difference vegetation index can also be replaced with the ratio vegetation index, the greenness vegetation index, etc., and the obtaining methods are all processes well-known to those skilled in the art and will not be limited and elaborated here.

[0068] Step S202: Analyze the difference between the absorptions of different sampling points at the same wavelength and the fluctuation of the absorptions of all sampling points at the same wavelength to obtain the density difference factor between different sampling points at this sampling moment.

[0069] To quantify the density difference between different sampling points, a wavelength can be arbitrarily selected as the wavelength to be measured first.

[0070] Then at this sampling moment, the difference between the absorptions of two different sampling points at the wavelength to be measured is used as the absorption difference value; the standard deviation of the absorptions of all sampling points at the wavelength to be measured is used as the absorption fluctuation value. Both the absorption difference value and the absorption fluctuation value characterize the change of the absorption of the sampling point at the same wavelength. Therefore, by combining the two, an absorption change index of the two sampling points at the wavelength to be measured is obtained according to the absorption fluctuation value and the absorption difference value, and both the absorption fluctuation value and the absorption difference value change in the same direction as the absorption change index.

[0071] Finally, the average value of the absorption change indexes of two different sampling points at all wavelengths is used as the density difference factor. The formula model of the density difference factor can be specifically, for example:

[0072]

[0073] Among them, represents at the At a sampling moment, the sampling points and the sampling point The density difference factor between them; Indicates the total number of wavelengths; Indicates the absorption rate of the wavelength; Indicates at the th wavelength, the absorption rate fluctuation value of all sampling points; Indicates at the th sampling moment, the sampling point At the th wavelength, the absorption rate; Indicates at the th sampling moment, the sampling point At the th wavelength, the absorption rate.

[0074] In the formula model of the density difference factor, at a certain sampling moment, from a global analysis, calculate the standard deviation of the absorption rates of all sampling points at the same wavelength , as the absorption rate fluctuation value, used to measure the change of the absorption rate at the same wavelength. The larger this value is, it means that in the overall situation, there are relatively significant differences in the absorption rates of all sampling points at this wavelength. Then, there is a greater possibility of differences in the density corresponding to this wavelength between the two sampling points. Then, from a local single wavelength analysis, calculate the difference in the absorption rates of the two sampling points at the same wavelength, as the absorption rate difference value . The larger this value is, it means that there are greater differences in the density corresponding to this wavelength between the two sampling points. Then, combine the absorption rate fluctuation value at each wavelength and the absorption rate difference value between the two sampling points, and take the product of the two as the absorption rate change index of the two sampling points at each wavelength . The larger this index is, it means that there are greater density differences between the two sampling points in the substance corresponding to this wavelength. Finally, take the average value of the absorption rate change indexes of the two sampling points at all wavelengths as the density difference factor between the two sampling points. The larger the density difference factor is, it means that there are greater differences in the density between the two sampling points.

[0075] Since both the absorption rate fluctuation value and the absorption rate difference value change in the same direction as the absorption rate change index, in other embodiments of the present invention, the absorption rate change index can also be expressed as , and the specific calculation method is not limited here.

[0076] Step S203: According to the vegetation condition difference value and the density difference factor between different sampling points, obtain the density difference parameter between different sampling points at this sampling moment.

[0077] Since the density difference factor only analyzes the differences in the absorption rates of different wavelengths in the spectral data of different sampling points, and there are obvious differences in the spectral data of sampling points with different vegetation densities, in order to reduce the impact of this situation, the vegetation condition difference value obtained in step S201 and the density difference factor obtained in step S202 should be combined to obtain the density difference parameter between different sampling points. The formula model of the density difference parameter can be specifically, for example:

[0078]

[0079] where, represents the density difference parameter between sampling point and sampling point at the th sampling moment; represents the vegetation condition difference value between sampling point and sampling point at the th sampling moment; represents the density difference factor between sampling point and sampling point at the th sampling moment; represents a preset first parameter; represents a normalization function.

[0080] In the formula model of the density difference parameter, when there are significant differences in the vegetation density at two sampling points, the vegetation condition difference value at these two sampling points will be larger, and then the difference in density between these two sampling points is more likely to be caused by the difference in vegetation coverage. That is, if the vegetation condition difference value between two sampling points is larger, the credibility of the density difference factor between these two sampling points should be reduced; conversely, if the vegetation density at two sampling points is relatively consistent, the vegetation condition difference value will be correspondingly smaller, and the credibility of the density difference factor between these two sampling points can be improved. Based on the foregoing logic, the vegetation condition difference value is used as the denominator part, and the density difference factor is used as the numerator part to obtain the density difference parameter between two sampling points.

[0081] It should be noted that the role of the preset first parameter is to prevent the denominator from being 0. Here, it can take the value of 0.001, and the specific value can be adjusted according to the actual situation and is not limited here.

[0082] Step S3: For any sampling point in the area to be measured, fuse and analyze the positional relationships between this sampling point and the remaining other sampling points at all sampling times, the differences between spectral data, and the density difference parameters to obtain the abnormal degree value of the building density at this sampling point; for any wavelength, fuse and analyze the fluctuations of the differences in the absorbance of this wavelength at different times for all sampling points and the abnormal degree value of the building density at each sampling point to obtain the retention degree value of this wavelength.

[0083] Based on the method in Step S2, at any sampling time, the density difference parameter between any two sampling points at this sampling time can be obtained. To obtain the wavelength that can represent the density of buildings, it is necessary to analyze the contribution degree of spectral data at different positions to the building density. Moreover, in the hyperspectral data obtained by aerial photography, the differences in spectral data of sampling points with similar positions may be caused by the sparsity of roads between buildings, and such differences cannot well show the difference relationship between the building density at this sampling point and the normal building density. Therefore, after obtaining the density difference parameter of the building density between sampling points, determine the local range of each sampling point, consider the distance metric between spectral data between sampling points within this local range, and combine the density difference parameter between sampling points, thereby calculating the abnormal degree value of the building density at each sampling point to measure the contribution degree of the spectral data at this sampling point to the analysis of the building density.

[0084] Preferably, in an embodiment of the present invention, the method for obtaining the abnormal degree value of the building density includes:

[0085] Please refer to Figure 3 , which shows the flowchart of a method for obtaining the abnormal degree value of the building density provided by an embodiment of the present invention. The method includes the following steps:

[0086] Step S301: Determine the comparison point of each sampling point according to the positional relationship between sampling points.

[0087] By selecting the comparison point of each sampling point through the positional relationship between sampling points, it can be ensured that the sampling point and the selected comparison point are relatively close in space, thereby increasing the spatial representativeness of the data. Such a comparison point is more likely to have similar environmental conditions or influencing factors with the corresponding sampling point, increasing the accuracy of the subsequent process.

[0088] Optionally select a sampling point as the target point, calculate the Euclidean distance between the target point and the positions of the remaining each sampling point as the distance factor, and then use the sampling points with the distance factor less than the preset distance threshold as the comparison points of the target point.

[0089] It should be noted that the calculation formula of the Euclidean distance is a process well-known to those skilled in the art and will not be elaborated here; the preset distance threshold is set to 7, and the specific value can be adjusted according to the implementation scenario and is not limited here; in other embodiments of the present invention, the value obtained by normalizing the Euclidean distance can also be used as the distance factor. At this time, the value of the distance factor ranges from 0 to 1, and the closer it is to 1, the farther the distance between the two sampling points. The distance threshold is used to screen out the comparison points closer to the sampling point. Therefore, the preset distance threshold can be set to 0.5, and the smaller the value, the stricter the screening conditions for the comparison points, and vice versa.

[0090] Step S302: Determine the abnormal degree value of the building housing density at each sampling point according to the difference between the spectral data of each sampling point and its comparison point at all sampling times and the building housing density difference parameter.

[0091] For any sampling time, the difference in the absorption rate of the same wavelength in the spectral data of the target point and each comparison point is used as the absorption rate distance value, and then the sum value of the absorption rate distance values of all wavelengths is calculated as the distance sum value. This distance sum value can be regarded as the distance between the two spectral data. Then, according to the density difference parameter and the distance sum value between the target point and each comparison point, the deviation value of the spectral data of the target point and each comparison point at this sampling time is obtained, and the deviation value is negatively correlated with the density difference parameter, and the distance sum value and the deviation value are positively correlated. Among them, the positive correlation means that the dependent variable will increase as the independent variable increases, and the dependent variable will decrease as the independent variable decreases. The specific relationship can be a multiplicative relationship, an additive relationship, the power of an exponential function, etc., which is determined by the actual application; the negative correlation means that the dependent variable will decrease as the independent variable increases, and the dependent variable will increase as the independent variable decreases, which can be a subtractive relationship, a divisive relationship, etc., and is determined by the actual application.

[0092] The average value of the deviation values of the spectral data of the target point and each comparison point at all sampling times is used as the abnormal factor of the building housing density of the target point; the average value of the abnormal factors of the building housing density of the target point and all comparison points is used as the abnormal degree value of the building housing density at the target point. The target point is taken as the sampling point as an example, the formula model of the abnormal degree value of the building housing density is:

[0093]

[0094] Among them, represents the abnormal degree value of the building housing density of the sampling point ; represents the total number of comparison points of the sampling point ; represents the total number of sampling times; represents at the At a sampling moment, the sampling point and the density difference parameter between the It represents at the sampling moment, the sampling point and the sum of distances between the spectral data of the reference points;

[0095] In the formula model of the abnormal degree value of the building housing density, when the distance between the spectral data of the sampling point and its reference point is large and the density difference parameter is small, the deviation value is larger, indicating that the difference in spectral data between the two sampling points is more likely to be caused by the different density of the roads between the buildings at the two sampling points. Such data changes cannot show the density of the buildings at the sampling point, so the contribution degree is low, and its abnormal degree is regarded as high. Then, the mean value of the deviation values between the target point and its reference points at all sampling moments is calculated as the abnormal factor of the building housing density at the target point , and finally, the mean value of the abnormal factors of the building housing density of the target point and all its corresponding reference points is used as the abnormal degree value of the building housing density at the target point. The larger this value is, the lower the contribution degree of the spectral data of the sampling point to the analysis of the building housing density, that is, the higher the abnormal degree of the building housing density.

[0096] It should be noted that the preset second parameter is used to prevent the denominator from being zero. Here, it can take the value of 0.001, and the specific value can be adjusted according to the implementation scenario and is not limited here.

[0097] So far, the abnormal degree value of the building housing density at each sampling point can be obtained, which is used to characterize the contribution degree of its spectral data to the building housing density.

[0098] Since the environmental factors change at different sampling moments, such as the weather changes, the water content in the air changes, resulting in the change of the environmental humidity. This environmental change will affect the absorption rate of some wavelengths in the spectral data. At the same time, because the influence caused by the environmental change is relatively uniform, the change of the wavelength absorption rate caused by the environmental factors at different moments will be different from the change of the wavelength absorption rate caused by the building housing density. Therefore, based on this feature, combined with the abnormal degree value of the building housing density at the sampling point obtained above, the research value of each wavelength for studying the building housing density is quantified, that is, the retention degree value of each wavelength is obtained.

[0099] Preferably, in an embodiment of the present invention, the method for obtaining the retention degree value includes:

[0100] All sampling moments are combined in pairs to obtain all non-repeating moment combinations. For example, if there are 4 sampling moments, denoted as sampling moment 1, sampling moment 2, sampling moment 3, and sampling moment 4 respectively, then all moment combinations include: (sampling moment 1, sampling moment 2), (sampling moment 1, sampling moment 3), (sampling moment 1, sampling moment 4), (sampling moment 2, sampling moment 3), (sampling moment 2, sampling moment 4), (sampling moment 3, sampling moment 4).

[0101] For any sampling point, in each moment combination, the abnormal degree value of the building density at this sampling point is fused with the difference in the absorption rate of the wavelength to be measured at this sampling point at two sampling moments:

[0102] The value obtained by performing negative correlation mapping and normalization on the abnormal degree value of the building density at each sampling point is used as the retention weight; for any sampling point, the difference in the absorption rate of the wavelength to be measured at this sampling point at two sampling moments is used as the absorption rate change factor; the product of the absorption rate change factor at this sampling point in each moment combination and the corresponding retention weight is used as the absorption rate change characteristic value of this sampling point. Taking the wavelength as an example, the formula model of the absorption rate change characteristic value can specifically be, for example:

[0103]

[0104] where represents the absorption rate change characteristic value of the wavelength to be measured at the th moment combination for the sampling point ; represents the abnormal degree value of the building density for the sampling point ; represents at the th moment combination, the absorption rate of the wavelength to be measured at the first sampling moment for the sampling point ; represents at the th moment combination, the absorption rate of the wavelength to be measured at the second sampling moment for the sampling point ; represents the exponential function with the natural constant as the base.

[0105] In the formula model of the absorption rate change eigenvalue, since the larger the abnormal degree value of the building housing density, the lower the contribution degree of the spectral data at the sampling point to the analysis of the building housing density. Therefore, the abnormal degree value of the building housing density at the sampling point is subjected to negative correlation mapping and normalization processing to achieve logical relationship correction, and used as the retention weight. ; At the same time, the absorption rate change factor characterizes the difference in the absorption rate of the wavelength to be measured at the same sampling point at different times. Then, the absorption rate change factor is combined with the retention weight of the sampling point, and the product of the two is used as the absorption rate change eigenvalue of the sampling point at the wavelength to be measured. The absorption rate change eigenvalue combines the contribution degree of the spectral data of the sampling point to the analysis of the building housing density and the change of the absorption rate of the wavelength to be measured at the sampling point at different times. Therefore, analyzing the absorption rate change eigenvalue in the subsequent process can better evaluate the retention situation of the wavelength to be measured.

[0106] Finally, the standard deviation of the absorption rate change eigenvalues of all sampling points under the same time combination is used as the retention factor of the wavelength to be measured under the same time combination; the mean value of the retention factors of the wavelength to be measured under all sampling time combinations is used as the retention degree value of the wavelength to be measured. Taking the wavelength as an example, the formula model of the retention degree value can be specifically, for example:

[0107]

[0108] Among them, represents the retention degree value of the wavelength to be measured ; represents the total number of time combinations; represents the sampling point at the th time combination, and the absorption rate change eigenvalue of the wavelength to be measured ; represents the retention factor of the wavelength to be measured at the th time combination; represents calculating the standard deviation.

[0109] Calculate the standard deviation of the absorption rate change eigenvalues of the wavelength to be measured for all sampling points under the same time combination to obtain the retention factor of the wavelength to be measured. When the retention factor is larger, it indicates that the absorption rate change fluctuations of the wavelength to be measured at different sampling points, that is, different positions, are larger, so it is more likely to be caused by changes in the state of the building housing itself; on the contrary, if the retention factor is smaller, it indicates that the absorption rate change fluctuations of the wavelength to be measured at different sampling points, that is, different positions, are smaller, indicating that the wavelength to be measured is more likely to be affected by environmental factors, so it has little reference value for studying the building housing density situation.

[0110] Based on the above process, retention degree values for each wavelength can be obtained.

[0111] Step S4: According to the hyperspectral data of all sampling points at all sampling times, calculate the initial covariance matrix corresponding to the spectral data; correct the element values in the initial covariance matrix according to the retention degree values of each wavelength to obtain an updated covariance matrix; perform principal component analysis on the updated covariance matrix to obtain a dimensionality reduction result.

[0112] The principal component analysis method is a commonly used dimensionality reduction technique, mainly used to discover the most important patterns and relationships in data. In the embodiments of the present invention, the principal component analysis method is used to reduce the dimensionality of all spectral data, so as to extract key features, and the dimensionality reduction result is used to reflect the density status of the building houses. And constructing a covariance matrix is a necessary step and process in the principal component analysis process, which describes the linear relationship between variables. Therefore, the initial covariance matrix can be calculated first according to the hyperspectral data.

[0113] Preferably, in an embodiment of the present invention, the method for obtaining the initial covariance matrix includes:

[0114] Since the hyperspectral data includes two position data and one spectral data; and in the embodiments of the present invention, the spectral data is mainly analyzed for dimensionality reduction, so each spectral data is used as a column vector, and the length of the column vector is equal to the number of wavelengths, and the element value is the wavelength absorption rate.

[0115] Then, the column vectors corresponding to the spectral data of all sampling points at all sampling times are used to obtain an initial matrix, and the initial matrix is calculated to obtain an initial covariance matrix, and the number of rows and columns of the initial covariance matrix are both equal to the number of wavelengths. It should be noted that the process of calculating the covariance matrix in the principal component analysis method is well known to those skilled in the art and will not be elaborated here.

[0116] The initial covariance matrix can provide an overall view of the relationships between the wavelengths in the dataset. Then, since the retention degree values of each wavelength are calculated in step S3, which reflects the importance degree of each wavelength, the initial covariance matrix can be corrected by the retention degree values of the wavelengths to obtain an updated covariance matrix, so as to emphasize the influence of important wavelengths and reduce the influence of unimportant wavelengths at the same time, which is beneficial to ensuring that the subsequent principal component analysis pays more attention to those wavelengths that are valuable for analyzing the density status of the building houses.

[0117] Preferably, in an embodiment of the present invention, the method for obtaining the updated covariance matrix includes:

[0118] Construct a correction matrix based on the retention degree values of each wavelength. The correction matrix is a diagonal matrix, and the value of each element on the diagonal is the retention degree value of each wavelength, and the values of the elements in the remaining positions are 0.

[0119] Multiply the matrix obtained by multiplying the correction matrix by the initial covariance matrix by the correction matrix again to obtain the updated covariance matrix. The formula model of the updated covariance matrix is:

[0120]

[0121] Where, represents the updated covariance matrix; represents the correction matrix; represents the initial covariance matrix.

[0122] In the formula model of the updated covariance matrix, multiplying the initial covariance matrix on the left by a correction matrix is equivalent to performing a row transformation on the initial covariance matrix, that is, scaling the element values of each row in the initial covariance matrix, and multiplying by a correction matrix on the right is equivalent to performing a column transformation, that is, scaling the element values of each column. And the element values on the diagonal of the correction matrix are the retention degree values of each wavelength. Therefore, after the above processing, the importance weights of each wavelength can be added to the original dimensions of each wavelength, so that the subsequent dimensionality reduction results can better represent the density information of the building houses and reduce the influence degree of unimportant wavelengths.

[0123] It should be noted that the number of rows and columns of the correction matrix is equal to the total number of wavelengths. Taking the total number of wavelengths as 3 as an example, the correction matrix is a 3×3 matrix, and the element value in the first row and the first column is the retention degree value of the first wavelength, the element value in the second row and the second column is the retention degree value of the second wavelength, and the element value in the third row and the third column is the retention degree value of the third wavelength, and the values of the elements in the remaining positions are 0.

[0124] After obtaining the updated covariance matrix, the relationship between each wavelength in the updated covariance matrix has been adjusted at this time. Therefore, the subsequent principal component analysis process can be continued for the updated covariance matrix to obtain the dimensionality reduction result. Since the influence of unimportant wavelengths is eliminated in the obtained dimensionality reduction result, it can more accurately reflect the density information of the building houses and can provide a more accurate basis for urban management.

[0125] It should be noted that for the convenience of calculation, all the index data participating in the calculation in the embodiments of the present invention have undergone data preprocessing to eliminate the influence of the dimension. The specific means of eliminating the dimension influence are well-known technical means to those skilled in the art and will not be limited here; the principal component analysis method used in the embodiments of the present invention is a well-known technical means to those skilled in the art, and the specific process will not be elaborated here.

[0126] In summary, the purpose of the embodiments of the present invention is to improve the importance degree of the wavelength reflecting the building houses, so as to improve the accuracy of the dimensionality reduction result representing the density status of the building houses in the subsequent dimensionality reduction process. First, hyperspectral data of the area to be measured at different sampling times is obtained. Since even at the same sampling time, if only the difference in wavelength absorption rates between sampling points is analyzed, it is impossible to better distinguish whether the difference is affected by the density of vegetation. Therefore, when analyzing the difference in the absorption rates of wavelengths between different sampling points and the fluctuation of the absorption rates of all sampling points at the same wavelength, the difference in the density of vegetation between different sampling points is combined to obtain the density difference parameter between different sampling points at each sampling time. This parameter more accurately reflects the density difference between different sampling points considering the influence of the density of vegetation. Further, in order to determine whether each wavelength determines the density of plants, it is necessary to analyze the contribution degree of the spectral data of different sampling points to the plant density. Since the difference in the spectral data of sampling points with similar positions may be caused by the different density of plants, based on the aforementioned density difference parameter, combined with the positional relationship between sampling points, the abnormal degree value of the building house density at each sampling point is obtained, which is used to quantify the contribution degree of the spectral data at this sampling point to the building house density. Further, by fusing and analyzing the difference in the absorption rates of each wavelength of different sampling points at different times and the abnormal degree value of the building house density, the fluctuation situation can be evaluated, and the retention degree value of each wavelength can be obtained, which is used to distinguish the wavelengths with high research value and low research value for the density status of the building houses. Then, based on all the hyperspectral data, the initial covariance matrix is calculated, and the element values in the initial covariance matrix are corrected based on the retention degree value of each wavelength, so that the wavelengths with high retention degree values have higher weights, and the updated covariance matrix is obtained. Finally, the updated covariance matrix is continued to be subjected to principal component analysis, and the obtained dimensionality reduction result will contain more key information and can more accurately represent the density status of the building houses.

[0127] This embodiment also provides a building house remote sensing mapping data processing system, which includes a processor 500, a memory 501, and a computer program. The memory 501 is used to store the corresponding computer program, and the processor 500 is used to run the corresponding computer program. When the computer program runs on the processor 500, it can implement the steps of any one of the building house remote sensing mapping data processing methods.

[0128] Please refer to Figure 4, which shows a system block diagram of a remote sensing mapping data processing system for building houses provided by an embodiment of the present invention, including: a data acquisition module 401: used to implement step S1 in the above method, a density difference analysis module 402: used to implement step S2 in the above method, a wavelength retention degree analysis module 403: used to implement step S3 in the above method, and a dimensionality reduction module 404: used to implement step S4 in the above method.

[0129] Please refer to Figure 5 , which shows a schematic diagram of the system provided by an embodiment of the present invention, including a processor 500, a memory 501, a bus 502, and a communication interface 503. The processor 500, the communication interface 503, and the memory 501 are connected through the bus 502; among them, the memory 501 may include a high-speed random access memory, and the bus 502 may be an ISA bus, a PCI bus, an EISA bus, etc. The processor 500 may be an integrated circuit chip with signal processing capabilities.

[0130] An embodiment of the present invention also provides a computer-readable storage medium corresponding to the method provided in the foregoing embodiment. Specifically, it may be, for example, an optical disc, on which a computer program (i.e., a program product) is stored. When the computer program is run by a processor, it will execute the method provided in any of the foregoing embodiments.

[0131] It should be noted that examples of the computer-readable storage medium may also include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), and other optical and magnetic storage media, which will not be elaborated here one by one.

[0132] It should be noted that the above sequence of embodiments of the present invention is only for description and does not represent the superiority or inferiority of the embodiments. The processes depicted in the drawings do not necessarily require the specific order or continuous order shown to achieve the desired result. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0133] Each embodiment in this specification is described in a progressive manner. The same or similar parts among the embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments.

Claims

1. A method for processing remote sensing mapping data of a building house, characterized in that, The method includes: Obtaining hyperspectral data of the area to be measured at different sampling times; For the hyperspectral data at any sampling time, according to the differences in vegetation density between different sampling points, the differences in absorption rates of wavelengths at different sampling points, and the fluctuation of absorption rates of all sampling points at the same wavelength, obtain the density difference parameter between different sampling points in the hyperspectral data at this sampling time; For any sampling point in the area to be measured, fuse and analyze the positional relationship, spectral data differences, and density difference parameters between this sampling point and the remaining other sampling points at all sampling times to obtain the abnormal degree value of building house density at this sampling point; for any wavelength, fuse and analyze the fluctuation of the differences in absorption rates of this wavelength at all sampling points at different times and the abnormal degree value of building house density of each sampling point to obtain the retention degree value of this wavelength; Calculate the initial covariance matrix corresponding to the spectral data according to the hyperspectral data of all sampling points at all sampling times; correct the element values in the initial covariance matrix according to the retention degree values of each wavelength to obtain an updated covariance matrix; perform principal component analysis on the updated covariance matrix to obtain a dimensionality reduction result; The method for obtaining the abnormal degree value of building house density includes: Optionally select a sampling point as the target point, and screen out the comparison points of the target point from the remaining sampling points according to the distances between the target point and each of the remaining sampling points; For any sampling time, take the difference in absorption rates of the same wavelength in the spectral data between the target point and each comparison point as the absorption rate distance value; take the sum value of the absorption rate distance values of all wavelengths as the distance sum value, and according to the density difference parameter between the target point and each comparison point and the distance sum value, obtain the deviation value of the spectral data between the target point and each comparison point at this sampling time, where the deviation value is negatively correlated with the density difference parameter, and the distance sum value and the deviation value are positively correlated; Take the mean value of the deviation values of the spectral data between the target point and each comparison point at all sampling times as the abnormal factor of building house density of the target point; take the mean value of the abnormal factors of building house density between the target point and all comparison points as the abnormal degree value of building house density at the target point.

2. The method for processing remote sensing mapping data of a building house according to claim 1, characterized in that, The method for obtaining the density difference parameter includes: Obtain the normalized difference vegetation index at each sampling point in the area to be measured; At any sampling time, take the difference in the normalized difference vegetation index between two different sampling points as the vegetation condition difference value; Analyze the difference in absorption rates of the same wavelength between two different sampling points and the fluctuation of absorption rates of all sampling points at the same wavelength to obtain the density difference factor between different sampling points at this sampling time; According to the vegetation condition difference value and the density difference factor between two different sampling points, obtain the density difference parameter between two different sampling points at this sampling time, where the density difference parameter is positively correlated with the density difference factor, the density difference parameter is negatively correlated with the vegetation condition difference value, and the value of the density difference parameter is a normalized value.

3. A method for processing remote sensing mapping data of a building house according to claim 2, characterized in that, The method for obtaining the density difference factor includes: Optionally select a wavelength as the wavelength to be measured. At this sampling moment, take the difference in absorbance at the wavelength to be measured between two different sampling points as the absorbance difference value; take the standard deviation of the absorbances at the wavelength to be measured of all sampling points as the absorbance fluctuation value; Obtain the absorbance change index at the wavelength to be measured for two sampling points based on the absorbance fluctuation value and the absorbance difference value. Both the absorbance fluctuation value and the absorbance difference value vary in the same direction as the absorbance change index; take the average value of the absorbance change indexes at all wavelengths for two different sampling points as the density difference factor.

4. A method for processing remote sensing mapping data of a building house according to claim 1, characterized in that, The method for obtaining the comparison points includes: Calculate the Euclidean distance between the target point and the positions of each of the remaining sampling points as the distance factor, and take the sampling points with the distance factor less than the preset distance threshold as the comparison points of the target point.

5. A method for processing remote sensing mapping data of a building house according to claim 3, characterized in that, The method for obtaining the retention degree value includes: Combine all sampling moments in pairs to obtain all non-repeating moment combinations; For any one sampling point, in each moment combination, fuse the abnormal degree value of the building house density at this sampling point with the absorbance difference at the wavelength to be measured at two sampling moments of this sampling point to obtain the absorbance change characteristic value of this sampling point; Take the standard deviation of the absorbance change characteristic values at the wavelength to be measured in the same moment combination for all sampling points as the retention factor at the wavelength to be measured in the same moment combination; take the average value of the retention factors at the wavelength to be measured in all moment combinations as the retention degree value of the wavelength to be measured.

6. A method for processing remote sensing mapping data of a building house according to claim 5, characterized in that, The method for obtaining the absorbance change characteristic value includes: Take the value obtained by performing negative correlation mapping and normalization on the abnormal degree value of the building house density at each sampling point as the retention weight; for any one sampling point, in each moment combination, take the absorbance difference at the wavelength to be measured at two sampling moments of this sampling point as the absorbance change factor; take the product of the absorbance change factor at this sampling point in each moment combination and the corresponding retention weight as the absorbance change characteristic value of this sampling point at the wavelength to be measured.

7. A method for processing remote sensing mapping data of a building house according to claim 1, characterized in that, The method for obtaining the initial covariance matrix includes: The hyperspectral data includes two position data and one spectral data; Take each spectral data as a column vector, the length of the column vector is equal to the number of wavelengths, and the element value is the wavelength absorbance; Form an initial matrix with the column vectors corresponding to the spectral data of all sampling points at all sampling moments, and perform calculations on the initial matrix to obtain the initial covariance matrix.

8. A method for processing remote sensing mapping data of a building house according to claim 1, characterized in that, The method for obtaining the updated covariance matrix includes: Construct a correction matrix based on the retention degree values of each wavelength. The correction matrix is a diagonal matrix, and each element value on the diagonal is the retention degree value of each wavelength, and the element values at the remaining positions are 0; Multiply the matrix obtained by multiplying the correction matrix by the initial covariance matrix by the correction matrix again to obtain the updated covariance matrix.

9. A remote sensing mapping data processing system for building houses, comprising a memory, a processor, and a computer program stored in the memory and capable of running on the processor, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 8.

Images