Rice crop monitoring in Vietnam using radar remote sensing data

Transcription

1 Earth Observation Technologies for Crop Monitoring: A Workshop to Promote Collaborations among GEOGLAM/JECAM/Asia-RiCE September 2018 Taichung, Taiwan Rice crop monitoring in Vietnam using radar remote sensing data Lam Dao Nguyen 1, Thuy Le-Toan 2, Hoang Phi Phung 1, Dang Diem Huong 1, Nguyen Van Anh Vu 1, Pham Duy Tien 3, Nguyen Duy Khang 1, Do Thanh Long 1 1 HCMC Space Technology Application Center (STAC) VNSC VAST, Vietnam 2 Centre d'etudes Spatiales de la BIOsphère (CESBIO), France 3 An Giang University (AGU), Vietnam

2 Contents 1. Introduction 2. Rice monitoring research works 3. Rice monitoring on Vietnam data cube 2

3 3 Introduction

4 Introduction Source: Parry, M.L. et al., 2007 Mekong Delta (Source: GSO, 2015) Area: 40,576 Km2 (1/8) Population: M (~1/5) MD accounts for more than half (25.7/45.1 Mt) of the country s rice production (>1/2) MD is one of the most affected regions in the world by global warming. Studies need to be conducted to quantify the changes observed by satellites in LULC, in cultural practices, etc. 4 Food security

5 Introduction Rice cropping system Rice season Single rice crop Traditional rice (rain-fed) Sowing-transplanting period Vegetative stage Double rice crop SA AW (rain-fed) WS SA (irrigated) Triple rice crop WS SA - AW Main rice-based cropping systems in the MD Reproductive stage Ripening stage Pictures of rice growing stages Rice crop Planting Harvesting English name Local name Winter Spring (WS) Dong Xuan Nov./Dec. Mar./Apr. Summer Autumn (SA) He Thu Apr./May Jul./Aug. Rainy season Thu Dong (Autumn Winter-AW) Jul./Sep. Oct./Dec. Mua (Traditional rice) Jul./Sep. Nov./Jan. 5 Main rice seasons in An Giang province, Mekong Delta

6 Introduction Why Satellite RS? Synoptic measurements on large areas; Near real time; Low cost; etc. Optical data Why use Radar RS for Rice? All weather capability; Day and night operation; No effects of atmospheric constituents. 6 Radar data 1 Nov. 2017

7 Introduction Previous projects: Data used: ERS2-SAR data, ; ENVISAT-ASAR data, ; TerraSAR-X & ENVISAT-ASAR data, ; COSMO-SkyMed, Rice mapping: Temporal change measurement, Single-date mapping algorithm Yield estimation model: Statistical model, Agrometeorological model (AMM) Ongoing projects: Asia-RiCE & GeoRice Data used: Sentinel-1, RADARSAT-2, ALOS-2. New project: VNRice Data used: Radar and/or optical RS data 7 Objectives To evaluate the use of remote sensing data in rice monitoring & yield estimation, towards an operational system for rice crop inventory in Vietnam.

8 Research works Asia-RiCE Technical Demonstrator Site Mekong Delta, VN An Giang province, Mekong Delta, Vietnam Geographic coordinates: UL: 10 58'47.38"N, '39.51"E, UR:10 58'35.84"N, '12.84"E LR: 10 05'24.65"N, '15.36"E, LL: 10 05'45.13"N, '23.41"E 8

9 Data set SAR data received: COSMO-SkyMed: Aug 2013 Feb 2014 RADARSAT-2: Aug 2013 present Sentinel-1: Aug 2014 present ALOS-2: Nov 2014 present 9 Sentinel-1 9 September 2018 COSMO-SkyMed data: - Band: X - Polarisation: HH&VV - Resolution: 20 m (StripMap PINGPONG) RADARSAT-2 data: - Band: C - Polarisation: VV&VH - Resolution: 10 m (Wide Fine) Sentinel-1 data: - Band: C - Polarisation: VV&VH - Resolution: 20 m (IW) ALOS-2 data: - Band: L - Polarisation: HH&HV - Resolution: 50 m (WS) & 12.5 m (Fine)

10 Methods ARD SAR data Ground-truth data σ o of sampling fields Classification method Thresholding 24/03/07 In situ rice yield Regression analysis Regression equation Classification map Estimated rice yield distribution map Accuracy assessment Reference data Ground data Estimated rice production (P5) Rice/Non-rice map (P1) Asia-RiCE products 10

11 Methods 11 RADARSAT-2 data: - Band: C - Polarisation: VV&VH - Resolution: 10 m (Wide Fine)

12 Data set at provincial scale District name Agency data (ha) Estimated area (ha) Percentage error (%) An Phu 13,640 13, Cho Moi 13,304 14, Chau Phu 34,940 30, Chau Thanh 28,630 26, Phu Tan 22,151 23, Tinh Bien 16,288 21, Chau Doc 6,315 6, Long Xuyen 4,518 3, Thoai Son 38,882 37, Tri Ton 42,210 38, Tan Chau 9,321 9, Total 230, , ,000 40,000 R² = 0.97 Estimated area (ha) 30,000 20,000 10,000 SA 2016 crop from RADARSAT-2 (15 Apr, 09 May, 02 Jun, 26Jun, 20 Jul & 13 Aug) ,000 20,000 30,000 40,000 50,000 Agency data (ha) 12

13 Data set at provincial scale 250,000 A distribution map of estimated rice yield of An Giang in SA 2016 crop using RADARSAT-2 data Estimated production (ton) 200, , ,000 50,000 R² = , , , , ,000 Agency data (ton) 13 District name Agency data (ton) Estimated production (ton) Percentage error (%) An Phú 73,656 68, Chợ Mới 77,296 77, Châu Phú 201, , Châu Thành 166, , Phú Tân 134, , Tịnh Biên 89, , Châu Đốc 37,890 32, Long Xuyên 29,503 16, Thoại Sơn 233, , Tri Tôn 229, , Tân Châu 54,994 53, Total 1,325,946 1,213,

14 Data set at provincial scale A distribution map of estimated rice yield of An Giang in SA 2016 crop using ALOS-2 data District name Agency data (ton) Estimated production (ton) Percentage error (%) An Phú 73,656 63, Chợ Mới 77,296 66, Châu Phú 201, , Châu Thành 166, , Phú Tân 134, , Tịnh Biên 89,584 99, Châu Đốc 37,890 34, Long Xuyên 29,503 17, Thoại Sơn 233, , Tri Tôn 229, , Tân Châu 54,994 53, Total 1,325,946 1,154, Statistical production (ton) 250, , , ,000 50,000 R² = , , , , ,000 Estimated production (ton) 14

15 Data set at provincial scale SA 2016 crop using ALOS-2 15 SA 2016 crop using S-1 SA 2016 crop using R-2

16 Data set at regional scale Map of WS Rice 2015 Map of WS Rice Reduced area in 2016 caused by shortage of water and saline water intrusion

17 Data set at regional scale 17 Map of rice age on Feb. 17 (WS 2018 crop) using Sentinel-1 data

18 New ongoing research works Applied research on the multi-temporal, multi-resolution optical and radar remote sensing data for rice planted area monitoring and rice yield, production estimation in the Mekong Delta and Red River Delta (VNRice) State level research project, 11/ /2020 (30 months) RS Data used: SAR data, optical data, fused data Outputs: Rice crop area maps; yield & production estimation; crop calendars/ crop growth status 18

19 Vietnam Data Cube Infrastructure: by help of IMSG (I. M. Systems Group, Inc.) Software: Supports from CSIRO/CEOS Satellite Data: Landsat from USGS ALOS from JAXA Sentinel-1&-2 from ESA Applications: Forest monitoring Rice monitoring Water monitoring Etc. 19

20 Vietnam Data Cube 20 1 head node, 1 login node, 7 compute nodes (each: 16 core 2.2GHz, 64GB RAM), 40Gb/s Infiniband Interconnect network Full redundant parallel storage LustreFS: 195TB, max read/write speed: 1.9GB/s Total storage capacity (including work, home, local): ~ 220TB Can be expanded to hundreds of compute nodes, hundreds of TBs of storage

21 Rice monitoring on Vietnam Data Cube Sentinel-1 data on Vietnam data cube 21

22 Rice monitoring on Vietnam Data Cube Tool Output 22

23 Rice monitoring on Vietnam Data Cube 23 Display classified image on QGIS

24 24 Rice monitoring on Vietnam Data Cube

25 Planned products of VDC: Rice/Non-rice maps Map of WS Rice 2017 Map of SA Rice 2017 Map of AW Rice 2017 Combination of 3 rice crop maps Rice cropping system map 25

26 Planned products of VDC: Rice/Non-rice maps Aug. 10, 2018 Aug. 15, Map of SA rice in the Mekong Delta and Red River Delta

27 Planned products of VDC: Days after sowing Feb. 17, 2018 July 11, Map of rice age in the Mekong Delta

28 Planned products of VDC: Days after sowing June 16, 2018 July 10, Map of rice age in the Red River Delta

29 Planned products of VDC: Rice yield maps A distribution map of estimated rice yield 29

30 30 Towards a National crop monitor

31 31 Dr. Lam Dao Nguyen