Index: /issm/trunk-jpl/src/m/utils/Exp/clicktoflowline.m =================================================================== --- /issm/trunk-jpl/src/m/utils/Exp/clicktoflowline.m (revision 11806) +++ /issm/trunk-jpl/src/m/utils/Exp/clicktoflowline.m (revision 11807) @@ -1,16 +1,24 @@ function clicktoflowline(index,x,y,u,v,filename) -%CLICKTOFLOWLINE create a flowline contour file (named 'filename') by clicking on a velocity field once. +%CLICKTOFLOWLINE - create a flowline ARGUS file % -% Usage: clicktoflowline(index,x,y,u,v,x0,y0,filename) +% create a flowline contour file (named 'filename') by clicking +% on a velocity field once (velocity must be plotted first) % -% Ex: clicktoflowline(md.mesh.elements,md.mesh.x,md.mesh.y,md.inversion.vx_obs,md.inversion.vy_obs,'flowline1.exp') +% Usage: +% clicktoflowline(index,x,y,u,v,x0,y0,filename) +% +% Example: +% clicktoflowline(md.mesh.elements,md.mesh.x,md.mesh.y,md.inversion.vx_obs,md.inversion.vy_obs,'flowline.exp') - +%Get click position [x0,y0]=ginput(1); +%Get flowline line=flowlines(index,x,y,u,v,x0,y0,200); +%plot hold on plot(line.x,line.y,'r-'); +%Write argus file expwrite(line,filename); Index: sm/trunk-jpl/src/m/utils/Exp/zinput.m =================================================================== --- /issm/trunk-jpl/src/m/utils/Exp/zinput.m (revision 11806) +++ (revision ) @@ -1,211 +1,0 @@ -function [out_regs] = zinput(arg1) -% ZINPUT - Graphical input from mouse with zoom -% -% [OUT_REGS] = ZINPUT(N) gets N points or regions from the -% current axes and returns the X- and Y-ranges in a length Nx4 -% matrix OUT_REGS. -% The cursor can be positioned using a mouse. Data points are entered by -% pressing the right mouse button, the region of the current axis -% are selected with the middle button, left button is for -% zooming, single cklick zooms in, click-and drag zooms to -% region (and doubble-click should zoom out - feature pending). -% Any key on the keyboard except carriage return zooms out to the -% orignal axis except carriage return, which terminates the input -% before N points are entered. -% -% [OUT_REGS] = ZINPUT gathers an unlimited number of points until the -% return key is pressed. -% -% Usage: -% [out_regs] = zinput(arg1) -% -% See also GINPUT - -% Copyright Bjorn Gustavsson 20050314 - -ax0 = axis; -out_regs = []; -c = computer; -if ~strcmp(c(1:2),'PC') - tp = get(0,'TerminalProtocol'); -else - tp = 'micro'; -end - -if ~strcmp(tp,'none') & ~strcmp(tp,'x') & ~strcmp(tp,'micro') & 0, - % I dont know about this so better make short-cut and blindly try - % what works for X in all environments - sorry about that. -else - - fig = gcf; - figure(gcf); - - if nargin == 0 - how_many = inf; - b = []; - else - how_many = arg1; - b = []; - if isstr(how_many) ... - | size(how_many,1) ~= 1 | size(how_many,2) ~= 1 ... - | ~(fix(how_many) == how_many) ... - | how_many < 0 - error('Requires a positive integer.'); - end - if how_many == 0 - ptr_fig = 0; - while(ptr_fig ~= fig) - ptr_fig = get(0,'PointerWindow'); - end - scrn_pt = get(0,'PointerLocation'); - loc = get(fig,'Position'); - pt = [scrn_pt(1) - loc(1), scrn_pt(2) - loc(2)]; - out1 = pt(1); y = pt(2); - elseif how_many < 0 - error('Argument must be a positive integer.'); - end - end - - % Remove figure button functions - state = uisuspend(fig); - pointer = get(gcf,'pointer'); - set(gcf,'pointer','fullcrosshair'); - - fig_units = get(fig,'units'); - char = 0; - while size(out_regs,1) < how_many - % Use no-side effect WAITFORBUTTONPRESS - waserr = 0; - try - keydown = wfbp; - catch - waserr = 1; - end - if(waserr == 1) - if(ishandle(fig)) - set(fig,'units',fig_units); - uirestore(state); - error('Interrupted'); - else - error('Interrupted by figure deletion'); - end - end - - ptr_fig = get(0,'CurrentFigure'); - if(ptr_fig == fig) - if keydown - axis(ax0); - char = get(fig, 'CurrentCharacter'); - else - char = get(fig, 'CurrentCharacter'); - button = abs(get(fig, 'CurrentCharacter')); - - pnt = get(gcf,'currentpoint'); - xy1 = get(gca,'currentpoint'); - rbbox([pnt 0 0],pnt) - xy2 = get(gca,'currentpoint'); - selection_type = get(gcf,'selectiontype'); - - if all(xy1==xy2) - ax1 = axis; - - dx = abs(ax1(2)-ax1(1)); - dy = abs(ax1(4)-ax1(3)); - - xmin = max(ax1(1),xy1(1)-dx/4); - xmax = min(ax1(2),xy1(1)+dx/4); - ymin = max(ax1(3),xy1(1,2)-dy/4); - ymax = min(ax1(4),xy1(1,2)+dy/4); - zoom2ax = [xmin xmax ymin ymax]; - else - %%% zoom to selected rectangle - zoom2ax = [sort([xy1(1,1) xy2(1,1)]) sort([xy1(1,2) xy2(1,2)])]; - end - - switch selection_type - case 'normal' - axis(zoom2ax) - reg = []; - case 'alt' - reg = axis; - axis(ax0) - case 'extend' - reg = xy1(1,[1 1 2 2]); - axis(ax0) - otherwise %%% open (doubleclick in linux) - axis(ax0); - reg = []; - end - end - pt = get(gca, 'CurrentPoint'); - - if (char == 'r') - rmi = size(out_regs,1); - if rmi > 0 - out_regs(rmi,:) = []; - end - set(fig, 'CurrentCharacter','q') - char = 'q'; - reg = []; - end - if(char == 13) % & how_many ~= 0) - % if the return key was pressed, char will == 13, - % and that's our signal to break out of here whether - % or not we have collected all the requested data - % points. - % If this was an early breakout, don't include - % the key info in the return arrays. - % We will no longer count it if it's the last input. - break; - end - - if ~isempty(reg) - out_regs = [out_regs;reg]; - end - end - %[size(out_regs,1), how_many, size(out_regs,1) < how_many] - end - - uirestore(state); - set(gcf,'pointer','arrow'); - set(fig,'units',fig_units); - set(fig, 'CurrentCharacter','q'); - -end - -%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -function key = wfbp -%WFBP Replacement for WAITFORBUTTONPRESS that has no side effects. - -fig = gcf; -current_char = []; - -% Now wait for that buttonpress, and check for error conditions -waserr = 0; -try - h=findall(fig,'type','uimenu','accel','C'); % Disabling ^C for edit menu so the only ^C is for - set(h,'accel',''); % interrupting the function. - keydown = waitforbuttonpress; - current_char = double(get(fig,'CurrentCharacter')); % Capturing the character. - if~isempty(current_char) & (keydown == 1) % If the character was generated by the - if(current_char == 3) % current keypress AND is ^C, set 'waserr'to 1 - waserr = 1; % so that it errors out. - end - end - - set(h,'accel','C'); % Set back the accelerator for edit menu. -catch - waserr = 1; -end -drawnow; -if(waserr == 1) - set(h,'accel','C'); % Set back the accelerator if it errored out. - error('Interrupted'); -end - -selection_type = get(gcf,'selectiontype'); -if strcmp(selection_type,'open') - axis(ax0) -end -if nargout>0, key = keydown; end -%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%