def get_pmf(self, point):
"""Produces an ODF from a SH bootstrap sample"""
where_dwi = self.where_dwi
bootdata = trilinear_interpolate4d(self.data, point)
dwidata = bootdata[where_dwi]
dwidata = shm.bootstrap_data_voxel(dwidata, self.H, self.R)
bootdata[where_dwi] = dwidata
fit = self.model.fit(bootdata)
return fit.odf(self.sphere)
def test_bootstrap_array():
B = np.array([[4, 5, 7, 4, 2.], [4, 6, 2, 3, 6.]])
H = hat(B.T)
R = np.zeros((5, 5))
d = np.arange(1, 6)
dhat = np.dot(H, d)
assert_array_almost_equal(bootstrap_data_voxel(dhat, H, R), dhat)
assert_array_almost_equal(bootstrap_data_array(dhat, H, R), dhat)
def test_bootstrap_array():
B = np.array([[4, 5, 7, 4, 2.],
[4, 6, 2, 3, 6.]])
H = hat(B.T)
R = np.zeros((5, 5))
d = np.arange(1, 6)
dhat = np.dot(H, d)
assert_array_almost_equal(bootstrap_data_voxel(dhat, H, R), dhat)
assert_array_almost_equal(bootstrap_data_array(dhat, H, R), dhat)
