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

Allow generation of .3ddose files for transformed EGS_XYZGeometry #1145

Open
blakewalters opened this issue Jun 5, 2024 · 3 comments · May be fixed by #1165
Open

Allow generation of .3ddose files for transformed EGS_XYZGeometry #1145

blakewalters opened this issue Jun 5, 2024 · 3 comments · May be fixed by #1165
Assignees
Milestone

Comments

@blakewalters
Copy link
Contributor

Is your feature request related to a problem? Please describe.
Currently, egs++ supports reading a CT .egsphant file into a rectilinear EGS_XYZGeometry. Using EGS_DoseScoring, a .3ddose file can potentially be generated for a simulation using this EGS_XYZGeometry. However, if, as is almost always the case, the EGS_XYZGeometry created using CT data must be rotated/translated from "CT space" into the simulation coordinate system, a .3ddose file for the transformed geometry can no longer be generated.

Describe the solution you'd like
It would be very useful to be able to generate a .3ddose file for CT data. By default, this data would be output in "CT space" and could be plotted atop the untransformed CT image. There may be a way to also transform the dose data into simulation space.

Alternatives
I suppose one could build their simulation geometry around the original CT data, but this seems unnecessarily complicated and inflexible.

@mchamberland
Copy link
Contributor

mchamberland commented Jun 6, 2024

egs_brachy did it first.

😉

Though, this situation does not come up as often in brachytherapy simulations. Also, I have not handled converting the 3ddose file to "simulation space". It seemed too daunting to me at the time, but I admit I only thought about the problem for about 5 minutes before coming to that conclusion!

@blakewalters
Copy link
Contributor Author

@mchamberland, I've been using a fix that undoes the transformations at the beginning of the run just to determine: 1) no. of voxels, 2) voxel boundaries, and 3) mapping of global dose region -> region no. in EGS_XYZGeometry. This seems to work okay...but it causes a crash in egs_view if the egs_dose ausgab object with the dose geometry isn't commented out before attempting to display.

@blakewalters blakewalters linked a pull request Aug 9, 2024 that will close this issue
@blakewalters
Copy link
Contributor Author

@rtownson and @ftessier: This issue is addressed by #1165, which might be an easy "add" to the 2024 release.

@ftessier ftessier added this to the Release 2024 milestone Aug 16, 2024
@rtownson rtownson linked a pull request Aug 19, 2024 that will close this issue
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging a pull request may close this issue.

3 participants