[x31,x39,x36,x56,x59]
fr
?micrograph image?
;<1> micx11
                                          
fr
?pw_avg?
;<2> power/pw_avgx11
                                          
fr
?pw_ro?
;<3>  power/pw_rox11

fr
?roo?
;<4>  power/roox11

fi x23,x24
<1>
12,2
;x24           Nrow
;x23           Nsam
;x87=((x23-2*x59)*(x24-2*x56))number of pieces in this micrograph

x36=100/(100-x36); normalization of % of the overlap in y 
x39=100/(100-x39); normalization of % of the overlap in x

x26=INT(x36*((x24-2*x56)/x31-1)+1)     ;number of pieces vertical dim.
x29=INT(x39*((x23-2*x59)/x31-1)+1)     ;number of pieces horizontal dim.
x93 = 0
x94 = 0
do lb2 j = 1, x26
x12 = x0
x21 = (x31/x36)*(x12-1) + x56   ;y-direction
do lb3 k = 1, x29
x13 = x0
x22 = (x31/x39)*(x13-1) + x59  ; x-direction
x93 = x93 + 1
wi
<1>       ;get small pieces of the input image 
_5
x31,x31           ;size of small pieces
x22,x21           ;coordinates of the upper left corner
ra                ;correct ramp effects 
_5
_6
pw                ;calc power spectrum of each small piece
_6
_7
sq               
_7  
_9
x94=x94+1
if (x94.ge.2) go to lb5    
cp
_9
_3
goto lb3
lb5
ad
_3
_9
_3
*
lb3
lb2
;
; average over power spectra of small pieces 
;
AR
_3
_8
((P1/x94)+0.) 
;power/pw_varx11
wu                ;calculate the square root 
_8
<2>
ro
<2>
<3>
li d
<3>
<4>
r
1
re