plot_sb_baselines.m

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

function []=plot_sb_baselines()
%PLOT_SB_BASELINES plot the small baselines in small_baselines.list
%
%   Andy Hooper, June 2007



currdir=pwd;
dirs=strread(currdir,'%s','delimiter','/');
if strcmp(dirs{end},'SMALL_BASELINES') 
    !\ls -d [1,2]* | sed 's/_/ /' > small_baselines.list
    load ../psver
    psname=['../ps',num2str(psver)];
else
    load psver
    psname=['ps',num2str(psver)];
end

sb=load('small_baselines.list');
ps=load(psname);

n_ifg=size(sb,1);
[yyyymmdd,I,J]=unique(sb);
ifg_ix=reshape(J,n_ifg,2);
x=ifg_ix(:,1);
y=ifg_ix(:,2);


day=str2num(datestr(ps.day,'yyyymmdd'));
[B,I]=intersect(day,yyyymmdd);

x=I(x);
y=I(y);

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}')