CN105737806B - Method and apparatus based on satellite-borne microwave radiometer inverting sea surface salinity - Google Patents
Method and apparatus based on satellite-borne microwave radiometer inverting sea surface salinity Download PDFInfo
- Publication number
- CN105737806B CN105737806B CN201610087880.XA CN201610087880A CN105737806B CN 105737806 B CN105737806 B CN 105737806B CN 201610087880 A CN201610087880 A CN 201610087880A CN 105737806 B CN105737806 B CN 105737806B
- Authority
- CN
- China
- Prior art keywords
- difference
- wave band
- vertical polarization
- observation
- reflectivity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C13/00—Surveying specially adapted to open water, e.g. sea, lake, river or canal
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Hydrology & Water Resources (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radiation Pyrometers (AREA)
- Image Processing (AREA)
Abstract
The present invention provides a kind of method and apparatus based on satellite-borne microwave radiometer inverting sea surface salinity, are related to marine exploration field, including：Obtain the observation bright temperature data auxiliary data corresponding with observation bright temperature data of target water area wave band c and x wave bands, and the difference of calculating observation c/x wave band vertical polarization reflectivity respectively by satellite-borne microwave radiometer；Difference to observing c/x wave band vertical polarization reflectivity carries out digital picture optimization processing, the difference of the observation c/x wave band vertical polarization reflectivity optimized；Inversion Calculation is carried out to the difference and auxiliary data of the observation c/x wave band vertical polarization reflectivity of optimization using default look-up table, obtain the sea surface salinity of target water area, it optimizes processing to the result of calculation for observing bright temperature data, make result of calculation closer to actual value, but also the distribution situation of the sea surface salinity obtained by the calculating structure inversion is closer to actual value, to improve the inversion accuracy of sea surface salinity.
Description
Technical field
The present invention relates to marine exploration fields, and satellite-borne microwave radiometer inverting sea salt is based in particular to one kind
The method and apparatus of degree.
Background technology
Sea surface salinity is to describe one of the important parameter of sea state, studies its variation and the regularity of distribution to analyzing ocean certainly
The effect of body characteristic and understanding ocean in sea-gas complication system is significant.Also, sea surface salinity be also meteorology,
The research object that the Other subjects such as ecological, hydrology and fishery are paid close attention to application field.
During studying sea surface salinity, it is necessary first to obtain sea surface salinity；Currently, obtaining sea surface salinity
Method is, using in orbit and producing the microwave radiometer of business product and obtaining corresponding data, then directly to obtaining
Data handled accordingly, to be finally inversed by corresponding sea surface salinity；Wherein, above-mentioned microwave radiometer can be Japan or
The AMSR series microwave radiometers in the U.S..
Inventor has found under study for action, is directly made in the prior art by the method for the data inversion sea surface salinity of acquisition
The obtained salt angle value of inverting and the deviation of practical salt angle value it is larger.
Invention content
The purpose of the present invention is to provide a kind of method and apparatus based on satellite-borne microwave radiometer inverting sea surface salinity, energy
The enough result of calculation by observation bright temperature data optimizes processing, to improve the inversion accuracy of sea surface salinity value.
In a first aspect, an embodiment of the present invention provides a kind of method based on satellite-borne microwave radiometer inverting sea surface salinity,
Method includes：
Obtain observation bright temperature data and the sight of target water area wave band c and x wave bands respectively by satellite-borne microwave radiometer
Survey the corresponding auxiliary data of bright temperature data；
According to the observation bright temperature data of c/x wave bands and corresponding auxiliary data, the reflection of calculating observation c/x wave band vertical polarizations
The difference of rate；
Difference to observing c/x wave band vertical polarization reflectivity carries out digital picture optimization processing, the observation c/x optimized
The difference of wave band vertical polarization reflectivity；
Inverting is carried out to the difference and auxiliary data of the observation c/x wave band vertical polarization reflectivity of optimization using default look-up table
It calculates, obtains the sea surface salinity of target water area.
With reference to first aspect, an embodiment of the present invention provides the first possible embodiments of first aspect, wherein right
The difference for observing c/x wave band vertical polarization reflectivity carries out digital picture optimization processing, the vertical pole of observation c/x wave bands optimized
Change reflectivity difference include：
According to dielectric constant model and auxiliary data, the difference of calculating simulation c/x wave band vertical polarization reflectivity；
The difference and the difference of simulation c/x wave band vertical polarization reflectivity of observing c/x wave band vertical polarization reflectivity are carried out pair
Than statistical disposition, obtain include observe c/x wave band vertical polarization reflectivity difference and simulation c/x wave band vertical polarizations reflectivity it
The two-dimensional digital image of the order of magnitude of difference；
Adjust the order of magnitude of the difference of observation c/x wave band vertical polarization reflectivity and simulation c/x wave bands in two-dimensional digital image
The order-of-magnitude agreement of the difference of vertical polarization reflectivity obtains revised two-dimensional digital image for the first time；
Second of correcting process is carried out to revised two-dimensional digital image for the first time using data rotation method, is obtained excellent
The difference of the observation c/x wave band vertical polarization reflectivity of change.
The possible embodiment of with reference to first aspect the first, an embodiment of the present invention provides second of first aspect
Possible embodiment, wherein adjustment two-dimensional digital image in observation c/x wave band vertical polarization reflectivity difference the order of magnitude with
The order-of-magnitude agreement of difference for simulating c/x wave band vertical polarization reflectivity includes：
The difference for simulating c/x wave band vertical polarization reflectivity is expressed as by observation c/ using the method for least squares line fitting
The function of the difference of x wave band vertical polarization reflectivity, obtains update equation Rmodel=0.0017Robs+0.0107；
According to update equation Rmodel=0.0017Robs+ 0.0107, it is vertical to adjust observation c/x wave bands in two-dimensional digital image
The order-of-magnitude agreement of the order of magnitude of the difference of polarization reflectivity and the difference of simulation c/x wave band vertical polarization reflectivity.
Second of possible embodiment with reference to first aspect, an embodiment of the present invention provides the third of first aspect
Possible embodiment, wherein revised two-dimensional digital image for the first time is repaiied for the second time using data rotation method
Positive processing, the difference of the observation c/x wave band vertical polarization reflectivity optimized include：
By the difference of the observation c/x wave band vertical polarization reflectivity in revised two-dimensional digital image for the first time and simulation c/
The difference of x wave band vertical polarization reflectivity carries out pairing statistical disposition, obtains the two-dimensional digital image with data value；Wherein, two
X-axis where the pictorial element of dimension word image indicates that model reflectivity, y-axis indicate observation reflectivity；
According to the center of mass point for the two dimensional image that all pictorial elements in two-dimensional digital image are constituted, two-dimensional digital figure is determined
The rotation axis of picture；
According to the APPROXIMATE DISTRIBUTION straight line of all pictorial elements in two-dimensional digital image, the rotation angle of two-dimensional digital image is determined
Degree；
According to rotation axis and rotation angle, rotation processing is carried out to revised two-dimensional digital image for the first time, is obtained
Second of revised two-dimensional digital image；
The difference of extraction observation c/x wave band vertical polarization reflectivity from second of revised two-dimensional digital image, as
The difference of the observation c/x wave band vertical polarization reflectivity of optimization.
The third possible embodiment with reference to first aspect, an embodiment of the present invention provides the 4th kind of first aspect
Possible embodiment, wherein according to the observation bright temperature data of c/x wave bands and corresponding auxiliary data, calculating observation c/x wave bands
The difference of vertical polarization reflectivity includes：
The difference for calculating the vertical polarization reflectivity of wave band c and the vertical polarization reflectivity of x wave bands obtains observation c/x waves
The difference of section vertical polarization reflectivity.
Second aspect, the embodiment of the present invention additionally provide a kind of dress based on satellite-borne microwave radiometer inverting sea surface salinity
It sets, device includes：
Acquisition module, the observation for obtaining target water area wave band c and x wave bands respectively by satellite-borne microwave radiometer
Bright temperature data auxiliary data corresponding with observation bright temperature data；
Computing module is used for the observation bright temperature data according to c/x wave bands and corresponding auxiliary data, calculating observation c/x waves
The difference of section vertical polarization reflectivity；
Digital picture optimization processing module, for carrying out digital picture to the difference for observing c/x wave band vertical polarization reflectivity
Optimization processing, the difference of the observation c/x wave band vertical polarization reflectivity optimized；
Inversion Calculation module, for the difference using default look-up table to the observation c/x wave band vertical polarization reflectivity of optimization
Inversion Calculation is carried out with auxiliary data, obtains the sea surface salinity of target water area.
In conjunction with second aspect, an embodiment of the present invention provides the first possible embodiments of second aspect, wherein number
Word image optimization processing module includes：
First computing unit, for according to dielectric constant model and auxiliary data, calculating simulation c/x wave band vertical polarizations to be anti-
Penetrate the difference of rate；
Statistical disposition unit is compared, for vertical with simulation c/x wave bands to the difference for observing c/x wave band vertical polarization reflectivity
The difference of polarization reflectivity carries out comparison statistical disposition, obtains including the difference and simulation c/x for observing c/x wave band vertical polarization reflectivity
The two-dimensional digital image of the order of magnitude of the difference of wave band vertical polarization reflectivity；
First time correcting process unit, for adjust observed in two-dimensional digital image c/x wave band vertical polarizations reflectivity it
The order-of-magnitude agreement of the order of magnitude of difference and the difference of simulation c/x wave band vertical polarization reflectivity, obtains revised two dimension for the first time
Digital picture；
Second of correcting process unit, for using data rotation method to revised two-dimensional digital image for the first time into
Second of correcting process of row, the difference of the observation c/x wave band vertical polarization reflectivity optimized.
In conjunction with the first possible embodiment of second aspect, an embodiment of the present invention provides second of second aspect
Possible embodiment, wherein first time correcting process unit includes：
Update equation obtains subelement, and the vertical pole of c/x wave bands will be simulated for the method using least squares line fitting
The difference for changing reflectivity is expressed as observing the function of the difference of c/x wave band vertical polarization reflectivity, obtains update equation Rmodel=
0.0017Robs+0.0107；
Subelement is adjusted, for according to update equation Rmodel=0.0017Robs+ 0.0107, it adjusts in two-dimensional digital image
It observes the order of magnitude of the difference of c/x wave band vertical polarization reflectivity and simulates the order of magnitude of the difference of c/x wave band vertical polarization reflectivity
Unanimously.
In conjunction with second of possible embodiment of second aspect, an embodiment of the present invention provides the third of second aspect
Possible embodiment, wherein second of correcting process unit include：
Statistical disposition subelement is matched, for the observation c/x wave bands in revised two-dimensional digital image for the first time to hang down
The difference of straight polarization reflectivity and the difference of simulation c/x wave band vertical polarization reflectivity carry out pairing statistical disposition, obtain with data
The two-dimensional digital image of value；Wherein, the x-axis where the pictorial element of two-dimensional digital image indicates that model reflectivity, y-axis indicate to see
Survey reflectivity；
Rotation axis determination subelement, the two dimensional image for being constituted according to all pictorial elements in two-dimensional digital image
Center of mass point, determine the rotation axis of two-dimensional digital image；
Rotation angle determination subelement is determined according to the APPROXIMATE DISTRIBUTION straight line of all pictorial elements in two-dimensional digital image
The rotation angle of two-dimensional digital image；
Rotation processing subelement is used for according to rotation axis and rotation angle, to revised two-dimensional digital figure for the first time
As carrying out rotation processing, second of revised two-dimensional digital image is obtained；
Subelement is extracted, for the extraction observation c/x wave band vertical polarizations from second of revised two-dimensional digital image
The difference of reflectivity, the difference of observation c/x wave band vertical polarization reflectivity as an optimization.
In conjunction with the third possible embodiment of second aspect, an embodiment of the present invention provides the 4th kind of second aspect
Possible embodiment, wherein computing module includes：
Third computing unit, the difference for calculating the vertical polarization reflectivity of wave band c and the vertical polarization reflectivity of x wave bands
Value obtains the difference of observation c/x wave band vertical polarization reflectivity.
A kind of method and apparatus based on satellite-borne microwave radiometer inverting sea surface salinity provided in an embodiment of the present invention, method
Including：Obtain the observation bright temperature data of target water area wave band c and x wave bands and corresponding respectively by satellite-borne microwave radiometer
Auxiliary data, and the difference of calculating observation c/x wave band vertical polarization reflectivity；Then, to observing c/x wave band vertical polarization reflectivity
Difference carry out digital picture optimization processing；Finally, the reflection using default look-up table to observation wave band c and x wave bands after optimization
The difference and auxiliary data of rate carry out Inversion Calculation, obtain the sea surface salinity of target water area, directly lead to in the prior art
Cross obtain data inversion sea surface salinity method so that the salt angle value that inverting obtains with practical salt angle value deviate it is larger compared with,
It optimizes processing to the result of calculation for observing bright temperature data so that result of calculation is closer to actual value, but also passing through this
The distribution situation for the sea surface salinity that structure inversion obtains is calculated closer to actual value, to improve the inverting essence of sea surface salinity
Degree.
To enable the above objects, features and advantages of the present invention to be clearer and more comprehensible, preferred embodiment cited below particularly, and coordinate
Appended attached drawing, is described in detail below.
Description of the drawings
In order to illustrate the technical solution of the embodiments of the present invention more clearly, below will be to needed in the embodiment attached
Figure is briefly described, it should be understood that the following drawings illustrates only certain embodiments of the present invention, therefore is not construed as pair
The restriction of range for those of ordinary skill in the art without creative efforts, can also be according to this
A little attached drawings obtain other relevant attached drawings.
Fig. 1 shows a kind of method based on satellite-borne microwave radiometer inverting sea surface salinity that the embodiment of the present invention is provided
Flow chart；
Fig. 2 shows another sides based on satellite-borne microwave radiometer inverting sea surface salinity that the embodiment of the present invention is provided
The flow chart of method；
Fig. 3 shows that the embodiment of the present invention includes observing the difference and simulation c/x wave bands of c/x wave band vertical polarization reflectivity
The schematic diagram of the two-dimentional dimension word image of the order of magnitude of the difference of vertical polarization reflectivity；
Fig. 4 shows the schematic diagram of the two-dimentional dimension word image after being corrected for the first time in the embodiment of the present invention；
Fig. 5 shows another side based on satellite-borne microwave radiometer inverting sea surface salinity that the embodiment of the present invention is provided
The flow chart of method；
Fig. 6 shows the schematic diagram of the two-dimensional digital image with data value in the embodiment of the present invention；
Fig. 7 shows the center-of-mass coordinate of 12 parts in the two-dimensional digital image with data value in the embodiment of the present invention
Schematic diagram；
Fig. 8 shows the schematic diagram of the two-dimensional digital image after being corrected by second in the embodiment of the present invention；
Fig. 9 shows a kind of device based on satellite-borne microwave radiometer inverting sea surface salinity that the embodiment of the present invention is provided
Structural schematic diagram；
Figure 10 shows a kind of dress based on satellite-borne microwave radiometer inverting sea surface salinity that the embodiment of the present invention is provided
Set the structural schematic diagram of middle digital picture optimization processing module and first time correcting process unit；
Figure 11 shows a kind of dress based on satellite-borne microwave radiometer inverting sea surface salinity that the embodiment of the present invention is provided
The structural schematic diagram of second correcting process unit and computing module in setting.
Specific implementation mode
Below in conjunction with attached drawing in the embodiment of the present invention, technical solution in the embodiment of the present invention carries out clear, complete
Ground describes, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Usually exist
The component of the embodiment of the present invention described and illustrated in attached drawing can be arranged and be designed with a variety of different configurations herein.Cause
This, the detailed description of the embodiment of the present invention to providing in the accompanying drawings is not intended to limit claimed invention below
Range, but it is merely representative of the selected embodiment of the present invention.Based on the embodiment of the present invention, those skilled in the art are not doing
The every other embodiment obtained under the premise of going out creative work, shall fall within the protection scope of the present invention.
For directly passing through the method for the data inversion sea surface salinity of acquisition in the prior art so that the salinity that inverting obtains
Value deviates larger problem with practical salt angle value, and an embodiment of the present invention provides one kind being based on satellite-borne microwave radiometer inverting sea
The method and apparatus of salinity optimize processing by the result of calculation to the observation bright temperature data directly acquired so that should
Result of calculation is closer to actual value, but also the distribution situation of the sea surface salinity obtained by the calculating structure inversion is closer in fact
Actual value, to improve the inversion accuracy of sea surface salinity.
Illustrate that an embodiment of the present invention provides a kind of based on satellite-borne microwave radiometer inverting sea surface salinity with reference to Fig. 1
Method, the method includes：
S101, target water area wave band c and the observation bright temperature data of x wave bands are obtained respectively by satellite-borne microwave radiometer
Auxiliary data corresponding with observation bright temperature data.
Specifically, microwave radiation from sea surface is mainly determined by factors such as sea surface temperature, sea surface salinity, the coarse situations in sea, this
Outside, microwave radiation from sea surface is also influenced by factors such as radiation frequency, observation angle and polarized states.Currently, distant using space-based
By SMOS and Aquarius completions, (SMOS and Aquarius herein is microwave radiation to the information of the means acquisition sea surface salinity of sense
Observation system), above-mentioned two observation system is all equipped with dedicated L-band and carries out radiation detection.In addition to above-mentioned two microwave spoke
It penetrates observation system, still in orbit and produces AMSR series and (the i.e. seas HY of China for having Japan/U.S. of business product
Ocean) series microwave radiometer (such as microwave radiometer of ocean two).The minimum work of sensor of the Seeds of First Post-flight of above-mentioned two series
Working frequency is located at C-band, wherein C-band and X-band are different for the sensibility of salinity, and all very low.Therefore utilize C waves
Section and the difference of X-band salinity susceptivity can be with inverting sea surface salinities；Based on this, in the embodiment of the present invention, pass through satellite-borne microwave
Radiometer obtains the observation bright temperature data of wave band c and x wave bands, to the sea surface salinity of follow-up inverting target area.
Experimental Area (i.e. above-mentioned target water area) is selected in China Seas and vicinity by the embodiment of the present invention
(Longitude：105 ° of E-135 ° of E, Latitude：0 ° of -40 ° of N), and pay close attention to Huanghai Sea East Sea sea area.This is because here
It is to change relative to thalassicum violent where Changjiang River into sea mouth.
The embodiment of the present invention is the observation bright temperature data and auxiliary that satellite-borne microwave radiometer first obtains wave band c and x wave bands
Then data carry out consistency treatment to the above three data of acquisition, that is, are chosen at the corresponding wave band c of same auxiliary data
Observe bright temperature data, the observation bright temperature data of x wave bands；Wherein, the satellite-borne microwave radiometer that the embodiment of the present invention uses it is preferred
For HY-2A microwave radiometers, i.e., first ocean dynamical environment satellite that China sends is obtained by the related data of the HY-2A
Take the observation bright temperature data of wave band c, the observation bright temperature data of x wave bands；Above-mentioned auxiliary data includes：It is compared for inversion result
WOA09 salinity Climatologicals data (being referred to as WOA09 Salinity Datas), for observing the bright matched AMSR-2 of temperature with HY-2A
L2b data.Wherein, above-mentioned l2b data, including sea surface temperature, ocean surface wind speed, moisture content, liquid water content；WOA09 salinity
Climatological data are salinity Climatological data, are provided by month.Above-mentioned auxiliary data can facilitate salinity inverting, and above-mentioned auxiliary
Help above-mentioned 5 physical quantitys (i.e. sea surface temperature, ocean surface wind speed, moisture content, liquid water content and the WOA09 salinity Climatologicals of data
Data) it the modes such as sounding balloon, buoy may be used obtains.
S102, the observation bright temperature data according to c/x wave bands and corresponding auxiliary data, the vertical pole of calculating observation c/x wave bands
Change the difference of reflectivity.
Specifically, according to formulaThen the observation bright temperature data for calculating separately c/x wave bands is
The difference of wave band c and the bright temperature data of x wave bands can be calculated；Wherein, R is the corresponding reflectivity of sea surface emissivity, R=1-E.M is one
It is a contain Downward atmospheric long-wave radiation, atmospheric transmittance, surface scattering correction, cosmic background radiation function；TBIndicate microwave spoke
Penetrate the bright temperature of meter detection, TBUIt is air uplink radiation (upwelling atmospheric brightness
Temperature), TBΩIt is surface scattering (the sky radiation scattered upward by Earth of atmospheric radiation
Surface), TsIt is sea surface temperature, τ is atmospheric transmittance (atmospheric transmission), and E is sea surface emissivity
(sea-surface emissivity)。
S103, digital picture optimization processing, the sight optimized are carried out to the difference for observing c/x wave band vertical polarization reflectivity
Survey the difference of c/x wave band vertical polarization reflectivity.
As described above, the observation bright temperature data of wave band c and x wave bands is obtained in the embodiment of the present invention using HY-2a radiometers,
In view of the observation data of above-mentioned HY-2a instruments may have systematic deviation, need to be modified the data of acquisition.
Specific correcting mode is as follows：The embodiment of the present invention：First according to dielectric constant model and auxiliary data, calculating simulation
The difference (i.e. the difference of model reflectivity) of c/x wave band vertical polarization reflectivity, then compares and observes bright temperature and supplementary number by HY-2a
According to (AMSR-2L2b products and WOA09 Salinity Datas) calculate gained C-band and X-band vertical polarization reflectivity difference and
The distribution of both differences of model reflectivity, and following data processing method is devised according to above-mentioned comparison result：The first step,
- translation processing is zoomed in and out, the order of magnitude that will calculate the difference of the observation c/x wave band vertical polarization reflectivity of gained is adjusted to and mould
The difference of type reflectivity it is consistent；Second step carries out roto-translatory processing, and the observation c/x wave bands corrected after the first successive step are vertical
The difference of polarization reflectivity.By the data processing of above-mentioned two step, the difference of C-band and X-band vertical polarization reflectivity is corrected,
And above-mentioned data correction mode is to maintain data point using data as two dimensional image processing and (observe the vertical pole of c/x wave bands
Change reflectivity difference and simulation c/x wave band vertical polarization reflectivity the corresponding data point of difference) between topological property it is constant.
S104, the difference using default look-up table to the observation c/x wave band vertical polarization reflectivity of optimization and corresponding auxiliary
Data carry out Inversion Calculation, obtain the sea surface salinity of target water area.
In the embodiment of the present invention, previously according to dielectric constant model and Fresnel reflection law, calculates sea surface salinity 0 and arrive
The difference of c/x wave band vertical polarization reflectivity in 40psu and sea surface temperature 290K to 315K ranges, and using the value as sea salt
The function of degree and sea surface temperature generates and presets look-up table, then utilizes the default look-up table to the sight after the optimization of above-mentioned calculating
The difference and corresponding auxiliary data for surveying c/x wave band vertical polarization reflectivity carry out inversion procedure, you can inverting obtains target sea
The sea surface salinity in region (is searched in default look-up table and the difference of above-mentioned observation c/x wave bands vertical polarization reflectivity and corresponding
The corresponding sea surface salinity of auxiliary data).
Specifically, the sea surface salinity of inverting is by identical sea surface temperature, model calculate and the vertical polarization reflectivity of observation it
The difference of difference, which minimizes, to be determined.
A kind of method based on satellite-borne microwave radiometer inverting sea surface salinity provided in an embodiment of the present invention, by spaceborne micro-
Wave radiation meter obtains the observation bright temperature data of target water area wave band c and x wave bands and corresponding auxiliary data respectively, and calculates
Observe the difference of c/x wave band vertical polarization reflectivity；Then, digitized map is carried out to the difference for observing c/x wave band vertical polarization reflectivity
As optimization processing；Finally, using default look-up table to the difference and supplementary number of the reflectivity of observation wave band c and x wave bands after optimization
According to carry out Inversion Calculation, obtain the sea surface salinity of target water area, with it is in the prior art directly by the data of acquisition it is anti-
The method for drilling sea surface salinity so that the salt angle value that inverting obtains is compared with the deviation of practical salt angle value is larger, to observing bright temperature number
According to result of calculation optimize processing so that result of calculation is closer to actual value, but also being obtained by the calculating structure inversion
The distribution situation of the sea surface salinity arrived is closer to actual value, to improve the inversion accuracy of sea surface salinity.
As described above, the observation bright temperature data of wave band c and x wave bands is obtained in the embodiment of the present invention using HY-2a radiometers,
In view of the observation data of above-mentioned HY-2a instruments may have systematic deviation, need to be modified the data of acquisition,
With reference to figure 2, the correcting mode of the embodiment of the present invention specifically comprises the following steps：
S201, according to dielectric constant model and auxiliary data, the difference of calculating simulation c/x wave band vertical polarization reflectivity.
The dielectric constant model that the embodiment of the present invention uses is indicated that Parameter Expression is published in AMSR-E by Debye equation
ATBD (2000).
S202, to observe c/x wave band vertical polarization reflectivity difference with simulation c/x wave band vertical polarization reflectivity difference into
Row comparison statistical disposition obtains including difference and the simulation c/x wave band vertical polarizations reflection for observing c/x wave band vertical polarization reflectivity
The two-dimensional digital image of the order of magnitude of the difference of rate.
Specifically, according to the difference of observation c/x wave band vertical polarization reflectivity and simulation c/x wave band vertical polarizations reflectivity it
The two-dimentional dimension word image (as shown in Figure 3) of difference structure, the x-axis in the two-dimensional digital image where pictorial element indicate that model reflects
Rate, y-axis indicate observation reflectivity, and the value of pixel is statistics number；And the difference and mould of c/x wave band vertical polarization reflectivity will be observed
The comparison that the difference of quasi- c/x wave band vertical polarization reflectivity carries out, the two-dimensional digital image include that observation c/x wave band vertical polarizations are anti-
Penetrate the corresponding order of magnitude of difference of the difference and simulation c/x wave band vertical polarization reflectivity of rate.
The order of magnitude with simulation c/x of the difference of observation c/x wave band vertical polarization reflectivity in S203, adjustment two-dimensional digital image
The order-of-magnitude agreement of the difference of wave band vertical polarization reflectivity obtains revised two-dimensional digital image for the first time.
From Fig. 3, it can be seen that observe the reflection of c/x wave band vertical polarizations in the two-dimensional digital image built in step 202
The difference of rate and the corresponding order of magnitude of difference of simulation c/x wave band vertical polarization reflectivity are not corresponding, therefore needing will be above-mentioned
The order of magnitude adjustment of the two is consistent, and specific method is, by the order of magnitude of the difference of the observation c/x wave band vertical polarization reflectivity in figure
It is adjusted to and simulates the order-of-magnitude agreement of the difference of c/x wave band vertical polarization reflectivity, you can obtain revised two dimension for the first time
Digital picture.
S204, second of correcting process is carried out to revised two-dimensional digital image for the first time using data rotation method,
The difference of the observation c/x wave band vertical polarization reflectivity optimized.
Specifically, in order to enable the closer actual value of c/x wave band vertical polarization reflectivity is observed, in the embodiment of the present invention
Above-mentioned two-dimensional digital image after being corrected first time using data rotation method is still carried out into second of correcting process.
The difference and simulation c/x waves of c/x wave band vertical polarization reflectivity are observed in the two-dimensional digital image built in step 202
The corresponding order of magnitude of difference of section vertical polarization reflectivity is not corresponding, therefore needs the order of magnitude of said two devices adjusting one
It causes, the specific method is as follows：
The difference for simulating c/x wave band vertical polarization reflectivity is expressed as by observation c/ using the method for least squares line fitting
The function of the difference of x wave band vertical polarization reflectivity, obtains update equation Rmodel=0.0017Robs+0.0107；Wherein, RmodelIt is
The difference (difference for simulating c/x wave band vertical polarization reflectivity) of the reflectivity of model；RobsIt is that the difference of reflectivity of observation (is seen
Survey the difference of c/x wave band vertical polarization reflectivity)；
According to update equation Rmodel=0.0017Robs+ 0.0107, it is vertical to adjust observation c/x wave bands in two-dimensional digital image
The order-of-magnitude agreement of the order of magnitude of the difference of polarization reflectivity and the difference of simulation c/x wave band vertical polarization reflectivity.Wherein, it is corrected
The two-dimensional digital image obtained afterwards is as shown in Figure 4.
In the embodiment of the present invention, in order to enable the closer actual value of observation c/x wave band vertical polarization reflectivity, the present invention
The above-mentioned two-dimensional digital image after being corrected first time using data rotation method is still carried out in embodiment and is repaiied for the second time
Positive processing, with reference to figure 5, the specific implementation of above-mentioned steps 204 is as follows：
S301, by the difference of the observation c/x wave band vertical polarization reflectivity in revised two-dimensional digital image for the first time and
The difference of simulation c/x wave band vertical polarization reflectivity carries out pairing statistical disposition, obtains the two-dimensional digital image with data value；Its
In, the x-axis where the pictorial element of two-dimensional digital image indicates that model reflectivity, y-axis indicate observation reflectivity.
The center of mass point of S302, the two dimensional image constituted according to all pictorial elements in two-dimensional digital image determine two dimension
The rotation axis of digital picture.
S303, according to the APPROXIMATE DISTRIBUTION straight line of all pictorial elements in two-dimensional digital image, determine two-dimensional digital image
Rotation angle.
S304, according to rotation axis and rotation angle, rotation processing is carried out to revised two-dimensional digital image for the first time,
Obtain second of revised two-dimensional digital image.
S305, the difference of c/x wave band vertical polarization reflectivity is observed in extraction from second of revised two-dimensional digital image,
The difference of observation c/x wave band vertical polarization reflectivity as an optimization.
With reference to figure 6, Fig. 7 and Fig. 8, in conjunction with above-mentioned steps 301-305, specifically, will first time modified two-dimensional digital figure
The difference (i.e. multiple data) of simulation c/x wave band vertical polarization reflectivity as in and observe c/x wave band vertical polarizations reflectivity it
Poor (i.e. multiple data) carry out pairing statistics, and the numerical digit of above-mentioned data is rounded to 0.000001.In this way, being counted in pairing
It just constitutes an abscissa later and ordinate all has data value two-dimensional digital image (as shown in Figure 6), wherein above-mentioned two dimension
The x-axis of image picture elements (i.e. pictorial element) is model reflectivity in digital picture, and y-axis is observation reflectivity, and the value of pixel is system
Count number.Then, the threshold value of above-mentioned pixel is taken as 0,1,5,10 respectively, the data that number is less than to threshold value remove, reservation
The individual percentage of former data is respectively 100%, 93.0%, 78.8%, 69.5% shared by data.And above-mentioned two-dimensional digital image master
The distribution of volume data be also be approximately straight line.Therefore the method for data rotation can be used to revised data for the first time
It is corrected again；Wherein, above-mentioned rotation key is the positioning of the selection and rotation axis (perpendicular to plane) of rotation angle.
Consider the barycenter (Centre of Mass) of the two dimensional image that all pixels are constituted as rotation axis by point,
The coordinate of the point is (0.0107577,0.010779), is located approximately on 45 ° of straight lines of origin；The angle of above-mentioned rotation is by main body
The near linear of data distribution determines, in order to obtain relevant parameter (the i.e. slope of linear equation and section according to ginseng of this straight line
Number), by all pixels of image according to model reflectivity be divided into from small to large 12 parts (preceding 11 part respectively include 1000
Pixel, last part include 784 pixels), find the center-of-mass coordinate of each part.The distribution of this 12 pairs of coordinates is approximately one
Straight line, equation are：Y=0.011057x+0.010660；
This straight line and x-axis angle are 0.63349 degree, therefore rotation angle is to rotate θ=44.46651 ° counterclockwise.Rotation
After turning, barycenter was moved on 45 ° of straight lines of origin (x coordinate is constant), second step is obtained and corrects data.
Then observation c/x wave band vertical polarization reflectivity can be obtained from above-mentioned second of modified two-dimensional digital image
Difference, the difference of observation c/x wave band vertical polarization reflectivity as an optimization, and participate in subsequent calculating.
In addition, in the embodiment of the present invention, the specific of the difference of c/x wave band vertical polarization reflectivity is observed in above-mentioned steps 102
Calculation is as follows：
The difference for calculating the vertical polarization reflectivity of wave band c and the vertical polarization reflectivity of x wave bands obtains observation c/x waves
The difference of section vertical polarization reflectivity.
Specifically, radiative transfer model (Radiative Transfer Model, RTM) is for describing microwave radiometer spy
The relationship of the physical factor arrived involved in the bright temperature surveyed and road radiation transmission process, can be expressed as following form：TB=TBU+τ
[ETs+TBΩ](1)
Wherein, TBIndicate the bright temperature of Based Microwave Radiometer Detection, TBUIt is air uplink radiation (upwelling
Atmospheric brightness temperature), TBΩIt is surface scattering (the sky radiation of atmospheric radiation
Scattered upward by Earth surface), TsIt is sea surface temperature, τ is atmospheric transmittance (atmospheric
Transmission), E is sea surface emissivity (sea-surface emissivity).TBΩIn the RTM of different parametrizations
Specific expression-form is different, can be unified for following equation：
TBΩ=RM (2)
Wherein, R is the corresponding reflectivity of sea surface emissivity, R=1-E.M is one and contains Downward atmospheric long-wave radiation, air
Transmitance, surface scattering correction, cosmic background radiation function.According to equation (1) and equation (2), sea surface reflection rate can be with table
Up to for following equation：
It can be obtained from the relational graph of vertical polarization reflectivity, sea surface temperature and sea surface salinity as drawn a conclusion：(1) X and C waves
The vertical polarization reflectivity of section is far below the response to sea surface temperature about the response of sea surface salinity.(2) vertical polarization of C-band
Reflectivity changes sea surface temperature and the variation of sea surface salinity is responded.(3) the vertical polarization reflectivity of X-band is to sea
Temperature change is responded, and is not responded to substantially to sea surface salinity variation in low temperature range.(4) vertical polarization of X and C-band
The difference of reflectivity changes sea surface temperature and sea surface salinity variation has apparent response.
Observe bright temperature using HY-2a scanning microwave radiometers in the embodiment of the present invention, and using the L2b products of AMSR-2 and
WOA09 Salinity Datas are as auxiliary data, the inverting sea surface salinity in CHINESE OFFSHORE sea area, by the experimental results showed that, pass through
The quality of HY-2a salinity retrieval products after data correction is better than the inverting of the L1a bright temperature datas directly using AMSR-E
As a result, its RMSE in entire CHINESE OFFSHORE sea area is 0.4086psu.
A kind of method based on satellite-borne microwave radiometer inverting sea surface salinity provided in an embodiment of the present invention, with the prior art
In the direct data inversion sea surface salinity by acquisition method so that salt angle value that inverting obtains and the deviation of practical salt angle value
Larger to compare, the result of calculation to observing bright temperature data optimizes processing so that result of calculation also makes closer to actual value
The distribution situation of the sea surface salinity obtained by the calculating structure inversion is obtained closer to actual value, to improve sea surface salinity
Inversion accuracy.
The present invention also provides a kind of device based on satellite-borne microwave radiometer inverting sea surface salinity, described device is for holding
Method of the row based on satellite-borne microwave radiometer inverting sea surface salinity, as shown in figure 9, described device includes：
Acquisition module 11, the sight for obtaining target water area wave band c and x wave bands respectively by satellite-borne microwave radiometer
Survey bright temperature data auxiliary data corresponding with observation bright temperature data；
Computing module 12 is used for the observation bright temperature data according to c/x wave bands and corresponding auxiliary data, calculating observation c/x
The difference of wave band vertical polarization reflectivity；
Digital picture optimization processing module 13, for carrying out digitized map to the difference for observing c/x wave band vertical polarization reflectivity
As optimization processing, the difference of the observation c/x wave band vertical polarization reflectivity optimized；
Inversion Calculation module 14, for using preset look-up table to the observation c/x wave band vertical polarizations reflectivity of optimization it
Difference and auxiliary data carry out Inversion Calculation, obtain the sea surface salinity of target water area.
As described above, the observation bright temperature data of wave band c and x wave bands is obtained in the embodiment of the present invention using HY-2a radiometers,
In view of the observation data of above-mentioned HY-2a instruments may have systematic deviation, need to be modified the data of acquisition,
With reference to figure 10, digital picture optimization processing module 13 includes：
First computing unit 131, for according to dielectric constant model and auxiliary data, the vertical pole of calculating simulation c/x wave bands
Change the difference of reflectivity；
Statistical disposition unit 132 is compared, for the difference and simulation c/x wave bands to observing c/x wave band vertical polarization reflectivity
The difference of vertical polarization reflectivity carries out comparison statistical disposition, obtains including the difference and mould for observing c/x wave band vertical polarization reflectivity
The two-dimensional digital image of the order of magnitude of the difference of quasi- c/x wave band vertical polarization reflectivity；
First time correcting process unit 133 observes c/x wave band vertical polarization reflectivity for adjusting in two-dimensional digital image
Difference the order of magnitude and simulate c/x wave band vertical polarization reflectivity difference order-of-magnitude agreement, obtain for the first time revised two
Dimension word image；
Second of correcting process unit 134, for utilizing data rotation method to revised two-dimensional digital figure for the first time
As carrying out second of correcting process, the difference of the observation c/x wave band vertical polarization reflectivity optimized.
Further, as shown in Figure 10, first time correcting process unit 133 includes：
Update equation obtains subelement 1331, is hung down for the method using least squares line fitting by c/x wave bands are simulated
The difference of straight polarization reflectivity is expressed as the function of the difference of observation c/x wave band vertical polarization reflectivity, obtains update equation Rmodel=
0.0017Robs+0.0107；
Subelement 1332 is adjusted, for according to update equation Rmodel=0.0017Robs+ 0.0107, adjust two-dimensional digital figure
The number of the difference of the order of magnitude of the difference of observation c/x wave band vertical polarization reflectivity and simulation c/x wave band vertical polarization reflectivity as in
Magnitude is consistent.
Further, as shown in figure 11, second of correcting process unit 134 includes：
Statistical disposition subelement 1341 is matched, is used for the observation c/x waves in revised two-dimensional digital image for the first time
The difference of section vertical polarization reflectivity and the difference of simulation c/x wave band vertical polarization reflectivity carry out pairing statistical disposition, are had
The two-dimensional digital image of data value；Wherein, the x-axis where the pictorial element of two-dimensional digital image indicates model reflectivity, y-axis table
Show observation reflectivity；
Rotation axis determination subelement 1342, the two dimension for being constituted according to all pictorial elements in two-dimensional digital image
The center of mass point of image determines the rotation axis of two-dimensional digital image；
Rotation angle determination subelement 1343, according to the APPROXIMATE DISTRIBUTION straight line of all pictorial elements in two-dimensional digital image,
Determine the rotation angle of two-dimensional digital image；
Rotation processing subelement 1344 is used for according to rotation axis and rotation angle, to revised two-dimemsional number for the first time
Word image carries out rotation processing, obtains second of revised two-dimensional digital image；
Subelement 1345 is extracted, it is vertical for extracting observation c/x wave bands from second of revised two-dimensional digital image
The difference of polarization reflectivity, the difference of observation c/x wave band vertical polarization reflectivity as an optimization.
In addition, in the embodiment of the present invention, as shown in figure 11, computing module 12 includes：
Third computing unit 122, for calculating the vertical polarization reflectivity of wave band c and the vertical polarization reflectivity of x wave bands
Difference, obtain observation c/x wave band vertical polarization reflectivity difference.
A kind of device based on satellite-borne microwave radiometer inverting sea surface salinity provided in an embodiment of the present invention, with the prior art
In the direct data inversion sea surface salinity by acquisition method so that salt angle value that inverting obtains and the deviation of practical salt angle value
Larger to compare, the result of calculation to observing bright temperature data optimizes processing so that result of calculation also makes closer to actual value
The distribution situation of the sea surface salinity obtained by the calculating structure inversion is obtained closer to actual value, to improve sea surface salinity
Inversion accuracy.
In several embodiments provided herein, it should be understood that disclosed systems, devices and methods, it can be with
It realizes by another way.The apparatus embodiments described above are merely exemplary, for example, the division of the unit,
Only a kind of division of logic function, formula that in actual implementation, there may be another division manner, in another example, multiple units or component can
To combine or be desirably integrated into another system, or some features can be ignored or not executed.Another point, it is shown or beg for
The mutual coupling, direct-coupling or communication connection of opinion can be by some communication interfaces, device or unit it is indirect
Coupling or communication connection can be electrical, machinery or other forms.
The unit illustrated as separating component may or may not be physically separated, aobvious as unit
The component shown may or may not be physical unit, you can be located at a place, or may be distributed over multiple
In network element.Some or all of unit therein can be selected according to the actual needs to realize the mesh of this embodiment scheme
's.
In addition, each functional unit in each embodiment of the present invention can be integrated in a processing unit, it can also
It is that each unit physically exists alone, it can also be during two or more units be integrated in one unit.
It, can be with if the function is realized in the form of SFU software functional unit and when sold or used as an independent product
It is stored in a computer read/write memory medium.Based on this understanding, technical scheme of the present invention is substantially in other words
The part of the part that contributes to existing technology or the technical solution can be expressed in the form of software products, the meter
Calculation machine software product is stored in a storage medium, including some instructions are used so that a computer equipment (can be
People's computer, server or network equipment etc.) it performs all or part of the steps of the method described in the various embodiments of the present invention.
And storage medium above-mentioned includes：USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited
The various media that can store program code such as reservoir (RAM, Random Access Memory), magnetic disc or CD.
The above description is merely a specific embodiment, but scope of protection of the present invention is not limited thereto, any
Those familiar with the art in the technical scope disclosed by the present invention, can easily think of the change or the replacement, and should all contain
Lid is within protection scope of the present invention.Therefore, protection scope of the present invention should be based on the protection scope of the described claims.
Claims (10)
1. a kind of method based on satellite-borne microwave radiometer inverting sea surface salinity, which is characterized in that the method includes：
Obtain the observation bright temperature data of target water area wave band c and x wave bands and the sight respectively by satellite-borne microwave radiometer
Survey the corresponding auxiliary data of bright temperature data；
According to the observation bright temperature data of c/x wave bands and corresponding auxiliary data, calculating observation c/x wave band vertical polarizations reflectivity it
Difference；
Digital picture optimization processing, the observation c/x optimized are carried out to the difference of the observation c/x wave band vertical polarization reflectivity
The difference of wave band vertical polarization reflectivity；
The difference of the observation c/x wave band vertical polarization reflectivity of the optimization and the auxiliary data are carried out using default look-up table
Inversion Calculation obtains the sea surface salinity of the target water area.
2. according to the method described in claim 1, it is characterized in that, described to the observation c/x wave band vertical polarization reflectivity
Difference carry out digital picture optimization processing, the difference of the observation c/x wave band vertical polarization reflectivity optimized includes：
According to dielectric constant model and the auxiliary data, the difference of calculating simulation c/x wave band vertical polarization reflectivity；
Difference to the observation c/x wave band vertical polarization reflectivity and the difference for simulating c/x wave band vertical polarization reflectivity carry out pair
Than statistical disposition, obtain including the difference of the observation c/x wave band vertical polarization reflectivity and the simulation c/x wave band vertical polarizations
The two-dimensional digital image of the order of magnitude of the difference of reflectivity；
Adjust the order of magnitude and simulation c/x wave bands of the difference that c/x wave band vertical polarization reflectivity is observed in the two-dimensional digital image
The order-of-magnitude agreement of the difference of vertical polarization reflectivity obtains revised two-dimensional digital image for the first time；
Second of correcting process is carried out to revised two-dimensional digital image for the first time using data rotation method, is optimized
Observe the difference of c/x wave band vertical polarization reflectivity.
3. according to the method described in claim 2, it is characterized in that, observing c/x waves in the adjustment two-dimensional digital image
The order of magnitude of the difference of section vertical polarization reflectivity and the order-of-magnitude agreement of the difference of simulation c/x wave band vertical polarization reflectivity include：
The difference of the simulation c/x wave band vertical polarization reflectivity is expressed as by observation c/ using the method for least squares line fitting
The function of the difference of x wave band vertical polarization reflectivity, obtains update equation Rmodel=0.0017Robs+0.0107；Wherein, RmodelIt is
Simulate the difference of c/x wave band vertical polarization reflectivity；RobsIt is the difference for observing c/x wave band vertical polarization reflectivity；
According to the update equation Rmodel=0.0017Robs+ 0.0107, adjust observation c/x wave bands in the two-dimensional digital image
The order-of-magnitude agreement of the order of magnitude of the difference of vertical polarization reflectivity and the difference of simulation c/x wave band vertical polarization reflectivity.
4. according to the method described in claim 3, it is characterized in that, described revised to first time using data rotation method
Two-dimensional digital image carries out second of correcting process, and the difference of the observation c/x wave band vertical polarization reflectivity optimized includes：
By the difference of the observation c/x wave band vertical polarization reflectivity in revised two-dimensional digital image for the first time and simulation c/x waves
The difference of section vertical polarization reflectivity carries out pairing statistical disposition, obtains the two-dimensional digital image with data value；Wherein, described two
X-axis where the pictorial element of dimension word image indicates that model reflectivity, y-axis indicate observation reflectivity；
According to the center of mass point for the two dimensional image that all pictorial elements in the two-dimensional digital image are constituted, the two-dimemsional number is determined
The rotation axis of word image；
According to the APPROXIMATE DISTRIBUTION straight line of all pictorial elements in the two-dimensional digital image, the rotation of the two-dimensional digital image is determined
Gyration；
According to the rotation axis and the rotation angle, rotation processing is carried out to revised two-dimensional digital image for the first time,
Obtain second of revised two-dimensional digital image；
The difference of extraction observation c/x wave band vertical polarization reflectivity from second of revised two-dimensional digital image, as
The difference of the observation c/x wave band vertical polarization reflectivity of optimization.
5. according to the method described in claim 4, it is characterized in that, the observation bright temperature data and correspondence according to c/x wave bands
Auxiliary data, the difference of calculating observation c/x wave band vertical polarization reflectivity includes：
The vertical polarization reflectivity of wave band c and x wave bands is calculated separately by following formulaWherein, TB
Indicate the bright temperature of Based Microwave Radiometer Detection, TBUIt is air uplink radiation, TBΩIt is the surface scattering of atmospheric radiation, TsIt is sea temperature
Degree, τ is atmospheric transmittance, and E is sea surface emissivity, and R is the corresponding reflectivity of sea surface emissivity, and R=1-E, M are one and contain
Downward atmospheric long-wave radiation, atmospheric transmittance, surface scattering correct and the function of cosmic background radiation；
The difference for calculating the vertical polarization reflectivity of wave band c and the vertical polarization reflectivity of x wave bands obtains observation c/x wave bands and hangs down
The difference of straight polarization reflectivity.
6. a kind of device based on satellite-borne microwave radiometer inverting sea surface salinity, which is characterized in that described device includes：
Acquisition module, the bright temperature of observation for obtaining target water area wave band c and x wave bands respectively by satellite-borne microwave radiometer
Data and the corresponding auxiliary data of the observation bright temperature data；
Computing module, is used for the observation bright temperature data according to c/x wave bands and corresponding auxiliary data, and calculating observation c/x wave bands hang down
The difference of straight polarization reflectivity；
Digital picture optimization processing module carries out digital picture for the difference to the observation c/x wave band vertical polarization reflectivity
Optimization processing, the difference of the observation c/x wave band vertical polarization reflectivity optimized；
Inversion Calculation module, for the difference using default look-up table to the observation c/x wave band vertical polarization reflectivity of the optimization
Inversion Calculation is carried out with the auxiliary data, obtains the sea surface salinity of the target water area.
7. device according to claim 6, which is characterized in that the digital picture optimization processing module includes：
First computing unit, for according to dielectric constant model and the auxiliary data, calculating simulation c/x wave band vertical polarizations to be anti-
Penetrate the difference of rate；
Statistical disposition unit is compared, it is vertical with simulation c/x wave bands for the difference to the observation c/x wave band vertical polarization reflectivity
Polarization reflectivity difference carry out comparison statistical disposition, obtain include the observation c/x wave band vertical polarization reflectivity difference and institute
State the two-dimensional digital image of the order of magnitude of the difference of simulation c/x wave band vertical polarization reflectivity；
First time correcting process unit, for adjust observed in the two-dimensional digital image c/x wave band vertical polarizations reflectivity it
The order-of-magnitude agreement of the order of magnitude of difference and the difference of simulation c/x wave band vertical polarization reflectivity, obtains revised two dimension for the first time
Digital picture；
Second of correcting process unit, for using data rotation method to revised two-dimensional digital image progress for the first time the
Second-order correction processing, the difference of the observation c/x wave band vertical polarization reflectivity optimized.
8. device according to claim 7, which is characterized in that the first time correcting process unit includes：
Update equation obtains subelement, for the method using least squares line fitting by the vertical pole of the simulation c/x wave bands
The difference for changing reflectivity is expressed as observing the function of the difference of c/x wave band vertical polarization reflectivity, obtains update equation Rmodel=
0.0017Robs+0.0107；Wherein, RmodelIt is the difference for simulating c/x wave band vertical polarization reflectivity；RobsIt is that observation c/x wave bands hang down
The difference of straight polarization reflectivity；
Subelement is adjusted, for according to the update equation Rmodel=0.0017Robs+ 0.0107, adjust the two-dimensional digital figure
The number of the difference of the order of magnitude of the difference of observation c/x wave band vertical polarization reflectivity and simulation c/x wave band vertical polarization reflectivity as in
Magnitude is consistent.
9. device according to claim 8, which is characterized in that second of correcting process unit includes：
Statistical disposition subelement is matched, is used for the vertical pole of observation c/x wave bands in revised two-dimensional digital image for the first time
The difference of the difference and simulation c/x wave band vertical polarization reflectivity of changing reflectivity carries out pairing statistical disposition, obtains having data value
Two-dimensional digital image；Wherein, the x-axis where the pictorial element of the two-dimensional digital image indicates that model reflectivity, y-axis indicate to see
Survey reflectivity；
Rotation axis determination subelement, the two dimensional image for being constituted according to all pictorial elements in the two-dimensional digital image
Center of mass point, determine the rotation axis of the two-dimensional digital image；
Rotation angle determination subelement is determined according to the APPROXIMATE DISTRIBUTION straight line of all pictorial elements in the two-dimensional digital image
The rotation angle of the two-dimensional digital image；
Rotation processing subelement is used for according to the rotation axis and the rotation angle, to revised two-dimemsional number for the first time
Word image carries out rotation processing, obtains second of revised two-dimensional digital image；
Subelement is extracted, for the extraction observation c/x wave band vertical polarizations from second of revised two-dimensional digital image
The difference of reflectivity, the difference of observation c/x wave band vertical polarization reflectivity as an optimization.
10. device according to claim 9, which is characterized in that the computing module includes：
Second computing unit, the vertical polarization reflectivity for calculating separately wave band c and x wave bands by following formulaWherein, TBIndicate the bright temperature of Based Microwave Radiometer Detection, TBUIt is air uplink radiation, TBΩIt is air
The surface scattering of radiation, TsIt is sea surface temperature, τ is atmospheric transmittance, and E is sea surface emissivity, and R is that sea surface emissivity is corresponding anti-
Rate, R=1-E are penetrated, M is one and contains Downward atmospheric long-wave radiation, atmospheric transmittance, surface scattering correction and cosmic background radiation
Function；
Third computing unit, the difference for calculating the vertical polarization reflectivity of wave band c and the vertical polarization reflectivity of x wave bands,
Obtain the difference of observation c/x wave band vertical polarization reflectivity.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610087880.XA CN105737806B (en) | 2016-02-16 | 2016-02-16 | Method and apparatus based on satellite-borne microwave radiometer inverting sea surface salinity |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610087880.XA CN105737806B (en) | 2016-02-16 | 2016-02-16 | Method and apparatus based on satellite-borne microwave radiometer inverting sea surface salinity |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105737806A CN105737806A (en) | 2016-07-06 |
| CN105737806B true CN105737806B (en) | 2018-08-17 |
Family
ID=56245221
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610087880.XA Expired - Fee Related CN105737806B (en) | 2016-02-16 | 2016-02-16 | Method and apparatus based on satellite-borne microwave radiometer inverting sea surface salinity |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105737806B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106092336B (en) * | 2016-07-20 | 2018-11-13 | 中国科学院国家空间科学中心 | A kind of one-dimensional interference type micro-wave Radiometer Image inversion method |
| CN109145494B (en) * | 2018-09-11 | 2020-04-07 | 北京师范大学 | Sea surface temperature inversion method and system |
| CN110411928B (en) * | 2019-08-02 | 2020-09-01 | 山东农业大学 | Method and device for obtaining concentration of suspended seawater sand based on airborne double-color laser point cloud |
| CN110793649B (en) * | 2019-10-14 | 2020-10-23 | 三亚中科遥感研究所 | Method for correcting brightness and temperature of rough sea surface by using backscattering cross section |
| CN112414554B (en) * | 2020-12-02 | 2021-12-07 | 国家卫星海洋应用中心 | Sea surface salinity obtaining method, device, equipment and medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104198052A (en) * | 2014-09-25 | 2014-12-10 | 国家卫星海洋应用中心 | Acquisition method for sea ice concentration on basis of ocean No. II satellite scanning microwave radiometer |
| CN105067120A (en) * | 2015-07-28 | 2015-11-18 | 国家卫星海洋应用中心 | Dynamic filtering and re-sampling method and dynamic filtering and re-sampling device for brightness temperature observation of space-borne microwave radiometer |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101163362B1 (en) * | 2011-11-16 | 2012-07-09 | 한국해양연구원 | Oceanographic waves observation system and method |
| KR101576167B1 (en) * | 2014-05-09 | 2015-12-09 | 한국해양과학기술원 | Apparatus for monitoring sea surface current |
-
2016
- 2016-02-16 CN CN201610087880.XA patent/CN105737806B/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104198052A (en) * | 2014-09-25 | 2014-12-10 | 国家卫星海洋应用中心 | Acquisition method for sea ice concentration on basis of ocean No. II satellite scanning microwave radiometer |
| CN105067120A (en) * | 2015-07-28 | 2015-11-18 | 国家卫星海洋应用中心 | Dynamic filtering and re-sampling method and dynamic filtering and re-sampling device for brightness temperature observation of space-borne microwave radiometer |
Non-Patent Citations (4)
| Title |
|---|
| "The WISE 2000 and 2001 field experiments in support of the SMOS mission: Sea surface L-band brightness temperature observations and their application to sea surface salinity retrieval";Camps A;《 IEEE Transactions on Geoscience and Remote Sensing》;20040430;第42卷(第2期);第804-823页 * |
| "一种用于微波辐射计遥感海表面盐度和温度的反演算法";殷晓斌;《中国科学 D辑 地球科学》;20061020;第36卷(第10期);第968-976页 * |
| "基于星载微波辐射低频双极化通道观测数据的微波辐射传输模型";卢勇夺;《中国海洋大学学报》;20130215;第43卷(第2期);第77-85页 * |
| "微波辐射计遥感海洋盐度的研究进展";李青侠;《海洋技术》;20070915;第26卷(第3期);第58-63页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105737806A (en) | 2016-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105737806B (en) | Method and apparatus based on satellite-borne microwave radiometer inverting sea surface salinity | |
| Sanderman et al. | A global map of mangrove forest soil carbon at 30 m spatial resolution | |
| Zhong et al. | Climate change trends and impacts on vegetation greening over the Tibetan Plateau | |
| Lay et al. | Ionospheric acoustic and gravity waves associated with midlatitude thunderstorms | |
| Elbelrhiti et al. | Barchan dune corridors: field characterization and investigation of control parameters | |
| Medvedev et al. | Relation of internal gravity wave anisotropy with neutral wind characteristics in the upper atmosphere | |
| Lang et al. | Transient luminous events above two mesoscale convective systems: Storm structure and evolution | |
| Le Pichon et al. | Infrasound monitoring of volcanoes to probe high‐altitude winds | |
| Ridley et al. | Toward resolution‐independent dust emissions in global models: Impacts on the seasonal and spatial distribution of dust | |
| CN104268423B (en) | Large scale dynamic evolution Sand-dust type aerosol inversion method | |
| Hare et al. | LOFAR lightning imaging: Mapping lightning with nanosecond precision | |
| Ajith et al. | On the fresh development of equatorial plasma bubbles around the midnight hours of June solstice | |
| Anguelova et al. | Whitecap fraction from satellite measurements: Algorithm description | |
| KR20100074264A (en) | Geospatial modeling system and related method using multiple sources of geographic information | |
| Ajith et al. | Explicit characteristics of evolutionary‐type plasma bubbles observed from Equatorial Atmosphere Radar during the low to moderate solar activity years 2010–2012 | |
| Delgado‐Fernandez et al. | Event‐scale dynamics of a parabolic dune and its relevance for mesoscale evolution | |
| CN108961199A (en) | Multi- source Remote Sensing Data data space-time fusion method and device | |
| Bór et al. | On the series of+ CG lightning strokes in dancing sprite events | |
| Goumehei et al. | Surface water body detection in polarimetric SAR data using contextual complex Wishart classification | |
| Curtarelli et al. | Evaluation of the Tropical Rainfall Measuring Mission 3B43 product over an inland area in Brazil and the effects of satellite boost on rainfall estimates | |
| Pardo‐igúzquiza et al. | Morphometric analysis of karst depressions on a Mediterranean karst massif | |
| Jonah et al. | Observation of TEC perturbation associated with medium‐scale traveling ionospheric disturbance and possible seeding mechanism of atmospheric gravity wave at a Brazilian sector | |
| Chybicki | Three-dimensional geographically weighted inverse regression (3GWR) model for satellite derived bathymetry using Sentinel-2 observations | |
| Favalli et al. | Digital elevation model construction from structured topographic data: the DEST algorithm | |
| CN105842259B (en) | The cold empty external calibration method and apparatus of spaceborne salinometer |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180817 Termination date: 20190216 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |