def test_dpy():
fd,fname = mkstemp()
dpw = Dpy(fname,'w')
A=np.ones((5,3))
B=2*A.copy()
C=3*A.copy()
dpw.write_track(A)
dpw.write_track(B)
dpw.write_track(C)
dpw.write_tracks([C,B,A])
dpw.close()
dpr = Dpy(fname,'r')
assert_equal(dpr.version()=='0.0.1',True)
T=dpr.read_tracksi([0,1,2,0,0,2])
print(T)
T2=dpr.read_tracks()
assert_equal(len(T2),6)
dpr.close()
assert_array_equal(A,T[0])
assert_array_equal(C,T[5])
os.remove(fname)
def roi_intersection(fdpy,fatlas,roi_no,froidpy):
dpr=Dpy(fdpy,'r')
T=dpr.read_tracksi(range(10000))
dpr.close()
Troi=[]
wI=get_roi(fatlas,roi_no,0)
for (i,t) in enumerate(T):
if i%1000==0:
print i
if track_roi_intersection_check(t,wI,.5):
Troi.append(t)
print(len(Troi))
dpw=Dpy(froidpy,'w')
dpw.write_tracks(Troi)
dpw.close()
'''
def test_dpy():
fname = 'test.bin'
dpw = Dpy(fname, 'w')
A = np.ones((5, 3))
B = 2 * A.copy()
C = 3 * A.copy()
dpw.write_track(A)
dpw.write_track(B)
dpw.write_track(C)
dpw.write_tracks([C, B, A])
dpw.close()
dpr = Dpy(fname, 'r')
assert_equal(dpr.version() == '0.0.1', True)
T = dpr.read_tracksi([0, 1, 2, 0, 0, 2])
T2 = dpr.read_tracks()
assert_equal(len(T2), 6)
dpr.close()
assert_array_equal(A, T[0])
assert_array_equal(C, T[5])
def test_dpy():
fname = 'test.bin'
with InTemporaryDirectory():
dpw = Dpy(fname, 'w')
A = np.ones((5, 3))
B = 2 * A.copy()
C = 3 * A.copy()
dpw.write_track(A)
dpw.write_track(B)
dpw.write_track(C)
dpw.write_tracks([C, B, A])
dpw.close()
dpr = Dpy(fname, 'r')
assert_equal(dpr.version() == '0.0.1', True)
T = dpr.read_tracksi([0, 1, 2, 0, 0, 2])
T2 = dpr.read_tracks()
assert_equal(len(T2), 6)
dpr.close()
assert_array_equal(A, T[0])
assert_array_equal(C, T[5])
def test_dpy():
fname = 'test.bin'
with InTemporaryDirectory():
dpw = Dpy(fname, 'w')
A = np.ones((5, 3))
B = 2 * A.copy()
C = 3 * A.copy()
dpw.write_track(A)
dpw.write_track(B)
dpw.write_track(C)
dpw.write_tracks(Streamlines([C, B, A]))
all_tracks = np.ascontiguousarray(np.vstack([A, B, C, C, B, A]))
npt.assert_array_equal(all_tracks, dpw.tracks[:])
dpw.close()
dpr = Dpy(fname, 'r')
npt.assert_equal(dpr.version() == u'0.0.1', True)
T = dpr.read_tracksi([0, 1, 2, 0, 0, 2])
T2 = dpr.read_tracks()
npt.assert_equal(len(T2), 6)
dpr.close()
npt.assert_array_equal(A, T[0])
npt.assert_array_equal(C, T[5])
dpw.write_tracks(streamlines)
"""
Write one track
"""
dpw.write_track(streamlines[0])
"""
or one track each time.
"""
for t in streamlines:
dpw.write_track(t)
dpw.close()
"""
Read streamlines directly from the disk using their indices
.. include:: ../links_names.inc
"""
dpr = Dpy('fornix.dpy', 'r')
some_streamlines = dpr.read_tracksi([0, 10, 20, 30, 100])
dpr.close()
print(len(streamlines))
print(len(some_streamlines))
"""
Write many streamlines at once.
"""
dpw.write_tracks(streamlines2)
"""
Write one track
"""
dpw.write_track(streamlines2[0])
"""
or one track each time.
"""
for t in streamlines:
dpw.write_track(t)
dpw.close()
"""
Read streamlines directly from the disk using their indices
.. include:: ../links_names.inc
"""
dpr = Dpy('fornix.dpy', 'r')
some_streamlines = dpr.read_tracksi([0, 10, 20, 30, 100])
dpr.close()
print(len(streamlines))
print(len(some_streamlines))
write many tracks at once
"""
dpw.write_tracks(tracks2)
"""
write one track
"""
dpw.write_track(tracks2[0] * 6)
"""
or one track each time
"""
for t in tracks:
dpw.write_track(t * 3)
dpw.close()
"""
read tracks directly from the disk using their indices
"""
dpr = Dpy('fornix.dpy', 'r')
some_tracks = dpr.read_tracksi([0, 10, 20, 30, 100])
dpr.close()
"""
Number of tracks in before and after
"""
print(len(tracks))
print(len(some_tracks))
"""
or one track each time
"""
for t in tracks:
dpw.write_track(t*3)
dpw.close()
"""
read tracks directly from the disk using their indices
"""
dpr=Dpy('fornix.dpy','r')
some_tracks=dpr.read_tracksi([0,10,20,30,100])
dpr.close()
"""
Number of tracks in before and after
"""
print(len(tracks))
print(len(some_tracks))
"""
.. include:: ../links_names.inc
"""
