X-Ray Phase Contrast Volume Reconstructions

Hi all,

I am doing a bit of research in x-ray phase contrast imaging. Part of the research is to demonstrate that this new imaging technique can be used for soft tissue volume reconstruction. The image slices demonstrate that glandular breast tissue can be clearly distinguished from the denser in situ ductal carcinoma. This is however not the case for the volume rendering and I believe the issue stems from the pixel values in the slices ranging from 1.2x10^-7 to 1.75x10^-7. I attached an image and will also attach a link to a drop box with the image slices. If anyone has any suggestions on how to improve the volume rendering please let me know.

Thanks,

Matthew

image

https://www.dropbox.com/sh/l80ywcjqd127i6h/AADpiQG6Te7wDx2Q6CgviIQ9a?dl=0

The data seems not to have exported right. I get lots of these warnings when I try to load:

TIFFReadDirectory: Warning, Unknown field with tag 50839 (0xc697) encountered.

Also when I spot check one image I see that the values are very small as you say. If you can export the data another way so it loads as a volume you can rescale the pixel values in Slicer pretty easily (a = array('Tiff*'); a*= 1e9 then save and reload the file to get a better window/level).

I’ve checked this now and the problem is that the folder contains two image series, each with different slice size. Since file name prefixes are almost the same, Slicer tries to load them as one stack and fails.

You can load the two series if you separate the files into two folders:

  • TiffSaver-tomo0111.nrrd … TiffSaver-tomo0476.tif
  • TiffSaver-tomo00000.tif … TiffSaver-tomo00180.tif

I would recommend to use the new utility functions for getting and updating the voxel array:

volNode = getNode('Tiff*')
a=arrayFromVolume(volNode)
a*= 1e10
arrayFromVolumeModified(volNode)

Thanks so much, apologies for the incorrect images being uploaded as well. I noticed that there also some images missing too. I can upload the correct tiff files when I get into work on Monday. The advice you’ve given me already has been a great help.

I have the link for the correct series of slices below https://www.dropbox.com/s/p9tq32mdcpkyaav/cropped3.tif?dl=0

I am new to 3-D slicer is there a document that describes how to access and use the utility functions? I can not find anything on it.

Kind regards,

Matthew

The functions are documented here: slicer — 3D Slicer documentation

1 Like