source: issm/trunk/src/m/utils/Exp/expmaster.m@ 33

Last change on this file since 33 was 33, checked in by seroussi, 16 years ago

initial input
File size: 7.2 KB
RevLine 
[1]1function expmaster(newfile,varargin)
2%EXPMASTER - allow to create, modify, add, cut, .. segments of domain outline together
3%
4% this routine is used to create, modify, cut,... an Argus file (.exp)
5%
6% expmaster(newprofile)
7% creation of an argus file newprofile
8%
9% expmaster(newprofile,oldprofile)
10% the modification of the file oldprofile will be saved in newprofile
11%
12% expmaster(newprofile,oldprofile1,oldprofile2,...)
13% the modification of the files oldprofile* will be saved in newprofile
14%
15% Usage:
16% expmaster(newfile,varargin)
17%
18% See also EXPDOC
19
20%Some checks
21if ~nargin | nargout
22 error('expmaster usage: expmaster(newfile,varargin)')
23elseif exist(newfile),
24 choice=input(['A file ' newfile ' already exists, do you want to modify it? (y/n)'],'s');
25 if ~strcmpi(choice,'y'),
26 error('no modification done ... exiting');
27 end
28end
29
30%put all the argus profiles given in input in one structure A
31A=struct([]);
32numprofiles=0;
33numpoints=0;
34closed=[];
[33]35for i=1:nargin-1
36 filename=varargin{i};
37 if ~exist(filename),
38 error(['expmaster error message:, ' filename ' does not exist. Exiting...']);
39 else
40 %read file
41 B=expread(filename,1);
42 %go through all profiles of B
43 for i=1:size(B,2)
44 %plug profile in A
45 if numprofiles
46 A(numprofiles+1)=B(i);
47 else
48 A=B(i);
[1]49 end
[33]50 %update numprofiles and numpoints
51 numpoints=numpoints+length(B(i).x);
52 numprofiles=numprofiles+1;
53 %figure out if the profile is closed or not
54 if (B(i).x(1)==B(i).x(end) & B(i).y(1)==B(i).y(end))
55 closed(numprofiles)=1;
56 else
57 closed(numprofiles)=0;
58 end
[1]59 end
60 end
61end
62
63%Get root of newfile
64[path root ext ver]=fileparts(newfile);
65
[33]66%get current figure
67set(gcf,'Renderer','zbuffer'); %fixes a bug on Mac OS X (not needed in future Matlab version
68F=getframe(gca);
69F=F.cdata;
70%get current axis
71xlim=get(gca,'Xlim');
72ylim=get(gca,'Ylim');
73%recreate x_m and y_m
74x_m=linspace(xlim(1),xlim(2),size(F,2));
75y_m=linspace(ylim(2),ylim(1),size(F,1)); %getframe reverse axis...
76%plot the data in another figure
77figure
78imagesc(x_m,y_m,F); set(gca,'Ydir','normal');
[1]79
80%plot existing profile if any
[33]81prevplot=1;
82prevplot2=1;
[1]83hold on
84if numprofiles
85 for i=1:numprofiles
86 plot(A(i).x,A(i).y,'-r','MarkerSize',10);
87 prevplot2=prevplot2+1;
88 end
89end
90
91%Build backup structre for do and redo
92backup=cell(1,3);
93backup{1,1}=A;
94backup{1,2}=numprofiles;
95backup{1,3}=numpoints;
96backup{1,4}=closed;
97
98loop=1;
99counter=1;
100while loop
101
[33]102 %Go through A and rule out the empty profiles
103 list=[];
104 for i=1:size(A,2);
105 if length(A(i).x)==0
106 list(end+1)=i;
107 numprofiles=numprofiles-1;
[1]108 end
[33]109 end
110 A(list)=[];
111 closed(list)=[];
[1]112
[33]113 %display menu
114 title('Main Menu','FontSize',14);
115 button=menu('Menu','add a profile (open)',...%1
116 'add a contour (closed)',... %2
117 'remove a profile',... %3
118 'modify the position of a point',... %4
119 'add points inside a profile',... %5
120 'add points at the end of a profile',... %6
121 'remove points',... %7
122 'remove several points',... %8
123 'cut a segment',... %9
124 'cut a large area',... %10
125 'merge profiles',... %11
126 'close profile',... %12
127 'undo',... %13
128 'redo',... %14
129 'quit'); %15
[1]130
131
[33]132 %UNDO??
133 if button==13;
134 if counter==1
135 disp('Already at oldest change');
136 else
137 counter=counter-1;
138 A=backup{counter,1};
139 numprofiles=backup{counter,2};
140 numpoints=backup{counter,3};
141 closed=backup{counter,4};
[1]142 end
[33]143 end
[1]144
[33]145 %REDO??
146 if button==14
147 if counter==size(backup,1)
148 disp('Already at newest change');
149 else
150 counter=counter+1;
151 A=backup{counter,1};
152 numprofiles=backup{counter,2};
153 numpoints=backup{counter,3};
154 closed=backup{counter,4};
[1]155 end
[33]156 end
[1]157
[33]158 switch button
[1]159
160 case 1
161
162 [A,numprofiles,numpoints,closed]=addprofile(A,numprofiles,numpoints,closed,prevplot2,root);
163 counter=counter+1;
164 backup{counter,1}=A;
165 backup{counter,2}=numprofiles;
166 backup{counter,3}=numpoints;
167 backup{counter,4}=closed;
168
169 case 2
170
[33]171 [A,numprofiles,numpoints,closed]=addcontour(A,numprofiles,numpoints,closed,prevplot2,root);
[1]172 counter=counter+1;
173 backup{counter,1}=A;
174 backup{counter,2}=numprofiles;
175 backup{counter,3}=numpoints;
176 backup{counter,4}=closed;
177
178 case 3
179
[33]180 [A,numprofiles,numpoints,closed]=removeprofile(A,numprofiles,numpoints,closed,prevplot2,root);
[1]181 counter=counter+1;
182 backup{counter,1}=A;
183 backup{counter,2}=numprofiles;
184 backup{counter,3}=numpoints;
185 backup{counter,4}=closed;
186
187 case 4
188
[33]189 [A,numprofiles,numpoints,closed]=modifyposition(A,numprofiles,numpoints,closed,prevplot,root);
[1]190 counter=counter+1;
191 backup{counter,1}=A;
192 backup{counter,2}=numprofiles;
193 backup{counter,3}=numpoints;
194 backup{counter,4}=closed;
195
196 case 5
197
[33]198 [A,numprofiles,numpoints,closed]=addinsideprofile(A,numprofiles,numpoints,closed,prevplot,root);
[1]199 counter=counter+1;
200 backup{counter,1}=A;
201 backup{counter,2}=numprofiles;
202 backup{counter,3}=numpoints;
203 backup{counter,4}=closed;
204
205 case 6
206
[33]207 [A,numprofiles,numpoints,closed]=addendprofile(A,numprofiles,numpoints,closed,prevplot2,root);
[1]208 counter=counter+1;
209 backup{counter,1}=A;
210 backup{counter,2}=numprofiles;
211 backup{counter,3}=numpoints;
212 backup{counter,4}=closed;
213
214 case 7
215
[33]216 [A,numprofiles,numpoints,closed]=removepoints(A,numprofiles,numpoints,closed,prevplot,root);
[1]217 counter=counter+1;
218 backup{counter,1}=A;
219 backup{counter,2}=numprofiles;
220 backup{counter,3}=numpoints;
221 backup{counter,4}=closed;
222
223 case 8
224
[33]225 [A,numprofiles,numpoints,closed]=removeseveralpoints(A,numprofiles,numpoints,closed,prevplot,root);
[1]226 counter=counter+1;
227 backup{counter,1}=A;
228 backup{counter,2}=numprofiles;
229 backup{counter,3}=numpoints;
230 backup{counter,4}=closed;
231
232 case 9
233
[33]234 [A,numprofiles,numpoints,closed]=cutprofile(A,numprofiles,numpoints,closed,prevplot,root);
[1]235 counter=counter+1;
236 backup{counter,1}=A;
237 backup{counter,2}=numprofiles;
238 backup{counter,3}=numpoints;
239 backup{counter,4}=closed;
240
241 case 10
242
[33]243 [A,numprofiles,numpoints,closed]=cutarea(A,numprofiles,numpoints,closed,prevplot,root);
244 counter=counter+1;
245 backup{counter,1}=A;
246 backup{counter,2}=numprofiles;
247 backup{counter,3}=numpoints;
248 backup{counter,4}=closed;
249
250 case 11
251
[1]252 [A,numprofiles,numpoints,closed]=mergeprofiles(A,numprofiles,numpoints,closed,prevplot,root);
253 counter=counter+1;
254 backup{counter,1}=A;
255 backup{counter,2}=numprofiles;
256 backup{counter,3}=numpoints;
257 backup{counter,4}=closed;
258
259
[33]260 case 12
[1]261
262 [A,numprofiles,numpoints,closed]=closeprofile(A,numprofiles,numpoints,closed,prevplot,root);
263 counter=counter+1;
264 backup{counter,1}=A;
265 backup{counter,2}=numprofiles;
266 backup{counter,3}=numpoints;
267 backup{counter,4}=closed;
268
[33]269 %QUIT
270 case 15
[1]271
272 loop=0;
273
274 otherwise
275
276 %do nothing
277
[33]278 end
[1]279
[33]280 %Now erase all that have been done and plot the new structure A as it is
281 undoplots(prevplot);
[1]282 if numprofiles
[33]283 prevplot2=1;
[1]284 for i=1:numprofiles
285 plot(A(i).x,A(i).y,'-r','MarkerSize',10);
286 prevplot2=prevplot2+1;
287 end
288 end
289end
290
291hold off
292
293%write contour using expwrite
[33]294title('New file written, exiting...','FontSize',14);
[1]295if isempty(A)
296 disp('Profile empty, no file written')
297else
298 expwrite(A,newfile);
299end
[33]300
301%close window
302close;
Note: See TracBrowser for help on using the repository browser.