CN107886473A - From the method for FY 3MWRI data inversion Sea Ice Model closenesses - Google Patents
From the method for FY 3MWRI data inversion Sea Ice Model closenesses Download PDFInfo
- Publication number
- CN107886473A CN107886473A CN201711101311.7A CN201711101311A CN107886473A CN 107886473 A CN107886473 A CN 107886473A CN 201711101311 A CN201711101311 A CN 201711101311A CN 107886473 A CN107886473 A CN 107886473A
- Authority
- CN
- China
- Prior art keywords
- mtd
- msub
- mtr
- data
- mrow
- 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
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000012937 correction Methods 0.000 claims abstract description 21
- 241001269238 Data Species 0.000 claims abstract description 11
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 11
- 238000001914 filtration Methods 0.000 claims abstract description 4
- 230000010287 polarization Effects 0.000 claims description 33
- 238000004364 calculation method Methods 0.000 claims description 5
- 102220486162 Alkaline ceramidase 1_T89V_mutation Human genes 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 claims description 3
- 235000013399 edible fruits Nutrition 0.000 claims 1
- 238000011160 research Methods 0.000 abstract description 9
- 239000013535 sea water Substances 0.000 description 2
- 230000001932 seasonal effect Effects 0.000 description 2
- 240000005636 Dryobalanops aromatica Species 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4038—Image mosaicing, e.g. composing plane images from plane sub-images
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V9/00—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2200/00—Indexing scheme for image data processing or generation, in general
- G06T2200/32—Indexing scheme for image data processing or generation, in general involving image mosaicing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Probability & Statistics with Applications (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Image Processing (AREA)
Abstract
The present invention relates to a kind of method from the MWRI data inversion Sea Ice Model closenesses of FY 3, step includes:The MWRI level one datas of FY 3 are obtained, data are carried out with radiant correction, data splicing, land mask process;Choose marine site and determine pure water, pure ice value, be fitted ice concentration, determine the expression formula of ice concentration, ice concentration result is finally given by filtering process twice;The present invention's can be by the domestic MWRI data applications of FY 3 into ASI algorithms, inverting research is carried out to Sea Ice Model closeness, laid a good foundation using the domestic MWRI data inversion Sea Ice Model closenesses of FY 3, research for arctic navigation channel provides technical guarantee and data are supported, improves the efficient application of domestic satellite remote sensing date.
Description
Technical field
The invention belongs to satellite remote sensing technology field, more particularly to one kind are close from FY-3MWRI data inversion Sea Ice Models
The method of intensity.
Background technology
FY-3MWRI (satellite of wind and cloud 3) is the polar orbiting meteorological satellite of new generation of China's independent research.With China in recent years
The rapid development of aerospace industry, the remote sensing satellite of China's transmitting is more and more, wherein than more typical such as wind and cloud series, high score
Series, environment series etc..But just at present, have one it is very prominent the drawbacks of be exactly existing these satellite datas in China profit
It is very insufficient.Trace it to its cause and be, it is also more ripe although the inversion algorithm at present on ice concentration is a lot,
Due to being had differences between satellite data, so original algorithm is each to need to redefine in the new data of utilization
The value of coefficient, and the value of each mooring points value, domestic satellite remote sensing date can not be applied efficiently.
The content of the invention
Present invention aims to overcome that the deficiencies in the prior art and provide a kind of from FY-3MWRI data inversions north
The method of pole ice concentration.
The object of the present invention is achieved like this:A kind of method from FY-3MWRI data inversion Sea Ice Model closenesses,
It is characterized in that:Comprise the following steps:
Step 1), obtain FY-3MWRI level one datas.
Step 2), the FY-3MWRI level one datas to acquisition pre-process;It is divided into radiant correction, data splicing, land
Mask;
Radiant correction:Using radiant correction formula L=Gain*DN+Bias, wherein L is the bright temperature value after radiant correction,
Gain and Bias is radiant correction parameter, and DN is pixel value in the level one data that FY-3MWRI is provided, and is carried according to FY-3MWRI
The data set of confession understands Gain=-0.01, Bias=0;It is public that the FY-3MWRI level one datas of acquisition are substituted into radiant correction
The data after correction are acquired in formula;
Data are spliced:The FY-3MWRI data of acquisition are single track imaging, image mosaic are carried out using ENVI softwares, by FY-
Longitude and latitude file in 3MWRI data sets is imported into ENVI, and data are carried out by building the GLT files with coordinate system
Splicing;
Land mask:Land mask is handled by the band math in ENVI softwares;
Step 3), two marine sites are chosen, choose the selection of the data progress sample point of 1 year, calculate two regions respectively
Polarization difference, carry out probability distribution statistical, choose the numerical value that daily maximum probability occurs, as the polarization difference value on the same day,
Then the average value of annual polarization difference value is taken;Determine the mooring points value of pure water and pure ice;P is designated as respectively0(pure water), P1It is (pure
Ice), P=T89V-T89H, wherein T89V be 89GHz frequency ranges vertical polarization, T89H be 89GHz frequency ranges horizontal polarization, P
Polarization for all kinds of pixels is poor;
Step 4), the P by acquisition0、P1Value, with interpolation calculation, the parameter of ASI ice concentration inversion formulas is determined,
ASI ice concentration inversion algorithms are specific as follows:
P=Tbv-Tbh (1)
Wherein, TbvIt is the bright temperature of vertical polarization, TbhIt is the bright temperature of horizontal polarization;
Select ice concentration of the three rank multinomials fitting from 0% to 100% as follows:
C=d3p3+d2p2+d1p+d0 (2)
Step 5), the P by acquisition0、P1Value is updated in formula (2), then again to formula (2) derivation, the polarization difference of ice face
Significantly less than the polarization difference of water, polarization difference P difference P during ice concentration C convergences 0 and 10And P1, draw and be used for
The quaternary once linear equation group of solution formula (2) coefficient, is shown in formula (3), and d is calculated using formula (3)0、d1、d2、d3
By d0、d1、d2、d3It is brought into formula (2) and draws ice concentration C.
Finally determination ice concentration C expression formula is:
C=6.45714 × 10-6P3-6.05256×10-4P2-9.22521×10-3P+1.10031 (4)
P is obtained by step 4)0、P1Value, and then obtain ice concentration formula provided by the invention:
C=1.29 × 10-5P3-1.28×10-3P2+1.01×10-2P+1.02 (5)
Ice concentration result based on the line number evidences of FY-3MWRI mono- is finally inversed by by formula (5);
Step 6), the method based on image entropy, it is determined that the weather filter threshold value based on FY-3MWRI data;Pass through
MATLAB majorized function fmin-search draws the threshold values of weather filter;Sea ice is finally given by filtering process twice
Closeness result.
The present invention's can be carried out anti-by domestic FY-3MWRI data applications into ASI algorithms to Sea Ice Model closeness
Drill research.But it can be varied from because ASI algorithms are directed to different data core parameters, so the present invention, which passes through, chooses typical case
The sample point in marine site, it is determined that the mooring points value of pure water and pure ice, and then the calculation formula of ice concentration is determined.And because
Often influenceed more seriously by weather using high-frequency data, so needing weather filter to be handled, because weather filters
The threshold value of device can also change with the change of data, so the method for the invention based on image entropy, it is determined that based on FY-3MWRI
The weather filter threshold value of data.The present invention is to have established base using domestic FY-3MWRI data inversions Sea Ice Model closeness
Plinth, the research of propulsion and arctic navigation channel for China polar region strategy provides technical guarantee and data are supported.Improve domestic
The efficient application of satellite remote sensing date.
Brief description of the drawings
Fig. 1 is FY-3MWRI Arctic Ocean ice concentration results on January 3rd, 2016.
Fig. 2 is SSM/I Arctic Ocean ice concentration results on January 3rd, 2016.
Embodiment
Embodiment 1, a kind of method from FY-3MWRI data inversion Sea Ice Model closenesses, it is characterised in that:Including with
Lower step:
Step 1), FY-3MWRI level one datas are obtained, FY-3MWRI TMI Major Systems parameters are as shown in table 1.
The FY-3MWRI TMI Major Systems parameters of table 1
Step 2), the FY-3MWRI level one datas to acquisition pre-process;It is divided into radiant correction, data splicing, land
Mask;
Radiant correction:The data that FY-3MWRI is provided are level one data, and data value is DN values, are anticipated without specific physics
Justice, therefore in inverting ice concentration, it is necessary to which the DN values of data to be changed into bright temperature value (with actual physical meaning
Thing reflectivity), using radiant correction formula L=Gain*DN+Bias, wherein L is the bright temperature value after radiant correction, Gain and
Bias is radiant correction parameter, and DN is pixel value in the level one data that FY-3MWRI is provided, according to the science of FY-3MWRI offers
Data set understands Gain=-0.01, Bias=0;The FY-3MWRI level one datas of acquisition are substituted into radiant correction formula and obtained
Data after to correction;
Data are spliced:The FY-3MWRI data of acquisition are single track imaging, if so wanting to carry out research needs to Arctic
Some width images are spliced, image mosaic are carried out using ENVI softwares, by the longitude and latitude file in FY-3MWRI data sets
It imported into ENVI, the splicing of data is carried out by building the GLT files with coordinate system;
Land mask:Because winter part land area has snow cover, in order to avoid this factor is intensive to sea ice
The identification of degree is impacted, and land mask is handled by the band math in ENVI softwares;
Step 3), two marine sites are chosen, i.e.,:To the north of Canadian archipelago for many years ice formation domain, on the south the sea ice edge of GRENLOND
Region, the selection of the data progress sample point of 1 year is chosen, the polarization difference in two regions is calculated respectively, carries out probability distribution
Statistics, the numerical value that daily maximum probability occurs is chosen, as the polarization difference value on the same day, then takes annual polarization difference value
Average value;Determine the mooring points value of pure water and pure ice;P is designated as respectively0(pure water), P1(pure ice), wherein P=T89V-T89H, T89V
For the vertical polarization of 89GHz frequency ranges, T89H is the horizontal polarization of 89GHz frequency ranges, and P is that the polarization of all kinds of pixels is poor;It is final to determine
P0=47.6K, P1=10.8K.
Step 5), the P by acquisition0=47.6K, P1=10.8K values, with interpolation calculation, determine that ASI ice concentrations are anti-
The parameter of formula is drilled, ASI ice concentration inversion algorithms are specific as follows:
P=Tbv-Tbh (1)
Wherein, TbvIt is the bright temperature of vertical polarization, TbhIt is the bright temperature of horizontal polarization;
Select ice concentration of the three rank multinomials fitting from 0% to 100% as follows:
C=d3p3+d2p2+d1p+d0 (2)
Step 6), the P by acquisition0=47.6K, P1=10.8K values are updated in formula (2), then again to formula (2) derivation,
The polarization difference of ice face significantly less than water polarization difference ice concentration C convergence 0 and 1 when polarization difference P difference P0
And P1, the quaternary once linear equation group for solving formula (2) coefficient is drawn, formula (3) is seen, d is calculated using formula (3)0、
d1、 d2、d3
By d0、d1、d2、d3It is brought into formula (2) and draws ice concentration C.
Finally determination ice concentration C expression formula is:
C=6.45714 × 10-6P3-6.05256×10-4P2-9.22521×10-3P+1.10031 (4)
P is obtained by step 3)0=47.6K, P1=10.8K values, and then it is public to obtain ice concentration provided by the invention
Formula:
C=1.29 × 10-5P3-1.28×10-3P2+1.01×10-2P+1.02 (5)
Ice concentration result based on the line number evidences of FY-3MWRI mono- is finally inversed by by formula (5);
Step 7), often influenceed by weather with high-frequency data it is more serious, so need weather filter to be handled,
Because the threshold value of weather filter can also change with the change of data, so the method for the invention based on image entropy, determines base
In the weather filter threshold value of FY-3MWRI data;Weather filter is drawn by MATLAB majorized function fmin-search
Threshold values;Ice concentration result is finally given by filtering process twice.The weather based on FY-3MWRI data is determined
Wave filter threshold value.The method of wherein image entropy is specific as follows:
Kapur methods are a kind of classical ways based on entropy, and for simple two classes partition problem, this method passes through excellent
Change following criterion function to select corresponding threshold value:
Topt=arg max [HI(T)+HW(T)]
Wherein:HIAnd H (T)W(T) size (entropy) that sea ice pixel and seawater pixel include information content is represented respectively;p(g)
Represent that gray level is the ratio shared by p pixel in image all pixels;P (T) represent accumulated probability, i.e., gray level from 0 to
Total pixel number purpose ratio shared by T pixel;G represents the gray level of maximum possible, i.e., when image has 256 gray levels,
G=255.The threshold value when object function maximizes is selected, as segmentation threshold.In specific method implementation process, adopt
With MATLAB, contained majorized function fmin-search is carried out in itself.
By MATLAB computings, the threshold value for drawing GR (37/19) weather filter is 0.05, i.e., when GR (37/19) is more than
It is 0 equal to 0.05 its seasonal ice concentration;The threshold value of GR (23/19) weather filter is 0.045, i.e., when GR (23/19) is big
In being 0 equal to 0.05 its seasonal ice concentration., will be all by being filtered twice to PRELIMINARY RESULTS obtained in the previous step
The weather sea ice that influences and cause to calculate mistake filter out, and then be corrected as seawater.After eventually passing through weather filter processing
Ice concentration result it is as shown in Figure 1.
Contrast verification:Arctic SSM/I data on January 3rd, 2016 are downloaded by ice and snow data center of the U.S., and adopted
With ASI ice concentration algorithms, ice concentration result (as shown in Figure 2) is finally finally inversed by, is carried out pair with the result of the present invention
Than checking.As a result show that both sea ice edges are consistent, ice concentration distribution is identical.It is provided by the present invention to be based on FY-
The Sea Ice Model closeness inverting research of 3MWRI data, for the first time by domestic satellite FY-3MWRI data, applies to ASI algorithms
In, it is innovative to be determined the core parameter based on ice concentration calculation formula in FY-3MWRI data ASI algorithms, and
The selection of effective threshold value of corresponding weather filter, therefore, research of the invention are to utilize domestic FY- to a certain extent
3MWRI data inversion Sea Ice Model closenesses have established certain basis, and technical guarantee sum is provided for the research in arctic navigation channel
According to support.
Examples detailed above is only the preferred embodiment of the present invention, is not intended to limit the invention, for the technology of this area
For personnel, the present invention can have various modifications and variations, and within the spirit and principles of the invention, that is made any repaiies
Change, equivalent substitution, improvement etc., should be included in the scope of the protection.
Claims (1)
- A kind of 1. method from FY-3MWRI data inversion Sea Ice Model closenesses, it is characterised in that:Comprise the following steps:Step 1), obtain FY-3MWRI level one datas.Step 2), the FY-3MWRI level one datas to acquisition pre-process;It is divided into radiant correction, data splicing, land mask;Radiant correction:Using radiant correction formula L=Gain*DN+Bias, wherein L is the bright temperature value after radiant correction, Gain with And Bias is radiant correction parameter, DN is pixel value in the level one data that FY-3MWRI is provided, according to the section of FY-3MWRI offers Learn data set and understand Gain=-0.01, Bias=0;The FY-3MWRI level one datas of acquisition are substituted into radiant correction formula and obtained Data after being corrected;Data are spliced:The FY-3MWRI data of acquisition are single track imaging, image mosaic are carried out using ENVI softwares, by FY-3MWRI Longitude and latitude file in data set is imported into ENVI, and the splicing of data is carried out by building the GLT files with coordinate system;Land mask:Land mask is handled by the band math in ENVI softwares;Step 3), two marine sites are chosen, choose the selection of the data progress sample point of 1 year, calculate the polarization in two regions respectively Difference, probability distribution statistical is carried out, the numerical value that daily maximum probability occurs is chosen, as the polarization difference value on the same day, Ran Houqu The average value of annual polarization difference value;Determine the mooring points value of pure water and pure ice;P is designated as respectively0(pure water), P1(pure ice), P= T89V-T89H, wherein T89V are the vertical polarization of 89GHz frequency ranges, and T89H is the horizontal polarization of 89GHz frequency ranges, and P is all kinds of pixels Polarization it is poor;Step 4), the P by acquisition0、P1Value, with interpolation calculation, determine the parameter of ASI ice concentration inversion formulas, ASI seas Ice concentration inversion algorithm is specific as follows:P=Tbv-Tbh (1)Wherein, TbvIt is the bright temperature of vertical polarization, TbhIt is the bright temperature of horizontal polarization;Select ice concentration of the three rank multinomials fitting from 0% to 100% as follows:C=d3p3+d2p2+d1p+d0 (2)Step 5), the P by acquisition0、P1Value is updated in formula (2), then significantly small to formula (2) derivation, the polarization difference of ice face again Polarization difference P distinguishes P when the polarization difference of water, ice concentration C convergences 0 and 10And P1, draw for solving formula (2) the quaternary once linear equation group of coefficient, is shown in formula (3), and d is calculated using formula (3)0、d1、d2、d3<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msup> <msub> <mi>P</mi> <mn>0</mn> </msub> <mn>3</mn> </msup> </mrow> </mtd> <mtd> <mrow> <msup> <msub> <mi>P</mi> <mn>0</mn> </msub> <mn>2</mn> </msup> </mrow> </mtd> <mtd> <msub> <mi>P</mi> <mn>0</mn> </msub> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <msub> <mi>P</mi> <mn>1</mn> </msub> <mn>3</mn> </msup> </mrow> </mtd> <mtd> <mrow> <msup> <msub> <mi>P</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> </mrow> </mtd> <mtd> <msub> <mi>P</mi> <mn>1</mn> </msub> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>3</mn> <msup> <msub> <mi>P</mi> <mn>0</mn> </msub> <mn>2</mn> </msup> </mrow> </mtd> <mtd> <msub> <mi>P</mi> <mn>1</mn> </msub> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>3</mn> <msup> <msub> <mi>P</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> </mrow> </mtd> <mtd> <mrow> <mn>2</mn> <msub> <mi>P</mi> <mn>1</mn> </msub> </mrow> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>&CenterDot;</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>d</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>d</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>d</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>d</mi> <mn>0</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <mn>1.14</mn> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <mn>0.14</mn> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>By d0、d1、d2、d3It is brought into formula (2) and draws ice concentration C.Finally determination ice concentration C expression formula is:C=6.45714 × 10-6P3-6.05256×10-4P2-9.22521×10-3P+1.10031 (4)P is obtained by step 4)0、P1Value, and then obtain ice concentration formula provided by the invention:C=1.29 × 10-5P3-1.28×10-3P2+1.01×10-2P+1.02 (5)Ice concentration result based on the line number evidences of FY-3MWRI mono- is finally inversed by by formula (5);Step 6), the method based on image entropy, it is determined that the weather filter threshold value based on FY-3MWRI data;Pass through MATLAB's Majorized function fmin-search draws the threshold values of weather filter;Ice concentration knot is finally given by filtering process twice Fruit.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711101311.7A CN107886473B (en) | 2017-11-09 | 2017-11-09 | Method for inverting north sea ice concentration from FY-3MWRI data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711101311.7A CN107886473B (en) | 2017-11-09 | 2017-11-09 | Method for inverting north sea ice concentration from FY-3MWRI data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107886473A true CN107886473A (en) | 2018-04-06 |
CN107886473B CN107886473B (en) | 2020-12-11 |
Family
ID=61779689
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711101311.7A Expired - Fee Related CN107886473B (en) | 2017-11-09 | 2017-11-09 | Method for inverting north sea ice concentration from FY-3MWRI data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107886473B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109765550A (en) * | 2019-01-17 | 2019-05-17 | 中国人民解放军61741部队 | Sea ice thickness inversion method, system and electronic equipment |
CN111104748A (en) * | 2019-12-23 | 2020-05-05 | 南京大学 | Method for calculating drift velocity of arctic sea ice |
CN111563318A (en) * | 2020-04-15 | 2020-08-21 | 中国科学院国家空间科学中心 | Method and system for inverting sea ice density by using 89GHz single-frequency multi-incidence-angle bright temperature difference |
CN112051221A (en) * | 2020-08-27 | 2020-12-08 | 武汉大学 | Sea ice density obtaining method based on space-time system point value |
CN112818851A (en) * | 2020-09-02 | 2021-05-18 | 河南工业大学 | Method for detecting icebound lake based on FY-3MWRI data |
CN113127794A (en) * | 2021-04-07 | 2021-07-16 | 中山大学 | Method for calculating density of arctic sea ice |
CN115856879A (en) * | 2022-11-30 | 2023-03-28 | 南京信息工程大学 | Sea ice melting period intensity inversion method |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101424741A (en) * | 2008-12-08 | 2009-05-06 | 中国海洋大学 | Real time extracting method for satellite remote sensing sea fog characteristic quantity |
CN102567730A (en) * | 2011-11-25 | 2012-07-11 | 中国海洋大学 | Method for automatically and accurately identifying sea ice edge |
CN102708369A (en) * | 2012-05-14 | 2012-10-03 | 大连理工大学 | Sea ice parameter extraction method on basis of satellite image |
CN102866171A (en) * | 2012-09-14 | 2013-01-09 | 江苏科技大学 | Backward electromagnetic scattering coefficient detection module of snow cover-covered sea ice |
CN103278083A (en) * | 2013-05-08 | 2013-09-04 | 南京信息工程大学 | Global navigation satellite signal reflectometry (GNSS-R) detection equipment for sea ice thickness and method for detecting sea ice thickness by utilizing equipment |
CN104198052A (en) * | 2014-09-25 | 2014-12-10 | 国家卫星海洋应用中心 | Acquisition method for sea ice concentration on basis of ocean No. II satellite scanning microwave radiometer |
CN106197383A (en) * | 2016-06-29 | 2016-12-07 | 南京大学 | A kind of remote sensing estimation method of sea ice volume |
US20170074793A1 (en) * | 2015-09-14 | 2017-03-16 | OptikTechnik LLC | Optical sensing device and method in a liquid treatment system |
-
2017
- 2017-11-09 CN CN201711101311.7A patent/CN107886473B/en not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101424741A (en) * | 2008-12-08 | 2009-05-06 | 中国海洋大学 | Real time extracting method for satellite remote sensing sea fog characteristic quantity |
CN102567730A (en) * | 2011-11-25 | 2012-07-11 | 中国海洋大学 | Method for automatically and accurately identifying sea ice edge |
CN102708369A (en) * | 2012-05-14 | 2012-10-03 | 大连理工大学 | Sea ice parameter extraction method on basis of satellite image |
CN102866171A (en) * | 2012-09-14 | 2013-01-09 | 江苏科技大学 | Backward electromagnetic scattering coefficient detection module of snow cover-covered sea ice |
CN103278083A (en) * | 2013-05-08 | 2013-09-04 | 南京信息工程大学 | Global navigation satellite signal reflectometry (GNSS-R) detection equipment for sea ice thickness and method for detecting sea ice thickness by utilizing equipment |
CN104198052A (en) * | 2014-09-25 | 2014-12-10 | 国家卫星海洋应用中心 | Acquisition method for sea ice concentration on basis of ocean No. II satellite scanning microwave radiometer |
US20170074793A1 (en) * | 2015-09-14 | 2017-03-16 | OptikTechnik LLC | Optical sensing device and method in a liquid treatment system |
CN106197383A (en) * | 2016-06-29 | 2016-12-07 | 南京大学 | A kind of remote sensing estimation method of sea ice volume |
Non-Patent Citations (3)
Title |
---|
TIMO VIHMA: "Effects of Arctic Sea Ice Decline on Weather and Climate: A Review", 《SURV GEOPHYS》 * |
刘志强等: "渤海AVHRR多通道海冰密集度反演算法试验研究", 《海洋学报》 * |
苏洁等: "极区海冰密集度AMSR-E数据反演算法的试验与验证", 《遥感学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109765550A (en) * | 2019-01-17 | 2019-05-17 | 中国人民解放军61741部队 | Sea ice thickness inversion method, system and electronic equipment |
CN111104748A (en) * | 2019-12-23 | 2020-05-05 | 南京大学 | Method for calculating drift velocity of arctic sea ice |
CN111104748B (en) * | 2019-12-23 | 2023-09-05 | 南京大学 | Calculation method of north pole sea ice drift velocity |
CN111563318A (en) * | 2020-04-15 | 2020-08-21 | 中国科学院国家空间科学中心 | Method and system for inverting sea ice density by using 89GHz single-frequency multi-incidence-angle bright temperature difference |
CN111563318B (en) * | 2020-04-15 | 2023-07-21 | 中国科学院国家空间科学中心 | Method and system for inverting sea ice concentration by utilizing 89GHz single-frequency multi-incidence-angle bright temperature difference |
CN112051221A (en) * | 2020-08-27 | 2020-12-08 | 武汉大学 | Sea ice density obtaining method based on space-time system point value |
CN112818851A (en) * | 2020-09-02 | 2021-05-18 | 河南工业大学 | Method for detecting icebound lake based on FY-3MWRI data |
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 |
CN115856879A (en) * | 2022-11-30 | 2023-03-28 | 南京信息工程大学 | Sea ice melting period intensity inversion method |
Also Published As
Publication number | Publication date |
---|---|
CN107886473B (en) | 2020-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107886473A (en) | From the method for FY 3MWRI data inversion Sea Ice Model closenesses | |
Xu et al. | Climate forcing and salinity variability in Chesapeake Bay, USA | |
CN101976429A (en) | Cruise image based imaging method of water-surface aerial view | |
CN104198052A (en) | Acquisition method for sea ice concentration on basis of ocean No. II satellite scanning microwave radiometer | |
Ding et al. | Effect of coastal-trapped waves on the synoptic variations of the Yellow Sea Warm Current during winter | |
CN112051221A (en) | Sea ice density obtaining method based on space-time system point value | |
Holland et al. | Ducted coastal ridging over SE Australia | |
Takagi et al. | Track analysis and storm surge investigation of 2017 Typhoon Hato: were the warning signals issued in Macau and Hong Kong timed appropriately? | |
Zhang et al. | Multi-temporal SAR image classification of coastal plain wetlands using a new feature selection method and random forests | |
Mbululo et al. | Climate characteristics over southern highlands Tanzania | |
Smeed et al. | Observed decline of the Atlantic Meridional Overturning Circulation 2004 to 2012. | |
Lorente et al. | Long-term skill assessment of SeaSonde radar-derived wave parameters in the Galician coast (NW Spain) | |
Dieng et al. | Trains of African easterly waves and their relationship to tropical cyclone genesis in the eastern Atlantic | |
Parra et al. | Hydrographic conditions during two austral summer situations (2015 and 2017) in the Gerlache and Bismarck straits, northern Antarctic Peninsula | |
CN111340668B (en) | Typhoon disaster assessment system | |
Liu et al. | A comparison of multiple approaches to study the modulation of ocean waves due to climate variability | |
CN107632328A (en) | A kind of coastal ocean weather information analysis and processing method | |
CN115586591A (en) | Automatic storm surge forecasting method and device | |
Lü et al. | Abnormal reverse intrusion of the Kuroshio Branch Current induced by super typhoon soudelor | |
Gartsman et al. | Analysis of the structure of river systems and the prospects for modeling hydrological processes | |
CN104408750A (en) | Method and system for calculating proportion of seawater to land | |
Zhou et al. | Cross‐Shelf Penetrating Fronts of Buoyant Coastal Currents Around the Headland | |
Xu | Automatic sampling of seawater quality based on electric propulsion unmanned ship | |
Berthet et al. | Subseasonal variability of mesoscale convective systems over the tropical northeastern Pacific | |
Strukov et al. | Relaxation Method for Assimilation of Sea Ice Concentration Data in the NEMO—LIM3 Multicategory Sea Ice Model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for 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: 20201211 Termination date: 20211109 |
|
CF01 | Termination of patent right due to non-payment of annual fee |