Annual Deforestation Mapping in Sumatera using multi temporal digital classification

Transcription

1 Annual Deforestation Mapping in Sumatera using multi temporal digital classification Kustiyo Remote Sensing Technology and Data Center LAPAN International Workshop on Land Use/Cover Change and Air Pollution in Asia, August 4 7th, 2015th, Bogor, Indonesia.

2 OUTLINE 1. Background 2. Objective 3. Data 4. Classification methodology 5. Result 6. Conclusion

3 Background 1. Needs: Annual deforestation Information from for FREL Submission 2. Available deforestation Information in Indonesia 1990, 1996, # Ministry of Environment and Forestry (Official) LAPAN (INCAS) MARGONO (SDSU) Classes 23 Class Forest/Non Forest Deforestation Periode 1990, 1996, 2000, 2003, 2006, 2009, 2011, 2012,2013 Metode Visual interpretation Digital Classification (semi automatic) Forest Definition Natural and secondary forest (annual) (annual) Tree cover ( >5 meter height, < 30% cover) Digital Classification Natural Forest Table: Available deforestation Information in Indonesia

4 Background 3. Landsat data availability of Indonesia Processed data for INCAS ( ) Number of Scene Sum 1399 Kal 1045 Sul 1062 Papua 1069 other 1614 Total 6189 Actual temporal density Landsat 5/7 Number of cloud free image (years : ) Raw Landsat data ( ) Tahun LS 5 Total TIFF_1T TIFF_1G Total

5 Objectives: Finding the algoritm for producing the annual deforestation information over Sumatera island using automatic processing Data: Annual Landsat 5/7 data 1990 to 2012 (main data) Forestmap 1990 (forest base map)

6 Processing Flow Landsat-5/7- L1G Product PREPROCESSING Annual Mosaic Reflectance Product INDEX CALCULATION and FILTERING Annual Deforestation Index CLASSIFICATION Annual Deforestation

7 Processing Flow (Preprocessing) Landsat-5/7- L1G Product Geometric and Radiometric Correction 1. Ortho Rectifications 2. Top of Atmesfer (ToA) and Bidirectional Reflectance Distribution Function (BRDF) Cloud masking Reflectance Product Cloud mask Masking Masked Masked Reflectance Reflectance Product Product Masked Reflectance Product Mosaik Annual Mosaiic Reflectance Product

8 Preprocessing results (Mosaic of Landsat data )(1/2)

9 Preprocessing results (Mosaic of Landsat data )(2/2)

10 Processing Flow (Classification) Annual Mosaiic Reflectance Annual Mosaiic Product Reflectance Annual Mosaiic Product Reflectance Annual Mosaiic Product Reflectance Product Deforestation Index calculation and filtering Annual Mosaiic Reflectance Annual Mosaiic Product Reflectance Annual Mosaiic Product Reflectance Annual Product Deforestation Index Alternative metodology: 1. DI= SWIR1 / NIR 2. NBR=(NIR SWIR2)/(NIR-SWIR2) 3. NDVI= (NIR-RED)/ (NIR+RED) Classification Annual Mosaiic Reflectance Annual Mosaiic Product Reflectance Annual Mosaiic Product Reflectance Annual Product Deforestation DI : Disturbance Index NBR : Normalized Burn Ratio NDVI : Normalized Difference Vegetation Index Forest base (1990) Def Index= (NIR - SWIR1) / (NIR + SWIR1) Annual Mosaiic Reflectance Annual Mosaiic Product Reflectance Annual Mosaiic Product Reflectance Final Annual Product Deforestation

11 Deforestation Index Index= (NIR - SWIR1) / (NIR + SWIR1) Caracteristics : Higher value for vegetation, less value for open area Less sensitive from Haze Sensitive with the logging Correct the terrain effect

12 Deforestation Index (Landsat8: PR , Acq.Date ) RGB (SWIR1, NIR, RED) INDEX(NIR-SWIR1)/(NIR+SWIR1) Higher value for vegetation, less value for open area

13 Deforestation Index (Landsat8: PR , Acq.Date ) INDEX: Open Area bush low-veg dence-veg

14 INDEX: Detect logging activities RGB (SWIR1, NIR, RED) INDEX(NIR-SWIR1)/(NIR+SWIR1)

15 INDEX: Less efected by haze RGB (SWIR1, NIR, RED) INDEX(NIR-SWIR1)/(NIR+SWIR1) RGB (RED, GREE, BLUE)

16 INDEX: Less impacted by terrain effect Band Blue Band Green Band Red Band NIR Band SWIR-1 INDEX

17 Deforestation Index (trend )

18 Deforestation Index (trend ) Landsat

19 Deforestation Classification methodology Y1 Y2... Yi.... Yn Index Stable forest i...n Deforestation on year-i Forest thershold Year Index Deforestation Images

20 Result Forest map in 1990 from MoEF Result of annual deforestation ( ) Forest in 1990 Forest in 2012 Deforestation

21 Result: Annual deforestation ( ) Note: The result was mask by forest in 1990 from MoEF Deforestation Forest in 2012

22 Result: Annual deforestation in Sumatera (km 2 ) The biggest deforestation in Sumatera is , It is related with strong Elnino 1997

23 Comparison of deforestation area in Sumatera between result and Ministry of Environment and Forestry KLHK DI Methode Interval year

24 Conclusions 1. This methode detected the biggest deforestation in Sumatera : In general, the result shows that deforestation trend is similar compare with deforestation from Ministry of environment and Forestry. 3. Improving the result using other information (Elevation and Slope) 4. Accuracy assesment is needed

25 THANK YOU