Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.apo.nmsu.edu/35m_operations/35m_manual/Instruments/CorMass/Files/magicsub2.pro Дата изменения: Thu Jul 8 07:52:44 2004 Дата индексирования: Sun Apr 10 09:20:50 2016 Кодировка: Поисковые слова: внешние планеты

pro magicsub2,oldim,procim

procim=fltarr(256,256)
unbias=fltarr(64)

for j=192,255 do begin

a=sort(oldim[*,j]) ;a contains indices of ascending ordered row
unbias[j-192]=oldim[a[64],j] ;choose mode

endfor

unbias=reverse(unbias)
mn=min(unbias)
unbias=unbias-mn

x=indgen(64)
result=poly_fit(x,unbias,2) ;fit to 2nd order polynomial
result=reform(result)

x=indgen(128)
y=result[0] + result[1]*x + result[2]*x^2

x=indgen(256)
ytot=[y,y]
ytot=reverse(ytot)
ytot=ytot+mn+result[0]

for h=0,255 do begin

procim[*,h]=oldim[*,h]-ytot[h] ;make reset-decay correction

endfor

return
end