CN113739788A - Geographical position correction method and device for brightness temperature data - Google Patents
Geographical position correction method and device for brightness temperature data Download PDFInfo
- Publication number
- CN113739788A CN113739788A CN202110884885.6A CN202110884885A CN113739788A CN 113739788 A CN113739788 A CN 113739788A CN 202110884885 A CN202110884885 A CN 202110884885A CN 113739788 A CN113739788 A CN 113739788A
- Authority
- CN
- China
- Prior art keywords
- coastline
- point
- brightness temperature
- points
- brightness
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000008859 change Effects 0.000 claims abstract description 37
- 238000012545 processing Methods 0.000 claims description 20
- 238000012888 cubic function Methods 0.000 claims description 12
- 230000003247 decreasing effect Effects 0.000 claims description 6
- 238000009795 derivation Methods 0.000 claims description 6
- 238000012163 sequencing technique Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 abstract description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 230000001737 promoting effect Effects 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 241001270131 Agaricus moelleri Species 0.000 description 1
- 238000002679 ablation Methods 0.000 description 1
- 238000012885 constant function Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Image Processing (AREA)
Abstract
本发明属于地球科学技术领域，公开了一种亮温数据的地理位置校正方法及装置，所述方法包括：获取亮温数据，根据亮温数据的亮温梯度变化得到海岸线点集；将海岸线点集中的第一海岸线点与第二海岸线点进行对比，得到第一海岸线点的地理定位误差集；统计所述地理定位误差集中的地理定位误差在沿轨和交轨方向上的分布特征；根据所述分布特征对亮温数据进行地理位置校正。有益效果：本发明提出以第一海岸线点为地面控制点的图像匹配法，在扫描坐标系中对观测角误差造成的地理定位误差进行订正，力图避免不同坐标系之间的转换，以较小的计算量和较简单的方法，实现卫星亮温数据的地理位置校正，为提高和促进国产卫星数据的应用提供有效的支持和保障。
The invention belongs to the field of earth science and technology, and discloses a method and device for correcting the geographic location of brightness temperature data. The method includes: acquiring brightness temperature data, obtaining a coastline point set according to the brightness temperature gradient change of the brightness temperature data; The concentrated first coastline points are compared with the second coastline points to obtain the geolocation error set of the first coastline point; the distribution characteristics of the geolocation errors in the geolocation error set in the along-track and rail-crossing directions are counted; The above distribution characteristics are used to correct the geographic location of the brightness temperature data. Beneficial effects: The present invention proposes an image matching method that takes the first coastline point as the ground control point, corrects the geographic positioning error caused by the observation angle error in the scanning coordinate system, and tries to avoid the conversion between different coordinate systems, so as to reduce the It can realize the geographic location correction of satellite brightness temperature data with a relatively simple method, and provide effective support and guarantee for improving and promoting the application of domestic satellite data.
Description
Technical Field
The invention relates to the technical field of geoscience, in particular to a geographical position correction method and a geographical position correction device for brightness and temperature data.
Background
The satellite-borne passive microwave data has the capability of carrying out earth observation all day long and all day long because the satellite-borne passive microwave data is not influenced by atmospheric conditions, and is widely applied to monitoring of variables such as sea ice density, ice cover/glacier ablation, soil humidity and the like. But the existing passive microwave data generally has geographical position deviation. If AMSR-E has a geometric geographic deviation of about 3.5 to 7km, FY-3B MWRI 89GHz data has an along-track position deviation of 3-4 km and an cross-track position deviation less than 1km, FY-3C MWRI 89GHz data has a geographic position deviation of 5-6 km, and early data of SSM/I has a geographic positioning error of up to 25 km. Accurate geographic positioning is exactly one of the basic requirements of quantitative inversion of ground and atmospheric parameters of passive microwave sensors, so that the correction of the geographic position of data is particularly important.
The brightness temperature data is one of satellite-borne passive microwave data. At present, satellite-borne passive microwave data correction methods comprise three types, namely an image correlation method, a lifting rail brightness temperature difference method and a coastline inflection point method. These methods mostly correct geometric geographical deviations by estimating the deviations of the satellite orbit and attitude parameters. The method has high correction precision, but complex conversion from an antenna coordinate system to a geographic coordinate is required in the correction process, and the calculation amount is huge. In recent years, with the continuous development of sensor technology and information acquisition technology, the remote sensing earth observation load with high resolution and high dynamic has the advantages of large number of wave bands, high spectral and spatial resolution, high data speed, short period and particularly large data volume. The method has the advantages that higher requirements are provided for satellite service processing capacity, the processing speed of remote sensing data is increased on the premise of ensuring product precision, and the problem of real-time or near-real-time remote sensing application requirements is very important.
Therefore, a new method and a device for correcting the geographical position of the brightness temperature data are needed, so that the geographical position of the brightness temperature data can be efficiently and accurately corrected, and the application requirements of the data can be met.
Disclosure of Invention
The purpose of the invention is: the geographical position correction method and device for the brightness temperature data are provided, so that the geographical position of the brightness temperature data can be efficiently and accurately corrected, and the application requirements of the data are met.
In order to achieve the above object, the present invention provides a geographical position correction method for brightness and temperature data, including:
the method comprises the steps of obtaining brightness temperature data of a first area including a sea and land boundary, and obtaining a coastline point set according to brightness temperature gradient changes of the brightness temperature data.
Comparing each first coastline point in the set of coastline points with a corresponding second coastline point in the set of real coastline points to obtain a set of geolocation errors including the geolocation error of each first coastline point.
And counting the distribution characteristics of each geographical positioning error in the geographical positioning error set in the along-track and cross-track directions.
And correcting the geographical position of the brightness temperature data according to the distribution characteristics.
Further, the obtaining of the coastline point set according to the brightness temperature gradient change of the brightness temperature data specifically includes:
acquiring bright temperature data of a first area in a scanning plane coordinate system, and marking a second coastline point according to a real coastline point set in the bright temperature data; the brightness temperature data is brightness temperature values of X rows and Y columns;
searching a plurality of brightness temperature points with the maximum brightness temperature gradient change in the X-axis direction or the Y-axis direction of each second coastline point; the brightness temperature gradient is the speed of the brightness temperature value change;
establishing a corresponding coordinate system according to the ordinate or abscissa of the plurality of bright temperature points corresponding to each second coastline point in the scanning plane coordinate system and the bright temperature values of the plurality of bright temperature points, and respectively performing function fitting;
and respectively carrying out second-order derivation on a plurality of first formulas obtained by function fitting, taking a point with a second-order derivative of zero as a first coastline point, and summarizing all the first coastline points into a coastline point set.
Further, the searching for a plurality of brightness temperature points with the largest brightness temperature gradient change specifically includes:
acquiring four brightness temperature points with the largest brightness temperature gradient change, and sequencing the four brightness temperature points according to the size of the abscissa or the ordinate of each brightness temperature point on a scanning plane;
when the extracted lighting temperatures of the four lighting temperature points are monotonically increased or decreased and the lighting temperature difference value between the first point and the fourth point of the four lighting temperature points is larger than the first threshold value, fitting a unitary cubic function to the coordinates of the four lighting temperature points in the scanning plane and taking the inflection point of the unitary cubic function as the coastline point.
Further, comparing each first coastline point in the set of coastline points with a corresponding second coastline point in the set of real coastline points to obtain a set of geolocation errors including the geolocation error of each first coastline point, specifically:
and each second coastline point corresponds to a first coastline point fitted by a function, the distance between each second coastline point and the corresponding first coastline point is obtained as a geographic positioning error, and all the geographic positioning errors are summarized into a geographic positioning error set.
Further, the counting of the distribution characteristics of each geographic positioning error in the geographic positioning error set in the along-rail and cross-rail directions specifically includes:
decomposing each geographical positioning error in the geographical positioning error set in the along track direction and the cross track direction, respectively performing function fitting on the geographical positioning errors decomposed in the along track direction and the cross track direction, and taking a second formula obtained by the function fitting as the distribution characteristics of the geographical positioning errors in the along track direction and the cross track direction.
Further, the geographic position correction of the brightness temperature data according to the distribution characteristics specifically includes:
substituting the coordinates of the first area brightness temperature data in the scanning plane coordinate system into a second formula to obtain a second geographical positioning error of the first area in the scanning plane coordinate system; and converting the scanning plane coordinate into a longitude and latitude coordinate according to a difference formula to obtain corrected brightness temperature data.
The invention also discloses a geographical position correction device of the brightness temperature data, which comprises the following components: the device comprises a first acquisition module, a first processing module, a second processing module and a third processing module.
The first obtaining module is used for obtaining brightness temperature data of a first area including a sea and land boundary, and obtaining a coastline point set according to brightness temperature gradient change of the brightness temperature data.
The first processing module is configured to compare each first coastline point in the set of coastline points with a corresponding second coastline point in the set of real coastline points to obtain a set of geographic positioning errors including a geographic positioning error of each first coastline point.
And the second processing module is used for counting the distribution characteristics of each geographical positioning error in the geographical positioning error set in the along-track and cross-track directions.
And the third processing module is used for correcting the geographical position of the brightness temperature data according to the distribution characteristics.
Further, the obtaining of the coastline point set according to the brightness temperature gradient change of the brightness temperature data specifically includes:
acquiring bright temperature data of a first area in a scanning plane coordinate system, and marking a second coastline point according to a real coastline point set in the bright temperature data; and the brightness temperature data is the brightness temperature values of X rows and Y columns.
Searching a plurality of brightness temperature points with the maximum brightness temperature gradient change in the X-axis direction or the Y-axis direction of each second coastline point; the brightness temperature gradient is the speed of the brightness temperature value change.
And establishing a corresponding coordinate system according to the ordinate or abscissa of the plurality of bright temperature points corresponding to each second coastline point in the scanning plane coordinate system and the bright temperature values of the plurality of bright temperature points, and respectively performing function fitting.
And respectively carrying out second-order derivation on a plurality of first formulas obtained by function fitting, taking a point with a second-order derivative of zero as a first coastline point, and summarizing all the first coastline points into a coastline point set.
Further, the searching for a plurality of brightness temperature points with the largest brightness temperature gradient change specifically includes:
and acquiring four brightness temperature points with the largest brightness temperature gradient change, and sequencing the four brightness temperature points according to the size of the abscissa or the ordinate of each brightness temperature point on the scanning plane.
When the extracted lighting temperatures of the four lighting temperature points are monotonically increased or decreased and the lighting temperature difference value between the first point and the fourth point of the four lighting temperature points is larger than the first threshold value, fitting a unitary cubic function to the coordinates of the four lighting temperature points in the scanning plane and taking the inflection point of the unitary cubic function as the coastline point.
Further, comparing each first coastline point in the set of coastline points with a corresponding second coastline point in the set of real coastline points to obtain a set of geolocation errors including the geolocation error of each first coastline point, specifically:
and each second coastline point corresponds to a first coastline point fitted by a function, the distance between each second coastline point and the corresponding first coastline point is obtained as a geographic positioning error, and all the geographic positioning errors are summarized into a geographic positioning error set.
Compared with the prior art, the geographic position correction method and the geographic position correction device for the brightness temperature data have the beneficial effects that: the invention provides an image matching method taking a first coastline point as a ground control point, corrects a geographical positioning error caused by an observation angle error in a scanning coordinate system, tries to avoid conversion among different coordinate systems, realizes geographical position correction of satellite brightness temperature data by a smaller calculation amount and a simpler method, and provides effective support and guarantee for improving and promoting application of domestic satellite data.
Drawings
FIG. 1 is a schematic flow chart of a geographical location correction method for brightness and temperature data according to the present invention;
FIG. 2 is a schematic diagram of a geographic position calibration apparatus for brightness and temperature data according to the present invention;
fig. 3 is a schematic reference diagram of a first coastline point acquired in the geographical location correction method of brightness and temperature data according to the present invention.
Detailed Description
The following detailed description of embodiments of the present invention is provided in connection with the accompanying drawings and examples. The following examples are intended to illustrate the invention but are not intended to limit the scope of the invention.
Name interpretation:
FY-3D (FengYun-3, FY-3) is the fourth satellite of the second generation polar orbit meteorological series satellite in China, and is loaded with 10 remote sensing instruments, thereby providing abundant information for land surface and atmosphere monitoring. The MWRI carried by FY-3D is an important radio sensor to monitor the earth in a cone scan with zenith angle of 53.1 ° and width of 1400 km. There are a total of 5 frequencies of 10.65GHz, 18.7GHz, 23.8GHz, 36.5GHz and 89 GHz. At 89GHz channel, the ground resolution of instantaneous field of view is 9km × 15 km.
MWRI (microwave Radiation imager) microwave imager.
Example 1:
the invention discloses a geographical position correction method of brightness temperature data, which is applied to geographical position correction of satellite brightness temperature data and mainly comprises the following steps:
and step S1, acquiring brightness temperature data of a first area including a boundary between the sea and the land, and obtaining a coastline point set according to brightness temperature gradient change of the brightness temperature data.
Step S2, comparing each first coastline point in the set of coastline points with a corresponding second coastline point in the set of real coastline points to obtain a set of geolocation errors including the geolocation error of each first coastline point.
And step S3, counting the distribution characteristics of each geographical positioning error in the geographical positioning error set in the along-track and cross-track directions.
And step S4, carrying out geographical position correction on the brightness temperature data according to the distribution characteristics.
In this embodiment, the brightness temperature data includes brightness temperature values of land and sea, and since the brightness temperature of land is higher than that of sea in the daytime, the brightness temperature gradient changes greatly in the brightness temperature transition region near the coastline. On this basis a first coastline point may be extracted.
In this embodiment, the brightness temperature data includes brightness temperature values in X rows and Y columns, and a corresponding coordinate system may be established according to rows and columns of the scanning plane coordinate system, and a corresponding coordinate system may be established by a row number and a column brightness temperature value or by a column number and a row brightness temperature value. Each light temperature value has a uniquely determined rank number. And each row number and one column number correspond to one longitude and latitude, and the longitude and latitude corresponding to each row number or column number are generated together when the brightness temperature data is generated when the satellite scans the earth surface.
In this embodiment, the method is preferably applied to the wind and cloud series satellites, and can effectively and accurately process and correct the brightness and temperature data of the wind and cloud series satellites.
In step S1, the obtaining a coastline point set according to the brightness temperature gradient change of the brightness temperature data specifically includes:
acquiring bright temperature data of a first area in a scanning plane coordinate system, and marking a second coastline point according to a real coastline point set in the bright temperature data; and the brightness temperature data is the brightness temperature values of X rows and Y columns.
Searching a plurality of brightness temperature points with the maximum brightness temperature gradient change in the X-axis direction or the Y-axis direction of each second coastline point; the brightness temperature gradient is the speed of the brightness temperature value change.
And establishing a corresponding coordinate system according to the ordinate or abscissa of the plurality of bright temperature points corresponding to each second coastline point in the scanning plane coordinate system and the bright temperature values of the plurality of bright temperature points, and respectively performing function fitting.
And respectively carrying out second-order derivation on a plurality of first formulas obtained by function fitting, taking a point with a second-order derivative of zero as a first coastline point, and summarizing all the first coastline points into a coastline point set.
In this embodiment, the searching for the plurality of brightness temperature points with the largest brightness temperature gradient change specifically includes:
acquiring four brightness temperature points with the largest brightness temperature gradient change, and sequencing the four brightness temperature points according to the size of the abscissa or the ordinate of each brightness temperature point on a scanning plane;
when the extracted lighting temperatures of the four lighting temperature points are monotonically increased or decreased and the lighting temperature difference value between the first point and the fourth point of the four lighting temperature points is larger than the first threshold value, fitting a unitary cubic function to the coordinates of the four lighting temperature points in the scanning plane and taking the inflection point of the unitary cubic function as the coastline point.
Referring to fig. 3(a), different gray scales represent different light temperature values, and the real coastline is represented by a specific gray scale. In the figure, the gray level of the rightmost black point is the area of the real coastline, and the center point of the gray level block is the real coastline point. Since the coastline in the drawing is in the column direction, four light temperature points having the largest change in luminance temperature gradient are searched in the row (X) direction of the real coastline point, such as four black points in fig. 3(a), the third point from left to right in the image is not the desired light temperature point, and the rightmost black point is the real coastline point (the second coastline point). Fig. 3(b) plots the light temperature of the four light temperature point positions and the coordinates in the scan plane coordinate system. These four bright temperature points are fitted with polynomial functions as follows:
wherein y is a brightness temperature value. x denotes the position in the scanning plane coordinate system. a, b and c are all constants. By calculating the second derivative of this function, the position in fig. 3(a) and (b) where the gradient change is the largest is obtained, and the third from the left in fig. 3(a) and (b) is the inflection point, which is the coastline inflection point and is also the first coastline point.
In step S2, the step of comparing each first coastline point in the set of coastline points with a corresponding second coastline point in the set of real coastline points to obtain a set of geographic positioning errors including a geographic positioning error of each first coastline point is specifically as follows:
and each second coastline point corresponds to a first coastline point fitted by a function, the distance between each second coastline point and the corresponding first coastline point is obtained as a geographic positioning error, and all the geographic positioning errors are summarized into a geographic positioning error set.
In step S3, the statistics of the distribution characteristics of each geolocation error in the set of geolocation errors in the along-track and cross-track directions specifically includes:
decomposing each geographical positioning error in the geographical positioning error set in the along track direction and the cross track direction, respectively performing function fitting on the geographical positioning errors decomposed in the along track direction and the cross track direction, and taking a second formula obtained by the function fitting as the distribution characteristics of the geographical positioning errors in the along track direction and the cross track direction.
In the present embodiment, for the geolocation error, fitting is performed with a constant function, a linear function, a polynomial function, and the like, in the two directions, respectively. Where y represents the geographical position error. x denotes the position in the scanning plane coordinate system. a, b, c are constants.
Respectively substituting the geographical positioning errors into the formulas to calculate the goodness of fit R2，R2The calculation formula of (a) is as follows:
wherein, yiFor the sample values, the values of the samples,
in order to predict the value of the target,
is the average of the sample values. R2The larger the fit, the better. According to R2And selecting a fitting formula with the best fitting effect from the three formulas to describe the distribution characteristics of the geographic positioning errors.
In step S4, the geographic location correction of the brightness temperature data according to the distribution characteristics specifically includes:
substituting the coordinates of the first area brightness temperature data in the scanning plane coordinate system into a second formula to obtain a second geographical positioning error of the first area in the scanning plane coordinate system; and converting the scanning plane coordinate into a longitude and latitude coordinate according to a difference formula to obtain corrected brightness temperature data.
Step S4 is to recalculate the geographic positioning error according to the distribution characteristics of the geographic positioning error, and correct the geographic positioning error. The step has the effects of reducing the geographical positioning error and improving the geographical positioning accuracy of satellite data. And substituting the coordinates of the MWRI data in the scanning plane coordinate system into the equation in the step S3 to obtain the geographical positioning error of the whole image in the scanning plane coordinate system. And then the scanning plane coordinates are converted into longitude and latitude by the following interpolation formula, and the corrected geographic position data can be obtained.
(lon′,lat′)＝f(x,y,x′,y′,lon,lat)
Where lon ', lat', x ', y' respectively represent longitude and latitude after correction, and x, y coordinates in the scanning plane coordinate system, and x, y, lon, lat respectively represent longitude and latitude before correction, and x, y coordinates in the scanning plane coordinate system. The value of the unknown function f at point P ═ x, y is solved as follows:
wherein f is in Q11＝(x1,y1)，Q12＝(x1,y2)，Q21＝(x2,y1)， Q22＝(x2,y2) The values of the four points are known quantities.
Example 2:
the invention also discloses a geographical position correction device of the brightness temperature data, which comprises the following components: a first acquisition module 101, a first processing module 102, a second processing module 103 and a third processing module 104.
The first obtaining module 101 is configured to obtain brightness temperature data of a first area including a boundary between an ocean and a land, and obtain a coastline point set according to brightness temperature gradient change of the brightness temperature data;
the first processing module 102 is configured to compare each first coastline point in the set of coastline points with a corresponding second coastline point in the set of real coastline points to obtain a set of geographic positioning errors including a geographic positioning error of each first coastline point;
the second processing module 103 is configured to count distribution characteristics of each geographic positioning error in the geographic positioning error set in along-track and cross-track directions;
the third processing module 104 is configured to perform geographic position correction on the brightness and temperature data according to the distribution characteristics.
In this embodiment, the obtaining of the coastline point set according to the brightness temperature gradient change of the brightness temperature data specifically includes:
acquiring bright temperature data of a first area in a scanning plane coordinate system, and marking a second coastline point according to a real coastline point set in the bright temperature data; and the brightness temperature data is the brightness temperature values of X rows and Y columns.
Searching a plurality of brightness temperature points with the maximum brightness temperature gradient change in the X-axis direction or the Y-axis direction of each second coastline point; the brightness temperature gradient is the speed of the brightness temperature value change.
And establishing a corresponding coordinate system according to the ordinate or abscissa of the plurality of bright temperature points corresponding to each second coastline point in the scanning plane coordinate system and the bright temperature values of the plurality of bright temperature points, and respectively performing function fitting.
And respectively carrying out second-order derivation on a plurality of first formulas obtained by function fitting, taking a point with a second-order derivative of zero as a first coastline point, and summarizing all the first coastline points into a coastline point set.
In this embodiment, the searching for the plurality of brightness temperature points with the largest brightness temperature gradient change specifically includes:
acquiring four brightness temperature points with the largest brightness temperature gradient change, and sequencing the four brightness temperature points according to the size of the abscissa or the ordinate of each brightness temperature point on a scanning plane;
when the extracted lighting temperatures of the four lighting temperature points are monotonically increased or decreased and the lighting temperature difference value between the first point and the fourth point of the four lighting temperature points is larger than the first threshold value, fitting a unitary cubic function to the coordinates of the four lighting temperature points in the scanning plane and taking the inflection point of the unitary cubic function as the coastline point.
In this embodiment, the comparing each first coastline point in the set of coastline points with a corresponding second coastline point in the set of real coastline points to obtain a set of geographic positioning errors including a geographic positioning error of each first coastline point specifically includes:
and each second coastline point corresponds to a first coastline point fitted by a function, the distance between each second coastline point and the corresponding first coastline point is obtained as a geographic positioning error, and all the geographic positioning errors are summarized into a geographic positioning error set.
In this embodiment, embodiment 2 is written based on embodiment 1, and therefore some of the same technical features are not described in detail.
To sum up, compared with the prior art, the geographic position correction method and device for brightness and temperature data provided by the embodiments of the present invention have the following beneficial effects: the invention provides an image matching method taking a first coastline point as a ground control point, corrects a geographical positioning error caused by an observation angle error in a scanning coordinate system, tries to avoid conversion among different coordinate systems, realizes geographical position correction of satellite brightness temperature data by a smaller calculation amount and a simpler method, and provides effective support and guarantee for improving and promoting application of domestic satellite data.
The above description is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, various modifications and substitutions can be made without departing from the technical principle of the present invention, and these modifications and substitutions should also be regarded as the protection scope of the present invention.
Claims (10)
1. A geographical position correction method of brightness temperature data is characterized by comprising the following steps:
obtaining brightness temperature data of a first area including a boundary between an ocean and a land, and obtaining a coastline point set according to brightness temperature gradient change of the brightness temperature data;
comparing each first coastline point in the set of coastline points with a corresponding second coastline point in the set of real coastline points to obtain a set of geographic positioning errors including the geographic positioning error of each first coastline point;
counting distribution characteristics of each geographical positioning error in the geographical positioning error set in the along-track and cross-track directions;
and correcting the geographical position of the brightness temperature data according to the distribution characteristics.
2. The geographical location correction method of brightness temperature data according to claim 1, wherein the obtaining of the coastline point set according to the brightness temperature gradient change of the brightness temperature data specifically comprises:
acquiring bright temperature data of a first area in a scanning plane coordinate system, and marking a second coastline point according to a real coastline point set in the bright temperature data; the brightness temperature data is brightness temperature values of X rows and Y columns;
searching a plurality of brightness temperature points with the maximum brightness temperature gradient change in the X-axis direction or the Y-axis direction of each second coastline point; the brightness temperature gradient is the speed of the brightness temperature value change;
establishing a corresponding coordinate system according to the ordinate or abscissa of the plurality of bright temperature points corresponding to each second coastline point in the scanning plane coordinate system and the bright temperature values of the plurality of bright temperature points, and respectively performing function fitting;
and respectively carrying out second-order derivation on a plurality of first formulas obtained by function fitting, taking a point with a second-order derivative of zero as a first coastline point, and summarizing all the first coastline points into a coastline point set.
3. The geographical location correction method of brightness temperature data according to claim 2, wherein the searching for a plurality of brightness temperature points with the largest brightness temperature gradient change specifically comprises:
acquiring four brightness temperature points with the largest brightness temperature gradient change, and sequencing the four brightness temperature points according to the size of the abscissa or the ordinate of each brightness temperature point on a scanning plane;
when the extracted lighting temperatures of the four lighting temperature points are monotonically increased or decreased and the lighting temperature difference value between the first point and the fourth point of the four lighting temperature points is larger than the first threshold value, fitting a unitary cubic function to the coordinates of the four lighting temperature points in the scanning plane and taking the inflection point of the unitary cubic function as the coastline point.
4. The method according to claim 2, wherein the step of comparing each first coastline point in the set of coastline points with a corresponding second coastline point in the set of real coastline points to obtain a set of geolocation errors including the geolocation error of each first coastline point comprises:
and each second coastline point corresponds to a first coastline point fitted by a function, the distance between each second coastline point and the corresponding first coastline point is obtained as a geographic positioning error, and all the geographic positioning errors are summarized into a geographic positioning error set.
5. The geographical location correction method of brightness temperature data according to claim 4, wherein the statistical distribution characteristics of each geographical positioning error in the geographical positioning error set in the along-track and cross-track directions are specifically:
decomposing each geographical positioning error in the geographical positioning error set in the along track direction and the cross track direction, respectively performing function fitting on the geographical positioning errors decomposed in the along track direction and the cross track direction, and taking a second formula obtained by the function fitting as the distribution characteristics of the geographical positioning errors in the along track direction and the cross track direction.
6. The geographical location correction method of brightness and temperature data according to claim 5, wherein the geographical location correction of the brightness and temperature data according to the distribution characteristics specifically includes:
substituting the coordinates of the first area brightness temperature data in the scanning plane coordinate system into a second formula to obtain a second geographical positioning error of the first area in the scanning plane coordinate system; and converting the scanning plane coordinate into a longitude and latitude coordinate according to a difference formula to obtain corrected brightness temperature data.
7. A geographical position correction device for brightness-temperature data, comprising: the device comprises a first acquisition module, a first processing module, a second processing module and a third processing module;
the first acquisition module is used for acquiring brightness temperature data of a first area including a sea and land boundary, and obtaining a coastline point set according to brightness temperature gradient change of the brightness temperature data;
the first processing module is configured to compare each first coastline point in the set of coastline points with a corresponding second coastline point in the set of real coastline points to obtain a set of geographic positioning errors including a geographic positioning error of each first coastline point;
the second processing module is used for counting distribution characteristics of each geographic positioning error in the geographic positioning error set in the along-track and cross-track directions;
and the third processing module is used for correcting the geographical position of the brightness temperature data according to the distribution characteristics.
8. The geographical location correction device of brightness temperature data according to claim 7, wherein the set of coastline points is obtained according to brightness temperature gradient changes of the brightness temperature data, and specifically comprises:
acquiring bright temperature data of a first area in a scanning plane coordinate system, and marking a second coastline point according to a real coastline point set in the bright temperature data; the brightness temperature data is brightness temperature values of X rows and Y columns;
searching a plurality of brightness temperature points with the maximum brightness temperature gradient change in the X-axis direction or the Y-axis direction of each second coastline point; the brightness temperature gradient is the speed of the brightness temperature value change;
establishing a corresponding coordinate system according to the ordinate or abscissa of the plurality of bright temperature points corresponding to each second coastline point in the scanning plane coordinate system and the bright temperature values of the plurality of bright temperature points, and respectively performing function fitting;
and respectively carrying out second-order derivation on a plurality of first formulas obtained by function fitting, taking a point with a second-order derivative of zero as a first coastline point, and summarizing all the first coastline points into a coastline point set.
9. The geographic position correction device of claim 8, wherein the searching for the plurality of brightness temperature points with the largest brightness temperature gradient change includes:
acquiring four brightness temperature points with the largest brightness temperature gradient change, and sequencing the four brightness temperature points according to the size of the abscissa or the ordinate of each brightness temperature point on a scanning plane;
when the extracted lighting temperatures of the four lighting temperature points are monotonically increased or decreased and the lighting temperature difference value between the first point and the fourth point of the four lighting temperature points is larger than the first threshold value, fitting a unitary cubic function to the coordinates of the four lighting temperature points in the scanning plane and taking the inflection point of the unitary cubic function as the coastline point.
10. The geographical location correction device of claim 8, wherein the comparing of each first coastline point in the set of coastline points with a corresponding second coastline point in the set of real coastline points results in a set of geographical positioning errors comprising a geographical positioning error of each first coastline point, specifically:
and each second coastline point corresponds to a first coastline point fitted by a function, the distance between each second coastline point and the corresponding first coastline point is obtained as a geographic positioning error, and all the geographic positioning errors are summarized into a geographic positioning error set.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110884885.6A CN113739788B (en) | 2021-08-03 | 2021-08-03 | Geographic position correction method and device for bright temperature data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110884885.6A CN113739788B (en) | 2021-08-03 | 2021-08-03 | Geographic position correction method and device for bright temperature data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113739788A true CN113739788A (en) | 2021-12-03 |
| CN113739788B CN113739788B (en) | 2023-05-23 |
Family
ID=78729865
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110884885.6A Active CN113739788B (en) | 2021-08-03 | 2021-08-03 | Geographic position correction method and device for bright temperature data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113739788B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115393438A (en) * | 2022-08-29 | 2022-11-25 | 北京华云星地通科技有限公司 | Satellite remote sensing image geographical positioning error detection method and system |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2000042385A1 (en) * | 1999-01-13 | 2000-07-20 | The Secretary Of State For Defence | Improved method of processing data and data processing apparatus |
| JP2005346665A (en) * | 2004-06-07 | 2005-12-15 | Nogiwa Sangyo Kk | Shoreline extraction method and shoreline-extracting system |
| JP2005346664A (en) * | 2004-06-07 | 2005-12-15 | Nogiwa Sangyo Kk | Coastline extraction method and coastline extraction system |
| US20080114546A1 (en) * | 2006-11-15 | 2008-05-15 | Space Systems/Loral, Inc. | Image navigation and registration accuracy improvement using parametric systematic error correction |
| CN104952080A (en) * | 2015-07-10 | 2015-09-30 | 中国科学院软件研究所 | Method for realizing coarse positioning of remote sensing images |
| CN106600589A (en) * | 2016-12-09 | 2017-04-26 | 中国科学院合肥物质科学研究院 | Coastline area remote sensing map-based satellite borne spectrograph pointing registration method |
| WO2017071160A1 (en) * | 2015-10-28 | 2017-05-04 | 深圳大学 | Sea-land segmentation method and system for large-size remote-sensing image |
| CN107870043A (en) * | 2017-10-25 | 2018-04-03 | 中国科学院国家空间科学中心 | A Synchronous Inversion Optimization Method for Sea Surface Parameters |
| CN108550174A (en) * | 2018-03-08 | 2018-09-18 | 中国地质大学（武汉） | A kind of coastline Super-resolution Mapping and system based on half global optimization |
| CN109741446A (en) * | 2018-12-12 | 2019-05-10 | 四川华控图形科技有限公司 | A kind of method of the fine coastal landform of dynamic generation in three-dimensional digital earth |
| CN109886988A (en) * | 2019-01-31 | 2019-06-14 | 湖北大学 | Method, system, device and medium for measuring positioning error of microwave imager |
| CN110991393A (en) * | 2019-12-17 | 2020-04-10 | 北京航天泰坦科技股份有限公司 | Method and device for remote sensing monitoring and analysis of coastline transition |
| CN112051221A (en) * | 2020-08-27 | 2020-12-08 | 武汉大学 | A method for obtaining sea ice density based on the point value of space-time system |
-
2021
- 2021-08-03 CN CN202110884885.6A patent/CN113739788B/en active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2000042385A1 (en) * | 1999-01-13 | 2000-07-20 | The Secretary Of State For Defence | Improved method of processing data and data processing apparatus |
| JP2005346665A (en) * | 2004-06-07 | 2005-12-15 | Nogiwa Sangyo Kk | Shoreline extraction method and shoreline-extracting system |
| JP2005346664A (en) * | 2004-06-07 | 2005-12-15 | Nogiwa Sangyo Kk | Coastline extraction method and coastline extraction system |
| US20080114546A1 (en) * | 2006-11-15 | 2008-05-15 | Space Systems/Loral, Inc. | Image navigation and registration accuracy improvement using parametric systematic error correction |
| CN104952080A (en) * | 2015-07-10 | 2015-09-30 | 中国科学院软件研究所 | Method for realizing coarse positioning of remote sensing images |
| WO2017071160A1 (en) * | 2015-10-28 | 2017-05-04 | 深圳大学 | Sea-land segmentation method and system for large-size remote-sensing image |
| CN106600589A (en) * | 2016-12-09 | 2017-04-26 | 中国科学院合肥物质科学研究院 | Coastline area remote sensing map-based satellite borne spectrograph pointing registration method |
| CN107870043A (en) * | 2017-10-25 | 2018-04-03 | 中国科学院国家空间科学中心 | A Synchronous Inversion Optimization Method for Sea Surface Parameters |
| CN108550174A (en) * | 2018-03-08 | 2018-09-18 | 中国地质大学（武汉） | A kind of coastline Super-resolution Mapping and system based on half global optimization |
| CN109741446A (en) * | 2018-12-12 | 2019-05-10 | 四川华控图形科技有限公司 | A kind of method of the fine coastal landform of dynamic generation in three-dimensional digital earth |
| CN109886988A (en) * | 2019-01-31 | 2019-06-14 | 湖北大学 | Method, system, device and medium for measuring positioning error of microwave imager |
| CN110991393A (en) * | 2019-12-17 | 2020-04-10 | 北京航天泰坦科技股份有限公司 | Method and device for remote sensing monitoring and analysis of coastline transition |
| CN112051221A (en) * | 2020-08-27 | 2020-12-08 | 武汉大学 | A method for obtaining sea ice density based on the point value of space-time system |
Non-Patent Citations (2)
| Title |
|---|
| 景振华;胡秀清;殷德奎;: "TG-2多角度偏振成像仪遥感影像地理定位与误差订正", 遥感技术与应用 * |
| 朱宁龙;: "基于海岸线状特征匹配的近海域卫星遥感影像定位", 海军航空工程学院学报 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115393438A (en) * | 2022-08-29 | 2022-11-25 | 北京华云星地通科技有限公司 | Satellite remote sensing image geographical positioning error detection method and system |
| CN115393438B (en) * | 2022-08-29 | 2023-03-10 | 北京华云星地通科技有限公司 | Method and system for detecting geographical positioning error of remote sensing image of satellite |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113739788B (en) | 2023-05-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wu et al. | Coarse scale in situ albedo observations over heterogeneous snow-free land surfaces and validation strategy: A case of MODIS albedo products preliminary validation over northern China | |
| CN110716185A (en) | An on-orbit radiation correction method for spaceborne microwave radiometers | |
| CN114755661B (en) | A method and device for calibrating parameters of a mobile laser scanning system | |
| CN110174359A (en) | A kind of Airborne Hyperspectral image heavy metal-polluted soil concentration evaluation method returned based on Gaussian process | |
| CN111144350B (en) | Remote sensing image positioning accuracy evaluation method based on reference base map | |
| CN110111274A (en) | A kind of spaceborne push-broom type optical sensor elements of exterior orientation calibrating method | |
| CN115407289B (en) | A geometric calibration method for spaceborne SAR based on control points provided by wind turbines | |
| CN104180781A (en) | Deformation monitoring data processing method for single frequency and double frequency GPS hybrid network | |
| CN116699609A (en) | High-precision airborne SAR image plane positioning method | |
| CN111123345B (en) | GNSS measurement-based empirical ionosphere model data driving method | |
| CN117994678A (en) | Positioning method and system for natural resource remote sensing mapping image | |
| CN113739788A (en) | Geographical position correction method and device for brightness temperature data | |
| CN108983172B (en) | SAR radiation calibration method based on random evolution constraint | |
| Yang et al. | Land contamination correction for passive microwave radiometer data: Demonstration of wind retrieval in the great lakes using SSM/I | |
| CN103777196B (en) | Based on terrain object distance single station measuring method and the measuring system thereof of geography information | |
| CN118606907A (en) | A random forest FY-4A water vapor correction method integrating elevation | |
| CN115688955A (en) | Vegetation growth prediction method, vegetation growth prediction device and storage medium | |
| CN109493298B (en) | A kind of airborne sweep type high-spectral data fast geometric bearing calibration | |
| CN111650128A (en) | A high-resolution atmospheric aerosol retrieval method based on the surface reflectance library | |
| Folger et al. | Lidar-based height correction for the assimilation of atmospheric motion vectors | |
| CN116299466A (en) | Geological deformation monitoring method and device for power transmission channel | |
| CN115830476A (en) | Terrain factor space downscaling method | |
| CN112068157B (en) | Method and device for realizing earth observation mode of stationary orbit multi-frequency terahertz detector | |
| Niu et al. | Development of a new algorithm to identify clear sky MSU data using AMSU-A data for verification | |
| Martínez et al. | Arctic+ Salinity: Algorithm Theoretical Baseline Document |
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 |