How did you get that geocoding? Or don’t you use Sentinel-1 data? All I can get is stretched image, without any position information.
Also, even after update, topo-phase removal doesn’t work. Problem is the same. During step 3 unnaturally huge amount of scatterers is found, which is completely filtered out during step 4.
It’s hard to explain how I did but I will try my best!
Yes, I have the same problem. I continue to have to modify the name of the bands i and q as @lveci has sugested above!
Do you rename them in file system, or using SNAP?
I have 20 slave images.
Here is my parameters. Right now almost everything except ‘unwrap_gold_n_win’ is left as default.
Created: '27-Dec-2016'
clap_alpha: 1
clap_beta: 0.3000
clap_low_pass_wavelength: 800
clap_win: 32
density_rand: 20
drop_ifg_index: []
filter_grid_size: 50
filter_weighting: 'P-square'
gamma_change_convergence: 0.0050
gamma_stdev_reject: 0
heading: 250.8561
insar_processor: 'gamma'
lambda: 0.0555
lonlat_offset: [0 0]
max_topo_err: 5
merge_resample_size: 0
merge_standard_dev: Inf
percent_rand: 20
plot_color_scheme: 'inflation'
plot_dem_posting: 90
plot_pixels_scatterer: 3
plot_scatterer_size: 120
quick_est_gamma_flag: 'y'
ref_centre_lonlat: [0 0]
ref_lat: [-Inf Inf]
ref_lon: [-Inf Inf]
ref_radius: Inf
scla_deramp: 'y'
scla_drop_index: []
scla_method: 'L2'
scn_deramp_ifg: []
scn_time_win: 365
scn_wavelength: 100
select_method: 'DENSITY'
shade_rel_angle: [90 45]
slc_osf: 1
small_baseline_flag: 'n'
subtr_tropo: 'n'
tropo_method: 'a_l'
unwrap_alpha: 8
unwrap_gold_alpha: 0.8000
unwrap_gold_n_win: 16
unwrap_grid_size: 200
unwrap_la_error_flag: 'y'
unwrap_method: '3D'
unwrap_patch_phase: 'n'
unwrap_prefilter_flag: 'y'
unwrap_spatial_cost_func_flag: 'n'
unwrap_time_win: 730
weed_max_noise: Inf
weed_neighbours: 'y'
weed_standard_dev: 1
weed_time_win: 730
weed_zero_elevation: 'n'Using SNAP. It’s the safest way, I think!
Like me. I also had to change that parameter.
I never studied that step. Do you have a reason to choose 16 or it was by chance?
Then SNAP already renames all i and q bands itself
I had this error message during Step 6 :
Minimum dimension of the resampled grid (17 pixels) is less than prefilter window size (32). So according to one mainsar forum topic I have changed the value to 16.
Dear Jaroslav and juanesburgo,
I managed to fulfill export to StaMPS when I re-run again Interferogram formation and TopoPhaseRemoval for the whole co-registered stack. It is not long.
But I still have the problem with the exported data. Below are the values which are imported by StaMPS. Heading is definitely wrong. Moreover look angles are not imported… Could you, please, send me your screen of Step 1 of StaMPS? Or probably you can give me any advise.
Thank you for your time
Katherine
jaroslav
December 28
juanesburgo:Using SNAP. It’s the safest way, I think!
Then SNAP already renames all i and q bands itself
juanesburgo:
Do you have a reason to choose 16 or it was by chance?
I had this error message during Step 6 :
Minimum dimension of the resampled grid (17 pixels) is less than prefilter window size (32). So according to one mainsar forum topic I have changed the value to 16.
Visit Topic or reply to this email to respond.
In Reply To
juanesburgo
December 28
Using SNAP. It’s the safest way, I think!
Like me. I also had to change that parameter. I never studied that step. Do you have a reason to choose 16 or it was by chance?
Visit Topic or reply to this email to respond.
To unsubscribe from these emails, click here .
All the best,
Katherine
How did you import hdr files into matlab? hdrread function didn’t work for me.
Did you change ijname=['pscands.1.ij']; line to ijname=['lonat.1.ij'];?
I didn’t get this kind of error, but for some reason your pscands.1.ij file is empty. I didn’t find, why can it come out like this though.
Dear Jaroslav,
Thank you very much for your comments, I have not noticed that this time pscands.1.ij are empty. Before updating it was OK. And the same question about heading - did you manage to get normal heading value for S-1 in Step 1 StaMPS?
Thank you for your time
Katherine
Hi all,
Maybe this helps, small Andy Hooper tutorial for Sentinel-1, he use different software for interferogram generation but also pointed out if we can export to StaMPS, this should work
http://seom.esa.int/landtraining2015/files/Day_4/D4P2a_LTC2015_Hooper.pdf
I hope it will help to correctly set the parameters
Best regards, Stanislav
You can read those files as ENVI files. There is matlab code to do that out there.
No, I didn’t.
I did this:
if exist(llname,'file')
fid=fopen(llname,'r');
lonlat=fread(fid,[2,inf],'float',endian);
lonlat=lonlat';
fclose(fid);
else
load('~/lonlat.mat');
%lonlat=local2llh(xy'/1000,[0;0])';
end
my modified code is between else and end.
Thank you very much Stanislav !
All the best
Katherine
Hi @ABraun ,
I started using the workflow between SNAP and stamps and I am currently facing the same problem as you did last month, namely that " Error opening file pscands.1.ij" message in the end of the mt_prep_gamma step. Did you by chance solve this problem?
Cheers,
Maryse
unfortunately not. Do you work on windows or linux?
I work on a Linux machine.
I think, that maybe the problem is due to this spelling error:
file name for zero amplitude PS: pscands.1.ij0
instead of pscands.1.ij
I think that the problem lies within the selpsc_patch script. I will try to get more information in the mainsar forum.
Maryse
alright, at least it is not platform-related.
You could be right about the variable pscands.1.ij0
But strange that it only happened to us… Please let me know if you receive an answer in the mainsar forum.
@ABraun, the solution is quite simple. Try to run mt_prep_gamma inside your “export” folder, like this:
mt_prep_gamma 20150806 /data2/…/Export_to_stamps/ 0.4 2 2 50 200
where:
mt_prep_gamma [your master] [full path to your working directory] [Amplitude disp.] [Number of patches in range] [Number of patches in azimuth] [overlapping pixels between patches in range] [overlapping pixels between patches in azimuth]
It worked for me, at least I didn’t get the error message concerning opening the pscands.1.ij file.
Unfortunately, I now got the same issue as @katherine before, namely that the pscands.1.ij file is empty…
you are right, thanks for clarification. It works for me now, too.
At least the error about opening pscands.1.ij is gone. But I can’t tell if the output makes sense, yet.
After SNAP export to STAMP, STAMP protsessing are possible Windows matlab?
you first need to run mt_prep_gamma
http://homepages.see.leeds.ac.uk/~earahoo/stamps/StaMPS_Manual_v3.3b1.pdf (beginning in chapter )
Please see also here:
The rest can be done in matlab.