equalize_limits.pro

IDL语言编写的用于天文自适应光学仿真的软件CAOS V6.0的第一部分。

; $Id: equalize_limits.pro,v 1.2 2002/03/14 11:49:12 riccardi Exp $

function equalize_limits, ima, LOW_THRESHOLD=low_t, HIGH_THRESHOLD=high_t, BINSIZE=binsize


if n_elements(low_t) eq 0 then low_t=0.01
if n_elements(high_t) eq 0 then high_t=0.99

n=1000
minv = min(ima, MAX=maxv)

if minv ne maxv then begin

	if n_elements(binsize) eq 0 then binsize = (maxv-minv)/float(n)

	p = double(histogram(ima, OMAX=maxv, OMIN=minv, BINSIZE=binsize))
	for i=1,n-1 do p[i]=p[i-1]+p[i]
	p = p/p[n-1]
	idx = where((p ge low_t) and (p le high_t), count)

	low_limit = minv+binsize*min(idx)
	high_limit = minv+binsize*(max(idx)+1)
endif else begin
	low_limit = minv
	high_limit = maxv
endelse

return, [low_limit, high_limit]

end