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