flatfield/off_flat_IDL/wo.pro

pro WO, reihe, nullen, sv=sv
;+
; NAME:
;       WO
; PURPOSE:
;       Finding invalid points in a time series
; EXPLANATION:
;       giving the first and last point of a sqequence of invalid points in a time series
;
;
; CALLING SEQUENCE:
;       WO, reihe, nullen [,sv=]
; INPUTS:
;       reihe: time series ((n) vector)
; OPTIONAL KEYWORS: sv: value of the invalid points (default criterium: finite(..)=0)
; OUTPUTS:
;       nullen: see explanaition ((2, m) vector)
; PROCEDURES USED:
;None
;
; MODIFICATION HISTORY:
;       Package 'SPM level 1.5' written between 2000 and 2003, Richard Wachter, PMOD/WRC
;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

 if n_elements(sv) eq 0 then w=where(finite(reihe) eq 0) else w=where(reihe eq sv)


 lw=(size(w))[1]

 seq=lonarr(2)
 nullen=[0L,0L]


 seq[0]=w[0]

 for k=1L, lw-1 do begin

  if ((w[k] - w[k-1]) gt 1) then begin

   seq[1]=w[k-1]
   nullen=[[nullen],[seq]]
   seq[0]=w[k]

  endif

 endfor

 seq[1]=w[lw-1]
 nullen=[[nullen],[seq]]
 nullen=nullen[*,1:*]

end


Revision:	1.1
Committed:	Fri Feb 18 00:26:40 2011 UTC (12 years, 7 months ago) by richard
Branch:	MAIN
CVS Tags:	Ver_6-0, Ver_6-1, Ver_6-2, Ver_6-3, Ver_6-4, Ver_9-1, Ver_5-14, Ver_5-13, Ver_LATEST, Ver_9-3, Ver_9-41, Ver_9-2, Ver_8-8, Ver_8-2, Ver_8-3, Ver_8-0, Ver_8-1, Ver_8-6, Ver_8-7, Ver_8-4, Ver_8-5, Ver_7-1, Ver_7-0, Ver_9-5, Ver_9-4, Ver_8-10, Ver_8-11, Ver_8-12, Ver_9-0, HEAD
Log Message:	IDL package to calculate offpoint flatfield, dark frames, and bad pixel lists from spacecraft offpoints 2011.02.17
#	Content
1	pro WO, reihe, nullen, sv=sv
2	;+
3	; NAME:
4	; WO
5	; PURPOSE:
6	; Finding invalid points in a time series
7	; EXPLANATION:
8	; giving the first and last point of a sqequence of invalid points in a time series
9	;
10	;
11	; CALLING SEQUENCE:
12	; WO, reihe, nullen [,sv=]
13	; INPUTS:
14	; reihe: time series ((n) vector)
15	; OPTIONAL KEYWORS: sv: value of the invalid points (default criterium: finite(..)=0)
16	; OUTPUTS:
17	; nullen: see explanaition ((2, m) vector)
18	; PROCEDURES USED:
19	;None
20	;
21	; MODIFICATION HISTORY:
22	; Package 'SPM level 1.5' written between 2000 and 2003, Richard Wachter, PMOD/WRC
23	;
24	;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
25
26	if n_elements(sv) eq 0 then w=where(finite(reihe) eq 0) else w=where(reihe eq sv)
27
28
29
30	lw=(size(w))[1]
31
32	seq=lonarr(2)
33	nullen=[0L,0L]
34
35
36	seq[0]=w[0]
37
38	for k=1L, lw-1 do begin
39
40	if ((w[k] - w[k-1]) gt 1) then begin
41
42	seq[1]=w[k-1]
43	nullen=[[nullen],[seq]]
44	seq[0]=w[k]
45
46	endif
47
48	endfor
49
50	seq[1]=w[lw-1]
51	nullen=[[nullen],[seq]]
52	nullen=nullen[,1:]
53
54	end
55
56
57