BMD computation

164 views
Skip to first unread message

Mahmoudi Sidi Ahmed

unread,
Mar 2, 2015, 4:43:39 AM3/2/15
to bonej-users-a...@googlegroups.com
Dear all,

I would like to computer the Bone Mineral Density (BMD) using DICOM CT scans. However, I don't find this functionnality within BoneJ.

I found the pQCT Density Distribtion but the latter requires pQCT snacs.

Do you have any idea for computing BMD in my case ?

Best regards,

Sidi, 

Michael Doube

unread,
Mar 2, 2015, 4:54:33 AM3/2/15
to bonej-users-a...@googlegroups.com
Dear Sidi,

This is usually achieved by scanning some phantoms of known 'BMD' then applying a calibration function to your image using ImageJ's Analyze > Calibrate... command. Be aware that the relation between CT pixel value and BMD may be a non-linear, and maybe even piecewise, function.

Then, the pixel values will be converted to BMD values for your measurements.

'BMD' measurements are highly scale dependent: low resolution techniques such as CT, pQCT and DXA are strongly influenced by BV/TV, i.e. how much marrow space there is & the size of the mineralised features in relation to the sample spacing (~= pixel 'size'). Meanwhile, high resolution techniques like qBSE SEM measure mineralisation of the bone extracellular matrix alone, unconfounded by unmineralised cell and marrow spaces because the sample spacing is on the order of 10s of nm (at least, submicron).

Michael

Michael Doube

unread,
Mar 2, 2015, 5:16:22 AM3/2/15
to bonej-users-a...@googlegroups.com
And I forgot to mention: ordinary CT images are typically full of all sorts of artefacts relating to X-ray physics (e.g. beam hardening) and reconstruction technique (e.g. edge sharpening filters), further reducing the linearity between pixel value and BMD. Be careful!

Michael
> --
> You received this message because you are subscribed to the Google Groups "BoneJ Users and Developers" group.
> To unsubscribe from this group and stop receiving emails from it, send an email to bonej-users-and-dev...@googlegroups.com <mailto:bonej-users-and-dev...@googlegroups.com>.
> To post to this group, send email to bonej-users-a...@googlegroups.com <mailto:bonej-users-a...@googlegroups.com>.
> Visit this group at http://groups.google.com/group/bonej-users-and-developers.
> For more options, visit https://groups.google.com/d/optout.

Ali Ghasem Zadeh

unread,
Mar 2, 2015, 5:29:59 AM3/2/15
to bonej-users-a...@googlegroups.com
Dear Sidi,
If you are using CTScaner , most of them have  got a specific phantom which called QC reference phantom  (with known volumetric density).
You can plot HU and real density on X and Y axis and then you will be able to calculate the function of HU and volumetric density.
Cheers
Ali

Sent from Samsung Mobile
To unsubscribe from this group and stop receiving emails from it, send an email to bonej-users-and-dev...@googlegroups.com.
To post to this group, send email to bonej-users-a...@googlegroups.com.

Mahmoudi Sidi Ahmed

unread,
Mar 2, 2015, 10:31:03 AM3/2/15
to bonej-users-a...@googlegroups.com
Thank you Michael for your answer,

I would like to know if Image J allows to calculate BMD for one image only ? In fact, I want to calculate BMD for the whole image sequence (DICOM files). Is it mandatory to select a ROI ?

Finally, I don't see how to include the phantom values within Image J.

Thanks in advance,

Sidi,

Michael Doube

unread,
Mar 2, 2015, 10:33:04 AM3/2/15
to bonej-users-a...@googlegroups.com
Dear Sidi,

Please take the time to read the ImageJ manual, which will help to answer these questions.

http://rsbweb.nih.gov/ij/docs/guide/

Michael

Mahmoudi Sidi Ahmed

unread,
Mar 3, 2015, 4:24:30 AM3/3/15
to bonej-users-a...@googlegroups.com, al...@unimelb.edu.au
Thank you Ali for your answer,


I would like to know if Image J allows to calculate BMD for one image only ? In fact, I want to calculate BMD for the whole image sequence (DICOM files). Is it mandatory to select a ROI ?

Finally, I don't see how to include the phantom values within Image J.

Thanks in advance,

Sidi,

> To unsubscribe from this group and stop receiving emails from it, send an email to bonej-users-and-developers+unsub...@googlegroups.com <mailto:bonej-users-and-developers+unsubscribe@googlegroups.com>.
> To post to this group, send email to bonej-users-a...@googlegroups.com <mailto:bonej-users-and-devel...@googlegroups.com>.

> Visit this group at http://groups.google.com/group/bonej-users-and-developers.
> For more options, visit https://groups.google.com/d/optout.

--
You received this message because you are subscribed to the Google Groups "BoneJ Users and Developers" group.
To unsubscribe from this group and stop receiving emails from it, send an email to bonej-users-and-developers+unsub...@googlegroups.com.

Mahmoudi Sidi Ahmed

unread,
Mar 6, 2015, 9:45:14 AM3/6/15
to bonej-users-a...@googlegroups.com
Dear all,

A suggested, I'm using The Image j calibration to compute BMD. In fact, after selecting the pjhantom values and launching "calibrate" I obtain a graph showing different pixel values with the corresponding densities.
 However, I would like to obtain the BMD values that corespond the input image sequence. 

Any idea ?

Best regards, 

Sidi,
Reply all
Reply to author
Forward
0 new messages