radar-science / Yang-2024-TropoTexCorr

Code and Notebooks for Yang et al. (2024) on InSAR tropospheric correction based on local texture correlation
https://doi.org/10.1109/TGRS.2024.3356749
Apache License 2.0
10 stars 4 forks source link

Question about running code #1

Open mohseniaref opened 8 months ago

mohseniaref commented 8 months ago

Hi,

Thanks so much for sharing the code. I am trying to run Texturecorrelation.m using my data. I would really appreciate if you could help me in this regard I am not sure about parameter but this is how i set it

N = 86;(number of date)
N_ref = 2;
Na_ref = 5113;azimuth number of pixel
Nr_ref = 5503;range number of pixel
lamda = 0.56;

and after this line

[k_LLF,d_LLF,k_HTC] = slope_estimation(phase_ts_deramp,dem,MASK,W,r,N_ref); % **6

I get this error. How can i fix this?

Warning: Polynomial is not unique; degree >= number of data points. 
> In polyfit (line 79)
In slope_estimation (line 57) 
Warning: Polynomial is not unique; degree >= number of data points. 
> In polyfit (line 79)
In slope_estimation (line 57) 
Warning: Polynomial is not unique; degree >= number of data points. 
> In polyfit (line 79)
In slope_estimation (line 57) 
Warning: Polynomial is not unique; degree >= number of data points. 
> In polyfit (line 79)
In slope_estimation (line 57) 
Warning: Polynomial is not unique; degree >= number of data points. 
> In polyfit (line 79)
In slope_estimation (line 57) 
Warning: Polynomial is not unique; degree >= number of data points. 
> In polyfit (line 79)
In slope_estimation (line 57) 
Warning: Polynomial is not unique; degree >= number of data points. 
> In polyfit (line 79)
In slope_estimation (line 57) 
Warning: Polynomial is not unique; degree >= number of data points. 
> In polyfit (line 79)
In slope_estimation (line 57) 
Warning: Polynomial is not unique; degree >= number of data points. 
> In polyfit (line 79)
In slope_estimation (line 57) 
Warning: Polynomial is not unique; degree >= number of data points. 
> In polyfit (line 79)
In slope_estimation (line 57) 
Unable to perform assignment because the size of the left side is 1-by-1 and the size of the
right side is 1-by-41.

Error in slope_estimation (line 122)
            Results_Compare(3,n) = k(index);
yunjunz commented 8 months ago

Hi @mohseniaref, thank you for your interest in trying our code! @yangqy916 has been busy writing her dissertation recently, so I would expect a slow response from her regarding your issue above.

I will try her Matlab code with the dataset used in the manuscript and let you know how it goes if I find some time.