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
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
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');
79
80%plot existing profile if any
81prevplot=1;
82prevplot2=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
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;
108 end
109 end
110 A(list)=[];
111 closed(list)=[];
112
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
130
131
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};
142 end
143 end
144
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};
155 end
156 end
157
158 switch button
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
171 [A,numprofiles,numpoints,closed]=addcontour(A,numprofiles,numpoints,closed,prevplot2,root);
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
180 [A,numprofiles,numpoints,closed]=removeprofile(A,numprofiles,numpoints,closed,prevplot2,root);
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
189 [A,numprofiles,numpoints,closed]=modifyposition(A,numprofiles,numpoints,closed,prevplot,root);
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
198 [A,numprofiles,numpoints,closed]=addinsideprofile(A,numprofiles,numpoints,closed,prevplot,root);
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
207 [A,numprofiles,numpoints,closed]=addendprofile(A,numprofiles,numpoints,closed,prevplot2,root);
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
216 [A,numprofiles,numpoints,closed]=removepoints(A,numprofiles,numpoints,closed,prevplot,root);
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
225 [A,numprofiles,numpoints,closed]=removeseveralpoints(A,numprofiles,numpoints,closed,prevplot,root);
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
234 [A,numprofiles,numpoints,closed]=cutprofile(A,numprofiles,numpoints,closed,prevplot,root);
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
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
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
260 case 12
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
269 %QUIT
270 case 15
271
272 loop=0;
273
274 otherwise
275
276 %do nothing
277
278 end
279
280 %Now erase all that have been done and plot the new structure A as it is
281 undoplots(prevplot);
282 if numprofiles
283 prevplot2=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
294title('New file written, exiting...','FontSize',14);
295if isempty(A)
296 disp('Profile empty, no file written')
297else
298 expwrite(A,newfile);
299end
300
301%close window
302close;
Note: See TracBrowser for help on using the repository browser.