CN114494864A - Method for extracting phenological information of snow based on remote sensing data - Google Patents
Method for extracting phenological information of snow based on remote sensing data Download PDFInfo
- Publication number
- CN114494864A CN114494864A CN202210049179.4A CN202210049179A CN114494864A CN 114494864 A CN114494864 A CN 114494864A CN 202210049179 A CN202210049179 A CN 202210049179A CN 114494864 A CN114494864 A CN 114494864A
- Authority
- CN
- China
- Prior art keywords
- snow
- remote sensing
- sensing image
- image data
- data set
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Image Processing (AREA)
Abstract
The invention relates to the technical field of remote sensing data, in particular to a snow phenology information extraction method based on remote sensing data, which comprises the following steps: acquiring an accumulated snow remote sensing image data set of a target area, wherein the accumulated snow remote sensing image data set comprises accumulated snow remote sensing image data of a plurality of days; acquiring an accumulated snow coverage rate curve corresponding to the accumulated snow remote sensing image data set according to the accumulated snow remote sensing image data set, acquiring an accumulated snow coverage rate corresponding to each point in the accumulated snow coverage rate curve, and acquiring an initial snow day position point and a final snow day position point in the accumulated snow coverage rate curve according to the accumulated snow coverage rate corresponding to each point and a preset error interval; acquiring snow remote sensing image data in a target time period of the target area according to the initial snow day position point and the final snow day position point, and filling the snow remote sensing image data in the target time period of the target area to construct a target snow remote sensing image data set.
Description
Technical Field
The invention relates to the technical field of remote sensing data, in particular to a method, a device, equipment and a storage medium for extracting snow phenological information based on remote sensing data.
Background
The climate of snow is an important sign of global climate system change, and the close relationship between snow and the terrestrial ecosystem determines that a series of changes of the climate of snow can have significant influence on the terrestrial ecosystem. For example, due to the advance of snow melting time and the reduction of snow coverage, most of the vegetation exhibits phenomena such as extension of growing season, advance of flowering phase, and the like: the change of the accumulated snow causes the reduction of the effective moisture of the soil and the increase of the temperature of the soil, and further causes the obvious change of the vegetation community composition and species diversity in a cold area; as snow cover continues to decrease, vegetation productivity declines and carbon absorption capacity also tends to decline. On the animal side, snow melting times are advanced and temperatures are raised, leading to a change in the life cycle of a large number of invertebrates, such as a reduction in hibernation; the reduction of the species of the Nippon anthozoan is caused by the advancing and shortening of the flowering period of the plants; some invertebrates such as spiders have obvious phenotypic variation; vertebrates also respond significantly to snow changes, such as changes in the food chain that result in changes in the biological cycle of some animals and a decrease in the number of species that are present in some animals. Therefore, accurately acquiring the snow climate information is extremely important for researches on global changes, snow hydrology, biodiversity and the like.
The remote sensing technology is a key means for extracting the information of the phenological signs of the accumulated snow. The accumulated snow coverage curve obtained by utilizing the remote sensing data can well identify the initial snow day and the final snow day, and further calculate the accumulated snow duration. However, under the influence of adverse factors such as cloud blocking, invalid observed values, sensor faults and track deviation, an accumulated snow coverage curve acquired by a remote sensing technology is often discontinuous, and abnormal high values and abnormal low values are easy to occur, so that accumulated snow phenological information is extracted wrongly, and the application of remote sensing data in accumulated snow phenological information extraction is greatly limited.
Disclosure of Invention
Based on the above, the present invention provides a method, an apparatus, a device, and a storage medium for extracting snow phenology information based on remote sensing data, wherein a high-precision snow remote sensing image dataset is constructed by generating a snow coverage curve and performing filling processing on the snow remote sensing image dataset based on the snow remote sensing image dataset, so as to meet the requirement of snow phenology extraction on the time-space continuity of input snow coverage data and improve the accuracy of snow phenology information extraction.
In a first aspect, an embodiment of the present application provides a method for extracting snow phenological information based on remote sensing data, including the following steps:
acquiring an accumulated snow remote sensing image data set of a target area, wherein the accumulated snow remote sensing image data set comprises accumulated snow remote sensing image data of a plurality of days;
acquiring an accumulated snow coverage rate curve corresponding to the accumulated snow remote sensing image data set according to the accumulated snow remote sensing image data set, wherein the accumulated snow coverage rate curve comprises an accumulated snow coverage rate curve and an ablated snow coverage rate curve;
acquiring the snow cover rate corresponding to each point in the snow cover rate curve, and acquiring an initial snow day position point and a final snow day position point in the snow cover rate curve according to the snow cover rate corresponding to each point and a preset error interval;
acquiring snow remote sensing image data in a target time period of the target area according to the initial snow day position point and the final snow day position point, and filling the snow remote sensing image data in the target time period of the target area to construct a target snow remote sensing image data set;
in response to the instruction is drawed to snow phenology information, instruction is drawed to snow phenology information includes the regional snow remote sensing image dataset that awaits measuring, according to the regional snow remote sensing image dataset that awaits measuring acquires the regional corresponding target snow remote sensing image dataset that awaits measuring, according to the regional corresponding target snow remote sensing image dataset that awaits measuring acquires the regional snow remote sensing image dataset that awaits measuring corresponds first snow day position point and end snow day position point, as snow phenology information.
In a second aspect, an embodiment of the present application provides an apparatus for extracting snow phenology information based on remote sensing data, including:
the snow remote sensing image data acquisition module is used for acquiring a snow remote sensing image data set of a target area, wherein the snow remote sensing image data set comprises snow remote sensing image data of a plurality of days;
the curve acquisition module is used for acquiring an accumulated snow coverage rate curve corresponding to the accumulated snow remote sensing image data set according to the accumulated snow remote sensing image data set, wherein the accumulated snow coverage rate curve comprises an accumulated snow coverage rate curve in an accumulation period and an ablated snow coverage rate curve;
the position point acquisition module is used for acquiring the snow cover rate corresponding to each point in the snow cover rate curve and acquiring a first snow day position point and a final snow day position point in the snow cover rate curve according to the snow cover rate corresponding to each point and a preset error interval;
the data set construction module is used for acquiring snow remote sensing image data in a target time period of the target area according to the initial snow day position point and the final snow day position point, and filling the snow remote sensing image data in the target time period of the target area to construct a target snow remote sensing image data set;
snow phenology information draws module for respond to snow phenology information and draw the instruction, snow phenology information draws the instruction and includes the regional snow remote sensing image data set that awaits measuring, according to the regional snow remote sensing image data set that awaits measuring acquires the regional corresponding target snow remote sensing image data set that awaits measuring, according to the regional corresponding target snow remote sensing image data set that awaits measuring acquires the regional snow remote sensing image data set that awaits measuring corresponds first snow day position point and end snow day position point as snow phenology information.
In a third aspect, an embodiment of the present application provides a computer device, including: a processor, a memory, and a computer program stored on the memory and executable on the processor; the computer program, when executed by the processor, implements the steps of the method for extracting snow climate information based on remote sensing data according to the first aspect.
In a fourth aspect, the present application provides a storage medium, where a computer program is stored, and the computer program, when executed by a processor, implements the steps of the method for extracting snow climate information based on remote sensing data according to the first aspect.
In the embodiment of the application, a method, a device, equipment and a storage medium for extracting snow phenology information based on remote sensing data are provided, a high-precision snow remote sensing image data set is constructed by generating a snow coverage curve and filling the snow remote sensing image data set based on a snow remote sensing image data set, the requirement of snow phenology extraction on time-space continuity of input snow coverage data is met, and accuracy of snow phenology information extraction is improved.
For a better understanding and practice, the invention is described in detail below with reference to the accompanying drawings.
Drawings
Fig. 1 is a schematic flowchart of a method for extracting snow climate information based on remote sensing data according to an embodiment of the present application;
fig. 2 is a schematic flowchart of S2 in the method for extracting snow climate information based on remote sensing data according to an embodiment of the present application;
FIG. 3 is a schematic flow chart illustrating a method for extracting snow climate information based on remote sensing data according to another embodiment of the present application;
fig. 4 is a schematic flowchart of S4 in the method for extracting snow climate information based on remote sensing data according to an embodiment of the present application;
fig. 5 is a schematic structural diagram of an apparatus for extracting snow climate information based on remote sensing data according to an embodiment of the present application;
fig. 6 is a schematic structural diagram of a computer device according to an embodiment of the present application.
Detailed Description
Reference will now be made in detail to the exemplary embodiments, examples of which are illustrated in the accompanying drawings. When the following description refers to the accompanying drawings, like numbers in different drawings represent the same or similar elements unless otherwise indicated. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the present application. Rather, they are merely examples of apparatus and methods consistent with certain aspects of the present application, as detailed in the appended claims.
The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the application. As used in this application and the appended claims, the singular forms "a", "an", and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It should also be understood that the term "and/or" as used herein refers to and encompasses any and all possible combinations of one or more of the associated listed items.
It is to be understood that although the terms first, second, third, etc. may be used herein to describe various information, such information should not be limited to these terms. These terms are only used to distinguish one type of information from another. For example, first information may also be referred to as second information, and similarly, second information may also be referred to as first information, without departing from the scope of the present application. The word "if/if" as used herein may be interpreted as "at … …" or "when … …" or "in response to a determination", depending on the context.
Referring to fig. 1, fig. 1 is a schematic flow chart of a method for extracting snow climate information based on remote sensing data according to an embodiment of the present application, where the method includes the following steps:
s1: acquiring an accumulated snow remote sensing image data set of a target area, wherein the accumulated snow remote sensing image data set comprises accumulated snow remote sensing image data of a plurality of days.
The main execution body of the method for extracting the snow phenology information based on the remote sensing data is extraction equipment (hereinafter referred to as extraction equipment for short) of the method for extracting the snow phenology information based on the remote sensing data, and in an optional embodiment, the extraction equipment can be one computer device, a server or a server cluster formed by combining a plurality of computer devices.
In this embodiment, the control device may establish a data connection with a preset network database, and obtain an accumulated snow remote sensing image dataset of a target area from the network database, where the accumulated snow remote sensing image dataset includes accumulated snow remote sensing image data for a plurality of days, and the pixel types of the accumulated snow remote sensing image data include an accumulated snow pixel, a non-accumulated snow pixel, and an unoccupied pixel.
S2: and acquiring an accumulated snow coverage rate curve corresponding to the accumulated snow remote sensing image data set according to the accumulated snow remote sensing image data set, wherein the accumulated snow coverage rate curve comprises an accumulated snow coverage rate curve and an ablated snow coverage rate curve.
In this embodiment, the extraction device performs time-series arrangement on the snow remote sensing image data in the snow remote sensing image data set according to hydrologic years, and for the northern hemisphere, the hydrologic years refer to 9 months to 8 months of the next year, where 9 months to 2 months of the next year are snow accumulation periods, and 3 months to 8 months are snow ablation periods. Acquiring an accumulated snow coverage rate curve corresponding to the accumulated snow remote sensing image data set according to the sorted accumulated snow remote sensing image data set, dividing the accumulated snow coverage rate curve according to hydrologic years, and acquiring an accumulated snow coverage rate curve and an ablated snow coverage rate curve.
Referring to fig. 2, fig. 2 is a schematic flowchart of S2 in the method for extracting snow climate information based on remote sensing data according to an embodiment of the present application, including step S201, which is as follows:
s201: and extracting the NDSI value of each pixel of the accumulated snow remote sensing image data in the accumulated snow remote sensing image data set, acquiring the accumulated snow coverage rate corresponding to each pixel according to the NDSI value and a preset accumulated snow coverage rate calculation algorithm, and constructing an accumulated snow coverage rate curve according to the accumulated snow coverage rate corresponding to each pixel.
The NDSI is a normalized snow cover index, and the covered part of the accumulated snow in the remote sensing image is highlighted by utilizing the combination of visible light and short wave infrared wave bands.
In this embodiment, the extraction device extracts an NDSI value of each pixel of the snow remote sensing image data in the snow remote sensing image data set, and obtains an snow coverage rate corresponding to each pixel according to the NDSI value and a preset snow coverage rate calculation algorithm, where the snow coverage rate calculation algorithm is:
FSC＝-0.01+(1.45*NDSI)
wherein FSC is the coverage rate of the accumulated snow, and NDSI is the value of NDSI.
And combining according to the snow coverage rate corresponding to each pixel and a time sequence to construct a snow coverage rate curve.
S3: and acquiring the snow cover rate corresponding to each point in the snow cover rate curve, and acquiring an initial snow day position point and a final snow day position point in the snow cover rate curve according to the snow cover rate corresponding to each point and a preset error interval.
In this embodiment, the extraction device obtains the snow coverage corresponding to each point in the snow coverage curve, and obtains the first snow day position point and the final snow day position point in the snow coverage curve according to the snow coverage corresponding to each point and a preset error interval, as follows:
and acquiring a point of the accumulation period snow coverage curve where the snow coverage is greater than a preset snow coverage on an early snow day as a predicted early snow day position point, wherein in an optional embodiment, the snow coverage on the early snow day may be set to 0.
According to the predicted first snow day position point and the error interval, when the value of the predicted first snow day position point is located in the error interval, taking the predicted first snow day position point as a first snow day position point; when the value of the predicted first-snow-day position point is larger than the maximum value of the error interval, acquiring a position point corresponding to the previous day of the predicted first-snow-day position point as a first-snow-day position point;
and acquiring a point, in the ablation stage snow coverage rate curve, where the snow coverage rate is equal to a preset snow coverage rate of the final snow day as a predicted position point of the final snow day, wherein in an alternative embodiment, the snow coverage rate of the final snow day may be set to 0.
Obtaining values of position points corresponding to one day before and after the predicted final snow day position point, and when the value of the position point corresponding to one day before the predicted final snow day position point is larger than the final snow day snow coverage rate and is positioned in the error interval, and the value of the position point corresponding to one day after the predicted final snow day position point is equal to the final snow day snow coverage rate, taking the predicted final snow day position point as the final snow day position point; and when the value of the position point corresponding to the day before the predicted final snow day position point is greater than the final snow day snow accumulation coverage rate and is greater than the maximum value of the error interval, and the value of the position point corresponding to the day after the predicted final snow day position point is equal to the final snow day snow accumulation coverage rate, taking the position point corresponding to the day after the predicted final snow day position point as the final snow day position point.
Referring to fig. 3, fig. 3 is a schematic flowchart of a method for extracting snow climate information based on remote sensing data according to another embodiment of the present application, including step S6, where before step S3, the step S6 specifically includes the following steps:
s6: obtaining MOD10A1F data of the target area and the snow coverage rate corresponding to each pixel in the MOD10A1F data, constructing a snow coverage rate curve corresponding to the MOD10A1F data according to the snow coverage rate corresponding to each pixel in the MOD10A1F data, obtaining a standard deviation and an average difference corresponding to the MOD10A1F data according to the snow coverage rate curve corresponding to the MOD10A1F data, and establishing an error interval according to the standard deviation and the average difference.
The MOD10A1F data are complete daily NDSI data with a global 500 m resolution space-time over the period of 2001 to 2020.
In this embodiment, the extraction device may obtain MOD10A1F data of a target area from a database, obtain the snow coverage corresponding to each pixel in the MOD10A1F data according to the MOD10A1F data, construct a snow coverage curve corresponding to the MOD10A1F data according to the snow coverage corresponding to each pixel in the MOD10A1F data, calculate a standard deviation and an average deviation of the snow coverage curve corresponding to the MOD10A1F data, and establish an error interval by taking a value of an origin of the snow coverage curve corresponding to the MOD10A1F data as a center and the standard deviation ± 2 times as an error range according to the standard deviation and the average deviation.
S4: acquiring snow remote sensing image data in a target time period of the target area according to the initial snow day position point and the final snow day position point, and filling the snow remote sensing image data in the target time period of the target area to construct a target snow remote sensing image data set.
In this embodiment, the extraction device obtains snow remote sensing image data in a target time period of the target area by using the initial snow day position point as a starting point and the final snow day position point as an end point, and constructs initial snow remote sensing image data;
in order to improve the space integrity of the preliminary snow remote sensing image data, on the basis of judgment and filling of the prior knowledge, according to the space correlation, filling processing is carried out on vacant pixels in the snow remote sensing image data in the target time period of the target area, and a target snow remote sensing image data set is constructed so as to improve the space integrity of the snow phenological data.
Referring to fig. 4, fig. 4 is a schematic flow chart of S4 in the method for extracting snow climate information based on remote sensing data according to an embodiment of the present application, including steps S401 to S402, which are as follows:
s401: and acquiring elevation data and a snow phenological value corresponding to each pixel in each snow remote sensing image data in the target snow remote sensing image data set.
The pixel of snow remote sensing image data has snow phenological value, snow phenological value is used for embodying first snow day position point and final snow day position point in the snow phenological information, the Elevation data does the Elevation data is Digital Elevation Model (DEM), and its Digital Elevation information that represents the ground topography.
In this embodiment, the extraction device acquires elevation image data close to the spatial resolution according to the spatial resolution of each snow remote sensing image data in the target snow remote sensing image data set, resamples the elevation image data to prevent spatial overlapping, and spatially matches the target snow remote sensing image data set. And acquiring elevation data corresponding to each pixel in each accumulated snow remote sensing image data in the target accumulated snow remote sensing image data set.
S402: and filling vacant pixels in the snow remote sensing image data according to the pixel types, the elevation data and the snow phenological values corresponding to the pixels in the snow remote sensing image data in the target snow remote sensing image data set, and acquiring a filled snow remote sensing image data set as a target snow remote sensing image data set.
In this embodiment, the extraction device fills the vacant pixels in each snow remote sensing image data according to the pixel types, the elevation data and the snow phenological values corresponding to each pixel in each snow remote sensing image data in the target snow remote sensing image data set, and obtains the filled snow remote sensing image data set as the target snow remote sensing image data set.
In an optional embodiment, the extracting device obtains neighboring pixels of the vacant pixel and a category of the neighboring pixels;
when all the adjacent pixels of the vacant pixel are the snow covered pixels, acquiring the adjacent pixel with the minimum difference value with the elevation data of the vacant pixel from the adjacent pixels according to the difference value of the elevation data of the vacant pixel and the adjacent pixels, taking the adjacent pixel as a target pixel, modifying the vacant pixel into the snow covered pixel, and filling the snow covered phenological value of the target pixel into the vacant pixel;
when the adjacent pixels of the vacant pixels comprise the snow covered pixels and the non-snow covered pixels, acquiring the adjacent pixels with the minimum difference value with the elevation data of the vacant pixels from the adjacent pixels as first target pixels according to the values of the elevation data of the vacant pixels and the adjacent pixels and when the value of the elevation data of the vacant pixels is larger than the value of the elevation data of any snow covered pixel in the adjacent pixels, modifying the vacant pixels into snow covered pixels, and filling the snow covered phenological value of the first target pixels into the vacant pixels;
when the value of the elevation data of the vacant pixel is smaller than the value of any one of the adjacent pixels which is not covered with snow, modifying the vacant pixel into the pixel which is not covered with snow;
when the value of the elevation data of the vacant pixel is smaller than the value of the elevation data of any one of the adjacent pixels and larger than the value of any one of the adjacent pixels without snow, calculating the mean value of the elevation data corresponding to the snow pixel and the non-snow pixel in the adjacent pixels, comparing the value of the elevation data of the vacant pixel with the mean value of the elevation data of the snow pixel in the adjacent pixels and the mean value of the elevation data corresponding to the non-snow pixel in the adjacent pixels respectively, acquiring a pixel with a value close to the value of the elevation data of the vacant pixel as a second target pixel, when the second target pixel is the non-snow pixel, non-modifying the vacant pixel into the snow pixel, when the second target pixel is the snow pixel, modifying the vacant pixel into the snow pixel, and acquiring the adjacent pixel with the minimum difference value of the elevation data of the vacant pixel from the adjacent pixels, and filling the snow phenological value of the third target pixel into the vacant pixel as a third target pixel.
S5: in response to the instruction is drawed to snow phenology information, instruction is drawed to snow phenology information includes the regional snow remote sensing image dataset that awaits measuring, according to the regional snow remote sensing image dataset that awaits measuring acquires the regional corresponding target snow remote sensing image dataset that awaits measuring, according to the regional corresponding target snow remote sensing image dataset that awaits measuring acquires the regional snow remote sensing image dataset that awaits measuring corresponds first snow day position point and end snow day position point, as snow phenology information.
The snow climate information is used for embodying an initial snow day, a final snow day and snow duration, the snow climate information extraction instruction is sent by a user, and the extraction equipment receives the snow climate information.
In this embodiment, the extraction device acquires the snow phenological information extraction instruction sent by the user, and responds to acquire the snow remote sensing image dataset of the area to be measured. The extraction equipment acquires the target snow remote sensing image data set corresponding to the area to be detected according to the snow remote sensing image data set of the area to be detected, acquires the initial snow day position point and the final snow day position point corresponding to the snow remote sensing image data set of the area to be detected according to the target snow remote sensing image data set corresponding to the area to be detected, and the initial snow day position point and the final snow day position point serve as snow phenological information and return to a display interface of the extraction equipment to be displayed and labeled.
Referring to fig. 5, fig. 5 is a schematic structural diagram of an apparatus for extracting snow climate information based on remote sensing data according to an embodiment of the present application, where the apparatus may implement all or a part of the apparatus for extracting snow climate information based on remote sensing data through software, hardware, or a combination of the two, and the apparatus 5 includes:
a data set obtaining module 51, configured to obtain an accumulated snow remote sensing image data set of a target area, where the accumulated snow remote sensing image data set includes accumulated snow remote sensing image data for a plurality of days;
the curve obtaining module 52 is configured to obtain an accumulated snow coverage curve corresponding to the accumulated snow remote sensing image data set according to the accumulated snow remote sensing image data set, where the accumulated snow coverage curve includes an accumulated snow coverage curve in an accumulation period and an ablated snow coverage curve in an ablation period;
a position point obtaining module 53, configured to obtain an accumulated snow coverage rate corresponding to each point in the accumulated snow coverage rate curve, and obtain a first snow day position point and a last snow day position point in the accumulated snow coverage rate curve according to the accumulated snow coverage rate corresponding to each point and a preset error interval;
the data set construction module 54 is configured to obtain snow remote sensing image data in the target time period of the target area according to the first snow day position point and the final snow day position point, and perform filling processing on the snow remote sensing image data in the target time period of the target area to construct a target snow remote sensing image data set;
snow phenology information draws module 55 for respond to snow phenology information and draw the instruction, snow phenology information draws the instruction and includes the regional snow remote sensing image dataset that awaits measuring, according to the regional snow remote sensing image dataset that awaits measuring acquires the regional corresponding target snow remote sensing image dataset that awaits measuring, according to the regional corresponding target snow remote sensing image dataset that awaits measuring acquires the regional snow remote sensing image dataset that awaits measuring corresponds first snow day position point and end snow day position point, as snow phenology information.
In the embodiment of the application, an accumulated snow remote sensing image data set of a target area is obtained through a data set obtaining module, wherein the accumulated snow remote sensing image data set comprises accumulated snow remote sensing image data of a plurality of days; acquiring an accumulated snow coverage rate curve corresponding to the accumulated snow remote sensing image data set through a curve acquisition module according to the accumulated snow remote sensing image data set, wherein the accumulated snow coverage rate curve comprises an accumulated snow coverage rate curve in an accumulation period and an ablated snow coverage rate curve; acquiring an accumulated snow coverage rate corresponding to each point in the accumulated snow coverage rate curve through a position point acquisition module, and acquiring an initial snow day position point and a final snow day position point in the accumulated snow coverage rate curve according to the accumulated snow coverage rate corresponding to each point and a preset error interval; acquiring snow remote sensing image data in a target time period of the target area according to the initial snow day position point and the final snow day position point through a data set construction module, and filling the snow remote sensing image data in the target time period of the target area to construct a target snow remote sensing image data set; through snow phenological information extraction module, respond to snow phenological information and draw the instruction, snow phenological information draws the instruction and includes the regional snow remote sensing image dataset that awaits measuring, according to the regional snow remote sensing image dataset that awaits measuring acquires the regional corresponding target snow remote sensing image dataset that awaits measuring, according to the regional corresponding target snow remote sensing image dataset that awaits measuring acquires the regional snow remote sensing image dataset that awaits measuring corresponds first snow day position point and end snow day position point as snow phenological information. Based on the accumulated snow remote sensing image data set, the accumulated snow coverage rate curve is generated, the accumulated snow remote sensing image data set is filled, so that the high-precision accumulated snow remote sensing image data set is constructed, the requirement of accumulated snow phenological extraction on the time-space continuity of the input accumulated snow coverage rate data is met, and the accuracy of accumulated snow phenological information extraction is improved.
Referring to fig. 6, fig. 6 is a schematic structural diagram of a computer device according to an embodiment of the present application, where the computer device 6 includes: a processor 61, a memory 62 and a computer program 63 stored on the memory 62 and executable on the processor 61; the computer device may store a plurality of instructions, where the instructions are suitable for being loaded by the processor 61 and executing the method steps in the embodiments shown in fig. 1 to fig. 4, and a specific execution process may refer to specific descriptions of the embodiments shown in fig. 1 to fig. 4, which is not described herein again.
The Memory 62 may include a Random Access Memory (RAM) 62, and may also include a Read-Only Memory (Read-Only Memory) 62. Optionally, the memory 62 includes a non-transitory computer-readable medium. The memory 62 may be used to store instructions, programs, code sets or instruction sets. The memory 62 may include a program storage area and a data storage area, wherein the program storage area may store instructions for implementing an operating system, instructions for at least one function (such as touch instructions, etc.), instructions for implementing the various method embodiments described above, and the like; the storage data area may store data and the like referred to in the above respective method embodiments. The memory 62 may optionally be at least one memory device located remotely from the aforementioned processor 61.
An embodiment of the present application further provides a storage medium, where the storage medium may store multiple instructions, and the instructions are suitable for being loaded by a processor and being executed by the method steps in the embodiments shown in fig. 1 to fig. 4, and a specific execution process may refer to specific descriptions of the embodiments shown in fig. 1 to fig. 4, which is not described herein again.
It should be clear to those skilled in the art that, for convenience and simplicity of description, the foregoing division of the functional units and modules is only used for illustration, and in practical applications, the above function distribution may be performed by different functional units and modules as needed, that is, the internal structure of the apparatus may be divided into different functional units or modules to perform all or part of the above described functions. Each functional unit and module in the embodiments may be integrated in one processing unit, or each unit may exist alone physically, or two or more units are integrated in one unit, and the integrated unit may be implemented in a form of hardware, or in a form of software functional unit. In addition, specific names of the functional units and modules are only for convenience of distinguishing from each other, and are not used for limiting the protection scope of the present application. The specific working processes of the units and modules in the system may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
In the above embodiments, the descriptions of the respective embodiments have respective emphasis, and reference may be made to the related descriptions of other embodiments for parts that are not described or illustrated in a certain embodiment.
Those of ordinary skill in the art will appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
In the embodiments provided in the present invention, it should be understood that the disclosed apparatus/terminal device and method may be implemented in other ways. For example, the above-described embodiments of the apparatus/terminal device are merely illustrative, and for example, the division of the modules or units is only one logical division, and there may be other divisions when actually implemented, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present invention may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit. The integrated unit can be realized in a form of hardware, and can also be realized in a form of a software functional unit.
The integrated modules/units, if implemented in the form of software functional units and sold or used as separate products, may be stored in a computer readable storage medium. Based on such understanding, all or part of the flow of the method according to the embodiments of the present invention may also be implemented by a computer program, which may be stored in a computer-readable storage medium, and when the computer program is executed by a processor, the steps of the method embodiments may be implemented. Wherein the computer program comprises computer program code, which may be in the form of source code, object code, an executable file or some intermediate form, etc.
The present invention is not limited to the above-described embodiments, and various modifications and variations of the present invention are intended to be included within the scope of the claims and the equivalent technology of the present invention if they do not depart from the spirit and scope of the present invention.
Claims (8)
1. A snow phenology information extraction method based on remote sensing data is characterized by comprising the following steps:
acquiring an accumulated snow remote sensing image data set of a target area, wherein the accumulated snow remote sensing image data set comprises accumulated snow remote sensing image data of a plurality of days;
acquiring an accumulated snow coverage rate curve corresponding to the accumulated snow remote sensing image data set according to the accumulated snow remote sensing image data set, wherein the accumulated snow coverage rate curve comprises an accumulated snow coverage rate curve in an accumulation period and an ablated snow coverage rate curve in an ablation period;
acquiring the snow cover rate corresponding to each point in the snow cover rate curve, and acquiring an initial snow day position point and a final snow day position point in the snow cover rate curve according to the snow cover rate corresponding to each point and a preset error interval;
acquiring snow remote sensing image data in a target time period of the target area according to the initial snow day position point and the final snow day position point, and filling the snow remote sensing image data in the target time period of the target area to construct a target snow remote sensing image data set;
in response to the instruction is drawed to snow phenology information, instruction is drawed to snow phenology information includes the regional snow remote sensing image dataset that awaits measuring, according to the regional snow remote sensing image dataset that awaits measuring acquires the regional corresponding target snow remote sensing image dataset that awaits measuring, according to the regional corresponding target snow remote sensing image dataset that awaits measuring acquires the regional snow remote sensing image dataset that awaits measuring corresponds first snow day position point and end snow day position point, as snow phenology information.
2. The method for extracting snow phenology information based on remote sensing data of claim 1, wherein the step of obtaining a snow coverage curve corresponding to the snow remote sensing image dataset according to the snow remote sensing image dataset comprises the steps of:
the method comprises the steps of extracting an NDSI value of each pixel of accumulated snow remote sensing image data in an accumulated snow remote sensing image data set, obtaining an accumulated snow coverage rate corresponding to each pixel according to the NDSI value and a preset accumulated snow coverage rate calculation algorithm, and constructing an accumulated snow coverage rate curve according to the accumulated snow coverage rate corresponding to each pixel, wherein the accumulated snow coverage rate calculation algorithm is as follows:
FSC＝-0.01+(1.45*NDSI)
wherein FSC is the coverage rate of the accumulated snow, and NDSI is the value of NDSI.
3. The method for extracting snow climate information based on remote sensing data according to claim 1, wherein before acquiring an initial snow day location point and a final snow day location point in the snow coverage curve according to the snow coverage corresponding to each point and a preset error interval, the method comprises:
obtaining MOD10A1F data of the target area and the snow coverage rate corresponding to each pixel in the MOD10A1F data, constructing a snow coverage rate curve corresponding to the MOD10A1F data according to the snow coverage rate corresponding to each pixel in the MOD10A1F data, obtaining a standard deviation and an average difference corresponding to the MOD10A1F data according to the snow coverage rate curve corresponding to the MOD10A1F data, and establishing an error interval according to the standard deviation and the average difference.
4. The method for extracting snow climate information based on remote sensing data according to claim 1, wherein: the pixel types of the snow accumulation remote sensing image data comprise snow accumulation pixels, non-snow accumulation pixels and vacant pixels.
5. The method for extracting snow phenology information based on remote sensing data according to claim 4, wherein the filling processing is performed on snow remote sensing image data in a target time period of the target area to construct a target snow remote sensing image data set, and the method comprises the following steps:
acquiring elevation data and snow phenological values corresponding to pixels in snow remote sensing image data in the target snow remote sensing image data set;
and filling vacant pixels in the snow remote sensing image data according to the pixel types, the elevation data and the snow phenological values corresponding to the pixels in the snow remote sensing image data in the target snow remote sensing image data set, and acquiring a filled snow remote sensing image data set as a target snow remote sensing image data set.
6. The utility model provides an extraction element of snow phenology information based on remote sensing data which characterized in that includes:
the snow remote sensing image data acquisition module is used for acquiring a snow remote sensing image data set of a target area, wherein the snow remote sensing image data set comprises snow remote sensing image data of a plurality of days;
the curve acquisition module is used for acquiring an accumulated snow coverage rate curve corresponding to the accumulated snow remote sensing image data set according to the accumulated snow remote sensing image data set, wherein the accumulated snow coverage rate curve comprises an accumulated snow coverage rate curve in an accumulation period and an ablated snow coverage rate curve;
the position point acquisition module is used for acquiring the snow cover rate corresponding to each point in the snow cover rate curve and acquiring a first snow day position point and a final snow day position point in the snow cover rate curve according to the snow cover rate corresponding to each point and a preset error interval;
the data set construction module is used for acquiring snow remote sensing image data in a target time period of the target area according to the initial snow day position point and the final snow day position point, and filling the snow remote sensing image data in the target time period of the target area to construct a target snow remote sensing image data set;
snow phenology information draws module for respond to snow phenology information and draw the instruction, snow phenology information draws the instruction and includes the regional snow remote sensing image data set that awaits measuring, according to the regional snow remote sensing image data set that awaits measuring acquires the regional corresponding target snow remote sensing image data set that awaits measuring, according to the regional corresponding target snow remote sensing image data set that awaits measuring acquires the regional snow remote sensing image data set that awaits measuring corresponds first snow day position point and end snow day position point as snow phenology information.
7. A computer device, comprising: a processor, a memory, and a computer program stored on the memory and executable on the processor; the computer program, when executed by the processor, implements the steps of the method for extracting snow climate information based on remote sensing data according to any of claims 1 to 5.
8. A storage medium, characterized by: the storage medium stores a computer program which, when executed by a processor, implements the steps of the method for remote sensing data-based extraction of snow climate information according to any of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210049179.4A CN114494864B (en) | 2022-01-17 | 2022-01-17 | Snow phenological information extraction method based on remote sensing data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210049179.4A CN114494864B (en) | 2022-01-17 | 2022-01-17 | Snow phenological information extraction method based on remote sensing data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114494864A true CN114494864A (en) | 2022-05-13 |
| CN114494864B CN114494864B (en) | 2022-12-27 |
Family
ID=81511042
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210049179.4A Active CN114494864B (en) | 2022-01-17 | 2022-01-17 | Snow phenological information extraction method based on remote sensing data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114494864B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115203594A (en) * | 2022-09-16 | 2022-10-18 | 成都国星宇航科技股份有限公司 | Multi-temporal remote sensing data display method, device, equipment and medium |
| CN116127273A (en) * | 2023-01-18 | 2023-05-16 | 中国科学院地理科学与资源研究所 | Snow index acquisition method, device, storage medium and equipment |
| CN117173578A (en) * | 2023-11-01 | 2023-12-05 | 长江水利委员会长江科学院 | Snow line elevation detection method and device, computer equipment and storage medium |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101799561A (en) * | 2010-02-05 | 2010-08-11 | 民政部国家减灾中心 | Snow disaster remote sensing monitoring simulation evaluation method based on disaster reduction small satellite |
| CN103984862A (en) * | 2014-05-15 | 2014-08-13 | 中国科学院遥感与数字地球研究所 | Multielement remote sensing information coordinated snow cover parameter inversion method |
| CN104021283A (en) * | 2014-05-23 | 2014-09-03 | 清华大学 | Prediction method and device of day runoff volume of snowmelt period |
| CN106780433A (en) * | 2016-11-16 | 2017-05-31 | 清华大学 | Accumulated snow phenology information fusion method and system based on multi- source Remote Sensing Data data |
| CN108765406A (en) * | 2018-05-31 | 2018-11-06 | 电子科技大学 | A kind of snow mountain detection method based on infrared remote sensing image |
| CN112052589A (en) * | 2020-09-03 | 2020-12-08 | 兰州交通大学 | Method and device for estimating snow cover proportion based on MODIS day-by-day snow product |
| US20210110157A1 (en) * | 2019-10-15 | 2021-04-15 | A.U.G. Signals Ltd. | Crop phenology estimation and tracking with remote sensing imagery |
| CN113936213A (en) * | 2021-10-14 | 2022-01-14 | 张定祥 | Glacier and permanent snow macroscopic monitoring method based on MODIS data |
-
2022
- 2022-01-17 CN CN202210049179.4A patent/CN114494864B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101799561A (en) * | 2010-02-05 | 2010-08-11 | 民政部国家减灾中心 | Snow disaster remote sensing monitoring simulation evaluation method based on disaster reduction small satellite |
| CN103984862A (en) * | 2014-05-15 | 2014-08-13 | 中国科学院遥感与数字地球研究所 | Multielement remote sensing information coordinated snow cover parameter inversion method |
| CN104021283A (en) * | 2014-05-23 | 2014-09-03 | 清华大学 | Prediction method and device of day runoff volume of snowmelt period |
| CN106780433A (en) * | 2016-11-16 | 2017-05-31 | 清华大学 | Accumulated snow phenology information fusion method and system based on multi- source Remote Sensing Data data |
| CN108765406A (en) * | 2018-05-31 | 2018-11-06 | 电子科技大学 | A kind of snow mountain detection method based on infrared remote sensing image |
| US20210110157A1 (en) * | 2019-10-15 | 2021-04-15 | A.U.G. Signals Ltd. | Crop phenology estimation and tracking with remote sensing imagery |
| CN112052589A (en) * | 2020-09-03 | 2020-12-08 | 兰州交通大学 | Method and device for estimating snow cover proportion based on MODIS day-by-day snow product |
| CN113936213A (en) * | 2021-10-14 | 2022-01-14 | 张定祥 | Glacier and permanent snow macroscopic monitoring method based on MODIS data |
Non-Patent Citations (2)
| Title |
|---|
| HUI GUO ET AL.: "Assessing Snow Phenology and Its Environmental Driving Factors in Northeast China", 《REMOTE SENSING》 * |
| 李晨昊等: "2000～2017年蒙古高原积雪时空变化及其对气候响应研究", 《中国草地学报》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115203594A (en) * | 2022-09-16 | 2022-10-18 | 成都国星宇航科技股份有限公司 | Multi-temporal remote sensing data display method, device, equipment and medium |
| CN115203594B (en) * | 2022-09-16 | 2022-11-29 | 成都国星宇航科技股份有限公司 | Multi-temporal remote sensing data display method, device, equipment and medium |
| CN116127273A (en) * | 2023-01-18 | 2023-05-16 | 中国科学院地理科学与资源研究所 | Snow index acquisition method, device, storage medium and equipment |
| CN116127273B (en) * | 2023-01-18 | 2023-12-12 | 中国科学院地理科学与资源研究所 | Snow index acquisition method, device, storage medium and equipment |
| CN117173578A (en) * | 2023-11-01 | 2023-12-05 | 长江水利委员会长江科学院 | Snow line elevation detection method and device, computer equipment and storage medium |
| CN117173578B (en) * | 2023-11-01 | 2024-02-06 | 长江水利委员会长江科学院 | Snow line elevation detection method and device, computer equipment and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114494864B (en) | 2022-12-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114494864B (en) | Snow phenological information extraction method based on remote sensing data | |
| Onojeghuo et al. | Rice crop phenology mapping at high spatial and temporal resolution using downscaled MODIS time-series | |
| Schmidtlein et al. | The spectral variability hypothesis does not hold across landscapes | |
| Mack et al. | A semi-automated approach for the generation of a new land use and land cover product for Germany based on Landsat time-series and Lucas in-situ data | |
| Guild et al. | Detection of deforestation and land conversion in Rondonia, Brazil using change detection techniques | |
| JP4873545B2 (en) | Field management support method and system | |
| Fernández-Manso et al. | Forecast of NDVI in coniferous areas using temporal ARIMA analysis and climatic data at a regional scale | |
| Pagani et al. | A high-resolution, integrated system for rice yield forecasting at district level | |
| Pennec et al. | Tropical forest phenology in French Guiana from MODIS time series | |
| Cao et al. | Mechanisms, monitoring and modeling of shrub encroachment into grassland: a review | |
| Ishii et al. | Early detection of the invasive alien plant Solidago altissima in moist tall grassland using hyperspectral imagery | |
| White et al. | Accuracy assessment of the vegetation continuous field tree cover product using 3954 ground plots in the south‐western USA | |
| Guo et al. | Forest cover dynamics from Landsat time-series data over Yan’an city on the Loess Plateau during the Grain for Green Project | |
| CN114202691B (en) | Rice identification method, device and equipment based on remote sensing image | |
| Ren et al. | Enhanced spatiotemporal heterogeneity and the climatic and biotic controls of autumn phenology in northern grasslands | |
| CN114612896A (en) | Rice yield prediction method, device and equipment based on remote sensing image | |
| Ferreira et al. | Equivalent water thickness in savanna ecosystems: MODIS estimates based on ground and EO-1 Hyperion data | |
| Silveira et al. | Modelling aboveground biomass in forest remnants of the Brazilian Atlantic Forest using remote sensing, environmental and terrain-related data | |
| Zhang et al. | Land–use and land-cover change detection using dynamic time warping–based time series clustering method | |
| Buxton et al. | Quantitatively monitoring the resilience of patterned vegetation in the Sahel | |
| Yang et al. | Estimation of tropical forest biomass using Landsat TM imagery and permanent plot data in Xishuangbanna, China | |
| Xie | Combining CERES-Wheat model, Sentinel-2 data, and deep learning method for winter wheat yield estimation | |
| JP6509887B2 (en) | Image change detection | |
| CN116579521B (en) | Yield prediction time window determining method, device, equipment and readable storage medium | |
| Laskin et al. | Evaluating the level of agreement between human and time-lapse camera observations of understory plant phenology at multiple scales |
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 |