plotclr.m

GPS多路径效应的谱分析工具

function plotclr(x,y,v,prn,marker)
% [e-mail]

if nargin <5
    marker='.';
end

map=colormap;
miv=min(min(v));
mav=max(max(v));

[row col]=size(x);

% Plot the points
hold on
for j=1:col,
    for i=1:row%length(x)
        in=round((v(i,j)-miv)*(length(map)-1)/(mav-miv));
        %--- Catch the out-of-range numbers
        if in==0;in=1;end
        if in > length(map);in=length(map);end
        if isnan(in);in=1;end
        plot3(x(i,j),y(i,j),v(i,j),marker,'color',map(in,:),'markerfacecolor',map(in,:))
        if i==row,
            if prn(j) < 10,
                ht=text(x(i,j),y(i,j),[' S0',num2str(prn(j))]);set(ht,'fontsize',8);
            else
                ht=text(x(i,j),y(i,j),[' S',num2str(prn(j))]);set(ht,'fontsize',8);
            end
        end
    end
end
hold off

% Re-format the colorbar
h=colorbar;

set(h,'fontsize',8);
set(get(h,'ylabel'),'string','[m]');
set(h,'ylim',[1 length(map)]);
yal=linspace(1,length(map),6);
set(h,'ytick',yal);
% Create the yticklabels
ytl=linspace(miv,mav,6);
s=char(6,4);
for i=1:6
    if min(abs(ytl)) >= 0.001
        B=sprintf('%4.2f',ytl(i));
    else
        B=sprintf('%4.2E',ytl(i));
    end
    s(i,1:length(B))=B;
end
set(h,'yticklabel',s);
grid on
view(2)