k-Wave Toolbox

hounsfield2density

Convert Hounsfield units to density

Syntax

density = hounsfield2density(ct_data)
density = hounsfield2density(ct_data, plot_fitting)
[density, sound_speed] = hounsfield2density(ct_data)
[density, sound_speed] = hounsfield2density(ct_data, plot_fitting)

Description

hounsfield2density converts Hounsfield units to units of density [kg/m^3] based on the experimental data given by Schneider, U., Pedroni, E., and Lomax A., "The calibration of CT Hounsfield units for radiotherapy treatment planning," Phys. Med. Biol., 41, pp. 111-124 (1996). The conversion is made using a four-part piece-wise linear fit to the data (see figure below). For soft-tissue, the approximate sound speed can also be returned using the empirical relationship given by Mast, T. D., "Empirical relationships between acoustic parameters in human soft tissues," Acoust. Res. Lett. Online, 1(2), pp. 37-42 (2000).

Inputs

ct_data

CT data in Hounsfield units to convert to density

Optional Inputs

plot_fitting

Boolean controlling whether the original data points and fitting is plotted (default = false)

Outputs

`density`	density in kg/m^3
`sound_speed`	sound speed in m/s

h5compare

interpCartData