CN114034651A - Method and device for generating global earth surface spectrum basic data - Google Patents
Method and device for generating global earth surface spectrum basic data Download PDFInfo
- Publication number
- CN114034651A CN114034651A CN202111325330.4A CN202111325330A CN114034651A CN 114034651 A CN114034651 A CN 114034651A CN 202111325330 A CN202111325330 A CN 202111325330A CN 114034651 A CN114034651 A CN 114034651A
- Authority
- CN
- China
- Prior art keywords
- basic data
- data
- spectrum
- earth surface
- global earth
- 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.)
- Pending
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 124
- 238000000034 method Methods 0.000 title claims abstract description 59
- 238000012545 processing Methods 0.000 claims abstract description 17
- 238000005070 sampling Methods 0.000 claims abstract description 10
- 230000001131 transforming effect Effects 0.000 claims abstract description 4
- 238000002310 reflectometry Methods 0.000 claims description 42
- 238000004590 computer program Methods 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 238000013500 data storage Methods 0.000 description 7
- 238000013461 design Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 4
- 230000005855 radiation Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/55—Specular reflectivity
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Image Processing (AREA)
Abstract
The invention relates to the technical field of computers, in particular to a method and a device for generating global earth surface spectrum basic data. The method comprises the following steps: selecting a preset MODIS data product as initial global earth surface spectrum basic data; repairing invalid pixels and cloud pixels in the initial global earth surface spectrum basic data to obtain repaired first basic data; extracting and projective transforming the first basic data in the HDF format by using MRT software to obtain second basic data projected by the geographic coordinates; the format of the second basic data is tif format; performing space scale consistency processing on second basic data of spring, summer, autumn and winter by adopting a nearest neighbor sampling method so as to enable the spatial resolution of the second basic data to be the same; and splicing the second basic data after the spatial scale unification treatment based on the wave bands of the second basic data to generate the target global earth surface spectrum basic data. The scheme can reduce the storage space of global earth surface spectrum basic data.
Description
Technical Field
The embodiment of the invention relates to the technical field of computers, in particular to a method and a device for generating global earth surface spectrum basic data.
Background
The global earth surface spectrum basic data is not only an important input for simulating a downward-looking detection scene of a hyperspectral satellite (namely, the first step of generating a background image), but also core data influencing the precision of background simulation; the global earth surface spectrum basic data comprise spectrum emissivity and spectrum reflectivity.
In the prior art, global earth surface spectrum basic data is stored in a manner of 500m spatial resolution. And a 500m by 500m image with 0.4-12.5 μm wave band occupies 2.25G of storage space, and if the global earth surface spectrum basic data is required to be stored, the storage space occupies 5.84T, which is not acceptable for computer hardware and retrieval speed.
Therefore, there is a need for a method and an apparatus for generating global earth surface spectrum basic data to solve the above technical problems.
Disclosure of Invention
The embodiment of the invention provides a method and a device for generating global earth surface spectrum basic data, which can reduce the storage space of the global earth surface spectrum basic data.
In a first aspect, an embodiment of the present invention provides a method for generating global earth surface spectrum basic data, including:
selecting a preset MODIS data product as initial global earth surface spectrum basic data; the initial global earth surface spectrum basic data comprise spectrum emissivity and spectrum reflectivity, the spectrum emissivity corresponds to three mid-infrared wave bands and three thermal infrared wave bands, and the spectrum reflectivity corresponds to seven wave bands;
repairing invalid pixels and cloud pixels in the initial global earth surface spectrum basic data to obtain repaired first basic data;
extracting and projective transforming the first basic data in the HDF format by using MRT software to obtain second basic data projected by the geographic coordinates; wherein the format of the second basic data is tif format;
performing spatial scale uniformization processing on the second basic data in spring, summer, autumn and winter by adopting a nearest neighbor sampling method so as to enable the spatial resolution of the second basic data to be the same;
and splicing the second basic data after the spatial scale conformance processing by adopting a ratio difference method based on the wave band of the second basic data to generate target global earth surface spectrum basic data.
In one possible design, the selecting of the preset MODIS data product as the initial global earth surface spectrum basic data includes:
selecting a seven-waveband earth surface reflectivity data product of 8-day synthetic data of MOD/MYD09A1 as the spectrum reflectivity of the initial global earth surface spectrum basic data;
the three-waveband surface emissivity data products of thermal infrared of 8-day synthetic data of MOD/MYD21A2 and the three-waveband surface emissivity data products of middle infrared of 8-day synthetic data of MOD/MYD11B2 are selected as the spectrum emissivity of the initial global surface spectrum basic data.
In one possible design, the first base data is obtained by the following formula:
wherein, imgy,t(i, j) is the pixel with the coordinate of (i, j) of the image to be restored, y and t are the year and date of obtaining the image to be restored respectively, imgy-1,t(i, j) is the pixel with coordinate (i, j) in the image of the ith date of the y-1 year in the same area, imgy+1,t(i, j) is the pixel with the coordinate (i, j) in the (y + 1) th date image of the same area, imgy,t-7(i, j) is the pixel with coordinates (i, j) in the t-7 th date image of the same area in the y year, imgy,t+7And (i, j) is a pixel with the coordinate (i, j) in the image of the ith +7 th date of the y year in the same area.
In one possible design, the spatial resolution of the second base data after the spatial scale unification process is 500 m.
In a second aspect, an embodiment of the present invention further provides an apparatus for generating global earth surface spectrum basic data, including:
the initialization module is used for selecting a preset MODIS data product as initial global earth surface spectrum basic data; the initial global earth surface spectrum basic data comprise spectrum emissivity and spectrum reflectivity, the spectrum emissivity corresponds to three mid-infrared wave bands and three thermal infrared wave bands, and the spectrum reflectivity corresponds to seven wave bands;
the restoration module is used for restoring the invalid pixel and the cloud pixel in the initial global earth surface spectrum basic data to obtain restored first basic data;
the projection conversion module is used for extracting and projection converting the first basic data in the HDF format by using MRT software to obtain second basic data projected by the geographic coordinates; wherein the format of the second basic data is tif format;
the consistency processing module is used for carrying out space scale consistency processing on the second basic data in spring, summer, autumn and winter by adopting a nearest neighbor sampling method so as to enable the spatial resolution of the second basic data to be the same;
and the generating module is used for splicing the second basic data after the spatial scale conformance processing by adopting a ratio difference method based on the wave band of the second basic data to generate the target global earth surface spectrum basic data.
In one possible design, the initialization module is configured to:
selecting a seven-waveband earth surface reflectivity data product of 8-day synthetic data of MOD/MYD09A1 as the spectrum reflectivity of the initial global earth surface spectrum basic data;
the three-waveband surface emissivity data products of thermal infrared of 8-day synthetic data of MOD/MYD21A2 and the three-waveband surface emissivity data products of middle infrared of 8-day synthetic data of MOD/MYD11B2 are selected as the spectrum emissivity of the initial global surface spectrum basic data.
In one possible design, the first base data is obtained by the following formula:
wherein, imgy,t(i, j) is the pixel with the coordinate of (i, j) of the image to be restored, y and t are the year and date of obtaining the image to be restored respectively, imgy-1,t(i, j) is the pixel with coordinate (i, j) in the image of the ith date of the y-1 year in the same area, imgy+1,t(i, j) is the pixel with the coordinate (i, j) in the (y + 1) th date image of the same area, imgy,t-7(i, j) is the pixel with coordinates (i, j) in the t-7 th date image of the same area in the y year, imgy,t+7And (i, j) is a pixel with the coordinate (i, j) in the image of the ith +7 th date of the y year in the same area.
In one possible design, the spatial resolution of the second base data after the spatial scale unification process is 500 m.
In a third aspect, an embodiment of the present invention further provides a computing device, including a memory and a processor, where the memory stores a computer program, and the processor, when executing the computer program, implements the method described in any one of the above.
In a fourth aspect, the present invention further provides a computer-readable storage medium, on which a computer program is stored, and when the computer program is executed in a computer, the computer program causes the computer to execute any one of the methods described above.
The embodiment of the invention provides a method and a device for generating global earth surface spectrum basic data, wherein the method for generating the global full-spectrum basic data covering 0.459-12.270 mu m is established by taking reflectivity and emissivity data of thirteen wave bands of a preset MODIS data product as the basis and sequentially carrying out methods such as preprocessing, spatial scale consistence, data splicing and the like. The method solves the problem of global earth surface spectrum basic data generation, and the data of thirteen wave bands are used for replacing the original high-spectrum storage mode, so that the global earth surface spectrum basic data storage space is reduced by 46 times, namely the global earth surface spectrum basic data storage space is reduced.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
FIG. 1 is a flow chart of a method for generating global earth surface spectrum basic data according to an embodiment of the present invention;
FIG. 2 is a diagram of a hardware architecture of a computing device according to an embodiment of the present invention;
fig. 3 is a schematic structural diagram of an apparatus for generating global earth surface spectrum basic data according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer and more complete, the technical solutions in the embodiments of the present invention will be described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention, and based on the embodiments of the present invention, all other embodiments obtained by a person of ordinary skill in the art without creative efforts belong to the scope of the present invention.
Referring to fig. 1, an embodiment of the present invention provides a method for generating global earth surface spectrum basic data, including:
102, repairing invalid pixels and cloud pixels in the initial global earth surface spectrum basic data to obtain repaired first basic data;
104, extracting and projective transforming the first basic data in the HDF format by using MRT software to obtain second basic data projected by the geographic coordinates; wherein the format of the second basic data is tif format;
106, performing spatial scale consistency processing on second basic data of spring, summer, autumn and winter by adopting a nearest neighbor sampling method so as to enable the spatial resolution of the second basic data to be the same;
and step 108, splicing the second basic data after the spatial scale unification treatment by adopting a ratio difference method based on the wave band of the second basic data to generate the target global earth surface spectrum basic data.
In the embodiment of the invention, a method for covering global full-spectrum basic data of 0.459-12.270 mu m is established by taking reflectivity and emissivity data of thirteen wave bands of a preset MODIS data product as a basis and sequentially carrying out methods such as preprocessing, spatial scale uniformization, data splicing and the like. The method solves the problem of global earth surface spectrum basic data generation, and the data of thirteen wave bands are used for replacing the original high-spectrum storage mode, so that the global earth surface spectrum basic data storage space is reduced by 46 times, namely the global earth surface spectrum basic data storage space is reduced.
The manner in which the various steps shown in fig. 1 are performed is described below.
For step 100, step 100 may specifically include:
selecting a preset MODIS data product as initial global earth surface spectrum basic data, wherein the method comprises the following steps:
selecting a seven-waveband earth surface reflectivity data product of 8-day synthetic data of MOD/MYD09A1 as the spectrum reflectivity of the initial global earth surface spectrum basic data;
the three-waveband surface emissivity data products of thermal infrared of 8-day synthetic data of MOD/MYD21A2 and the three-waveband surface emissivity data products of middle infrared of 8-day synthetic data of MOD/MYD11B2 are selected as the spectrum emissivity of the initial global surface spectrum basic data.
Aiming at the real requirements of full spectrum simulation and simulation of earth surface scenes in a large global area, the wave band coverage of basic data is determined to be 0.4-12.5 mu m, the spatial resolution is 500m, the scene change caused by the weather factors is considered, and the basic data of global background radiation covers four seasons of spring, summer, autumn and winter.
And aiming at the relevant requirements during the selection of the basic data, on the basis of the analysis of the index parameters and the data products of the existing satellite data, selecting an MODIS data product as the basic data to realize the generation and the storage of the global earth surface spectrum basic data. In the aspect of surface reflectivity, 1-7 waveband surface reflectivity data products of 8-day synthetic data of MOD/MYD09A1 and the like are selected, the spatial resolution of the data is 500m, and the reflectivity product data are shown in Table 1.
TABLE 1 MODIS surface reflectivity data case table
In the aspect of surface emissivity, the MODIS remote sensor has 6 channels to obtain surface emissivity data of mid-infrared and thermal-infrared bands, and the corresponding band information of the obtained emissivity data is shown in table 2.
Table 2 wavelength information of obtaining band by MODIS emissivity data
MOD/MYD21A2 and MOD/MYD11B2 surface emissivity products were selected that were 8-day synthetic data, including surface emissivity data for six bands of mid-infrared (B20, B22, B23), thermal infrared (B29, B31, B32), with spatial resolutions ranging from 1km to 6km, and MODIS surface emissivity data as shown in Table 3.
TABLE 3 MODIS surface emissivity data condition table
In the aspect of data phase and season phase selection, four time nodes of 4 months and 15 days, 7 months and 15 days, 10 months and 15 days and 1 month and 15 days respectively represent four seasons of spring, summer, autumn and winter, and data MODIS reflectivity and emissivity of corresponding time nodes are respectively obtained.
Aiming at the step 102, part of the pixels in the MODIS surface reflectivity and emissivity product are invalid pixels and only have filling values. Meanwhile, under the influence of cloud, part of pixels have no real values of surface reflectance values and emissivity values. Therefore, through filling value retrieval, the QC data layer and the state data layer (the QC data layer and the state data layer are respectively two data layers in MODIS data) are positioned as the cloud pixels, and data in the same period of the previous year and adjacent time before and after are used for repairing the invalid pixels and the cloud pixels.
The first basic data is obtained by the following formula:
wherein, imgy,t(i, j) is the pixel with the coordinate of (i, j) of the image to be restored, y and t are the year and date of obtaining the image to be restored respectively, imgy-1,t(i, j) is the pixel with coordinate (i, j) in the image of the ith date of the y-1 year in the same area, imgy+1,t(i, j) is the pixel with the coordinate (i, j) in the (y + 1) th date image of the same area, imgy,t-7(i, j) is the pixel with coordinates (i, j) in the t-7 th date image of the same area in the y year, imgy,t+7And (i, j) is a pixel with the coordinate (i, j) in the image of the ith +7 th date of the y year in the same area.
In this embodiment, an effective pixel can be obtained by repairing the invalid pixel and the cloud pixel.
In step 104, the MODIS reflectivity and emissivity product generally includes a plurality of data layers, such as data quality, cloud coverage, observation geometry, transit time, and the like, in addition to the data product layer, and in order to obtain the global earth surface spectrum basic data, the corresponding data layer needs to be extracted from the corresponding data product file, and the MODIS data product is generally a sinusoidal projection, and unlike the conventional geographic coordinate projection mode, data layer extraction and projection conversion need to be performed.
And (2) adopting MODIS Reproduction Tool (MRT) software to extract a reflectivity and emissivity data layer for an MODIS product under HDF-EOS format and sinusoidal projection, outputting the reflectivity and emissivity data layer to a common tif and other data format, and simultaneously carrying out re-projection and re-sampling, wherein projection information is a WGS84 reference plane, Geographic longitude and latitude, and strip data with the nearest re-sampling type and the spatial resolution of 0.005 DEG are adopted.
For step 106, the MOD/MYD09a1 surface reflectivity product contains reflectivity data of seven bands including B01-07, and the spatial resolution is 500m, while the MOD/MYD11B2 surface reflectivity product contains medium wave emissivity data of three bands including Emis _20, Emis _22 and Emis _23, and the spatial resolution is 6km, and the MOD/MYD21a2 surface reflectivity product contains long wave emissivity data of three bands including Emis _29, Emis _31 and Emis _32, and the spatial resolution is 1 km.
In order to meet the requirement of global data spatial resolution, a nearest neighbor sampling method is adopted to carry out spatial scale conversion on reflectivity, emissivity, classified products and the like, so that various data products are unified into 500m spatial resolution. Meanwhile, the reflectivity and emissivity data of four seasons of spring, summer, autumn and winter are processed by referring to the same method.
Aiming at step 108, 13 waveband images including B01-07, Emis _20, Emis _22, Emis _23, Emis _29, Emis _31 and Emis _32 after the spatial scale is uniformized are respectively subjected to image splicing to form 13 independent data layers in the global range. And longitudinally associating the 13 independent data layers according to the geographic coordinates to synthesize global background radiation basic data. The basic data comprises the basic radiation data with the global scope of 0.46-12.01 mu m, and the seasons are applicable to four seasons of spring, summer, autumn and winter.
In subsequent application, basic data, surface coverage classification and a reference spectrum library of 13 corresponding wave bands are extracted according to the selected longitude and latitude, and hyperspectral data of corresponding spatial resolution is generated through ratio-interpolation calculation, so that big data storage and retrieval are avoided. The data storage format is shown in table 4, for example:
table 4 data storage format
In conclusion, the method for generating and storing the global earth surface spectrum basic data can generate the global earth surface spectrum basic data and obviously reduce the storage capacity of global earth surface radiation mass basic data. The invention can provide global earth surface spectrum data for satellite downward-looking infrared camera detection scene simulation software and support satellite downward-looking infrared camera detection scene earth background simulation.
As shown in fig. 2 and fig. 3, the embodiment of the invention provides a device for generating global earth surface spectrum basic data. The device embodiments may be implemented by software, or by hardware, or by a combination of hardware and software. From a hardware aspect, as shown in fig. 2, a hardware architecture diagram of a computing device in which a global earth surface spectrum basic data generation apparatus provided for the embodiment of the invention is located is shown, in addition to the processor, the memory, the network interface, and the nonvolatile memory shown in fig. 2, the computing device in which the apparatus is located in the embodiment may also generally include other hardware, such as a forwarding chip responsible for processing a message, and the like. Taking a software implementation as an example, as shown in fig. 3, as a logical apparatus, a CPU of a computing device in which the apparatus is located reads a corresponding computer program in a non-volatile memory into a memory to run. The embodiment provides a generation device of global earth surface spectrum basic data, wherein the spatial target geometric model comprises at least two first geometric bodies for characterizing a spatial target main body and two second geometric bodies for characterizing a spatial target cell, and the device comprises:
the initialization module is used for selecting a preset MODIS data product as initial global earth surface spectrum basic data; the initial global earth surface spectrum basic data comprise spectrum emissivity and spectrum reflectivity, the spectrum emissivity corresponds to three mid-infrared wave bands and three thermal infrared wave bands, and the spectrum reflectivity corresponds to seven wave bands;
the restoration module is used for restoring the invalid pixel and the cloud pixel in the initial global earth surface spectrum basic data to obtain restored first basic data;
the projection conversion module is used for extracting and projection converting the first basic data in the HDF format by using MRT software to obtain second basic data projected by the geographic coordinates; wherein the format of the second basic data is tif format;
the consistency processing module is used for carrying out space scale consistency processing on the second basic data in spring, summer, autumn and winter by adopting a nearest neighbor sampling method so as to enable the spatial resolution of the second basic data to be the same;
and the generating module is used for splicing the second basic data after the spatial scale conformance processing by adopting a ratio difference method based on the wave band of the second basic data to generate the target global earth surface spectrum basic data.
In an embodiment of the present invention, the initialization module is configured to perform the following operations:
selecting a seven-waveband earth surface reflectivity data product of 8-day synthetic data of MOD/MYD09A1 as the spectrum reflectivity of the initial global earth surface spectrum basic data;
the three-waveband surface emissivity data products of thermal infrared of 8-day synthetic data of MOD/MYD21A2 and the three-waveband surface emissivity data products of middle infrared of 8-day synthetic data of MOD/MYD11B2 are selected as the spectrum emissivity of the initial global surface spectrum basic data.
In an embodiment of the present invention, the first basic data is obtained by the following formula:
wherein, imgy,t(i, j) is the pixel with the coordinate of (i, j) of the image to be repaired, and y and t are respectively the pixelTaking year and date of image to be restored, imgy-1,t(i, j) is the pixel with coordinate (i, j) in the image of the ith date of the y-1 year in the same area, imgy+1,t(i, j) is the pixel with the coordinate (i, j) in the (y + 1) th date image of the same area, imgy,t-7(i, j) is the pixel with coordinates (i, j) in the t-7 th date image of the same area in the y year, imgy,t+7And (i, j) is a pixel with the coordinate (i, j) in the image of the ith +7 th date of the y year in the same area.
In one embodiment of the present invention, the spatial resolution of the second base data after the spatial scale unification process is 500 m.
It is to be understood that the illustrated structure of the embodiment of the present invention does not constitute a specific limitation to a global earth surface spectrum basic data generation device. In other embodiments of the present invention, a global earth surface spectrum basis data generation apparatus may include more or fewer components than those shown, or combine certain components, or split certain components, or a different arrangement of components. The illustrated components may be implemented in hardware, software, or a combination of software and hardware.
Because the content of information interaction, execution process, and the like among the modules in the device is based on the same concept as the method embodiment of the present invention, specific content can be referred to the description in the method embodiment of the present invention, and is not described herein again.
The embodiment of the invention also provides a computing device, which comprises a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to realize the generation method of the global earth surface spectrum basic data in any embodiment of the invention.
Embodiments of the present invention further provide a computer-readable storage medium, on which a computer program is stored, where the computer program, when executed by a processor, causes the processor to execute a method for generating global earth surface spectrum basic data according to any embodiment of the present invention.
Specifically, a system or an apparatus equipped with a storage medium on which software program codes that realize the functions of any of the above-described embodiments are stored may be provided, and a computer (or a CPU or MPU) of the system or the apparatus is caused to read out and execute the program codes stored in the storage medium.
In this case, the program code itself read from the storage medium can realize the functions of any of the above-described embodiments, and thus the program code and the storage medium storing the program code constitute a part of the present invention.
Examples of the storage medium for supplying the program code include a floppy disk, a hard disk, a magneto-optical disk, an optical disk (e.g., CD-ROM, CD-R, CD-RW, DVD-ROM, DVD-RAM, DVD-RW, DVD + RW), a magnetic tape, a nonvolatile memory card, and a ROM. Alternatively, the program code may be downloaded from a server computer via a communications network.
Further, it should be clear that the functions of any one of the above-described embodiments may be implemented not only by executing the program code read out by the computer, but also by causing an operating system or the like operating on the computer to perform a part or all of the actual operations based on instructions of the program code.
Further, it is to be understood that the program code read out from the storage medium is written to a memory provided in an expansion board inserted into the computer or to a memory provided in an expansion module connected to the computer, and then causes a CPU or the like mounted on the expansion board or the expansion module to perform part or all of the actual operations based on instructions of the program code, thereby realizing the functions of any of the above-described embodiments.
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 similar elements in a process, method, article, or apparatus that comprises the element.
Those of ordinary skill in the art will understand that: all or part of the steps for realizing the method embodiments can be completed by hardware related to program instructions, the program can be stored in a computer readable storage medium, and the program executes the steps comprising the method embodiments when executed; and the aforementioned storage medium includes: various media that can store program codes, such as ROM, RAM, magnetic or optical disks.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (10)
1. A method for generating global earth surface spectrum basic data is characterized by comprising the following steps:
selecting a preset MODIS data product as initial global earth surface spectrum basic data; the initial global earth surface spectrum basic data comprise spectrum emissivity and spectrum reflectivity, the spectrum emissivity corresponds to three mid-infrared wave bands and three thermal infrared wave bands, and the spectrum reflectivity corresponds to seven wave bands;
repairing invalid pixels and cloud pixels in the initial global earth surface spectrum basic data to obtain repaired first basic data;
extracting and projective transforming the first basic data in the HDF format by using MRT software to obtain second basic data projected by the geographic coordinates; wherein the format of the second basic data is tif format;
performing spatial scale uniformization processing on the second basic data in spring, summer, autumn and winter by adopting a nearest neighbor sampling method so as to enable the spatial resolution of the second basic data to be the same;
and splicing the second basic data after the spatial scale conformance processing by adopting a ratio difference method based on the wave band of the second basic data to generate target global earth surface spectrum basic data.
2. The method according to claim 1, wherein the selecting of the preset MODIS data product as the initial global earth surface spectrum basic data comprises:
selecting a seven-waveband earth surface reflectivity data product of 8-day synthetic data of MOD/MYD09A1 as the spectrum reflectivity of the initial global earth surface spectrum basic data;
the three-waveband surface emissivity data products of thermal infrared of 8-day synthetic data of MOD/MYD21A2 and the three-waveband surface emissivity data products of middle infrared of 8-day synthetic data of MOD/MYD11B2 are selected as the spectrum emissivity of the initial global surface spectrum basic data.
3. The method of claim 1, wherein the first base data is obtained by the following formula:
wherein, imgy,t(i, j) is the pixel with the coordinate of (i, j) of the image to be restored, y and t are the year and date of obtaining the image to be restored respectively, imgy-1,t(i, j) is the pixel with coordinate (i, j) in the image of the ith date of the y-1 year in the same area, imgy+1,t(i, j) is the pixel with the coordinate (i, j) in the (y + 1) th date image of the same area, imgy,t-7(i, j) is the pixel with coordinates (i, j) in the t-7 th date image of the same area in the y year, imgy,t+7And (i, j) is a pixel with the coordinate (i, j) in the image of the ith +7 th date of the y year in the same area.
4. The method according to claim 1, wherein the spatial resolution of the second base data after the spatial scale unification processing is 500 m.
5. An apparatus for generating global earth surface spectrum basic data, comprising:
the initialization module is used for selecting a preset MODIS data product as initial global earth surface spectrum basic data; the initial global earth surface spectrum basic data comprise spectrum emissivity and spectrum reflectivity, the spectrum emissivity corresponds to three mid-infrared wave bands and three thermal infrared wave bands, and the spectrum reflectivity corresponds to seven wave bands;
the restoration module is used for restoring the invalid pixel and the cloud pixel in the initial global earth surface spectrum basic data to obtain restored first basic data;
the projection conversion module is used for extracting and projection converting the first basic data in the HDF format by using MRT software to obtain second basic data projected by the geographic coordinates; wherein the format of the second basic data is tif format;
the consistency processing module is used for carrying out space scale consistency processing on the second basic data in spring, summer, autumn and winter by adopting a nearest neighbor sampling method so as to enable the spatial resolution of the second basic data to be the same;
and the generating module is used for splicing the second basic data after the spatial scale conformance processing by adopting a ratio difference method based on the wave band of the second basic data to generate the target global earth surface spectrum basic data.
6. The apparatus of claim 5, wherein the initialization module is configured to:
selecting a seven-waveband earth surface reflectivity data product of 8-day synthetic data of MOD/MYD09A1 as the spectrum reflectivity of the initial global earth surface spectrum basic data;
the three-waveband surface emissivity data products of thermal infrared of 8-day synthetic data of MOD/MYD21A2 and the three-waveband surface emissivity data products of middle infrared of 8-day synthetic data of MOD/MYD11B2 are selected as the spectrum emissivity of the initial global surface spectrum basic data.
7. The apparatus of claim 1, wherein the first basic data is obtained by the following formula:
wherein, imgy,t(i, j) is the pixel with the coordinate of (i, j) of the image to be restored, y and t are the year and date of obtaining the image to be restored respectively, imgy-1,t(i, j) is the pixel with coordinate (i, j) in the image of the ith date of the y-1 year in the same area, imgy+1,t(i, j) is the pixel with the coordinate (i, j) in the (y + 1) th date image of the same area, imgy,t-7(i, j) is the pixel with coordinates (i, j) in the t-7 th date image of the same area in the y year, imgy,t+7And (i, j) is a pixel with the coordinate (i, j) in the image of the ith +7 th date of the y year in the same area.
8. The apparatus according to claim 1, wherein the spatial resolution of the second base data after the spatial scale unification process is 500 m.
9. A computing device comprising a memory having stored therein a computer program and a processor that, when executing the computer program, implements the method of any of claims 1-4.
10. A computer-readable storage medium, on which a computer program is stored which, when executed in a computer, causes the computer to carry out the method of any one of claims 1-4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111325330.4A CN114034651A (en) | 2021-11-10 | 2021-11-10 | Method and device for generating global earth surface spectrum basic data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111325330.4A CN114034651A (en) | 2021-11-10 | 2021-11-10 | Method and device for generating global earth surface spectrum basic data |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114034651A true CN114034651A (en) | 2022-02-11 |
Family
ID=80143870
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111325330.4A Pending CN114034651A (en) | 2021-11-10 | 2021-11-10 | Method and device for generating global earth surface spectrum basic data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114034651A (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040164235A1 (en) * | 2003-01-21 | 2004-08-26 | Miller Steven D. | System and method for significant dust detection and enhancement of dust images over land and ocean |
| US20180000339A1 (en) * | 2016-06-29 | 2018-01-04 | Ace Vision Group, Inc. | System and methods using real-time predictive virtual 3d eye finite element modeling for simulation of ocular structure biomechanics |
| CN110348107A (en) * | 2019-07-08 | 2019-10-18 | 中国农业科学院农业资源与农业区划研究所 | A kind of pixel real surface temperature rebuilding method under cloud |
| US20200141877A1 (en) * | 2018-11-06 | 2020-05-07 | Nanjing Agricultural University | Method for estimating aboveground biomass of rice based on multi-spectral images of unmanned aerial vehicle |
| CN112163523A (en) * | 2020-09-29 | 2021-01-01 | 北京环境特性研究所 | Abnormal target detection method and device and computer readable medium |
| AU2021102705A4 (en) * | 2021-05-20 | 2021-07-08 | Dongguan University Of Technology | THz-TDS image defocus processing method based on deep learning |
| CN113256493A (en) * | 2021-05-28 | 2021-08-13 | 北京环境特性研究所 | Thermal infrared remote sensing image reconstruction method and device |
-
2021
- 2021-11-10 CN CN202111325330.4A patent/CN114034651A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040164235A1 (en) * | 2003-01-21 | 2004-08-26 | Miller Steven D. | System and method for significant dust detection and enhancement of dust images over land and ocean |
| US20180000339A1 (en) * | 2016-06-29 | 2018-01-04 | Ace Vision Group, Inc. | System and methods using real-time predictive virtual 3d eye finite element modeling for simulation of ocular structure biomechanics |
| US20200141877A1 (en) * | 2018-11-06 | 2020-05-07 | Nanjing Agricultural University | Method for estimating aboveground biomass of rice based on multi-spectral images of unmanned aerial vehicle |
| CN110348107A (en) * | 2019-07-08 | 2019-10-18 | 中国农业科学院农业资源与农业区划研究所 | A kind of pixel real surface temperature rebuilding method under cloud |
| CN112163523A (en) * | 2020-09-29 | 2021-01-01 | 北京环境特性研究所 | Abnormal target detection method and device and computer readable medium |
| AU2021102705A4 (en) * | 2021-05-20 | 2021-07-08 | Dongguan University Of Technology | THz-TDS image defocus processing method based on deep learning |
| CN113256493A (en) * | 2021-05-28 | 2021-08-13 | 北京环境特性研究所 | Thermal infrared remote sensing image reconstruction method and device |
Non-Patent Citations (11)
| Title |
|---|
| JING CONG: "A Pilot Study: Changes of Intestinal Microbiota of Patients With Non-small Cell Lung Cancer in Response to Osimertinib Therapy", 《SYSTEMS MICROBIOLOGY》, vol. 11, pages 1 - 4 * |
| XIDONG CHEN: "An Assessment of Water Color for Inland Water in China Using a Landsat 8-Derived Forel–Ule Index and the Google Earth Engine Platform", 《 IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING》, vol. 14, pages 5773, XP011863242, DOI: 10.1109/JSTARS.2021.3085411 * |
| 刘浩: "卫星观测火箭尾喷焰红外动态场景生成研究", 《红外与激光工程》, vol. 50, no. 7, pages 124 - 133 * |
| 史超: "基于IDL的高光谱数据异常像元快速检测与修复", 《地理空间信息》, vol. 18, no. 8, pages 50 - 55 * |
| 叶昱含: "基于MODIS的中红外全球背景辐射数据生成方法研究", 《中国优秀硕士学位论文全文数据库基础科学辑》, pages 9 - 19 * |
| 曹鑫;陈学泓;张委伟;廖安平;陈利军;陈志刚;陈晋;: "全球30m空间分辨率耕地遥感制图研究", 中国科学:地球科学, no. 11, pages 18 - 27 * |
| 梁顺林;白瑞;陈晓娜;程洁;范闻捷;何涛;贾坤;江波;蒋玲梅;焦子锑;刘元波;倪文俭;邱凤;宋柳霖;孙林;唐伯惠;闻建光;吴桂平;谢东辉;姚云军;袁文平;张永光;张玉珍;张云腾;张晓通;赵天杰;赵祥;: "2019年中国陆表定量遥感发展综述", 遥感学报, no. 06, pages 13 - 66 * |
| 王东广: "结合纹理信息的高分辨率遥感图像变化检测方法", 《国土资源遥感》, no. 4, pages 76 - 81 * |
| 蔡禾: "Kamers-Kronig关系在反射太赫兹时域光谱测量中的应用", 《光谱学与光谱分析》, vol. 41, no. 1, pages 105 - 110 * |
| 黄宇: "自适应多点定标非均匀性校正算法", 《红外技术》, vol. 42, no. 7, pages 637 - 643 * |
| 龙耀威: "镀膜型光谱成像数据提取与作物叶绿素分布探测研究", 《光谱学与光谱分析》, vol. 40, no. 5, pages 1581 - 1587 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wan | MODIS land surface temperature products users’ guide | |
| TW200813887A (en) | Method and system for generating an image-textured digital surface model (DSM) for a geographical area of interest | |
| CN114758337B (en) | Semantic instance reconstruction method, device, equipment and medium | |
| Demetrescu et al. | Digital replica of cultural landscapes: An experimental reality-based workflow to create realistic, interactive open world experiences | |
| CN115358486B (en) | Port freight volume prediction method, system and application based on stereoscopic satellite images | |
| James et al. | An innovative photogrammetry color segmentation based technique as an alternative approach to 3D scanning for reverse engineering design | |
| CN114387289B (en) | Semantic segmentation method and device for three-dimensional point cloud of power transmission and distribution overhead line | |
| Aati et al. | Comparative study of photogrammetry software in industrial field | |
| Zhang et al. | An improved multi‐task pointwise network for segmentation of building roofs in airborne laser scanning point clouds | |
| CN114565768A (en) | Image segmentation method and device | |
| Syrris et al. | Mosaicking Copernicus Sentinel-1 data at global scale | |
| CN116912143A (en) | Image synthesis method, apparatus and computer readable medium | |
| CN114034651A (en) | Method and device for generating global earth surface spectrum basic data | |
| KR102390433B1 (en) | Global convergence video production device and production method | |
| CN115496654A (en) | Image super-resolution reconstruction method, device and medium based on self-attention mechanism | |
| CN115170962A (en) | Forest land change detection method based on unmanned aerial vehicle and satellite film | |
| Di Staso et al. | Heterogeneous-Resolution and multi-source terrain builder for CesiumJS WebGL virtual globe | |
| Maiwald et al. | Semi-automatic generation of historical urban 3D models at a larger scale using Structure-from-Motion, Neural Rendering and historical maps | |
| Koutaki et al. | Ishigaki retrieval system using 3D shape matching and combinatorial optimization | |
| Erwin et al. | The national geospatial digital archives—Collection development: Lessons learned | |
| CN116883575B (en) | Building group rendering method, device, computer equipment and storage medium | |
| THIRUCHELVAM et al. | Crime scene reconstruction based on a suitable software: A comparison study | |
| Collaro et al. | Detection of landscape features with visible and thermal imaging at the Castle of Puerta Arenas | |
| Wiser | Managing data for comprehensive environmental projects: case studies and recommendations | |
| CN117853664B (en) | Three-dimensional face reconstruction method based on double-branch feature fusion |
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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Zhu Xijuan Inventor after: Wu Jie Inventor after: Liu Dong Inventor after: Wu Kaifeng Inventor before: Zhu Xijuan Inventor before: Wu Jie Inventor before: Chen Shanjing Inventor before: Liu Dong Inventor before: Wu Kaifeng |
|
| CB03 | Change of inventor or designer information | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20220211 |
|
| RJ01 | Rejection of invention patent application after publication |