Home Explore Help
Sign In
IGT-Madison
/
WiscPlan_olderMV
2
0
Fork 0
Files Issues 0 Pull Requests 0 Wiki
Tree: 71798cc853
Branches Tags
WiscPlan_old...
/
opt
/
linlsqOpt
/
input
/
dosePaintSetup.asv

dosePaintSetup.asv 4.9 KB
History Raw

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485 
  1. %% script to apply non-uniform prescription in IMRT dose optimization
    2. 

  2. %  prescription matrix 'PTVDosePlus.bin' and 'PTVDoseMinus.bin'
    3. 

  3. 

  4. % addpath('C:\Documents and Settings\Steve\Desktop\daveTPS\photonDoseCalc\scripts_m');
    5. 

  5. % addpath('C:\Documents and Settings\Steve\Desktop\daveTPS\opt\linlsqUtil');
    6. 

  6. 

  7. % open prescription in sparse format matlab structure and create PTV mask
    8. 

  8. %
    9. 

  9. %  presc.x_count (int)
    10. 

  10. %  presc.y_count (int)
    11. 

  11. %  presc.z_count (int)
    12. 

  12. %  presc.Nind    (int)
    13. 

  13. %  presc.non_zero_indices  (int32)
    14. 

  14. %  presc.non_zero_values   (single)
    15. 

  15. 

  16. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV003f\Doses\AV003f_Continuous_cropped.am');
    17. 

  17. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV003f\Doses\AV003f_PVC_Continuous_cropped.am');
    18. 

  18. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV003f\Doses\AV003f_Threshold_cropped.am');
    19. 

  19. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV003f\Doses\AV003f_PVC_Threshold_cropped.am');
    20. 

  20. 

  21. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV004c\Doses\AV004c_Continuous_cropped.am');
    22. 

  22. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV004c\Doses\AV004c_PVC_Continuous_cropped.am');
    23. 

  23. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV004c\Doses\AV004c_Threshold_cropped.am');
    24. 

  24. pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV004c\Doses\AV004c_PVC_Threshold_cropped.am');
    25. 

  25. 

  26. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\HN002b\Doses\HN002b_Continuous_cropped.am');
    27. 

  27. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\HN002b\Doses\HN002b_PVC_Continuous_cropped.am');
    28. 

  28. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\HN002b\Doses\HN002b_Threshold_cropped.am');
    29. 

  29. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\HN002b\Doses\HN002b_PVC_Threshold_cropped.am');
    30. 

  30. 

  31. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV001c\Doses\AV001c_Continuous_cropped.am');
    32. 

  32. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV001c\Doses\AV001c_PVC_Continuous_cropped.am');
    33. 

  33. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV001c\Doses\AV001c_Threshold_cropped.am');
    34. 

  34. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV001c\Doses\AV001c_PVC_Threshold_cropped.am');
    35. 

  35. 

  36. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV001f\Doses\AV001f_Continuous_cropped.am');
    37. 

  37. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV001f\Doses\AV001f_PVC_Continuous_cropped.am');
    38. 

  38. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV001f\Doses\AV001f_Threshold_cropped.am');
    39. 

  39. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV001f\Doses\AV001f_PVC_Threshold_cropped.am');
    40. 

  40. 

  41. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\HN003c\Doses\HN003c_Continuous_cropped.am');
    42. 

  42. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\HN003c\Doses\HN003c_PVC_Continuous_cropped.am');
    43. 

  43. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\HN003c\Doses\HN003c_Threshold_cropped.am');
    44. 

  44. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\HN003c\Doses\HN003c_PVC_Threshold_cropped.am');
    45. 

  45. 

  46. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV004f\Doses\AV004f_Continuous_cropped.am');
    47. 

  47. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV004f\Doses\AV004f_PVC_Continuous_cropped.am');
    48. 

  48. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV004f\Doses\AV004f_Threshold_cropped.am');
    49. 

  49. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV004f\Doses\AV004f_PVC_Threshold_cropped.am');
    50. 

  50. 

  51. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV002f\Doses\AV002f_Continuous_cropped.am');
    52. 

  52. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV002f\Doses\AV002f_PVC_Continuous_cropped.am');
    53. 

  53. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV002f\Doses\AV002f_Threshold_cropped.am');
    54. 

  54. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV002f\Doses\AV002f_PVC_Threshold_cropped.am');
    55. 

  55. 

  56. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV005c\Doses\AV005c_Continuous_cropped.am');
    57. 

  57. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV005c\Doses\AV005c_PVC_Continuous_cropped.am');
    58. 

  58. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV005c\Doses\AV005c_Threshold_cropped.am');
    59. 

  59. % pet = am2geom('D:\FTP\Dose Painting\AAPM09\AV005c\Doses\AV005c_PVC_Threshold_cropped.am');
    60. 

  60. 

  61. % load('D:\FTP\Dose Painting\AAPM09\HN003c\Doses\HN003c_U01_3MM_continuous.mat');
    62. 

  62. % load('D:\FTP\Dose Painting\AAPM09\HN003c\Doses\HN003c_U01_5MM_continuous.mat');
    63. 

  63. % load('D:\FTP\Dose Painting\AAPM09\HN003c\Doses\HN003c_U01_7MM_continuous.mat');
    64. 

  64. % 
    65. 

  65. % load('D:\FTP\Dose Painting\AAPM09\HN003c\Doses\HN003c_U01_3MM_SUV70_threshold.mat');
    66. 

  66. % load('D:\FTP\Dose Painting\AAPM09\HN003c\Doses\HN003c_U01_5MM_SUV70_threshold.mat');
    67. 

  67. % load('D:\FTP\Dose Painting\AAPM09\HN003c\Doses\HN003c_U01_7MM_SUV70_threshold.mat');
    68. 

  68. % 
    69. 

  69. % load('D:\FTP\Dose Painting\AAPM09\HN003c\Doses\HN003c_U01_3MM_SUV40_CTV54_GTV72.mat');
    70. 

  70. load('D:\FTP\Dose Painting\AAPM09\HN003c\Doses\HN003c_U01_5MM_SUV40_CTV54_GTV72.mat');
    71. 

  71. % load('D:\FTP\Dose Painting\AAPM09\HN003c\Doses\HN003c_U01_7MM_SUV40_CTV54_GTV72.mat');
    72. 

  72. % 
    73. 

  73. % load('D:\FTP\Dose Painting\AAPM09\HN003c\Doses\HN003c_U01_5MM_SUV40_CTV54_GTV94.mat');
    74. 

  74. % load('D:\FTP\Dose Painting\AAPM09\HN003c\Doses\HN003c_U01_5MM_SUV40_CTV54_GTV99.mat');
    75. 

  75. 

  76. petSparse = sparse3D(presc);
    77. 

  77. 

  78. % petSparse = sparse3D(pet.data);
    79. 

  79. petSparse.Nind = length(petSparse.non_zero_indices);
    80. 

  80. 

  81. % write prescription to file
    82. 

  82.  
    83. 

  83. write_presc(petSparse,'GTV');
    84. 

  84. 

  85. fprintf('Uniform prescription sucessfully replaced with non-uniform prescription.\n');
    86.