Introduc)on*to*Asia/RiCE* Technical*Demonstra)on*Sites*II.** - - PowerPoint PPT Presentation

introduc on to asia rice technical demonstra on sites ii
SMART_READER_LITE
LIVE PREVIEW

Introduc)on*to*Asia/RiCE* Technical*Demonstra)on*Sites*II.** - - PowerPoint PPT Presentation

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


slide-1
SLIDE 1


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

Kei$Oyoshi1,$Nobuhiro$Tomiyama2,$Toshio$Okumura2,$Shinichi$Sobue1,2,$
 Wahyu$Supriatna3,$Rizatus$Shofiya?3,$Lam$Dao$Nguyen4$and$Thuy$Le$Toan5*

!

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)
 5Centre$d'Etudes$Spa?ales$de$la$BIOsphère$(CESBIO)



 7th$GEOSSYAP$Symposium@Tokyo,$Japan
 27$May$2014

slide-2
SLIDE 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.$

2

[AsiaYRiCE$GEOGLAM$Phase$1$Space$Data$ Requirements$Technical$Demonstrator$Site$ Characteriza?on,$2013]$

slide-3
SLIDE 3

Target$Products$for$AsiaYRiCE

3

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 Rice*Crop*Area*Map

@Cul?vated$area$(every$year)$
 @Inventory$of$agricultural$facili?es

P2 Crop*Calendars

@Timing$of$sowing,$plan?ng,$growing$and$ harves?ng/growing$status
 @Iden?fica?on$of$growth$stages$and$

P3 Crop*Damage*Assessment

@Detec?on$of$flooding$and$other$disaster$ impacted$area
 @Detec?on$of$drought$area


P4 Agro/met*informa)on

@Nowcast$and$anomaly$of$radia?on,$land$ surface$temperature,$precipita?on,$soil$ moisture$and$so$on.

P5 Rice*Produc)on*Es)ma)on*and* Forecas)ng

@EmpiricalYsta?s?cal$model$es?mate$ with$satellite$derived$agroYmet$


[AsiaYRiCE$Team,$2012,$$hdp://www.asiaYrice.org/files/workplan.pdf]

slide-4
SLIDE 4

4 4

India Thailand Indonesia Philippines China Japan Laos Vietnam*(North) Malaysia Phase*1B:*Apr*2014*/ Phase*1A:*Jun*2013*/

Technical*Demonstra)on*Sites

Taiwan Vietnam*(South)

slide-5
SLIDE 5

Top*20*Rice*Produc)on*Countries*in*2011 [FAOSTAT,$2012]

China India Indonesia

Bangladesh

Vietnam Thailand Myanmar Philippines Brazil Cambodia Japan

World$Rice$Produc?on$(2011)

! Asian$countries$produce$and$consume$90%$of$global$rice.$

5

slide-6
SLIDE 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

slide-7
SLIDE 7

Japan: Yamagata Pref.

slide-8
SLIDE 8

Loca?on$of$Technical$Demonstra?on$Site

8

90x150km

Month Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec

Calendar

Plan?ng Harves?ng

Rice*Crop*Calendar*in*Yamagata*Prefecture* Top*10
 Rice*Produc)on*Prefecture*(t)*

0"" 100,000"" 200,000"" 300,000"" 400,000"" 500,000"" 600,000"" 700,000"" Nigata" Hokkaido" Akita" Ibaraki" Yamagata" Miyagi" Fukushima" Tochigi" Chiba" Iwate" Tokyo"

2012 Transplan?ng,$Irrigated$cropping,$Yield$=$approx$6.0$t/ha

slide-9
SLIDE 9

Rice$Cul?vated$Area$Mapping$Sojware

1.$Main$Window

! INAHOR$:$INterna?onal$Asian$Harvest$mOnitoring$system$for$Rice

9

This*so]ware*enable*us*to*easily*map*rice*cul)vated*area*from*)me/series*SAR*data.

2.$Select$SAR$data$for$Processing 3.$Rice$Cul?vated$Area$Mapping$ 4.$Calculate$Mapping$Area

slide-10
SLIDE 10

Paddy$Field$Mapping$by$ALOS$PALSAR$over$Yamagata

Google$Eearth

Tsuruoaka*City

slide-11
SLIDE 11

Paddy$Field$Mapping$by$ALOS$PALSAR$over$Yamagata

ALOS$PALSAR$FBS$(2009):$Planting$Season$(Min)

Tsuruoaka*City

slide-12
SLIDE 12

Paddy$Field$Mapping$by$ALOS$PALSAR$over$Yamagata

ALOS$PALSAR$FBS$(2009):$Vegetative$Season$(Max)

Tsuruoaka*City

slide-13
SLIDE 13

Paddy$Field$Mapping$by$ALOS$PALSAR$over$Yamagata

ALOS$PALSAR$FBS$(2009):$Range$(Max$Y$Min)

Tsuruoaka*City

slide-14
SLIDE 14

ALOS$PALSAR$FBS$(2009)

Paddy$Field$Mapping$by$ALOS$PALSAR$over$Yamagata

Tsuruoaka*City

slide-15
SLIDE 15

Field$Survey$for$Valida?on

! Landcover$PreYsurvey$(Apr$2014)$

  • Paddy$field,$soybeans,$agricultural$land$abandonment$etc.

15

Landcover*survey*will*be*conducted*in*the*beginning*of*Jun.

slide-16
SLIDE 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)

16

[Developed$by
 Prof.$Mizoguchi,$U.Tokyo]

slide-17
SLIDE 17

Near$realY?me$inYsitu$and$Satellite$Monitoring$System

17

! Tsuruoka$City$in$Yamagata$Prefecture$

Near*real/)me*photo*and*satellite*monitoring*system

  • NDVI$(Vegeta?on)$
  • NDWI$(Water)$
  • NDSI$(Snow)

AWS*with*Camera

slide-18
SLIDE 18

NDVI$(MOD13A1$500m)$vs$FR$Photo

18

DOY = 162! (installing syatem) DOY = 195! (Max NDVI) DOY = 244! DOY = 264! (Harvest) DOY = 348! (Snowing) DOY = 296! (tiller) N38˚ 41’ 47”! E139˚ 49’ 16” SAR*observa)on*data*(ALOS/2*etc.)*will*be*included*in*this*comparison*in*2014.

slide-19
SLIDE 19

Current$TDS$Tsuruoka$Site$:$20th$May$2014

19

NDVI NDWI NDSI NDWI NDVI Planting Stage

slide-20
SLIDE 20

RADARSATY2$Acquisi?on$Plan$Over$Tsuruoka

20

Date UTC Site Beam Pass 16YAprY2014 $08:36:08.331 Japan Wide$Fine$ Ascending 23YAprY2014 $20:46:23.634 Japan Wide$Fine$ Descending 10YMayY2014 $08:36:08.331 Japan Wide$Fine$ Ascending 17YMayY2014 $20:46:23.634 Japan Wide$Fine$ Descending 3YJunY2014 $08:36:08.331 Japan Wide$Fine$ Ascending 10YJunY2014 $20:46:23.633 Japan Wide$Fine$ Descending 27YJunY2014 $08:36:08.331 Japan Wide$Fine$ Ascending 4YJulY2014 $20:46:23.633 Japan Wide$Fine$ Descending 21YJulY2014 $08:36:08.331 Japan Wide$Fine$ Ascending 28YJulY2014 $20:46:23.633 Japan Wide$Fine$ Descending 14YAugY2014 $08:36:08.331 Japan Wide$Fine$ Ascending 21YAugY2014 $20:46:23.633 Japan Wide$Fine$ Descending 7YSepY2014 $08:36:08.331 Japan Wide$Fine$ Ascending 14YSepY2014 $20:46:23.633 Japan Wide$Fine$ Descending 1YOctY2014 $08:36:08.331 Japan Wide$Fine$ Ascending 8YOctY2014 $20:46:23.633 Japan Wide$Fine$ Descending 25YOctY2014 $08:36:08.331 Japan Wide$Fine$ Ascending 1YNovY2014 $20:46:23.633 Japan Wide$Fine$ Descending 18YNovY2014 $08:36:08.331 Japan Wide$Fine$ Ascending 25YNovY2014 $20:46:23.633 Japan Wide$Fine$ Descending

Wide Fine Mode!

  • 5.2x7.7 m !
  • 170x150 km!
  • Inc. 20-40 deg!
  • VV+VH
slide-21
SLIDE 21

Indonesia: West Java

slide-22
SLIDE 22

Loca?on$of$Technical$Demonstra?on$Site

22

West$Java$


(Subang,$Karawang,$Indramayu)

slide-23
SLIDE 23

Airborne$SAR$Observa?on$Campaign$(PiYSARYL2)

23

90*km 20*km Azimuth Range 4*August*2012

slide-24
SLIDE 24

Rice Phenological Mapping by Airborne SAR

Airborne SAR (Pi-SAR-L2) data analysis


  • Pi-SAR-L2 carries similar SAR sensor as ALOS-2

slide-25
SLIDE 25

25

Phenological Stages

Each sampling area was classified into 6 phenological stages based on the visual interpretation of photos.

Harvesting Planting Vegetated2 Vegetated1 Harvested Drought

slide-26
SLIDE 26

26

SAR Images of Phenological Stages

  • R:G:B = HH:HV:VV

Harvesting Planting Vegetated2 Vegetated1 Harvested Drought

HH VV HV

slide-27
SLIDE 27

27

slide-28
SLIDE 28

28

18

Rice Phenology Mapping by Airborne SAR

Planting Vegetated$1 Vegetated$2 RIpening

slide-29
SLIDE 29

29

Class$distribution$for$selected$area

No Class Percent 1 rice150days_kb3 8.818 2 rice105days_kb5 23.071 3 rice140days_kb2 8.127 4 water 1.024 5 mixedtreecrops 16.864 6

  • thers

17.120 7 nodata 24.976 Total 100.000

Classification Result in Batola, South Kalimantan

Overall$Class$Performance$:$ (13953$/$16954$)$=$$82.3%$ Kappa$Statistic$(X100)$=$$73.4%

Courtesy:$Wahyu$Supriatna$ and$$Rizatus$Shofiya?$

slide-30
SLIDE 30

Dry$Biomass$vs$Backscadering$Coef.

30

Dry*Biomass Backsca`ering*Coef.*(dB)

n$=$37 r=0.85 r=0.74 r=0.62

HH HV VV

HH*showed*highest*correla)on*coefficient.
 ALOS/2*ScanSAR*mode*has*HH*mode.

slide-31
SLIDE 31

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 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 RADARSAT-2 Data (Ascending) RADARSAT-2 Data (Descending)

31

SAR Data Collection : Radarsat-2

Rice cultivated area mapping

Prelim inary result Processed by I NAHOR

slide-32
SLIDE 32

Vietnam: An Giang Prov.

slide-33
SLIDE 33
  • 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.

33

Mekong*Delta,*Vietnam

❑ 13*provinces*and*city; ❑ Area:*40,553*Km2*

(12%*or*1/8);

❑ Population:*17.39*M*

(20%*or*1/5);

❑ Rice*production:**

24.293*Mton**************** (>*50%*or*1/2).

******Source:*GSO,*2012

6th SAFE Workshop– Kuala Lumpur, 23rd April 2014 [Dr.*Lam*Dao*Nguyen*et*al.]

slide-34
SLIDE 34

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 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 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 6th SAFE Workshop– Kuala Lumpur, 23rd April 2014 [Dr.*Lam*Dao*Nguyen*et*al.]

slide-35
SLIDE 35

9

  • Field samples: 40
  • Check points: 40 rice

points, 24 non-rice points 6th SAFE Workshop– Kuala Lumpur, 23rd April 2014 [Dr.*Lam*Dao*Nguyen*et*al.]

slide-36
SLIDE 36

10

Cosmo SkyMed data: StripMap Pingpong: HHVV 10 dates (19 Aug – 23 Nov)

!

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

Sample rice field MPD3 sown on 17 Aug 2013 19 Aug 2013 (2) 4 Sep 2013 (18) 13 Sep 2013 (27) 19 Aug 2013 4 Sep 2013 20 Sep 2013

6th SAFE Workshop– Kuala Lumpur, 23rd April 2014 [Dr.*Lam*Dao*Nguyen*et*al.]

slide-37
SLIDE 37

15

HH/VV

HH/VV*(dB) \10.5 \7.0 \3.5 0.0 3.5 7.0 10.5 14.0 Ngày*sau*khi*sạ 10 20 30 40 50 60 70 80 90 100 110 CSK_PP_20130819 CSK_PP_20130904 CSK_PP_20130920 CSK_PP_20131006 CSK_PP_20131014 CSK_PP_20131030 CSK_PP_20131107 CSK_PP_20131115

Days after sowing

6th SAFE Workshop– Kuala Lumpur, 23rd April 2014 [Dr.*Lam*Dao*Nguyen*et*al.]

slide-38
SLIDE 38

16

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

slide-39
SLIDE 39

18

6th SAFE Workshop– Kuala Lumpur, 23rd April 2014

Commune* Name Estimated*area* (ha) Agency*data* (ha) Percentage*error* (%)

An*Bình

2,274 2,287 \0.55

Bình*Thành

2,313 2,325 \0.52

Định*Mỹ

3,178 2,965 6.70

Định*Thành

2,613 2,454 6.09

Mỹ Phú*Đông

2,676 2,740 \2.39

Phú*Thuận

1,131 835 26.19

Thoại*Giang

2,203 2,248 \2.03

Núi*Sập

455 442 2.82

Óc*Eo

522 592 \13.46

Phú*Hòa

294 310 \5.48

Vĩnh*Khánh

2,500 2,810 \12.39

Vĩnh*Phú

3,069 3,153 \2.73

Vĩnh*Trạch

1,269 1,424 \12.19

Vĩnh*Chánh

2,357 1,996 15.31

Vọng*Đông

2,601 2,096 19.42

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

Agency*data*(ha) 1,000 2,000 3,000 4,000 Es)mated*area*(ha) 1,000 2,000 3,000 4,000 y*=*0.9639x*+*18.956 R²*=*0.9511

[Dr.*Lam*Dao*Nguyen*et*al.]

slide-40
SLIDE 40

25

YRa = 6.4575 - 2.6085 *Ra1+ 0.9477 *Ra2 - 0.5232 *Ra3 r2 = 0.802, sey = 0.282 ton/ha Where: YRa: rice yield (ton/ha), Ra1 : polarisation ratio of first date image,

Ra2 : polarisation ratio of second date image, Ra3 : polarisation ratio of third date image, r2 : the coefficient of determination, sey : 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

6th SAFE Workshop– Kuala Lumpur, 23rd April 2014 [Dr.*Lam*Day*Nguyen*et*al.]

slide-41
SLIDE 41

ALOS\2*:*New*SAR*Satellite

41

ALOS/2*satellite*

*Launch** :*24th*May*2014* *Orbit*type** :*Sun\synchronous Altitude** :*628*km*+/\*500*m*(for*reference*orbit) *Revisit*time** :*14*days LSDN* * :*12:00*+/\*15*min

SAR*antenna

X Y Z

Solar*paddles 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.

12: 05 p.m. on May 24

slide-42
SLIDE 42

ALOS\2*Observaton*Mode

42

Spotlight Ultra Fine High sensitive Fine ScanSAR nominal ScanSAR wide Bandwidth 84MHz 84MHz 42MHz 28MHz 14MHz 28MHz 14MHz Resolution Rg×Az 3×1m 3m 6m 10m 100m 60m Swath Rg×Az 25×25km 50km 50km 70km 350km (5-scan) 490km (7-scan) Polarization SP SP/DP SP/DP/FP/CP SP/DP NESZ

  • 24dB
  • 24dB
  • 28dB
  • 26dB
  • 26dB
  • 23dB
  • 23dB

S/A Rg 25dB 25dB 23dB 25dB 25dB 20dB Az 20dB 25dB 20dB 23dB 20dB 20dB

SP : HH or VV or HV , 
 DP : HH+HV or VV+VH , 
 FP : HH+HV+VH+VV ,
 CP : Compact pol

FP
 1 cov/yr
 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.

! ALOS*has*several*observaton*modes*different*in*polarizaton,** spatal*and*actual*temporal*resolutons.*

HH+HV
 6 cov/yr
 (SE Asia)

slide-43
SLIDE 43

Summary

! Progress* ! 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

slide-44
SLIDE 44

Thank You for Your Attention !

!

  • hyoshi.kei@jaxa.jp