def test_quant_cross_section_ac(self):
s = self.signal
method = 'cross_section'
zfactors = [20, 50]
factors = utils_eds.zeta_to_edx_cross_section(zfactors, ['Al', 'Zn'])
intensities = s.get_lines_intensity()
res = s.quantification(intensities, method, factors,
absorption_correction=True)
res2 = s.quantification(intensities, method='zeta',
factors=zfactors,
absorption_correction=True)
np.testing.assert_allclose(res[0][1].data, np.ones((2, 2)) * 44.02534,
atol=1e-3)
np.testing.assert_allclose(res2[0][0].data, res[0][0].data, atol=1e-3)
def test_quant_cross_section_ac(self):
s = self.signal
method = 'cross_section'
factors = [3, 5]
intensities = s.get_lines_intensity()
res = s.quantification(intensities, method, factors,
absorption_correction=True)
zfactors = utils_eds.zeta_to_edx_cross_section(factors, ['Al', 'Zn'])
res2 = s.quantification(intensities, method='zeta', factors=[22.402, 21.7132],
absorption_correction=True)
np.testing.assert_allclose(res[0][0].data, np.array(
[[49.4889, 49.4889],
[49.4889, 49.4889]]), atol=1e-3)
np.testing.assert_allclose(res[0][0].data, np.array(
[[49.4889, 49.4889],
[49.4889, 49.4889]]), atol=1e-3)
def test_zeta_vs_cross_section(self):
s = self.signal
factors = [3, 5]
method = 'zeta'
intensities = s.get_lines_intensity()
zfactors = utils_eds.edx_cross_section_to_zeta([3, 5], ['Al', 'Zn'])
factors2 = utils_eds.zeta_to_edx_cross_section(zfactors, ['Al', 'Zn'])
np.testing.assert_allclose(factors, factors2, atol=1e-3)
res = s.quantification(intensities,
method,
factors=utils_eds.edx_cross_section_to_zeta(
[22.402, 21.7132], ['Al', 'Zn']))
res2 = s.quantification(intensities,
method='cross_section',
factors=[22.402, 21.7132])
np.testing.assert_allclose(res[0][0].data, res2[0][0].data, atol=1e-3)
np.testing.assert_allclose(res[0][0].data,
np.ones((2, 2)) * 36.2969,
atol=1e-3)
def test_zeta_to_edx_cross_section(self):
factors = [1079.815272, 162.4378035]
elements = ['Pt', 'Ni']
res = utils_eds.zeta_to_edx_cross_section(factors, elements)
np.testing.assert_allclose(res, [3, 6], atol=1e-3)
def test_zeta_to_edx_cross_section(self):
factors = [1079.815272, 162.4378035]
elements = ['Pt', 'Ni']
res = utils_eds.zeta_to_edx_cross_section(factors, elements)
np.testing.assert_allclose(res, [3, 6], atol=1e-3)
