sb_find.m

来自「StaMps最新测试版」· M 代码 · 共 79 行

function []=sb_find(rho_min,ddiff_max,bdiff_max)
%SB_FIND find small baselines
%   SB_FIND(RHO_MIN,DDIFF_MAX,BDIFF_MAX)
%
%   Andy Hooper, September 2006
%
%   ======================================================================
%   04/2008 AH: Criteria for selection changed to coherence basis 
%   ======================================================================


if nargin<1
    rho_min=0.5
end

if nargin<2
    ddiff_max=1500
end

if nargin<3
    bdiff_max=1070
end

load psver
psname=['ps',num2str(psver)];
ps=load(psname);


[X,Y]=meshgrid(ps.bperp',ps.bperp);
bdiff=abs(X-Y);
[X,Y]=meshgrid(ps.day',ps.day);
ddiff=abs(X-Y);
N=size(X,1);

rho=(1-bdiff/bdiff_max).*(1-ddiff/ddiff_max); % correlation
rho(bdiff>bdiff_max|ddiff>ddiff_max)=0;
rho=rho-eye(N);
ix=(rho>rho_min&~tril(ones(ps.n_ifg)));

[dummy,I]=max(rho);
best_ix=(sort([1:N;I]));
for i=1:N
    ix(best_ix(1,i),best_ix(2,i))=1;
end
[x,y]=find(ix);


%[x,y]=find((rho>rho_min|I|I2)&~tril(ones(ps.n_ifg)));
%[x,y]=find(bdiff<bdiff_max&ddiff<ddiff_max&~tril(ones(ps.n_ifg)));

no_conn_ix=find(sum(rho>rho_min)==0);



sortxy=sortrows([x,y]);
x=sortxy(:,1);
y=sortxy(:,2);

clf

for i=1:length(x)
    l=line([ps.day(x(i)),ps.day(y(i))],[ps.bperp(x(i)),ps.bperp(y(i))]);
    text((ps.day(x(i))+ps.day(y(i)))/2,(ps.bperp(x(i))+ps.bperp(y(i)))/2,num2str(i));
    set(l,'color',[0 1 0],'linewidth',2)
end
hold on
p=plot(ps.day,ps.bperp,'ro');
set(p,'markersize',12,'linewidth',2)
hold off
datetick('x',12)
ylabel('B_{perp}')

fid=fopen('small_baselines.list','w');
for i=1:length(x)
    fprintf(fid,'%s %s\n',datestr(ps.day(x(i)),'yyyymmdd'),datestr(ps.day(y(i)),'yyyymmdd'));
end
fclose(fid);