CN112612916B - Method and device for generating inspection error space distribution diagram of marine satellite data - Google Patents
Method and device for generating inspection error space distribution diagram of marine satellite data Download PDFInfo
- Publication number
- CN112612916B CN112612916B CN202011596962.XA CN202011596962A CN112612916B CN 112612916 B CN112612916 B CN 112612916B CN 202011596962 A CN202011596962 A CN 202011596962A CN 112612916 B CN112612916 B CN 112612916B
- Authority
- CN
- China
- Prior art keywords
- longitude
- latitude
- marine satellite
- satellite data
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000007689 inspection Methods 0.000 title claims abstract description 50
- 238000010586 diagram Methods 0.000 title claims description 27
- 238000012545 processing Methods 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims description 25
- 238000012360 testing method Methods 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 11
- 238000003908 quality control method Methods 0.000 claims description 9
- 238000012795 verification Methods 0.000 claims description 8
- 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
- 230000007812 deficiency Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000006467 substitution reaction 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
- 238000007796 conventional method 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
- 230000000007 visual effect Effects 0.000 description 1
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
Abstract
The embodiment of the application provides a method and a device for generating an inspection error space distribution map of ocean satellite data, and relates to the technical field of ocean engineering. Performing 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 an observation error of the marine satellite data by using the reference data; projecting the observation error on a space equal longitude and latitude grid generated by dividing a preset equal longitude and latitude network based on marine satellite data so as to obtain a projection result; calculating a scaling factor between the spatial distance and the longitude grid at different latitudes based on the curvature of the earth; and carrying out adjacent point interpolation by taking the proportional coefficient as an interpolation window according to the projection result so as to fill in uncovered equal longitude and latitude grids in the projection result and generate a marine satellite inspection error space distribution map, thereby solving the problem that the prior method cannot intuitively display and accurately display the space distribution of inspection errors.
Description
Technical Field
The application relates to the technical field of ocean engineering, in particular to a method and a device for generating an inspection error space distribution diagram of ocean satellite data.
Background
With the rapid development of aerospace technology, remote sensing technology and computer technology, marine satellite remote sensing is in a new stage of updating and rapid development. The ocean satellite has the characteristics of large-scale and high-precision global ocean observation, and the application fields of the ocean satellite comprise ocean ecological environment monitoring, ocean power environment forecasting, ocean disaster prevention and reduction, ocean rights and interests maintenance and the like. The precision of the marine satellite product is the basis of marine satellite remote sensing application, and the good product quality directly determines the data application effect, but the existing method cannot intuitively and accurately display the spatial distribution of the inspection errors.
Disclosure of Invention
The embodiment of the application aims to provide a method and a device for generating a space distribution diagram of an inspection error of marine satellite data, which are used for accurately calculating and visually displaying the inspection error in a form of the space distribution diagram of the inspection error, and solve the problem that the existing method cannot visually display and accurately display the space distribution of the inspection error.
The embodiment of the application provides a method for generating a test error space distribution diagram of marine satellite data, which comprises the following steps:
performing space-time three-dimensional spline interpolation processing on the 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 using the reference data;
projecting the observation errors 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;
calculating a scaling factor between the spatial distance and the longitude grid at different latitudes based on the curvature of the earth;
and carrying out adjacent point interpolation by taking the proportionality coefficient as an interpolation window according to the projection result so as to fill in uncovered equal longitude and latitude grids in the projection result and generate a marine satellite inspection error space distribution diagram.
In the implementation process, an automatic generation method of the marine satellite data inspection error space distribution diagram is provided, the calculation of the observation error of the marine satellite data is carried out 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 the positioning of the marine satellite data, and the display precision of the marine satellite data inspection error space distribution diagram is improved; the method has the advantages that the proportional coefficient between the space distance and the longitude grid in different dimensions is the same as the dimension of the space coverage deficiency, the optimization of the setting of an interpolation window is realized, the searching 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 insufficient coverage of equal longitude and latitude grid projection is effectively solved by adopting a method of adjacent point interpolation, 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 improved, and the problem that the space distribution of inspection errors cannot be intuitively 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 includes:
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 so as to acquire the reference data.
In the implementation process, a cubic spline interpolation method is adopted to conduct space-time three-dimensional spline interpolation processing on the inspection source data, and reference data corresponding to the marine satellite data is obtained through calculation.
Further, the calculating the observation error of the marine satellite data using the reference data includes:
obtaining an ocean parameter inversion result corresponding to the ocean satellite data;
acquiring a quality identifier corresponding to the ocean parameter inversion result so as to control the quality of the ocean parameter inversion result;
and calculating the observation error of the marine satellite data by using the marine parameter inversion result after quality control and the reference data.
In the implementation process, the observation errors of the marine satellite data are calculated based on the reference data, and the marine satellite observation data which need to be removed in the marine parameter inversion result quality identification are set as invalid values, so that the accuracy of the observation error calculation result is improved.
Further, the projecting the observation error on a spatial equal longitude and latitude grid generated by dividing the equal longitude and latitude network based on the marine satellite data to obtain a projection result includes:
performing equal longitude and latitude grid division according to the orbital spatial resolution of the marine satellite data to generate a spatial equal longitude and latitude grid, wherein the equal longitude and latitude grid size is expressed as:
Grid_Size＝Spatial_Resolution/100；
wherein grid_size represents the Size of the equal warp and weft grids, and the unit is an angle; spatial_resolution represents the orbital Spatial Resolution of the marine satellite data in kilometers;
and calculating row-column index values of the space equal-longitude-latitude grids by using longitude and latitude of the ocean satellite data so as to assign the observation errors to the corresponding space equal-longitude-latitude grids.
In the implementation process, longitude and latitude grid division such as global (-90 degrees to 90 degrees, 0 degrees to 360 degrees) and the like is carried out according to the spatial resolution of the marine satellite data, and space and other longitude and latitude grid projection is carried out on the observation error data.
Further, the calculating the scaling factor between the spatial distance and the longitude grid on different latitudes based on the curvature of the earth includes:
calculating the proportional relation between the spatial distance and the longitude grid on different latitudes through the curvature of the earth, wherein the proportional coefficient is expressed as follows:
wherein R represents the proportionality coefficient, ceiling represents an upward rounding, and Latitude represents a Latitude.
In the implementation process, the curvature of the earth is approximately spherical, the proportional relation between the space distance and the longitude grid on different latitudes is calculated, and the proportional coefficient directly calculated by the latitudes is rounded to obtain the integral proportional coefficient.
Further, the performing interpolation of adjacent points according to the projection result by using the scaling factor as an interpolation window to fill in uncovered equal longitude and latitude grids in the projection result and generate a marine satellite inspection error space distribution diagram, including:
circularly judging whether each equal warp and weft grid is assigned;
if not, searching by taking the proportional coefficient corresponding to the equal warp and weft grids as an interpolation window to obtain effective projection closest to the interpolation window;
assigning the effective projection to the equal warp and weft grid.
In the implementation process, the equal longitude and latitude grids are traversed, scaling coefficients on different latitudes are used as interpolation windows, adjacent point interpolation is carried out, and areas with incomplete coverage of the equal longitude and latitude grid projections are filled.
The embodiment of the application also provides a device for generating the test error space distribution map of the marine satellite data, which comprises:
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;
an error calculation module for calculating an observation error of the marine satellite data 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 scaling factor calculation module is used for calculating scaling factors between the space distances and longitude grids on different latitudes based on the curvature of the earth;
and the interpolation module is used for carrying out adjacent point interpolation by taking the proportionality coefficient as an interpolation window according to the projection result so as to fill in uncovered equal longitude and latitude grids in the projection result and generate a marine satellite inspection error space distribution diagram.
In the implementation process, the visual display of the space distribution of the inspection errors is realized by generating the space distribution diagram of the inspection errors of the marine satellite, and the problem that the space distribution of the inspection errors cannot be intuitively and accurately displayed in the conventional method is solved.
Further, the projection module includes:
the grid division module is used for carrying out equal-longitude and latitude grid division according to the orbital spatial resolution of the marine satellite data to generate space equal-longitude and latitude grids, wherein the equal-longitude and latitude grids are expressed as:
Grid_Size＝Spatial_Resolution/100；
wherein grid_size represents the Size of the equal warp and weft grids, and the unit is an angle; spatial_resolution represents the orbital Spatial Resolution of the marine satellite data in kilometers;
and the assignment module is used for calculating row-column index values of the space equal-longitude-latitude grids by using the longitude and latitude of the ocean satellite data so as to assign the observation errors to the corresponding space equal-longitude-latitude grids.
In the implementation process, the equal longitude and latitude grid arrangement is carried out according to the orbital resolution of the ocean observation satellite, the equal longitude and latitude division is carried out on the global scope (-90 degrees to 90 degrees and 0 degrees to 360 degrees), and the ocean observation satellite error is projected into the corresponding grid according to the equal longitude and latitude grid.
The embodiment of the application also provides electronic equipment, which comprises a memory and a processor, wherein the memory is used for storing a computer program, and the processor runs the computer program to enable the electronic equipment to execute the method for generating the verification error space distribution map of the marine satellite data.
The embodiment of the application also provides a readable storage medium, wherein the readable storage medium stores computer program instructions, and when the computer program instructions are read and executed by a processor, the method for generating the verification error space distribution map of the marine satellite data is executed.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are needed 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 should not be considered as limiting the scope, and other related drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flowchart of a method for generating a test error spatial distribution map of marine satellite data according to an embodiment of the present application;
FIG. 2 is a flowchart of reference data acquisition provided in an embodiment of the present application;
FIG. 3 is a flow chart of calculating the observed error of marine satellite data according to an embodiment of the present application;
FIG. 4 is a flow chart of projection provided by an embodiment of the present application;
FIG. 5 is a flowchart of filling an area with insufficient coverage of a grid projection of equal warp and weft according to an embodiment of the present application;
FIG. 6 is a block diagram of a device for generating a test error spatial distribution map of marine satellite data according to an embodiment of the present application;
fig. 7 is an overall block diagram of a device for generating a test error space distribution map of marine satellite data according to an embodiment of the present application.
Icon:
100-a reference data acquisition module; 101-a parameter acquisition module; 102-a reference data calculation module; 200-an error calculation module; 201, an inversion result acquisition module; 202-a quality control module; 203, an observation error calculation module; 300-a projection module; 301-a meshing module; 302-a first assignment module; 400-a scaling factor calculation module; 500-an interpolation module; 501-a judging module; 502-an effective projection acquisition module; 503-a second assignment 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 numerals and letters denote like items in the following figures, and thus once an item is defined in one figure, no further definition or explanation thereof is necessary in the following figures. Meanwhile, in the description of the present application, the terms "first", "second", and the like are used only to distinguish 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 verification error spatial distribution map of marine satellite data according to an embodiment of the present application. The method is used for automatically generating a test error space distribution diagram and specifically comprises the following steps of:
step S100: performing space-time three-dimensional spline interpolation processing on the 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, a flowchart of reference data acquisition 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 so as to acquire the reference data.
For example, pattern data such as ECMWF or NCEP may be used as test 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 to perform space-time three-dimensional spline interpolation, and reference data corresponding to the marine satellite data is calculated.
The method specifically comprises the steps of reading mode data such as ECMWF or NCEP, 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 0-360 degrees; the observation time, the longitude and latitude of the observation data, and the like in the marine satellite data are read, and the same is true. Converting the observation time into julian days, and converting the longitude and latitude of the observation data into a range of 0-360 degrees; the time, longitude and latitude grids, physical parameters, observation time and observation longitude and latitude of the pattern data are taken as inputs, a cubic spline interpolation method is adopted to conduct space-time three-dimensional spline interpolation processing on the pattern data, and reference data corresponding to the ocean satellite data are obtained through calculation.
Further, the reference data corresponding to the marine satellite data exceeding the time of checking the source data and the latitude and longitude grid range is set to an invalid value.
By adopting the three-dimensional spline interpolation method to calculate the reference data, the three-dimensional spline interpolation accords with the space-time change rule of the 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 the positioning 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 using the reference data;
as shown in fig. 3, for a flowchart for calculating the observed error of the marine satellite data, the steps may include:
step S201: obtaining an ocean parameter inversion result corresponding to the ocean satellite data;
step S202: acquiring a quality identifier corresponding to the ocean parameter inversion result so as to control the quality of 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 quality control and the reference data.
The sea parameter inversion results corresponding to the sea satellite data are read, such as sea surface temperature, sea surface wind speed, sea surface wind direction, effective wave height, wavelength, wave direction, wave period, atmosphere water vapor content, cloud liquid water content and the like; reading a quality identifier corresponding to the 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 the space-time three-dimensional spline interpolation processing.
In addition, the ocean satellite observation data to be removed in the quality identification of the ocean parameter inversion result is set to be an invalid value.
Step S300: projecting the observation errors 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;
as shown in fig. 4, which is a projection flowchart, the steps may specifically include:
step S301: performing equal longitude and latitude grid division according to the orbital spatial resolution of the marine satellite data to generate a spatial equal longitude and latitude grid, wherein the equal longitude and latitude grid size is expressed as:
Grid_Size＝Spatial_Resolution/100；
wherein grid_size represents the Size of the equal warp and weft grids, and the unit is an angle; spatial_resolution represents the orbital Spatial Resolution of the marine satellite data in kilometers;
step S302: and calculating row-column index values of the space equal-longitude-latitude grids by using longitude and latitude of the ocean satellite data so as to assign the observation errors to the corresponding space equal-longitude-latitude grids.
Illustratively, the equal longitude and latitude grid size is set according to the orbital spatial resolution of the marine satellite data, for example, the orbital spatial resolution is 25km, and the equal longitude and latitude grid is set to 0.25 °; dividing the longitude and latitude of the world (-90 degrees to 90 degrees, 0 degrees to 360 degrees) into 720 latitude grids and 1440 longitude grids; and calculating row and column index values of the equal longitude and latitude grids by using longitude and latitude of the ocean satellite data, and assigning the ocean satellite data observation errors to the corresponding row and column grids.
Step S400: calculating a scaling factor between the spatial distance and the longitude grid at different latitudes based on the curvature of the earth;
specifically, the curvature of the earth is approximated to be a sphere, and the proportional relation between the spatial distances and the longitude grids on different latitudes is calculated through the curvature of the earth, wherein the proportional coefficient is expressed as follows:
wherein R represents the proportionality coefficient, ceiling represents an upward rounding, and Latitude represents a Latitude.
Specifically, according to the approximate earth curvature, cosine calculation is carried out on different latitudes, absolute values are taken, and then the calculation result is rounded to obtain the integral proportionality coefficient.
Step S500: and carrying out adjacent point interpolation by taking the proportionality coefficient as an interpolation window according to the projection result so as to fill in uncovered equal longitude and latitude grids in the projection result and generate a marine satellite inspection error space distribution diagram.
As shown in fig. 5, to fill up the area with insufficient coverage of the warp and weft grid projection, the steps may specifically include:
step S501: circularly judging whether each equal warp and weft grid is assigned;
step S502: if not, searching by taking the proportional coefficient corresponding to the equal warp and weft grids as an interpolation window to obtain effective projection closest to the interpolation window;
step S503: assigning the effective projection to the equal warp and weft grid.
Circularly traversing the 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 grids are not assigned, searching by taking a proportional coefficient corresponding to the longitude and latitude grids as an interpolation window, namely searching for covered longitude and latitude grids in the interpolation window, recording the grid point distance from the longitude and latitude grids to uncovered grids, and obtaining the longitude and latitude grids with the nearest grid point distance from the covered grid point to be effective projection; searching the effective projection closest to the longitude and latitude grid in the interpolation window; the effective projection value is assigned to the longitude and latitude grid.
And calculating the proportional relation between the spatial distances and the longitude grids on different latitudes by adopting the earth curvature, wherein the proportional coefficient is the same as the scale of the space coverage deficiency, so that the optimization of the setting of an interpolation window is realized, the searching range of the interpolation window is effectively reduced, the space coverage interpolation efficiency is improved, and the automatic generation time of the marine satellite data inspection error space distribution map is shortened.
The proportional relation between the space distances and the longitude grids on different latitudes is used as an interpolation window, a method of adjacent point interpolation is adopted, the problem of incomplete coverage of equal longitude and latitude grid projections is effectively solved, a 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, and the readability of an inspection report is improved.
Example 2
An embodiment of the present application provides a device for generating a test error spatial distribution map of marine satellite data, which is applied to the method for generating a test error spatial distribution map of marine satellite data in embodiment 1, as shown in fig. 6, and is a structural block diagram of the device for generating a test error spatial distribution map of marine satellite data, and the device includes:
the reference data acquisition module 100 is used for performing 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;
an error calculation module 200 for calculating an observation error of the marine satellite data using the reference data;
the projection module 300 is configured to project the observation error on a spatial equal longitude and latitude grid generated by dividing the space based on the marine satellite data by equal longitude and latitude network, so as to obtain a projection result;
a scaling factor calculation module 400 for calculating scaling factors between spatial distances at different latitudes and longitude grids based on the curvature of the earth;
and the interpolation module 500 is used for carrying out adjacent point interpolation by taking the proportionality coefficient as an interpolation window according to the projection result so as to fill in the uncovered equal longitude and latitude grids in the projection result and generate a marine satellite inspection error space distribution diagram.
As shown in fig. 7, an overall block diagram of the apparatus for generating a test error spatial distribution map of marine satellite data, wherein the reference data acquisition module 100 includes:
a parameter acquisition module 101, configured to acquire physical parameters of inspection source data;
and the reference data calculation module 102 is used for 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 so as to acquire the reference data.
The error calculation module 200 includes:
an inversion result obtaining module 201, configured to obtain an inversion result of the marine parameter corresponding to the marine satellite data;
the quality control module 202 is configured to obtain a quality identifier corresponding to the marine parameter inversion result, so as to perform quality control on the marine parameter inversion result;
and the observation error calculation module 203 is configured to calculate an observation error of the marine satellite data using the quality-controlled marine parameter inversion result and the reference data.
Wherein the projection module 300 comprises:
the meshing module 301 is configured to perform equal-longitude and latitude meshing according to the orbital spatial resolution of the marine satellite data, and generate a spatial equal-longitude and latitude mesh, where the equal-longitude and latitude mesh size is expressed as:
Grid_Size＝Spatial_Resolution/100；
wherein grid_size represents the Size of the equal warp and weft grids, and the unit is an angle; spatial_resolution represents the orbital Spatial Resolution of the marine satellite data in kilometers;
and the first assignment module 302 is configured to calculate a rank index value of the space equal longitude and latitude grid by using 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 relation between the spatial distance and the longitude grid at different latitudes is calculated through the curvature of the earth, and the proportional coefficient is expressed as:
wherein R represents the proportionality coefficient, ceiling represents an upward rounding, and Latitude represents a Latitude.
The interpolation module 500 includes:
a judging module 501, configured to circularly judge whether each of the equal warp and weft grids performs assignment;
the effective projection obtaining module 502 is configured to search for an interpolation window by using the scaling coefficient corresponding to the equal longitude and latitude grid if no assignment is performed, so as to obtain an effective projection closest to the interpolation window;
a second assignment module 503, configured to assign the effective projection to the equal warp and weft grid.
The embodiment of the application further provides an electronic device, which includes a memory and a processor, where the memory is configured to store a computer program, and the processor is configured to execute the computer program to cause the electronic device to execute the method for generating the verification error space distribution map of the marine satellite data according to the embodiment 1.
The present embodiment also provides a readable storage medium having stored therein computer program instructions which, when read and executed by a processor, perform the method for generating a verification error spatial distribution map of marine satellite data according to embodiment 1.
In the several embodiments provided in this application, it should be understood that the disclosed apparatus and method may be implemented in other manners as well. The apparatus embodiments described above are merely illustrative, for example, flow diagrams 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, the functional modules in the embodiments of the present application may be integrated together to form a single part, or each module may exist alone, or two or more modules may be integrated to form a single 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 may be embodied essentially or in a part contributing to the prior art or in a part of the technical solution, in the form of a software product stored in a storage medium, including several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods described in 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, random Access Memory), a magnetic disk, or an optical disk, or other various media capable of storing program codes.
The foregoing is merely exemplary embodiments of the present application and is not intended to limit the scope of the present application, and various modifications and variations may be suggested to one skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the present application should be included in the protection scope of the present application. It should be noted that: like reference numerals and letters denote like items in the following figures, and thus once an item is defined in one figure, no further definition or explanation thereof is necessary in the following figures.
The foregoing is merely 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 think about changes or substitutions within the technical scope of the present application, and the changes and substitutions are intended to 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 relational terms such as first and second, and the like are 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. Moreover, 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 one … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
Claims (8)
1. A method for generating a test error spatial distribution map of marine satellite data, the method comprising:
performing space-time three-dimensional spline interpolation processing on the inspection source data based on parameter information of the marine satellite data to obtain reference data of the marine satellite data, specifically, performing space-time three-dimensional spline interpolation processing on the pattern data by using time, longitude and latitude grids, physical parameters, observation time and observation longitude and latitude of the pattern data as inputs and calculating to obtain reference data corresponding to the marine satellite data, wherein the pattern data is the inspection source data;
calculating an observation error of the marine satellite data by using the reference data, and specifically, acquiring a marine parameter inversion result corresponding to the marine satellite data; acquiring a quality identifier corresponding to the ocean parameter inversion result so as to control the quality of the ocean parameter inversion result; calculating the observation error of the marine satellite data by using the marine parameter inversion result after quality control and the reference data;
projecting the observation errors 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;
calculating a scaling factor between the spatial distance and the longitude grid at different latitudes based on the curvature of the earth;
and carrying out adjacent point interpolation by taking the proportionality coefficient as an interpolation window according to the projection result so as to fill in uncovered equal longitude and latitude grids in the projection result and generate a marine satellite inspection error space distribution diagram.
2. The method for generating a test error spatial distribution map of marine satellite data according to claim 1, wherein projecting the observation error on a spatial 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:
performing equal longitude and latitude grid division according to the orbital spatial resolution of the marine satellite data to generate a spatial equal longitude and latitude grid, wherein the equal longitude and latitude grid size is expressed as:
Grid_Size=Spatial_Resolution/100；
wherein grid_size represents the Size of the equal warp and weft grids, and the unit is an angle; spatial_resolution represents the orbital Spatial Resolution of the marine satellite data in kilometers;
and calculating row-column index values of the space equal-longitude-latitude grids by using longitude and latitude of the ocean satellite data so as to assign the observation errors to the corresponding space equal-longitude-latitude grids.
3. The method for generating a test error spatial distribution map of marine satellite data according to claim 1, wherein calculating a scaling factor between spatial distances at different latitudes and a longitude grid based on the curvature of the earth comprises:
calculating the proportional relation between the spatial distance and the longitude grid on different latitudes through the curvature of the earth, wherein the proportional coefficient is expressed as follows:
；
wherein,Rthe scale factor is indicated as such and,ceilingthe representation is rounded up and down to the top,Latituderepresenting the latitude.
4. The method for generating a test error spatial distribution map of marine satellite data according to claim 2, wherein said performing, based on the projection result, a neighboring interpolation with the scaling factor as an interpolation window to fill in uncovered equal longitude and latitude grids in the projection result and generating the test error spatial distribution map of marine satellite comprises:
circularly judging whether each equal warp and weft grid is assigned;
if not, searching by taking the proportional coefficient corresponding to the equal warp and weft grids as an interpolation window to obtain effective projection closest to the interpolation window;
assigning the effective projection to the equal warp and weft grid.
5. A device for generating a test error spatial distribution map of marine satellite data, said device 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 parameter information of the marine satellite data to acquire reference data of the marine satellite data, specifically, taking time, longitude and latitude grids, physical parameters of pattern data, observation time and observation longitude and latitude of the marine satellite data as inputs, carrying out space-time three-dimensional spline interpolation processing on the pattern data by adopting a cubic spline interpolation method, and calculating to obtain reference data corresponding to the marine satellite data, wherein the pattern data is the inspection source data;
an error calculation module for calculating an observation error of the marine satellite data using the reference data, the error calculation module comprising: the inversion result acquisition module is used for acquiring an ocean parameter inversion result corresponding to the ocean satellite data; the quality control module is used for acquiring a quality identifier corresponding to the ocean parameter inversion result so as to control the quality of the ocean parameter inversion result; the observation error calculation module is used for calculating the observation error of the marine satellite data by using the marine parameter inversion result after quality control and 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 scaling factor calculation module is used for calculating scaling factors between the space distances and longitude grids on different latitudes based on the curvature of the earth;
and the interpolation module is used for carrying out adjacent point interpolation by taking the proportionality coefficient as an interpolation window according to the projection result so as to fill in uncovered equal longitude and latitude grids in the projection result and generate a marine satellite inspection error space distribution diagram.
6. The apparatus for generating a verification error spatial distribution map of marine satellite data according to claim 5, wherein said projection module comprises:
the grid division module is used for carrying out equal-longitude and latitude grid division according to the orbital spatial resolution of the marine satellite data to generate space equal-longitude and latitude grids, wherein the equal-longitude and latitude grids are expressed as:
Grid_Size=Spatial_Resolution/100；
wherein grid_size represents the Size of the equal warp and weft grids, and the unit is an angle; spatial_resolution represents the orbital Spatial Resolution of the marine satellite data in kilometers;
and the assignment module is used for calculating row-column index values of the space equal-longitude-latitude grids by using the longitude and latitude of the ocean satellite data so as to assign the observation errors to the corresponding space equal-longitude-latitude grids.
7. An electronic device comprising a memory for storing a computer program and a processor that runs the computer program to cause the electronic device to perform the method of generating a verification error spatial profile of marine satellite data according to any one of claims 1 to 4.
8. A readable storage medium having stored therein computer program instructions which, when read and executed by a processor, perform the method of generating a verification error spatial profile of marine satellite data as claimed in any one of claims 1 to 4.
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 CN112612916A (en) | 2021-04-06 |
| CN112612916B true 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) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114494811B (en) * | 2022-02-07 | 2023-03-24 | 国家海洋环境预报中心 | Method and system for fusing abnormal height data of satellite along sea level |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| 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 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7966301B2 (en) * | 2003-05-09 | 2011-06-21 | Planeteye Company Ulc | System and method for employing a grid index for location and precision encoding |
| FR2888023B1 (en) * | 2005-06-29 | 2008-02-29 | Webraska Mobile Technologies S | LAND GLOBE MAPPING PROJECTION METHOD AND MAP PRODUCED ACCORDING TO THE METHOD |
| US10448838B2 (en) * | 2015-10-19 | 2019-10-22 | Biosense Webster (Israel) Ltd. | Illustrating error in a temperature distribution map |
-
2020
- 2020-12-29 CN CN202011596962.XA patent/CN112612916B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| 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 (7)
| Title |
|---|
| "滑动窗区空间相关系数加权集合方法及其在IPCC-AR4多模式集合模拟和预测中的应用;杜振彩;《大气科学》;第34卷(第8期);第1168-1186页 * |
| 一种基于卫星遥感与地面测站数据融合技术的雪深动态反演方法;赵亮等;气象学报;第71卷(第4期);第769-782页 * |
| 卫星测高重力数据反演海底地形的理论和方法研究;范雕;中国优秀硕士学位论文全文数据库 基础科学辑;第A008-77页 * |
| 基于Himawari-8静止气象卫星的输电线路山火监测与告警技术;周游等;《高电压技术》;第46卷(第7期);第2561-2569页 * |
| 基于多源遥感数据的高时空分辨率地表温度产品生成方法研究;王雨生;中国优秀硕士学位论文全文数据库 基础科学辑;第A009-46页 * |
| 廖玉芳等.《湖南省现代农业气候区划》.湖南大学出版社,2010,第135-136页. * |
| 自主海洋测高卫星串飞模式的设计与重力场反演精度分析;翟振和;孙中苗;肖云;李芳;;武汉大学学报(信息科学版)(07);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112612916A (en) | 2021-04-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Fisher | Improved modeling of elevation error with geostatistics | |
| US10134152B2 (en) | Method and system for determining cells traversed by a measuring or visualization axis | |
| US9378462B2 (en) | Probability mapping system | |
| CN102346923B (en) | Data grading organization method based on longitude and latitude grid | |
| CN106709127B (en) | Extraction of pollutant emission source using CFD and satellite data | |
| Cannata et al. | Two-dimensional dam break flooding simulation: a GIS-embedded approach | |
| US20060228019A1 (en) | System and method for three dimensional change detection and measurement of a scene using change analysis | |
| Deshpande | Improved floodplain delineation method using high‐density LiDAR data | |
| CN112612916B (en) | Method and device for generating inspection error space distribution diagram of marine satellite data | |
| CN116205541B (en) | Method and device for evaluating influence of local pollution source on environmental air quality | |
| Varela-González et al. | Performance testing of LiDAR exploitation software | |
| Su et al. | Bounding the residual tropospheric error by interval analysis | |
| CN111524220B (en) | Method, system, electronic device, and medium for sectioning three-dimensional model | |
| Zhu et al. | Feature-constrained automatic geometric deformation analysis method of bridge models toward digital twin | |
| RU2796370C1 (en) | System and method for dynamic visualization of environmental pollution | |
| Chang et al. | Modeling spatial variations of sea level rise and corresponding inundation impacts: A case study for Florida, USA | |
| Krdžalić et al. | A precise geoid model of Bosnia and Herzegovina by the KTH method and its validation | |
| JP2021089158A (en) | Wind direction and velocity estimation device, wind direction and velocity estimation method, and wind direction and velocity estimation program | |
| Elmore et al. | Aggregation Methods Using Bathymetry Sources of Differing Subjective Reliabilities for Navigation Mapping | |
| Sultanov et al. | Development of a processing method of digital maps and satellite images for solving problems of emergencies | |
| US7869951B1 (en) | Method of determining local electromagnetic sensor refractivity | |
| Lytvyn et al. | Aerospace pictures processing by means of interstripation of functions of two variables | |
| Pérez et al. | Efficient methods to convert LiDAR-derived ellipsoid heights to orthometric heights | |
| Lambrechts et al. | Geographic information systems and remote sensing | |
| Sayed et al. | Deep Learning Approach for Detecting Digital Elevation Model (DEM) Uncertainty to Enhance Assessment of Water Resources |
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 |