Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.eso.org/~chummel/oyster/lab/params.fringevissq.html Дата изменения: Mon Nov 26 08:49:04 2007 Дата индексирования: Fri Feb 28 14:18:35 2014 Кодировка: Поисковые слова: ngc 5128

function fringevissq,c,k,n,normal=normal
;
; Return the squared normalized visibility amplitude. Zero if no counts.
; The first index in the frames must correspond to the phase bins.
; Use k to select a baseline, if multiple baseline are present.
;
common Constants,c_light,pi_circle,e_euler,i_complex,a_disp,b_disp
;
if n_elements(k) eq 0 then k=1
if n_elements(n) eq 0 then n=1
;
if n_elements(normal) eq 0 then normal=0 else normal=normal gt 0
;
c=reform(c)
s=fringenphot(c)
r=size(c)
n0=r(1)
nominator=add(abs(fringevis(c-avg(c),k))^2-s,n,r(0)-2)/n
denominator=(add(s,n,r(0)-2)/n)^2
index=where(denominator eq 0,count)
if count gt 0 then denominator(index)=1
v=nominator/denominator
if count gt 0 then v(index)=0
;
if normal then norm_factor=4.0/sinc(float(k)/n0)^2 $
else norm_factor=1.0
return,v*norm_factor
;
end