function norm3d,img
image=reform(img)
n=size(image)
if n(0) eq 3 then begin
out=float(image)
for i=0,n(3)-1 do begin
 mx=max(float(out(*,*,i)),min=mn)
 if mx eq 0 then mx=1.
 if mn eq 0 then mn=1.
 out(*,*,i)=(out(*,*,i)>0)/mx+(out(*,*,i)<0)/abs(mn)
end
return,out
endif else begin
 mx=max(float(image),min=mn)
 if mx eq 0 then mx=1.
 if mn eq 0 then mn=1.
 return, (image>0)/mx+(image<0)/abs(mn)
end
end