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

Last change on this file since 285 was 285, checked in by Mathieu Morlighem, 16 years ago

Improved expmaster (if no current figure)
File size: 7.3 KB
Line 
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=[];
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);
49 end
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
59 end
60 end
61end
62
63%Get root of newfile
64[path root ext ver]=fileparts(newfile);
65
66%get current figure
67if ~isempty(get(0,'children')),%if there is already a figure (return the number of opened figures)
68 set(gcf,'Renderer','zbuffer'); %fixes a bug on Mac OS X (not needed in future Matlab version)
69 F=getframe(gca);
70 F=F.cdata;
71 %get current axis
72 xlim=get(gca,'Xlim');
73 ylim=get(gca,'Ylim');
74 %recreate x_m and y_m
75 x_m=linspace(xlim(1),xlim(2),size(F,2));
76 y_m=linspace(ylim(2),ylim(1),size(F,1)); %getframe reverse axis...
77 %plot the data in another figure
78 figure
79 imagesc(x_m,y_m,F); set(gca,'Ydir','normal');
80else
81 figure
82end
83
84%plot existing profile if any
85prevplot=1;
86prevplot2=1;
87hold on
88if numprofiles
89 for i=1:numprofiles
90 plot(A(i).x,A(i).y,'-r','MarkerSize',10);
91 max(A(i).x)
92 end
93end
94
95%Build backup structre for do and redo
96backup=cell(1,3);
97backup{1,1}=A;
98backup{1,2}=numprofiles;
99backup{1,3}=numpoints;
100backup{1,4}=closed;
101
102loop=1;
103counter=1;
104while loop
105
106 %Go through A and rule out the empty profiles
107 list=[];
108 for i=1:size(A,2);
109 if length(A(i).x)==0
110 list(end+1)=i;
111 numprofiles=numprofiles-1;
112 end
113 end
114 A(list)=[];
115 closed(list)=[];
116
117 %display menu
118 title('Main Menu','FontSize',14);
119 button=menu('Menu','add a profile (open)',...%1
120 'add a contour (closed)',... %2
121 'remove a profile',... %3
122 'modify the position of a point',... %4
123 'add points inside a profile',... %5
124 'add points at the end of a profile',... %6
125 'remove points',... %7
126 'remove several points',... %8
127 'cut a segment',... %9
128 'cut a large area',... %10
129 'merge profiles',... %11
130 'close profile',... %12
131 'undo',... %13
132 'redo',... %14
133 'quit'); %15
134
135
136 %UNDO??
137 if button==13;
138 if counter==1
139 disp('Already at oldest change');
140 else
141 counter=counter-1;
142 A=backup{counter,1};
143 numprofiles=backup{counter,2};
144 numpoints=backup{counter,3};
145 closed=backup{counter,4};
146 end
147 end
148
149 %REDO??
150 if button==14
151 if counter==size(backup,1)
152 disp('Already at newest change');
153 else
154 counter=counter+1;
155 A=backup{counter,1};
156 numprofiles=backup{counter,2};
157 numpoints=backup{counter,3};
158 closed=backup{counter,4};
159 end
160 end
161
162 switch button
163
164 case 1
165
166 [A,numprofiles,numpoints,closed]=addprofile(A,numprofiles,numpoints,closed,prevplot2,root);
167 counter=counter+1;
168 backup{counter,1}=A;
169 backup{counter,2}=numprofiles;
170 backup{counter,3}=numpoints;
171 backup{counter,4}=closed;
172
173 case 2
174
175 [A,numprofiles,numpoints,closed]=addcontour(A,numprofiles,numpoints,closed,prevplot2,root);
176 counter=counter+1;
177 backup{counter,1}=A;
178 backup{counter,2}=numprofiles;
179 backup{counter,3}=numpoints;
180 backup{counter,4}=closed;
181
182 case 3
183
184 [A,numprofiles,numpoints,closed]=removeprofile(A,numprofiles,numpoints,closed,prevplot2,root);
185 counter=counter+1;
186 backup{counter,1}=A;
187 backup{counter,2}=numprofiles;
188 backup{counter,3}=numpoints;
189 backup{counter,4}=closed;
190
191 case 4
192
193 [A,numprofiles,numpoints,closed]=modifyposition(A,numprofiles,numpoints,closed,prevplot,root);
194 counter=counter+1;
195 backup{counter,1}=A;
196 backup{counter,2}=numprofiles;
197 backup{counter,3}=numpoints;
198 backup{counter,4}=closed;
199
200 case 5
201
202 [A,numprofiles,numpoints,closed]=addinsideprofile(A,numprofiles,numpoints,closed,prevplot,root);
203 counter=counter+1;
204 backup{counter,1}=A;
205 backup{counter,2}=numprofiles;
206 backup{counter,3}=numpoints;
207 backup{counter,4}=closed;
208
209 case 6
210
211 [A,numprofiles,numpoints,closed]=addendprofile(A,numprofiles,numpoints,closed,prevplot2,root);
212 counter=counter+1;
213 backup{counter,1}=A;
214 backup{counter,2}=numprofiles;
215 backup{counter,3}=numpoints;
216 backup{counter,4}=closed;
217
218 case 7
219
220 [A,numprofiles,numpoints,closed]=removepoints(A,numprofiles,numpoints,closed,prevplot,root);
221 counter=counter+1;
222 backup{counter,1}=A;
223 backup{counter,2}=numprofiles;
224 backup{counter,3}=numpoints;
225 backup{counter,4}=closed;
226
227 case 8
228
229 [A,numprofiles,numpoints,closed]=removeseveralpoints(A,numprofiles,numpoints,closed,prevplot,root);
230 counter=counter+1;
231 backup{counter,1}=A;
232 backup{counter,2}=numprofiles;
233 backup{counter,3}=numpoints;
234 backup{counter,4}=closed;
235
236 case 9
237
238 [A,numprofiles,numpoints,closed]=cutprofile(A,numprofiles,numpoints,closed,prevplot,root);
239 counter=counter+1;
240 backup{counter,1}=A;
241 backup{counter,2}=numprofiles;
242 backup{counter,3}=numpoints;
243 backup{counter,4}=closed;
244
245 case 10
246
247 [A,numprofiles,numpoints,closed]=cutarea(A,numprofiles,numpoints,closed,prevplot,root);
248 counter=counter+1;
249 backup{counter,1}=A;
250 backup{counter,2}=numprofiles;
251 backup{counter,3}=numpoints;
252 backup{counter,4}=closed;
253
254 case 11
255
256 [A,numprofiles,numpoints,closed]=mergeprofiles(A,numprofiles,numpoints,closed,prevplot,root);
257 counter=counter+1;
258 backup{counter,1}=A;
259 backup{counter,2}=numprofiles;
260 backup{counter,3}=numpoints;
261 backup{counter,4}=closed;
262
263
264 case 12
265
266 [A,numprofiles,numpoints,closed]=closeprofile(A,numprofiles,numpoints,closed,prevplot,root);
267 counter=counter+1;
268 backup{counter,1}=A;
269 backup{counter,2}=numprofiles;
270 backup{counter,3}=numpoints;
271 backup{counter,4}=closed;
272
273 %QUIT
274 case 15
275
276 loop=0;
277
278 otherwise
279
280 %do nothing
281
282 end
283
284 %Now erase all that have been done and plot the new structure A as it is
285 undoplots(prevplot);
286 if numprofiles
287 prevplot2=1;
288 for i=1:numprofiles
289 plot(A(i).x,A(i).y,'-r','MarkerSize',10);
290 prevplot2=prevplot2+1;
291 end
292 end
293end
294
295hold off
296
297%write contour using expwrite
298title('New file written, exiting...','FontSize',14);
299if isempty(A)
300 disp('Profile empty, no file written')
301else
302 expwrite(A,newfile);
303end
304
305%close window
306close;
Note: See TracBrowser for help on using the repository browser.