Feature request to reproduce "Synthetic CT example" from MatRAD

[mlapaeva]

Dear all,

First of all, thank you so much for making matRad in Python a reality!
I would like to reproduce the Synthetic CT (sCT) example from this pull request in Python(e0404/matRad#827).

For that, in the PyRadPlan example, a few functionalities are missing:

1. DICOM import → This is required to import CT and sCT .dcm images.
2. DVH and QI extraction → The main goal of the use case is to compare plans and their DVH parameters calculated on CT and sCT. So, it would be great to easily extract DVH.
   In matRad, we have:

resultGUI = matRad_planAnalysis(resultGUI, ct, cst, stf, pln);
% Get the plan parameters
dvh = resultGUI.dvh;
qi = resultGUI.qi;

3. Calculate dose direct function → In our workflow, we first calculate and optimize the plan for the real CT image, then extract the optimized weights, replace the image with the Synthetic CT image, and calculate the dose using the previous weight parameters. It would be great to do this in PyRadPlan.
   In matRad, we have:

weights = resultGUI.w;
resultGUI = matRad_calcDoseDirect(ct, stf, pln, cst, weights);

Could you please let us know when this functionality will be available?

Regards,
Mariia

[wahln]

Thanks for the request:

1. we could start without dicom import by just modifying the CT HUs directly, however @g-stanic is working on loading dicom at the moment.
2. yes, not in there yet, but @g-stanic is looking into that as well so we hopefully have an update there soon. So we just wait for this (probably my review of the code is the bottleneck).
3. this function exists with a new name and signature. In pyRadPlan (as well as in matRad >= 3.1, however it also accepts the older one), it is called calc_dose_forward(ct,cst,stf,pln,w)