Forest Retrievals. using SAR Polarimetry. (Practical Session D3P2a)

Transcription

1 Forest Retrievals using SAR Polarimetry (Practical Session D3P2a) Laurent FERRO-FAMIL - Eric POTTIER University of Rennes 1

2 Pol-InSAR Practical Forest Application

3

4 PolSARpro SIM PolSARproSim is a rapid, coherent, fully polarimetric and interferometric SAR simulation of forest.

5 PolSARpro SIM The SAR image is evaluated as a coherent sum of scattering events from small elements of the scene Direct-Ground, Direct-Volume and Ground-Volume contributions are included, with both trees and short vegetation comprising Volume terms. DECIDUOUS GV DG DV PINE RANDOM HEDGE Given the map of tree locations and dimensions a grid of points is used to sample the attenuation of the coherent wave in 3D

6 B = 10m PolSARpro SIM Pol-InSAR Data = 45 o h=3km r 2 r 1 h v =10m y Bragg Surface Scattering Geometric configuration Platform altitude : 3000m Incidence angle: 45 Horizontal Baseline : 10m Vertical Baseline : 0m System Configuration Frequency : 1.5 GHz Azimuth resolution : m Range resolution : m Ground Surface Configuration Surface properties : 0 (smoothest) Ground moisture Content : 0 (driest) Azimuth / Range ground slope : 0 % Forest configuration Tree Species : 0 (hedge) Tree Height: 10m Forest stand density : 0.2 Forest Stand Circular Area : 1 Ha

7 PolSARpro SIM DATA_MASTERDIR config.txt s11.bin, s12.bin s21.bin, s22.bin DATA_SLAVEDIR config.txt s11.bin, s12.bin s21.bin, s22.bin flat_earth.bin kz.bin

8 MAIN MENU

9 ENVIRONMENT Configure Data Main Directories location Input Master Directory: C:/Pol-InSAR_Training_Course/Master_Track Input Slave Directory: C:/Pol-InSAR_Training_Course/Slave_Track

10 PROCESS DATA

11 PROCESS DATA

12 ELEMENTS Do it Yourself: Select some elements, set the parameters and view the corresponding BMP files (select BMP).

13 ELEMENTS DATA_MASTERDIR config.txt s11.bin, s12.bin s21.bin, s22.bin Axy.bin, Ixy.bin Ixy_db.bin sxy_pha.bin Axy.bmp, Ixy.bmp Ixy_db.bmp sxy_pha.bmp

14 PROCESS DATA

15 RAW INTERFEROGRAM Do it Yourself: Select polarization channels, set the parameters and view the corresponding BMP files. Note: The Output Directory is automatically set to: MasterDir_SlaveDir

16 RAW INTERFEROGRAM DATA_MASTERDIR_SLAVEDIR config.txt interferogram_xx_xx.bin interferogram_xx_xx.bmp

17 PROCESS DATA

18 FLAT EARTH REMOVAL DATA_SLAVEDIR config.txt s11.bin, s12.bin s21.bin, s22.bin DATA_SLAVEDIR_FER config.txt s11.bin, s12.bin s21.bin, s22.bin Do it Yourself: Enter Flat Earth file name, set the parameters and run the function. Note: The Input Slave Directory is automatically set to: SlaveDir_FER

19 RAW INTERFEROGRAM Do it Yourself: Select polarization channels, set the parameters and view the corresponding BMP files. Note: The Output Directory is automatically set to: MasterDir_SlaveDir_FER

20 RAW INTERFEROGRAM DATA_MASTERDIR_SLAVEDIR_FER config.txt interferogram_xx_xx.bin interferogram_xx_xx.bmp

21 PROCESS DATA

22 COHERENCE ESTIMATION Do it Yourself: Select polarization channels (linear, circular, pauli), set the parameters (Box Car = 11x11) and view the corresponding BMP files (select BMP).

23 COHERENCE ESTIMATION DATA_MASTERDIR_SLAVEDIR_FER config.txt cmplx_coh_xx.bin cmplx_coh_xx_mod.bmp cmplx_coh_xx_pha.bmp

24 PROCESS DATA

25 HEIGHT ESTIMATION

26 HEIGHT ESTIMATION INVERSION PROCEDURES DEM Differencing Algorithm Coherence Amplitude Inversion Procedure Ground Phase Estimation RVOG Inversion Procedure

27 VOLUME COHERENCE MODEL Modeling z Parameter Estimation z h v h v 0 0 Simplifications : Only 2 significant mechanisms Low density medium No refraction

28 VOLUME COHERENCE MODEL VOL j 0 e h v 0 h f(ze) v 0 jkz z dz f(z) dz 0 k z Topographic Phase 4 sin( ) 0 Vertical Wavenumber POLARIZATION INDEPENDENT

29 VOLUME COHERENCE MODEL VOL j 0 e h v 0 h f(ze) v 0 jkz z dz f(z) dz Vertical Structure function z cos( ) f(z) e 0 Case of Uniform Random Layer 0 Incidence Angle Extinction Coefficient POLARIZATION INDEPENDENT

30 RVOG COHERENCE MODEL RVOG = Random Volume Over Ground 2 Layer Combined Surface and random Volume Scattering w 1w w w e j 0 VOL Surface Scattering Contribution Volume Scattering Contribution B. Treuhaft (2000), S.R. Cloude (2003) POLARIZATION DEPENDENT G / V ratio

31 FOREST HEIGHT ESTIMATION w v Polarisation Channel corresponding to Volume Scattering j v e 0 VOL 0 w 2HV w s Polarisation Channel corresponding to Surface Scattering j we s ws e 0 1w 0VOL j s HH-VV

32 FOREST HEIGHT ESTIMATION 0 j0 we v VOL h v 1 w v j VOL w s Model 0 we s w s 1w s 4 Parameters 4 Observables INVERSION

33 v z 0 0 h jk s VOL s v j s VOL j v e w w w e w e w DEM Differencing Algorithm z s v v k ) w ( arg ) w ( arg h FOREST HEIGHT ESTIMATION

34 FOREST HEIGHT ESTIMATION Coherence Amplitude Inversion Procedure Assumption: Only Volume Scattering is present j0 wv e VOL wv VOL min h v ( w v ) p p 1 e e p 1 h ph v v D Search Procedure Look Up Table (LUT)

35 L 1 w w arg ˆ L L 1 w w e w 1 w e w e w v s 0 v s j s s VOL j s VOL j v Topographic Phase Estimation s s w 1 w L With: Estimation of L 0 C BL AL 1 L L 1 w w 2 2 v s 2 v s s * v s 2 v w w C w w w 2 B 1 w A 2A 4 AC B B L 2 FOREST HEIGHT ESTIMATION

36 FOREST HEIGHT ESTIMATION RVOG Inversion Procedure min h v, ( w v ) e j ˆ 0 p p 1 e e p 1 h ph v v 1 1 Expensive 2-D Search Procedure! h v arg ˆ 1 ( w ) 2 sinc ( w ) k z v 0 k z v Suitable Compromise DEM Differencing Inversion Coherence Amplitude Inversion

37 HEIGHT ESTIMATION DATA_MASTERDIR_SLAVEDIR_FER config.txt DEM_diff_heights.bin, Coh_heights.bin Ground_phase.bin, Ground_phase_median.bin RVOG_phase_heights.bin, RVOG_heights.bin DEM_diff_heights.bmp, Coh_heights.bmp Ground_phase.bmp, Ground_phase_median.bmp RVOG_phase_heights.bmp, RVOG_heights.bmp 2HV Do it Yourself: Set the parameters (Median Size = 21, Factor = 0.4) and view the corresponding BMP files. HH-VV

38 HEIGHT ESTIMATION DEM_diff_heights Coh_heights Ground_phase Ground_phase_median RVOG_phase_heights RVOG_heights

39 PROCESS DATA

40 Do it Yourself: Select a BMP file Select a BIN file Select Input Data Format Select Show Select Area (line or rect) SAVE PLOT HEIGHT ESTIMATION

41 HEIGHT ESTIMATION DEM_diff_heights Coh_heights RVOG_heights

42 PROCESS DATA

43 Do it Yourself: Select a BMP file Select a BIN file Select Input Data Format Select Pixel Select Show Select Representation X Range / Y Range = 200pix XY Range = 30 pix (3D) Set Min / Max Values PLOT HEIGHT ESTIMATION

44 HEIGHT ESTIMATION DEM_diff_heights Coh_heights RVOG_heights

45 Pol-InSAR Practical Forest Application

46 MAIN MENU

47 ENVIRONMENT Configure Data Main Directories location Input Master Directory: C:/Taunstein_ESAR/master_slc Input Slave Directory: C:/Traunstein_ESAR/slave_slc

48 PROCESS DATA

49 PROCESS DATA

50 RAW INTERFEROGRAM Do it Yourself: Select polarization channels, set the parameters and view the corresponding BMP files. Note: The Output Directory is automatically set to: MasterDir_SlaveDir

51 RAW INTERFEROGRAM DATA_MASTERDIR_SLAVEDIR config.txt interferogram_xx_xx.bin interferogram_xx_xx.bmp

52 PROCESS DATA

53 FLAT EARTH REMOVAL DATA_SLAVEDIR config.txt s11.bin, s12.bin s21.bin, s22.bin DATA_SLAVEDIR_FER config.txt s11.bin, s12.bin s21.bin, s22.bin Do it Yourself: Enter Flat Earth file name, set the parameters and run the function. Note: The Input Slave Directory is automatically set to: SlaveDir_FER

54 RAW INTERFEROGRAM Do it Yourself: Select polarization channels, set the parameters and view the corresponding BMP files. Note: The Output Directory is automatically set to: MasterDir_SlaveDir_FER

55 RAW INTERFEROGRAM DATA_MASTERDIR_SLAVEDIR_FER config.txt interferogram_xx_xx.bin interferogram_xx_xx.bmp

56 Convert S2 - T6 : Multilook

57 Convert S2 - T6 : Multilook Do it Yourself: Select Multi Look : Row = 6 and Col = 2 Select Output Data Format : 2 x [S2] >> [T6] Note: The Output Directory is automatically set to: MasterDir_SlaveDir_FER_MLK

58 PROCESS DATA

59 PROCESS DATA Do it Yourself: Set the parameters : Num Looks = 3 ; Window Size = 3. Note: The Output Directory is automatically set to: MasterDir_SlaveDir_FER_MLK_LEE

60 Create RGB file PROCESS DATA

61 PROCESS DATA

62 PROCESS DATA

63 PROCESS DATA Do it Yourself: Select the correlation coefficients, set the parameters (Box Car= 5x5) and view the corresponding BMP files.

64 PROCESS DATA

65 PROCESS DATA

66 PROCESS DATA

67 COHERENCE ESTIMATION Do it Yourself: Select polarization channels (linear, circular, pauli), set the parameters (Box Car = 7x7) and view the corresponding BMP files (select BMP).

68 COHERENCE ESTIMATION

69 COHERENCE ESTIMATION

70 PROCESS DATA

71 HEIGHT ESTIMATION

72 HEIGHT ESTIMATION 2HV HH-VV Do it Yourself: Set the parameters (Median Size = 21, Factor = 0.4) and view the corresponding BMP files. 2D Kz File : DataDirectory / kz-ph-mlk / kz.bin

73 HEIGHT ESTIMATION

74 HEIGHT ESTIMATION Min = 0m Max = 50m

75 HEIGHT ESTIMATION Data Analysis : Histogram

76 HEIGHT ESTIMATION Do it Yourself: Select a BMP file Select a BIN file Select Input Data Format Select Show Select Area SAVE PLOT

77 PROCESS DATA

78 Do it Yourself: Select a BMP file Select a BIN file Select Input Data Format Select Pixel Select Show Select Representation X Range / Y Range = 30pix XY Range = 30 pix (3D) Set Min / Max Values PLOT HEIGHT ESTIMATION

79 Questions?