User Tools

Site Tools


code:work:lvisf:processing:idl:start

This is an old revision of the document!


IDL Processing

2017.04.07

Stationary RX

Roll 25000

@read_lvisf_merged.pro
; select file #1
v=where(c.system_status GT 100 and c.system_status LT 500)
pa = peak(c.txawaveform)
pb = peak(c.txbwaveform)
pc = peak(c.txcwaveform)
mtitle='2017.03.31 CHA Roll=25000 RX=12104 File=131830'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pa(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHB Roll=25000 RX=12104 File=131830'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pb(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHC Roll=25000 RX=12104 File=131830'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pc(v)),psym=2,symsize=0.3

Roll 35000

@read_lvisf_merged.pro
; select file #2
v=where(c.system_status GT 100 and c.system_status LT 500)
pa = peak(c.txawaveform)
pb = peak(c.txbwaveform)
pc = peak(c.txcwaveform)
mtitle='2017.03.31 CHB Roll=35000 RX=24712 File=131930'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pa(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHB Roll=35000 RX=24712 File=131930'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pb(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHB Roll=35000 RX=24712 File=131930'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pc(v)),psym=2,symsize=0.3

Roll 45000

@read_lvisf_merged.pro
; select file #3
i=where(c.time_gps_sod GT 1)
plot,c(i).time_gps_sod,c(i).system_status
v=where(c.system_status GT 100 and c.system_status LT 500 and c.time_gps_sod LT 48080)
pa = peak(c.txawaveform)
pb = peak(c.txbwaveform)
pc = peak(c.txcwaveform)
mtitle='2017.03.31 CHA Roll=45000 RX=35836 File=132030'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pa(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHB Roll=45000 RX=35836 File=132030'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pb(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHC Roll=45000 RX=35836 File=132030'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pc(v)),psym=2,symsize=0.3

Roll 55000

@read_lvisf_merged.pro
; select file #3
c1=c
@read_lvisf_merged.pro
; select file #4
c2=c
c=[c1,c2] & c1=0 & c2=0
i=where(c.time_gps_sod GT 1)
plot,c(i).time_gps_sod,c(i).system_status
v=where(c.system_status GT 100 and c.system_status LT 500 and c.time_gps_sod GT 48080 and c.time_gps_sod LT 48120)
pa = peak(c.txawaveform)
pb = peak(c.txbwaveform)
pc = peak(c.txcwaveform)
mtitle='2017.03.31 CHA Roll=55000 RX=47340 File=132130'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pa(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHB Roll=55000 RX=47340 File=132130'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pb(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHC Roll=55000 RX=47340 File=132130'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pc(v)),psym=2,symsize=0.3

60000

@read_lvisf_merged.pro
; select file #3
c1=c
@read_lvisf_merged.pro
; select file #4
c2=c
c=[c1,c2] & c1=0 & c2=0
i=where(c.time_gps_sod GT 1)
plot,c(i).time_gps_sod,c(i).system_status
v=where(c.system_status GT 100 and c.system_status LT 500 and c.time_gps_sod GT 48120)
pa = peak(c.txawaveform)
pb = peak(c.txbwaveform)
pc = peak(c.txcwaveform)
mtitle='2017.03.31 CHA Roll=60000 RX=54054 File=132130'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pa(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHB Roll=60000 RX=54054 File=132130'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pb(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHC Roll=60000 RX=54054 File=132130'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pc(v)),psym=2,symsize=0.3
code/work/lvisf/processing/idl/start.1491601704.txt.gz · Last modified: 2017/04/07 21:48 by david