Thanks a lot for replying.
I here attached a spreadsheet containing a group of Monte Carlo calculated dose data and film measurement data. Film data have a higher resolution.
Any suggestions to do gamma analysis for these data?
So Nice!
Thank you Simmon and @nko41.
Yes, I have read the posts by you and nko41.
I have also successfully implemented the example you provided. https://docs.pymedphys.com/lib/howto/gamma/from-dicom.html
It’s a really helpful tool
BTW, I have some comments and suggestions about this example.
I will post them later.
From what I have done, I would suggest that you keep your coordinates for both film data and MC dose data the same, I assume that you are evaluating the film data against the MC dose data? Maybe you can try expanding your MC dose data so that each dose data matches each point in the film data. So for example, for the coordinate -49.5 in the MC dose data, you could expand that from -49.5 to -48.6 and use the original value (5.68E-06) in the MC dose data for all those points. You could do this for all the points in your MC dose data and evaluate your film dose data against the reference dose data or vice versa.
Actually, I would recommend to not do the interpolation yourself, instead, individually provide the different coordinate systems into their respective coordinate argument values.
Internally, the gamma code interpolates between the values you provide it in order to do its own searching. To get the most accurate gamma result I would recommend leaving the interpolation to the gamma internals.
Yup, that makes sense. I suspect I was thinking in “world coordinates” when I did those axes. But patient coordinates makes far more sense in that case.
Might you be willing to make the fix yourself? You can open up the file in GitHub over at pymedphys/from-dicom.ipynb at main · pymedphys/pymedphys · GitHub and then it is possible to download the notebook file, make the changes, and then re-upload the file within a pull request for approval.
It’s a bit of work to give it a shot, but the more hands on deck building PyMedPhys together the better PyMedPhys becomes for everyone
I also made some other changes.
For example, the current script always shows the title of the gamma figures with “local gamma” disregarding the value of gamma_options[‘local_gamma’]
I fixed the figure titles to make them consistent with this gamma_options[‘local_gamma’]
I also fixed the labels overlapping issue between the two lower sub plots.