CN112612916A - Method and device for generating inspection error spatial distribution map of ocean satellite data - Google Patents
Method and device for generating inspection error spatial distribution map of ocean satellite data Download PDFInfo
- Publication number
- CN112612916A CN112612916A CN202011596962.XA CN202011596962A CN112612916A CN 112612916 A CN112612916 A CN 112612916A CN 202011596962 A CN202011596962 A CN 202011596962A CN 112612916 A CN112612916 A CN 112612916A
- Authority
- CN
- China
- Prior art keywords
- latitude
- longitude
- satellite data
- marine satellite
- data
- 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
- 238000007689 inspection Methods 0.000 title claims abstract description 74
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000012545 processing Methods 0.000 claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims description 21
- 238000004590 computer program Methods 0.000 claims description 13
- 238000003908 quality control method Methods 0.000 claims description 11
- 238000012795 verification Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 abstract description 28
- 238000007796 conventional method Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 12
- 230000006870 function Effects 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/50—Information retrieval; Database structures therefor; File system structures therefor of still image data
- G06F16/58—Retrieval characterised by using metadata, e.g. metadata not derived from the content or metadata generated manually
- G06F16/587—Retrieval characterised by using metadata, e.g. metadata not derived from the content or metadata generated manually using geographical or spatial information, e.g. location
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/50—Information retrieval; Database structures therefor; File system structures therefor of still image data
- G06F16/54—Browsing; Visualisation therefor
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computer Graphics (AREA)
- Remote Sensing (AREA)
- Library & Information Science (AREA)
- Image Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The embodiment of the application provides a method and a device for generating a spatial distribution diagram of inspection errors of ocean satellite data, and relates to the technical field of ocean engineering. The method comprises the steps of carrying out space-time three-dimensional spline interpolation processing on inspection source data based on parameter information of marine satellite data to obtain reference data of the marine satellite data; calculating the observation error of the marine satellite data by using the reference data; projecting the observation error on a preset space equal longitude and latitude grid generated by equal longitude and latitude network division based on marine satellite data to obtain a projection result; calculating a proportionality coefficient between space distances and longitude grids at different latitudes based on the curvature of the earth; and performing near point interpolation by taking the proportional coefficient as an interpolation window according to the projection result to fill uncovered equal longitude and latitude grids in the projection result and generate a marine satellite inspection error spatial distribution map, so that the problem that the spatial distribution of the inspection error cannot be visually and accurately displayed by the conventional method is solved.
Description
Technical Field
The application relates to the technical field of ocean engineering, in particular to a method and a device for generating a spatial distribution diagram of inspection errors of ocean satellite data.
Background
With the rapid development of aerospace technology, remote sensing technology and computer technology, ocean satellite remote sensing is in a new stage of updating and rapid development. The marine satellite has the characteristics of wide-range and high-precision global marine observation, and the application fields of the marine satellite comprise marine ecological environment monitoring, marine dynamic environment forecasting, marine disaster prevention and reduction, marine rights and interests maintenance and the like. The precision of a marine satellite product is the basis of marine satellite remote sensing application, good product quality directly determines the data application effect, but the existing method cannot visually display and accurately display the spatial distribution of inspection errors.
Disclosure of Invention
The embodiment of the application aims to provide a method and a device for generating a spatial distribution diagram of inspection errors of ocean satellite data, which can accurately calculate and visually display the inspection errors in the form of the spatial distribution diagram of the inspection errors, and solve the problem that the spatial distribution of the inspection errors cannot be visually displayed and accurately displayed in the existing method.
The embodiment of the application provides a method for generating a spatial distribution diagram of inspection errors of marine satellite data, which comprises the following steps:
performing space-time three-dimensional spline interpolation processing on inspection source data based on parameter information of the marine satellite data to obtain reference data of the marine satellite data;
calculating an observation error of the marine satellite data by using the reference data;
projecting the observation error on a preset space equal longitude and latitude grid generated by equal longitude and latitude network division based on the marine satellite data to obtain a projection result;
calculating a proportionality coefficient between space distances and longitude grids at different latitudes based on the curvature of the earth;
and performing near point interpolation by taking the proportion coefficient as an interpolation window according to the projection result so as to fill the uncovered equal longitude and latitude grids in the projection result and generate a marine satellite inspection error space distribution map.
In the implementation process, the automatic generation method of the marine satellite data inspection error spatial distribution map is provided, the observation error of the marine satellite data is calculated through space-time three-dimensional spline interpolation, so that the calculation result is more accurate, the positioning of the marine satellite observation error is consistent with that of the marine satellite data, and the display precision of the marine satellite data inspection error spatial distribution map is improved; the method has the advantages that the proportion coefficient between the space distance and the longitude grids in different dimensions is adopted, the proportion coefficient is the same as the dimension of space coverage loss, optimization of interpolation window setting is achieved, the search range of the interpolation window is effectively reduced, the space coverage interpolation efficiency is improved, the automatic generation time of the marine satellite data inspection error space distribution map is shortened, the problem of incomplete coverage of equal longitude and latitude grid projection is effectively solved by adopting a near point interpolation method, the marine satellite inspection error space distribution map with continuous and consistent space is generated, the visualization effect of the marine satellite inspection error space distribution map is effectively improved, the readability of an inspection report is increased, and the problem that the space distribution of inspection errors cannot be visually displayed and accurately displayed by the existing method is solved.
Further, the performing space-time three-dimensional spline interpolation processing on the inspection source data based on the parameter information of the marine satellite data to obtain the reference data of the marine satellite data comprises:
acquiring physical parameters of inspection source data;
and performing space-time three-dimensional spline interpolation processing on the inspection source data based on the physical parameters and the parameter information of the marine satellite data to acquire the reference data.
In the implementation process, the time-space three-dimensional spline interpolation processing is carried out on the inspection source data by adopting a cubic spline interpolation method, and the reference data corresponding to the ocean satellite data is obtained through calculation.
Further, said calculating an observation error of said marine satellite data using said reference data comprises:
acquiring an ocean parameter inversion result corresponding to the ocean satellite data;
acquiring a quality identifier corresponding to the ocean parameter inversion result to perform quality control on the ocean parameter inversion result;
and calculating the observation error of the marine satellite data by using the marine parameter inversion result after the quality control and the reference data.
In the implementation process, the observation error of the marine satellite data is calculated based on the reference data, and the marine satellite observation data needing to be eliminated in the marine parameter inversion result quality identification is set to be an invalid value, so that the accuracy of the observation error calculation result is improved.
Further, the projecting the observation error on a preset space equal longitude and latitude grid generated by equal longitude and latitude network division based on the marine satellite data to obtain a projection result includes:
and performing equal-longitude-latitude grid division according to the spatial resolution of the ocean satellite data along the orbit to generate a spatial equal-longitude-latitude grid, wherein the size of the equal-longitude-latitude grid is expressed as:
Grid_Size＝Spatial_Resolution/100；
wherein, Grid _ Size represents the Size of equal longitude and latitude grids, and the unit is an angle; spatial _ Resolution represents the along-track Spatial Resolution of the marine satellite data, in kilometers;
and calculating the row and column index values of the space equal longitude and latitude grids by utilizing the longitude and latitude of the marine satellite data so as to assign the observation error to the corresponding space equal longitude and latitude grids.
In the implementation process, global (-90 degrees, 0-360 degrees) equal longitude and latitude grid division is carried out according to the spatial resolution of the marine satellite data, and spatial equal longitude and latitude grid projection is carried out on the observation error data.
Further, the computing a scaling factor between spatial distances and longitude grids at different latitudes based on the curvature of the earth includes:
calculating the proportional relation between space distances and longitude grids at different latitudes through the curvature of the earth, wherein the proportional coefficient is expressed as:
wherein, R represents the proportionality coefficient, ceiling represents rounding up, and Latitude represents Latitude.
In the implementation process, the curvature of the earth is approximate to a spherical shape, the proportional relation between the space distance and the longitude grid at different latitudes is calculated, and the proportional coefficient directly calculated by the latitude is rounded to obtain the proportional coefficient of an integer.
Further, the performing near point interpolation by using the proportionality coefficient as an interpolation window according to the projection result to fill the uncovered equal longitude and latitude grids in the projection result and generate a marine satellite inspection error spatial distribution map includes:
circularly judging whether each equal longitude and latitude grid is assigned or not;
if not, searching by taking the proportional coefficients corresponding to the equal longitude and latitude grids as an interpolation window to obtain the nearest effective projection;
and assigning the effective projection to the equal longitude and latitude grids.
In the implementation process, the equal longitude and latitude grids are traversed, the proportional coefficients at different latitudes are used as interpolation windows, the adjacent point interpolation is carried out, and the area with incomplete projection coverage of the equal longitude and latitude grids is filled.
The embodiment of the present application further provides a device for generating a spatial distribution diagram of inspection errors of marine satellite data, where the device includes:
the reference data acquisition module is used for carrying out space-time three-dimensional spline interpolation processing on the inspection source data based on the parameter information of the marine satellite data so as to acquire the reference data of the marine satellite data;
the error calculation module is used for calculating the observation error of the marine satellite data by using the reference data;
the projection module is used for projecting the observation error on a preset space equal longitude and latitude grid generated by equal longitude and latitude network division based on the marine satellite data so as to obtain a projection result;
the scale factor calculation module is used for calculating scale factors between space distances and longitude grids at different latitudes based on the curvature of the earth;
and the interpolation module is used for performing near point interpolation by taking the proportionality coefficient as an interpolation window according to the projection result so as to fill the uncovered equal longitude and latitude grids in the projection result and generate a marine satellite inspection error space distribution map.
In the implementation process, the visual display of the spatial distribution of the inspection errors is realized by generating the spatial distribution diagram of the inspection errors of the marine satellite, and the problem that the spatial distribution of the inspection errors cannot be visually displayed and accurately displayed in the conventional method is solved.
Further, the projection module includes:
the grid division module is used for performing equal longitude and latitude grid division according to the spatial resolution of the ocean satellite data along the track to generate spatial equal longitude and latitude grids, wherein the size of the equal longitude and latitude grids is represented as follows:
Grid_Size＝Spatial_Resolution/100；
wherein, Grid _ Size represents the Size of equal longitude and latitude grids, and the unit is an angle; spatial _ Resolution represents the along-track Spatial Resolution of the marine satellite data, in kilometers;
and the assignment module is used for calculating the row and column index values of the space equal longitude and latitude grids by utilizing the longitude and latitude of the marine satellite data so as to assign the observation error to the corresponding space equal longitude and latitude grids.
In the implementation process, equal longitude and latitude grid setting is carried out according to the resolution of the marine observation satellite along the track, equal longitude and latitude division is carried out on the global range (-90 degrees to 90 degrees and 0 degree to 360 degrees), and the error of the marine observation satellite is projected into the corresponding grid according to the equal longitude and latitude grid.
An embodiment of the present application further provides an electronic device, where the electronic device includes a memory and a processor, where the memory is used to store a computer program, and the processor runs the computer program to make the electronic device execute the method for generating the spatial distribution diagram of inspection error of the marine satellite data described in any one of the above.
The embodiment of the present application further provides a readable storage medium, in which computer program instructions are stored, and when the computer program instructions are read and executed by a processor, the method for generating the spatial distribution diagram of inspection error of the marine satellite data described in any one of the above is performed.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are required to be used in the embodiments of the present application will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered as limiting the scope, and that those skilled in the art can also obtain other related drawings based on the drawings without inventive efforts.
Fig. 1 is a flowchart of a method for generating a spatial distribution diagram of inspection errors of marine satellite data according to an embodiment of the present disclosure;
FIG. 2 is a flow chart of reference data acquisition provided by an embodiment of the present application;
FIG. 3 is a flowchart for calculating an observation error of marine satellite data according to an embodiment of the present disclosure;
FIG. 4 is a flow chart of projection provided by an embodiment of the present application;
fig. 5 is a flowchart for filling up an area where the equal longitude and latitude grid projection does not completely cover, according to an embodiment of the present application;
fig. 6 is a block diagram of a structure of a device for generating a spatial distribution diagram of inspection errors of marine satellite data according to an embodiment of the present disclosure;
fig. 7 is a block diagram of an overall structure of a device for generating a spatial distribution diagram of inspection errors of marine satellite data according to an embodiment of the present disclosure.
Icon:
100-a reference data acquisition module; 101-a parameter acquisition module; 102-a reference data calculation module; 200-an error calculation module; 201-inversion result obtaining module; 202-quality control module; 203-observation error calculation module; 300-a projection module; 301-a meshing module; 302-a first valuation module; 400-a scaling factor calculation module; 500-an interpolation module; 501-a judging module; 502-effective projection acquisition module; 503-second valuation module.
Detailed Description
The technical solutions in the embodiments of the present application will be described below with reference to the drawings in the embodiments of the present application.
It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures. Meanwhile, in the description of the present application, the terms "first", "second", and the like are used only for distinguishing the description, and are not to be construed as indicating or implying relative importance.
Example 1
Referring to fig. 1, fig. 1 is a flowchart of a method for generating a spatial distribution diagram of inspection errors of marine satellite data according to an embodiment of the present disclosure. The method is used for automatically generating the inspection error spatial distribution map, and specifically comprises the following steps:
step S100: performing space-time three-dimensional spline interpolation processing on inspection source data based on parameter information of the marine satellite data to obtain reference data of the marine satellite data;
specifically, as shown in fig. 2, which is a flowchart of reference data acquisition, the step may specifically include:
step S101: acquiring physical parameters of inspection source data;
step S102: and performing space-time three-dimensional spline interpolation processing on the inspection source data based on the physical parameters and the parameter information of the marine satellite data to acquire the reference data.
For example, the mode data such as ECMWF or NCEP can be used as inspection source data, parameter information such as observation time and observation position of the marine satellite data is used as input, a cubic spline interpolation method is used for space-time three-dimensional spline interpolation, and reference data corresponding to the marine satellite data is calculated.
Specifically, reading mode data such as ECMWF or NCEP and the like, specifically comprising time data, space longitude and latitude data, physical parameter data such as sea surface temperature, wind speed U component, wind speed V component, effective wave height, wavelength, wave direction, wave period, atmosphere water vapor content, cloud liquid water content and the like of the mode data, converting the time data into julian days, converting the space longitude and latitude data into a range of 0-360 degrees, and performing matrix conversion on the physical parameters according to the range of 0-360 degrees; and reading the observation time, the observation data longitude and latitude and the like in the marine satellite data, and the like. Converting the observation time into julian days, and converting the longitude and latitude of the observation data into the range of 0-360 degrees; and inputting the time of the mode data, the latitude and longitude grids, the physical parameters, the observation time of the marine satellite data and the observation latitude and longitude, performing space-time three-dimensional spline interpolation processing on the mode data by adopting a cubic spline interpolation method, and calculating to obtain reference data corresponding to the marine satellite data.
In addition, the reference data corresponding to the marine satellite data beyond the time and longitude and latitude grid range of the inspection source data is set to be invalid.
The reference data is calculated by adopting a three-dimensional spline interpolation method, the three-dimensional spline interpolation accords with the space-time change rule of observation elements such as marine dynamic environment, marine ecological environment and the like, so that the calculation of the marine satellite observation reference data is more accurate, the positioning of the reference data is consistent with that of the marine satellite observation data, and the display precision of the marine satellite data inspection error space distribution diagram is improved.
Step S200: calculating an observation error of the marine satellite data by using the reference data;
as shown in fig. 3, to calculate the observation error of the marine satellite data, the steps may include:
step S201: acquiring an ocean parameter inversion result corresponding to the ocean satellite data;
step S202: acquiring a quality identifier corresponding to the ocean parameter inversion result to perform quality control on the ocean parameter inversion result;
step S203: and calculating the observation error of the marine satellite data by using the marine parameter inversion result after the quality control and the reference data.
Reading an ocean parameter inversion result corresponding to ocean satellite data, such as sea surface temperature, sea surface wind speed, sea surface wind direction, effective wave height, wavelength, wave direction, wave period, atmospheric water vapor content, cloud liquid water content and the like; reading a quality identifier corresponding to an ocean parameter inversion result; performing quality control on the ocean parameter inversion result by using the quality identification of the ocean parameter inversion result; and calculating the observation error of the marine satellite data by using the reference data obtained by space-time three-dimensional spline interpolation processing.
In addition, marine satellite observation data needing to be eliminated in the quality identification of the marine parameter inversion result is set to be an invalid value.
Step S300: projecting the observation error on a preset space equal longitude and latitude grid generated by equal longitude and latitude network division based on the marine satellite data to obtain a projection result;
as shown in fig. 4, which is a flowchart of projection, the step may specifically include:
step S301: and performing equal-longitude-latitude grid division according to the spatial resolution of the ocean satellite data along the orbit to generate a spatial equal-longitude-latitude grid, wherein the size of the equal-longitude-latitude grid is expressed as:
Grid_Size＝Spatial_Resolution/100；
wherein, Grid _ Size represents the Size of equal longitude and latitude grids, and the unit is an angle; spatial _ Resolution represents the along-track Spatial Resolution of the marine satellite data, in kilometers;
step S302: and calculating the row and column index values of the space equal longitude and latitude grids by utilizing the longitude and latitude of the marine satellite data so as to assign the observation error to the corresponding space equal longitude and latitude grids.
Exemplarily, the equal-longitude-latitude grid size setting is carried out according to the spatial resolution along the rail of the marine satellite data, for example, the spatial resolution along the rail is 25km, and the equal-longitude-latitude grid is set to be 0.25 degrees; carrying out equal longitude and latitude division on the global longitude and latitude (-90 degrees and 0-360 degrees), wherein the number of latitude grids is 720, and the number of longitude grids is 1440; and calculating row and column index values of the equal longitude and latitude grids by utilizing the longitude and latitude of the marine satellite data, and assigning the marine satellite data observation error to the corresponding row and column grids.
Step S400: calculating a proportionality coefficient between space distances and longitude grids at different latitudes based on the curvature of the earth;
specifically, the curvature of the earth is approximated to be a spherical shape, and the proportional relation between the spatial distance and the longitude grid at different latitudes is calculated through the curvature of the earth, wherein the proportional coefficient is expressed as:
wherein, R represents the proportionality coefficient, ceiling represents rounding up, and Latitude represents Latitude.
Specifically, according to the approximate curvature of the earth, cosine calculation is carried out on different latitudes, absolute values are obtained, and then the calculation results are rounded, so that the proportional coefficient of the integer is obtained.
Step S500: and performing near point interpolation by taking the proportion coefficient as an interpolation window according to the projection result so as to fill the uncovered equal longitude and latitude grids in the projection result and generate a marine satellite inspection error space distribution map.
As shown in fig. 5, for a flowchart of filling up an area where the equal longitude and latitude grid projection coverage is not complete, the step may specifically include:
step S501: circularly judging whether each equal longitude and latitude grid is assigned or not;
step S502: if not, searching by taking the proportional coefficients corresponding to the equal longitude and latitude grids as an interpolation window to obtain the nearest effective projection;
step S503: and assigning the effective projection to the equal longitude and latitude grids.
Circularly traversing the equal longitude and latitude grids, and judging whether the longitude and latitude grids are assigned or not; if the longitude and latitude grid is assigned, jumping to the next longitude and latitude grid; if the longitude and latitude grid is not assigned, searching by taking a proportional coefficient corresponding to the longitude and latitude grid as an interpolation window, namely searching the covered longitude and latitude grid in the interpolation window, recording the grid point distance from the longitude and latitude grid to the uncovered grid, and taking the longitude and latitude grid with the nearest covered grid point distance as an effective projection; searching an effective projection which is closest to the longitude and latitude grid in the interpolation window; and assigning the effective projection value to the longitude and latitude grid.
The method adopts the earth curvature to calculate the proportional relation between the space distance and the longitude grid at different latitudes, the proportional coefficient is the same as the scale of the space coverage loss, the optimization of the interpolation window setting is realized, the search range of the interpolation window is effectively reduced, the space coverage interpolation efficiency is improved, and the automatic generation time of the space distribution diagram of the ocean satellite data inspection error is shortened.
The proportional relation between the spatial distance and the longitude grid at different latitudes is used as an interpolation window, the problem of incomplete coverage of equal longitude and latitude grid projection is effectively solved by adopting a near point interpolation method, a spatially continuous and consistent marine satellite inspection error spatial distribution map is generated, the visualization effect of the marine satellite inspection error spatial distribution map is effectively improved, and the readability of an inspection report is increased.
Example 2
The embodiment of the present application provides a device for generating a spatial distribution diagram of inspection errors of marine satellite data, which is applied to the method for generating a spatial distribution diagram of inspection errors of marine satellite data in embodiment 1, and as shown in fig. 6, is a structural block diagram of the device for generating a spatial distribution diagram of inspection errors of marine satellite data, and the device includes:
the reference data acquisition module 100 is configured to perform space-time three-dimensional spline interpolation processing on inspection source data based on parameter information of marine satellite data to acquire reference data of the marine satellite data;
an error calculation module 200, configured to calculate an observation error of the marine satellite data by using the reference data;
the projection module 300 is configured to project the observation error on a preset space equal-longitude and latitude grid generated by equal-longitude and latitude network division based on the marine satellite data to obtain a projection result;
a scale factor calculation module 400 for calculating scale factors between spatial distances and longitude grids at different latitudes based on the curvature of the earth;
and the interpolation module 500 is configured to perform near point interpolation according to the projection result by using the proportionality coefficient as an interpolation window, so as to fill the uncovered equal longitude and latitude grids in the projection result and generate a marine satellite inspection error spatial distribution map.
As shown in fig. 7, it is a block diagram of the overall structure of the device for generating the spatial distribution diagram of the inspection error of the marine satellite data, wherein the reference data acquiring module 100 includes:
a parameter obtaining module 101, configured to obtain physical parameters of inspection source data;
a reference data calculation module 102, configured to perform space-time three-dimensional spline interpolation processing on the inspection source data based on the physical parameters and the parameter information of the marine satellite data to obtain the reference data.
The error calculation module 200 includes:
an inversion result obtaining module 201, configured to obtain an ocean parameter inversion result corresponding to the ocean satellite data;
the quality control module 202 is configured to obtain a quality identifier corresponding to the ocean parameter inversion result, so as to perform quality control on the ocean parameter inversion result;
and the observation error calculation module 203 is configured to calculate an observation error of the marine satellite data by using the inversion result of the marine parameters after the quality control and the reference data.
Wherein the projection module 300 includes:
the grid division module 301 is configured to perform equal-longitude-latitude grid division according to the spatial resolution of the ocean satellite data along the track, and generate equal-longitude-latitude grids in space, where the equal-longitude-latitude grid size is expressed as:
Grid_Size＝Spatial_Resolution/100；
wherein, Grid _ Size represents the Size of equal longitude and latitude grids, and the unit is an angle; spatial _ Resolution represents the along-track Spatial Resolution of the marine satellite data, in kilometers;
the first assignment module 302 is configured to calculate a rank index value of the space equal longitude and latitude grid by using the longitude and latitude of the marine satellite data, so as to assign the observation error to the corresponding space equal longitude and latitude grid.
The calculation process of the scaling factor calculation module 400 is as follows:
specifically, the proportional relationship between the spatial distance and the longitude grid at different latitudes is calculated through the curvature of the earth, and the proportionality coefficient is expressed as:
wherein, R represents the proportionality coefficient, ceiling represents rounding up, and Latitude represents Latitude.
The interpolation module 500 includes:
a judging module 501, configured to circularly judge whether each equal longitude and latitude grid is assigned;
an effective projection obtaining module 502, configured to search for an interpolation window using the scaling factor corresponding to the equal longitude and latitude grids if no assignment is performed, so as to obtain an effective projection closest to the equal longitude and latitude grid;
a second assigning module 503, configured to assign the effective projection to the equal longitude and latitude grid.
An embodiment of the present application further provides an electronic device, where the electronic device includes a memory and a processor, the memory is used to store a computer program, and the processor runs the computer program to make the electronic device execute the method for generating the spatial distribution diagram of inspection error of marine satellite data according to embodiment 1.
The embodiment of the present application further provides a readable storage medium, in which computer program instructions are stored, and when the computer program instructions are read and executed by a processor, the method for generating the inspection error space distribution map of the marine satellite data according to embodiment 1 is performed.
In the embodiments provided in the present application, it should be understood that the disclosed apparatus and method can be implemented in other ways. The apparatus embodiments described above are merely illustrative, and for example, the flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of apparatus, methods and computer program products according to various embodiments of the present application. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
In addition, functional modules in the embodiments of the present application may be integrated together to form an independent part, or each module may exist separately, or two or more modules may be integrated to form an independent part.
The functions, if implemented in the form of software functional modules and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present application or portions thereof that substantially contribute to the prior art may be embodied in the form of a software product stored in a storage medium and including instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present application. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
The above description is only an example of the present application and is not intended to limit the scope of the present application, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, improvement and the like made within the spirit and principle of the present application shall be included in the protection scope of the present application. It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures.
The above description is only for the specific embodiments of the present application, but the scope of the present application is not limited thereto, and any person skilled in the art can easily conceive of the changes or substitutions within the technical scope of the present application, and shall be covered by the scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.
It is noted that, herein, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other identical elements in a process, method, article, or apparatus that comprises the element.
Claims (10)
1. A method for generating a spatial distribution map of inspection errors of marine satellite data, the method comprising:
performing space-time three-dimensional spline interpolation processing on inspection source data based on parameter information of the marine satellite data to obtain reference data of the marine satellite data;
calculating an observation error of the marine satellite data by using the reference data;
projecting the observation error on a preset space equal longitude and latitude grid generated by equal longitude and latitude network division based on the marine satellite data to obtain a projection result;
calculating a proportionality coefficient between space distances and longitude grids at different latitudes based on the curvature of the earth;
and performing near point interpolation by taking the proportion coefficient as an interpolation window according to the projection result so as to fill the uncovered equal longitude and latitude grids in the projection result and generate a marine satellite inspection error space distribution map.
2. The method as claimed in claim 1, wherein the step of performing spatio-temporal three-dimensional spline interpolation on the inspection source data based on the parameter information of the marine satellite data to obtain the reference data of the marine satellite data comprises:
acquiring physical parameters of inspection source data;
and performing space-time three-dimensional spline interpolation processing on the inspection source data based on the physical parameters and the parameter information of the marine satellite data to acquire the reference data.
3. The method as claimed in claim 1, wherein said calculating the observation error of said marine satellite data using said reference data comprises:
acquiring an ocean parameter inversion result corresponding to the ocean satellite data;
acquiring a quality identifier corresponding to the ocean parameter inversion result to perform quality control on the ocean parameter inversion result;
and calculating the observation error of the marine satellite data by using the marine parameter inversion result after the quality control and the reference data.
4. The method as claimed in claim 1, wherein the step of projecting the observation error on a preset equal longitude and latitude grid generated by equal longitude and latitude grid division based on the marine satellite data to obtain a projection result comprises:
and performing equal-longitude-latitude grid division according to the spatial resolution of the ocean satellite data along the orbit to generate a spatial equal-longitude-latitude grid, wherein the size of the equal-longitude-latitude grid is expressed as:
Grid_Size＝Spatial_Resolution/100；
wherein, Grid _ Size represents the Size of equal longitude and latitude grids, and the unit is an angle; spatial _ Resolution represents the along-track Spatial Resolution of the marine satellite data, in kilometers;
and calculating the row and column index values of the space equal longitude and latitude grids by utilizing the longitude and latitude of the marine satellite data so as to assign the observation error to the corresponding space equal longitude and latitude grids.
5. A method for generating a spatial distribution map of inspection errors for marine satellite data as claimed in claim 1, wherein said computing a scaling factor between spatial distance and longitude grids at different latitudes based on earth curvature comprises:
calculating the proportional relation between space distances and longitude grids at different latitudes through the curvature of the earth, wherein the proportional coefficient is expressed as:
wherein, R represents the proportionality coefficient, ceiling represents rounding up, and Latitude represents Latitude.
6. The method as claimed in claim 4, wherein the performing a near point interpolation according to the projection result with the scaling factor as an interpolation window to fill in the uncovered equal longitude and latitude grids in the projection result and generate the marine satellite inspection error spatial distribution map comprises:
circularly judging whether each equal longitude and latitude grid is assigned or not;
if not, searching by taking the proportional coefficients corresponding to the equal longitude and latitude grids as an interpolation window to obtain the nearest effective projection;
and assigning the effective projection to the equal longitude and latitude grids.
7. An apparatus for generating a spatial distribution map of inspection errors of marine satellite data, said apparatus comprising:
the reference data acquisition module is used for carrying out space-time three-dimensional spline interpolation processing on the inspection source data based on the parameter information of the marine satellite data so as to acquire the reference data of the marine satellite data;
the error calculation module is used for calculating the observation error of the marine satellite data by using the reference data;
the projection module is used for projecting the observation error on a preset space equal longitude and latitude grid generated by equal longitude and latitude network division based on the marine satellite data so as to obtain a projection result;
the scale factor calculation module is used for calculating scale factors between space distances and longitude grids at different latitudes based on the curvature of the earth;
and the interpolation module is used for performing near point interpolation by taking the proportionality coefficient as an interpolation window according to the projection result so as to fill the uncovered equal longitude and latitude grids in the projection result and generate a marine satellite inspection error space distribution map.
8. The device for generating the spatial distribution map of inspection errors of marine satellite data as claimed in claim 7, wherein the projection module comprises:
the grid division module is used for performing equal longitude and latitude grid division according to the spatial resolution of the ocean satellite data along the track to generate spatial equal longitude and latitude grids, wherein the size of the equal longitude and latitude grids is represented as follows:
Grid_Size＝Spatial_Resolution/100；
wherein, Grid _ Size represents the Size of equal longitude and latitude grids, and the unit is an angle; spatial _ Resolution represents the along-track Spatial Resolution of the marine satellite data, in kilometers;
and the assignment module is used for calculating the row and column index values of the space equal longitude and latitude grids by utilizing the longitude and latitude of the marine satellite data so as to assign the observation error to the corresponding space equal longitude and latitude grids.
9. An electronic device, characterized in that the electronic device comprises a memory for storing a computer program and a processor for executing the computer program to cause the electronic device to perform the method of verification of the spatial distribution of error of marine satellite data according to any of claims 1 to 6.
10. A readable storage medium, wherein computer program instructions are stored in the readable storage medium, and when the computer program instructions are read and executed by a processor, the method for generating the spatial distribution map of inspection error of the marine satellite data according to any one of claims 1 to 6 is performed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011596962.XA CN112612916B (en) | 2020-12-29 | 2020-12-29 | Method and device for generating inspection error space distribution diagram of marine satellite data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011596962.XA CN112612916B (en) | 2020-12-29 | 2020-12-29 | Method and device for generating inspection error space distribution diagram of marine satellite data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112612916A true CN112612916A (en) | 2021-04-06 |
| CN112612916B CN112612916B (en) | 2024-02-06 |
Family
ID=75249084
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011596962.XA Active CN112612916B (en) | 2020-12-29 | 2020-12-29 | Method and device for generating inspection error space distribution diagram of marine satellite data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112612916B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114494811A (en) * | 2022-02-07 | 2022-05-13 | 国家海洋环境预报中心 | Method and system for fusing abnormal height data of satellite along sea level |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040225665A1 (en) * | 2003-05-09 | 2004-11-11 | Microsoft Corporation | System and method for employing a grid index for location and precision encoding |
| US20070003911A1 (en) * | 2005-06-29 | 2007-01-04 | Julien Serre | Method and system for cartographic projection of the terrestrial globe and map produced by this method |
| JP2008058109A (en) * | 2006-08-30 | 2008-03-13 | Central Res Inst Of Electric Power Ind | Observation data estimation method and observation data estimation program |
| CN102033898A (en) * | 2010-09-27 | 2011-04-27 | 华东师范大学 | Extraction method for local cloud cover information metadata of moderate resolution imaging spectral image |
| CN103324941A (en) * | 2013-06-19 | 2013-09-25 | 鲁东大学 | Remote sensing classification pattern spot boundary precision evaluation method based on close distance |
| CN103942737A (en) * | 2014-05-09 | 2014-07-23 | 国家电网公司 | Drawing method of historical forest fire distribution of power transmission line |
| CN104268429A (en) * | 2014-10-15 | 2015-01-07 | 湖北大学 | Satellite-borne SAR (Synthetic Aperture Radar) based offshore wind energy resource remote sensing method and system |
| US20170105627A1 (en) * | 2015-10-19 | 2017-04-20 | Biosense Webster (Israel) Ltd. | Illustrating error in a temperature distribution map |
| JP2019020260A (en) * | 2017-07-18 | 2019-02-07 | アイサンテクノロジー株式会社 | Parameter distribution system |
| JP2019104432A (en) * | 2017-12-14 | 2019-06-27 | 日本電気株式会社 | Correction device, system, correction method, and program |
| CN110531444A (en) * | 2019-08-29 | 2019-12-03 | 中国气象局广州热带海洋气象研究所（广东省气象科学研究所） | The error source of numerical weather prediction model determines method and device |
| KR20200059085A (en) * | 2018-11-20 | 2020-05-28 | 서울대학교산학협력단 | A Method for Sea Surface Temperature Retrieval using Surface Drifter Temperature Data and Satellite Infrared Images |
| CN111639149A (en) * | 2020-05-29 | 2020-09-08 | 山东浪潮通软信息科技有限公司 | Ocean data visualization method and device |
-
2020
- 2020-12-29 CN CN202011596962.XA patent/CN112612916B/en active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040225665A1 (en) * | 2003-05-09 | 2004-11-11 | Microsoft Corporation | System and method for employing a grid index for location and precision encoding |
| US20070003911A1 (en) * | 2005-06-29 | 2007-01-04 | Julien Serre | Method and system for cartographic projection of the terrestrial globe and map produced by this method |
| JP2008058109A (en) * | 2006-08-30 | 2008-03-13 | Central Res Inst Of Electric Power Ind | Observation data estimation method and observation data estimation program |
| CN102033898A (en) * | 2010-09-27 | 2011-04-27 | 华东师范大学 | Extraction method for local cloud cover information metadata of moderate resolution imaging spectral image |
| CN103324941A (en) * | 2013-06-19 | 2013-09-25 | 鲁东大学 | Remote sensing classification pattern spot boundary precision evaluation method based on close distance |
| CN103942737A (en) * | 2014-05-09 | 2014-07-23 | 国家电网公司 | Drawing method of historical forest fire distribution of power transmission line |
| CN104268429A (en) * | 2014-10-15 | 2015-01-07 | 湖北大学 | Satellite-borne SAR (Synthetic Aperture Radar) based offshore wind energy resource remote sensing method and system |
| US20170105627A1 (en) * | 2015-10-19 | 2017-04-20 | Biosense Webster (Israel) Ltd. | Illustrating error in a temperature distribution map |
| JP2019020260A (en) * | 2017-07-18 | 2019-02-07 | アイサンテクノロジー株式会社 | Parameter distribution system |
| JP2019104432A (en) * | 2017-12-14 | 2019-06-27 | 日本電気株式会社 | Correction device, system, correction method, and program |
| KR20200059085A (en) * | 2018-11-20 | 2020-05-28 | 서울대학교산학협력단 | A Method for Sea Surface Temperature Retrieval using Surface Drifter Temperature Data and Satellite Infrared Images |
| CN110531444A (en) * | 2019-08-29 | 2019-12-03 | 中国气象局广州热带海洋气象研究所（广东省气象科学研究所） | The error source of numerical weather prediction model determines method and device |
| CN111639149A (en) * | 2020-05-29 | 2020-09-08 | 山东浪潮通软信息科技有限公司 | Ocean data visualization method and device |
Non-Patent Citations (6)
| Title |
|---|
| 周游等: "基于Himawari-8静止气象卫星的输电线路山火监测与告警技术", 《高电压技术》, vol. 46, no. 7, pages 2561 - 2569 * |
| 杜振彩: ""滑动窗区空间相关系数加权集合方法及其在IPCC-AR4多模式集合模拟和预测中的应用", 《大气科学》, vol. 34, no. 8, pages 1168 - 1186 * |
| 王雨生: "基于多源遥感数据的高时空分辨率地表温度产品生成方法研究", 中国优秀硕士学位论文全文数据库 基础科学辑, pages 009 - 46 * |
| 翟振和;孙中苗;肖云;李芳;: "自主海洋测高卫星串飞模式的设计与重力场反演精度分析", 武汉大学学报(信息科学版), no. 07 * |
| 范雕: "卫星测高重力数据反演海底地形的理论和方法研究", 中国优秀硕士学位论文全文数据库 基础科学辑, pages 008 - 77 * |
| 赵亮等: "一种基于卫星遥感与地面测站数据融合技术的雪深动态反演方法", 气象学报, vol. 71, no. 4, pages 769 - 782 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114494811A (en) * | 2022-02-07 | 2022-05-13 | 国家海洋环境预报中心 | Method and system for fusing abnormal height data of satellite along sea level |
| CN114494811B (en) * | 2022-02-07 | 2023-03-24 | 国家海洋环境预报中心 | Method and system for fusing abnormal height data of satellite along sea level |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112612916B (en) | 2024-02-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | Google Earth elevation data extraction and accuracy assessment for transportation applications | |
| Böhm et al. | Atmospheric effects in space geodesy | |
| Frigg et al. | An assessment of the foundational assumptions in high-resolution climate projections: the case of UKCP09 | |
| Dąbrowski et al. | Integration of multi-source geospatial data from GNSS receivers, terrestrial laser scanners, and unmanned aerial vehicles | |
| Suyunov et al. | Field studies of electronic total stations in a special reference satellite geodetic basis | |
| US9626796B2 (en) | Method to optimize the visualization of a map's projection based on data and tasks | |
| Alcaras et al. | The influence of interpolated point location and density on 3D bathymetric models generated by kriging methods: an application on the Giglio Island Seabed (Italy) | |
| CN111611540A (en) | Image control point elevation precise cloud computing conversion method based on thousand seeking positions | |
| Wilson et al. | Seamless bathymetry and topography datasets for New South Wales, Australia | |
| Bures et al. | River bathymetry model based on floodplain topography | |
| CN116205541A (en) | Method and device for evaluating influence of local pollution source on environmental air quality | |
| CN112612916B (en) | Method and device for generating inspection error space distribution diagram of marine satellite data | |
| Alcaras et al. | The importance of the coordinate transformation process in using heterogeneous data in coastal and marine geographic information system | |
| Kastrisios et al. | Increasing Efficiency of Nautical Chart Production and Accessibility to Marine Environment Data through an Open-Science Compilation Workflow | |
| Schweitzer et al. | A method for analysis of spatial uncertainty in image based surface velocimetry | |
| Su et al. | Bounding the residual tropospheric error by interval analysis | |
| Lee et al. | Analyzing precision and efficiency of global navigation satellite system-derived height determination for coastal and island areas | |
| Kalaitzis et al. | Interactive web mapping applications for 2D and 3D geo-visualization of persistent scatterer interferometry SAR data | |
| Masetti et al. | Automated identification of discrepancies between nautical charts and survey soundings | |
| Liu et al. | On the study of influences of different factors on the rapid tropospheric tomography | |
| Krdžalić et al. | A precise geoid model of Bosnia and Herzegovina by the KTH method and its validation | |
| Wang et al. | Local oceanic vertical deflection determination with gravity data along a profile | |
| CN112987010A (en) | System and method for multi-radar mapping of robot | |
| Seo et al. | Quality assessment of linear data | |
| Lewicka et al. | Analysis of Transformation Methods of Hydroacoustic and Optoelectronic Data Based on the Tombolo Measurement Campaign in Sopot |
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 |