pi_pick.m

基于matlab的反演程序,用于地球物理勘探中射线追踪及偏移成像程序.

function pick(arg1,arg2,arg3);
global PICKS PICKCOLOR XAXISTOP ZOOM_VALUE ZOOM_LOCKS
hax=findobj(gcf,'type','axes','tag','MAINAXES');
clickcheck=get(gcf,'selectiontype');
h=get(gcf,'userdata');
hmsg=h(2);
hi=h(5);
xdat=get(hi,'xdata');
ydat=get(hi,'ydata');
hzoompick=h(9);
value=get(hzoompick,'value');
if(strcmp(clickcheck,'normal'))
    h=get(gcf,'userdata');
    hzoompick=h(9);
    if(value~=3)
        pick=drawlinefini;
        
    else
        % user has called auto picker, need to get picks info from handle
        % where picks have been stored
        picksout=picksle('EXPORT');
        % need to set data so this thingy can plot properly... thingy is
        % the technical term of course
        % must be [x1 y1 x2 y2] style
        if(size(picksout,1)<1)
            stringinfo='Automatic picking did not find any picks';
            set(hmsg,'string',stringinfo,'backgroundcolor',[1 1 1]);
            % not enough points to make a line
            return
        end
        pick=[];
        for ii=1:size(picksout,1)-1
            pick01=[picksout(ii,2) picksout(ii,1) picksout(ii+1,2) picksout(ii+1,1)];
            pick=[pick;pick01];
        end
    end
    if(size(pick,2)~=5)
        
    elseif(length(pick)==5)
        if(iscell(pick))
            return
        elseif(pick(5)>0)
            delete(pick(5));
        end
        pick=[pick(1) pick(2) pick(3) pick(4)];
        test=pick(1)-pick(3)+pick(4)-pick(2);%zero for a double click
        ck1=sort([pick(1)-xdat(1) xdat(end)-pick(1) pick(3)-xdat(1) xdat(end)-pick(3),...
                pick(2)-ydat(1) ydat(end)-pick(2) pick(4)-ydat(1) ydat(end)-pick(4)]);
        if(test==0)
            stringinfo='Please create Pick with greater then zero length';
            col=[1 1 0];
        elseif(ck1(1)<=0)
            stringinfo='Please keep Pick within the image area';
            col=[1 1 0];
        end
    end
    for ii=1:size(pick,1)
        % first two are markes that will appear at the ends
        % of each pick line when moving is instigated
        hpick=line([pick(ii,1) pick(ii,3)],[pick(ii,2) pick(ii,4)],[1 1],'linewidth',2,...
            'color',PICKCOLOR,'buttondownfcn','plotimage(''picklinemenu'')',...
            'userdata','','tag','PICKS');
        if(isempty(PICKS))
            PICKS=cell(1,3);
            PICKS{1}=gcf;
            PICKS{2}=[pick(ii,1) pick(ii,2) pick(ii,3) pick(ii,4)];
            PICKS{3}=hpick;
        else
            checkpicks=[];
            for jj=1:size(PICKS,1)
                CheckMasterFigure=PICKS{jj,1};
                if(CheckMasterFigure==gcf)
                    PicksPositions=PICKS{jj,2};
                    PicksHandles=PICKS{jj,3};
                    PicksPositions=[PicksPositions;pick(ii,1) pick(ii,2) pick(ii,3) pick(ii,4)];
                    PicksHandles=[PicksHandles;hpick];
                    PICKS{jj,2}=PicksPositions;
                    PICKS{jj,3}=PicksHandles;
                    checkpicks='Pick Has Been Set';
                end
            end
            if(isempty(checkpicks))
                % for some reason, need to add picks.
                NewPICK={gcf [pick(ii,1) pick(ii,2) pick(ii,3) pick(ii,4)] hpick};
                PICKS(size(PICKS,1)+1,:)=NewPICK;
            end
        end
    end
    stringinfo=['Pick has been added'];
    col=[1 1 1];
    set(hmsg,'string',stringinfo,'backgroundcolor',col);
end