orthogonal array projections [message #7545] |
Wed, 20 November 1996 00:00 |
deb
Messages: 12 Registered: April 1996
|
Junior Member |
|
|
I'm trying to make a 2-d projection of max amplitude in a 3-d image
array. The following bit of Fortran-esque code does what i want,
but it's excruciatingly slow:
for k=0,n-1 do begin
for l=0,z-1 do begin
ayz(k,l)=max(image(low:hi,k,l)
endfor
endfor
Is there a way to do this without using loops? If i simply say
ayz(*,*)=max(image(low:hi,*,*))
i end up with the max value over the whole array and not the max
value in each depth bin in the array. (which makes sense looking
at the code, but i can't figure out how to isolate each depth bin
w/o using the loop business). Voxel_proj and Project_vol are
overkill for this application (and even slower than the loop stuff
anyway).
Any thoughts would be greatly appreciated:)
-Deb Summa
summa@lanl.gov
|
|
|