/
mage_findslits.pro
35 lines (28 loc) · 966 Bytes
/
mage_findslits.pro
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
function mage_findslits, traceflat, orderfile = orderfile
; files = ['mage0114.fits','mage0115.fits','mage0116.fits','mage0117.fits','mage0118.fits']
; for i=0, n_elements(files)-1 do begin
; mage_proc, '../Raw/'+files[i], tmp, hdr=hdr
; if i EQ 0 then begin
; trcimg = tmp
; endif else begin
; trcimg += tmp
; endelse;
;
; endfor
mage_proc, traceflat, trcimg, hdr=hdr
tset_slits = mage_traceorders(trcimg, /chk)
ordermask=mage_ordermask(tset_slits)
;;slitmask = long_slits2mask(tset_slits)
;;ordermask = 0 * slitmask
;;ordermask[WHERE(slitmask GT 0)] = -slitmask[WHERE(slitmask GT 0)] + 21L
if keyword_set(orderfile) then begin
mwrfits, ordermask, orderfile, /create
mwrfits, tset_slits, orderfile
endif else begin
file = "Orders.fits"
mwrfits, ordermask, file, /create
mwrfits, tset_slits, file
endelse
return, tset_slits
print, "mage_findslits: all done!"
end