PolSARpro for Agriculture Forest

Transcription

1 PolSARpro for Agriculture Forest Eric POTTIER Wednesday 5 September, Lecture D3PB-Part 3 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 1

2 EMISAR DCRS (DK) G3 Aircraft L, C-Band (Quad) 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 2

3 Foulum (L,C-Band) - Land Cover Monitoring LandSat TM Classification 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 3

4 MAIN MENU PolSARpro Full Software Single Data Set Tutorial on POLSAR and POLinSAR Help Files Viewer Display Tools Version for the EO Scientific Investigator Spaceborne Sensors: ALOS, ENVISAT, RADARSAT2, TerraSar, SIR-C Airborne Sensors: AIRSAR, Convair, EMISAR, ESAR, PISAR, RAMSES 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 4

5 PROCESSING CHAIN Configuration Data Import Input Data File: Raw Binary Data Airborne Sensors QuickLook Extract Raw Data Data Process Polarimetric Elements Speckle Filetring Data Analysis Statistics / Histograms / Profiles Polarimetric Segmentation H / A / alpha classification Wishart - H / A / alpha classification Supervised Wishart classification Rules-Based Hierarchical Classification 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 5

6 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 6

7 Processing Configuration Configure Data Main Directory location Input Directory: C:/EMISAR/Foulum_L Note: Copy the data directory to the root directory to shorten the path name length 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 7

8 Data IMPORT 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 8

9 Data IMPORT EMISAR Data File Do it Yourself: Data Set: EMISAR / FOULUM_L_Format_S2 Input Data File configuration: 1750 rows x 1000 col Convert Input IEEE binary Format 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 9

10 Data IMPORT Quick Look Create a QuickLook (24-bit RGB Output BMP File) color image: Pauli Color coded Image T22 = HH-VV T33 = HV+VH T11 = HH+VV Do it Yourself: Run Quick look 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 10

11 Data IMPORT Extract Binary Data Convert Fully Polarimetric Data Files to different Polarimetric Output Format INPUT DIR OUTPUT DIR Fully Polarimetric Data Files config.txt s11.bin, s12.bin, s21.bin, s22.bin 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 11

12 Data IMPORT Extract Binary Data Do it Yourself: Select an area, then Extract the Data : Full Resolution -> [S2] View the Pauli-RGB File automatically created 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 12

13 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 13

14 Data PROCESS [S2] Elements MENU 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 14

15 Data PROCESS [S2] Elements MENU 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 15

16 Data PROCESS [S2] Elements MENU DATADIR 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 Do it Yourself: Select some elements, set the parameters and view the corresponding BMP files (select BMP). 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 16

17 Data PROCESS [S2] Elements SPECKLE FILTER 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 17

18 Data PROCESS [S2] Elements SPECKLE FILTER DATADIR config.txt s11.bin, s12.bin s21.bin, s22.bin Lee_Refined_Filter BoxCar_Filter DATADIR_BOX DATADIR_LEE T3 config.txt [T3x3] Elements Do it Yourself: Apply a Box Speckle Filter then Create a RGB File ([T3]) (Don tforgetto change the Main directory) Apply a Lee Speckle Filter then Create a RGB File ([T3]) 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 18

19 Data PROCESS [S2] Elements DATA ANALYSIS - STATISTICS 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 19

20 Data PROCESS [S2] Elements DATA ANALYSIS - STATISTICS Channel: Polarimetric Elements [S2] Elements: Real part, Imaginary part, Amplitude, Phase Theoretical PDF: Gaussian, Exponential, Rayleigh, Uniform 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 20

21 Data PROCESS [S2] Elements DATA ANALYSIS - HISTOGRAMS 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 21

22 Data PROCESS [S2] Elements DATA ANALYSIS - HISTOGRAMS 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 22

23 Data PROCESS [S2] Elements DATA ANALYSIS - PROFILES 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 23

24 Data PROCESS [S2] Elements DATA ANALYSIS - PROFILES X / Y Range : 200 pix XY Range (3D) : 30 pix 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 24

25 Data PROCESS [S2] Elements POLARIMETRIC SEGMENTATION 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 25

26 POLARIMETRIC SEGMENTATION H / A / Alpha Classification DATADIR config.txt s11.bin, s12.bin s21.bin, s22.bin entropy.bin, anisotropy.bin, alpha.bin combination_ha.bin, combination_1mha.bin, combination_h1ma.bin, combination_1mh1ma.bin H_A_class.bin, H_Alpha_class.bin, A_Alpha_class.bin entropy.bmp, anisotropy.bmp, alpha.bmp combination_ha.bmp, combination_1mha.bmp, combination_h1ma.bmp, combination_1mh1ma.bmp H_A_class.bmp, H_Alpha_class.bmp, A_Alpha_class.bmp H_A_occurence.bmp, H_Alpha_ occurence.bmp, A_Alpha_ occurence.bmp, H_A_segmented.bmp, H_Alpha_ segmented.bmp, A_Alpha_ segmented.bmp X_Y_occurrence_plane X_Y_segmented_plane HAlphaLambda_RGB.bmp, HAAlpha_RGB.bmp HACombinations_RGB.bmp Do it Yourself: Select some elements, set the parameters (Nwin = 1)and view 05/09/07 Lecture D3PB-Part 3 the corresponding BMP files. PolSARpro for Agriculture Forest Eric POTTIER 26

27 POLARIMETRIC SEGMENTATION H / A / Alpha Classification 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 27

28 POLARIMETRIC SEGMENTATION H / A / Alpha Classification 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 28

29 POLARIMETRIC SEGMENTATION H / A / Alpha Classification X_Y_occurrence_plane X_Y_segmented_plane 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 29

30 POLARIMETRIC SEGMENTATION H / A / Alpha Classification DATADIR config.txt s11.bin, s12.bin s21.bin, s22.bin H_alpha_lambda_class1(2,3).bin, H_alpha_lambda_occurence_class1(2,3).bin, H_alpha_lambda_segmented_class1(2,3).bin, H_alpha_lambda_class.bin, H_alpha_lambda_class1(2,3).bmp, H_alpha_lambda_occurence_class1(2,3).bmp, H_alpha_lambda_segmented_class1(2,3).bmp, H_alpha_lambda_class.bmp, X_Y_occurrence_plane X_Y_segmented_plane Do it Yourself: Select some elements, set the parameters (Nwin = 1)and view the corresponding BMP files. 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 30

31 POLARIMETRIC SEGMENTATION H / A / Alpha Classification 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 31

32 Data PROCESS [S2] Elements POLARIMETRIC SEGMENTATION 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 32

33 POLARIMETRIC SEGMENTATION H / A / Alpha - Wishart Classification DATADIR config.txt s11.bin, s12.bin s21.bin, s22.bin Wishart_H_alpha_class_X.bin Wishart_H_A_alpha_class_X.bin Wishart_H_alpha_class_X.bmp Wishart_H_A_alpha_class_X.bmp X = window size Do it Yourself: Set the parameters, run and view the corresponding BMP files. 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 33

34 Data PROCESS [S2] Elements POLARIMETRIC SEGMENTATION 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 34

35 POLARIMETRIC SEGMENTATION Supervised Wishart Classification 1: Select the Output Directory 2 : Select the Output Image Number of Rows and Columns 3 : Select BMP creation files, Reject Class and Confusion Matrix. Note: BMP and Confusion Matrix are selected by default 4 : Enter the Window Size (Nwin = 5 ) and the Reject Ratio Coefficient.(if Reject Class option is selected) 5 : Edit the ColorMap and choose the colours associated to the different classes 6 : Select (or not) the Coded ColorMap option 7 : Enter the name of the Training Areas list text file The default output file name is set to: Config / training_areas.txt 8 : Open Graphic Editor to define graphically the Areas of Interest (AoI) or run Text Editor to enter the coordinates of the Areas of Interest (AoI)in the Training Areas text file 9 : Run Training Process. This program will define the training clusters centres from the selected Areas of Interest (AoI). 10 : Run the Wishart Supervised Classification procedure 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 35

36 POLARIMETRIC SEGMENTATION Supervised Wishart Classification 1 : Create a new Training Class and add it to the list. 2 : Delete the current Training Class from the list. 3 : Create a new Training Area and add it to the list of the current Training Class. 4 : Delete the current Training Area from the list. 5 : Clear the different contours on the chosen image. 6 : Delete all the Training Classes and associated Training Areas from the list. 7 : Save the Training Class and associated Training Areas list....the default output file name is set to MD / training_areas.txt. 8 : Rectangular selection of Area of Interest (AoI). 9 : Polygonal selection of Area of Interest (AoI).... Note: The contour is automatically closed by clicking on the Mouse Right Button. 10 : Toggle selected area contour color (black / white). Class : Redraw all the Training Areas of the current Training Class. Class : Move in the up/down direction in the Training Class list. Area : Redraw the current Training Area of the current Training Class. Area : Move in the up/down direction in the Training Area list of the current Training Class. 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 36

37 POLARIMETRIC SEGMENTATION Supervised Wishart Classification 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 37

38 POLARIMETRIC SEGMENTATION Supervised Wishart Classification DATADIR config.txt s11.bin, s12.bin s21.bin, s22.bin Run Training Process Training_areas.txt Training_areas.bin Training_cluster_set.bmp Run Classification Do it Yourself: Set the parameters, select different classes, run and view the corresponding BMP files. Supervised_class_X.bin Supervised_class_rej_X.bin Confusion_matrix_X.txt Confusion_matrix_rej_X.txt Supervised_class_X.bmp Supervised_class_rej_X.bmp Classified_cluster_set.bmp Classified_cluster_set_rej.bmp X = window size 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 38

39 POLARIMETRIC SEGMENTATION Supervised Wishart Classification 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 39

40 POLARIMETRIC SEGMENTATION Supervised Wishart Classification 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 40

41 POLARIMETRIC SEGMENTATION Supervised Wishart Classification 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 41

42 Data PROCESS [S2] Elements POLARIMETRIC SEGMENTATION 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 42

43 POLARIMETRIC SEGMENTATION Rule-Based Hierarchical Classification NOTE: Change Input Directory to: ~/EMISAR/Foulum_C 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 43

44 POLARIMETRIC SEGMENTATION Rule-Based Hierarchical Classification 1 : New 2 : Enter Parameter Label 3 : Select the Parameter File 4 : Enter 5 : go to step 1 6 : Save & Exit Parameter Label Parameter File sigma0_hv_l_april ~/EMISAR/Foulum_C/sigma0_hv_L_April.bin sigma0_hh_c_april ~/EMISAR/Foulum_C/sigma0_hh_C_April.bin rho_hhvv_c_april ~/EMISAR/Foulum_C/rho_hhvv_C_April.bin sigma0_hv_c_april ~/EMISAR/Foulum_C/sigma0_hv_C_April.bin sigma0_rr_ov_rl_c_july ~/EMISAR/Foulum_C/sigma0_rr_Ov_sigma0_rl_C_July.bin sigma0_hv_c_june ~/EMISAR/Foulum_C/sigma0_hv_C_June.bin rho_hhvv_c_june ~/EMISAR/Foulum_C/rho_hhvv_C_June.bin 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 44

45 POLARIMETRIC SEGMENTATION Rule-Based Hierarchical Classification Node Class 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 45

46 POLARIMETRIC SEGMENTATION Rule-Based Hierarchical Classification Node Class Active Node Class Definition 1 : Enter Class Number 2 : Enter Class Definition 1 : Select Para 1 2 : Select Para 2 3 : Enter Coeff 1 4 : Enter Coeff 2 5 : Select Operator 6 : Enter Threshold 7 : Enter Para1 x Coeff1 + Para2 x Coeff2 < Threshold 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 46

47 POLARIMETRIC SEGMENTATION Rule-Based Hierarchical Classification sigma0_hv_l_april > -18dB 8 sigma0_hh_c_april > -30dB rho_hhvv_c_april 0.02 x sigma0_hv_c_april > 1 9 sigma0_hv_c_june > -18dB sigma0_hv_c_june > -12dB 7 sigma0_rr_ov_rl_c_july > sigma0_hv_c_june > -30dB rho_hhvv_c_june > /09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 47

48 POLARIMETRIC SEGMENTATION Rule-Based Hierarchical Classification 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 48

49 Questions? 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 49

50 Single Data Set Package POLinSAR Tutorial Training Course 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 50

51 Direct access to the Tutorial while using PolSARpro facilities The Tutorial is made available in PDF format. 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 51

52 Direct access to the Tutorial while using PolSARpro facilities The Tutorial is made available in PDF format. 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 52

53 MAIN MENU PolSARpro Full Software Multi Data Sets Dual POLinSAR Tutorial on POLSAR and POLinSAR Help Files Viewer Display Tools Version for the EO Scientific Investigator Spaceborne Sensors: ALOS, ENVISAT, RADARSAT2, TerraSar, SIR-C Airborne Sensors: AIRSAR, Convair, EMISAR, ESAR, PISAR, RAMSES 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 53

54 Processing Configuration Configure Data Main Input Master Directory location Configure Data Main Input Slave Directory location Input Master Directory: C:/POLinSAR_Training_Course/Master_Track Input Slave Directory: C:/POLinSAR_Training_Course/Slave_Track Note: Copy the data directory to the root directory to shorten the path name length 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 54

55 MENU 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 55

56 ELEMENTS 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 56

57 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 Do it Yourself: Select some elements, set the parameters and view the corresponding BMP files (select BMP). 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 57

58 RAW INTERFEROGRAM 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 58

59 RAW INTERFEROGRAM DATA_MASTERDIR_SLAVEDIR config.txt interferogram_xx_xx.bin interferogram_xx_xx.bmp 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 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 59

60 FLAT EARTH REMOVAL 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 60

61 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 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 61

62 RAW INTERFEROGRAM DATA_MASTERDIR_SLAVEDIR_FER config.txt interferogram_xx_xx.bin interferogram_xx_xx.bmp 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 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 62

63 COHERENCE ESTIMATION 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 63

64 COHERENCE ESTIMATION DATA_MASTERDIR_SLAVEDIR_FER config.txt cmplx_coh_xx.bin cmplx_coh_xx_mod.bmp cmplx_coh_xx_pha.bmp Do it Yourself: Select polarization channels (linear, circular, pauli), set the parameters (Box Car = 11x11) and view the corresponding BMP files (select BMP). 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 64

65 HEIGHT ESTIMATION INVERSION PROCEDURES 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 65

66 HEIGHT ESTIMATION INVERSION PROCEDURES INVERSION PROCEDURES DEM Differencing Algorithm Coherence Amplitude Inversion Procedure Ground Phase Estimation RVOG Inversion Procedure 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 66

67 VOLUME COHERENCE MODEL Modeling z Parameter Estimation z h v h v 0 0 Simplifications : Only 2 significant mechanisms Low density medium No refraction 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 67

68 VOLUME COHERENCE MODEL φ0 f ( z ) γ VOL = j 0 e φ Topographic Phase Vertical Structure function h v 0 f ( z )e h v 0 jk z z f ( z )dz k dz z = 4πΔθ λ sin( θ Vertical Wavenumber 0 ) 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 68

69 VOLUME COHERENCE MODEL Case of Uniform Random Layer f ( z ) = e σ z cos( θ 0 ) θ 0 σ Incidence Angle Extinction Coefficient 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 69

70 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 70 VOLUME COHERENCE MODEL 1 e 1 e e p p e 1 e 1 e p p e v v z v 0 v v 1 0 ph h jk ph 1 j ph h p 1 j VOL = = φ φ γ z 1 0 jk p p ) cos( p + = = θ σ ) sin( 4 k 0 z θ λ πδθ = POLARIZATION INDEPENDENT

71 RVOG COHERENCE MODEL (Random Volume Over Ground) 2 Layer Combined Surface and random Volume Scattering μ( w) = γ ( w) jφ VOL e 0 γ + = 1 + μμ ( w) ( w) Surface Scattering Contribution Volume Scattering Contribution B. Treuhaft (2000), S.R. Cloude (2003) POLARIZATION DEPENDENT G / V ratio 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 71

72 FOREST HEIGHT ESTIMATION w v Polarisation Channel corresponding to Volume Scattering γ jφ0 ( w v ) a = e γ VOL μa0 2HV w s Polarisation Channel corresponding to Surface Scattering γ ( w ) s = e jφ γ μ μ( w s ) ( w ) 0 VOL j s a μa e φ 0 HH-VV 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 72

73 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 73 FOREST HEIGHT ESTIMATION ( ) ( ) ( ) ( ) [ ] ( ) ( ) = + + = = s v 1 v 0 s s VOL j s VOL j v w w Model h w 1 w e w e w 0 0 γ γ μ σ φ μ μ γ γ γ γ φ φ a 4 Parameters 4 Observables INVERSION

74 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 74 FOREST HEIGHT ESTIMATION ( ) ( ) ( ) ( ) ( ) v z 0 0 h jk s VOL s v j s VOL j v e w w w e w e w α γ γ γ γ γ γ γ φ φ = = a a DEM Differencing Algorithm [ ] [ ] z s v v k ) w ( arg ) w ( arg h γ γ

75 FOREST HEIGHT ESTIMATION DEM Differencing Algorithm 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 75

76 FOREST HEIGHT ESTIMATION Coherence Amplitude Inversion Procedure Assumption: Only Volume Scattering is present γ jφ0 ( w v ) = e γ VOL a γ ( w v ) = γ VOL min h v γ ( w v ) p p 1 e e p 1 h ph v v D Search Procedure Look Up Table (LUT) 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 76

77 FOREST HEIGHT ESTIMATION Coherence Amplitude Inversion Procedure 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 77

78 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 78 FOREST HEIGHT ESTIMATION ( ) ( ) ( ) ( ) ( ) ( )( ) ( ) ( )( ) [ ] 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 = = + + = = γ γ φ γ γ μ μ γ γ γ γ φ φ φ a 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 γ γ γ γ γ γ = R = = 2 A 4 AC B B L 2 =

79 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 0.3 ε 0.5 Suitable Compromise DEM Differencing Inversion Coherence Amplitude Inversion 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 79

80 FOREST HEIGHT ESTIMATION RVOG Inversion Procedure 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 80

81 HEIGHT ESTIMATION INVERSION PROCEDURES 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 = 4) and view the corresponding BMP files. HH-VV 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 81

82 HEIGHT ESTIMATION INVERSION PROCEDURES DEM_diff_heights Coh_heights Ground_phase Ground_phase_median RVOG_phase_heights RVOG_heights 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 82

83 PHASE CENTER HEIGHT HISTOGRAMS Do it Yourself: Select a BMP file Select a BIN file Select Input Data Format Select Show Select Area SAVE PLOT 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 83

84 HEIGHT ESTIMATION INVERSION PROCEDURES DEM_diff_heights Coh_heights RVOG_heights 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 84

85 DATA ANALYSIS DATA PROFILES 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 85

86 DATA ANALYSIS DATA PROFILES 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 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 86

87 DATA ANALYSIS DATA PROFILES DEM_diff_heights Coh_heights RVOG_heights 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 87

88 Questions? 05/09/07 Lecture D3PB-Part 3 PolSARpro for Agriculture Forest Eric POTTIER 88