Hi,
I am trying to use wrap solidify extension. And I am following these steps
But when I am using apply button I am getting this error. Why I am getting this error?
You need to set some content in the segment that you would like to solidify. In CT images, you can segment the skin surface or all the bones using simple thresholding (using Threshold effect).
See these segmentation recipes for step-by-step instructions:
Hi,
I follower the instruction your suggested. When I tried with the outer surface option, the wrap solidify worked. But I sill can see the holes in the bone labels. Screenshopt attached.
Does this functionality helps to fill the holes?
I am sharing my data file here. Could you please look into it. That will be really helpful.
https://drive.google.com/drive/folders/1s1W06mB6EIyTfMsKf86K3-l6zuTimxFx?usp=sharing
If you see a “hole” on a slice then it does not mean that it is actually a hole in 3D. It may be just a concavity that may show up in a hole in a certain intersection, but if you look at it on other slices (hold down Shift key while moving the mouse pointer to see the same point in all image slices, you can also enable crosshair display by clicking on the crosshair icon in the toolbar) then you’ll see what it is.
You may still have holes in the final result at places where the surface becomes very thin. To fix those, you can use Crop volume module to resample the input volume to have finer and isotropic resolution by setting spacing scale to a value <1 (e.g. 0.5).
Note that there are a few issues with your data set:
It seems that your image has been imported from DICOM incorrectly. What software did you use to import it? If it comes from a CBCT then you may need to run the image through DICOM patcher and/or set the correct spacing in Volumes module / Volume information / Image spacing section. Currently, LR and AP spacing seems 2-3x larger than reality (if it is a human skull) and IS spacing is relatively lower (by a factor of 1.5-2x). Fix this before doing anything else with your volume.
The volume seems to be pre-processed by some strange algorithm and made very noisy (lots of quantization noise; CTs are normally 16-bit images, but this one is just 8). After thresholding, you can reduce effect of this noise by removing small speckles by using Smoothing effect with Median method.
Thank you so much for your help. I really appreciate your time. It was not in dicom format. I downloaded the png sequences from below link(you need to open it in Microsoft edge). Which is isotropic. In this link they also have the label images as well of different (bone , Holes). I have labeled with these two labeled images.
I have availed the original ct images downloaded from the website in the shared drive folder.(axct folder)
https://drive.google.com/drive/folders/1s1W06mB6EIyTfMsKf86K3-l6zuTimxFx?usp=sharing
I did a resampling. Okay, I will do the resampling with spacing values 0.5 on your suggestion. I did not do any other pre-processing accept doing resampling. I did use smoothing effect but I have seen a lot of islands and noises in the surface extraction. (The screenshot)
I am getting back to you after doing a threshold and smoothing.
What is the end goal of your project?
To create a deformable surface model of cranium foramina which can be warped onto patient skull ct image to identify the boundary of the cranial foramina. Before this the steps are:
- Multi-material 2-simplex is going to applied for the surface extraction to create triangular bone and hole mesh.
- Decimation is needed to control resolution and to solve problem of producing tractable nonrigid registration–without a coarse-to-fine approach, the fine mesh can become easily entangled in false minima. ( I am going to use ACVD for the decimation by keeping the shared boundary(hole and bone ) intact:: that is why I need a smooth simple bone mesh)
- A deformable model registration at first with bone mesh and later hole mesh with patient ct data.
Images that you downloaded from http://headneckbrainspine.com are quite bad quality screenshots. You can get proper CTs from various medical image databases. For example, check out head CT images in TCIA’s head and neck cancer (CPTAC-HNSCC) colletion. They were acquired 20 years ago, still better quality than the image you got from that training website. Due to better dose efficiency and improvement of image reconstruction methods, most likely you can get even better quality images (lower noise, maybe even higher resolution) from current CT scanners. I would highly recommend to ask for images from your clinical collaborators that are acquired with today’s technology.
If you need the foramina then you need to enable “carve holes” and set the value lower than the smallest holes you want to preserve. If you are interested in segmenting very delicate structures then you may need to oversample your volume. For example, to reliable segment a 2mm hole, you need 0.4mm but preferable about 0.2mm resolution in your binary labelmap. If your input CT does not have this high resolution then you can oversample it using Crop volume module. To handle such high resolution images, you may either need to reduce your region of interest size or use strong computer, with lots of RAM (about 10x more than your cropped volume size).
Thank you so much for your insights. But if I use other ct images, I don’t have the segmented label images for that those ct images. Only the headneckbrainspine has segmented label images for bone and holes. That is the limitation that I cannot use other ct images.
I have added one folder axatlas (segmented label bone image) and axforamenmagnum(segmented hole label image) for example.
https://drive.google.com/drive/folders/1oq9F9exKBsq-Mn7JFFY0ojwaLxUGq4ZB?usp=sharing
You can segment any CT image using Segment Editor module.