Missing cum_filt.h5 file and only cum.h5 exits. about licsbas HOT 10 OPEN

I used the batch processing script for processing. after it was finished I could only find cum.h5 file inside the folder. what process did i do wrong?

Please show me the full log.

when doing batch processing , does it take ascending or descending interferograms?

It depends on the frame you choose. Your result looks ascending.

Also how to find the LOS and flight direction ?

GEOCml*/[ENU].geo contain LOS information as unit vectors.

How to conver the velocity to EW, NS and vertical displacement. Also i cannot move the black dot which is already in map. whenever i click anywhere i cannot generate the time series on second image. Only when i draw a region, it gives me time series. what could be wrong?

We cannot convert a single LOS to 3D.
Please double-click to select a point as written in the left top of the window.

from licsbas.

Thank you for your kind reply. I wanted to compare GNSS time series and InSAR time series. I obtained the GNSS time series from NGSS website with 'NS', 'EW' and 'UD' displacement and I found out that I need to convert the GNSS displacement to LOS using

LOS = D(NS)sin(θ) cos(α) + D(EW)sin(θ) sin(α) + D(UD)cos(θ)

where
D(NS) = displacement in NS of GNSS
D(EW) = displacement in EW of GNSS
D(UD) = displacement in UD of GNSS

and θ is the incidence angle and α is the azimuthal angle. where can i find α and θ ? My aim is to find the subsidence rate at my study area....what is best way to find subsidence rate from LOS velocity ?

from licsbas.

GEOCml*/E.geo, N.geo, and U.geo are the unit vectors. You can compute LOS by D(NS)*N.geo+D(EW)*E.geo+D(UD)*U.geo.

from licsbas.

It should be very basics and i might sound very stupid. But how can I read N.geo, E.geo and U.geo file? the GPS data file is in csv format with four column 'Date' , 'NS', 'EW' and 'UD'. I tried to open .geo format with QGIS both vector and raster but it is is not opening giving me error as not supported raster data source. I also used notepad to open any of .geo file but it is showing weird characters.

from licsbas.

You need to convert the float file (*.geo) to GeoTIFF using LiCSBAS_flt2geotiff.py.

from licsbas.

I am currently using sample to learn more about LiCSBAS before i actually do processing and the sample doesn't appear to contain N.geo, E.geo and U.geo file in it LiCSBAS_sample_CF.tar.gz , so it would be helpful if these could be provided with sample.

and next is that i earlier i used batch script to process one of my study area and when i shared the output with others, i receive a comment that 'obital ramps' might not have been removed. Is it true ? isn't the batch script should remove it? if not how can i remove it? following is my output

from licsbas.

You can create the E/N/U.geo for the sample by following this instruction: https://github.com/yumorishita/LiCSBAS/wiki/2_0_workflow#2-running-licsbas--quick-start

By default, LiCSBAS does not apply deramping because the orbital and ionospheric ramp is generally very small in Sentinel-1. You can activate deramping here:

from licsbas.

You can create the E/N/U.geo for the sample by following this instruction: https://github.com/yumorishita/LiCSBAS/wiki/2_0_workflow#2-running-licsbas--quick-start

By default, LiCSBAS does not apply deramping because the orbital and ionospheric ramp is generally very small in Sentinel-1. You can activate deramping here:

LiCSBAS/batch_LiCSBAS.sh

Line 62 in 5db6516

p16_deg_deramp="" # 1, bl, or 2. default: no deramp

It appears that to create E/N/U.geo file, instead of downloading the results provided into zip format , i will need to run the batchscript from and start from first step. It will take some time...I'll comment if i get any errors. Thank you for all your help.

from licsbas.

LiCSBAS_flt2geotiff.py.

Thank you for your comments Prof. Yu Morishita, I applied LiCSBAS_flt2geotiff.py to convert N.geo file to raster file and used QGIS to open the file. If I have a GPS station in study area, do I have to use the pixel value of raster file using same coordinates of GNSS station and multiply that with D(NS) in following formula?

D(NS)*N.geo+D(EW)*E.geo+D(UD)*U.geo.

Here is illustration, Using LiCSBAS_flt2geotiff.py. to convert N.geo >N.geo.tif and importing it into QGIS

Then using pixel value of N.geo.tif to multiply with North displacement of GNSS station

and repeat same with U.geo and E.geo and finally using summation 'D(NS)*N.geo+D(EW)*E.geo+D(UD)*U.geo.' to get the LOS of 3D GNSS displacement as per this formula 'D(NS)*N.geo+D(EW)*E.geo+D(UD)*U.geo.'. Am I in the right direction Professor ?

from licsbas.

Yes, I think you are right.

from licsbas.