Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

cbct post processing #22

Closed
hyaihjq opened this issue Feb 28, 2024 · 8 comments
Closed

cbct post processing #22

hyaihjq opened this issue Feb 28, 2024 · 8 comments

Comments

@hyaihjq
Copy link

hyaihjq commented Feb 28, 2024

9ebaea98c6e5d5d9b2f04d451cbcbe8a

What post-processing can be used to convert the right image into the left image?thanks

@kylechampley
Copy link
Collaborator

It's hard to say without knowing the source of the data and how the reconstruction was performed, but the artifacts I am seeing are most likely caused by scatter and/ or beam hardening. LEAP does have beam hardening correction algorithms and these are applied to the data prior to reconstruction. If you can provide me more information about this image I might be able to help out more.

@hyaihjq
Copy link
Author

hyaihjq commented Feb 28, 2024

It's hard to say without knowing the source of the data and how the reconstruction was performed, but the artifacts I am seeing are most likely caused by scatter and/ or beam hardening. LEAP does have beam hardening correction algorithms and these are applied to the data prior to reconstruction. If you can provide me more information about this image I might be able to help out more.

thanks for your reply. The right image is the result of reconstruction, and the left image is the visualization result after the software is opened. It can be determined that it is not related to reconstruction, but only the difference in post-processing.

@kylechampley
Copy link
Collaborator

Hmm, interesting. I guess you could try some denoising method. I’d check the visualization software manual.

@hyaihjq
Copy link
Author

hyaihjq commented Feb 29, 2024

Hmm, interesting. I guess you could try some denoising method. I’d check the visualization software manual.

thanks,Is there any way to enhance the contrast of the teeth area while reducing the contrast of the soft tissue?

@kylechampley
Copy link
Collaborator

Ya, TV denoising will do this. Try this:
f_denoised = leapct.diffuse(f, delta, numIter)

where f is the reconstruction volume and numIter is the number of iterations (use 10-50 iterations). You'll have to find the value of delta yourself. The delta parameter specifies a transition value where
if the difference of two neighboring voxels are less than delta, then they will be blurred together, reducing contrast and noise
on the other hand if the difference of two neighboring voxels are larger than delta, then they will have little influence on each other and the edge between these will be preserved.

Thus you should set the delta parameter to be less than the difference between the teeth and the soft tissue. I usually manually tune delta, but as a rule of thumb, I start with this guess:
delta = ((typical voxel value of material A) - (typical voxel value of material B)) / 20
where material A and B are materials that you want to have clear separation.

Good luck!

@hyaihjq
Copy link
Author

hyaihjq commented Feb 29, 2024

Ya, TV denoising will do this. Try this: f_denoised = leapct.diffuse(f, delta, numIter)

where f is the reconstruction volume and numIter is the number of iterations (use 10-50 iterations). You'll have to find the value of delta yourself. The delta parameter specifies a transition value where if the difference of two neighboring voxels are less than delta, then they will be blurred together, reducing contrast and noise on the other hand if the difference of two neighboring voxels are larger than delta, then they will have little influence on each other and the edge between these will be preserved.

Thus you should set the delta parameter to be less than the difference between the teeth and the soft tissue. I usually manually tune delta, but as a rule of thumb, I start with this guess: delta = ((typical voxel value of material A) - (typical voxel value of material B)) / 20 where material A and B are materials that you want to have clear separation.

Good luck!

thanks for your help! I will try your methed.

@kylechampley
Copy link
Collaborator

Any luck? Do you still want to keep this issue open?

@hyaihjq
Copy link
Author

hyaihjq commented Mar 18, 2024

Any luck? Do you still want to keep this issue open?

no, thanks

@hyaihjq hyaihjq closed this as completed Mar 18, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants