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

AttenuationConv data should be in units of cm-1 #44

Open
robbietuk opened this issue Jul 30, 2020 · 1 comment
Open

AttenuationConv data should be in units of cm-1 #44

robbietuk opened this issue Jul 30, 2020 · 1 comment
Labels
good first issue Good for newcomers help wanted Extra attention is needed
Milestone

Comments

@robbietuk
Copy link
Collaborator

robbietuk commented Jul 30, 2020

In

10
0 0 Air
114 114 Lung
115 383 Body
384 384 Adipose
385 394 Body
395 395 Muscle
396 402 Body
403 403 Body
404 420 Body
421 1000 SpineBone

the attenuation coefficient ranges are not intuitive. STIR uses attenuation coefficient units of cm^-1.
The input attenuation voxelised phantom for this project should use this attenuation coefficient unit for future STIR processing.

Body tissue has an approximate attenuation coefficient of 0.096cm^-1 (i.e. water).
In

stir_math --including-first --times-scalar 10000 $AttenuationFilenameGATE".hv" $AttenuationFilename

the interfile attenuation image is multipled by 10,000 to create a GATE interfile image, which would corespond to a AttenuationConv of 960 for body tissue/water. This would corresponds to SpineBone in GATE.

Two options:

  1. change the multiplier in In

    stir_math --including-first --times-scalar 10000 $AttenuationFilenameGATE".hv" $AttenuationFilename
    to not be 10,000 but some scaling value that maps between the values in an XCAT and the expected GATE materials.

  2. Change the values in

    10
    0 0 Air
    114 114 Lung
    115 383 Body
    384 384 Adipose
    385 394 Body
    395 395 Muscle
    396 402 Body
    403 403 Body
    404 420 Body
    421 1000 SpineBone

    to better align the multiplied coefficients.

In my opinion, 2. is the better method.
To do this, we need to use something like XCAT to assertain the attenuation coefficients ranges of different body materials, multiply by 10,000, and update the AttenuationConv.dat file. Should be done for #37

@robbietuk robbietuk added help wanted Extra attention is needed good first issue Good for newcomers labels Jul 30, 2020
@robbietuk robbietuk added this to the v1.0 milestone Jul 30, 2020
@robbietuk robbietuk mentioned this issue Jul 30, 2020
9 tasks
@KrisThielemans
Copy link
Collaborator

I agree. conversion from XCAT to STIR is one thing. Conversion from STIR to GATE is another. They should be separate. This table should work for attenuation images generated by STIR really.

I think handling XCAT specifics belongs in https://github.com/UCL/STIR-GATE-Connection/tree/master/VoxelisedSimulation/SubScripts/CreateXCATImages (although that directory is in a bit a funny place as it doesn't contain just scripts., it might be better to move it one level up or so)

@robbietuk robbietuk mentioned this issue Nov 16, 2020
8 tasks
@robbietuk robbietuk modified the milestones: v1.0, v1.1 Nov 16, 2020
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
good first issue Good for newcomers help wanted Extra attention is needed
Projects
None yet
Development

No branches or pull requests

2 participants