In [1]:
import ifg.shis.rdr.valid as val
from pyfbf import Workspace
import numpy as np
W = Workspace('scatter_test/rdr20170323T170847end20170323T233859sdr20170324T194155')
In [2]:
%matplotlib inline
from pylab import *
In [3]:
debug = {}
ok, seg = val.segments(W, debug=debug)
In [58]:
figure(figsize=(20,10))
plot(seg)
grid()
title('segments unit test for sh170323')
Out[58]:
In [62]:
figure(figsize=(20,10))
# plot(seg)
smp = W['sceneMirrorPosition'][:]
segments = W['segments'][:]
gudrex = np.argwhere(smp==ord('E'))
plot(gudrex, segments[gudrex], color='green')
grid()
title('segments.int4 test for E views sh170323')
Out[62]:
In [5]:
min(seg)
Out[5]:
In [6]:
' '.join(str(x) for x in sorted(set(seg)))
Out[6]:
In [7]:
W['Roll'][np.argwhere(seg==-1).squeeze()]
Out[7]:
In [8]:
figure(figsize=(20,10))
plot(debug['flat'])
# plot(W['Roll'][:] > 2.0)
Out[8]:
In [16]:
figure(figsize=(20,10))
plot(debug['turning'])
# plot(W['Roll'][:] > 2.0)
Out[16]:
In [18]:
figure(figsize=(20,10))
plot(W['Roll'][:])
plot(W['Pitch'][:])
legend(['Roll', 'Pitch'])
grid()
In [19]:
figure(figsize=(20,10))
s,e = 10000,25000
plot(W['Roll'][s:e])
plot(W['Pitch'][s:e])
plot(np.sign(seg[s:e]) * 20)
plot(np.require(ok[s:e], dtype=np.int32) * 25)
legend(['Roll', 'Pitch', 'segments sign', 'segments mask'])
grid()
title('Attitude excerpt sh170323'); xlabel('aircraftTime UTC hours'); ylabel('degrees')
Out[19]:
In [37]:
figure(figsize=(20,10))
s,e = 11000,23000
plot(W['Roll'][s:e])
plot(W['Pitch'][s:e])
plot(np.sign(seg[s:e]) * 20)
plot(np.require(ok[s:e], dtype=np.int32) * 25)
legend(['Roll', 'Pitch', 'segments sign', 'segments mask'])
grid()
title('Attitude excerpt sh170323'); xlabel('record offset'); ylabel('degrees')
Out[37]:
In [23]:
from ifg.shis.rdr.valid import diffmask
limits = val.DEFAULT_SEGMENT_LIMITS
In [24]:
attitude_changing = diffmask(abs(np.diff(pitch)) > limits.max_delta_pitch) | diffmask(abs(np.diff(roll)) > limits.max_delta_roll)
attitude_stable = ~attitude_changing
steady_not_level = ~level & ~attitude_changing
In [72]:
figure(figsize=(20,10))
s,e = 11000,23000
plot(W['Roll'][s:e])
plot(W['Pitch'][s:e])
plot(np.sign(W['segments'][s:e]) * 20)
plot(np.require(ok[s:e], dtype=np.int32) * 25)
legend(['Roll', 'Pitch', 'segments.int4 sign', 'segments mask'])
grid()
title('Attitude excerpt sh170323'); xlabel('record offset'); ylabel('degrees')
Out[72]:
In [25]:
figure(figsize=(20,10))
plot(pitch[:3000])
plot(roll[:3000])
plot(attitude_changing[:3000]*25)
plot(steady_not_level[:3000]*25)
legend(['pitch', 'roll', 'att_change', 'steady_not_level'])
Out[25]:
In [26]:
plot(W['intactMirrorSequences'][:])
Out[26]:
In [45]:
lat,lon,segments = W['instrumentLatitude'][:], W['instrumentLongitude'][:], W['segments'][:]
s,e = 11000, 23000
rex = np.argwhere(lon[s:e]!=0).squeeze() + s
gudrex = np.argwhere((lon[s:e]!=0) & (segments[s:e]!=0)).squeeze() + s
badrex = np.argwhere((lon[s:e]!=0) & (segments[s:e]==0)).squeeze() + s
figure(figsize=(20,15))
plot(lon[gudrex], lat[gudrex], '.', color='blue')
plot(lon[badrex], lat[badrex], '.', color='red')
Out[45]:
In [46]:
lat,lon,segments = W['instrumentLatitude'][:], W['instrumentLongitude'][:], W['segments'][:]
s,e = 11000, 23000
rex = np.argwhere(lon[s:e]!=0).squeeze() + s
gudrex = np.argwhere((lon[s:e]!=0) & (seg[s:e]!=0)).squeeze() + s
badrex = np.argwhere((lon[s:e]!=0) & (seg[s:e]==0)).squeeze() + s
figure(figsize=(20,15))
plot(lon[gudrex], lat[gudrex], '.', color='blue')
plot(lon[badrex], lat[badrex], '.', color='red')
Out[46]:
In [70]:
cal = val.calibration(W)
lat,lon,segments,smp = W['instrumentLatitude'][:], W['instrumentLongitude'][:], W['segments'][:], W['sceneMirrorPosition'][:]
s,e = 11000, 23000
rex = np.argwhere(lon[s:e]!=0).squeeze() + s
gudrex = np.argwhere((lon[s:e]!=0) & (seg[s:e]!=0) & (cal[s:e]) & (smp[s:e]==ord('E'))).squeeze() + s
badrex = np.argwhere((lon[s:e]!=0) & (seg[s:e]!=0) & (~cal[s:e]) & (smp[s:e]==ord('E'))).squeeze() + s
gonrex = np.argwhere((lon[s:e]!=0) & (segments[s:e]==0) & (smp[s:e]==ord('E'))).squeeze() + s
figure(figsize=(20,15))
plot(lon[gonrex], lat[gonrex], '.', color='orange')
plot(lon[gudrex], lat[gudrex], '.', color='blue')
plot(lon[badrex], lat[badrex], '.', color='red')
grid()
title('segmentize!=0 & calibration for earth views %d:%d %s' % (s,e,W.path))
legend(['segments==0', 'Earth segmentize!=0 & cal OK', 'earth segmentize!=0 & bad cal'])
Out[70]:
In [55]:
cal = val.calibration(W)
lat,lon,segments,smp = W['instrumentLatitude'][:], W['instrumentLongitude'][:], W['segments'][:], W['sceneMirrorPosition'][:]
s,e = 11000, 23000
rex = np.argwhere(lon[s:e]!=0).squeeze() + s
gudrex = np.argwhere((lon[s:e]!=0) & (segments[s:e]!=0) & (cal[s:e]) & (smp[s:e]==ord('E'))).squeeze() + s
badrex = np.argwhere((lon[s:e]!=0) & (segments[s:e]!=0) & (~cal[s:e]) & (smp[s:e]==ord('E'))).squeeze() + s
figure(figsize=(20,15))
plot(lon[gudrex], lat[gudrex], '.', color='blue')
plot(lon[badrex], lat[badrex], '.', color='red')
title('segments.int4!=0 & calibration for earth views')
legend(['good cal', 'bad cal'])
Out[55]:
In [ ]: