Thesis-PDF - IAP/TU Wien

function join_amptampr (w) // joining amplitude trace and retrace
Wave w
Variable laynr
Variable i
Variable j
Variable offset=0
Variable used_offset
Variable test
Variable test1
Variable separationline=0.00000005
// used to distinguish from which area information comes
Variable sep=0 // if 1: active
Variable fuzz=0.2 // 0: offset to be used: mean-offset,
// 1: dot-per-dot-offset, or use value in between
print ("separation line height if present is:"), separationline
print ("fuzz parameter is:"), fuzz
laynr = DimSize(w, 2) // number of layers
Redimension/N=(-1,-1,laynr+1) w // create a new layer
for(i=0;i< DimSize(w, 0);i+=1) // initialize variables;continue test
for(j=0;j< DimSize(w, 1);j+=1)
offset = offset+w[i][j][2]+w[i][j][3]
endfor // condition;update loop variables
endfor
//calculating mean offset for amplitude traces
offset=offset/(2*DimSize(w, 0)*DimSize(w,1))
print ("mean offset is:"), offset
for(i=0;i< DimSize(w, 0);i+=1)
for(j=0;j< DimSize(w, 1);j+=1)
// initialize variables;continue test
used_offset = offset*(1-fuzz) +fuzz*(w[i][j][2]+w[i][j][3])/2
if (w[i][j][0]==min(w[i][j][0],w[i][j][1]))
w[i][j][laynr]=w[i][j][2]-used_offset
//taken from amplTRACE if corresponding heightTRACE is min
test=0
else
106