CN104198052A - Acquisition method for sea ice concentration on basis of ocean No. II satellite scanning microwave radiometer - Google Patents
Acquisition method for sea ice concentration on basis of ocean No. II satellite scanning microwave radiometer Download PDFInfo
- Publication number
- CN104198052A CN104198052A CN201410498638.2A CN201410498638A CN104198052A CN 104198052 A CN104198052 A CN 104198052A CN 201410498638 A CN201410498638 A CN 201410498638A CN 104198052 A CN104198052 A CN 104198052A
- Authority
- CN
- China
- Prior art keywords
- bright
- temperature data
- ice
- frequency range
- wyntet
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Abstract
The invention relates to the technical field of satellite communication and particularly relates to an acquisition method for sea ice concentration on the basis of the Ocean No. II satellite scanning microwave radiometer. The acquisition method for the sea ice concentration on the basis of the Ocean II satellite scanning microwave radiometer comprises the following steps: according to a first brightness temperature data of a typically characteristic area, a brightness temperature characteristic value is obtained; the typical characteristic area comprises a characteristic area of one-year ice, a characteristic area of multiyear ice and characteristic area of sea water; the first brightness temperature data is acquired through the Ocean II satellite scanning microwave radiometer; according to the brightness temperature characteristic value and a second brightness temperature data of a monitored area, dynamic brightness temperature characteristic value of the monitored area is acquired; the second brightness temperature data is acquired through the Ocean No. II satellite scanning microwave radiometer; according to the dynamic brightness temperature characteristic value of the monitored area, the dynamic sea ice concentration of the monitored area is acquired. The acquisition method for sea ice concentration on the basis of the Ocean II satellite scanning microwave radiometer overcomes the technical problem that the sea ice concentration of a correspondingly monitored area cannot be acquired on the basis of the data acquired by the Ocean II satellite in the prior art.
Description
Technical field
The present invention relates to technical field of satellite communication, in particular to the ice concentration acquisition methods based on No. two satellite scanning microwave radiometers in ocean.
Background technology
Ice concentration is to describe the major parameter of polar region sea ice, is defined as sea ice in unit area and covers shared number percent.At present be mainly derived from satellite-borne microwave radiometer for the data of ice concentration inverting, this radiometer, as a kind of passive microwave sensor, in conjunction with the bright temperature of observation under the different polarization modes of different-waveband, can be distinguished sea ice and seawater.Ice and snow data center of the U.S. provides the two poles of the earth ice concentration data that the spatial resolution of having utilized satellite-borne microwave radiometer to obtain since 1978 is 25km, and the data of multiple microwave radiometers are for the product of this long-time sequence of inverting, as the SMMR on U.S. Nimbus-7 and Seasat-A (Scanning Multi-frequency Microwave Radiometer), SSMI (Special Sensor Microwave/Imager) on U.S. DMSP (Defense Meteorological Satellite Program), AMSR-E (Advanced Microwave Scanning Radiometer for EOS) on U.S. Aqua satellite etc.Now DMSP-F17 SSMIS, DMSP-F18 SSMIS in orbit, GCOM-W1 AMSR2 provide data source continuing as this product.
No. two, ocean HY-2 satellite is the satellite for ocean dynamical environment monitoring of China's independent research.In HY-2 satellite, be equipped with radar altimeter, microwave scatterometer, scanning microwave radiometer and proofread and correct the load such as radiometer, wherein scanning microwave radiometer frequency of operation is 6.6GHz, 10.7GHz, 18.7GHz, 23.8GHz and 27.0GHz, except 23.8GHz only has vertical polarization, other 4 working frequency range all have two kinds of polarization modes of horizontal and vertical, sensor scan swath is greater than 1600km, sensitivity is better than 0.8K in 37.0GHz frequency range, and all the other frequency ranges are better than 0.5K.Utilize this sensing data, can inverting obtain the parameter such as liquid water content, sea ice and rainfall amount in Global Sea-level temperature, Ocean Wind-field, atmospheric water steam content, cloud.
Existing method in the world adopts the bright temperature data acquisition of spaceborne radiometer 89GHz monitored area ice concentration, but not this wave band of HY-2 satellite microwave scanning radiometer; The ice concentration acquisition methods of the existing bright temperature data based on lower frequency wave band is all based on external spaceborne radiometer, can not be directly applied for HY-2 satellite.
Summary of the invention
The object of the present invention is to provide the ice concentration acquisition methods based on No. two satellite scanning microwave radiometers in ocean, to solve the above problems.
Ice concentration acquisition methods based on No. two satellite scanning microwave radiometers in ocean is provided in an embodiment of the present invention, comprise: according to first of characteristic feature region the bright temperature data, obtain bright wyntet's sign value, wherein said characteristic feature region comprises one year ice characteristic area, for many years ice characteristic area and seawater characteristics region, and described the first bright temperature data are obtained by No. two, ocean satellite scanning microwave radiometer; According to the second bright temperature data of described bright wyntet's sign value and monitored area, obtain the dynamic bright wyntet's sign value of described monitored area, wherein said the second bright temperature data are obtained by No. two, ocean satellite scanning microwave radiometer; According to the dynamic bright wyntet's sign value of described monitored area, obtain the dynamic ice concentration of described monitored area.
Preferably, describedly obtain bright wyntet's sign value according to first of characteristic feature region the bright temperature data, comprising: according to the bright temperature data of the 18.7GHz frequency range vertical polarization directions of described one year ice characteristic area, obtain bright wyntet's sign value T
b, FY, 18.7V; According to the bright temperature data of the 18.7GHz frequency range horizontal polarization directions of described one year ice characteristic area, obtain bright wyntet's sign value T
b, FY, 18.7H; According to the bright temperature data of the 37GHz frequency range vertical polarization directions of described one year ice characteristic area, obtain bright wyntet's sign value T
b, FY, 37V; According to the bright temperature data of the 18.7GHz frequency range vertical polarization directions of the described characteristic area of ice for many years, obtain bright wyntet's sign value T
b, MY, 18.7V; According to the bright temperature data of the 18.7GHz frequency range horizontal polarization directions of the described characteristic area of ice for many years, obtain bright wyntet's sign value T
b, MY, 18.7H; According to the bright temperature data of the 37GHz frequency range vertical polarization directions of the described characteristic area of ice for many years, obtain bright wyntet's sign value T
b, MY, 37V; According to the bright temperature data of the 18.7GHz frequency range vertical polarization directions in described seawater characteristics region, obtain bright wyntet's sign value T
b, OW, 18.7V; According to the bright temperature data of the 18.7GHz frequency range horizontal polarization directions in described seawater characteristics region, obtain bright wyntet's sign value T
b, OW, 18.7H; According to the bright temperature data of the 37GHz frequency range vertical polarization directions in described seawater characteristics region, obtain bright wyntet's sign value T
b, OW, 37V.
Preferably, according to the second bright temperature data of described bright wyntet's sign value and monitored area, obtain the dynamic bright wyntet's sign value of described monitored area, comprise: by the pixel of dividing in advance, obtain the bright temperature data of each described pixel in monitored area, and utilize the bright wyntet's sign value in characteristic feature region and the bright temperature data of each pixel to calculate respectively the corresponding one year ice closeness of each described pixel and ice concentration for many years; According to the described one year ice closeness of each pixel and the result of calculation of ice concentration for many years, described monitored area is divided into one year ice region, for many years territory, ice formation and seawater region; Utilize the bright temperature data of the default monitoring time in described one year ice region, described territory, ice formation for many years and described seawater region, obtain the corresponding bright wyntet's sign value of each described monitoring time, as the dynamic bright wyntet's sign value of described monitored area.
Preferably, describedly obtain the dynamic ice concentration of described monitored area according to the dynamic bright wyntet's sign value of described monitored area, comprising: utilize the corresponding bright wyntet's sign value of each described monitoring time, obtain the corresponding bright temperature coefficient M of each described monitoring time
i, F
iand D
i, i=0～3; Utilize the corresponding described bright wyntet's sign value of each described monitoring time and described bright temperature coefficient, obtain the ice concentration of monitored area described in each described monitoring time.
Preferably, the method also comprises: the described ice concentration of utilizing the 37GHz frequency range of described monitored area and the bright temperature Data correction of 18.7GHz frequency range to obtain.
Preferably, the described described ice concentration of utilizing the 37GHz frequency range of described monitored area and the bright temperature Data correction of 18.7GHz frequency range to obtain, comprising: the bright temperature data acquisition spectrum gradient rate GR (37/18.7) that utilizes 37GHz frequency range and the 18.7GHz frequency range of described monitored area; Judge whether described GR (37/18.7) >=first preset value is set up, if set up, the order sea ice density corresponding with it is zero.
Preferably, the method also comprises: the described ice concentration of utilizing the 23.8GHz frequency range of described monitored area and the bright temperature Data correction of 18.7GHz frequency range to obtain.
Preferably, the described described ice concentration of utilizing the 23.8GHz frequency range of described monitored area and the bright temperature Data correction of 18.7GHz frequency range to obtain, comprising: the bright temperature data acquisition spectrum gradient rate GR (23.8/18.7) that utilizes 23.8GHz frequency range and the 18.7GHz frequency range of described monitored area; Judge whether described GR (23.8/18.7) >=second preset value is set up, if set up, the order sea ice density corresponding with it is zero.
Preferably, the bright temperature data acquisition spectrum gradient rate GR (37/18.7) of described 37GHz frequency range and the 18.7GHz frequency range of utilizing described monitored area, comprising: utilize formula GR (37/18.7)=(T
b, 37v-T
b, 18.7v)/(T
b, 37v+ T
b, 18.7v) calculate described spectrum gradient rate GR (37/18.7), wherein T
b, 37Vfor the bright temperature data of described monitored area 37GHz frequency range vertical polarization directions; T
b, 18.7Vfor the bright temperature data of described monitored area 18.7GHz frequency range vertical polarization directions.
Preferably, the bright temperature data acquisition spectrum gradient rate GR (23.8/18.7) of described 23.8GHz frequency range and the 18.7GHz frequency range of utilizing described monitored area, comprising: utilize formula GR (23.8/18.7)=(T
b, 23.8v-T
b, 18.7v)/(T
b, 23.8v+ T
b, 18.7v) calculate described spectrum gradient rate GR (23.8/18.7), wherein T
b, 23.8vfor the bright temperature data of described monitored area 23.8GHz frequency range vertical polarization directions; T
b, 18.7Vfor the bright temperature data of described monitored area 18.7GHz frequency range vertical polarization directions.
The ice concentration acquisition methods based on No. two satellite scanning microwave radiometers in ocean that the embodiment of the present invention provides, utilize No. two, ocean satellite scanning microwave radiometer to obtain the first bright temperature data in characteristic feature region, and obtain the second bright temperature data of monitored area; Obtain bright wyntet's sign value according to these first bright temperature data; Further utilize the second bright temperature data and this bright wyntet's sign value to obtain dynamic bright wyntet's sign value, obtain thereby utilize behavioral characteristics value and bright temperature data to carry out monitored area ice concentration, overcome the technical matters of the ice concentration of the corresponding monitored area of data acquisition that cannot obtain based on No. two, ocean satellite in correlation technique.
Brief description of the drawings
Fig. 1 shows the process flow diagram of the ice concentration acquisition methods based on No. two satellite scanning microwave radiometers in ocean in the embodiment of the present invention;
Fig. 2 shows the schematic diagram of characteristic feature region chosen position in the embodiment of the present invention;
Fig. 3 shows the schematic diagram of the polarization gradient rate in the characteristic feature region of choosing in the embodiment of the present invention;
Fig. 4 shows the schematic diagram of the spectrum gradient rate in the characteristic feature region of choosing in the embodiment of the present invention;
Fig. 5 shows the scatter diagram of polarization gradient rate and the spectrum gradient rate in the characteristic feature region of choosing in the embodiment of the present invention;
Fig. 6 is the partial enlarged drawing of Fig. 5;
Fig. 7 shows the process flow diagram that obtains monitored area ice concentration in the embodiment of the present invention.
Embodiment
Also by reference to the accompanying drawings the present invention is described in further detail below by specific embodiment.
The embodiment of the present invention provides a kind of ice concentration acquisition methods based on No. two satellite scanning microwave radiometers in ocean, and as shown in Figure 1, main treatment step comprises:
Step S11: according to first of characteristic feature region the bright temperature data, obtain bright wyntet's sign value, wherein characteristic feature region comprises one year ice characteristic area, for many years ice characteristic area and seawater characteristics region, and the first bright temperature data are obtained by No. two, ocean satellite scanning microwave radiometer;
Step S12: according to second of bright wyntet's sign value and monitored area the bright temperature data, obtain the dynamic bright wyntet's sign value of monitored area, wherein the second bright temperature data are obtained by No. two, ocean satellite scanning microwave radiometer;
Step S13: according to the dynamic bright wyntet's sign value of monitored area, obtain the dynamic ice concentration of monitored area.
The ice concentration acquisition methods of the embodiment of the present invention can be directly applied for No. two, ocean satellite.
In existing ice concentration acquisition methods, be that bright wyntet's sign value based on fixing is carried out the inverting of ice concentration, do not consider the seasonal variety of the radiation characteristic of different Sea Ice Types.
In this method, adopt dynamic bright wyntet's sign value, by one year ice and for many years the seasonal variations of ice radiation characteristic take into account, make the ice concentration degree of accuracy that obtains higher, the user demand of more satisfied reality.
In the present invention, according to the concrete grammar of the bright wyntet's sign value of bright temperature data acquisition in selected characteristic feature region be: according to the bright temperature data of the 18.7GHz frequency range vertical polarization directions of one year ice characteristic area, obtain bright wyntet's sign value T
b, FY, 18.7V; According to the bright temperature data of the 18.7GHz frequency range horizontal polarization directions of one year ice characteristic area, obtain bright wyntet's sign value T
b, FY, 18.7H; According to the bright temperature data of the 37GHz frequency range vertical polarization directions of one year ice characteristic area, obtain bright wyntet's sign value T
b, FY, 37V; According to the bright temperature data of the 18.7GHz frequency range vertical polarization directions of ice characteristic area for many years, obtain bright wyntet's sign value T
b, MY, 18.7V; According to the bright temperature data of the 18.7GHz frequency range horizontal polarization directions of ice characteristic area for many years, obtain bright wyntet's sign value T
b, MY, 18.7H; According to the bright temperature data of the 37GHz frequency range vertical polarization directions of ice characteristic area for many years, obtain bright wyntet's sign value T
b, MY, 37V; According to the bright temperature data of the 18.7GHz frequency range vertical polarization directions in seawater characteristics region, obtain bright wyntet's sign value T
b, OW, 18.7V; According to the bright temperature data of the 18.7GHz frequency range horizontal polarization directions in seawater characteristics region, obtain bright wyntet's sign value T
b, OW, 18.7H; According to the bright temperature data of the 37GHz frequency range vertical polarization directions in seawater characteristics region, obtain bright wyntet's sign value T
b, OW, 37V.
Obtain after the bright wyntet's sign value in characteristic feature region, can utilize T
b, FY, 18.7V, T
b, FY, 18.7H, T
b, FY, 37V, T
b, MY, 18.7V, T
b, MY, 18.7H, T
b, MY, 37V, T
b, OW, 18.7V, T
b, OW, 18.7Hand T
b, OW, 37Vobtain bright temperature coefficient M
i, F
iand D
i, i=0～3; Wherein M
i, F
iand D
icomputing method be:
Wherein, 18.7GHz vertical polarization directions can be simplified shown as 18.7V; 18.7GHz horizontal polarization directions can be simplified shown as 18.7H; 37VGHz vertical polarization directions can be simplified shown as 37V.
In this method, choose the sea ice region with characteristic feature and observe, wherein, one year ice characteristic area, for many years ice characteristic area and seawater characteristics region can need to be chosen according to observation.
As shown in Figure 2, one year ice characteristic area can be chosen Chukchi Sea Partial Sea Area, ice characteristic area can be chosen marine site to the north of Greenland for many years, chooses Norwegian Sea Partial Sea Area without ice sea, and the boxed area that wherein in Fig. 2, MY (Multi-year) is corresponding refers to ice characteristic area for many years; The boxed area that FY (First year) is corresponding refers to one year ice characteristic area; The boxed area that Water is corresponding refers to seawater characteristics region.
Determine behind characteristic feature region, obtain the bright temperature data in characteristic feature region, if Fig. 3 is the annual variations in 2012 of above-mentioned three regional polarization gradient rates; As 4 being the annual variations of above-mentioned three region spectrum gradient rates 2012.Can find out winter and spring (May November-next year) on the Northern Hemisphere by Fig. 3 and Fig. 4, the polarization gradient rate in seawater characteristics region and spectrum gradient rate are than one year ice characteristic area and the relevant parameter of ice characteristic area is high for many years, and the polarization gradient rate of one year ice characteristic area and for many years ice characteristic area is not obviously distinguished; The spectrum gradient rate of one year ice characteristic area will be higher than the spectrum gradient rate of ice characteristic area for many years; In two seasons of autumn in summer, along with the thawing of one year ice, its polarization gradient rate and spectrum gradient rate raise gradually, until 8～September, one year ice all melts as seawater, and for many years these two parameters of ice also because sea ice surface temperature raises and sea ice thawing presents growth in various degree.
For 3 marine sites choosing, draw the polarization gradient rate of monthly average every month and the scatter diagram of spectrum gradient rate, as shown in Figure 5, amplify the part that Fig. 6 is Fig. 5.As can be seen from Figure 5, the gradient rate parameter in seawater characteristics region is more concentrated, so choose month for annual; One year ice characteristic area is since June, two gradient rate parameters all start to rise, reach maximal value to September, then start to reduce, it is mainly because one year ice thawing in summer causes gradient rate parameter generation respective change, so the feature month of one year ice choose other month of removing for 6～October; As can be seen from Figure 5 6～August for many years the parameter of ice characteristic area and the parameter differences of one year ice characteristic area less, together with one year ice feature month, point " mixed ", thus choose in the feature month of ice characteristic area for many years remove 6,7, other month in August.According to the feature month of choosing, calculate 18.7V, the bright temperature mean value of 18.7H and 37V, concrete bright wyntet's sign value and choose month as table 1.
Three kinds of typical sea ice regions of table 1 are in the bright wyntet's sign value of 18.7V, 18.7H and 37V wave band
One year ice characteristic area, for many years ice characteristic area and the seawater characteristics region that utilization is chosen obtains the above-mentioned formula (1) of bright wyntet's sign value substitution and formula (2) obtains M
i, F
iand D
i, i=0～3.
Obtaining after the bright wyntet's sign value in characteristic feature region, can, according to the second bright temperature data of the bright wyntet's sign value in the characteristic feature region obtaining and monitored area, obtain the dynamic bright wyntet's sign value of monitored area, concrete grammar, as shown in Figure 7, key step comprises:
Step S21: by the pixel of dividing in advance, obtain the bright temperature data of each pixel in monitored area, and utilize the bright wyntet's sign value in characteristic feature region and the bright temperature data of each pixel to calculate respectively the corresponding one year ice closeness of each pixel and ice concentration for many years;
Step S22: according to the one year ice closeness of each pixel and the result of calculation of ice concentration for many years, monitored area is divided into one year ice region, for many years territory, ice formation and seawater region;
Step S23: utilize the bright temperature data of the default monitoring time in one year ice region, for many years territory, ice formation and seawater region, obtain the corresponding bright wyntet's sign value of each monitoring time, as the dynamic bright wyntet's sign value of monitored area;
Step S24: utilize the corresponding bright wyntet's sign value of each monitoring time, obtain the corresponding bright temperature coefficient M of each monitoring time
i, F
iand D
i, i=0～3;
Step S25: utilize the corresponding bright wyntet's sign value of each monitoring time and bright temperature coefficient, obtain the ice concentration of each monitoring time monitored area.
In step S21, utilize the bright wyntet's sign value in characteristic feature region and the bright temperature data of each pixel calculate respectively the corresponding one year ice closeness of each pixel and for many years the method for ice concentration be: utilize the bright wyntet's sign value in characteristic feature region to calculate bright temperature coefficient according to formula (1), (2); The bright temperature data of this bright temperature coefficient and each pixel are updated in the computing formula of ice concentration, calculate the one year ice closeness of each pixel and ice concentration for many years.
In this method, the computing method of ice concentration are: utilize formula
Calculate one year ice closeness;
Utilize formula
Calculate ice concentration for many years;
Ice concentration is C
t=C
fY+ C
mY, wherein D=D
0+ D
1* PR+D
2* GR+D
3* PR*GR, in the computing formula of above-mentioned ice concentration, the polarization gradient rate that PR is monitored area, the spectrum gradient rate that GR is monitored area.
Utilize formula PR (18.7)=(T
b, 18.7V-T
b, 18.7H)/(T
b, 18.7V+ T
b, 18.7H) calculate the polarization gradient rate of monitored area, wherein T
b, 18.7Vfor the bright temperature data of monitored area 18.7GHz frequency range vertical polarization directions, T
b, 18.7Hfor the bright temperature data of monitored area 18.7GHz frequency range horizontal polarization directions.
Utilize formula GR (37/18.7)=(T
b, 37v-T
b, 18.7v)/(T
b, 37v+ T
b, 18.7v) calculate the spectrum gradient rate of monitored area, wherein T
b, 37Vfor the bright temperature data of monitored area 37GHz frequency range vertical polarization directions; T
b, 18.7Vfor the bright temperature data of monitored area 18.7GHz frequency range vertical polarization directions.
Can find out according to the computing formula of above-mentioned ice concentration, obtain can obtaining after the bright wyntet's sign value in characteristic feature region and the second bright temperature data of monitored area the initial ice concentration of monitored area; But in the actual acquisition process of ice concentration, can find out that from Fig. 3 and Fig. 4 one year ice and spectrum gradient rate and the meeting of polarization gradient rate of ice for many years change because of the variation in season, make thus the ice concentration of the monitored area obtaining have larger error, for reducing the ice concentration error causing due to seasonal effect, the initial ice concentration also further calculating according to monitored area in this method is carried out the division in one year ice region, for many years territory, ice formation and seawater region to monitored area; For one year ice region, for many years territory, ice formation and seawater region after dividing, obtain the bright temperature data in each region according to monitoring time, obtain the corresponding bright wyntet's sign value of each monitoring time, obtain the time dependent dynamic bright wyntet's sign value in monitored area.
The above-mentioned concrete mode of obtaining dynamic bright wyntet's sign value is: obtain the one year ice closeness of the each pixel in monitored area and for many years after ice concentration according to the bright wyntet's sign value in characteristic feature region, based on the initial result of this closeness, one year ice threshold range and ice threshold range are for many years set, for example one year ice closeness is greater than to 90% region and is divided into one year ice region, ice concentration is for many years greater than to 90% region and is divided into territory, ice formation for many years, monitored area except one year ice region and for many years the region territory, ice formation be divided into seawater region.For the bright temperature data of repartitioning the one year ice region, for many years territory, ice formation and the seawater region that obtain and can obtain respectively each monitoring time, can recalculate and obtain the bright wyntet's sign value corresponding with each monitoring time section with the bright temperature data under 37V according to the 18.7V of each monitoring time, 18.7H, be i.e. the dynamic bright wyntet's sign value corresponding with each monitoring time section.
It should be noted that, above-mentioned each monitoring time can refer to monitoring time point, also can refer to monitoring time section.
For monitored area, obtain after dynamic bright temperature value, can utilize the bright temperature data of obtaining in each monitoring time section again to obtain polarization gradient rate and the spectrum gradient rate in this monitoring time section; Utilize the polarization gradient rate corresponding with each monitoring time section, spectrum gradient rate and dynamic bright wyntet's sign value to obtain the one year ice closeness in the one year ice region in this monitoring time section and the ice concentration for many years in territory, ice formation for many years; Utilize one year ice closeness and for many years ice concentration obtain corresponding ice concentration.
In the acquisition methods of the ice concentration of the embodiment of the present invention, take into full account one year ice, there is the feature of seasonal variety in ice radiation characteristic for many years, obtain the ice concentration of monitored area based on dynamic bright wyntet's sign value, the precision that improves ice concentration, provides technical support for business provides ice concentration.
Utilize the ice concentration that said method obtains to have wrong situation, for example, can calculate the sea ice of low closeness at exposed waters, this kind of situation is mainly because the phenomenons such as water vapor, cloud liquid water, rainfall in atmosphere cause.In addition, at fringe of land, be subject to that land ratio radiance is higher to be affected, bright temperature is higher, easily causes that ice concentration result makes a mistake in fringe of land region.
For the mistake of effectively avoiding ice concentration result to occur, the ice concentration of utilizing the 37GHz frequency range of monitored area and the bright temperature Data correction of 18.7GHz frequency range to obtain in this method.
Concrete methods of realizing is: the ice concentration of utilizing the 37GHz frequency range of monitored area and the bright temperature Data correction of 18.7GHz frequency range to obtain, comprising: the bright temperature data acquisition spectrum gradient rate GR (37/18.7) that utilizes 37GHz frequency range and the 18.7GHz frequency range of monitored area; Judge whether GR (37/18.7) >=first preset value is set up, if set up, the order sea ice density corresponding with it is zero, utilizes this kind of method can effectively remove the situation due to the higher caused ice concentration result erroneous judgement of cloud liquid water content.
Wherein the first above-mentioned preset value can arrange according to actual needs, preferably, provides example in this method, and making the first preset value is 0.13, when GR (37/18.7)>=0.13, makes ice concentration C
t=0.
The ice concentration of utilizing the 23.8GHz frequency range of monitored area and the bright temperature Data correction of 18.7GHz frequency range to obtain in this method.
Concrete methods of realizing is: the ice concentration of utilizing the 23.8GHz frequency range of monitored area and the bright temperature Data correction of 18.7GHz frequency range to obtain, comprising: the bright temperature data acquisition spectrum gradient rate GR (23.8/18.7) that utilizes 23.8GHz frequency range and the 18.7GHz frequency range of monitored area; Judge whether GR (23.8/18.7) >=second preset value is set up, if set up, the order sea ice density corresponding with it is zero, utilizes the method can effectively remove the situation due to the caused ice concentration result erroneous judgement of the phenomenon such as water vapor and rainfall in atmosphere; And can effectively avoid at fringe of land, owing to being subject to land ratio radiance to affect compared with high the situation that the ice concentration result that causes makes a mistake.
In the method, the second preset value can arrange according to actual needs, preferably, provides example in this method, and making the second preset value is 0.085, when GR (23.8/18.7)>=0.085, makes ice concentration C
t=0.
Wherein GR (23.8/18.7)=(T
b, 23.8V-T
b, 18.7V)/(T
b, 23.8V+ T
b, 18.7V), for the bright temperature of the vertical polarization mode of 23.8GHz and 18.7GHz is calculated spectrum gradient rate.
In the acquisition methods of ice concentration of the present invention, utilize dynamic bright wyntet's sign value, by one year ice and for many years the seasonal variations of ice radiation characteristic take into account; And calculating in the process of dynamic bright wyntet's sign value, rejected the impact of the phenomenons such as water vapor, cloud liquid water, rainfall in atmosphere, make the bright wyntet's sign statistics of different extra large table section more accurate.
These are only the preferred embodiments of the present invention, be not limited to the present invention, for a person skilled in the art, the present invention can have various modifications and variations.Within the spirit and principles in the present invention all, any amendment of doing, be equal to replacement, improvement etc., within all should being included in protection scope of the present invention.
Claims (10)
1. the ice concentration acquisition methods based on No. two satellite scanning microwave radiometers in ocean, is characterized in that, comprising:
According to first of characteristic feature region the bright temperature data, obtain bright wyntet's sign value, wherein said characteristic feature region comprises one year ice characteristic area, for many years ice characteristic area and seawater characteristics region, and described the first bright temperature data are obtained by No. two, ocean satellite scanning microwave radiometer;
According to the second bright temperature data of described bright wyntet's sign value and monitored area, obtain the dynamic bright wyntet's sign value of described monitored area, wherein said the second bright temperature data are obtained by No. two, ocean satellite scanning microwave radiometer;
According to the dynamic bright wyntet's sign value of described monitored area, obtain the dynamic ice concentration of described monitored area.
2. method according to claim 1, is characterized in that, describedly obtains bright wyntet's sign value according to first of characteristic feature region the bright temperature data, comprising:
According to the bright temperature data of the 18.7GHz frequency range vertical polarization directions of described one year ice characteristic area, obtain bright wyntet's sign value T
b, FY, 18.7V;
According to the bright temperature data of the 18.7GHz frequency range horizontal polarization directions of described one year ice characteristic area, obtain bright wyntet's sign value T
b, FY, 18.7H;
According to the bright temperature data of the 37GHz frequency range vertical polarization directions of described one year ice characteristic area, obtain bright wyntet's sign value T
b, FY, 37V;
According to the bright temperature data of the 18.7GHz frequency range vertical polarization directions of the described characteristic area of ice for many years, obtain bright wyntet's sign value T
b, MY, 18.7V;
According to the bright temperature data of the 18.7GHz frequency range horizontal polarization directions of the described characteristic area of ice for many years, obtain bright wyntet's sign value T
b, MY, 18.7H;
According to the bright temperature data of the 37GHz frequency range vertical polarization directions of the described characteristic area of ice for many years, obtain bright wyntet's sign value T
b, MY, 37V;
According to the bright temperature data of the 18.7GHz frequency range vertical polarization directions in described seawater characteristics region, obtain bright wyntet's sign value T
b, OW, 18.7V; According to the bright temperature data of the 18.7GHz frequency range horizontal polarization directions in described seawater characteristics region, obtain bright wyntet's sign value T
b, OW, 18.7H;
According to the bright temperature data of the 37GHz frequency range vertical polarization directions in described seawater characteristics region, obtain bright wyntet's sign value T
b, OW, 37V.
3. method according to claim 1 and 2, is characterized in that, according to the second bright temperature data of described bright wyntet's sign value and monitored area, obtains the dynamic bright wyntet's sign value of described monitored area, comprising:
By the pixel of dividing in advance, obtain the bright temperature data of each described pixel in monitored area, and utilize the bright wyntet's sign value in characteristic feature region and the bright temperature data of each pixel to calculate respectively the corresponding one year ice closeness of each described pixel and ice concentration for many years;
According to the described one year ice closeness of each pixel and the result of calculation of ice concentration for many years, described monitored area is divided into one year ice region, for many years territory, ice formation and seawater region;
Utilize the bright temperature data of the default monitoring time in described one year ice region, described territory, ice formation for many years and described seawater region, obtain the corresponding bright wyntet's sign value of each described monitoring time, as the dynamic bright wyntet's sign value of described monitored area.
4. method according to claim 3, is characterized in that, describedly obtains the dynamic ice concentration of described monitored area according to the dynamic bright wyntet's sign value of described monitored area, comprising:
Utilize the corresponding bright wyntet's sign value of each described monitoring time, obtain the corresponding bright temperature coefficient M of each described monitoring time
i, F
iand D
i, i=0～3;
Utilize the corresponding described bright wyntet's sign value of each described monitoring time and described bright temperature coefficient, obtain the ice concentration of monitored area described in each described monitoring time.
5. method according to claim 4, is characterized in that, the method also comprises: the described ice concentration of utilizing the 37GHz frequency range of described monitored area and the bright temperature Data correction of 18.7GHz frequency range to obtain.
6. method according to claim 5, is characterized in that, the described described ice concentration of utilizing the 37GHz frequency range of described monitored area and the bright temperature Data correction of 18.7GHz frequency range to obtain, comprising:
Utilize the bright temperature data acquisition spectrum gradient rate GR (37/18.7) of 37GHz frequency range and the 18.7GHz frequency range of described monitored area;
Judge whether described GR (37/18.7) >=first preset value is set up, if set up, the order sea ice density corresponding with it is zero.
7. method according to claim 4, is characterized in that, the method also comprises: the described ice concentration of utilizing the 23.8GHz frequency range of described monitored area and the bright temperature Data correction of 18.7GHz frequency range to obtain.
8. method according to claim 7, is characterized in that, the described described ice concentration of utilizing the 23.8GHz frequency range of described monitored area and the bright temperature Data correction of 18.7GHz frequency range to obtain, comprising:
Utilize the bright temperature data acquisition spectrum gradient rate GR (23.8/18.7) of 23.8GHz frequency range and the 18.7GHz frequency range of described monitored area;
Judge whether described GR (23.8/18.7) >=second preset value is set up, if set up, the order sea ice density corresponding with it is zero.
9. method according to claim 6, is characterized in that, the bright temperature data acquisition spectrum gradient rate GR (37/18.7) of described 37GHz frequency range and the 18.7GHz frequency range of utilizing described monitored area, comprising:
Utilize formula GR (37/18.7)=(T
b, 37v-T
b, 18.7v)/(T
b, 37v+ T
b, 18.7v) calculate described spectrum gradient rate GR (37/18.7), wherein T
b, 37Vfor the bright temperature data of described monitored area 37GHz frequency range vertical polarization directions; T
b, 18.7Vfor the bright temperature data of described monitored area 18.7GHz frequency range vertical polarization directions.
10. method according to claim 8, is characterized in that, the bright temperature data acquisition spectrum gradient rate GR (23.8/18.7) of described 23.8GHz frequency range and the 18.7GHz frequency range of utilizing described monitored area, comprising:
Utilize formula GR (23.8/18.7)=(T
b, 23.8v-T
b, 18.7v)/(T
b, 23.8v+ T
b, 18.7v) calculate described spectrum gradient rate GR (23.8/18.7), wherein T
b, 23.8vfor the bright temperature data of described monitored area 23.8GHz frequency range vertical polarization directions; T
b, 18.7Vfor the bright temperature data of described monitored area 18.7GHz frequency range vertical polarization directions.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410498638.2A CN104198052B (en) | 2014-09-25 | 2014-09-25 | Ice concentration acquisition methods based on the satellite scanning microwave radiometer of ocean two |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410498638.2A CN104198052B (en) | 2014-09-25 | 2014-09-25 | Ice concentration acquisition methods based on the satellite scanning microwave radiometer of ocean two |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104198052A true CN104198052A (en) | 2014-12-10 |
| CN104198052B CN104198052B (en) | 2017-07-14 |
Family
ID=52083380
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410498638.2A Active CN104198052B (en) | 2014-09-25 | 2014-09-25 | Ice concentration acquisition methods based on the satellite scanning microwave radiometer of ocean two |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104198052B (en) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105737806A (en) * | 2016-02-16 | 2016-07-06 | 国家卫星海洋应用中心 | Method and device for inverting sea surface salinity based on satellite-borne microwave radiometer |
| CN106777246A (en) * | 2016-12-27 | 2017-05-31 | 大连海事大学 | There is the extracting method of frozen water domain course line navigation window |
| CN106872466A (en) * | 2016-12-31 | 2017-06-20 | 中国科学院遥感与数字地球研究所 | A kind of lake ice freeze thawing monitoring method and system based on dynamic Decomposition of Mixed Pixels method |
| CN106871877A (en) * | 2017-02-13 | 2017-06-20 | 国家卫星海洋应用中心 | Sea ice mark determines method and device |
| CN107886473A (en) * | 2017-11-09 | 2018-04-06 | 河南工业大学 | From the method for FY 3MWRI data inversion Sea Ice Model closenesses |
| CN108106734A (en) * | 2017-12-07 | 2018-06-01 | 中国科学院海洋研究所 | One kind seeks calculation optimization algorithm towards IDL language ice concentrations |
| CN109520451A (en) * | 2018-11-13 | 2019-03-26 | 南京大学 | A kind of remote sensing estimation method of sea ice surface roughness |
| CN112966656A (en) * | 2021-03-29 | 2021-06-15 | 国家卫星海洋应用中心 | Data processing method and device |
| CN113093183A (en) * | 2021-04-02 | 2021-07-09 | 国家卫星海洋应用中心 | Threshold determination method, sea ice thickness inversion method, device, equipment and medium |
| CN113127794A (en) * | 2021-04-07 | 2021-07-16 | 中山大学 | Method for calculating density of arctic sea ice |
| CN113610189A (en) * | 2021-08-23 | 2021-11-05 | 国家卫星海洋应用中心 | Remote sensing classification method and device for sea ice and sea water and electronic equipment |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6137437A (en) * | 1999-03-24 | 2000-10-24 | Agence Spatiale Europeenne | Spaceborne scatterometer |
| CN1908695A (en) * | 2005-08-01 | 2007-02-07 | 中国海洋大学 | Dynamic inserted technique for sea ice satellite remote sensing data in Bohai cloud coverage area |
| CN101424741A (en) * | 2008-12-08 | 2009-05-06 | 中国海洋大学 | Real time extracting method for satellite remote sensing sea fog characteristic quantity |
| CN101587190A (en) * | 2009-06-25 | 2009-11-25 | 国家海洋局第二海洋研究所 | Satellite remote-sensing monitoring method of daytime marine fog |
| US20100063733A1 (en) * | 2008-09-09 | 2010-03-11 | Thomas Patrick Yunck | Cellular Interferometer for Continuous Earth Remote Observation (CICERO) |
| CN101799561A (en) * | 2010-02-05 | 2010-08-11 | 民政部国家减灾中心 | Snow disaster remote sensing monitoring simulation evaluation method based on disaster reduction small satellite |
-
2014
- 2014-09-25 CN CN201410498638.2A patent/CN104198052B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6137437A (en) * | 1999-03-24 | 2000-10-24 | Agence Spatiale Europeenne | Spaceborne scatterometer |
| CN1908695A (en) * | 2005-08-01 | 2007-02-07 | 中国海洋大学 | Dynamic inserted technique for sea ice satellite remote sensing data in Bohai cloud coverage area |
| US20100063733A1 (en) * | 2008-09-09 | 2010-03-11 | Thomas Patrick Yunck | Cellular Interferometer for Continuous Earth Remote Observation (CICERO) |
| CN101424741A (en) * | 2008-12-08 | 2009-05-06 | 中国海洋大学 | Real time extracting method for satellite remote sensing sea fog characteristic quantity |
| CN101587190A (en) * | 2009-06-25 | 2009-11-25 | 国家海洋局第二海洋研究所 | Satellite remote-sensing monitoring method of daytime marine fog |
| CN101799561A (en) * | 2010-02-05 | 2010-08-11 | 民政部国家减灾中心 | Snow disaster remote sensing monitoring simulation evaluation method based on disaster reduction small satellite |
Non-Patent Citations (2)
| Title |
|---|
| LIJIAN SHI: "SEA ICE THICKNESS RETRIEVAL FROM SAR IMAGRY OVER BOHAI SEA", 《GEOSCIENCE AND REMOTE SENSING SYMPOSIUM,2014 IEEE INTERNATIONAL》 * |
| 郝春江 等: "利用卫星遥感监测南极海冰的研究和进展", 《海洋预报》 * |
Cited By (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105737806A (en) * | 2016-02-16 | 2016-07-06 | 国家卫星海洋应用中心 | Method and device for inverting sea surface salinity based on satellite-borne microwave radiometer |
| CN105737806B (en) * | 2016-02-16 | 2018-08-17 | 国家卫星海洋应用中心 | Method and apparatus based on satellite-borne microwave radiometer inverting sea surface salinity |
| CN106777246A (en) * | 2016-12-27 | 2017-05-31 | 大连海事大学 | There is the extracting method of frozen water domain course line navigation window |
| CN106777246B (en) * | 2016-12-27 | 2020-12-11 | 大连海事大学 | Method for extracting navigation window of air route with ice water area |
| CN106872466B (en) * | 2016-12-31 | 2019-08-13 | 中国科学院遥感与数字地球研究所 | A kind of lake ice freeze thawing monitoring method and system based on dynamic Decomposition of Mixed Pixels method |
| CN106872466A (en) * | 2016-12-31 | 2017-06-20 | 中国科学院遥感与数字地球研究所 | A kind of lake ice freeze thawing monitoring method and system based on dynamic Decomposition of Mixed Pixels method |
| CN106871877A (en) * | 2017-02-13 | 2017-06-20 | 国家卫星海洋应用中心 | Sea ice mark determines method and device |
| CN106871877B (en) * | 2017-02-13 | 2018-08-28 | 国家卫星海洋应用中心 | Sea ice mark determines method and device |
| CN107886473A (en) * | 2017-11-09 | 2018-04-06 | 河南工业大学 | From the method for FY 3MWRI data inversion Sea Ice Model closenesses |
| CN107886473B (en) * | 2017-11-09 | 2020-12-11 | 河南工业大学 | Method for inverting north sea ice concentration from FY-3MWRI data |
| CN108106734B (en) * | 2017-12-07 | 2019-08-02 | 中国科学院海洋研究所 | One kind seeking calculation optimization algorithm towards IDL language ice concentration |
| CN108106734A (en) * | 2017-12-07 | 2018-06-01 | 中国科学院海洋研究所 | One kind seeks calculation optimization algorithm towards IDL language ice concentrations |
| CN109520451A (en) * | 2018-11-13 | 2019-03-26 | 南京大学 | A kind of remote sensing estimation method of sea ice surface roughness |
| CN112966656A (en) * | 2021-03-29 | 2021-06-15 | 国家卫星海洋应用中心 | Data processing method and device |
| CN113093183A (en) * | 2021-04-02 | 2021-07-09 | 国家卫星海洋应用中心 | Threshold determination method, sea ice thickness inversion method, device, equipment and medium |
| CN113093183B (en) * | 2021-04-02 | 2022-03-08 | 国家卫星海洋应用中心 | Threshold determination method, sea ice thickness inversion method, device, equipment and medium |
| CN113127794A (en) * | 2021-04-07 | 2021-07-16 | 中山大学 | Method for calculating density of arctic sea ice |
| CN113127794B (en) * | 2021-04-07 | 2023-02-03 | 中山大学 | Method for calculating density of arctic sea ice |
| CN113610189A (en) * | 2021-08-23 | 2021-11-05 | 国家卫星海洋应用中心 | Remote sensing classification method and device for sea ice and sea water and electronic equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104198052B (en) | 2017-07-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104198052A (en) | Acquisition method for sea ice concentration on basis of ocean No. II satellite scanning microwave radiometer | |
| Gómez-Gesteira et al. | The state of climate in NW Iberia | |
| Du et al. | Satellite and Argo observed surface salinity variations in the tropical Indian Ocean and their association with the Indian Ocean dipole mode | |
| Siddique-E-Akbor et al. | Satellite precipitation data–driven hydrological modeling for water resources management in the Ganges, Brahmaputra, and Meghna Basins | |
| Moawad | Analysis of the flash flood occurred on 18 January 2010 in wadi El Arish, Egypt (a case study) | |
| CN106908415A (en) | A kind of big region crops time of infertility Soil Moisture Monitoring method based on amendment NDVI time serieses | |
| CN102722766B (en) | Wheat output predication method based on revised regional climate mode data | |
| Hao et al. | Oasis cold island effect and its influence on air temperature: a case study of Tarim Basin, Northwest China | |
| Qi et al. | Climate change on the southern slope of Mt. Qomolangma (Everest) Region in Nepal since 1971 | |
| Zecchetto et al. | Diagnostics of T1279 ECMWF analysis winds in the Mediterranean Basin by comparison with ASCAT 12.5 km winds | |
| CN107886473A (en) | From the method for FY 3MWRI data inversion Sea Ice Model closenesses | |
| Guo et al. | Impact of lake/reservoir expansion and shrinkage on energy and water vapor fluxes in the surrounding area | |
| Ashkenazy et al. | On the probability and spatial distribution of ocean surface currents | |
| Singh et al. | Flood monitoring using microwave passive remote sensing (AMSR-E) in part of the Brahmaputra basin, India | |
| Heathfield et al. | Erosive water level regime and climatic variability forcing of beach–dune systems on south‐western Vancouver Island, British Columbia, Canada | |
| Kayser et al. | Assessing geeSEBAL automated calibration and meteorological reanalysis uncertainties to estimate evapotranspiration in subtropical humid climates | |
| Nunez et al. | Rainfall estimation in south-west Tasmania using satellite images and phytosociological calibration | |
| CN111597692A (en) | Surface net radiation estimation method, system, electronic equipment and storage medium | |
| Beaton et al. | Investigating snowpack across scale in the northern Great Lakes–St. Lawrence forest region of Central Ontario, Canada | |
| Hanslow et al. | Rising tides: tidal inundation in south east Australian estuaries | |
| Wang et al. | The relationship between sea level change of China's coast and ENSO | |
| Siyal et al. | Study of soil, water, and cropping pattern in Danastar Wah (Manchar Lake) command area using geospatial tools | |
| Koffi et al. | The water vapour flux above Switzerland and its role in the August 2005 extreme precipitation and flooding | |
| Ashraf et al. | Remote sensing of the glacial environ‐ment influenced by climate change | |
| Masina et al. | Variations in monthly maximum gust speed at St Mary's, Isles of Scilly (UK) |
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 |