Revisión | e2e49dc4bd4327c897d866888f985af01c36e57e (tree) |
---|---|
Tiempo | 2014-11-04 04:39:28 |
Autor | Takuo Yasunaga <yasunaga@bio....> |
Commiter | Takuo Yasunaga |
Bug fixed:
Eos_env
install-eos.sh
@@ -1,498 +0,0 @@ | ||
1 | -.SUFFIXES: .tif .dat .prectf .fft .ctf .ctfinfoauto .ctfinfo .mrc .divideroi .divideroiinfo \ | |
2 | - .box .roi \ | |
3 | - .roiinfo-coarse .roiall-coarse .roiinforef-coarse .roiallref-coarse .roiinfoprev-coarse \ | |
4 | - .roiinfo .roiall .roiinforef .roiallref .roiinfoprev \ | |
5 | - .box2roiinfo .roiinfo2box .corinfo .divideroiinfo \ | |
6 | - .roilst0 .roilst .mrc3d .ds6 .mrc3dcentremed .mrc3dcentre .mrc3dcentremedpad .orimrccentre .oricorinfo .coarroi .coarroiinfo .roi2 .pad .mul | |
7 | - | |
8 | - | |
9 | -#comment#Forall | |
10 | -TILTAXISNUMBER=2 | |
11 | -ORIFILE1=a | |
12 | -ORIFILE2=b | |
13 | -ZEROTILTNUM=60 | |
14 | -MINTILTNUM=0 | |
15 | -MAXTILTNUM=120 | |
16 | -STEPTILTNUM=4 | |
17 | -T=43 | |
18 | - | |
19 | -#comment#Fordata2mrc | |
20 | -RESOLUTION=1.7 | |
21 | - | |
22 | -#comment#For2ndAxisRotate | |
23 | -ROTATIONANGLE=-90 | |
24 | - | |
25 | -#comment#ForCOARCreate | |
26 | -ORISIZEX=2048 | |
27 | -ORISIZEY=2048 | |
28 | -CORMODE=16 | |
29 | -ORINX=1024 | |
30 | -ORINY=1024 | |
31 | -ORINZ=1 | |
32 | - | |
33 | -#comment#Forroiinfo | |
34 | -ORIBASE=$(ORIFILE1)0$(ZEROTILTNUM) | |
35 | - | |
36 | -#comment#ForPadAndMul | |
37 | -TILTAXISANGLE=0 | |
38 | -IMAGEPADm=15 | |
39 | -IMAGEPADV=0 | |
40 | -IMAGEMULTIPLYINGv=-1 | |
41 | - | |
42 | -#comment#ForroilstCreate | |
43 | -EULARANGLE=ZENS | |
44 | -MINTILTANGLE=-60 | |
45 | -STEPTILTANGLE=4 | |
46 | - | |
47 | -#comment#Fords6Create | |
48 | -MRC2MAPmul=1 | |
49 | -SMOOTHr=5 | |
50 | -SMOOTHm=1 | |
51 | -SMOOTHM=1 | |
52 | -CENTERNX=200 | |
53 | -CENTERNY=200 | |
54 | -CENTERNZ=200 | |
55 | - | |
56 | -all: data2mrc Smoothing 2ndAxisRotate COARCreate $(ORIBASE).roiinfo BoxCreate PadAndMul roilstCreate ds6Create | |
57 | -data2mrc:: | |
58 | - for (( j=1 ; $$j <= $(TILTAXISNUMBER) ; j=$$(($$j + 1 )))) \ | |
59 | - do \ | |
60 | - if [ $$j -eq 1 ]; then \ | |
61 | - orifile=$(ORIFILE1); \ | |
62 | - else \ | |
63 | - orifile=$(ORIFILE2); \ | |
64 | - fi ; \ | |
65 | - for (( i=0 ; $$i <= $(MAXTILTNUM) ; i=$$(($$i + $(STEPTILTNUM))))) \ | |
66 | - do \ | |
67 | - ori=`printf "%s%03d" $$orifile $$i`; \ | |
68 | - make $$ori.prectf; \ | |
69 | - make $$ori.divideroiinfo; \ | |
70 | - awk '{ printf("%s\n", $$1); }' $$ori.divideroiinfo >> dvdroilst ; \ | |
71 | - for divideroivaliable in `ls -1 $$ori*.divideroi` ; do \ | |
72 | - dividefilename=`basename $$divideroivaliable .divideroi`; \ | |
73 | - make $$dividefilename.fft; \ | |
74 | - make $$dividefilename.ctfinfo; \ | |
75 | - make $$dividefilename.ctf; \ | |
76 | - done ; \ | |
77 | - make $$ori.mrc ; \ | |
78 | - done ; \ | |
79 | - done | |
80 | - | |
81 | -.dat.prectf: | |
82 | - emData2mrc -i $*.dat -o $*.prectf | |
83 | - | |
84 | -.tif.prectf: | |
85 | - tiff2mrc -i $*.tif -o $*.prectf -m 1 -r $(RESOLUTION) | |
86 | - | |
87 | -.prectf.divideroiinfo: | |
88 | - mrcImageDivideInfoForTomography -i $*.prectf -o $*.divideroiinfo -m 1 -d 4 | |
89 | - sed 's/roi/divideroi/' $*.divideroiinfo >> $*.tmproiinfo | |
90 | - mv $*.tmproiinfo $*.divideroiinfo | |
91 | - mrcImageROIs -i $*.prectf -o $*.mon -I $*.divideroiinfo -S 1 | |
92 | - | |
93 | -.divideroi.fft: | |
94 | - mrcImagePad -i $*.divideroi -o $*.dividepad -W 1024 -H 1024 -m 3 ; \ | |
95 | - mrcImageAbnormalValueRemove -i $*.dividepad -o $*.prectfab -m 1 ; \ | |
96 | - mrcImageFFT -i $*.prectfab -o $*.fft ; \ | |
97 | - | |
98 | -.fft.ctfinfo: | |
99 | - ctfDetermine -i $*.fft -o $*.ctfinfo -D 10000 -m 128 -CutLow 0.05 -CutHigh 0.2 -Ac 200 -d 3 -Cc 0.7 -Cs 1.7 -ratio 0.2 | grep defocus > CTFD ; \ | |
100 | - awk -v BASE=$* '{ printf("%s %s\n", BASE, $$2); }' CTFD >> $*.defocuslst ; \ | |
101 | - | |
102 | -.fft.ctf: | |
103 | - mrcImageCTFCompensation -i $*.fft -o $*.ctf -info2 $*.ctfinfo -m 1 ; \ | |
104 | - mrcImageCenterGet -i $*.ctf -o $*.ctfcentre -Nx 512 -Ny 512 ; \ | |
105 | - | |
106 | -.prectf.mrc: | |
107 | - cat $*.divideroiinfo | sed s/divideroi/ctfcentre/g > $*.divideroiinfo2 | |
108 | - mrcImageCTFCompensationForTiltImage -i $*.prectf -I $*.divideroiinfo2 -o $*.mrc | |
109 | - | |
110 | -#.fft.ctfinfoauto: | |
111 | -# ctfDetermine -i $*.fft -o $*.ctfinfoauto -D 50000 -m 163 -CutLow 0.04 -CutHigh 0.15 -Ac 200 -d 3 -Cc 0.7 -Cs 2.1 | |
112 | -## awk '/^[^:]/ {print $$0};' $*.ctfinfoauto > $*.ctfinfo | |
113 | -## Display2 -i $*.fft -H 1e7 | |
114 | - | |
115 | -#.ctfinfoauto.ctfinfo: | |
116 | -# if [ ! -f $*.ctfinfo ] ; then \ | |
117 | - if [ -f default.ctfinfo ] ; then \ | |
118 | - cat default.ctfinfo >> $*.ctfinfo ; \ | |
119 | - else \ | |
120 | - cat $*.ctfinfoauto >> $*.ctfinfo ; \ | |
121 | - fi \ | |
122 | - fi | |
123 | -# ctfDisplay -i $*.fft | |
124 | -# mv $@ $@.tmp | |
125 | -# grep .: $@.tmp > $@ | |
126 | -# cp $@ default.ctfinfo | |
127 | - | |
128 | -#.ctfinfo.mrc: | |
129 | -# mv $*.ctfinfo $*.ctfinfotmp | |
130 | -# awk '/^[^:]/ {print $$0};' $*.ctfinfotmp > $*.ctfinfo | |
131 | -# mrcImageCTFCompensation -i $*.fft -o $*.mrc -info2 $*.ctfinfo -m 1 | |
132 | -# Display2 -i $*.mrc -geometry 1024x1024+0+0 -Zoom 0.5 & | |
133 | - | |
134 | -mrcImageSmoothing:: | |
135 | - for (( j=1 ; $$j <= $(TILTAXISNUMBER) ; j=$$(($$j + 1 )))) \ | |
136 | - do \ | |
137 | - if [ $$j -eq 1 ]; then \ | |
138 | - orifile=$(ORIFILE1); \ | |
139 | - else \ | |
140 | - orifile=$(ORIFILE2); \ | |
141 | - fi ; \ | |
142 | - for (( i=0 ; $$i <= $(MAXTILTNUM) ; i=$$(($$i + $(STEPTILTNUM))))) \ | |
143 | - do \ | |
144 | - ori=`printf "%s%03d" $$orifile $$i`; \ | |
145 | - mrcImageSmoothing -i $$ori.mrc -o $$ori.smth -sx 5 -sy 5 -m 1 ; \ | |
146 | - mv $$ori.smth $$ori.mrc ; \ | |
147 | - done ; \ | |
148 | - done | |
149 | - | |
150 | -mrcImageWindowing:: | |
151 | - for (( j=1 ; $$j <= $(TILTAXISNUMBER) ; j=$$(($$j + 1 )))) \ | |
152 | - do \ | |
153 | - if [ $$j -eq 1 ]; then \ | |
154 | - orifile=$(ORIFILE1); \ | |
155 | - else \ | |
156 | - orifile=$(ORIFILE2); \ | |
157 | - fi ; \ | |
158 | - for (( i=0 ; $$i <= $(MAXTILTNUM) ; i=$$(($$i + $(STEPTILTNUM))))) \ | |
159 | - do \ | |
160 | - ori=`printf "%s%03d" $$orifile $$i`; \ | |
161 | - mrcImageWindowing -i $$ori.mrc -o $$ori.wnd ; \ | |
162 | - mv $$ori.wnd $$ori.mrc ; \ | |
163 | - done ; \ | |
164 | - done | |
165 | - | |
166 | -2ndAxisRotate:: | |
167 | - for (( i=0 ; $$i <= $(MAXTILTNUM) ; i=$$(( $$i + $(STEPTILTNUM) )) )) \ | |
168 | - do \ | |
169 | - ori=`printf "%s%03d" $(ORIFILE2) $$i`; \ | |
170 | - mrcImageRotation -i $$ori.mrc -o $$ori.mrc2 -a $(ROTATIONANGLE) ; \ | |
171 | - mv $$ori.mrc2 $$ori.mrc ; \ | |
172 | - done | |
173 | - | |
174 | -##sep | |
175 | - | |
176 | -COARCreate:: | |
177 | - for (( j=1 ; $$j <= $(TILTAXISNUMBER) ; j=$$(($$j + 1 )))) \ | |
178 | - do \ | |
179 | - if [ $$j -eq 1 ]; then \ | |
180 | - orifile=$(ORIFILE1); \ | |
181 | - oribase=$(ORIFILE1)0$(ZEROTILTNUM); \ | |
182 | - else \ | |
183 | - orifile=$(ORIFILE2); \ | |
184 | - fi ; \ | |
185 | - ORIZEROFILE=$$oribase.mrc; \ | |
186 | - previous=$$oribase ; \ | |
187 | - echo "Create the zero-tilt Image ROI (same as mrc, just copy)"; \ | |
188 | - make $$oribase.roiinfo-coarse; echo "$$oribase.roi will be created."; \ | |
189 | - make $$oribase.roiall-coarse ; echo "$$oribase.roi was created." ; \ | |
190 | - make $$oribase.orimrccentre ; echo "The centre of $$oribase.roi was created." ; \ | |
191 | - for delta in +$(STEPTILTNUM) -$(STEPTILTNUM) ; do \ | |
192 | - echo "$$delta"; \ | |
193 | - previous=$$oribase; \ | |
194 | - TOTALX=0 ; \ | |
195 | - TOTALY=0 ; \ | |
196 | - echo "Total: $$TOTALX $$TOTALY" > $$previous.oricorinfo; \ | |
197 | - for (( i=$$(($(ZEROTILTNUM) + $$delta )) ; (( $(MINTILTNUM) <= $$i && $$i <= $(MAXTILTNUM) )) ; i=$$(($$i + $$delta)) )) \ | |
198 | - do \ | |
199 | - ori=`printf "%s%03d" $$orifile $$i`; \ | |
200 | - echo $$i, $(MAXTILTNUM), $(MINTILTNUM); \ | |
201 | - echo $$previous > $$ori.prevname ; \ | |
202 | - ln -sf $$previous.roiinfo-coarse $$ori.roiinfoprev-coarse ; make $$ori.roiinforef-coarse ; \ | |
203 | - make $$ori.roiallref-coarse ; \ | |
204 | - make $$ori.orimrccentre ; \ | |
205 | - TOTALX=`awk '/Total:/ { print $$2 }' $$previous.oricorinfo`; \ | |
206 | - TOTALY=`awk '/Total:/ { print $$3 }' $$previous.oricorinfo`; \ | |
207 | - ln -sf $$previous.orimrccentre $$ori.oriref ; make $$ori.oricorinfo ; \ | |
208 | - echo "TotalPrev: $$TOTALX $$TOTALY" >> $$ori.oricorinfo; \ | |
209 | - make $$ori.roiinfo-coarse; \ | |
210 | - make $$ori.roiall-coarse; \ | |
211 | - make $$ori.orimrccentre ; sleep 1; \ | |
212 | - touch $$ori.oricorinfo ; sleep 1; \ | |
213 | - touch $$ori.roiinfo-coarse; sleep 1; \ | |
214 | - touch $$ori.roiall-coarse; sleep 1; \ | |
215 | - previous=$$ori; \ | |
216 | - done ; \ | |
217 | - echo $$i; \ | |
218 | - done ; \ | |
219 | - done | |
220 | - | |
221 | -CleanCOARCreate:: | |
222 | - rm *.roiinfo-coarse | |
223 | - | |
224 | -$(ORIFILE1)0$(ZEROTILTNUM).roiinfo-coarse: | |
225 | - printf "%s%03d.roi Rect2 %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f\n" $(ORIFILE1) $(ZEROTILTNUM) $$(($(ORISIZEX)/2)) $$(($(ORISIZEY)/2)) $(ORISIZEX) $(ORISIZEY) > $(ORIFILE1)0$(ZEROTILTNUM).roiinfo-coarse | |
226 | - | |
227 | -$(ORIFILE2)0$(ZEROTILTNUM).roiinfo-coarse: | |
228 | - printf "%s%03d.roi Rect2 %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f\n" $(ORIFILE2) $(ZEROTILTNUM) $$(($(ORISIZEX)/2)) $$(($(ORISIZEY)/2)) $(ORISIZEX) $(ORISIZEY) > $(ORIFILE2)0$(ZEROTILTNUM).roiinfo-coarse | |
229 | - | |
230 | -.orimrccentre.oricorinfo: | |
231 | - echo "orimrccentre -> oricorinfo: SHIFT:= current.orimrccentre for prev.oirmrccentre" | |
232 | - mrcImageCorrelation -i $*.orimrccentre -r $*.oriref -c $*.oricorinfo -o $*.oricormap -m $(CORMODE) -s $*.orishift | |
233 | - | |
234 | -.oricorinfo.roiinfo-coarse: | |
235 | - TOTALX=`awk '/TotalPrev:/ { print $$2 }' $*.oricorinfo`; \ | |
236 | - TOTALY=`awk '/TotalPrev:/ { print $$3 }' $*.oricorinfo`; \ | |
237 | - SHIFTX=`awk '/Shift:/ { print $$2 }' $*.oricorinfo`; \ | |
238 | - SHIFTY=`awk '/Shift:/ { print $$3 }' $*.oricorinfo`; \ | |
239 | - echo "$$SHIFTX,$$SHIFTY"; \ | |
240 | - if [ $$(($(ORINX)/2)) -lt $$SHIFTX ]; then \ | |
241 | - SHIFTX=$$(($$SHIFTX - $(ORINX))) ; \ | |
242 | - fi ; \ | |
243 | - if [ $$(($(ORINY)/2)) -lt $$SHIFTY ]; then \ | |
244 | - SHIFTY=$$(($$SHIFTY - $(ORINY))) ; \ | |
245 | - fi ; \ | |
246 | - TOTALX=$$(($$TOTALX + $$SHIFTX)) ; \ | |
247 | - TOTALY=$$(($$TOTALY + $$SHIFTY)) ; \ | |
248 | - echo "Total: $$TOTALX $$TOTALY" >> $*.oricorinfo; \ | |
249 | - printf "%s.roi Rect2 %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f\n" $* $$(($(ORISIZEX)/2+$$TOTALX)) $$(($(ORISIZEY)/2+$$TOTALY)) $(ORISIZEX) $(ORISIZEY) $$TOTALX $$TOTALY > $*.roiinfo-coarse ; | |
250 | - | |
251 | -.roiinfo-coarse.roiall-coarse: | |
252 | - make $*.mrc | |
253 | - mrcImageROIs -i $*.mrc -I $*.roiinfo-coarse | |
254 | - touch $*.roiall-coarse | |
255 | - | |
256 | -.roiinforef-coarse.roiallref-coarse: | |
257 | - mrcImageROIs -i $*.mrc -I $*.roiinforef-coarse | |
258 | - touch $*.roiallref-coarse | |
259 | - | |
260 | -.roiinfoprev-coarse.roiinforef-coarse: | |
261 | - previous=`cat $*.prevname`; \ | |
262 | - sed -s s/$$previous/$*/ $*.roiinfoprev-coarse > $*.roiinforef-coarse | |
263 | - | |
264 | -.roi.orimrccentre: | |
265 | - echo "roi -> orimrccentre" | |
266 | - mrcImageCenterGet -i $*.roi -o $*.orimrccentre -Nx $(ORINX) -Ny $(ORINY) -Nz $(ORINZ) | |
267 | - mrcImageWindowing -i $*.orimrccentre -o $*.wnd -m 16 | |
268 | - mv $*.wnd $*.orimrccentre | |
269 | - | |
270 | -$(ORIBASE).roiinfo: | |
271 | - Display2 -i $*.mrc ; | |
272 | -# if [ $(TILTAXISNUMBER) -eq 2 ]; then \ | |
273 | -# Display2 -i $(ORIFILE2)0$(ZEROTILTNUM).mrc ;\ | |
274 | -# fi; \ | |
275 | - | |
276 | -##sep | |
277 | - | |
278 | -BoxCreate:: | |
279 | - for (( j=1 ; $$j <= $(TILTAXISNUMBER) ; j=$$(($$j + 1 )))) \ | |
280 | - do \ | |
281 | - if [ $$j -eq 1 ]; then \ | |
282 | - orifile=$(ORIFILE1); \ | |
283 | - oribase=$(ORIFILE1)0$(ZEROTILTNUM); \ | |
284 | - else \ | |
285 | - orifile=$(ORIFILE2); \ | |
286 | - fi ; \ | |
287 | - BASE=$$oribase; \ | |
288 | - make $$BASE.roiall ; \ | |
289 | - for delta in +$(STEPTILTNUM) -$(STEPTILTNUM) ; \ | |
290 | - do \ | |
291 | - echo $$delta,$$BASE ; \ | |
292 | - previous=$$BASE ; \ | |
293 | - for (( i=$$(($(ZEROTILTNUM)+$$delta)) ; (( $(MINTILTNUM) <= $$i && $$i <= $(MAXTILTNUM) )) ; i=$$(($$i + $$delta)) )) \ | |
294 | - do \ | |
295 | - name=`printf "%s%03d" $$orifile $$i`; \ | |
296 | - echo ">>>> $$name $$i start: $(ZEROTILTNUM) $$(($(ZEROTILTNUM)+$$delta))"; \ | |
297 | - echo $$previous > $$name.prevname; \ | |
298 | - ln -sf $$previous.roiinfo $$name.roiinfoprev ; make $$name.roiinforef ; \ | |
299 | - make $$name.roiinfo ; \ | |
300 | - make $$name.roiall ; \ | |
301 | - previous=$$name ; \ | |
302 | - done ; \ | |
303 | - done ; \ | |
304 | - done | |
305 | - | |
306 | -CleanBoxCreate:: | |
307 | - rm *.roiall *.roiinfoprev *.roiinforef | |
308 | - | |
309 | -.roiinfo.roiall: | |
310 | - make $*.roi | |
311 | - mrcImageROIs -i $*.roi -I $*.roiinfo | |
312 | - touch $*.roiall | |
313 | - | |
314 | -.roiinforef.roiallref: | |
315 | - mrcImageROIs -i $*.roi -I $*.roiinforef | |
316 | - touch $*.roiallref | |
317 | - | |
318 | -.roi.corinfo: | |
319 | - mrcImageCorrelation -i $*.roi -r $*.ref -c $*.corinfo | |
320 | - | |
321 | -.roiinforef.roiinfo: | |
322 | - # Extract ROIs using roiinfo of previous image (roiinforef) | |
323 | - previous=`cat $*.prevname`; \ | |
324 | - echo $$previous; \ | |
325 | - echo "start initial ROI"; \ | |
326 | - make $*.roiallref ; \ | |
327 | - ROINUM=`wc -l $*.roiinforef | awk '{print $$1}'` ; \ | |
328 | - echo $$ROINUM ; \ | |
329 | - rm -f $*.roiinfo ; \ | |
330 | - for (( j=0 ; $$j < $$ROINUM ; j=$$(($$j+1)))) \ | |
331 | - do \ | |
332 | - input=`printf "%s-%04d" $* $$j`; \ | |
333 | - refer=`printf "%s-%04d" $$previous $$j`; \ | |
334 | - echo $$input, $$refer ; \ | |
335 | - ln -sf $$refer.roi $$input.ref ; make $$input.corinfo ; \ | |
336 | - SHIFTX=`awk '/Shift:/ { print $$2 }' $$input.corinfo`; \ | |
337 | - SHIFTY=`awk '/Shift:/ { print $$3 }' $$input.corinfo`; \ | |
338 | - echo "$$SHIFTX,$$SHIFTY,$$input"; \ | |
339 | - awk -v SHIFTXAWK0=$$SHIFTX -v SHIFTYAWK0=$$SHIFTY -v NAMEAWK=$$input -v BASEAWK=$* \ | |
340 | - ' /'"$$input"'/ { \ | |
341 | - SizeX = int(sqrt(($$5-$$3)*($$5-$$3)+( $$6-$$4)*( $$6-$$4))+0.5)+1; \ | |
342 | - SizeY = int(sqrt(($$9-$$3)*($$9-$$3)+($$10-$$4)*($$10-$$4))+0.5)+1; \ | |
343 | - if(SizeX/2<SHIFTXAWK0) { \ | |
344 | - SHIFTXAWK0 -= SizeX ; \ | |
345 | - } \ | |
346 | - if(SizeY/2<SHIFTYAWK0) { \ | |
347 | - SHIFTYAWK0 -= SizeY ; \ | |
348 | - } \ | |
349 | - Angle = atan2($$6-$$4, $$5-$$3) ; \ | |
350 | - SHIFTXAWK = +1*( SHIFTXAWK0*cos(Angle) - SHIFTYAWK0*sin(Angle)); \ | |
351 | - SHIFTYAWK = +1*( SHIFTXAWK0*sin(Angle) + SHIFTYAWK0*cos(Angle)); \ | |
352 | - printf("%s-%04d.roi Rect %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f\n", BASEAWK, FNR-1, $$3+SHIFTXAWK, $$4+SHIFTYAWK, $$5+SHIFTXAWK, $$6+SHIFTYAWK, $$7+SHIFTXAWK, $$8+SHIFTYAWK, $$9+SHIFTXAWK, $$10+SHIFTYAWK, SHIFTXAWK, SHIFTYAWK, Angle*180/3.141592653589793238, SizeX, SizeY);}' \ | |
353 | - $*.roiinforef | tee -a $*.roiinfo ; \ | |
354 | - done ; | |
355 | - | |
356 | -.roiinfoprev.roiinforef: | |
357 | - previous=`cat $*.prevname`; \ | |
358 | - sed -s s/$$previous/$*/ $*.roiinfoprev > $*.roiinforef | |
359 | - | |
360 | -##sep | |
361 | - | |
362 | -PadAndMul:: | |
363 | - for (( k=1 ; $$k <= $(TILTAXISNUMBER) ; k=$$(($$k + 1 )))) \ | |
364 | - do \ | |
365 | - if [ $$k -eq 1 ]; then \ | |
366 | - orifile=$(ORIFILE1); \ | |
367 | - oribase=$(ORIFILE1)0$(ZEROTILTNUM); \ | |
368 | - else \ | |
369 | - orifile=$(ORIFILE2); \ | |
370 | - oribase=$(ORIFILE2)0$(ZEROTILTNUM); \ | |
371 | - fi ; \ | |
372 | - roiNum=`wc -l $$oribase.roiinfo | awk '{print $$1}' `; \ | |
373 | - for ((i=$(MINTILTNUM);i<=$(MAXTILTNUM);i=$$(($$i+$(STEPTILTNUM))) )) \ | |
374 | - do \ | |
375 | - for((j=0 ; j < $$roiNum ; j = $$(($$j+1)))) \ | |
376 | - do \ | |
377 | - name=`printf "%s%03d-%04d" $$orifile $$i $$j` ; \ | |
378 | - echo $$(($$i-$(ZEROTILTNUM))) > $$name.tilt ; \ | |
379 | - ln -sf $$name.roi $$name.roi2 ; make $$name.pad ; \ | |
380 | - make $$name.mul ; \ | |
381 | - done ; \ | |
382 | - done ; \ | |
383 | - done | |
384 | - | |
385 | -CleanPadAndMul:: | |
386 | - rm *.pad | |
387 | - | |
388 | -.roi2.pad: | |
389 | - tilt=`cat $*.tilt` ; \ | |
390 | - echo $$tilt ; \ | |
391 | - mrcImagePad -i $*.roi2 -o $*.pad -Tilt $(TILTAXISANGLE) $$tilt -m $(IMAGEPADm) -V $(IMAGEPADV) | |
392 | - | |
393 | -.pad.mul: | |
394 | - mrcImageMultiplying -i $*.pad -o $*.mul -v $(IMAGEMULTIPLYINGv) | |
395 | - | |
396 | -roilstCreate:: | |
397 | - for (( j=1 ; $$j <= $(TILTAXISNUMBER) ; j=$$(($$j + 1 )))) \ | |
398 | - do \ | |
399 | - if [ $$j -eq 1 ]; then \ | |
400 | - orifile=$(ORIFILE1); \ | |
401 | - oribase=$(ORIFILE1)0$(ZEROTILTNUM); \ | |
402 | - echo "$(ORIFILE1)" > roilstorifile.txt; \ | |
403 | - else \ | |
404 | - orifile=$(ORIFILE2); \ | |
405 | - oribase=$(ORIFILE2)0$(ZEROTILTNUM); \ | |
406 | - echo "$(ORIFILE2)" > roilstorifile.txt; \ | |
407 | - fi ; \ | |
408 | - roiNum=`wc -l $$oribase.roiinfo | awk '{print $$1}' `;\ | |
409 | - echo $$roiNum; \ | |
410 | - for (( i=0 ; $$i < $$roiNum ; i=$$(($$i + 1)))) \ | |
411 | - do \ | |
412 | - tiltnum=`printf "%s%04d" $$orifile $$i`; \ | |
413 | - echo "$$i" > $$tiltnum.roilst0; \ | |
414 | - make $$tiltnum.roilst ; \ | |
415 | - done ; \ | |
416 | - done ; \ | |
417 | - roiNum=`wc -l $$oribase.roiinfo | awk '{print $$1}' `;\ | |
418 | - echo $$roiNum; \ | |
419 | - for (( i=0 ; $$i < $$roiNum ; i=$$(($$i + 1)))) \ | |
420 | - do \ | |
421 | - tiltnum=`printf "%04d" $$i`; \ | |
422 | - echo "$$i" > $$tiltnum.roilst0; \ | |
423 | - cat $(ORIFILE1)$$tiltnum.roilst >> $$tiltnum-d.roilst ; \ | |
424 | - cat $(ORIFILE2)$$tiltnum.roilst >> $$tiltnum-d.roilst ; \ | |
425 | - done | |
426 | - | |
427 | -Cleanroilst:: | |
428 | - rm *.roilst | |
429 | - | |
430 | -.roilst0.roilst: | |
431 | - i=`cat $*.roilst0`; \ | |
432 | - orifile=`cat roilstorifile.txt`; \ | |
433 | - oribase=`printf "%s%03d" $$orifile $(ZEROTILTNUM)`; \ | |
434 | - num=`printf "%04d" $$i`; \ | |
435 | - roiname=`printf "%s-%04d" $$oribase $$i`; \ | |
436 | - echo "$$roiname in $$oribase.roiinfo"; \ | |
437 | - if [ $$orifile = a ]; then \ | |
438 | - echo "$*.roilst creating" ; \ | |
439 | - ls -1 $$orifile*-$$num.mul | awk -v eularangle=$(EULARANGLE) ' BEGIN {angle=$(MINTILTANGLE)} {printf("%s %s %f %f %f\n",$$1, eularangle, 0, angle, 0); angle+=$(STEPTILTANGLE); if (angle == 0 ) { angle+=$(STEPTILTANGLE)} }' > $*.roilst ; \ | |
440 | - else \ | |
441 | - echo "$*.roilst creating" ; \ | |
442 | - ls -1 $$orifile*-$$num.mul | awk -v eularangle=$(EULARANGLE) ' BEGIN {angle=$(MINTILTANGLE)} {printf("%s %s %f %f %f\n",$$1, eularangle, angle, 0, 0); angle+=$(STEPTILTANGLE); if (angle == 0 ){ angle+=$(STEPTILTANGLE)} }' > $*.roilst ; \ | |
443 | - fi | |
444 | -##sep | |
445 | - | |
446 | -ds6Create:: | |
447 | - roiNum=`ls -1 *-d.roilst | wc -l `; \ | |
448 | - echo $$roiNum; \ | |
449 | - for (( i=0 ; $$i < $$roiNum ; i=$$(($$i + 1)) )) \ | |
450 | - do \ | |
451 | - tiltnum=`printf "%04d-d" $$i`; \ | |
452 | - make $$tiltnum.mrc3d ; \ | |
453 | - make $$tiltnum.mrc3dcentre ; \ | |
454 | - make $$tiltnum.mrc3dcentremed ; \ | |
455 | - make $$tiltnum.mrc3dcentremedpad ; \ | |
456 | - make $$tiltnum.ds6 ; \ | |
457 | - done ; | |
458 | - for (( j=1 ; $$j <= $(TILTAXISNUMBER) ; j=$$(($$j + 1 )))) \ | |
459 | - do \ | |
460 | - if [ $$j -eq 1 ]; then \ | |
461 | - orifile=$(ORIFILE1); \ | |
462 | - else \ | |
463 | - orifile=$(ORIFILE2); \ | |
464 | - fi ; \ | |
465 | - roiNum=`ls -1 $$orifile*.roilst | wc -l `; \ | |
466 | - echo $$roiNum; \ | |
467 | - for (( i=0 ; $$i < $$roiNum ; i=$$(($$i + 1)) )) \ | |
468 | - do \ | |
469 | - tiltnum=`printf "%s%04d" $$orifile $$i`; \ | |
470 | - make $$tiltnum.mrc3d ; \ | |
471 | - make $$tiltnum.mrc3dcentre ; \ | |
472 | - make $$tiltnum.mrc3dcentremed ; \ | |
473 | - make $$tiltnum.mrc3dcentremedpad ; \ | |
474 | - make $$tiltnum.ds6 ; \ | |
475 | - done ; \ | |
476 | - done | |
477 | - | |
478 | -Cleands6:: | |
479 | - rm *.mrc3d | |
480 | - | |
481 | -.roilst.mrc3d: | |
482 | - mrc2Dto3D -I $*.roilst -o $*.mrc3d -single 0 -InterpolationMode 2 -m 1 | |
483 | - | |
484 | -.mrc3dcentremedpad.ds6: | |
485 | - mrc2map -i $*.mrc3d -o $*.ds6 -mul $(MRC2MAPmul) | |
486 | - mrc2map -i $*.mrc3dcentre -o $*.centre.ds6 -mul $(MRC2MAPmul) | |
487 | - mrc2map -i $*.mrc3dcentremed -o $*.centremed.ds6 -mul $(MRC2MAPmul) | |
488 | - mrc2map -i $*.mrc3dcentremedpad -o $*.centremedpad.ds6 -mul $(MRC2MAPmul) | |
489 | - | |
490 | -.mrc3dcentre.mrc3dcentremed: | |
491 | - mrcImageSmoothing -i $*.mrc3dcentre -o $*.mrc3dcentremed -r $(SMOOTHr) -m $(SMOOTHm) -M $(SMOOTHM) | |
492 | - | |
493 | -.mrc3d.mrc3dcentre: | |
494 | - mrcImageCenterGet -i $*.mrc3d -o $*.mrc3dcentre -Nx $(CENTERNX) -Ny $(CENTERNY) -Nz $(CENTERNZ) | |
495 | - | |
496 | -.mrc3dcentremed.mrc3dcentremedpad: | |
497 | - mrcImage3DPad -i $*.mrc3dcentremed -o $*.mrc3dcentremedpad -r 90 -m 1 -M 1 -v 0 | |
498 | - |
@@ -1,377 +0,0 @@ | ||
1 | -.SUFFIXES: .tif .dat .prectf .fft .ctfinfo .mrc \ | |
2 | - .box .roi \ | |
3 | - .roiinfo-coarse .roiall-coarse .roiinforef-coarse .roiallref-coarse .roiinfoprev-coarse \ | |
4 | - .roiinfo .roiall .roiinforef .roiallref .roiinfoprev \ | |
5 | - .box2roiinfo .roiinfo2box .corinfo \ | |
6 | - .roilst0 .roilst .mrc3d .ds6 .mrc3dmed .mrc3dmedcentre .orimrccentre .oricorinfo .coarroi .coarroiinfo .roi2 .pad .mul | |
7 | - | |
8 | -RESOLUTION=2.0 | |
9 | -ORIFILE1=a | |
10 | -ORIFILE2=b | |
11 | -EULARANGLE1=XOYS | |
12 | -EULARANGLE2=YOYS | |
13 | -ZEROTILTNUM=60 | |
14 | -ORIBASE=$(ORIFILE1)0$(ZEROTILTNUM) | |
15 | -MINTILTNUM=0 | |
16 | -MAXTILTNUM=121 | |
17 | -STEPTILTNUM=4 | |
18 | -MINTILTANGLE=-60 | |
19 | -STEPTILTANGLE=4 | |
20 | -TILTAXISNUMBER=2 | |
21 | - | |
22 | -SMOOTHr=5 | |
23 | -SMOOTHm=1 | |
24 | -SMOOTHM=1 | |
25 | -CENTERNX=100 | |
26 | -CENTERNY=100 | |
27 | -CENTERNZ=100 | |
28 | - | |
29 | -ORINX=1024 | |
30 | -ORINY=1024 | |
31 | - | |
32 | -ORINZ=1 | |
33 | -ORISIZEX=2048 | |
34 | -ORISIZEY=2048 | |
35 | -TILTAXISANGLE=0 | |
36 | -IMAGEPADm=15 | |
37 | -IMAGEPADV=0 | |
38 | -IMAGEMULTIPLYINGv=-1 | |
39 | -MRC2MAPmul=1 | |
40 | - | |
41 | -all: data2mrc COARCreate $(ORIBASE).roiinfo BoxCreate PadAndMul roilstCreate ds6Create | |
42 | - | |
43 | -data2mrc:: | |
44 | - for (( j=1 ; $$j <= $(TILTAXISNUMBER) ; j=$$(($$j + 1 )))) \ | |
45 | - do \ | |
46 | - if [ $$j -eq 1 ]; then \ | |
47 | - orifile=$(ORIFILE1); \ | |
48 | - else \ | |
49 | - orifile=$(ORIFILE2); \ | |
50 | - fi ; \ | |
51 | - for (( i=0 ; $$i <= $(MAXTILTNUM) ; i=$$(($$i + $(STEPTILTNUM))))) \ | |
52 | - do \ | |
53 | - ori=`printf "%s%03d" $$orifile $$i`; \ | |
54 | - make $$ori.prectf; \ | |
55 | - make $$ori.fft; \ | |
56 | - make $$ori.ctfinfo; \ | |
57 | - make $$ori.mrc ; \ | |
58 | - done ; \ | |
59 | - done | |
60 | - | |
61 | -.dat.prectf: | |
62 | - emData2mrc -i $*.dat -o $*.prectf | |
63 | - | |
64 | -.tif.prectf: | |
65 | - tiff2mrc -i $*.tif -o $*.prectf -m 1 -r $(RESOLUTION) | |
66 | - | |
67 | -.prectf.fft: | |
68 | - mrcImageFFT -i $*.prectf -o $*.fft | |
69 | - | |
70 | -.fft.ctfinfo: | |
71 | - ctfDetermine -i $*.fft -o $*.ctfinfoauto -D 50000 -m 163 -CutLow 0.07 -CutHigh 0.15 -Ac 200 -d 3 -Cc 0.7 -Cs 2.1 | |
72 | - awk '/^[^:]/ {print $$0};' $*.ctfinfoauto > $*.ctfinfo | |
73 | - Display2 -i $*.fft -H 1e7 | |
74 | - | |
75 | -.ctfinfo.mrc: | |
76 | - mv $*.ctfinfo $*.ctfinfotmp | |
77 | - awk '/^[^:]/ {print $$0};' $*.ctfinfotmp > $*.ctfinfo | |
78 | - mrcImageCTFCompensation -i $*.fft -o $*.mrc -info2 $*.ctfinfo -m 1 | |
79 | - Display2 -i $*.mrc -geometry 1024x1024+0+0 -Zoom 0.5 & | |
80 | - | |
81 | -##sep | |
82 | - | |
83 | -COARCreate:: | |
84 | - for (( j=1 ; $$j <= $(TILTAXISNUMBER) ; j=$$(($$j + 1 )))) \ | |
85 | - do \ | |
86 | - if [ $$j -eq 1 ]; then \ | |
87 | - orifile=$(ORIFILE1); \ | |
88 | - oribase=$(ORIFILE1)0$(ZEROTILTNUM); \ | |
89 | - else \ | |
90 | - orifile=$(ORIFILE2); \ | |
91 | - oribase=$(ORIFILE2)0$(ZEROTILTNUM); \ | |
92 | - fi ; \ | |
93 | - ORIZEROFILE=$$oribase.mrc; \ | |
94 | - previous=$$oribase ; \ | |
95 | - echo "Create the zero-tilt Image ROI (same as mrc, just copy)"; \ | |
96 | - make $$oribase.roiinfo-coarse; echo "$$oribase.roi will be created."; \ | |
97 | - make $$oribase.roiall-coarse ; echo "$$oribase.roi was created." ; \ | |
98 | - make $$oribase.orimrccentre ; echo "The centre of $$oribase.roi was created." ; \ | |
99 | - for delta in -$(STEPTILTNUM) +$(STEPTILTNUM) ; do \ | |
100 | - echo "$$delta"; \ | |
101 | - previous=$$oribase; \ | |
102 | - TOTALX=0 ; \ | |
103 | - TOTALY=0 ; \ | |
104 | - echo "Total: $$TOTALX $$TOTALY" > $$previous.oricorinfo; \ | |
105 | - for (( i=$$(($(ZEROTILTNUM) + $$delta )) ; (( $(MINTILTNUM) <= $$i && $$i <= $(MAXTILTNUM) )) ; i=$$(($$i + $$delta)) )) \ | |
106 | - do \ | |
107 | - ori=`printf "%s%03d" $$orifile $$i`; \ | |
108 | - echo $$i, $(MAXTILTNUM), $(MINTILTNUM); \ | |
109 | - echo $$previous > $$ori.prevname ; \ | |
110 | - ln -sf $$previous.roiinfo-coarse $$ori.roiinfoprev-coarse ; make $$ori.roiinforef-coarse ; \ | |
111 | - make $$ori.roiallref-coarse ; \ | |
112 | - make $$ori.orimrccentre ; \ | |
113 | - TOTALX=`awk '/Total:/ { print $$2 }' $$previous.oricorinfo`; \ | |
114 | - TOTALY=`awk '/Total:/ { print $$3 }' $$previous.oricorinfo`; \ | |
115 | - ln -sf $$previous.orimrccentre $$ori.oriref ; make $$ori.oricorinfo ; \ | |
116 | - echo "TotalPrev: $$TOTALX $$TOTALY" >> $$ori.oricorinfo; \ | |
117 | - make $$ori.roiinfo-coarse; \ | |
118 | - make $$ori.roiall-coarse; \ | |
119 | - touch $$ori.orimrccentre ; sleep 1; \ | |
120 | - touch $$ori.oricorinfo ; sleep 1; \ | |
121 | - touch $$ori.roiinfo-coarse; sleep 1; \ | |
122 | - touch $$ori.roiall-coarse; sleep 1; \ | |
123 | - previous=$$ori; \ | |
124 | - done ; \ | |
125 | - echo $$i; \ | |
126 | - done ; \ | |
127 | - done | |
128 | - | |
129 | -CleanCOARCreate:: | |
130 | - rm *.roiinfo-coarse | |
131 | - | |
132 | -$(ORIFILE1)0$(ZEROTILTNUM).roiinfo-coarse: | |
133 | - printf "%s%03d.roi Rect2 %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f\n" $(ORIFILE1) $(ZEROTILTNUM) $$(($(ORISIZEX)/2)) $$(($(ORISIZEY)/2)) $(ORISIZEX) $(ORISIZEY) > $(ORIFILE1)0$(ZEROTILTNUM).roiinfo-coarse | |
134 | - | |
135 | -$(ORIFILE2)0$(ZEROTILTNUM).roiinfo-coarse: | |
136 | - printf "%s%03d.roi Rect2 %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f\n" $(ORIFILE2) $(ZEROTILTNUM) $$(($(ORISIZEX)/2)) $$(($(ORISIZEY)/2)) $(ORISIZEX) $(ORISIZEY) > $(ORIFILE2)0$(ZEROTILTNUM).roiinfo-coarse | |
137 | - | |
138 | -.orimrccentre.oricorinfo: | |
139 | - echo "orimrccentre -> oricorinfo: SHIFT:= current.orimrccentre for prev.oirmrccentre" | |
140 | - mrcImageCorrelation -i $*.orimrccentre -r $*.oriref -c $*.oricorinfo | |
141 | - | |
142 | -.oricorinfo.roiinfo-coarse: | |
143 | - TOTALX=`awk '/TotalPrev:/ { print $$2 }' $*.oricorinfo`; \ | |
144 | - TOTALY=`awk '/TotalPrev:/ { print $$3 }' $*.oricorinfo`; \ | |
145 | - SHIFTX=`awk '/Shift:/ { print $$2 }' $*.oricorinfo`; \ | |
146 | - SHIFTY=`awk '/Shift:/ { print $$3 }' $*.oricorinfo`; \ | |
147 | - echo "$$SHIFTX,$$SHIFTY"; \ | |
148 | - if [ $$(($(ORINX)/2)) -lt $$SHIFTX ]; then \ | |
149 | - SHIFTX=$$(($$SHIFTX - $(ORINX))) ; \ | |
150 | - fi ; \ | |
151 | - if [ $$(($(ORINY)/2)) -lt $$SHIFTY ]; then \ | |
152 | - SHIFTY=$$(($$SHIFTY - $(ORINY))) ; \ | |
153 | - fi ; \ | |
154 | - TOTALX=$$(($$TOTALX + $$SHIFTX)) ; \ | |
155 | - TOTALY=$$(($$TOTALY + $$SHIFTY)) ; \ | |
156 | - echo "Total: $$TOTALX $$TOTALY" >> $*.oricorinfo; \ | |
157 | - printf "%s.roi Rect2 %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f\n" $* $$(($(ORISIZEX)/2+$$TOTALX)) $$(($(ORISIZEY)/2+$$TOTALY)) $(ORISIZEX) $(ORISIZEY) $$TOTALX $$TOTALY > $*.roiinfo-coarse ; | |
158 | - | |
159 | -.roiinfo-coarse.roiall-coarse: | |
160 | - make $*.mrc | |
161 | - mrcImageROIs -i $*.mrc -I $*.roiinfo-coarse | |
162 | - touch $*.roiall-coarse | |
163 | - | |
164 | -.roiinforef-coarse.roiallref-coarse: | |
165 | - mrcImageROIs -i $*.mrc -I $*.roiinforef-coarse | |
166 | - touch $*.roiallref-coarse | |
167 | - | |
168 | -.roiinfoprev-coarse.roiinforef-coarse: | |
169 | - previous=`cat $*.prevname`; \ | |
170 | - sed -s s/$$previous/$*/ $*.roiinfoprev-coarse > $*.roiinforef-coarse | |
171 | - | |
172 | -.roi.orimrccentre: | |
173 | - echo "roi -> orimrccentre" | |
174 | - mrcImageCenterGet -i $*.roi -o $*.orimrccentre -Nx $(ORINX) -Ny $(ORINY) -Nz $(ORINZ) | |
175 | - | |
176 | - | |
177 | -$(ORIBASE).roiinfo: | |
178 | - Display2 -i $*.mrc ; \ | |
179 | - if [ $(TILTAXISNUMBER) -eq 2 ]; then \ | |
180 | - Display2 -i $(ORIFILE2)0$(ZEROTILTNUM).mrc ;\ | |
181 | - fi; \ | |
182 | - | |
183 | -##sep | |
184 | - | |
185 | -BoxCreate:: | |
186 | - for (( j=1 ; $$j <= $(TILTAXISNUMBER) ; j=$$(($$j + 1 )))) \ | |
187 | - do \ | |
188 | - if [ $$j -eq 1 ]; then \ | |
189 | - orifile=$(ORIFILE1); \ | |
190 | - oribase=$(ORIFILE1)0$(ZEROTILTNUM); \ | |
191 | - else \ | |
192 | - orifile=$(ORIFILE2); \ | |
193 | - oribase=$(ORIFILE2)0$(ZEROTILTNUM); \ | |
194 | - fi ; \ | |
195 | - BASE=$$oribase; \ | |
196 | - make $$BASE.roiall ; \ | |
197 | - for delta in +$(STEPTILTNUM) -$(STEPTILTNUM) ; \ | |
198 | - do \ | |
199 | - echo $$delta,$$BASE ; \ | |
200 | - previous=$$BASE ; \ | |
201 | - for (( i=$$(($(ZEROTILTNUM)+$$delta)) ; (( $(MINTILTNUM) <= $$i && $$i <= $(MAXTILTNUM) )) ; i=$$(($$i + $$delta)) )) \ | |
202 | - do \ | |
203 | - name=`printf "%s%03d" $$orifile $$i`; \ | |
204 | - echo ">>>> $$name $$i start: $(ZEROTILTNUM) $$(($(ZEROTILTNUM)+$$delta))"; \ | |
205 | - echo $$previous > $$name.prevname; \ | |
206 | - ln -sf $$previous.roiinfo $$name.roiinfoprev ; make $$name.roiinforef ; \ | |
207 | - make $$name.roiinfo ; \ | |
208 | - make $$name.roiall ; \ | |
209 | - previous=$$name ; \ | |
210 | - done ; \ | |
211 | - done ; \ | |
212 | - done | |
213 | - | |
214 | -CleanBoxCreate:: | |
215 | - rm *.roiall *.roiinfoprev *.roiinforef | |
216 | - | |
217 | -.roiinfo.roiall: | |
218 | - make $*.roi | |
219 | - mrcImageROIs -i $*.roi -I $*.roiinfo | |
220 | - touch $*.roiall | |
221 | - | |
222 | -.roiinforef.roiallref: | |
223 | - mrcImageROIs -i $*.roi -I $*.roiinforef | |
224 | - touch $*.roiallref | |
225 | - | |
226 | -.roi.corinfo: | |
227 | - mrcImageCorrelation -i $*.roi -r $*.ref -c $*.corinfo | |
228 | - | |
229 | -.roiinforef.roiinfo: | |
230 | - # Extract ROIs using roiinfo of previous image (roiinforef) | |
231 | - previous=`cat $*.prevname`; \ | |
232 | - echo $$previous; \ | |
233 | - echo "start initial ROI"; \ | |
234 | - make $*.roiallref ; \ | |
235 | - ROINUM=`wc -l $*.roiinforef | awk '{print $$1}'` ; \ | |
236 | - echo $$ROINUM ; \ | |
237 | - rm -f $*.roiinfo ; \ | |
238 | - for (( j=0 ; $$j < $$ROINUM ; j=$$(($$j+1)))) \ | |
239 | - do \ | |
240 | - input=`printf "%s-%04d" $* $$j`; \ | |
241 | - refer=`printf "%s-%04d" $$previous $$j`; \ | |
242 | - echo $$input, $$refer ; \ | |
243 | - ln -sf $$refer.roi $$input.ref ; make $$input.corinfo ; \ | |
244 | - SHIFTX=`awk '/Shift:/ { print $$2 }' $$input.corinfo`; \ | |
245 | - SHIFTY=`awk '/Shift:/ { print $$3 }' $$input.corinfo`; \ | |
246 | - echo "$$SHIFTX,$$SHIFTY,$$input"; \ | |
247 | - awk -v SHIFTXAWK0=$$SHIFTX -v SHIFTYAWK0=$$SHIFTY -v NAMEAWK=$$input -v BASEAWK=$* \ | |
248 | - ' /'"$$input"'/ { \ | |
249 | - SizeX = int(sqrt(($$5-$$3)*($$5-$$3)+( $$6-$$4)*( $$6-$$4))+0.5)+1; \ | |
250 | - SizeY = int(sqrt(($$9-$$3)*($$9-$$3)+($$10-$$4)*($$10-$$4))+0.5)+1; \ | |
251 | - if(SizeX/2<SHIFTXAWK0) { \ | |
252 | - SHIFTXAWK0 -= SizeX ; \ | |
253 | - } \ | |
254 | - if(SizeY/2<SHIFTYAWK0) { \ | |
255 | - SHIFTYAWK0 -= SizeY ; \ | |
256 | - } \ | |
257 | - Angle = atan2($$6-$$4, $$5-$$3) ; \ | |
258 | - SHIFTXAWK = +1*( SHIFTXAWK0*cos(Angle) - SHIFTYAWK0*sin(Angle)); \ | |
259 | - SHIFTYAWK = +1*( SHIFTXAWK0*sin(Angle) + SHIFTYAWK0*cos(Angle)); \ | |
260 | - printf("%s-%04d.roi Rect %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f %15.6f\n", BASEAWK, FNR-1, $$3+SHIFTXAWK, $$4+SHIFTYAWK, $$5+SHIFTXAWK, $$6+SHIFTYAWK, $$7+SHIFTXAWK, $$8+SHIFTYAWK, $$9+SHIFTXAWK, $$10+SHIFTYAWK, SHIFTXAWK, SHIFTYAWK, Angle*180/3.141592653589793238, SizeX, SizeY);}' \ | |
261 | - $*.roiinforef | tee -a $*.roiinfo ; \ | |
262 | - done ; | |
263 | - | |
264 | -.roiinfoprev.roiinforef: | |
265 | - previous=`cat $*.prevname`; \ | |
266 | - sed -s s/$$previous/$*/ $*.roiinfoprev > $*.roiinforef | |
267 | - | |
268 | -PadAndMul:: | |
269 | - for (( j=1 ; $$j <= $(TILTAXISNUMBER) ; j=$$(($$j + 1 )))) \ | |
270 | - do \ | |
271 | - if [ $$j -eq 1 ]; then \ | |
272 | - orifile=$(ORIFILE1); \ | |
273 | - oribase=$(ORIFILE1)0$(ZEROTILTNUM); \ | |
274 | - else \ | |
275 | - orifile=$(ORIFILE2); \ | |
276 | - oribase=$(ORIFILE2)0$(ZEROTILTNUM); \ | |
277 | - fi ; \ | |
278 | - roiNum=`wc -l $$oribase.roiinfo | awk '{print $$1}' `; \ | |
279 | - for ((i=$(MINTILTNUM);i<=$(MAXTILTNUM);i=$$(($$i+$(STEPTILTNUM))) )) \ | |
280 | - do \ | |
281 | - for((j=0;j<$$roiNum;j=$$(($$j+1)))) \ | |
282 | - do \ | |
283 | - name=`printf "%s%03d-%04d" $$orifile $$i $$j` ; \ | |
284 | - echo $$(($$i-$(ZEROTILTNUM))) > $$name.tilt ; \ | |
285 | - ln -sf $$name.roi $$name.roi2 ; make $$name.pad ; \ | |
286 | - make $$name.mul ; \ | |
287 | - done ; \ | |
288 | - done ; \ | |
289 | - done | |
290 | - | |
291 | -CleanPadAndMul:: | |
292 | - rm *.pad | |
293 | - | |
294 | -.roi2.pad: | |
295 | - tilt=`cat $*.tilt` ; \ | |
296 | - echo $$tilt ; \ | |
297 | - mrcImagePad -i $*.roi2 -o $*.pad -Tilt $(TILTAXISANGLE) $$tilt -m $(IMAGEPADm) -V $(IMAGEPADV) | |
298 | - | |
299 | -.pad.mul: | |
300 | - mrcImageMultiplying -i $*.pad -o $*.mul -v $(IMAGEMULTIPLYINGv) | |
301 | - | |
302 | -roilstCreate:: | |
303 | - for (( j=1 ; $$j <= $(TILTAXISNUMBER) ; j=$$(($$j + 1 )))) \ | |
304 | - do \ | |
305 | - if [ $$j -eq 1 ]; then \ | |
306 | - orifile=$(ORIFILE1); \ | |
307 | - oribase=$(ORIFILE1)0$(ZEROTILTNUM); \ | |
308 | - echo "$(ORIFILE1)" > roilstorifile.txt; \ | |
309 | - else \ | |
310 | - orifile=$(ORIFILE2); \ | |
311 | - oribase=$(ORIFILE2)0$(ZEROTILTNUM); \ | |
312 | - echo "$(ORIFILE2)" > roilstorifile.txt; \ | |
313 | - fi ; \ | |
314 | - roiNum=`wc -l $$oribase.roiinfo | awk '{print $$1}' `; \ | |
315 | - echo $$roiNum; \ | |
316 | - for (( i=0 ; $$i < $$roiNum ; i=$$(($$i + 1)))) \ | |
317 | - do \ | |
318 | - tiltnum=`printf "%s%04d" $$orifile $$i`; \ | |
319 | - echo "$$i" > $$tiltnum.roilst0; \ | |
320 | - make $$tiltnum.roilst ; \ | |
321 | - done ; \ | |
322 | - done | |
323 | - | |
324 | -Cleanroilst:: | |
325 | - rm *.roilst | |
326 | - | |
327 | -.roilst0.roilst: | |
328 | - i=`cat $*.roilst0`; \ | |
329 | - orifile=`cat roilstorifile.txt`; \ | |
330 | - oribase=`printf "%s%03d" $$orifile $(ZEROTILTNUM)`; \ | |
331 | - num=`printf "%04d" $$i`; \ | |
332 | - roiname=`printf "%s-%04d" $$oribase $$i`; \ | |
333 | - echo "$$roiname in $$oribase.roiinfo"; \ | |
334 | - tiltangle=`awk -v ORIBASE=$$oribase '/'"$$roiname"'/ { printf("%15.6f", atan2($$6-$$4, $$5-$$3)*180/3.141592653589793238); }' ORIBASE.roiinfo`; \ | |
335 | - if [ $orifile = a ]; then \ | |
336 | - eularangle=$(EULARANGLE1); \ | |
337 | - else \ | |
338 | - eularangle=$(EULARANGLE2); \ | |
339 | - fi ; \ | |
340 | - echo "$*.roilst creating" ; \ | |
341 | - ls -1 *-$$num.mul | awk ' BEGIN {angle=$(MINTILTANGLE)} {printf("%s %s %f %f %f\n",$$1, '"$$eularangle"', -1*angle, -1*'"$$tiltangle"', 0); angle+=$(STEPTILTANGLE)}' > $*.roilst ; \ | |
342 | - | |
343 | -##sep | |
344 | - | |
345 | -##noshow | |
346 | -ds6Create:: | |
347 | - roiNum=`ls -1 *.roilst | wc -l `; \ | |
348 | - echo $$roiNum; \ | |
349 | - for (( i=0 ; $$i < $$roiNum ; i=$$(($$i + 1)) )) \ | |
350 | - do \ | |
351 | - tiltnum=`printf "%s-%04d" $(ORIFILE) $$i`; \ | |
352 | - make $$tiltnum.mrc3d ; \ | |
353 | - make $$tiltnum.mrc3dmed ; \ | |
354 | - make $$tiltnum.mrc3dmedcentre ; \ | |
355 | - make $$tiltnum.ds6 ; \ | |
356 | - done ; | |
357 | - | |
358 | -##noshow | |
359 | - | |
360 | -Cleands6:: | |
361 | - rm *.mrc3d | |
362 | - | |
363 | -.roilst.mrc3d: | |
364 | - mrc2Dto3D -I $*.roilst -o $*.mrc3d -single 0 -InterpolationMode 2 -m 1 | |
365 | - | |
366 | -.mrc3dmedcentre.ds6: | |
367 | - mrc2map -i $*.mrc3d -o $*.ds6 -mul $(MRC2MAPmul) | |
368 | - mrc2map -i $*.mrc3dmed -o $*.med.ds6 -mul $(MRC2MAPmul) | |
369 | - mrc2map -i $*.mrc3dmedcentre -o $*.medcentre.ds6 -mul $(MRC2MAPmul) | |
370 | - | |
371 | -.mrc3d.mrc3dmed: | |
372 | - mrcImageSmoothing -i $*.mrc3d -o $*.mrc3dmed -r $(SMOOTHr) -m $(SMOOTHm) -M $(SMOOTHM) | |
373 | - | |
374 | -.mrc3dmed.mrc3dmedcentre: | |
375 | - mrcImageCenterGet -i $*.mrc3dmed -o $*.mrc3dmedcentre -Nx $(CENTERNX) -Ny $(CENTERNY) -Nz $(CENTERNZ) | |
376 | - | |
377 | - |
@@ -1 +0,0 @@ | ||
1 | -test |
@@ -1 +0,0 @@ | ||
1 | -test |
@@ -1,6 +0,0 @@ | ||
1 | -data2mrc | |
2 | -COARCreate | |
3 | -roilstCreate | |
4 | -mrcImagePadAndMultiplying | |
5 | -ds6Create | |
6 | - |
@@ -347,10 +347,11 @@ then | ||
347 | 347 | source $ATISTREAMROOT/.bash_profile |
348 | 348 | fi |
349 | 349 | else |
350 | - OIFS=$IFS | |
351 | - IFS=":" | |
352 | - read k v < <(echo "$LD_LIBRARY_PATH" | egrep "$ATISTREAMSDKROOT" ) | |
353 | - IFS=$OIFS | |
350 | +# OIFS=$IFS | |
351 | +# IFS=":" | |
352 | +# read k v < <(echo "$LD_LIBRARY_PATH" | egrep "$ATISTREAMSDKROOT" ) | |
353 | + k=`echo "$LD_LIBRARY_PATH" | egrep "$ATISTREAMSDKROOT"` | |
354 | +# IFS=$OIFS | |
354 | 355 | # echo "()k=$k" |
355 | 356 | # echo "v=$v" |
356 | 357 | if [ -z "$k" ]; |
@@ -385,10 +386,11 @@ then | ||
385 | 386 | done |
386 | 387 | fi |
387 | 388 | else |
388 | - OIFS=$IFS | |
389 | - IFS=":" | |
390 | - read k v < <(echo "$LD_LIBRARY_PATH" | egrep "$LD_LIBGFORTRAN" ) | |
391 | - IFS=$OIFS | |
389 | +# OIFS=$IFS | |
390 | +# IFS=":" | |
391 | +# read k v < <(echo "$LD_LIBRARY_PATH" | egrep "$LD_LIBGFORTRAN" ) | |
392 | + k=`echo "$LD_LIBRARY_PATH" | egrep "$LD_LIBGFORTRAN` | |
393 | +# IFS=$OIFS | |
392 | 394 | # echo "()k=$k" |
393 | 395 | # echo "v=$v" |
394 | 396 | if [ -z "$k" ]; |
@@ -16,7 +16,9 @@ if [ ! -d $EOS_HOME0 ] ; then | ||
16 | 16 | mkdir $EOS_HOME0 |
17 | 17 | fi |
18 | 18 | |
19 | -curl -L $ROOT_EOS_SOURCEFORGE/env/Eos_env?export=raw | sh; git clone --depth 1 git://git.sourceforge.jp/gitroot/eos/base.git $EOS_HOME0 | |
19 | +#curl -L $ROOT_EOS_SOURCEFORGE/env/Eos_env?export=raw | sh - ; | |
20 | + | |
21 | +git clone --depth 1 git://git.sourceforge.jp/gitroot/eos/base.git $EOS_HOME0 | |
20 | 22 | |
21 | 23 | if [ -z $EOS_HOME0 ] ; then |
22 | 24 | export EOS_HOME0=$HOME/Eos |