CN1556381A - Aviation high spectrum remote sensing flight ground synchronous scaling and reflectivity conversion method - Google Patents
Aviation high spectrum remote sensing flight ground synchronous scaling and reflectivity conversion method Download PDFInfo
- Publication number
- CN1556381A CN1556381A CNA2003101229433A CN200310122943A CN1556381A CN 1556381 A CN1556381 A CN 1556381A CN A2003101229433 A CNA2003101229433 A CN A2003101229433A CN 200310122943 A CN200310122943 A CN 200310122943A CN 1556381 A CN1556381 A CN 1556381A
- Authority
- CN
- China
- Prior art keywords
- ground
- pixel
- image
- reflectivity
- calibration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Image Processing (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
The method includes three sub procedures: general layout of ground synchronous calibration experiment, capture of ground synchronous calibration data and indoor reflectivity conversion. Based on ground synchronous calibration data captured from ground synchronous calibration experiment, and selected sun elevation angle as basic reference, reflectivity conversion is carried out for aviation high spectrum remote sensing data. The method obtains accurate and reliable ground calibration data, and considers influence on radiance of ground object from time so as to raise precision of conversion. The invention also solves issues of too long time for capturing aerial line images or images not containing ground calibration points.
Description
Technical field
The present invention relates to the sensor information process field, particularly a kind of method of utilizing ground synchronous calibration the aviation high-spectrum remote sensing data to be carried out reflectivity conversion.
Background technology
High-spectrum remote sensing data has higher resolution on the spectrum dimension, the comparatively detailed spectral characteristic of complete sum of ground object target can be provided.Because be subjected to the influence of atmosphere, spectral signal that remote sensor receives and atural object real spectrum characteristic have certain deviation; The reflectance spectrum characteristic is intrinsic physics of atural object and chemical feature index in addition, the intrinsic propesties that more can reflect atural object than radiation wave spectrum.Therefore, have only high-spectrum remote sensing data is converted to reflectivity, could more effectively carry out follow-up applied analysis.For high-spectrum remote-sensing was used, reflectivity conversion was a very important element task.
Reflectivity conversion method commonly used at present comprises utilizes radiation transfer equation to carry out reflectivity conversion, utilize image itself to carry out reflectivity conversion such as flat field model, inherent average relative reflectance model, logarithm residual error model etc., and the method for carrying out reflectivity conversion by ground actual measurement reflectance spectrum data.These three class methods respectively have characteristics: radiation transfer theory is comparatively ripe, but needs to know atmospheric parameter and relevant data in concrete the application, and attitude instability in the airborne remote sensing system operational process does not have real-time atmospheric parameter record, is difficult to use.It is fairly simple to utilize image itself to carry out the method for reflectivity conversion, but what obtain is relative reflectance.Therefore carrying out reflectivity conversion by the ground synchronous calibration measured data is easy, the most effective and the most practical method that operates.Analyze to have now and test the method for carrying out reflectivity conversion, in concrete the application, generally also have following problem by ground synchronous calibration:
1) ground synchronous calibration experimental arrangement error or unreasonable
The ground synchronous calibration experiment is the important content that carries out reflectivity conversion by ground data, does not have ground synchronous calibration data accurately and reliably, just can not obtain high-precision reflectivity.This is a very numerous and diverse job, comprises artificial and the determining, manually calibrate the layout of atural object, calibrate layouting of atural object naturally of calibration type of ground objects naturally, and the obtaining etc. of ground synchronous calibration data.But in real world applications, exist such as causing the selection of calibration type of ground objects improper without wave spectrum experiment fully and analysis, calibration atural object size design is improper to make the influence that can't effectively eliminate mixed pixel in application, and the ground data synchronism of obtaining is crossed problem such as difference and made the ground calibration data of obtaining not accurate enough reliably.Therefore press for that a cover is complete, system arrange the ground synchronous calibration experiment carrying out ground synchronous calibration according to aviation flight, and then obtain the method for accurately reliable ground scalar number certificate.
2) the reflectivity conversion method is not considered the influence of time factor
Having now utilizes the ground synchronous data to carry out the method for reflectivity conversion, generally be that entire image is done as a whole the processing, do not consider the influence of time factor to atural object target emanation characteristic, this disposal route is comparatively suitable for the little small size test site of aerial remotely-sensed data acquisition time span.Aviation high-spectrum remote sensing data spatial resolution is generally higher, can not be to the large area region instantaneous imaging, and common aviation flight experiment can cover many course lines; And because the restriction of experiment condition, ground synchronous calibration can only concentrate on the nearer several representative points of certain bar airline distance.If therefore course line, ground calibration point place is long, when perhaps utilizing the ground calibration data to look like to carry out reflectivity conversion, consider that or else the influence of time factor will cause bigger deviation to not containing ground calibration point enroute chart.
Summary of the invention
The objective of the invention is to overcome above-mentioned the deficiencies in the prior art, a kind of aviation high-spectrum remote-sensing flight ground synchronous calibration and reflectivity conversion method are provided, to obtain accurately reliable ground scalar number certificate, and when reflectivity conversion, will the time influence of ground object target spoke brightness be taken into account, improve the precision of reflectivity conversion.
For achieving the above object, the invention provides a kind of aviation high-spectrum remote-sensing flight ground synchronous calibration and reflectivity conversion method, its characteristics are to comprise the following steps:
(1) determines the latitude and longitude coordinates of flight path and major control point;
(2) determine the scope and the ground calibration thing type of ground scaling point according to flight path;
(3) time of determining ground experiment is the moment that aircraft is crossed ground calibration point;
(4) carry out the ground synchronous calibration experiment, obtain ground synchronous calibration data, this synchronous calibration data comprises ground spectral data and ground calibration auxiliary data, and wherein the ground calibration auxiliary data comprises the longitude and latitude of time that each spectral data obtains and each calibration thing position;
(5) reflectivity conversion step:, make benchmark with the sun altitude of choosing the aviation high-spectrum remote sensing data is carried out reflectivity conversion according to the ground synchronous calibration data of obtaining.
Described ground calibration thing comprises artificial and natural ground calibration thing, and wherein artificial ground calibration thing comprises as the byssus of bright body with as the black cotton of black matrix; Naturally the calibration thing comprises concrete floor as bright body, water body and as the dump of black matrix, its area all is not less than 10 * 10 pixel areas, and is positioned at point under the machine.
Described reflectivity conversion step comprises the following steps:
When (1) determining that the aviation high-spectrum remote sensing carried out reflectivity conversion, be elected to be the sun altitude of benchmark;
(2) basis is the reflectivity when entire image is converted to same sun altitude, still image pixel value is converted to the reflectivity constantly of imaging separately, determines the spoke brightness value L of each pixel of image
S λ
(3) determine the inverting coefficient k of reflectivity
λAnd b
λ
(4) with the inverting coefficient k
λAnd b
λBe applied to the spoke brightness value L of step (2) gained according to the form of ρ=kL+b
S λ, can obtain reflectivity data ρ
G λ
Described inverting coefficient k
λAnd b
λBe the spoke brightness value of the respective regions that extracts according to ground actual measurement synchronous calibration data and from image, utilize least square method to calculate.
Described aviation high-spectrum remote sensing data is the enroute chart picture that comprises ground calibration point, and when corresponding sun altitude was benchmark during with the imaging of ground calibration point, described reflectivity conversion step comprised the following steps:
(1) determine the sun altitude of each pixel of enroute chart picture, selecting the corresponding pixel area of ground scaling point on image is the interesting image district, determines and with the pairing average sun altitude h of all images region of interest
tBe benchmark;
(2), determine that entire image calculating entire image is at average sun altitude h according to following formula
tSpoke brightness value L under the image-forming condition
S λ t:
Wherein, h
iBe the pairing sun altitude of pixel i; L
S λ iIt is the original spoke brightness value of entire image correspondence;
(3) the spoke brightness calculation inverting coefficient that utilizes ground synchronous calibration data and extract from the interesting image district comprises the following steps:
1. extracting each interesting image district sun altitude is h
tThe time each passage pairing spoke brightness value;
2. channel parameters each corresponding ground, the interesting image district spectral data that experiment is obtained to ground synchronous calibration according to imaging spectrometer resamples, and makes ground wave spectrum and image wave spectrum have identical wave band number, the ground wave spectrum after obtaining to resample;
3. pairing image wave spectrum in each interesting image district and ground wave spectrum are matched, and set up system of equations, and try to achieve the inverting coefficient k of λ wave band correspondence by following equation
λAnd b
λ:
ρ
gλR=k
λ·L
sλR+b
λ
ρ wherein
G λ RBe the ground wave spectrum after resampling; L
S λ RThe sun altitude that extracts for image respective image region of interest is h
tThe spoke brightness value;
(4) the inverting coefficient k
λAnd b
λEquation shown in (3) is used for the simulation spoke brightness value L that step (2) obtains set by step
S λ t, can obtain entire image is h at sun altitude
tReflectivity ρ
G λ t
Described aviation high-spectrum remote sensing data is the course line image pixel that comprises ground calibration point, and with each pixel imaging of image when corresponding sun altitude is benchmark constantly, described reflectivity conversion step comprises the following steps:
(1) determine the sun altitude of each pixel of enroute chart picture, selecting the corresponding pixel area of ground scaling point on image is the interesting image district, extracts the spoke brightness value L in each interesting image district from image
S λ t, determine and with the average sun altitude h of its imaging
tBe benchmark;
(2) calculate each interesting image district at the corresponding sun altitude h of pixel i according to following formula
tSimulation spoke brightness value L under the condition
S λ 0i:
(3) utilize ground synchronous calibration data and the spoke brightness value that extracts from the interesting image district calculates the inverting coefficient of each pixel, comprise the following steps:
1. channel parameters each corresponding ground, the interesting image district spectral data that experiment is obtained to ground synchronous calibration according to imaging spectrometer resamples, and makes ground wave spectrum and image wave spectrum have identical wave band number, the ground wave spectrum after obtaining to resample;
2. for pixel i, the simulation spoke brightness value and resampling ground, the ground wave spectrum in each interesting image district matched, and set up system of equations by following relation, try to achieve the corresponding inverting coefficient of each pixel:
ρ
gλ0=k
i·L
sλ0i+b
i
Wherein, ρ
G λ 0Be the ground wave spectrum after resampling; L
S λ 0iSimulation spoke brightness value from the extraction of respective image region of interest; k
iAnd b
iBe the pairing reflectivity inverting of the pixel i that is asked coefficient.
(4) trying to achieve inverting coefficient k
iAnd b
iBe applied to the original spoke brightness value of pixel i as follows, promptly obtain the reflectivity ρ of pixel i
G λ i:
ρ
gλi=k
i·L
sλi+b
i
L wherein
S λ iThe spoke brightness value that receives for pixel i position remote sensor.
Described aviation high-spectrum remote sensing data is not for containing the course line image pixel of ground calibration point, and with image each when pixel imaging moment, corresponding sun altitude was benchmark, described reflectivity conversion step may further comprise the steps:
(1) determines the sun altitude of each pixel of enroute chart picture, from containing the spoke brightness value L in each interesting image district that the ground calibration dot image reads
S λ tAnd corresponding sun altitude h
t
(2) for pixel i, establishing its sun altitude is h
i, calculate each interesting image district as follows at this sun altitude h
iCorresponding down simulation radiance L
S λ 0t:
(3) utilize the spoke brightness value in ground synchronous calibration data and interesting image district to calculate the inverting coefficient, comprise the following steps:
1. channel parameters each corresponding ground, the interesting image district spectral data that experiment is obtained to ground synchronous calibration according to imaging spectrometer resamples, and makes ground wave spectrum and image wave spectrum have identical wave band number, the ground wave spectrum after obtaining to resample;
2. for pixel i, the simulation spoke brightness value and resampling ground, the ground wave spectrum in each interesting image district matched, and set up system of equations by following relation, try to achieve the corresponding reflectivity inverting of each pixel i coefficient k
iAnd b
i:
ρ
gλ0=k
i·L
sλ0i+b
i;
Wherein, ρ
G λ 0Be the ground wave spectrum after resampling; L
S λ 0iSimulation spoke brightness value for the respective image region of interest;
(4) with the inverting coefficient k that obtains
iAnd b
iBe applied to the original spoke brightness value of pixel i as follows, promptly obtain the reflectivity ρ of pixel i
G λ i:
ρ
gλi=k
i·L
sλi+b
i
L wherein
S λ iThe spoke brightness value that receives for pixel i position remote sensor.
Described aviation high-spectrum remote sensing data is not for containing the enroute chart picture of ground calibration point, and when corresponding sun altitude was benchmark during with the imaging of ground calibration point, described reflectivity conversion step may further comprise the steps:
(1) determines the sun altitude of each pixel of enroute chart picture, read existing inverting coefficient k
λAnd b
λAnd pairing imaging sun altitude h
t
(2), entire image is converted to sun altitude h according to following formula
tSimulation spoke brightness value L under the image-forming condition
S λ t:
H wherein
iSun altitude for pixel i; L
S λ iBe the original spoke brightness value of entire image;
(3) with the inverting coefficient k
λAnd b
λBe applied to the simulation spoke luminance picture of above-mentioned steps (2) gained according to the form of ρ=kL+b, can obtain entire image is h at altitude of the sun
tReflectivity.
Described sun altitude is that promptly the imaging time of each pixel correspondence and latitude and longitude information calculate according to remote sensing flight auxiliary data.
The present invention tests ground synchronous calibration before ground synchronous calibration experiment and has carried out scientific and reasonable general layout, select suitable calibration type of ground objects through wave spectrum experiment and analysis, suitably the atural object size is calibrated in design, effectively subdue the influence of mixed pixel, assurance obtains the synchronism of ground data, thereby guaranteed the accurately reliable of the ground synchronous calibration data obtained in the ground synchronous calibration experiment; The ground synchronous calibration data that the present invention simultaneously utilizes the ground synchronous calibration realization to obtain is carried out reflectivity conversion to the aviation high-spectrum remote sensing data, corresponding sun altitude is that benchmark carries out reflectivity conversion to the aviation high-spectrum remote sensing data during promptly with the imaging of ground calibration point, to the time influence of ground object target spoke brightness be taken into account, it is long or when not containing ground calibration point enroute chart and look like to carry out reflectivity conversion to have overcome course line, ground calibration point place, do not consider that time factor will bring the problem than large deviation, improve the precision of reflectivity conversion.
The present invention is further illustrated below in conjunction with accompanying drawing and embodiment.
Description of drawings
Fig. 1 is a workflow diagram of the present invention.
Fig. 2 is the workflow diagram of airborne remote sensing flight ground synchronous calibration experimental set-up of the present invention.
Fig. 3 is the process flow diagram that the aviation high-spectrum remote sensing data is carried out the embodiment 1 of reflectivity conversion of the present invention.
Fig. 4 is the process flow diagram that the aviation high-spectrum remote sensing data is carried out the embodiment 2 of reflectivity conversion of the present invention.
Fig. 5 is the method flow diagram of calculating sun altitude of the present invention.
Fig. 6 is the process flow diagram that the aviation high-spectrum remote sensing data is carried out the embodiment 3 of reflectivity conversion of the present invention.
Fig. 7 is the process flow diagram that the aviation high-spectrum remote sensing data is carried out the embodiment 4 of reflectivity conversion of the present invention.
Embodiment
Aviation high-spectrum remote-sensing flight ground synchronous calibration of the present invention and reflectivity conversion method utilize ground synchronous calibration that the aviation high-spectrum remote sensing data is carried out reflectivity conversion, as shown in Figure 1, this method comprises obtaining and three subprocess of indoor reflection rate conversion of ground synchronous calibration experimental set-up, ground synchronous calibration data.
Wherein ground synchronous calibration experimental set-up process is earlier determined flight path according to application demand as shown in Figure 2, and provides the latitude and longitude coordinates of major control point, particularly turning point or the latitude and longitude coordinates of interesting target particularly.
Secondly determine the scope of ground scaling point and suitable ground calibration thing type (abbreviation type of ground objects), wherein the ground calibration thing comprises artificial and calibrates two kinds of atural objects naturally.In conjunction with the line of flight, choose the suitable position of layouting, as the crossing on the not too crowded harbour of traffic, road, park etc. so that artificial calibration thing is put in existing enough spaces, comprise again satisfy condition calibrate atural object naturally.The atural object that can be used as scaling point requires to meet the following conditions: significantly do not absorb or reflection, have good lambert's bulk properties; Spectral characteristic is vary stable in time; Size is suitable, and the area that is not less than 10 * 10 pixels is arranged, to reduce the area of mixed pixel; The scaling point number is 5~7, and its reflectivity covers high, medium and low reflectivity range; Scaling point is positioned at point under the machine, to avoid how much and radiometric distortion; The distance of scaling point distance will guarantee synchronism to greatest extent testing in the acceptable scope synchronously.Choose artificial calibration atural object through the preparation of indoor wave spectrum comparative analysis and man-made features; Through field reconnaissance and ground-object spectrum comparative analysis, choose nature calibration atural object, determine the position that it is concrete, and the longitude and latitude of this position is in time fed back to aviation flight organization unit cover this zone to guarantee the course line.
Determine the time of ground experiment at last: flight was accurately obtained the concrete time of taking off on same day, spent the moment of ground calibration point to determine aircraft.
Then carry out the ground synchronous calibration experiment, enter the acquisition process of ground synchronous calibration data, this is the prerequisite of the aviation high-spectrum remote sensing data being carried out reflectivity conversion.This synchronous calibration data comprises ground spectral data and ground calibration auxiliary data, and wherein the ground calibration auxiliary data comprises time that each spectral data obtains and the longitude and latitude that each calibrates the thing position etc.
Obtaining of above ground synchronous calibration experimental set-up and ground synchronous calibration data can describe in detail by following examples.
Carry out the high-spectrum remote-sensing flight experiment of OMIS (practical modularization imaging spectrometer) and PHI (pull-broom type area array CCD hyperspectral imager) in the area, Shanghai, before the flight, determine following several contents:
1) determine flight path according to application demand, Dianshan Lake → Xu Pu bridge → Nanpu bridge → Yangpu Bridge (along Pudong) → Jiang Wan airport (area) → Yangpu Bridge → Nanpu bridge (Yan Puxi), and provide the latitude and longitude coordinates of major control point.
2) determine the approximate range of scaling point and suitable type of ground objects
In order to select ground synchronous calibration experiment man-made features, white, the black cloth of differing textures such as polyester-cotton blend, cotton yarn card, Nylon Taffeta, velveteen, cotton textiles, Lycra, silk, silk fabric carried out indoor spectrum comparative analysis experiment, 12 * 12m is processed in selection at last
2Black cotton and byssus respectively as calibration with black matrix and bright body.
Through field reconnaissance and ground-object spectrum comparative analysis, select concrete floor, vegetation, floor tile, marble, dump, soil, water body as calibrating atural object naturally, determine the position that it is concrete, and the longitude and latitude of position is in time fed back to aviation flight organization unit.
3) step is determined the time of ground experiment
Flight was accurately obtained the concrete time of taking off on same day, spent the moment of ground calibration point to determine aircraft.
4) carry out the ground synchronous calibration experiment, obtain ground synchronous calibration data
Aviation flight same day, utilize black, bright body and various natural feature on a map, ground such as section, Jiang Wan airport and park, square, Nanpu have carried out the ground synchronous calibration experiment in the thing outbeach, on the south the Xu Pu bridge, obtain the wave spectrum and the associated auxiliary data of various natural feature on a map targets such as concrete floor, vegetation, floor tile, marble, dump, soil, water body, for the processing of aviation high-spectrum remote sensing data provides data source.
After obtaining ground synchronous calibration data,, make benchmark with the sun altitude of choosing the aviation high-spectrum remote sensing data is carried out reflectivity conversion according to this ground synchronous calibration data.An aviation high-spectrum remote-sensing flight generally can cover several course lines, and because the restriction of experiment condition, ground synchronous calibration often can only concentrate on the nearer several representative points of certain bar airline distance.Therefore utilize reflectivity conversion method proposed by the invention to look like to contain or do not contain the ground calibration point with regard to the enroute chart of proofreading and correct respectively below and provide concrete example.
1, the enroute chart picture that contains ground calibration point
For the reflectivity conversion that contains ground calibration point enroute chart picture, corresponding sun altitude is that benchmark carries out reflectivity conversion to entire image in the time of can the imaging of ground calibration point, and the specific implementation step is seen Fig. 3; Image pixel can also be converted to the reflectivity under the image-forming condition separately, the specific implementation step is seen Fig. 4.Details are as follows for it:
1) corresponding sun altitude is benchmark (embodiment 1) during the imaging of ground calibration point
Each pixel sun altitude of computed image at first; Next select the ground scaling point, calculate the average sun altitude of ground calibration point, as the benchmark of entire image reflection conversion; Calculate the spoke brightness value of entire image under this benchmark altitude of the sun corner condition; Utilize ground synchronous calibration data and (Region of Interest, be called for short ROI) extracts from the interesting image district spoke brightness value calculates the inverting coefficient by least square method; At last the inverting coefficient is applied to entire image and obtains reflectivity.As shown in Figure 3, it specifically may further comprise the steps:
Step 1: read the enroute chart picture that contains ground calibration point that will carry out reflectivity conversion, requirement is the spoke brightness value.
Step 2: according to the corresponding remote sensing of reading images flight auxiliary data, obtain the imaging time and the latitude and longitude information of each pixel correspondence, and calculate the pairing sun altitude of every pixel, its calculation procedure as shown in Figure 5:
Step 21: the year, month, day according to imaging is calculated a day angle J;
Step 22; Calculate declination angle δ by day angle J;
Step 23: calculate solar hour angle τ, performing step is:
Step 231: according to longitude information, S changes local time S into during with Beijing of imaging
d,
Step 232: calculate the time difference by day angle J, and local time is carried out the time difference correct, get very
Solar time,
Step 233: calculate solar hour angle τ;
Step 24: calculate too according to latitude information and gained declination angle δ, solar hour angle τ
Sun elevation angle h.
Generally speaking, the auxiliary data precision of airborne remote sensing flight is limited, but interpolation is obtained the supplementary of every row or every pixel necessary the time.
Step 3: display image, select the corresponding pixel area of ground scaling point on image, i.e. interesting image district (ROI district).
Step 4: the pairing average sun altitude h in all ROI districts of statistical computation
t
Step 5:, entire image is converted to average sun altitude h according to following formula
tSpoke brightness value L under the image-forming condition
S λ t:
Wherein, h
iBe the pairing sun altitude of pixel i; L
S λ tIt is the original spoke brightness value of entire image correspondence; L
S λ tFor entire image is h at sun altitude
tThe time pairing spoke brightness value.
Step 6: extracting each ROI district sun altitude is h
tThe time each passage correspondence the spoke brightness value.
Step 7: read each ROI district corresponding ground return wave spectrum data that the ground synchronous calibration experiment is obtained.
Step 8: the channel parameters according to imaging spectrometer resamples to the ground return wave spectrum data that it obtained, and makes ground wave spectrum and image wave spectrum have identical wave band number.
Step 9: the ground return wave spectrum after obtaining to resample.
Step 10: pairing image wave spectrum in each ROI district and ground wave spectrum are matched.
Step 11: according to step 10, set up system of equations, and try to achieve the inverting coefficient k by following equation
λAnd b
λ:
ρ
gλR=k
λ·L
sλR+b
λ
Wherein, ρ
G λ RBe the ground actual measurement reflectance spectrum after resampling; L
S λ RThe sun altitude that extracts for the corresponding ROI of image district is h
tThe spoke brightness value; k
λAnd b
λInverting coefficient for the λ wave band correspondence that requires.
Step 12: at last again k
λAnd b
λEquation shown in 11 is used for the spoke brightness value that obtains of step 5 set by step, and can obtain entire image is h at sun altitude
tReflectivity.
2) image pixel is converted to reflectivity (embodiment 2) under each pixel image-forming condition
At first calculate each pixel sun altitude; Select the ground scaling point, calculate the spoke brightness value and the average sun altitude in each ROI district; For every pixel, calculate the simulation spoke brightness value of each ROI district under this sun altitude; Utilize ground calibration data and simulation spoke brightness value to calculate the inverting coefficient of each pixel; At last the inverting coefficient is used for each pixel and promptly gets reflectivity.As shown in Figure 4, it specifically may further comprise the steps:
Step 1: read the enroute chart picture that contains ground calibration point that will carry out reflectivity conversion, requirement is the spoke brightness value.
Step 2: by the sun altitude of the every pixel of step computed image shown in Figure 5.
Step 3: display image, select the corresponding pixel area of ground scaling point on image, i.e. interesting image district (ROI district).
Step 4:, extract its spoke brightness value L from image for each ROI district
S λ t, and the average sun altitude h of its imaging of statistical computation
t,
Step 5: calculate it as follows at the corresponding sun altitude h of pixel i
tSimulation spoke brightness value L under the condition
S λ 0t:
Step 6: read each ROI district corresponding ground return wave spectrum data that the ground synchronous calibration experiment is obtained.
Step 7: the channel parameters according to imaging spectrometer resamples to the ground return wave spectrum data that it obtained, and makes ground wave spectrum and image wave spectrum have identical wave band number.
Step 8: the ground return wave spectrum after obtaining to resample.
Step 9: for pixel i, resampling ground, the ground wave spectrum that simulation spoke brightness value that the step 5 in each ROI district is obtained and step 8 obtain matches.
Step 10: according to step 9, set up system of equations, try to achieve the corresponding inverting coefficient of each pixel by following relation:
ρ
gλ0=k
i·L
sλ0i+b
i
Wherein, ρ
G λ 0Be the ground calibration point actual measurement reflectance spectrum after resampling; L
S λ 0iThe simulation spoke brightness value that extracts from corresponding ROI district; k
iAnd b
iBe the pairing reflectivity inverting of the pixel i that is asked coefficient.
Step 11: the last k that tries to achieve
iAnd b
iBe applied to the original spoke brightness value of pixel i as follows, promptly obtain its reflectivity:
ρ
gλt=k
i·L
sλt+b
i
Wherein, ρ
G λ iBe the pairing reflectivity of pixel i; L
S λ iThe spoke brightness value that receives for pixel i position remote sensor.
2, the enroute chart picture that does not comprise ground calibration point
For the reflectivity conversion that does not contain ground calibration point enroute chart picture that flight experiment obtained the same day, can import according to the user, entire image is converted to reflectivity under certain sun altitude image-forming condition, its specific implementation step is seen Fig. 6; Image pixel can also be converted to the reflectivity under the image-forming condition separately, its specific implementation step is seen Fig. 7.
1), entire image is converted to reflectivity (embodiment 3) under certain sun altitude image-forming condition according to user input
At first read existing reflectivity inverting coefficient and the pairing sun altitude of this inverting coefficient; The radiance value of computed image under this sun altitude image-forming condition; Utilize linear model, obtain the reflectivity of new enroute chart picture under this altitude of the sun corner condition that view picture does not contain ground calibration point.As shown in Figure 6, it specifically may further comprise the steps:
Step 1: read the enroute chart picture that contains ground calibration point that will carry out reflectivity conversion, requirement is the spoke brightness value.
Step 2: by the sun altitude of the every pixel of step computed image shown in Figure 5.
Step 3: read existing reflectivity conversion coefficient k
λAnd b
λ, and pairing imaging sun altitude h
t
Step 4:, entire image is converted to sun altitude h according to following formula
tSimulation spoke brightness value under the image-forming condition:
Wherein, h
iSun altitude for pixel i; L
S λ iBe the original spoke brightness value of its image; L
S λ tFor it is h at sun altitude
tThe time pairing spoke brightness value.
Step 5: with conversion coefficient k
λAnd b
λ, be applied to the simulation spoke brightness value that step 4 obtains according to the form of ρ=kL+b, can obtain entire image is h at sun altitude
tReflectivity.
2) image pixel is converted to reflectivity (embodiment 4) under each pixel image-forming condition
At first calculate each pixel sun altitude; Read then by containing each ROI district spoke brightness value and sun altitude that ground calibration point enroute chart picture extracts; For each pixel, calculate the simulation spoke brightness value of each ROI district under this sun altitude; Utilize ground calibration data and simulation spoke brightness value to calculate the inverting coefficient of each pixel; At last the inverting coefficient is used for each pixel and promptly gets reflectivity.As shown in Figure 7, it specifically may further comprise the steps:
Step 1: read the enroute chart picture that contains ground calibration point that will carry out reflectivity conversion, requirement is the spoke brightness value.
Step 2: by the sun altitude of the every pixel of step computed image shown in Figure 5.
Step 3: read from containing the radiation wave spectrum and the corresponding sun altitude in each ROI district that the ground calibration dot image obtains.
Step 4: for pixel i, establishing its sun altitude is h
t, calculate each ROI district as follows at this sun altitude h
tCorresponding down simulation spoke brightness value:
Wherein, L
S λ tFor being h at sun altitude from containing the ROI district that the ground calibration dot image extracts
tThe time the spoke brightness value, L
S λ 0tFor this ROI district is h at sun altitude
iThe time the spoke brightness value.
Step 5: read each corresponding ground, ROI district spectral data that the ground synchronous calibration experiment is obtained.
Step 6: the channel parameters according to imaging spectrometer resamples to it, makes ground wave spectrum and image wave spectrum have identical wave band number.
Step 7: the ground return wave spectrum after obtaining to resample.
Step 8: for pixel i, resampling ground, the ground wave spectrum that simulation spoke brightness value that the step 4 in each ROI district is obtained and step 7 obtain matches.
Step 9: according to step 8, set up system of equations, try to achieve the corresponding inverting coefficient of each pixel by following relation:
ρ
gλ0=k
i·L
sλ0i+b
i
Wherein, ρ
G λ 0Be the ground actual measurement reflectance spectrum after resampling; L
S λ 0iSimulation spoke brightness value for corresponding ROI district; k
tAnd b
iBe the pairing reflectivity inverting of pixel i coefficient.
Step 10: last k
iAnd b
iBe applied to the original spoke brightness value of pixel i as follows, promptly obtain its reflectivity:
ρ
gλi=k
i·L
sλi+b
i
Wherein, ρ
G λ iBe the pairing reflectivity of pixel i; L
S λ iThe spoke brightness value that receives for pixel i position remote sensor.
The present invention before ground synchronous calibration experiment to ground synchronous calibration experiment carried out scientific and reasonable general layout, guaranteed the ground synchronous calibration data obtained in the ground synchronous calibration experiment accurately not reliably; The ground synchronous calibration data that the present invention simultaneously utilizes the ground synchronous calibration experiment to obtain is carried out reflectivity conversion to the aviation high-spectrum remote sensing data, to the time influence of ground object target spoke brightness be taken into account, it is long or when not containing ground calibration point enroute chart and look like to carry out reflectivity conversion to have overcome course line, ground calibration point place, do not consider that time factor will bring the problem than large deviation, improve the precision of reflectivity conversion.When the present invention has proposed with the imaging of ground calibration point when reflectivity conversion in addition corresponding sun altitude to be benchmark carry out reflectivity conversion and image pixel be converted to the method for reflectivity under the image-forming condition separately entire image, and look like to contain ground calibration point according to desire conversion enroute chart and handle respectively.
Claims (9)
1, a kind of aviation high-spectrum remote-sensing flight ground synchronous calibration and reflectivity conversion method is characterized in that comprising the following steps:
(1) determines the latitude and longitude coordinates of flight path and major control point;
(2) determine the scope and the ground calibration thing type of ground scaling point according to flight path;
(3) time of determining ground experiment is the moment that aircraft is crossed ground calibration point;
(4) carry out the ground synchronous calibration experiment, obtain ground synchronous calibration data, this synchronous calibration data comprises ground spectral data and ground calibration auxiliary data, and wherein the ground calibration auxiliary data comprises the longitude and latitude of time that each spectral data obtains and each calibration thing position;
(5) reflectivity conversion step; According to the ground synchronous calibration data of obtaining, make benchmark with the sun altitude of choosing the aviation high-spectrum remote sensing data is carried out reflectivity conversion.
2, aviation high-spectrum remote-sensing flight ground synchronous calibration according to claim 1 and reflectivity conversion method, it is characterized in that described ground calibration thing comprises artificial and natural ground calibration thing, wherein artificial ground calibration thing comprises as the byssus of bright body with as the black cotton of black matrix; Naturally the calibration thing comprises concrete floor as bright body, water body and as the dump of black matrix, its area all is not less than 10 * 10 pixel areas, and is positioned at point under the machine.
3, aviation high-spectrum remote-sensing flight ground synchronous calibration according to claim 1 and 2 and reflectivity conversion method is characterized in that described reflectivity conversion step comprises the following steps:
When (1) determining that the aviation high-spectrum remote sensing carried out reflectivity conversion, be elected to be the sun altitude of benchmark;
(2) basis is the reflectivity when entire image is converted to same sun altitude, still image pixel value is converted to the reflectivity constantly of imaging separately, determines the spoke brightness value L of each pixel of image
S λ
(3) determine the inverting coefficient k of reflectivity
λAnd b
λ
(4) with the inverting coefficient k
λAnd b
λBe applied to the spoke brightness value L of step (2) gained according to the form of ρ=kL+b
S λ, can obtain reflectivity data ρ
G λ
4, aviation high-spectrum remote-sensing flight ground synchronous calibration according to claim 3 and reflectivity conversion method is characterized in that described inverting coefficient k
λAnd b
λBe the spoke brightness value of the respective regions that extracts according to ground actual measurement synchronous calibration data and from image, utilize least square method to calculate.
5, aviation high-spectrum remote-sensing flight ground synchronous calibration according to claim 4 and reflectivity conversion method, it is characterized in that described aviation high-spectrum remote sensing data is the enroute chart picture that comprises ground calibration point, and when corresponding sun altitude was benchmark during with the imaging of ground calibration point, described reflectivity conversion step comprised the following steps:
(1) determine the sun altitude of each pixel of enroute chart picture, selecting the corresponding pixel area of ground scaling point on image is the interesting image district, determines and with the pairing average sun altitude h of all images region of interest
tBe benchmark;
(2), determine that entire image calculating entire image is at average sun altitude h according to following formula
tSpoke brightness value L under the image-forming condition
S λ t:
Wherein, h
iBe the pairing sun altitude of pixel i; L
S λ iIt is the original spoke brightness value of entire image correspondence;
(3) the spoke brightness calculation inverting coefficient that utilizes ground synchronous calibration data and extract from the interesting image district comprises the following steps:
1. extracting each interesting image district sun altitude is h
tThe time each passage pairing spoke brightness value;
2. channel parameters each corresponding ground, the interesting image district spectral data that experiment is obtained to ground synchronous calibration according to imaging spectrometer resamples, and makes ground wave spectrum and image wave spectrum have identical wave band number, the ground wave spectrum after obtaining to resample;
3. pairing image wave spectrum in each interesting image district and ground wave spectrum are matched, and set up system of equations, and try to achieve the inverting coefficient k of λ wave band correspondence by following equation
λAnd b
λ:
ρ
gλR=k
λ·L
sλR+b
λ
ρ wherein
G λ RBe the ground wave spectrum after resampling; L
S λ RThe sun altitude that extracts for image respective image region of interest is h
tThe spoke brightness value;
(4) the inverting coefficient k
λAnd b
λEquation shown in (3) is used for the simulation spoke brightness value L that step (2) obtains set by step
S λ t, can obtain entire image is h at sun altitude
tReflectivity ρ
G λ t
6, aviation high-spectrum remote-sensing flight ground synchronous calibration according to claim 4 and reflectivity conversion method, it is characterized in that described aviation high-spectrum remote sensing data is the course line image pixel that comprises ground calibration point, and when corresponding sun altitude was benchmark constantly with each pixel imaging of image, described reflectivity conversion step comprised the following steps:
(1) determine the sun altitude of each pixel of enroute chart picture, selecting the corresponding pixel area of ground scaling point on image is the interesting image district, extracts the spoke brightness value L in each interesting image district from image
S λ t, determine and with the average sun altitude h of its imaging
tBe benchmark;
(2) calculate each interesting image district at the corresponding sun altitude h of pixel i according to following formula
iSimulation spoke brightness value L under the condition
S λ 0i:
(3) utilize ground synchronous calibration data and the spoke brightness value that extracts from the interesting image district calculates the inverting coefficient of each pixel, comprise the following steps:
1. channel parameters each corresponding ground, the interesting image district spectral data that experiment is obtained to ground synchronous calibration according to imaging spectrometer resamples, and makes ground wave spectrum and image wave spectrum have identical wave band number, the ground wave spectrum after obtaining to resample;
2. for pixel i, the simulation spoke brightness value and resampling ground, the ground wave spectrum in each interesting image district matched, and set up system of equations by following relation, try to achieve the corresponding inverting coefficient of each pixel:
ρ
gλ0=k
i·L
sλ0i+b
i
Wherein, ρ
G λ 0Be the ground wave spectrum after resampling; L
S λ 0iSimulation spoke brightness value from the extraction of respective image region of interest; k
tAnd b
iBe the pairing reflectivity inverting of the pixel i that is asked coefficient.
(4) trying to achieve inverting coefficient k
iAnd b
iBe applied to the original spoke brightness value of pixel i as follows, promptly obtain the reflectivity ρ of pixel i
G λ i:
ρ
G λ i=k
iL
S λ i+ b
iL wherein
S λ iThe spoke brightness value that receives for pixel i position remote sensor.
7, aviation high-spectrum remote-sensing flight ground synchronous calibration according to claim 4 and reflectivity conversion method, it is characterized in that described aviation high-spectrum remote sensing data is not for containing the course line image pixel of ground calibration point, and when corresponding sun altitude was benchmark constantly with each pixel imaging of image, described reflectivity conversion step may further comprise the steps:
(1) determines the sun altitude of each pixel of enroute chart picture, from containing the spoke brightness value L in each interesting image district that the ground calibration dot image reads
S λ tAnd corresponding sun altitude h
t
(2) for pixel i, establishing its sun altitude is h
i, calculate each interesting image district as follows at this sun altitude h
iCorresponding down simulation radiance L
S λ 0i:
(3) utilize the spoke brightness value in ground synchronous calibration data and interesting image district to calculate the inverting coefficient, comprise the following steps:
1. channel parameters each corresponding ground, the interesting image district spectral data that experiment is obtained to ground synchronous calibration according to imaging spectrometer resamples, and makes ground wave spectrum and image wave spectrum have identical wave band number, the ground wave spectrum after obtaining to resample;
2. for pixel i, the simulation spoke brightness value and resampling ground, the ground wave spectrum in each interesting image district matched, and set up system of equations by following relation, try to achieve the corresponding reflectivity inverting of each pixel i coefficient k
iAnd b
i:
ρ
G λ 0=k
iL
S λ 0i+ b
iWherein, ρ
G λ 0Be the ground wave spectrum after resampling; L
S λ 0iSimulation spoke brightness value for the respective image region of interest;
(4) with the inverting coefficient k that obtains
iAnd b
iBe applied to the original spoke brightness value of pixel i as follows, promptly obtain the reflectivity ρ of pixel i
G λ t:
ρ
gλi=k
i·L
sλt+b
i
L wherein
S λ tThe spoke brightness value that receives for pixel i position remote sensor.
8, aviation high-spectrum remote-sensing flight ground synchronous calibration according to claim 3 and reflectivity conversion method, it is characterized in that described aviation high-spectrum remote sensing data is not for containing the enroute chart picture of ground calibration point, and when corresponding sun altitude was benchmark during with the imaging of ground calibration point, described reflectivity conversion step may further comprise the steps:
(1) determines the sun altitude of each pixel of enroute chart picture, read existing inverting coefficient k
λAnd b
λAnd pairing imaging sun altitude h
t
(2), entire image is converted to sun altitude h according to following formula
tSimulation spoke brightness value L under the image-forming condition
S λ t:
H wherein
iSun altitude for pixel i; L
S λ iBe the original spoke brightness value of entire image;
(3) with the inverting coefficient k
λAnd b
λBe applied to the simulation spoke luminance picture of above-mentioned steps (2) gained according to the form of ρ=kL+b, can obtain entire image is h at altitude of the sun
tReflectivity.
9, aviation high-spectrum remote-sensing flight ground synchronous calibration according to claim 1 and 2 and reflectivity conversion method, it is characterized in that described sun altitude is that promptly the imaging time of each pixel correspondence and latitude and longitude information calculate according to remote sensing flight auxiliary data.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2003101229433A CN100344948C (en) | 2003-12-30 | 2003-12-30 | Aviation high spectrum remote sensing flight ground synchronous scaling and reflectivity conversion method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2003101229433A CN100344948C (en) | 2003-12-30 | 2003-12-30 | Aviation high spectrum remote sensing flight ground synchronous scaling and reflectivity conversion method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1556381A true CN1556381A (en) | 2004-12-22 |
CN100344948C CN100344948C (en) | 2007-10-24 |
Family
ID=34338842
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2003101229433A Expired - Fee Related CN100344948C (en) | 2003-12-30 | 2003-12-30 | Aviation high spectrum remote sensing flight ground synchronous scaling and reflectivity conversion method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100344948C (en) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101887012A (en) * | 2010-06-28 | 2010-11-17 | 中国国土资源航空物探遥感中心 | Spectral reflectance peak decomposition based quantitative inversion method of hyperspectral remote sensing mineral content |
CN102346247A (en) * | 2010-07-30 | 2012-02-08 | 中国科学院遥感应用研究所 | Computing method and calibration method for HJ-1B B08 effective waveband width |
CN102346070A (en) * | 2010-07-30 | 2012-02-08 | 中国科学院遥感应用研究所 | HJ-1A satellite hyper-spectral imaging instrument orbit radiation scaling method |
CN102435309A (en) * | 2011-09-20 | 2012-05-02 | 北京农业信息技术研究中心 | Field reflectance calibration method and system of agricultural imaging hyperspectral spectrometer |
CN102538820A (en) * | 2011-12-13 | 2012-07-04 | 中国测绘科学研究院 | Calibration method of aerial remote sensing integrated system |
CN102607592A (en) * | 2012-02-24 | 2012-07-25 | 北京大学 | Remote sensing calibration comprehensive method and calibration equipment vehicle |
CN102788643A (en) * | 2012-07-13 | 2012-11-21 | 中国科学院长春光学精密机械与物理研究所 | Method for calibrating ontrack high-precision optical spectrum of space remote sensing optical spectrum instrument |
CN102830392A (en) * | 2012-08-29 | 2012-12-19 | 南京大学 | Relative radiation and luminosity correcting method for optical remote sensing data of moon |
CN104246640A (en) * | 2011-09-30 | 2014-12-24 | 洛克希德马丁公司 | Method and apparatus for dynamic air traffic trajectory synchronization |
CN104866585A (en) * | 2015-05-28 | 2015-08-26 | 华北计算技术研究所 | Comprehensive processing system for experimental flight-test data |
CN105241428A (en) * | 2015-08-27 | 2016-01-13 | 北京市遥感信息研究所 | Water depth retrieval method with hyper-spectrum |
CN106644081A (en) * | 2017-01-10 | 2017-05-10 | 首都师范大学 | Multispectral image calibration method |
CN107687900A (en) * | 2017-09-01 | 2018-02-13 | 中国科学院上海技术物理研究所 | One kind is applied to airborne thermal infrared imaging spectrometer atmospheric correction system and method |
CN108120510A (en) * | 2017-12-08 | 2018-06-05 | 中国科学院合肥物质科学研究院 | A kind of in-orbit absolute radiation calibration method of optical sensor based on reflection mirror array |
CN114113000A (en) * | 2021-12-14 | 2022-03-01 | 江苏省地质勘查技术院 | Method, system and equipment for improving spectrum reflectivity inversion based on empirical linear model |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000047949A1 (en) * | 1999-02-09 | 2000-08-17 | Leica Geosystems Ag | Method and system for taking air-borne photogrammetric pictures (photogrammetry prism) |
JP2000258154A (en) * | 1999-03-12 | 2000-09-22 | Nec Corp | Image pickup device |
EP1048928A1 (en) * | 1999-04-27 | 2000-11-02 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Spaceborne hot temperature event (HTE) detection arrangement and multi-satellite HTE detection system |
-
2003
- 2003-12-30 CN CNB2003101229433A patent/CN100344948C/en not_active Expired - Fee Related
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101887012A (en) * | 2010-06-28 | 2010-11-17 | 中国国土资源航空物探遥感中心 | Spectral reflectance peak decomposition based quantitative inversion method of hyperspectral remote sensing mineral content |
CN102346247A (en) * | 2010-07-30 | 2012-02-08 | 中国科学院遥感应用研究所 | Computing method and calibration method for HJ-1B B08 effective waveband width |
CN102346070A (en) * | 2010-07-30 | 2012-02-08 | 中国科学院遥感应用研究所 | HJ-1A satellite hyper-spectral imaging instrument orbit radiation scaling method |
CN102346247B (en) * | 2010-07-30 | 2013-08-07 | 中国科学院遥感应用研究所 | Computing method and calibration method for HJ-1B B08 effective waveband width |
CN102435309A (en) * | 2011-09-20 | 2012-05-02 | 北京农业信息技术研究中心 | Field reflectance calibration method and system of agricultural imaging hyperspectral spectrometer |
CN104246640B (en) * | 2011-09-30 | 2017-06-09 | 洛克希德马丁公司 | For the method and apparatus of dynamic air traffic trajectory synchronization |
CN104246640A (en) * | 2011-09-30 | 2014-12-24 | 洛克希德马丁公司 | Method and apparatus for dynamic air traffic trajectory synchronization |
CN102538820A (en) * | 2011-12-13 | 2012-07-04 | 中国测绘科学研究院 | Calibration method of aerial remote sensing integrated system |
CN102538820B (en) * | 2011-12-13 | 2015-05-20 | 中国测绘科学研究院 | Calibration method of aerial remote sensing integrated system |
CN102607592B (en) * | 2012-02-24 | 2014-11-26 | 北京大学 | Remote sensing calibration comprehensive method and calibration equipment vehicle |
CN102607592A (en) * | 2012-02-24 | 2012-07-25 | 北京大学 | Remote sensing calibration comprehensive method and calibration equipment vehicle |
CN102788643A (en) * | 2012-07-13 | 2012-11-21 | 中国科学院长春光学精密机械与物理研究所 | Method for calibrating ontrack high-precision optical spectrum of space remote sensing optical spectrum instrument |
CN102788643B (en) * | 2012-07-13 | 2014-03-12 | 中国科学院长春光学精密机械与物理研究所 | Method for calibrating ontrack high-precision optical spectrum of space remote sensing optical spectrum instrument |
CN102830392B (en) * | 2012-08-29 | 2014-09-17 | 南京大学 | Relative radiation and luminosity correcting method for optical remote sensing data of moon |
CN102830392A (en) * | 2012-08-29 | 2012-12-19 | 南京大学 | Relative radiation and luminosity correcting method for optical remote sensing data of moon |
CN104866585B (en) * | 2015-05-28 | 2018-05-22 | 华北计算技术研究所 | A kind of experiment test flight data total system |
CN104866585A (en) * | 2015-05-28 | 2015-08-26 | 华北计算技术研究所 | Comprehensive processing system for experimental flight-test data |
CN105241428A (en) * | 2015-08-27 | 2016-01-13 | 北京市遥感信息研究所 | Water depth retrieval method with hyper-spectrum |
CN106644081B (en) * | 2017-01-10 | 2018-12-28 | 首都师范大学 | A kind of multispectral image scaling method |
CN106644081A (en) * | 2017-01-10 | 2017-05-10 | 首都师范大学 | Multispectral image calibration method |
CN107687900A (en) * | 2017-09-01 | 2018-02-13 | 中国科学院上海技术物理研究所 | One kind is applied to airborne thermal infrared imaging spectrometer atmospheric correction system and method |
CN107687900B (en) * | 2017-09-01 | 2023-09-12 | 中国科学院上海技术物理研究所 | Atmospheric correction system and method suitable for airborne thermal infrared imaging spectrometer |
CN108120510A (en) * | 2017-12-08 | 2018-06-05 | 中国科学院合肥物质科学研究院 | A kind of in-orbit absolute radiation calibration method of optical sensor based on reflection mirror array |
CN108120510B (en) * | 2017-12-08 | 2021-01-12 | 中国科学院合肥物质科学研究院 | On-orbit absolute radiation calibration method of optical remote sensor based on reflector array |
CN114113000A (en) * | 2021-12-14 | 2022-03-01 | 江苏省地质勘查技术院 | Method, system and equipment for improving spectrum reflectivity inversion based on empirical linear model |
Also Published As
Publication number | Publication date |
---|---|
CN100344948C (en) | 2007-10-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100344948C (en) | Aviation high spectrum remote sensing flight ground synchronous scaling and reflectivity conversion method | |
CN1191462C (en) | Individual element photoelectric measuring device for pane object | |
CN1302277C (en) | Field crops nutrition diagnostic method | |
CN1188672C (en) | Individual element photoelectric measuring device for plane object | |
CN110599548A (en) | Camera calibration method and device, camera and computer readable storage medium | |
CN109816680B (en) | High-throughput calculation method for crop plant height | |
CN1898697A (en) | Techniques for modifying image field data | |
CN110987821B (en) | Hyperspectral rapid atmospheric correction parameterization method | |
RU2009140389A (en) | DEVICE AND METHODS OF EXAMINATION, EVALUATION AND CLASSIFICATION OF PRECIOUS STONES | |
CN101039439A (en) | Method and apparatus for realizing correction of white balance | |
CN1378643A (en) | Method for colour matching by means of electronic imaging device | |
CN101210876A (en) | Rice nutrient information measuring method based on visible/near infrared multispectral imaging | |
CN1756374A (en) | Have color temperature conversion method and device based on the color correction functions of pixel brightness | |
CN114384031A (en) | Satellite-air-ground hyperspectral remote sensing water body heavy metal pollution three-dimensional monitoring method | |
CN106468536A (en) | A kind of leaf area measurement method | |
CN101034475A (en) | Method for generating non-shadow satellite remote sensing orthographic digital image by computer | |
CN101034472A (en) | Landform transformation of satellite remote sensing digital image supported by GIS | |
CN1677085A (en) | Agricultural application integrating system for earth observation technique and its method | |
CN110689567B (en) | Method for measuring and calculating total leaf area of whole arbor plant | |
CN101034471A (en) | Landform transformation method for sonar remote sensing digital image of underwater digit land model | |
CN1885345A (en) | Analysis method of digital image color analysis system | |
Zhang et al. | Analysis of vegetation indices derived from aerial multispectral and ground hyperspectral data | |
McCarthy et al. | Digital analysis of lichen cover: a technique for use in lichenometry and licnenology | |
CN113063741B (en) | Urban building material extraction method and device based on spectral characteristics | |
CN101034476A (en) | Method for generating underwater non-shadow sonar remote sensing orthographic digital image by computer |
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: 20071024 Termination date: 20100201 |