CN108230310B - Method for extracting non-fire spatio-temporal data based on semi-variogram - Google Patents

Method for extracting non-fire spatio-temporal data based on semi-variogram Download PDF

Info

Publication number
CN108230310B
CN108230310B CN201810003693.8A CN201810003693A CN108230310B CN 108230310 B CN108230310 B CN 108230310B CN 201810003693 A CN201810003693 A CN 201810003693A CN 108230310 B CN108230310 B CN 108230310B
Authority
CN
China
Prior art keywords
fire
data
pixels
pixel
buffer
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810003693.8A
Other languages
Chinese (zh)
Other versions
CN108230310A (en
Inventor
何彬彬
文崇波
全兴文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201810003693.8A priority Critical patent/CN108230310B/en
Publication of CN108230310A publication Critical patent/CN108230310A/en
Application granted granted Critical
Publication of CN108230310B publication Critical patent/CN108230310B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/62Analysis of geometric attributes of area, perimeter, diameter or volume
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10032Satellite or aerial image; Remote sensing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Probability & Statistics with Applications (AREA)
  • Operations Research (AREA)
  • Evolutionary Biology (AREA)
  • Algebra (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Geometry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Quality & Reliability (AREA)
  • Image Processing (AREA)

Abstract

The invention relates to a method for extracting non-fire spatio-temporal data based on a half-variation function, which comprises the steps of firstly determining the positions and the occurrence moments of all fire pixels in combustion area data, judging the vegetation types of all the fire pixels, then extracting the data of all the fire pixels and the surrounding non-fire pixel data with the same vegetation type as the fire pixels, carrying out spatial correlation analysis on the data by using the half-variation function, establishing a corresponding fire buffer zone by taking a variation process as a buffer distance, then extracting the data of the fire pixels at the relevant time phase, repeating the process of establishing the fire buffer zone, and finally selecting the data of the pixels which do not fall into all the buffer zones as the non-fire data. The method can eliminate subjective influence caused by human factors, eliminate differences caused by factors such as geographical positions, vegetation types and the like, eliminate accidental factors of fire occurrence, improve the contrast of non-fire data and be beneficial to improving the precision of a simulation result of a training model.

Description

Method for extracting non-fire spatio-temporal data based on semi-variogram
Technical Field
The invention particularly relates to a method for extracting non-fire spatio-temporal data based on a half-variogram.
Background
The forest fire is a very common and destructive natural disaster, the forest fire happens more than 20 ten thousand times every year in the world on average, and the burning area of the forest occupies more than 1 per thousand of the total area of the forest in the world. The forest fire happens more than 1 ten thousand times per year in China, hundreds of thousands of forests to millions of hectares are burnt, and the forest area occupies 5-8 per thousand of the whole country. Forest fires not only burn trees to directly reduce the forest area, but also seriously destroy the forest structure and the forest environment, so that the forest ecological system is out of balance, the forest biomass is reduced, and even people and livestock are killed.
The forest fire danger level prediction method is very meaningful and necessary work, can realize the forest fire danger level prediction with large range and high space-time resolution by a remote sensing technology, and has certain scientific guiding significance for the allocation of forest fire protection resources and the personnel movement. At present, the research on fire risk level prediction models is relatively deep, and the training models generally need to input climate parameters, terrain parameters and vegetation parameters, wherein the vegetation parameters can reflect the physical characteristics (such as type, coverage and the like) of vegetation. The training model requires comparison data, i.e. the input data to the model includes not only the data of a fire (fire data) but also the data of a non-fire (non-fire data). The quality of the input data has a significant impact on the final simulation effect of the training model.
At present, a method for acquiring fire data and non-fire data based on a remote sensing technology mainly comprises the following steps: firstly, carrying out fire burning area mapping according to a fire point extraction algorithm or directly utilizing related fire products, and then extracting data of fire pixel according to a mapping result or the fire products, thereby acquiring fire data; the obtained non-fire data usually takes a fire pixel as a circle center, the buffer radius is determined according to experience to further establish a buffer area, and data of the non-fire pixel is randomly extracted from an area outside the buffer area to serve as comparison data of the fire pixel. Although the method for extracting the non-fire data is simple in operation, the method has the following problems: firstly, the determination of the fire pixel buffer radius is divided according to experience, subjective influence caused by human factors exists, meanwhile, the size of the buffer radius changes along with the difference of a research area, a vegetation type and the like, that is, the buffer radius of each fire pixel is different, if the buffer radius is too small, space correlation exists between extracted non-fire data and corresponding fire data, so that the contrast is lacked, the simulation effect of a training model is influenced, and the prediction purpose cannot be well achieved; if the buffer radius is set too large, the area for non-fire data extraction will be too small. Secondly, the extraction method of non-fire data only considers spatial correlation and does not consider temporal correlation, such as: a certain picture element has no fire in a certain time phase, but a fire in the next time phase occurs, which indicates that the data of the picture element is close to the fire data, so that the non-fire data extracted by the method still lacks the contrast. In summary, how to accurately, conveniently and efficiently extract non-fire data becomes a technical problem to be solved urgently in the field.
Disclosure of Invention
In view of the above, the present invention aims to: aiming at the problems of poor contrast and the like of the existing non-fire data extraction method, the invention provides the method for extracting the non-fire spatio-temporal data by combining the fitting half-variation function value of the spherical model and the multi-temporal characteristics of the remote sensing image.
The technical scheme provided by the invention is as follows:
a method for extracting non-fire data based on a semi-variogram is characterized by comprising the following steps:
step 1: determining the position and the occurrence moment of a fire pixel in the combustion area data;
step 2: determining the vegetation type of the fire pixel;
and step 3: extracting data of fire pixels and surrounding non-fire pixel data with the same type as the vegetation type of the fire pixels;
and 4, step 4: performing spatial correlation analysis on the fire pixels by using a half-variation function, taking a variation range as a buffer distance of the fire pixels and establishing a corresponding fire buffer area;
and 5: extracting data of a relevant time phase of a fire pixel;
step 6: repeating the step 4 to establish a fire buffer zone of a time phase related to the fire pixel;
and 7: and (4) selecting pixels which do not fall into the fire buffer areas obtained in the step (4) and the step (6) and extracting data of the pixels to serve as non-fire data.
Furthermore, the variable range is obtained by fitting a spherical model to the half-variation function value.
Furthermore, the buffer area is specifically established by taking the fire pixel as the center of a circle and taking the variable range as the radius to form the buffer area.
Furthermore, the relevant time phases of the fire pixel in the invention are specifically as follows: and taking a single fire pixel as a research object, and taking the moment of fire occurrence of the pixel as a target time phase, wherein the adjacent time phases of the target time phase and the time phases corresponding to the target time phase in the rest years are related time phases.
Further, the combustion area data is specifically remote sensing image data.
The inventive concept is explained in detail below: in order to consider the spatial correlation between the fire pixel and the non-fire pixel, a buffer area is often required to be established for the fire pixel. And a plurality of fire pixels are often found in the fire distribution remote sensing image, so that buffer areas of all the fire pixels need to be obtained, and the comparison data of the fire pixels cannot fall into any one buffer area. Because the vegetation types have difference, the invention introduces the vegetation type factor in the process of establishing the buffer zone, and utilizes the spherical model to fit the variation trend of the half variation function value, takes the variation as the threshold distance of the spatial correlation, and establishes the buffer zone of the fire pixel; meanwhile, the time correlation of the fire pixel at different time phases is considered, and the extraction of non-fire pixel data close to the fire pixel data is avoided, so that the difference between the fire data and the non-fire data is improved, and the aim of improving the simulation precision of the fire risk grade prediction model is fulfilled.
Compared with the prior art, the invention has the beneficial effects that:
according to the invention, vegetation type factors are introduced, and the variation range of a half variation function is used as a buffer radius to establish a buffer area of a fire pixel, so that not only are subjective influences caused by human factors eliminated, but also differences caused by different geographic positions, vegetation types and the like are eliminated; meanwhile, the time relativity of the fire pixel at different time phases is considered, so that non-fire data with insufficient difference with fire data are eliminated; therefore, the method and the device can enhance the contrast between fire data and non-fire data, and further improve the simulation precision of the fire risk grade prediction model.
Drawings
FIG. 1 is a schematic view of a fire distribution image and fire pixel buffer radius in Michelle county, Yunnan province.
FIG. 2 is a graph showing the variation of the half-value of the variation function of a Normalized Difference Vegetation Index (NDVI).
FIG. 3 is a graph showing the variation of the half-value of the variation function of the water Content (FMC) of combustible.
FIG. 4 is a diagram illustrating the superposition of buffers for different index parameters.
FIG. 5 is a schematic diagram of the superposition of buffers for different fire pixels.
Fig. 6 is a diagram of a multi-temporal buffer overlay.
Detailed Description
The principles and features of the present invention will be further explained with reference to the drawings and specific embodiments of the specification:
the occurrence and development of fire are very complex processes, and many influencing factors are involved, so that different representative fire risk evaluation indexes are selected according to different time, space scales and research purposes in fire risk grade evaluation. At present, the fire risk indexes commonly used can be roughly divided into the aspects of terrain conditions, weather environment, vegetation factors and the like, wherein the vegetation factors are one of the main influence factors influencing the fire maturity of vegetation.
The invention selects two important factors, namely combustible Moisture Content (FMC) and Normalized Difference Vegetation Index (NDVI), as training data of a fire risk grade prediction model. The FMC is a ratio obtained by dividing the difference between the wet weight and the dry weight of the plants in the sampling unit by the dry weight, reflects the water content of the leaves of the unit vegetation, directly influences the fire risk grade, is an important fire induction factor, and can be obtained through inversion of a vegetation radiation transmission model in the prior art. The NDVI is a common vegetation index, can reflect the growth condition and the vegetation coverage of vegetation, and can be obtained by the first and second wave band calculation of a remote sensing satellite reflectivity product (MOD09A 1). Embodiments of the present invention will be described in detail below with reference to these two vegetation parameters as examples:
example (b):
a method for extracting non-fire data based on a semi-variogram comprises the following steps:
the method comprises the following steps: preparing data;
this example is detailed in 1 month 2016 in Yunnan province, Maitreya:
the present embodiment provides a fire distribution image using a remote sensing satellite product (MCD64a1) that records the time and place of each month the fire occurred; because of the difference between the vegetation types, different vegetation types should be treated separately, and the vegetation type of each fire pixel is judged by using a land cover product (MCD12Q1) provided by a remote sensing satellite.
Step two: spatial decorrelation;
statistics on MCD64A1 fire distribution products shows that 46 pixels in Mueller county 1 month in 2016 have fires, and in the embodiment, a fire buffer area is established by calculating the buffer radius of each pixel according to FMC and NDVI. Referring to fig. 1, the left image is a fire distribution image in the mile county of Yunnan province, and the right image is a schematic buffer radius diagram of a certain fire pixel therein. Dark squares in the circles indicate fire pixels, light squares in the circles and light squares outside the circles indicate pixels without fire, and areas occupied by the circles indicate buffer areas. The pixels within the circle are within the buffer area and therefore cannot be selected as non-fire pixels, but only pixels outside the buffer area are selected as non-fire pixels. The operation of creating a buffer for only one of the fire pixels will be described in detail below, and the processing for the remaining fire pixels is the same.
The invention mainly carries out spatial decorrelation processing based on a half-mutation function:
the semivariation function is also called semivariance function, is a key function for researching soil variability in geostatistics, and is a continuous function for describing spatial continuous variation of soil properties, and can reflect the variation between different distance observed values of soil properties, and the basic formula of the semivariance function is shown as the following formula (1).
Figure BDA0001537882440000051
Wherein R is(h)Is a half-variation function, N is the number of pixels at h from the fire point, Z(x)Is the pixel value, Z, of the fire pixel(xi+h)Is the pixel value at the pixel h away from the fire point.
The present embodiment uses a spherical model to fit half-value variation function values. The spherical model is also referred to as a Materon model. It is known in the art that the distance at which a model first assumes a horizontal state is a variable distance, and a sample position separated by a distance closer than the variable distance is spatially autocorrelation, while a sample position farther than the variable distance is not spatially autocorrelation, and the basic formula is shown in equation (2).
Figure BDA0001537882440000052
Wherein, YhIs an independent variable, h is a dependent variable, C0C, a are model parameters, and parameter a represents a variation.
Calculating half variation function values of FMC and NDVI with distance according to formula (1), wherein the half variation function values of NDVI with distance are shown in FIG. 2, and the half variation function values of FMC with distance are shown in FIG. 3; then, the model parameter a (namely the variation) is obtained by fitting the formula (2), so that the variation of the NDVI is 33 pixel distances, and the variation of the FMC is 35 pixel distances. At this time, the fire pixel has two buffer radiuses, and a larger buffer radius is selected to establish a fire buffer area of the fire pixel, so that a superposition schematic diagram of the fire buffer area shown in fig. 4 is obtained, wherein for the same fire pixel, the inner circle is an NDVI buffer area, and the outer circle is an FMC buffer area. Similarly, the above processing procedure is repeated for the rest fire pixels, and then the corresponding fire buffer areas are established, so that the fire buffer area superposition diagram shown in fig. 5 can be obtained, and the selection of the non-fire pixels cannot fall into any one of the buffer areas.
Step three: time decorrelation;
in practical situations, there are strong contingencies in the occurrence or non-occurrence of a fire. For a picture element, neither a fire has occurred in that phase nor is it within the buffer of any fire picture element in that phase, but its data indicates that the fire risk level for that picture element is high, most likely because a fire has occurred in its adjacent phase or the corresponding phase of another year. If the fire occurrence condition of the picture element is not considered in the adjacent time phase and the corresponding time phase of other years, the picture element is possibly selected as a non-fire picture element, so that the difference between fire and non-fire data is reduced.
Based on the multi-temporal characteristics of remote sensing images, with months as a time scale, 2016 (year 1) as a target time phase, 2015 (year 12), 2016 (year 2) and 2001-2017 (MCD64A1 begins at year 11) and yearly 1 as related time phases, extracting fire distribution images of the related time phases, and then repeating the operation of the step two to establish a fire buffer area of a fire pixel in each related time phase. As shown in fig. 6, the two-step fire buffer and the three-step fire buffer are superimposed to obtain an image, where (a) is a schematic diagram of the buffer at the target time phase, (b) is a schematic diagram of the buffer at the previous time phase (i.e., 12 months 2015), fig. (c) is a schematic diagram of the buffer at the next time phase (i.e., 2016 months 2), fig. (d) is a schematic diagram of the buffer at month 2001, fig. (e) is a schematic diagram of the buffer at month 2002, fig. (f) is a schematic diagram of the buffer at month n, and fig. (g) is a schematic diagram of the superimposition of the target time phase and all the buffers at the time phases related to the target time phase.
Step four: extracting non-fire data;
in fig. 6, 46 pixels which do not fall into the fire buffer are randomly selected as the comparison data of the fire data to ensure the extraction quality of the non-fire data.
While the present invention has been described with reference to the embodiments illustrated in the drawings, the present invention is not limited to the embodiments, which are illustrative rather than restrictive, and it will be apparent to those skilled in the art that many more modifications and variations can be made without departing from the spirit of the invention and the scope of the appended claims.

Claims (3)

1. A method for extracting non-fire data based on a semi-variogram is characterized by comprising the following steps:
step 1: determining the position and the occurrence moment of a fire pixel in the combustion area data;
step 2: determining the vegetation type of the fire pixel;
and step 3: extracting data of fire pixels and surrounding non-fire pixel data with the same type as the vegetation type of the fire pixels;
and 4, step 4: performing spatial correlation analysis on the fire pixels by using a half-variation function, taking a variation range as a buffer distance of the fire pixels and establishing a corresponding fire buffer area; the specific method comprises the following steps: acquiring a variable range by adopting a spherical model to fit a half-variation function value, and calculating a radius of each fire pixel according to FMC and NDVI by taking the fire pixel as a circle center to establish a fire buffer area;
and 5: extracting data of a relevant time phase of a fire pixel;
step 6: repeating the step 4 to establish a fire buffer zone of a time phase related to the fire pixel;
and 7: and (4) selecting pixels which do not fall into the fire buffer areas obtained in the step (4) and the step (6) and extracting data of the pixels to serve as non-fire data.
2. The method for extracting non-fire data based on the semi-variogram according to claim 1, wherein the combustion region data in step 1 is specifically remote sensing image data.
3. The method for extracting non-fire data based on the semi-variogram according to claim 1 or 2, wherein the relevant time phases of the fire pixels in step 6 are: a single fire pixel is used as a research object, the moment of fire of the pixel is used as a target time phase, and the adjacent time phases of the target time phase and the time phases corresponding to the target time phase in the rest years are related time phases.
CN201810003693.8A 2018-01-03 2018-01-03 Method for extracting non-fire spatio-temporal data based on semi-variogram Active CN108230310B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810003693.8A CN108230310B (en) 2018-01-03 2018-01-03 Method for extracting non-fire spatio-temporal data based on semi-variogram

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810003693.8A CN108230310B (en) 2018-01-03 2018-01-03 Method for extracting non-fire spatio-temporal data based on semi-variogram

Publications (2)

Publication Number Publication Date
CN108230310A CN108230310A (en) 2018-06-29
CN108230310B true CN108230310B (en) 2021-12-17

Family

ID=62645236

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810003693.8A Active CN108230310B (en) 2018-01-03 2018-01-03 Method for extracting non-fire spatio-temporal data based on semi-variogram

Country Status (1)

Country Link
CN (1) CN108230310B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1763560A (en) * 2005-10-20 2006-04-26 中国农业科学院农业资源与农业区划研究所 Method based on MODIS data automatic steppe burn scar detection
CN102819829A (en) * 2012-07-17 2012-12-12 中国科学院地理科学与资源研究所 Rebuilding algorithm for super-resolution remote sensing image based on fractal theory
CN103632065A (en) * 2013-12-19 2014-03-12 中国科学院地理科学与资源研究所 Simulated-annealing-based method for validating precision of pixel value of remote-sensing product of complex ground surface space
CN103700092A (en) * 2013-12-04 2014-04-02 中国科学院遥感与数字地球研究所 Forest burned area automatic extraction method based on time sequence remote sensing image
CN103902802A (en) * 2012-12-29 2014-07-02 中国科学院深圳先进技术研究院 Method for rebuilding spatial-information-considered vegetation index time series data
CN104851087A (en) * 2015-04-17 2015-08-19 华中农业大学 Multi-scale forest dynamic change monitoring method
CN105279738A (en) * 2015-07-15 2016-01-27 电子科技大学 Coping method of shadow problem in vegetation parameter estimation, based on remote sensing images
CN106296630A (en) * 2015-05-18 2017-01-04 中国科学院遥感与数字地球研究所 A kind of website LAI observes spatial representative evaluation methodology in Remote Sensing Products grid cell size

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9952042B2 (en) * 2013-07-12 2018-04-24 Magic Leap, Inc. Method and system for identifying a user location

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1763560A (en) * 2005-10-20 2006-04-26 中国农业科学院农业资源与农业区划研究所 Method based on MODIS data automatic steppe burn scar detection
CN102819829A (en) * 2012-07-17 2012-12-12 中国科学院地理科学与资源研究所 Rebuilding algorithm for super-resolution remote sensing image based on fractal theory
CN103902802A (en) * 2012-12-29 2014-07-02 中国科学院深圳先进技术研究院 Method for rebuilding spatial-information-considered vegetation index time series data
CN103700092A (en) * 2013-12-04 2014-04-02 中国科学院遥感与数字地球研究所 Forest burned area automatic extraction method based on time sequence remote sensing image
CN103632065A (en) * 2013-12-19 2014-03-12 中国科学院地理科学与资源研究所 Simulated-annealing-based method for validating precision of pixel value of remote-sensing product of complex ground surface space
CN104851087A (en) * 2015-04-17 2015-08-19 华中农业大学 Multi-scale forest dynamic change monitoring method
CN106296630A (en) * 2015-05-18 2017-01-04 中国科学院遥感与数字地球研究所 A kind of website LAI observes spatial representative evaluation methodology in Remote Sensing Products grid cell size
CN105279738A (en) * 2015-07-15 2016-01-27 电子科技大学 Coping method of shadow problem in vegetation parameter estimation, based on remote sensing images

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Combining AVHRR and meteorological data for estimating live fuel moisture content;MarianoGarcía et al;《Remote Sensing of Environment》;20080915;第112卷(第9期);3618-3627 *
The wildland fuel cell concept: an approach to characterize fine-scale variation in fuels and fire in frequently burned longleaf pine forests;J. Kevin Hiers等;《International Journal of Wildland Fire》;20070529;第18卷(第3期);315-325 *
三峡库区森林植被生物多样性遥感定量监测方法研究;徐文婷;《万方数据知识服务平台》;20050425;1-141 *
基于半变异函数的多极化SAR图像地表淹没程度分析;沈国状等;《遥感技术与应用》;20060316;第20卷(第6期);569-573 *
空间统计学及其在森林图形与图像处理中应用的研究;冯益明;《中国博士学位论文全文数据库 (农业科技辑)》;20050215(第2期);D049-1 *

Also Published As

Publication number Publication date
CN108230310A (en) 2018-06-29

Similar Documents

Publication Publication Date Title
Collins et al. Training data requirements for fire severity mapping using Landsat imagery and random forest
Page et al. A harmonized image processing workflow using Sentinel-2/MSI and Landsat-8/OLI for mapping water clarity in optically variable lake systems
Gizachew et al. Mapping and estimating the total living biomass and carbon in low-biomass woodlands using Landsat 8 CDR data
Barakat et al. Monitoring of forest cover dynamics in eastern area of Béni-Mellal Province using ASTER and Sentinel-2A multispectral data
Achard et al. Pan-tropical monitoring of deforestation
De Luca et al. A workflow based on Sentinel-1 SAR data and open-source algorithms for unsupervised burned area detection in Mediterranean ecosystems
Ndalila et al. Geographic patterns of fire severity following an extreme eucalyptus forest fire in southern Australia: 2013 Forcett-Dunalley fire
Levanoni et al. Can we predict butterfly diversity along an elevation gradient from space?
Grzegozewski et al. Mapping soya bean and corn crops in the State of Paraná, Brazil, using EVI images from the MODIS sensor
Meng et al. A linear mixed-effects model of biomass and volume of trees using Landsat ETM+ images
Pandit et al. Exploring the inclusion of Sentinel-2 MSI texture metrics in above-ground biomass estimation in the community forest of Nepal
Chen et al. A simple method for detecting phenological change from time series of vegetation index
Kolarik et al. A multi-plot assessment of vegetation structure using a micro-unmanned aerial system (UAS) in a semi-arid savanna environment
Safanelli et al. Leveraging the application of Earth observation data for mapping cropland soils in Brazil
Reimer et al. Advancing reference emission levels in subnational and national REDD+ initiatives: a CLASlite approach
Castanho et al. A close look at above ground biomass of a large and heterogeneous Seasonally Dry Tropical Forest-Caatinga in North East of Brazil
Melgar-García et al. A new big data triclustering approach for extracting three-dimensional patterns in precision agriculture
Melo et al. Striking divergences in Earth Observation products may limit their use for REDD+
Aynekulu et al. Review of methodologies for land degradation neutrality baselines: Sub-national case studies from Costa Rica and Namibia
De Grandi et al. Wavelet based analysis of TanDEM-X and LiDAR DEMs across a tropical vegetation heterogeneity gradient driven by fire disturbance in Indonesia
Nguyen et al. Mapping Land use/land cover using a combination of Radar Sentinel-1A and Sentinel-2A optical images
Migolet et al. Estimation of the total dry aboveground biomass in the tropical forests of Congo Basin using optical, LiDAR, and radar data
Du et al. Mapping annual global forest gain from 1983 to 2021 with landsat imagery
Aszalós et al. Accurate prediction of ice disturbance in European deciduous forests with generalized linear models: a comparison of field-based and airborne-based approaches
Kranz et al. Earth observation based multi-scale assessment of logging activities in the Democratic Republic of the Congo

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