CN101871884B - Atmospheric correction and regional mineral map spotting method utilizing multi-scene ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) remote sensing data - Google Patents
Atmospheric correction and regional mineral map spotting method utilizing multi-scene ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) remote sensing data Download PDFInfo
- Publication number
- CN101871884B CN101871884B CN2010101968685A CN201010196868A CN101871884B CN 101871884 B CN101871884 B CN 101871884B CN 2010101968685 A CN2010101968685 A CN 2010101968685A CN 201010196868 A CN201010196868 A CN 201010196868A CN 101871884 B CN101871884 B CN 101871884B
- Authority
- CN
- China
- Prior art keywords
- data
- atmospheric
- atmospheric correction
- aster
- reflectivity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention relates to an atmospheric correction and regional mineral map spotting method utilizing multi-scene ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) remote sensing data, comprising the following nine steps of: (1) reading in ASTER data; (2) selecting reference data, wherein one-scene data is selected from multi-scene data as a reference for atmospheric correction; (3) carrying out the atmospheric correction on the reference data to obtain reference reflectivity data; (4) carrying out atmospheric correction on non-reference data to obtain non-reference reflectivity data; (5) calculating the overlay region spectral characteristic difference of the reference reflectivity data and the non-reference reflectivity data; (6) calculating the atmospheric water vapor content and the atmospheric visibility of the non-reference data; (7) newly carrying out the atmospheric correction on the non-reference data by utilizing the calculated atmospheric water vapor content and the atmospheric visibility to obtain non-reference reflectivity data; (8) carrying out inlaying treatment on the multi-scene reflectivity data including the reference data and the non-reference data; and (9) carrying out mineral identification and map spotting by utilizing spectral indexes. The method effectively eliminates inlaid gaps among the multi-scene data, thereby realizing the regional mineral map spotting.
Description
Technical field
The invention belongs to the geologic remote sensing technology, the method for especially a kind of many scapes ASTER remotely-sensed data atmospheric correction and regional mineral map plotting.
Background technology
The ASTER data are ETM/ETM
+The earth observation satellite remotely-sensed data of a new generation after the data.These data can be discerned the altered mineral relevant with metalliferous deposit mineralization, and such as sericite, smalite, chlorite, epidote etc., application potential is huge in resource exploration.Yet, this data processing complex, traditional ETM/ETM
+Processing method of data has been difficult to meet the demands.Data processing mainly is made up of two big steps: 1) atmospheric correction; 2) mineral identification and charting.
Official's algorithm by the ASTER data atmospheric correction of ASTER sensor design team issue is; At first utilize the atmospheric parameters such as MODIS data extract atmospheric water body burden, aerosol optical depth of phase simultaneously; Utilize the atmospheric parameter that extracts that data are carried out atmospheric correction afterwards, eliminate the influence of atmosphere spectrum.Utilize the complex algorithm of MODIS data inversion atmosphere profile parameter in this method, and need obtain the MODIS data of phase simultaneously, for domestic consumer, the difficulty of realization is very big.In addition, when utilizing many scapes data to carry out regional mineral map plotting, the error of atmospheric parameter inverting can cause many scapes data to be inlayed and the difficulty of splicing.
The ASTER data are compared ETM/ETM
+Data have more wave band, when mineral identification and charting, and traditional ETM/ETM
+The method of data two waveband computing can't be given full play to the advantage of ASTER data, needs exploitation to be applicable to the mineral identification spectrum index of ASTER data characteristics.
Summary of the invention
The objective of the invention is, inlay the difficulty with splicing, and the two waveband operational method can't give full play to the ASTER data edge, compare ETM/ETM to the complicacy of ASTER data official atmospheric correction algorithm and many scapes data
+Data; Problem with more wave band; A kind of many scapes ASTER data atmospheric correction and the regional mineral map plotting method that propose; At the atmospheric correction that need not only can accomplish under the situation of atmospheric parameter the ASTER data, can also effectively eliminate the mineral map plotting results' that cause because of the atmospheric parameter difference between many scapes data difference.
Technical solution of the present invention is to comprise following nine steps:
(1) the ASTER data reads in;
(2) reference data is selected, from many scapes data, chooses the benchmark of scape data as atmospheric correction;
(3) atmospheric correction of reference data obtains the benchmark reflectivity data;
(4) atmospheric correction of non-reference data obtains non-benchmark reflectivity data;
(5) calculating benchmark reflectivity data and non-benchmark reflectivity data overlay region spectral signature are poor;
(6) calculate non-reference data atmosphere vapour content and atmospheric visibility;
(7) utilize the atmosphere vapour content and the atmospheric visibility that calculate that non-reference data is carried out atmospheric correction again, obtain non-benchmark reflectivity data;
Many scapes reflectivity data that (8) will comprise benchmark reflectivity data and non-benchmark reflectivity data carries out damascene;
(9) utilize spectrum index to carry out mineral identification and charting, obtain the mineral distribution plan.
At first carry out reading in of step (1) described ASTER data;
Selecting of the described reference data of step (2), it is chosen to be,
The data of choosing near the regional center of charting are reference data, and other data comprise non-reference data, all are that benchmark carries out atmospheric correction with the reference data.
The method of the atmospheric correction of the non-reference data described in the atmospheric correction of the reference data described in the step (3) and the step (4) is:
To big pneumatic jack 0-100km, atmosphere is divided into 36 layers from the face of land, the condition of each atmosphere comprises physics, chemical parameters such as temperature, pressure, various atmospheric molecule concentration.Atmospheric correction method adopts 8 stream radiation delivery models.Atmosphere profile parameter is from the change curve of 0-36 atmosphere parameter, by the normal atmosphere mode-definition.
After atmospheric correction is accomplished; Utilize some object spectrums in the image; Whether atmosphere vapour content that defines in blank ground, vegetation, the water body criterion atmospherical model and aerosol optical depth and true atmospheric parameter be approaching; Difference according to the two resets atmosphere vapour content and aerosol optical depth in the normal atmosphere pattern, and utilizes 8 radiation delivery models that data are carried out atmospheric correction again.
Step (5) described calculating benchmark reflectivity data and non-benchmark reflectivity data overlay region spectral signature are poor, and its method is:
At first, calculate the absorption degree of depth of benchmark reflectivity data and non-benchmark reflectivity data overlay region wave band 8 respectively, wave band 7 and the ratio of wave band 9 reflectivity sums, and the ratio of wave band 4 and wave band 3 with wave band 8 reflectivity; Then, calculate benchmark reflectivity data and non-benchmark reflectivity data overlay region difference.
Step (6) described calculating non-reference data atmosphere vapour content and atmospheric visibility, its method is:
According to statistical relationship and the ratio of wave band 4 and 3 and the statistical relationship of atmospheric visibility of the wave band 8 absorption degree of depth and atmosphere vapour content, calculate the atmosphere vapour content and the atmospheric visibility of non-reference data.
Atmosphere vapour content and atmospheric visibility that the described utilization of step (7) is calculated carry out atmospheric correction again to non-reference data, obtain non-benchmark reflectivity data, and its method is:
With atmosphere vapour content and the atmospheric visibility data input normal atmosphere pattern that step (6) is calculated, utilize 8 stream radiation delivery models that non-reference data is carried out atmospheric correction again, obtain non-benchmark reflectivity data.
The described many scapes reflectivity data of benchmark reflectivity data and non-benchmark reflectivity data that will comprise of step (8) carries out damascene, and its method is:
In order to obtain the unified mosaic map mosaic of tone, at first carry out the gray scale coupling of each band image.Partly have the gray-scale value of each image of principle adjustment of identical average gray and variance according to the doubling of the image, can make mosaic map mosaic not have tangible splicing seams like this.
The described spectrum index of utilizing of step (9) is carried out mineral identification and charting, obtains the mineral distribution plan, and its method is:
Based on the mineral identification of spectrum index and the basis of charting is the foundation of mineral identification spectrum index.Spectrum according to the common altered mineral of metalliferous deposit.Set up the spectrum index table of mineral identifications.Enforcement utilizes spectrum index to carry out mineral identification and charting, obtains the mineral distribution plan.
Advantage of the present invention is at the atmospheric correction that need not only can accomplish under the situation of atmospheric parameter the ASTER data, can also effectively eliminate the mineral map plotting results' that cause because of the atmospheric parameter difference between many scapes data difference.
Description of drawings
Fig. 1 is a process flow diagram of the present invention;
Fig. 2 is the influence figure of ASTER data atmosphere vapour content to the atmospheric correction result;
Wherein atmospheric visibility is 40km, and reflectance spectrum is corresponding: 1. atmosphere vapour content 1.4637g/cm
2, 2. atmosphere vapour content 1.171g/cm
2, 3. atmosphere vapour content 0.8182g/cm
2, 4. atmosphere vapour content 0.5855g/cm
2, 5. atmosphere vapour content 0.2927g/cm
2
Fig. 3 is the influence figure of ASTER data atmospheric visibility to the atmospheric correction result;
Wherein atmosphere vapour content is 0.9514g/cm
2, the atmosphere that reflectance spectrum is corresponding: 6. visibility 60km, 7. visibility 40km, 8. visibility is respectively 20km;
Fig. 4 is the graph of a relation that ASTER data wave band 8 absorbs the degree of depth and atmosphere vapour content;
Wherein atmospheric visibility is 40km, and fitting formula is y=0.059x+0.9496, and related coefficient is 0.9992;
Fig. 5 is the graph of a relation of ASTER data wave band 3,4 ratios and visibility;
Wherein the gasoloid pattern is rural pattern;
Fig. 6 is the spectrogram of the common altered mineral of metalliferous deposit;
Mineral from top to bottom are respectively alunite, white (thin,tough silk) mica, jarosite, smalite, smectite, rauhkalk, epidote, chlorite, kalzit; Wherein dotted line is that spectral resolution is the spectrum of 10 nanometers, and solid line is the spectrum that is resampled to behind the ASTER data wave band.
Embodiment
For many scapes ASTER remotely-sensed data atmospheric correction and the regional mineral map plotting method that the present invention relates to better is described, utilize East Tianshan of Xinjiang area three scape ASTER to carry out mineral map plotting, realization flow is as shown in Figure 1, and concrete steps are following:
(1) high-spectral data reads in
Read in East Tianshan of Xinjiang area three scape ASTER, former size of data is for being respectively 2827 (OK) * 2513 (row) * 9 (wave band) (first scape), 2797 (OK) * 2475 (row) * 9 (wave band) (second scape), 2826 (OK) * 2512 (row) * 9 (wave band) (the 3rd scape).
(2) reference data is selected
First scape is chosen to be reference data.
(3) atmospheric correction of reference data
To big pneumatic jack 0-100km, atmosphere is divided into 36 layers from the face of land, the condition of each atmosphere comprises physics, chemical parameters such as temperature, pressure, various atmospheric molecule concentration.Atmospheric correction method adopts 8 stream radiation delivery models.Atmosphere profile parameter is from the change curve of 0-36 atmosphere parameter, by the normal atmosphere mode-definition.
After atmospheric correction is accomplished; Utilize some object spectrums in the image; Blank ground, vegetation, water body; Whether atmosphere vapour content that defines in the criterion atmospherical model and atmospheric visibility and true atmospheric parameter be approaching, resets atmosphere vapour content and atmospheric visibility in the normal atmosphere pattern according to the two difference, and utilize the radiation delivery model that data are carried out atmospheric correction again.
Analyze the influence rule of different atmospheric parameters to face of land reflectivity inverting, be the atmosphere vapour content that defines in the criterion atmospherical model and atmospheric visibility and true atmospheric parameter whether near and reset atmosphere vapour content and the key of atmospheric visibility in the normal atmosphere pattern.Atmosphere vapour content mainly impacts ASTER data 3,5,6,8,9 wave band reflectivity, wherein 8,9 wave bands more so (accompanying drawing 2).Accompanying drawing 2 is the test site, even blank ground reflectivity inversion result when different atmosphere vapour content in the East Tianshan of Xinjiang area, and the reducing of the atmosphere vapour content of importing during with atmospheric correction, wave band 3,5,6,8,9 reflectivity are reduction in various degree.Can know that by the ground measured spectra blank ground spectrum does not have tangible absorption characteristic at wave band 8 places, with atmosphere vapour content in the accompanying drawing 2 be 0.8182g/cm
2Reflectivity when (the atmosphere vapour pressure on every sq) is approaching, can think that thus the atmosphere vapour content when this scape ASTER data are obtained is 0.8182g/cm
2The atmosphere vapour content value just can be imported 0.8182g/cm during atmospheric correction
2Atmospheric visibility mainly influences the reflectivity (accompanying drawing 3) of ASTER the 3rd, 4 wave bands, and the reflectance varies of wave band 3,4 can not influence follow-up mineral identification and charting, as long as farthest eliminate the difference in reflectivity between many scapes data.For the reference data of ceiling unlimited mist, atmospheric visibility generally is made as 40-60km.After having confirmed the atmosphere vapour content and atmospheric visibility when reference data is obtained, 8 stream radiation delivery models promptly capable of using are realized atmospheric correction.
Atmosphere vapour molecular conecentration value is input as 0.8182g/cm
2, the atmospheric visibility value is input as 40km.
(4) atmospheric correction of non-reference data
To big pneumatic jack 0-100km, atmosphere is divided into 36 layers from the face of land, the condition of each atmosphere comprises physics, chemical parameters such as temperature, pressure, various atmospheric molecule concentration.Atmospheric correction method adopts 8 stream radiation delivery models.Atmosphere profile parameter, the change curve of from 0 to 36 atmosphere parameter is by the normal atmosphere mode-definition.
(5) calculating reference data and non-reference data overlay region spectral signature are poor
At first, calculate the absorption degree of depth of benchmark reflectivity data and non-benchmark reflectivity data overlay region wave band 8 respectively, wave band 7 and the ratio of wave band 9 reflectivity sums, and the ratio of wave band 4 and wave band 3 with wave band 8 reflectivity; Then, calculate benchmark reflectivity data and non-benchmark reflectivity data overlay region difference.
The difference that the difference that second scape, the 3rd scape and the first scape overlay region wave band 8 absorb the degree of depth is respectively 0.01,0.03, the second scape, the 3rd scape and the first scape overlay region wave band 4 and wave band 3 ratios is respectively 0.8565,0.8976.
(6) calculate non-reference data atmosphere vapour content and atmospheric visibility
According to relation (accompanying drawing 4) and the ratio of wave band 4 and 3 and the statistical relationship (accompanying drawing 5) of atmospheric visibility of the wave band 8 absorption degree of depth and atmosphere vapour content, calculate the atmosphere vapour content and the atmospheric visibility of non-reference data.
The atmosphere vapour content of second scape, the 3rd scape data is respectively 1.0g/cm
2, 1.38g/cm
2The atmospheric visibility of second scape, the 3rd scape data is respectively 43km, 57km.
(7) utilize the atmosphere vapour content and the atmospheric visibility that calculate that non-reference data is carried out atmospheric correction again, obtain non-benchmark reflectivity data, its method is:
With atmosphere vapour content and the atmospheric visibility data input normal atmosphere pattern that step (6) is calculated, utilize 8 stream radiation delivery models that non-reference data is carried out atmospheric correction again, obtain non-benchmark reflectivity data.
Many scapes reflectivity data that (8) will comprise benchmark reflectivity data and non-benchmark reflectivity data carries out damascene, and its method is:
In order to obtain the unified mosaic map mosaic of tone, at first carry out the gray scale coupling of each band image.Partly have the gray-scale value of each image of principle adjustment of identical average gray and variance according to the doubling of the image, can make mosaic map mosaic not have tangible splicing seams like this.
(9) utilize spectrum index to carry out mineral identification and charting, obtain the mineral distribution plan, its method is:
Based on the mineral identification of spectrum index and the basis of charting is the foundation of mineral identification spectrum index.According to the spectrum (accompanying drawing 6) of the common altered mineral of metalliferous deposit, the spectrum index table of these mineral identifications below having set up.
The spectrum index table of common mineral identification
Explain: numeral is all represented the reflectivity of corresponding wave band in the table, like the ratio of spectrum index 2/1 expression wave band 2 with wave band 1 reflectivity; The U presentation logic " with ".
Ferriferous oxide refers to oxides-containing irons such as limonite, haematite, and ASTER data wave band 2,1 and ETM+ wave band 3,1 position are approaching, and wave band 2 is identical with the method that the ratio that utilizes ETM+ wave band 3 with wave band 1 reflectivity extracts ferriferous oxide with the ratio of wave band 1 reflectivity.Silicate iron refers to compose the iron that is stored in the lattice silicate; Be mainly used in the higher altered minerals of iron content such as dark-coloured rock forming minerals such as extracting pyroxene, hornblende and epidote, chlorite; Can know by accompanying drawing 6; The ratio of chlorite, epidote wave band 5 and wave band 4 reflectivity is all greater than 1, and other mineral are then less than, the spectrum index of therefore ratio of wave band 5 and 4 being discerned as silicate iron.The spectrum index of the identification of jarosite does, if wave band 6,8 reflectivity all greater than 7, promptly available wave band 7 is discerned with respect to the absorption degree of depth ((6+8) * 0.5/7) of wave band 6,8, it is big more to absorb the degree of depth, jarosite content is high more.Kaolinic identification is mainly according to the weak absorption characteristic at wave band 5,6 places; The spectrum index of its identification does; If (If of table 1) wave band 5,7 reflectivity are all greater than wave band 6 reflectivity, there is poor absorption in wave band 5 with respect to wave band 4,6 simultaneously, extracts smalite with regard to available band 6 with respect to the absorption degree of depth ((5+7) * 0.5/6) of wave band 5,7; It is big more to absorb the degree of depth; Smalite content is high more, and wave band 5 is discerned with respect to the poor absorption available spectrum index (1.082192*4-0.082192*6)＞5 of wave band 4,6, (1.082192*4-0.082192*6) for utilizing the reflectivity of inserting the wave band 5 that obtains in the wave band 4,6.White mica is very similar in ASTER wave band polishing wax characteristic with smectite; Be difficult to distinguish identification; Identification spectrum index be similar to smalite, unique difference be wave band 5 with respect to wave band 4,6, do not have the absorption characteristic; I.e. (1.082192*4-0.082192*6)＜5, the same smalite of physical significance (1.082192*4-0.082192*6).The spectrum index of alunite identification does, if wave band 4,6 reflectivity all greater than 5, promptly available wave band 5 is discerned with respect to the absorption degree of depth ((4+6) * 0.5/5) of wave band 4,6, it is big more to absorb the degree of depth, alunite content is high more.Carbonate, chlorite, epidote are very similar in ASTER wave band polishing wax characteristic; Be difficult to distinguish; The spectrum index of its identification does, if wave band 7,9 reflectivity all greater than 8, promptly available wave band 5 is discerned with respect to the absorption degree of depth ((7+9) * 0.5/8) of wave band 4,6; It is big more to absorb the degree of depth, and its content is high more.
According to the spectral signature of altered mineral, set the spectrum index of identification and charting, accomplish mineral identification and charting according to spectrum index.
Claims (8)
1. remotely-sensed data atmospheric correction of scape ASTER more than a kind and regional mineral map plotting method, it is characterized in that: it comprises following nine steps:
(1) the ASTER data are read in;
(2) reference data is selected, from many scapes data, chooses the benchmark of scape data as atmospheric correction;
(3) atmospheric correction of reference data obtains the benchmark reflectivity data;
(4) atmospheric correction of non-reference data obtains non-benchmark reflectivity data;
(5) calculating benchmark reflectivity data and non-benchmark reflectivity data overlay region spectral signature are poor;
(6) calculate non-reference data atmosphere vapour content and atmospheric visibility;
(7) utilize the atmosphere vapour content and the atmospheric visibility that calculate that non-reference data is carried out atmospheric correction again, obtain non-benchmark reflectivity data;
(8) the many scapes reflectivity data that comprises reference data and non-reference data is carried out damascene;
(9) utilize spectrum index to carry out mineral identification and charting, obtain the mineral distribution plan.
2. by described a kind of many scapes ASTER remotely-sensed data atmospheric correction of claim 1 and regional mineral map plotting method, it is characterized in that: the selected method of the reference data described in the step (2) is: selected data near the charting regional center are reference data.
3. by described a kind of many scapes ASTER remotely-sensed data atmospheric correction of claim 1 and regional mineral map plotting method, it is characterized in that: the method for the atmospheric correction of the non-reference data described in the atmospheric correction of the reference data described in the step (3) and the step (4) is:
The extremely big pneumatic jack 0-100km from the face of land; Atmosphere is divided into 36 layers; The condition of each atmosphere comprises temperature, pressure, various atmospheric molecule concentration parameter, and atmospheric correction method adopts 8 stream radiation delivery models, atmosphere profile parameter; From the change curve of 0-36 atmosphere parameter, by the normal atmosphere mode-definition;
After atmospheric correction is accomplished; Utilize some atural objects in the image; Whether atmosphere vapour content that defines in the standard spectrum criterion atmospherical model of blank ground, vegetation, water body and aerosol optical depth and true atmospheric parameter be approaching; Difference according to the two resets atmosphere vapour content and aerosol optical depth in the normal atmosphere pattern, and utilizes 8 radiation delivery models that data are carried out atmospheric correction again.
4. by described a kind of many scapes ASTER remotely-sensed data atmospheric correction of claim 1 and regional mineral map plotting method, it is characterized in that: the method for calculating benchmark reflectivity data described in the step (5) and non-benchmark reflectivity data overlay region spectral signature difference is:
At first, calculate the absorption degree of depth of benchmark reflectivity data and non-benchmark reflectivity data overlay region wave band 8 respectively, wave band 7 and the ratio of wave band 9 reflectivity sums, and the ratio of wave band 4 and wave band 3 with wave band 8 reflectivity; Then, calculate benchmark reflectivity data and non-benchmark reflectivity data overlay region difference.
5. by described a kind of many scapes ASTER remotely-sensed data atmospheric correction of claim 1 and regional mineral map plotting method, it is characterized in that: the calculating non-reference data atmosphere vapour content described in the step (6) and the method for atmospheric visibility are:
Absorb the statistical relationship of the degree of depth and atmosphere vapour content according to wave band 8, and the statistical relationship of the ratio of wave band 4 and 3 and atmospheric visibility, calculate the atmosphere vapour content and the atmospheric visibility of non-reference data.
6. by described a kind of many scapes ASTER remotely-sensed data atmospheric correction of claim 1 and regional mineral map plotting method, it is characterized in that: the method that atmosphere vapour content that the utilization described in the step (7) is calculated and atmospheric visibility carry out atmospheric correction again to non-reference data is:
With atmosphere vapour content and the atmospheric visibility data input normal atmosphere pattern that step (6) is calculated, utilize 8 stream radiation delivery models that non-reference data is carried out atmospheric correction again, obtain non-benchmark reflectivity data.
7. by described a kind of many scapes ASTER remotely-sensed data atmospheric correction of claim 1 and regional mineral map plotting method, it is characterized in that: the method that the many scapes reflectivity data to comprising reference data and non-reference data described in the step (8) carries out damascene is:
At first carry out the gray scale coupling of each band image, adjust the gray-scale value of each image, make mosaic map mosaic not have tangible splicing seams.
8. by described a kind of many scapes ASTER remotely-sensed data atmospheric correction of claim 1 and regional mineral map plotting method, it is characterized in that: the method for utilizing spectrum index to carry out mineral identification and charting described in the step (9) is:
According to the spectral signature of altered mineral, set the spectrum index of identification and charting, accomplish mineral identification and charting according to spectrum index.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2010101968685A CN101871884B (en) | 2010-06-02 | 2010-06-02 | Atmospheric correction and regional mineral map spotting method utilizing multi-scene ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) remote sensing data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2010101968685A CN101871884B (en) | 2010-06-02 | 2010-06-02 | Atmospheric correction and regional mineral map spotting method utilizing multi-scene ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) remote sensing data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101871884A CN101871884A (en) | 2010-10-27 |
| CN101871884B true CN101871884B (en) | 2012-06-27 |
Family
ID=42996878
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2010101968685A Expired - Fee Related CN101871884B (en) | 2010-06-02 | 2010-06-02 | Atmospheric correction and regional mineral map spotting method utilizing multi-scene ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) remote sensing data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101871884B (en) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102032903B (en) * | 2010-12-01 | 2012-11-21 | 南京大学 | Coral island reef remote sensing information automatic extraction method based on Landsat data source |
| CN102426625B (en) * | 2011-09-21 | 2014-06-18 | 中国国土资源航空物探遥感中心 | Method for extracting porphyry copper typical altered zone mineral information based on ASTER satellite data |
| CN103267533B (en) * | 2013-05-08 | 2015-09-30 | 中国科学院光电研究院 | A kind of practical high-spectrum remote sensing air automatic correcting method |
| CN103984940B (en) * | 2014-06-03 | 2017-12-26 | 核工业北京地质研究院 | A kind of method based on high-spectral data identification hematization |
| CN104990880B (en) * | 2015-07-09 | 2017-10-27 | 核工业北京地质研究院 | Sericite mineral relative abundance computational methods based on Airborne Hyperspectral remotely-sensed data |
| CN105651710A (en) * | 2015-12-30 | 2016-06-08 | 淮北师范大学 | Optical remote sensing method for acquiring absolute light intensity of aerosol |
| CN109300133B (en) * | 2018-11-19 | 2020-10-23 | 珠江水利委员会珠江水利科学研究院 | Urban river network area water body extraction method |
| CN111044477A (en) * | 2019-12-20 | 2020-04-21 | 核工业北京地质研究院 | Remote sensing method for rapidly distinguishing silicification zone from acidic dike |
| CN113406041B (en) * | 2021-05-31 | 2023-03-17 | 核工业北京地质研究院 | Method for obtaining key altered mineral combination of sodium-substituted rock type uranium ore |
| CN114076955B (en) * | 2021-09-27 | 2022-07-19 | 农业农村部规划设计研究院 | Method for identifying lodging corn based on high-score first satellite data |
| CN114332589B (en) * | 2021-12-31 | 2023-07-11 | 中国科学院紫金山天文台 | Fine detection method for surface water or hydroxyl of atmospheric celestial body |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100458562C (en) * | 2004-04-28 | 2009-02-04 | 北京矿产地质研究院 | Method for making remote sensing shadowgraph of positive landform |
| CN1317569C (en) * | 2004-06-29 | 2007-05-23 | 中国国土资源航空物探遥感中心 | Layered lineage identification method for high spectrum mineral |
| CN1317551C (en) * | 2004-06-29 | 2007-05-23 | 中国国土资源航空物探遥感中心 | High spectrum minerals maximum correlation identification method based on spectrum hybrid composition |
| CN1595110A (en) * | 2004-06-29 | 2005-03-16 | 中国国土资源航空物探遥感中心 | High spectrum multiple end element linear fitting mineral identification and mineralizing indication method |
| CN101551471B (en) * | 2009-05-19 | 2012-01-18 | 中国国土资源航空物探遥感中心 | High-spectrum remote-sensing mineral content quantitative inversion method |
-
2010
- 2010-06-02 CN CN2010101968685A patent/CN101871884B/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN101871884A (en) | 2010-10-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101871884B (en) | Atmospheric correction and regional mineral map spotting method utilizing multi-scene ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) remote sensing data | |
| CN102324106B (en) | SFS (Shape From Shading) three-dimensional reconstruction sparse-DEM (Digital Elevation Model) encrypting method considering surface spectral information | |
| Degerickx et al. | Enhancing the performance of Multiple Endmember Spectral Mixture Analysis (MESMA) for urban land cover mapping using airborne lidar data and band selection | |
| Jones et al. | Assessing the utility of airborne hyperspectral and LiDAR data for species distribution mapping in the coastal Pacific Northwest, Canada | |
| Small et al. | Multi-scale standardized spectral mixture models | |
| Xin et al. | Using remote sensing data to quantify changes in glacial lakes in the Chinese Himalaya | |
| CN103984862B (en) | A kind of snow parameter inversion method of polynary remote sensing information collaboration | |
| CN104123559B (en) | Multi-source remote sensing discrimination method and system for underground potassium-containing brine resources in salt lake region | |
| EP3750104A1 (en) | A system and method to fuse multiple sources of optical data to generate a high-resolution, frequent and cloud-/gap-free surface reflectance product | |
| CN105651263B (en) | Shallow water depth multi-source remote sensing merges inversion method | |
| CN102938066A (en) | Method for reconstructing outer outline polygon of building based on multivariate data | |
| CN101916377B (en) | Method for intelligently identifying mineral information based on spectral information | |
| CN104808255A (en) | Fractal theory-based mineralization anomaly information mining method | |
| CN106600574A (en) | Landslide extraction method based on remote-sensing image and altitude data | |
| CN103136736B (en) | Hyperspectral remote sensing data non-linear dimension descending method | |
| CN104463164A (en) | Tree canopy structure information extraction method based on rib method and crown height ratio | |
| Kassner et al. | Analysis of the solar potential of roofs by using official lidar data | |
| CN103592235B (en) | A kind of spectral correction method of imaging spectral reflectivity data | |
| CN101673339A (en) | Target identification method of high spectroscopic data rearranged based on spectral absorption characteristics | |
| Fu et al. | Evaluation of LAI estimation of mangrove communities using DLR and ELR algorithms with UAV, hyperspectral, and SAR images | |
| CN103926203B (en) | A kind of for the probabilistic Spectral Angle Mapping method of object spectrum | |
| CN103413289B (en) | Image processing method and device for mining area | |
| Sun et al. | Identifying terraces in the hilly and gully regions of the Loess Plateau in China | |
| Oliveira et al. | The use of remote sensing tools for accurate charcoal kilns’ inventory and distribution analysis: Comparative assessment and prospective | |
| CN109839356A (en) | Alteration information extraction method based on WorldView-3 satellite data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| C17 | Cessation of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120627 Termination date: 20130602 |