CN103295251B - Infrared strong absorbing path image simulation method and device in one - Google Patents

Abstract

Embodiments provide infrared strong absorbing path image simulation method and device in one, spectroscopic data in the wavelength band that spectroscopic data in the wavelength band that second channel is corresponding is corresponding with first passage has linear relationship, therefore, by the wave band response function of first passage and the wave band response function of second channel, select target wavelength from sample library of spectra, and calculate response function corresponding to each target wavelength, and then determine the equivalent spectrum data of first passage and the equivalent spectrum data of second channel, finally can determine the linearly dependent coefficient of first passage and second channel, thus, known linear related coefficient and obtained by second channel remote sensing images time, the remote sensing images obtained by first passage can be simulated, therefore, the remote sensing images of non-strong absorbing path can be converted to the remote sensing images of strong absorbing path, because the remote sensing images of non-strong absorbing path relatively easily obtain, therefore, facilitate the simulation of the remote sensing images to strong absorbing path.

Description

Infrared strong absorbing path image simulation method and device in one

Technical field

The present invention relates to remote sensing fields, particularly relate to infrared strong absorbing path image simulation method and device in one.

Background technology

Remote sensing image simulation refers to on the analysis and understanding basis of remote sensor imaging process, is obtained the process of image, thus obtain the technology of simulation remote sensing images by digital model simulation remote sensor.Remote sensing image simulation can be implemented in virtual environment and builds and test novel remote sensor, proves the technical feasibility of the earth observation plan based on novel remote sensor, thus reduces implement general plan and cost, shorten the satellite lead time.+

And the strong absorbing path of air existed in middle infra-red range, because atmospheric transmittance is low, so the remote sensor observed these passages is less, so lack relevant measured data when carrying out the remote sensing image simulation for strong absorbing path, therefore, difficulty is caused to the remote sensing image simulation of the strong absorbing path in middle infra-red range.

Summary of the invention

In view of this, embodiments provide infrared strong absorbing path image simulation method and device in one, object is the problem of the remote sensing image simulation difficulty of the strong absorbing path in solution in infra-red range.

To achieve these goals, the embodiment of the invention discloses following technical scheme:

In one, infrared strong absorbing path image simulation method, is applied to the first passage of remote sensor, and described first passage is the strong absorbing path in middle infra-red range, comprising:

According to the wave band response function of the wave band response function of described first passage and the second channel of described remote sensor, select target wavelength from the wavelength of sample library of spectra, and calculate response function corresponding to each target wavelength respectively, described second channel is non-strong absorbing path, and the spectroscopic data in spectroscopic data in wavelength band corresponding to the described second channel wavelength band corresponding with described first passage has linear relationship, described sample library of spectra comprises the spectroscopic data of N number of atural object, and wherein N is integer;

The spectroscopic data of the response function corresponding according to described each target wavelength and described sample library of spectra, determine the equivalent spectrum data of the equivalent spectrum data of M atural object of described first passage and M atural object of described second channel respectively, wherein, described M is integer, and is less than or equal to N;

According to the equivalent spectrum data of the equivalent spectrum data of M atural object of described first passage and M atural object of described second channel, determine the linearly dependent coefficient of described first passage and described second channel;

According to described linearly dependent coefficient and by the Remote Spectra image that described second channel obtains, simulate the remote sensing images obtained by described first passage, described remote sensing images comprise emissivity image or albedo image.

Preferably, described according to the wave band response function of described first passage and the wave band response function of second channel, from the wavelength of sample library of spectra, select target wavelength comprises:

According to the wave band response function of described first passage and the wave band response function of second channel, determine the union of the range of wavelengths of described first passage and described second channel;

Described and concentrated wavelength will be positioned at as target wavelength in sample library of spectra.

Preferably, the response function that each target wavelength of described calculating is corresponding comprises:

Travel through each target wavelength in accordance with the following methods, until obtain response function corresponding to each target wavelength:

From described and concentratedly to determine and the immediate first wave length of current goal wavelength and second wave length, make described target wavelength be more than or equal to described first wave length and be less than or equal to described second wave length;

According to the response function of described first wave length and the response function of described second wave length, according to the first rule, determine the response function of described current goal wavelength.

Preferably, described first rule comprises:

$f\left(j\right)+\frac{f(j+1)-f\left(j\right)}{{\mathrm{\λ}}_F(j+1)-{\mathrm{\λ}}_F\left(j\right)}*({\mathrm{\λ}}_B\left(i\right)-{\mathrm{\λ}}_F\left(j\right)),$ Wherein, described λ _bi () is described current goal wavelength, described λ _fj () is described first wave length, described λ _f(j+1) be described second wave length, described f (j) is the response function of described first wave length, the response function that described f (j+1) is described second wave length.

Preferably, the spectroscopic data of the response function corresponding according to described each target wavelength and described sample library of spectra, determine that the equivalent spectrum data of the equivalent spectrum data of M atural object of described first passage and M atural object of described second channel comprise respectively:

For each atural object, according to the response function of the wavelength in described sample library of spectra in the range of wavelengths of described first passage, according to Second Rule, determine the equivalent spectrum data of described first passage, travel through M atural object, until obtain the equivalent spectrum data of M atural object of described first passage;

For each atural object, according to the response function of the wavelength in described Sample Storehouse in the range of wavelengths of described second channel, according to described Second Rule, determine the equivalent spectrum data of described second channel, travel through M atural object, until obtain the equivalent spectrum data of M atural object of described second channel.

Preferably, described Second Rule comprises:

wherein, described in be the equivalent reflectivity (or equivalent emissivity) of first (or second channel), described x _ifor the emissivity (or reflectivity) that i-th target wavelength in the target wavelength in described sample library of spectra is corresponding, described R _ifor the response function that the target wavelength in described sample library of spectra is corresponding, described n is the number of the target wavelength in the range of wavelengths of first passage (or second channel).

Preferably, the equivalent spectrum data of the equivalent spectrum data of the described M according to a described first passage atural object and M atural object of described second channel, determine that the linearly dependent coefficient of described first passage and described second channel comprises:

Foundation $a=\frac{M\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{\mathrm{\ϵ}}_d\left(i\right)*{\mathrm{\ϵ}}_t\left(i\right)-\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{\mathrm{\ϵ}}_d\left(i\right)*\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{\mathrm{\ϵ}}_t\left(i\right)}{M\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{\mathrm{\ϵ}}_d^2\left(i\right)-{\left(\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{\mathrm{\ϵ}}_d\left(i\right)\right)}^2}$ Calculate linear multiplication coefficient;

Foundation calculate side-play amount, wherein ε _dthe equivalent spectrum data of i i-th atural object that () is second channel, ε _tthe equivalent spectrum data of i i-th atural object that () is first passage, a is linear multiplication coefficient, and b is side-play amount.

In one, infrared strong absorbing path image simulation device, is applied to the first passage of remote sensor, and described first passage is the strong absorbing path in middle infra-red range, comprising:

Response function computing module, for the wave band response function according to the wave band response function of described first passage and the second channel of described remote sensor, select target wavelength from the wavelength of sample library of spectra, and calculate response function corresponding to each target wavelength respectively, described second channel is non-strong absorbing path, and the spectroscopic data in spectroscopic data in wavelength band corresponding to the described second channel wavelength band corresponding with described first passage has linear relationship, described sample library of spectra comprises the spectroscopic data of N number of atural object, and wherein N is integer;

Equivalent spectrum data determining module, for the spectroscopic data according to response function corresponding to described each target wavelength and described sample library of spectra, determine the equivalent spectrum data of the equivalent spectrum data of M atural object of described first passage and M atural object of described second channel respectively, wherein, described M is integer, and is less than or equal to N;

Linearly dependent coefficient determination module, for the equivalent spectrum data of the equivalent spectrum data of M the atural object according to described first passage and M atural object of described second channel, determines the linearly dependent coefficient of described first passage and described second channel;

Remote sensing image simulation module, for according to described linearly dependent coefficient and the remote sensing images that obtained by described second channel, simulate the remote sensing images obtained by described first passage, described remote sensing images comprise emissivity image or albedo image.

Preferably, described response function computing module comprises:

Wavelength selection unit, for according to the wave band response function of described first passage and the wave band response function of second channel, determine the union of the range of wavelengths of described first passage and described second channel, and described and concentrated wavelength will be positioned in sample library of spectra as target wavelength;

Computing unit, for from described and concentratedly to determine and the immediate first wave length of current goal wavelength and second wave length, make described target wavelength be more than or equal to described first wave length and be less than or equal to described second wave length, according to the response function of described first wave length and the response function of described second wave length, according to the first rule, determine the response function of described current goal wavelength;

Control module, travels through each target wavelength for controlling described computing unit, until obtain response function corresponding to each target wavelength.

Preferably, described equivalent spectrum data determining module comprises:

First computing unit, for the response function according to the wavelength in described sample library of spectra in the range of wavelengths of described first passage, according to Second Rule, determines the equivalent spectrum data of arbitrary atural object of described first passage;

Second computing unit, for the response function according to the wavelength in described Sample Storehouse in the range of wavelengths of described second channel, according to described Second Rule, determines the equivalent spectrum data of arbitrary atural object of described second channel;

Control module, for controlling described first computing unit and described second unit traversal M atural object, until obtain the equivalent spectrum data of M atural object of described first passage and described second channel.

Infrared strong absorbing path image image simulation method and device during the embodiment of the present invention provides, spectroscopic data in the wavelength band that spectroscopic data in the wavelength band that second channel is corresponding is corresponding with first passage has linear relationship, therefore, by the wave band response function of first passage and the wave band response function of second channel, select target wavelength from sample library of spectra, and calculate response function corresponding to each target wavelength, and then determine the equivalent spectrum data of first passage and the equivalent spectrum data of second channel, finally can determine the linearly dependent coefficient of first passage and second channel, thus, known linear related coefficient and obtained by second channel remote sensing images time, the remote sensing images obtained by first passage can be simulated, described remote sensing images comprise emissivity image or albedo image, therefore, method described in the present embodiment and device, the remote sensing images of non-strong absorbing path can be converted to the remote sensing images of strong absorbing path, because the remote sensing images of non-strong absorbing path relatively easily obtain, therefore, facilitate the simulation of the remote sensing images to strong absorbing path.

Accompanying drawing explanation

In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.

The process flow diagram of infrared strong absorbing path image simulation method during Fig. 1 is a kind of disclosed in the embodiment of the present invention;

Fig. 2 is the process flow diagram of the disclosed infrared strong absorbing path image simulation method in another of the embodiment of the present invention;

The structural representation of infrared strong absorbing path image simulation device during Fig. 3 is a kind of disclosed in the embodiment of the present invention.

Embodiment

The embodiment of the invention discloses infrared strong absorbing path image simulation method and device in one, its central inventive point is, linear relationship between the spectroscopic data in the wavelength band that the spectroscopic data in the wavelength band that in utilization, infrared strong absorption first passage is corresponding is corresponding with non-strong absorption second channel, by the equivalent spectrum data of first passage and second channel, obtain the linearly dependent coefficient of first passage and second channel, thus according to passing through acquisition remote sensing images and the linearly dependent coefficient of second channel, simulate the remote sensing images that can be obtained by first passage, described remote sensing images comprise emissivity image or albedo image.

Below in conjunction with the accompanying drawing in the embodiment of the present invention, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.

Infrared strong absorbing path image simulation method in one disclosed in the embodiment of the present invention, the first passage of application remote sensor, described first passage is the strong absorbing path in middle infra-red range, described passage, refers to the wavelength band of sensor imaging, the image of certain passage, refer to the response element of this passage by sensor, the image only the spoke luminosity response of the wavelength band of this passage formed, method described in the present embodiment as shown in Figure 1, comprises the following steps:

S101: according to the wave band response function of the wave band response function of described first passage and the second channel of described remote sensor, select target wavelength from the wavelength of sample library of spectra, and calculate response function corresponding to each target wavelength respectively;

Described second channel is non-strong absorbing path, and the spectroscopic data in spectroscopic data in the wavelength band of its correspondence wavelength band corresponding with described first passage has linear relationship, and described sample library of spectra comprises the spectroscopic data of N number of atural object.

Usually, first passage and the wave band response function of second channel known.

Sample library of spectra for the library of spectra pre-set, wherein can comprise the spectroscopic data of N number of atural object usually, and described spectroscopic data can be emissivity data, also can be reflectivity data.Usually, sample library of spectra is the two-dimensional array of the capable W row of N, and wherein, often row represents a kind of atural object, and be often classified as emissivity data or the reflectivity data at a wavelength place, wavelength and spectroscopic data one_to_one corresponding, wherein N, W are integer.

S102: the spectroscopic data of the response function corresponding according to described each target wavelength and described sample library of spectra, determines the equivalent spectrum data of the equivalent spectrum data of M atural object of described first passage and M atural object of described second channel respectively;

Wherein, described M is integer, and is less than or equal to N;

That is, in this step, the whole atural objects in sample library of spectra can be used, a part also can be selected to use.

S103: according to the equivalent spectrum data of the equivalent spectrum data of M atural object of described first passage and M atural object of described second channel, determine the linearly dependent coefficient of described first passage and described second channel;

S104: according to described linearly dependent coefficient and the remote sensing images that obtained by described second channel, simulate the remote sensing images obtained by described first passage, described remote sensing images comprise emissivity image or albedo image.

In prior art, due to the analog data source of infrared strong absorbing path in shortage, therefore, the simulation of the remote sensing images of the infrared strong absorbing path of centering is caused to cause difficulty.

And in the present embodiment, make use of non-strong absorbing path and in linear relationship between infrared strong absorbing path, by wave band response function and the sample library of spectra of infrared strong absorbing path in known and non-strong absorbing path, determine the linearly dependent coefficient between infrared strong absorbing path and non-strong absorbing path, final according to described linearly dependent coefficient, the remote sensing images of infrared strong absorbing path during the remote sensing images that non-strong absorbing path obtains are converted to, thus realize the simulation of the remote sensing images of the infrared strong absorbing path of centering, with in compared with infrared strong absorbing path, the remote sensing images of non-strong absorbing path easily obtain, therefore, method described in the present embodiment be convenient to obtain in the remote sensing images of infrared strong absorbing path.

The disclosed infrared strong absorbing path image simulation method in another of the embodiment of the present invention, be applied to the first passage of remote sensor, described first passage is the strong absorbing path in middle infra-red range, as shown in Figure 2, comprising:

S201: according to the wave band response function of described first passage and the wave band response function of second channel, determine the union of the range of wavelengths of described first passage and described second channel;

In practical application, the passage of 3.9-4.5 mum wavelength scope has linear dependence usually, such as, 23 passages of MODIS sensor, wavelength band is 4.020-4.080 micron, and the wavelength band needing the passage of simulation is 4.1-4.3 micron, and both have linear dependence, then can using 23 passages as second channel, using passage to be simulated as first passage.

Usually, the wave band response function of passage is the two-dimensional array that 2 row L arrange, and wherein, a behavior wavelength, the response that another behavior is corresponding with wavelength, L is the number of the wavelength that wave band response function comprises.

Therefore, according to the wave band response function of first passage and the wave band response function of second channel, the range of wavelengths of first passage and second channel can be determined, further, determine the union of range of wavelengths.

S202: described and concentrated wavelength will be positioned in sample library of spectra as target wavelength;

That is, the wavelength of the wavelength coverage of first passage and second channel will be positioned in sample library of spectra as target wavelength.

S203: calculate the response function that each target wavelength is corresponding;

Namely each target wavelength is traveled through in accordance with the following methods, until obtain response function corresponding to each target wavelength:

From described and concentratedly to determine and current goal wavelength X _b(i) immediate first wave length λ _f(j) and second wave length λ _f(j+1), make described target wavelength be more than or equal to described first wave length and be less than or equal to described second wave length, be i.e. λ _f(j)≤λ _b(i)≤λ _f(j+1);

According to response function f (j) of described first wave length and the response function f (j+1) of described second wave length, according to the first rule, determine the response function of described current goal wavelength;

Described first rule is as the formula (1):

$f\left(j\right)+\frac{f(j+1)-f\left(j\right)}{{\mathrm{\λ}}_F(j+1)-{\mathrm{\λ}}_F\left(j\right)}*({\mathrm{\λ}}_B\left(i\right)-{\mathrm{\λ}}_F\left(j\right))---\left(1\right)$

Namely the wave band response function of foundation first passage and second channel, determines the response function that target wavelength is corresponding.

Such as, sample library of spectra is as shown in table 1:

Table 1

Wavelength (unit: micron)	Reflectivity
		4	0.362148
4.02	0.362944
		4.04	0.362174
4.06	0.35888
		4.08	0.356122
4.1	0.357292
		4.12	0.351872
4.14	0.350434
		4.16	0.346154
4.18	0.338216
		4.2	0.332216
4.22	0.324168
		4.24	0.314548
4.26	0.307466
		4.28	0.30641
4.3	0.307656
		4.32	0.30565
4.34	0.300018
		4.36	0.294546
4.38	0.288622
		4.4	0.281962

First passage spectral response functions is as shown in table 2:

Table 2

Wavelength (unit: micron)	Reflectivity
		4.208	0.96
4.217	0.97
		4.229	0.98
4.241	0.99
		4.251	0.99
4.259	0.99
		4.271	0.99
4.282	0.99
		4.289	0.99
4.3	0.99
		4.311	0.98
4.321	0.97
		4.329	0.96
4.338	0.95

4.351	0.94
		4.361	0.93

Second channel spectral response functions is as shown in table 3:

Table 3

Wavelength (unit: micron)	Reflectivity
		4.0183	0.6458
4.0218	0.750556
		4.0253	0.817849
4.0288	0.896326
		4.0323	0.960123
4.0358	0.97399
		4.0393	0.945746
4.0428	0.959025
		4.0463	0.943537
4.0498	0.99655
		4.0533	0.937629
4.0568	0.925526
		4.0603	0.920733
4.0638	0.951187
		4.0673	0.948747
4.0708	0.96493
		4.0743	0.98478
4.0778	0.998212
		4.0813	0.997535
4.0848	0.975225
		4.0883	0.925132
4.0918	0.838034
		4.0953	0.716928
4.0988	0.58802

The range of wavelengths of first passage and second channel is respectively 4.2-4.36 and 4.020-4.080, then the wavelength can selecting the range of wavelengths being positioned at first passage from table 1 is 4.2,4.22,4.24,4.26,4.28,4.3,4.32,4.34,4.36, the wavelength being positioned at the range of wavelengths of second channel is 4.02,4.04,4.06,4.08, then in sample library of spectra, above-mentioned wavelength is target wavelength.

For any one target wavelength, the wavelength 4.22 of such as first passage, select wavelength 4.217 and 4.229 immediate with it in the wavelength that first passage response function comprises, according to the response function of 4.217 and 4.229, determine the response function of current goal wavelength.The wavelength 4.06 of such as second channel again, selects wavelength 4.0568 and 4.0603 immediate with it in the wavelength that second channel response function comprises, and according to the response function of 4.0568 and 4.0603, determines the response function of current goal wavelength.

It should be noted that first passage is different with the range of wavelengths of second channel.So any one target wavelength, corresponding wave band may be only first or second channel both one, instead of to be present within the scope of the first and second channel wavelengths simultaneously.Here, illustrate for second channel wavelength.

S204: the equivalent spectrum data calculating first passage and second channel respectively:

For each atural object, according to the response function of the wavelength in described sample library of spectra in the range of wavelengths of described first passage, according to Second Rule, determine the equivalent spectrum data of described first passage, travel through M atural object, until obtain the equivalent spectrum data of M atural object of described first passage;

For each atural object, according to the response function of the wavelength in described Sample Storehouse in the range of wavelengths of described second channel, according to described Second Rule, determine the equivalent spectrum data of described second channel, travel through M atural object, until obtain the equivalent spectrum data of M atural object of described second channel.

That is, the equivalent spectrum data of the response function determination first passage utilizing the target wavelength in the range of wavelengths of first passage corresponding, the equivalent spectrum data of the response function determination second channel utilizing again the target wavelength of second channel corresponding.

In the present embodiment, Second Rule as the formula (2):

$\stackrel{\‾}{y}=\underset{i=1}{\overset{i=n}{\mathrm{\Σ}}}{x}_i{R}_i/\underset{i=1}{\overset{i=n}{\mathrm{\Σ}}}{R}_i,---\left(2\right)$

Wherein, be the equivalent reflectivity (or equivalent emissivity) of first (or second channel), x _ifor the emissivity (or reflectivity) that i-th target wavelength in the target wavelength in sample library of spectra is corresponding, R _ifor the response function that i-th target wavelength in the target wavelength in sample library of spectra is corresponding, n is the number of the target wavelength in the range of wavelengths of first passage (or second channel).

S205: the equivalent spectrum data of foundation first passage and second channel, utilize formula (3) and formula (4) to determine the linearly dependent coefficient of described first passage and described second channel:

$a=\frac{M\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{\mathrm{\ϵ}}_d\left(i\right)*{\mathrm{\ϵ}}_t\left(i\right)-\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{\mathrm{\ϵ}}_d\left(i\right)*\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{\mathrm{\ϵ}}_t\left(i\right)}{M\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{\mathrm{\ϵ}}_d^2\left(i\right)-{\left(\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{\mathrm{\ϵ}}_d\left(i\right)\right)}^2}---\left(3\right)$

$b=\frac{\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{\mathrm{\ϵ}}_t\left(i\right)-a*\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{\mathrm{\ϵ}}_d\left(i\right)}{M}---\left(4\right)$

Wherein, ε _dthe equivalent emissivity of i i-th atural object that () is second channel, ε _tthe equivalent emissivity of i i-th atural object that () is first passage, a is linear multiplication coefficient, and b is side-play amount.It should be noted that, above formula is for equivalent emissivity, if sample library of spectra provides reflectivity data, also equivalent emissivity can be replaced with equivalent reflectivity, application of formula (3) and (4) calculate linear multiplication coefficient and side-play amount.

S206: according to described linearly dependent coefficient and the remote sensing images that obtained by described second channel, simulate the remote sensing images obtained by described first passage, described remote sensing images comprise emissivity image or albedo image.

Particularly, suppose that the remote sensing images obtained by second channel are D, then the remote sensing images T obtained by first passage is:

T=a*D+b。

Remote sensing image simulation method described in the present embodiment, as as described in above steps, calculate the linearly dependent coefficient between first passage and second channel, according to described linearly dependent coefficient, the remote sensing images of second channel are converted to the remote sensing images of first passage, because second channel is non-strong absorbing path, so, its remote sensing images easily obtain, and therefore, are convenient to the remote sensing images obtaining first passage.

Corresponding with said method embodiment, the embodiment of the invention also discloses infrared strong absorbing path image simulation device in one, as shown in Figure 3, comprising:

Response function computing module 301, for the wave band response function according to the wave band response function of described first passage and the second channel of described remote sensor, select target wavelength from the wavelength of sample library of spectra, and calculate response function corresponding to each target wavelength respectively, described second channel is non-strong absorbing path, and the spectroscopic data in spectroscopic data in wavelength band corresponding to the described second channel wavelength band corresponding with described first passage has linear relationship, described sample library of spectra comprises the spectroscopic data of N number of atural object, and wherein N is integer;

Equivalent spectrum data determining module 302, for the spectroscopic data according to response function corresponding to described each target wavelength and described sample library of spectra, determine the equivalent spectrum data of the equivalent spectrum data of M atural object of described first passage and M atural object of described second channel respectively, wherein, described M is integer, and is less than or equal to N;

Linearly dependent coefficient determination module 303, for the equivalent spectrum data of the equivalent spectrum data of M the atural object according to described first passage and M atural object of described second channel, determines the linearly dependent coefficient of described first passage and described second channel;

Remote sensing image simulation module 304, for according to described linearly dependent coefficient and the remote sensing images that obtained by described second channel, simulate the remote sensing images obtained by described first passage, described remote sensing images comprise emissivity image or albedo image.

In the present embodiment, the concrete grammar of the remote sensing images that calculating linearly dependent coefficient and simulation are obtained by described first passage is as described in above-described embodiment.

Further, described in the present embodiment, response function computing module comprises:

Wavelength selection unit, for according to the wave band response function of described first passage and the wave band response function of second channel, determine the union of the range of wavelengths of described first passage and described second channel, and described and concentrated wavelength will be positioned in sample library of spectra as target wavelength;

Computing unit, for from described and concentratedly to determine and the immediate first wave length of current goal wavelength and second wave length, make described target wavelength be more than or equal to described first wave length and be less than or equal to described second wave length, according to the response function of described first wave length and the response function of described second wave length, according to the first rule, determine the response function of described current goal wavelength;

Control module, travels through each target wavelength for controlling described computing unit, until obtain response function corresponding to each target wavelength.

Wherein, the first rule, Ru shown in above-mentioned formula (1), repeats no more here.

Further, described in the present embodiment, equivalent spectrum data determining module comprises:

First computing unit, for the response function according to the wavelength in described sample library of spectra in the range of wavelengths of described first passage, according to Second Rule, determines the equivalent spectrum data of arbitrary atural object of described first passage;

Second computing unit, for the response function according to the wavelength in described Sample Storehouse in the range of wavelengths of described second channel, according to described Second Rule, determines the equivalent spectrum data of arbitrary atural object of described second channel;

Control module, for controlling described first computing unit and described second computing unit traversal M atural object, until obtain the equivalent spectrum data of M atural object of described first passage and described second channel.

Wherein, Second Rule, as shown in above-mentioned formula (2), repeats no more here.

If the function described in the present embodiment method using the form of SFU software functional unit realize and as independently production marketing or use time, can be stored in a computing equipment read/write memory medium.Based on such understanding, the part of the part that the embodiment of the present invention contributes to prior art or this technical scheme can embody with the form of software product, this software product is stored in a storage medium, comprising some instructions in order to make a computing equipment (can be personal computer, server, mobile computing device or the network equipment etc.) perform all or part of step of method described in each embodiment of the present invention.And aforesaid storage medium comprises: USB flash disk, portable hard drive, ROM (read-only memory) (ROM, Read-OnlyMemory), random access memory (RAM, RandomAccessMemory), magnetic disc or CD etc. various can be program code stored medium.

In this instructions, each embodiment adopts the mode of going forward one by one to describe, and what each embodiment stressed is the difference with other embodiment, between each embodiment same or similar part mutually see.

To the above-mentioned explanation of the disclosed embodiments, professional and technical personnel in the field are realized or uses the present invention.To be apparent for those skilled in the art to the multiple amendment of these embodiments, General Principle as defined herein can without departing from the spirit or scope of the present invention, realize in other embodiments.Therefore, the present invention can not be restricted to these embodiments shown in this article, but will meet the widest scope consistent with principle disclosed herein and features of novelty.

Claims (2)

1. in, an infrared strong absorbing path image simulation method, is applied to the first passage of remote sensor, and described first passage is the strong absorbing path in middle infra-red range, it is characterized in that, comprising:

According to the wave band response function of the wave band response function of described first passage and the second channel of described remote sensor, select target wavelength from the wavelength of sample library of spectra, and calculate response function corresponding to each target wavelength respectively, described second channel is non-strong absorbing path, and the spectroscopic data in spectroscopic data in wavelength band corresponding to the described second channel wavelength band corresponding with described first passage has linear relationship, described sample library of spectra comprises the spectroscopic data of N number of atural object, and wherein N is integer;

The spectroscopic data of the response function corresponding according to described each target wavelength and described sample library of spectra, determine the equivalent spectrum data of the equivalent spectrum data of M atural object of described first passage and M atural object of described second channel respectively, wherein, described M is integer, and is less than or equal to N;

According to the equivalent spectrum data of the equivalent spectrum data of M atural object of described first passage and M atural object of described second channel, determine the linearly dependent coefficient of described first passage and described second channel;

According to described linearly dependent coefficient and the remote sensing images that obtained by described second channel, simulate the remote sensing images obtained by described first passage, described remote sensing images comprise emissivity image or albedo image;

Described according to the wave band response function of described first passage and the wave band response function of second channel, from the wavelength of sample library of spectra, select target wavelength comprises:

According to the wave band response function of described first passage and the wave band response function of second channel, determine the union of the range of wavelengths of described first passage and described second channel;

Described and concentrated wavelength will be positioned at as target wavelength in sample library of spectra;

Response function corresponding to each target wavelength of described calculating comprises:

Travel through each target wavelength in accordance with the following methods, until obtain response function corresponding to each target wavelength:

From described and concentratedly to determine and the immediate first wave length of current goal wavelength and second wave length, make described target wavelength be more than or equal to described first wave length and be less than or equal to described second wave length;

According to the response function of described first wave length and the response function of described second wave length, according to the first rule, determine the response function of described current goal wavelength;

Described first rule comprises:

wherein, described λ _bi () is described current goal wavelength, described λ _fj () is described first wave length, described λ _f(j+1) be described second wave length, described f (j) is the response function of described first wave length, the response function that described f (j+1) is described second wave length;

The spectroscopic data of the described response function corresponding according to described each target wavelength and described sample library of spectra, determine that the equivalent spectrum data of the equivalent spectrum data of M atural object of described first passage and M atural object of described second channel comprise respectively:

For each atural object, according to the response function of the wavelength in described sample library of spectra in the range of wavelengths of described first passage, according to Second Rule, determine the equivalent spectrum data of described first passage, travel through M atural object, until obtain the equivalent spectrum data of M atural object of described first passage;

For each atural object, according to the response function of the wavelength in described Sample Storehouse in the range of wavelengths of described second channel, according to described Second Rule, determine the equivalent spectrum data of described second channel, travel through M atural object, until obtain the equivalent spectrum data of M atural object of described second channel;

Described Second Rule comprises:

wherein, described in be first or the equivalent reflectivity of second channel or equivalent emissivity, described x _ithe emissivity corresponding for i-th target wavelength in the target wavelength in described sample library of spectra or reflectivity, described R _ifor the response function that the target wavelength in described sample library of spectra is corresponding, described n is the number of the target wavelength in the range of wavelengths of first passage or second channel;

The equivalent spectrum data of the equivalent spectrum data of the described M according to a described first passage atural object and M atural object of described second channel, determine that the linearly dependent coefficient of described first passage and described second channel comprises:

Foundation $a=\frac{M\underset{i=1}{\overset{M}{\Σ}}{\mathrm{\ϵ}}_d\left(i\right)*{\mathrm{\ϵ}}_t\left(i\right)-\underset{i=1}{\overset{M}{\Σ}}{\mathrm{\ϵ}}_d\left(i\right)*\underset{i=1}{\overset{M}{\Σ}}{\mathrm{\ϵ}}_t\left(i\right)}{M\underset{i=1}{\overset{M}{\Σ}}{\mathrm{\ϵ}}_d^2\left(i\right)-{\left(\underset{i=1}{\overset{M}{\Σ}}{\mathrm{\ϵ}}_d\left(i\right)\right)}^2}$ Calculate linear multiplication coefficient;

Foundation calculate side-play amount, wherein ε _dthe equivalent spectrum data of i i-th atural object that () is second channel, ε _tthe equivalent spectrum data of i i-th atural object that () is first passage, a is linear multiplication coefficient, and b is side-play amount.

2. in, an infrared strong absorbing path image simulation device, is applied to the first passage of remote sensor, and described first passage is the strong absorbing path in middle infra-red range, it is characterized in that, comprising:

Response function computing module, for the wave band response function according to the wave band response function of described first passage and the second channel of described remote sensor, select target wavelength from the wavelength of sample library of spectra, and calculate response function corresponding to each target wavelength respectively, described second channel is non-strong absorbing path, and the spectroscopic data in spectroscopic data in wavelength band corresponding to the described second channel wavelength band corresponding with described first passage has linear relationship, described sample library of spectra comprises the spectroscopic data of N number of atural object, and wherein N is integer;

Equivalent spectrum data determining module, for the spectroscopic data according to response function corresponding to described each target wavelength and described sample library of spectra, determine the equivalent spectrum data of the equivalent spectrum data of M atural object of described first passage and M atural object of described second channel respectively, wherein, described M is integer, and is less than or equal to N;

Linearly dependent coefficient determination module, for the equivalent spectrum data of the equivalent spectrum data of M the atural object according to described first passage and M atural object of described second channel, determines the linearly dependent coefficient of described first passage and described second channel;

Remote sensing image simulation module, for according to described linearly dependent coefficient and the remote sensing images that obtained by described second channel, simulate the remote sensing images obtained by described first passage, described remote sensing images comprise emissivity image or albedo image;

Described response function computing module comprises:

Wavelength selection unit, for according to the wave band response function of described first passage and the wave band response function of second channel, determine the union of the range of wavelengths of described first passage and described second channel, and described and concentrated wavelength will be positioned in sample library of spectra as target wavelength;

Computing unit, for from described and concentratedly to determine and the immediate first wave length of current goal wavelength and second wave length, make described target wavelength be more than or equal to described first wave length and be less than or equal to described second wave length, according to the response function of described first wave length and the response function of described second wave length, according to the first rule, determine the response function of described current goal wavelength;

Control module, travels through each target wavelength for controlling described computing unit, until obtain response function corresponding to each target wavelength;

Described equivalent spectrum data determining module comprises:

First computing unit, for the response function according to the wavelength in described sample library of spectra in the range of wavelengths of described first passage, according to Second Rule, determines the equivalent spectrum data of arbitrary atural object of described first passage;

Second computing unit, for the response function according to the wavelength in described Sample Storehouse in the range of wavelengths of described second channel, according to described Second Rule, determines the equivalent spectrum data of arbitrary atural object of described second channel;

Control module, for controlling described first computing unit and described second unit traversal M atural object, until obtain the equivalent spectrum data of M atural object of described first passage and described second channel.