https://wiki-pnb.eri.ucsb.edu/index.php?action=history&feed=atom&title=Bt_ac9.m 2026-05-22T09:33:58Z Revision history for this page on the wiki MediaWiki 1.34.0 https://wiki-pnb.eri.ucsb.edu/index.php?title=Bt_ac9.m&diff=181&oldid=prev 2010-08-02T23:55:01Z

<p>Created page with &#039;&lt;pre&gt; % command file for matlab bottom top determination: Daniel C. Konnoff % 10/08/92 % modified: JBU Oct 97 - fits new 1997 ac9 processing scheme % - yadda yadda yadda …&#039;</p> <p><b>New page</b></p><div>&lt;pre&gt;<br /> % command file for matlab bottom top determination: Daniel C. Konnoff<br /> % 10/08/92<br /> % modified: JBU Oct 97 - fits new 1997 ac9 processing scheme<br /> % - yadda yadda yadda<br /> % modified: SRW 98 - now a function that returns [d]<br /> function [d]=bt_ac9(mat_in);<br /> clf;<br /> % set graph geometry<br /> set (gcf,'units','inches', 'pos',[2 1 9 5])<br /> % declare variables<br /> u=[]; % holds the points<br /> d=[]; % holds the index's only<br /> cast=0;<br /> plot(mat_in(:,1),mat_in(:,2),'.');grid; <br /> while(1)<br /> cast=cast+1<br /> which=['Selecting points for cast: ',num2str(cast)];<br /> % first point<br /> disp('select the first region')<br /> title(which);<br /> [x,y]=ginput(2);<br /> % zoom to start point of the cast<br /> if y(1)&lt;y(2)<br /> if x(1)&lt;x(2)<br /> axis([x(1),x(2),y(1),y(2)])<br /> else<br /> axis([x(2),x(1),y(1),y(2)])<br /> end<br /> else<br /> if x(1)&lt;x(2)<br /> axis([x(1),x(2),y(2),y(1)])<br /> else<br /> axis([x(2),x(1),y(2),y(1)])<br /> end % end if<br /> end % end if<br /> % pretty the graph<br /> hold on<br /> plot(mat_in(:,1),mat_in(:,2),'ro'); grid; title(which);<br /> plot([x(1), x(2)], [0, 0], 'g-.'); <br /> hold off;<br /> % find specific start point<br /> disp('click on start of cast')<br /> [x,y]=ginput(1);<br /> % add the points to cast_id array<br /> u=[u;x,y]; <br /> % last point<br /> plot(mat_in(:,1),mat_in(:,2),'.');grid; title(which);<br /> hold on;<br /> %d=ceil(u(:,1))<br /> d=[];<br /> for id=1:size(u,1),<br /> dd=find(mat_in(:,1)&lt;u(id,1));%___________________________<br /> d=[d; dd(end)];%-----------------------------------------<br /> end<br /> for ip=1:2:size(d,1),<br /> plot(u(:,1),u(:,2),'ro');<br /> end<br /> % $$$ plot(mat_in(d(1):d(end),1),...<br /> % $$$ mat_in(d(1):d(end),2),'g',u(:,1),u(:,2),'ro');<br /> disp('select the end region')<br /> [x,y]=ginput(2);<br /> % zoom to end point of the cast<br /> if y(1)&lt;y(2)<br /> if x(1)&lt;x(2)<br /> axis([x(1),x(2),y(1),y(2)])<br /> else<br /> axis([x(2),x(1),y(1),y(2)])<br /> end<br /> else<br /> if x(1)&lt;x(2)<br /> axis([x(1),x(2),y(2),y(1)])<br /> else<br /> axis([x(2),x(1),y(2),y(1)])<br /> end % end if<br /> end % end if<br /> % pretty the graph<br /> hold on<br /> plot(mat_in(:,1),mat_in(:,2),'ro');grid; title(which);<br /> plot([x(1), x(2)], [0, 0], 'g-.'); <br /> hold off;<br /> % find specific end point<br /> disp('click on end of cast')<br /> [x,y]=ginput(1);<br /> % add the points to cast_id array<br /> u=[u;x,y]; <br /> % are we done yet?<br /> axis('auto');<br /> plot(mat_in(:,1),mat_in(:,2),'.');grid; title(which);<br /> hold on;<br /> %d=ceil(u(:,1))<br /> d=[];<br /> for id=1:size(u,1),<br /> dd=find(mat_in(:,1)&lt;u(id,1));<br /> d=[d; dd(end)];<br /> end<br /> for ip=1:2:size(d,1),<br /> plot(mat_in(d(ip):d(ip+1),1),...<br /> mat_in(d(ip):d(ip+1),2),'g',u(:,1),u(:,2),'ro');<br /> end <br /> % $$$ plot(mat_in(d(1):d(end),1),...<br /> % $$$ mat_in(d(1):d(end),2),'g',u(:,1),u(:,2),'ro');<br /> title('Selected points'); grid;<br /> disp(' ');<br /> %disp('Want to save these points?');<br /> %done=input('y or n: ','s');<br /> %disp(' ')<br /> % $$$ if (done=='y'|done=='Y')<br /> % $$$ d(1,1)=ceil(u(1,1)); % copy index's to d vector<br /> % $$$ d(2,1)=floor(u(2,1));<br /> % $$$ d=[d, mat_in(d,1), -mat_in(d,2)];<br /> % $$$ disp('Your final cast_id array:')<br /> % $$$ disp(d)<br /> % $$$ save bt.pts d /ascii;<br /> % $$$ else<br /> % $$$ disp('castID not saved, re-run by typing &quot;bt_ac9&quot; at the prompt');<br /> % $$$ end % end if-else<br /> ButtonName=questdlg('Save these points to CastID matrix?',...<br /> '','Yes','No','No');<br /> switch ButtonName,<br /> case 'Yes'<br /> %d=ceil(u(:,1));<br /> d=[];<br /> for id=1:size(u,1),<br /> dd=find(mat_in(:,1)&lt;u(id,1));<br /> d=[d; dd(end)];<br /> end<br /> %d(1,1)=ceil(u(1,1)) % copy index's to d vector<br /> %d(2,1)=floor(u(2,1))<br /> case 'No'<br /> disp('Points not saved to castID, ');<br /> for ip=1:2:size(d,1),<br /> plot(mat_in(d(ip):d(ip+1),1),...<br /> mat_in(d(ip):d(ip+1),2));<br /> end <br /> u(length(u)-1:length(u),:)=[];<br /> %d=ceil(u(:,1));<br /> d=[];<br /> for id=1:size(u,1),<br /> dd=find(mat_in(:,1)&lt;u(id,1));<br /> d=[d; dd(end)];<br /> end<br /> end % switch<br /> ButtonName=questdlg('Pick points for another cast?',...<br /> '','Yes','No','No');<br /> switch ButtonName,<br /> case 'Yes'<br /> case 'No'<br /> if ~isempty(d),<br /> d=[d, mat_in(d,1), -mat_in(d,2)];<br /> disp('Your final cast_id array:')<br /> disp(d)<br /> %cd(pb) %--------------------------<br /> save bt.pts d /ascii;<br /> % cd ..<br /> end<br /> return;<br /> end % switch<br /> <br /> end % while <br /> return<br /> &lt;/pre&gt;<br /> <br /> &lt;big&gt;[[bt_ac9]]&lt;big&gt;</div>

128.111.101.185