def test_rhoAsh_equ_correct2(self):
hu = uniform(0,2000)
self.param['calibrationCorrect'] = True
self.param['numCTparam'] = 'single'
self.param['rhoAsha1'] = uniform(0,2000)
self.param['rhoAshb1'] = uniform(0,2000)
rhoQCT, rhoAsh, modulus = _define_equations(self.param)
self.assertEqual(rhoAsh(rhoQCT(hu)), self.param['rhoAsha1'] + (hu*self.param['rhoAshb1']))
def setUp(self):
self.vtk = create_vtk()
self.voxels = [[1, 2], [6, 10, 12], [13, 16, 18, 22], [26]]
self.param = {'integration': 'E',
'rhoQCTa': 0.,
'rhoQCTb': 1.,
'calibrationCorrect': False,
'numEparam': 'single',
'minVal': 0.0001,
'Ea1': 0.,
'Eb1': 1.,
'Ec1': 1.}
self.rhoQCT_equ, self.rhoAsh_equ, self.modulus_equ = _define_equations(self.param)
self.moduli = [[_apply_equations(self.vtk.lookup[i],
self.rhoQCT_equ,
self.rhoAsh_equ,
self.modulus_equ) for i in v] for v in self.voxels]
self.mean_moduli = [mean(m) for m in self.moduli]
def test_rhoAsh_equ_correct3(self):
hu = uniform(0,2000)
self.param['calibrationCorrect'] = True
self.param['numCTparam'] = 'triple'
self.param['rhoAsha1'] = uniform(0,2000)
self.param['rhoAshb1'] = uniform(0,2000)
self.param['rhoAsha2'] = uniform(0,2000)
self.param['rhoAshb2'] = uniform(0,2000)
self.param['rhoAsha3'] = uniform(0,2000)
self.param['rhoAshb3'] = uniform(0,2000)
tmp = [uniform(0,2000), uniform(0,2000)]
self.param['rhoThresh1'] = min(tmp)
self.param['rhoThresh2'] = max(tmp)
self.param['rhoQCTa'] = 0.
self.param['rhoQCTb'] = 1.
rhoQCT, rhoAsh, modulus = _define_equations(self.param)
if rhoQCT(hu) < self.param['rhoThresh1']:
res = self.param['rhoAsha1'] + (self.param['rhoAshb1'] * rhoQCT(hu))
elif (rhoQCT(hu) >= self.param['rhoThresh1']) & (rhoQCT(hu) <= self.param['rhoThresh2']):
res = self.param['rhoAsha2'] + (self.param['rhoAshb2'] * rhoQCT(hu))
else:
res = self.param['rhoAsha3'] + (self.param['rhoAshb3'] * rhoQCT(hu))
self.assertEqual(rhoAsh(rhoQCT(hu)), res)
def test_rhoAsh_equ_correct(self):
hu = uniform(0,2000)
self.param['rhoQCTa'] = uniform(0,2000)
self.param['rhoQCTb'] = uniform(0,2000)
rhoQCT, rhoAsh, modulus = _define_equations(self.param)
self.assertEqual(rhoAsh(rhoQCT(hu)), self.param['rhoQCTa'] + (hu*self.param['rhoQCTb']))
