add_lidar_elevation.py

[falkamelung]

That would be another useful script to have. It reads the lat_geocorr,lon_geocorr from the csv file and adds another column lidar_dem. This will then be displayed in insarmaps as well. To be more precise, we should set a radius (or rectangle half-length) of 5 meters and calculate the mean and variance of the LIDAR DEM within this area. If the variance is high we may be looking at a scatterer on the facade. Scatterers on the rooftop and ground will have low variance.

It should check whether the x_geocorr column exist and exit if it does not.

There could be an option --radius to select the radius for the mean and variance calculation (Default: 5 meters).

add_lidar_elevation.py  MiamiSenA28_MDCBeaches_20191001-20231031.csv ----dem   DEM/MiamiBeach.tif
add_lidar_elevation.py  MiamiSenA28_MDCBeaches_20191001-20231031.csv ----dem   DEM/MiamiBeach.tif --radius 10

This will then also be displayed in insarmaps.

• We actually have to carefully look at the geolocation of our pixels. I am not sure that we take the geoid calculation properly into account. That was something on my list to check. I have previously seen that our pixel locations may be off by 5-20 meters or so: Geolocation problem examples. Cape Canaveral and Banda Api #557