-
Notifications
You must be signed in to change notification settings - Fork 11
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
Documentation of pseudoCT processing #43
Comments
I am trying to debug a pseudoCT run and noticed that my max. supported freq is below the trageted 500 kHz. Anyone run into this before? The relevant computation is:
The This is the usual spec I get (which in this case also produces NaNs however):
I presume it's because the minimum sound speed is unreasonable (below that of water), so a problem with the CT-to-tissueproperty transform. |
Update: determining the lower bound of soundspeed as the one of water indeed fixes at least the max. supported frequency issue. The relevant fix is bdf0ff9. Given that soundspeed is directly translated from the estimated density this may only be a superficial fix for a deeper problem. |
The pCT-to-tissue conversion algorithm documented in Yaakub et al. is now implemented in Ref: Yaakub, S. N. et al. Pseudo-CTs from T1-Weighted MRI for Planning of Low-Intensity Transcranial Focused Ultrasound Neuromodulation: An Open-Source Tool. Brain Stimulation. 16. 75–78 (2023). |
For All following steps are applied to the skull mask only.
For both variants, 𝛼 is bounded based on estimates made at 500 kHz (i.e., 𝛼(f); see Aubry, J.-F., 2022 for prior benchmark simulations). We actually require Reference: Adapted from Carpino et al. (2024). Transcranial ultrasonic stimulation of the human amygdala to modulate threat learning. MSc thesis. |
For I have set it to HU_min for now with commit 39e453e. |
Now that we have pseudoCT pseudo-integrated, it would be worthwhile to document how we are processing the related tissue mask (e.g., smoothing choices) and the conversion into
pseudo-Hounsfield
,sound_speed
andalpha
units.The relevant functions are:
The text was updated successfully, but these errors were encountered: