8000 discrepancy in voxel size/FOV · Issue #147 · LLNL/LEAP · GitHub
[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Skip to content

discrepancy in voxel size/FOV #147

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

Open
karimi-2024 opened this issue Feb 13, 2025 · 2 comments
Open

discrepancy in voxel size/FOV #147

karimi-2024 opened this issue Feb 13, 2025 · 2 comments

Comments

@karimi-2024
Copy link

When I try to set leap image volume parameters using an initialState.sct file, or when I set the parameters using setDefaultVolume(), I get a slight difference. Consider an image that is S x S, and the pixelsize is P

When using initialState.sct, the FOV = S x P, but when I use setDefaultVolume(), the FOV =(S-1) x P.

@kylechampley
Copy link
Collaborator

Are you saying that there is a discrepancy between LTT and LEAP? Is it that the default volume gives different number of voxels or is it that the FOV mask is different?

The code that calculates the default volume is here:
https://github.com/LLNL/LEAP/blob/main/src/parameters.cpp#L786

And the code that calculates the FOV is here:
https://github.com/LLNL/LEAP/blob/main/src/parameters.cpp#L397

@karimi-2024
Copy link
Author
karimi-2024 commented Feb 14, 2025

Hi Kyle, sorry, yes I should have been much clearer. In one method, I am using LTT's loadsct, followed by setLEAPfromLTT. In the second method, I am setting the values with leapct.set_modularbeam(). Then, leapct.print_print() prints out 1 pixel-width different FOVs from the two methods (outputs are below).

The 1-pixel difference changes the projections. If define an object (f) with Leap primitives, and forward project it to projections (g). f is different relative to g in the two methods. That's the 6149 reason I stumbled onto this discrepancy. I will start looking into the links and let you know what I find. Thanks for the response and links! --Seemeen

leapct2.print_parameters()

======== CT Modular-Beam Geometry ========
number of angles: 5
number of detector elements (rows, cols): 600 x 600
detector pixel size: 0.000120 mm x 0.000120 mm
center detector pixel: 299.500000, 299.500000
mean sod = 10000.017578 mm
mean sdd = 10010.017578 mm

======== CT Volume ========
number of voxels (x, y, z): 600 x 600 x 600
voxel size: 0.000120 mm x 0.000120 mm x 0.000120 mm
FOV: [-0.035964, 0.035964] x [-0.035964, 0.035964] x [-0.035964, 0.035964]

======== Processing Settings ========
GPU processing on device 0
GPU with least amount of memory: 6.278028 GB

True

leapct.print_parameters()

======== CT Modular-Beam Geometry ========
number of angles: 5
number of detector elements (rows, cols): 600 x 600
detector pixel size: 0.000120 mm x 0.000120 mm
center detector pixel: 299.500000, 299.500000
mean sod = 10000.017578 mm
mean sdd = 10010.017578 mm

======== CT Volume ========
number of voxels (x, y, z): 600 x 600 x 600
voxel size: 0.000120 mm x 0.000120 mm x 0.000120 mm
FOV: [-0.036000, 0.036000] x [-0.036000, 0.036000] x [-0.036000, 0.036000]

======== Processing Settings ========
GPU processing on device 0
GPU with least amount of memory: 6.278028 GB

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
0