Introduc)ontoAsia/RiCE* TechnicalDemonstra)onSites*II.** Japan,Indonesia,Vietnam*

Size: px

Start display at page:

Download "Introduc)on*to*Asia/RiCE* Technical*Demonstra)on*Sites*II.** Japan,*Indonesia,*Vietnam*"

Terence Parks
5 years ago
Views:

Introduc)on*to*Asia/RiCE* Technical*Demonstra)on*Sites*II.** Japan,*Indonesia,*Vietnam* Kei$Oyoshi 1,$Nobuhiro$Tomiyama 2,$Toshio$Okumura 2,$Shinichi$Sobue 1,2,$ Wahyu$Supriatna 3,$Rizatus$Shofiya?

1 Introduc)on*to*Asia/RiCE* Technical*Demonstra)on*Sites*II.** Japan,*Indonesia,*Vietnam* Kei$Oyoshi 1,$Nobuhiro$Tomiyama 2,$Toshio$Okumura 2,$Shinichi$Sobue 1,2,$ Wahyu$Supriatna 3,$Rizatus$Shofiya? 3,$Lam$Dao$Nguyen 4 $and$thuy$le$toan 5 *! 1$ Japan$Aerospace$Explora?on$Agency$(JAXA) 2$ Remote$Sensing$Technology$Center$(RESTEC) 2$ Indonesian$Center$for$Agricultural$Land$Resources$Research$and$Development$(ICALRD) 4$ Vietnam$Southern$Satellite$Technology$Applica?on$Center$(VAST/STAC) 5 Centre$d'Etudes$Spa?ales$de$la$BIOsphère$(CESBIO) 7th$GEOSSYAP$Symposium@Tokyo,$Japan 27$May$2014

2 AsiaYRiCE$:$Aims$of$TDS$Ac?vi?es! Each$technical$demonstra?on$site$covers$approximately$100$ square$kilometers$or$a$single$province$and$represents$one$of$the$ key$components$of$asiayrice.$! $They$seek$to$develop,$exercise$and$refine$the$process$ $from$ stakeholder$engagement,$through$resource$marshalling$and$ deployment,$and$on$to$trial$implementa?on$and$product$ genera?on$using$sar$and$other$earth$observa?on$data$for$rice$ crop$monitoring.$ [AsiaYRiCE$GEOGLAM$Phase$1$Space$Data$ Requirements$Technical$Demonstrator$Site$ Characteriza?on,$2013]$ 2

3 Target$Products$for$AsiaYRiCE Enhance$our$rice$monitoring$ability$through$the$use$of$blending$ various$earth$observa?on$data,$especially$the$u?liza?on$of$sar$data$ ID Product Descrip)on P1 P2 P3 P4 P5 Rice*Crop*Area*Map Crop*Calendars Crop*Damage*Assessment Agro/met*informa)on @Nowcast$and$anomaly$of$radia?on,$land$ surface$temperature,$precipita?on,$soil$ with$satellite$derived$agroymet$ [AsiaYRiCE$Team,$2012,$$hdp:// 3

4 Technical*Demonstra)on*Sites Japan India China Taiwan Thailand Laos Vietnam*(North) Philippines Vietnam*(South) Malaysia Indonesia Phase*1A:*Jun*2013*/ Phase*1B:*Apr*2014*/ 4 4

5 World$Rice$Produc?on$(2011)! Asian$countries$produce$and$consume$90%$of$global$rice.$ Top*20*Rice*Produc)on*Countries*in*2011 China India Indonesia Bangladesh Vietnam Thailand Myanmar Philippines Brazil Cambodia Japan [FAOSTAT,$2012] 5

6 Overview$of$Progress$in$TDS$Ac?vi?es! Japan*(Yamagata*Pref.) Y$Field$survey$and$installing$AWS$system$with$camera Y$Acquisi?on$of$SAR$data$(RadarsatY2)$ Y$Development$of$paddy$field$mapping$sojware! Indonesia*(West*Java,*South*Kalimantan) Y$Field$Survey Y$Acquisi?on$of$SAR$data$(RadarsatY2) Y$Airborne$SAR$observa?on$and$phenology$classifica?on Y$Preparing$SAR$data$processing$for$paddy$field$mapping$!! Vietnam/South*(An*Giang*Prov.) Y$Field$Survey Y$Acquisi?on$of$SAR$data$(CSK,$RadarsatY2) Y$Paddy$field$mapping$by$SAR$(CSK)$and$valida?on Y$Yield$es?ma?on$by$SAR$(CSK)$and$valida?on 6

7 Japan: Yamagata Pref.

600,000"" 2012 500,000"" 400,000"" 300,000"" 200,000"" 100,000"" 0""

Chiba" Iwate" Tokyo" Rice*Crop*Calendar*in*Yamagata*Prefecture* Month Jan

8 Loca?on$of$Technical$Demonstra?on$Site 90x150km Top*10 Rice*Produc)on*Prefecture*(t)* 700,000"" 600,000"" ,000"" 400,000"" 300,000"" 200,000"" 100,000"" 0"" Nigata" Hokkaido" Akita" Ibaraki" Yamagata" Miyagi" Fukushima" Tochigi" Chiba" Iwate" Tokyo" Rice*Crop*Calendar*in*Yamagata*Prefecture* Month Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec Calendar Plan?ng Harves?ng Transplan?ng,$Irrigated$cropping,$Yield$=$approx$6.0$t/ha 8

9 Rice$Cul?vated$Area$Mapping$Sojware! INAHOR$:$INterna?onal$Asian$Harvest$mOnitoring$system$for$Rice 1.$Main$Window 2.$Select$SAR$data$for$Processing 3.$Rice$Cul?vated$Area$Mapping$ 4.$Calculate$Mapping$Area This*so]ware*enable*us*to*easily*map*rice*cul)vated*area*from*)me/series*SAR*data. 9

10 Paddy$Field$Mapping$by$ALOS$PALSAR$over$Yamagata Google$Eearth Tsuruoaka*City

11 Paddy$Field$Mapping$by$ALOS$PALSAR$over$Yamagata ALOS$PALSAR$FBS$ 2009 :$Planting$Season$(Min) Tsuruoaka*City

12 Paddy$Field$Mapping$by$ALOS$PALSAR$over$Yamagata ALOS$PALSAR$FBS$ 2009 :$Vegetative$Season$(Max) Tsuruoaka*City

13 Paddy$Field$Mapping$by$ALOS$PALSAR$over$Yamagata ALOS$PALSAR$FBS$ 2009 :$Range$(Max$Y$Min) Tsuruoaka*City

14 Paddy$Field$Mapping$by$ALOS$PALSAR$over$Yamagata ALOS$PALSAR$FBS$ 2009 Tsuruoaka*City

15 Field$Survey$for$Valida?on! Landcover$PreYsurvey$(Apr$2014)$ Paddy$field,$soybeans,$agricultural$land$abandonment$etc. Landcover*survey*will*be*conducted*in*the*beginning*of*Jun. 15

Meteorological$sensor Y$Air$Temperature Y$Precipita?on$(rain$gauge) Y$Radia?

16 Field$Router$(Automa?c$Weather$Sta?on)! Characteris?cs Y$Automa?c$data$collec?on Y$Wireless$communica?on$(daily) Y$Powered$by$solar$panel! Meteorological$sensor Y$Air$Temperature Y$Precipita?on$(rain$gauge) Y$Radia?on Y$Wind$direc?on$and$velocity! Image$sensor Y$CMOS$camera$(120$M) [Developed$by Prof.$Mizoguchi,$U.Tokyo] 16

17 Near$realY?me$inYsitu$and$Satellite$Monitoring$System! Tsuruoka$City$in$Yamagata$Prefecture$ NDVI$(Vegeta?on)$ NDWI$(Water)$ NDSI$(Snow) AWS*with*Camera Near*real/)me*photo*and*satellite*monitoring*system 17

18 NDVI$(MOD13A1$500m)$vs$FR$Photo N ! E DOY = 348! (Snowing) DOY = 296! (tiller) DOY = 162! (installing syatem) DOY = 195! (Max NDVI) DOY = 244! DOY = 264! (Harvest) SAR*observa)on*data*(ALOS/2*etc.)*will*be*included*in*this*comparison*in*

19 Current$TDS$Tsuruoka$Site$:$20th$May$2014 NDVI Planting Stage NDWI NDSI NDVI NDWI 19

20 RADARSATY2$Acquisi?on$Plan$Over$Tsuruoka Date UTC Site Beam Pass 16YAprY2014 $08:36: Japan Wide$Fine$ Ascending 23YAprY2014 $20:46: Japan Wide$Fine$ Descending 10YMayY2014 $08:36: Japan Wide$Fine$ Ascending 17YMayY2014 $20:46: Japan Wide$Fine$ Descending 3YJunY2014 $08:36: Japan Wide$Fine$ Ascending 10YJunY2014 $20:46: Japan Wide$Fine$ Descending 27YJunY2014 $08:36: Japan Wide$Fine$ Ascending 4YJulY2014 $20:46: Japan Wide$Fine$ Descending 21YJulY2014 $08:36: Japan Wide$Fine$ Ascending 28YJulY2014 $20:46: Japan Wide$Fine$ Descending 14YAugY2014 $08:36: Japan Wide$Fine$ Ascending 21YAugY2014 $20:46: Japan Wide$Fine$ Descending 7YSepY2014 $08:36: Japan Wide$Fine$ Ascending 14YSepY2014 $20:46: Japan Wide$Fine$ Descending 1YOctY2014 $08:36: Japan Wide$Fine$ Ascending 8YOctY2014 $20:46: Japan Wide$Fine$ Descending 25YOctY2014 $08:36: Japan Wide$Fine$ Ascending 1YNovY2014 $20:46: Japan Wide$Fine$ Descending 18YNovY2014 $08:36: Japan Wide$Fine$ Ascending 25YNovY2014 $20:46: Japan Wide$Fine$ Descending Wide Fine Mode! - 5.2x7.7 m! - 170x150 km! - Inc deg! - VV+VH 20

21 Indonesia: West Java

22 Loca?on$of$Technical$Demonstra?on$Site West$Java$ (Subang,$Karawang,$Indramayu) 22

23 Airborne$SAR$Observa?on$Campaign$(PiYSARYL2) 4*August*2012 Azimuth Range 20*km 90*km 23

24 Rice Phenological Mapping by Airborne SAR Airborne SAR (Pi-SAR-L2) data analysis - Pi-SAR-L2 carries similar SAR sensor as ALOS-2

25 Phenological Stages Each sampling area was classified into 6 phenological stages based on the visual interpretation of photos. Planting Vegetated1 Vegetated2 Harvesting Harvested Drought 25

26 SAR Images of Phenological Stages R:G:B = HH:HV:VV VV HH HV Planting Vegetated1 Vegetated2 Harvesting Harvested Drought 26

27 27

28 Rice Phenology Mapping by Airborne SAR 18 Planting Vegetated$1 Vegetated$2 RIpening 28

Classification Result in Batola, South Kalimantan Courtesy:$Wahyu$Supriatna$ and$$rizatus$shofiya?$ Class$distribution$for$selected$area No Class Percent 1 rice150days_kb3 8.818 2 rice105days_kb5 23.

29 Classification Result in Batola, South Kalimantan Courtesy:$Wahyu$Supriatna$ and$$rizatus$shofiya?$ Class$distribution$for$selected$area No Class Percent 1 rice150days_kb rice105days_kb rice140days_kb water mixedtreecrops others nodata Total Overall$Class$Performance$:$ (13953$/$16954$)$=$$82.3%$ Kappa$Statistic$(X100)$=$$73.4% 29

30 Dry$Biomass$vs$Backscadering$Coef. HH HV r=0.85 r=0.62 VV Dry*Biomass n$=$37 Backsca`ering*Coef.*(dB) r=0.74 HH*showed*highest*correla)on*coefficient. ALOS/2*ScanSAR*mode*has*HH*mode. 30

Fine W2 RADARSAT-2 Data (Descending) Obs.

31 SAR Data Collection : Radarsat-2 RADARSAT-2 Data (Ascending) Obs.*date Orbit Mode Beam 01#Dec#2013 Asc Fine W2 15#Dec#2013 Asc Fine W2 18#Jan#2014 Asc Fine W2 11#Feb#2014 Asc Fine W2 RADARSAT-2 Data (Descending) Obs.*date Orbit Mode Beam 23#Aug#2013 Des Fine W2 03#Nov#2013 Des Fine W2 27#Nov#2013 Des Fine W2 21#Dec#2013 Des Fine W2 14#Jan#2014 Des Fine W2 07#Feb#2014 Des Fine W2 Preliminary result Processed by INAHOR Rice cultivated area mapping 31

32 Vietnam: An Giang Prov.

*remote*sensing/based* technology*has*not*been*applied*in*practice* to*replace*traditional*ground/based*survey* methods.

33 6 th SAFE Workshop Kuala Lumpur, 23 rd April 2014 Background: Vietnam*is*one*of*the*leading* rice*producers*and*exporters*in*the*world.* Monitoring,*distribution*mapping,*and*yield* estimating*of*rice*crop,*hence,*have*become* the*imperative*tasks.*remote*sensing/based* technology*has*not*been*applied*in*practice* to*replace*traditional*ground/based*survey* methods.*this*is*because*of**limitations*of* image*processing*techniques,*and*available* EO*data*obtained,**which*have*a*remarkable* influence*on*results*achieved.* Objectives: Monitoring*of*rice*cropping*area* and*rice*growth;*and*estimating*rice*yield* and*production,*based*on*realistic*eo*data* acquisition*plan. Mekong*Delta,*Vietnam 13*provinces*and*city; Area:*40,553*Km 2 * (12%*or*1/8); Population:*17.39*M* (20%*or*1/5); Rice*production:** *Mton**************** (>*50%*or*1/2). ******Source:*GSO,* [Dr.*Lam*Dao*Nguyen*et*al.]

34 6 th SAFE Workshop Kuala Lumpur, 23 rd April 2014 AW 2013 and WS 2014 rice crop SAR data Dates of satellite data acquired: Cosmo-Skymed: 13 dates Radarsat-2: 12 dates Sensor Date*of*required Pass Cosmod\Skymed 2013\Aug\19 *DES Cosmod\Skymed 2013\Sep\04 *DES Cosmod\Skymed 2013\Sep\20 *DES Cosmod\Skymed 2013\Oct\06 *DES Cosmod\Skymed 2013\Oct\14 *DES Cosmod\Skymed 2013\Oct\22 *DES Cosmod\Skymed 2013\Oct\30 *DES Cosmod\Skymed 2013\Nov\07 *DES Cosmod\Skymed 2013\Nov\15 *DES Cosmod\Skymed 2013\Nov\23 *DES Cosmod\Skymed 2014\Jan\26 *DES Cosmod\Skymed 2014\Feb\11 *DES Cosmod\Skymed 2014\Feb\27 *DES Sensor Date*of*required Pass RADARSAT\2 2013\Aug\29 *DES RADARSAT\2 2013\Sep\22 *DES RADARSAT\2 2013\Oct\16 *DES RADARSAT\2 2013\Nov\09 *DES RADARSAT\2 2013\Dec\03 *DES RADARSAT\2 2013\Dec\08 *ASC RADARSAT\2 2013\Dec\27 *DES RADARSAT\2 2014\Jan\01 *ASC RADARSAT\2 2014\Jan\20 *DES RADARSAT\2 2014\Jan\25 *ASC RADARSAT\2 2014\Feb\13 *DES RADARSAT\2 2014\Apr\07 ASC [Dr.*Lam*Dao*Nguyen*et*al.]

6 th SAFE Workshop Kuala Lumpur, 23 rd April 2014 - Field samples: 40 -

35 6 th SAFE Workshop Kuala Lumpur, 23 rd April Field samples: 40 - Check points: 40 rice points, 24 non-rice points 9 [Dr.*Lam*Dao*Nguyen*et*al.]

6 th SAFE Workshop Kuala Lumpur, 23 rd April 2014 19 Aug

An Giang (Thoai Son & Chau Thanh districts): samples in

36 6 th SAFE Workshop Kuala Lumpur, 23 rd April Aug 2013 (2) Sample rice field MPD3 sown on 17 Aug Aug 2013 Cosmo SkyMed data: StripMap Pingpong: HHVV 10 dates (19 Aug 23 Nov)! An Giang (Thoai Son & Chau Thanh districts): samples in red. 4 Sep 2013 (18) 4 Sep [Dr.*Lam*Dao*Nguyen*et*al.] 13 Sep 2013 (27) 20 Sep 2013

37 6 th SAFE Workshop Kuala Lumpur, 23 rd April 2014 HH/VV HH/VV*(dB) \3.5 CSK_PP_ CSK_PP_ CSK_PP_ CSK_PP_ CSK_PP_ CSK_PP_ CSK_PP_ CSK_PP_ \7.0 \ Ngày*sau*khi*sạ Days after sowing 15 [Dr.*Lam*Dao*Nguyen*et*al.]

38 6 th SAFE Workshop Kuala Lumpur, 23 rd April 2014 Methods used for rice/non-rice mapping using multi-temporal CSK radar imagery 16 [Dr.*Lam*Dao*Nguyen*et*al.]

39 6 th SAFE Workshop Kuala Lumpur, 23 rd April 2014 Commune* Estimated*area* Agency*data* Percentage*error* Name (ha) (ha) (%) An*Bình 2,274 2,287 \0.55 Bình*Thành 2,313 2,325 \0.52 Định*Mỹ 3,178 2, Định*Thành 2,613 2, Mỹ Phú*Đông 2,676 2,740 \2.39 Phú*Thuận 1, Thoại*Giang 2,203 2,248 \2.03 Núi*Sập ,000 Óc*Eo \13.46 Phú*Hòa \5.48 Vĩnh*Khánh 2,500 2,810 \12.39 Vĩnh*Phú 3,069 3,153 \2.73 3,000 Vĩnh*Trạch 1,269 1,424 \12.19 Vĩnh*Chánh 2,357 1, Vọng*Đông 2,601 2, Agency*data*(ha) 2,000 Accuracy assessment of rice/non-rice mapping method for AW 2013 crop in Thoai Son district using two-date CSK data (20 Sep and 14 Nov) based on statistical data y*=*0.9639x*+* R²*=* , [Dr.*Lam*Dao*Nguyen*et*al.] 0 0 1,000 2,000 3,000 4,000 Es)mated*area*(ha)

6 th SAFE Workshop Kuala Lumpur, 23 rd April 2014 Y Ra = 6.4575-2.

282 ton/ha Where: Y Ra : rice yield (ton/ha), Ra 1 : polarisation ratio of

polarisation ratio of third date image, 25 r 2 : the coefficient of

A*distribution*map*of*estimated*rice*yield*in*AW*

40 6 th SAFE Workshop Kuala Lumpur, 23 rd April 2014 Y Ra = *Ra *Ra *Ra 3 r 2 = 0.802, se y = ton/ha Where: Y Ra : rice yield (ton/ha), Ra 1 : polarisation ratio of first date image, Ra 2 : polarisation ratio of second date image, Ra 3 : polarisation ratio of third date image, 25 r 2 : the coefficient of determination, se y : the standard error for the y estimate. A*distribution*map*of*estimated*rice*yield*in*AW* 2013*crop*at*Chau*Thanh*and*Thoai*Son*district* using*three\date*polarisation*ratios: *CSK*04*Sep*013* *CSK*20*Sep*2013 *CSK*06*Oct*2013 [Dr.*Lam*Day*Nguyen*et*al.]

$ALOS\2*:*New*SAR*Satellite ALOS/2*satellite* *Launch** *Orbit*type** Altitude** *Revisit*time** LSDN* * :*24th*May*2014* :*Sun\synchronous$ $+/\*500*m*$

41 ALOS\2*:*New*SAR*Satellite ALOS/2*satellite* *Launch** *Orbit*type** Altitude** *Revisit*time** LSDN* * :*24th*May*2014* :*Sun\synchronous :*628*km*+/\*500*m*(for*reference*orbit) :*14*days :*12:00*+/\*15*min Solar*paddles SAR*antenna Z Y X 12:05 p.m. on May 24 PALSAR/2** *L/band*Synthetic*Aperture*Radar *Active*Phased*Array*Antenna*type ****two*dimensions*scan*(range*and*azimuth) *Antenna*size*:*3m(El)*x*10m(Az) *Bandwidth*:*14* *84MHz *Peak*transmit*Power**:*5100W *Observation*swath*:*25* *490km *Resolution*:*Range:*3*m*to*100*m Azimuth:*1*m*to*100*m Successfully*launched*on*24th*May*2014*and*antenna*deployment*was*completed. 41

42 ALOS\2*Observaton*Mode! ALOS*has*several*observaton*modes*different*in*polarizaton,** spatal*and*actual*temporal*resolutons.* Spotlight Ultra Fine High sensitive Fine ScanSAR nominal ScanSAR wide Bandwidth 84MHz 84MHz 42MHz 28MHz 14MHz 28MHz 14MHz Resolution Swath Rg Az 3 1m Rg Az 25 25km 3m 6m 10m 100m 60m 50km 50km 70km 350km (5-scan) 490km (7-scan) Polarization SP SP/DP SP/DP/FP/CP SP/DP S/A NESZ -24dB -24dB -28dB -26dB -26dB -23dB -23dB SP : HH or VV or HV, DP : HH+HV or VV+VH, FP : HH+HV+VH+VV, CP : Compact pol Rg 25dB 25dB 23dB 25dB 25dB 20dB Az 20dB 25dB 20dB 23dB 20dB 20dB FP 1 cov/yr HH+HV 6 cov/yr (SE Asia) HH+HV 9 cov/yr (SE Asia) ScanSAR*mode*can*provide*data*every*42/day,* *ScanSAR*is*useful*for*periodical**monitoring*with*provincial/country/level. 42

43 ! Progress* Summary! Field*survey*is*conductng*periodically*in*three*TDS.*! Preliminary*paddy*field*map*were*generated*and*in*validaton.*! Rice*yield*estmaton*by*SAR*was*conducted*(Vietnam).*! Radarsat\2*(wide*fine*mode)*data*are*acquired.*!! Future*Works*! Contnuaton*of*acquiring*satellite*and*in\situ*data*! Refinement*of*the*algorithms*to*estmate*paddy*field*area* and*yield*estmaton*based*on*the*in\situ*data*! Demonstraton*of*INAHOR*for*paddy*field*mapping*! Preparaton*of*the*utlizing*new*SAR*data*including*ALOS\2* Sentnel\1*and*RISAT\1 43

44 Thank You for Your Attention!!

The Current Status and Brief Results of Engineering Model for PALSAR-2 onboard ALOS-2 and Science Project

The Current Status and Brief Results of Engineering Model for PALSAR-2 onboard ALOS-2 and Science Project + The 16 th KC meeting Japan Aerospace Exploration Agency Masanobu Shimada, Yukihiro KANKAKU The