1 / 43

BARE EARTH DEM FROM SRTM DATA

BARE EARTH DEM FROM SRTM DATA. CG. Objective of the project. Improve Digital Elevation Models (DEMs) from SRTM mission for hydrodynamic modeling and other applications. Obtain spatial and temporal structure of vegetation biomass.

auryon
Télécharger la présentation

BARE EARTH DEM FROM SRTM DATA

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. BARE EARTH DEM FROM SRTM DATA CG

  2. Objective of the project • Improve Digital Elevation Models (DEMs) from SRTM mission for hydrodynamic modeling and other applications. • Obtain spatial and temporal structure of vegetation biomass. • Obtain a high resolution bare-earth DEM (90m) for areas where the elevation data is not available.

  3. LiDAR 1m vegetation appears as point clouds

  4. Sensor dependence and data smoothing: LiDAR at 30 m

  5. Sensor dependence: Vegetation appears as uniformly distributed noise at some scales at SRTM 30 m NED data @30 m SRTM Data @ 30 m

  6. SRTM VEGETATION RESPONSE • A Fundamental law of science applies to the SRTM data • “ One scientist’s noise is another scientist’s signal”

  7. Fourier Analysis • Vege. Topo. = bare earth Topo.+ veg. ht Unknown, but at large scale bare-earth topo. Amplitude is approximately equals to the amplitude of the SRTM signal. Phase is unknown and will be resolved iteratively Unknown, but we have shown that phase of the veg. height can be obtained using SRTM signal. Amplitude can will be resolved iteratively Known from SRTM Signal We have two equations in the Fourier domain and two unknowns- can be resolved. Can be extended to the 2D

  8. Spatial domain frequency domain Real part Imaginary part

  9. Fourier based approach • Linearity of the Fourier transform is used. • If: x1[n] + x2[n] = x3[n], • then: ReX1[f] + ReX2[f] = ReX3[f] and ImX1[f] + ImX2[f] = ImX3[f].

  10. Additively in Spectral Domain: proposed model Fourier Amplitude Fourier Phase

  11. model…. For bare-earth

  12. Study Area: SF Eel River

  13. Fourier Power spectrum : vegetation topography, bare-earth, and vegetation height- SF Eel LiDAR data

  14. Study area 2

  15. Study area 3: Tenderfoot Creek Vertical profiles Vegetation fractional cover ~ 90

  16. Tendorfoot Creek:

  17. Vegetation effect on SRTM:

  18. SRTM and ASTER

  19. A Simple approach: • Low frequency content of the vegetation topography is coming from bare-earth topography • High frequency content of the vegetation topography is coming from the vegetation height. • Two regions are separated by a scale break at frequency fc

  20. Approach…. low frequency content~ bare-earth topo. High frequency content~ vegetation height

  21. Model parameters-Amplitudes of the bare-earth and vegetation topography: LiDAR data at low frequencies

  22. Model parameters-phase of the bare-earth and vegetation topography: at low frequencies

  23. Canopy topography phase VS Canopy height phase CANOPY HEIGHT PHASE CAN BE OBTAINED! AMPLITUDE CAN BE RESOLVED ITERATIVELY

  24. Back to the model equations.. For bare-earth = and = The low frequency content is obtained using model eqn. The high frequency information of the bare-earth is obtained by Fourier interpolations

  25. model…. For vegetation height = and = High frequency contact is obtained from the model eqn. The low frequency information of the bare-earth is obtained by Fourier interpolations

  26. Results..

  27. Model vegetation height

  28. Vegetation height: model vs. LiDAR

  29. Results: A comparison of Model topo. and LiDAR bare-earth: It preserves the second order moments

  30. Conclusions: • Vegetation topography leaves a clear statistical signature about the vegetation height and the bare-earth ( can be extracted using minimally used ancillary data.) • A clear signature of the vegetation height data – two power law scaling regimes (i.e., low frequency and high frequency) with a scaling break at a intermediate characteristic frequency. • The characteristic frequency depends on the vegetation density and the grid resolution • vegetation mainly distort bare-earth power law scaling regime at the high frequency range • If the objective is the canopy model (hydrodynamic model) model resolution should be selected from low frequency (high frequency) range.

  31. Can we improve our solution! • Fourier approach can not localized spatially! • Solution obtained from iterative optimization may not be always realistic. • Can we integrate remotely sensed ground truth data to obtain realistic solutions • We want to zoom down to vegetation patch scales and implement scale dependant interpolating approach to remove the vegetation effect/ obtain the vegetation height while preserving observed statistical properties of the bare earth

  32. SCALE 1 SACLE 2 SCALE 3 SCALE 4 SCALE 5 Chandana Gangodagamage

  33. FUSION OF DIFFERENT SENSOR DATA ASTER and SRTM

  34. Comparison of ASTER-SRTM-LiDAR data elevation

  35. TfTgTvAfAgAv eq1= 'Af/(1+(tan(Tf))^2)=Ag/(1+(tan(Tg))^2)+Av/(1+(tan(Tv))^2)' eq2= 'Af/(1+(tan(Tf))^2)*tan(Tf)=Ag/(1+(tan(Tg))^2)*tan(Tg)+Av/(1+(tan(Tv))^2)*tan(Tv)' [Tg Av]=solve(eq1eq2,Tg,Av) • [Tg Av]=solve(eq1, eq2,Tg, Av) • Tg= • (tan(Tf)+tan(Tf)*tan(Tv)^2-1/2/(-Af*tan(Tv)+Af*tan(Tf))*(Ag+Ag*tan(Tf)^2+(Ag^2+2*Ag^2*tan(Tf)^2+Ag^2*tan(Tf)^4+4*Af*tan(Tv)^2*Ag*tan(Tf)^2-4*Af^2*tan(Tv)^2+4*Af*tan(Tv)^2*Ag+8*Af^2*tan(Tv)*tan(Tf)-4*Af*tan(Tf)^3*Ag*tan(Tv)-4*Af*tan(Tf)*Ag*tan(Tv)-4*Af^2*tan(Tf)^2)^(1/2))-1/2/(-Af*tan(Tv)+Af*tan(Tf))*(Ag+Ag*tan(Tf)^2+(Ag^2+2*Ag^2*tan(Tf)^2+Ag^2*tan(Tf)^4+4*Af*tan(Tv)^2*Ag*tan(Tf)^2-4*Af^2*tan(Tv)^2+4*Af*tan(Tv)^2*Ag+8*Af^2*tan(Tv)*tan(Tf)-4*Af*tan(Tf)^3*Ag*tan(Tv)-4*Af*tan(Tf)*Ag*tan(Tv)-4*Af^2*tan(Tf)^2)^(1/2))*tan(Tv)^2)*Af/(-1/2/(-Af*tan(Tv)+Af*tan(Tf))*(Ag+Ag*tan(Tf)^2+(Ag^2+2*Ag^2*tan(Tf)^2+Ag^2*tan(Tf)^4+4*Af*tan(Tv)^2*Ag*tan(Tf)^2-4*Af^2*tan(Tv)^2+4*Af*tan(Tv)^2*Ag+8*Af^2*tan(Tv)*tan(Tf)-4*Af*tan(Tf)^3*Ag*tan(Tv)-4*Af*tan(Tf)*Ag*tan(Tv)-4*Af^2*tan(Tf)^2)^(1/2))-1/2/(-Af*tan(Tv)+Af*tan(Tf))*(Ag+Ag*tan(Tf)^2+(Ag^2+2*Ag^2*tan(Tf)^2+Ag^2*tan(Tf)^4+4*Af*tan(Tv)^2*Ag*tan(Tf)^2-4*Af^2*tan(Tv)^2+4*Af*tan(Tv)^2*Ag+8*Af^2*tan(Tv)*tan(Tf)-4*Af*tan(Tf)^3*Ag*tan(Tv)-4*Af*tan(Tf)*Ag*tan(Tv)-4*Af^2*tan(Tf)^2)^(1/2))*tan(Tf)^2+tan(Tv)*tan(Tf)^2+tan(Tv)) • (tan(Tf)+tan(Tf)*tan(Tv)^2-1/2/(-Af*tan(Tv)+Af*tan(Tf))*(Ag+Ag*tan(Tf)^2-(Ag^2+2*Ag^2*tan(Tf)^2+Ag^2*tan(Tf)^4+4*Af*tan(Tv)^2*Ag*tan(Tf)^2-4*Af^2*tan(Tv)^2+4*Af*tan(Tv)^2*Ag+8*Af^2*tan(Tv)*tan(Tf)-4*Af*tan(Tf)^3*Ag*tan(Tv)-4*Af*tan(Tf)*Ag*tan(Tv)-4*Af^2*tan(Tf)^2)^(1/2))-1/2/(-Af*tan(Tv)+Af*tan(Tf))*(Ag+Ag*tan(Tf)^2-(Ag^2+2*Ag^2*tan(Tf)^2+Ag^2*tan(Tf)^4+4*Af*tan(Tv)^2*Ag*tan(Tf)^2-4*Af^2*tan(Tv)^2+4*Af*tan(Tv)^2*Ag+8*Af^2*tan(Tv)*tan(Tf)-4*Af*tan(Tf)^3*Ag*tan(Tv)-4*Af*tan(Tf)*Ag*tan(Tv)-4*Af^2*tan(Tf)^2)^(1/2))*tan(Tv)^2)*Af/(-1/2/(-Af*tan(Tv)+Af*tan(Tf))*(Ag+Ag*tan(Tf)^2-(Ag^2+2*Ag^2*tan(Tf)^2+Ag^2*tan(Tf)^4+4*Af*tan(Tv)^2*Ag*tan(Tf)^2-4*Af^2*tan(Tv)^2+4*Af*tan(Tv)^2*Ag+8*Af^2*tan(Tv)*tan(Tf)-4*Af*tan(Tf)^3*Ag*tan(Tv)-4*Af*tan(Tf)*Ag*tan(Tv)-4*Af^2*tan(Tf)^2)^(1/2))-1/2/(-Af*tan(Tv)+Af*tan(Tf))*(Ag+Ag*tan(Tf)^2-(Ag^2+2*Ag^2*tan(Tf)^2+Ag^2*tan(Tf)^4+4*Af*tan(Tv)^2*Ag*tan(Tf)^2-4*Af^2*tan(Tv)^2+4*Af*tan(Tv)^2*Ag+8*Af^2*tan(Tv)*tan(Tf)-4*Af*tan(Tf)^3*Ag*tan(Tv)-4*Af*tan(Tf)*Ag*tan(Tv)-4*Af^2*tan(Tf)^2)^(1/2))*tan(Tf)^2+tan(Tv)*tan(Tf)^2+tan(Tv))

  36. Av = atan(1/2/(-Af*tan(Tv)+Af*tan(Tf))*(Ag+Ag*tan(Tf)^2+(Ag^2+2*Ag^2*tan(Tf)^2+Ag^2*tan(Tf)^4+4*Af*tan(Tv)^2*Ag*tan(Tf)^2-4*Af^2*tan(Tv)^2+4*Af*tan(Tv)^2*Ag+8*Af^2*tan(Tv)*tan(Tf)-4*Af*tan(Tf)^3*Ag*tan(Tv)-4*Af*tan(Tf)*Ag*tan(Tv)-4*Af^2*tan(Tf)^2)^(1/2))) atan(1/2/(-Af*tan(Tv)+Af*tan(Tf))*(Ag+Ag*tan(Tf)^2-(Ag^2+2*Ag^2*tan(Tf)^2+Ag^2*tan(Tf)^4+4*Af*tan(Tv)^2*Ag*tan(Tf)^2-4*Af^2*tan(Tv)^2+4*Af*tan(Tv)^2*Ag+8*Af^2*tan(Tv)*tan(Tf)-4*Af*tan(Tf)^3*Ag*tan(Tv)-4*Af*tan(Tf)*Ag*tan(Tv)-4*Af^2*tan(Tf)^2)^(1/2)))

  37. You waited all this time! I thought I get enough time to play with you now !!!!!!!!!!!!!!!!! ? ChandanaG

  38. Phase topography vs. canopy surface VEG HT PHASE OF THE BARE EARTH IS NOT CORRELATED WITH THAT OF IN THE SRTM SIGNAL. SOLVE ITERATIVELY BUT THE AMPLITUDE OF THE FOURIER TRANSFORM OF BARE EARTH CAN BE FOUND AT LARGE SCALE AND CAN BE SYNTHESIZED AT SMALL SCALES.

  39. Results.. Results are obtained before solving the equation iteratively. Iterative solutions results a narrow Pro. Dis Function.

More Related