Pub Date : 2024-11-01Epub Date: 2022-12-10DOI: 10.1016/j.zemedi.2022.11.005
Nicole Fehler, Christian Lingenfelder, Sebastian Kupferschmid, Martin Hessling
Purpose: With diaphanoscopic illumination of the eye, the intensity of light entering its interior depends on the transmission properties of the eyewall. Light that passes through the eyewall can cause damage to the retina. Therefore, in this study, the intraocular irradiances are determined at different positions on the retina, directly behind the illuminated eyewall, the opposite eyewall and near the macula of ex-vivo porcine eyes. These irradiances are examined for their dependence on the pressure applied on the eyewall with the illuminating fiber and for the influence of the pigmentation of the eye.
Methods: In total 221 ex-vivo porcine eyes were investigated. For transscleral illumination an illumination fiber with a diffusing adapter cap is pressed against the equatorial eyewall. The illumination fiber is pressed onto the eye and the pressure is measured in the anterior chamber. Three different pressures are applied, 23, 78 and 132 mmHg. A detection fiber with diffusing fiber tip is inserted into the eye at the desired position. The eyes were divided in groups with high and less pigmentation to investigate the influence of the pigmentation on the intraocular irradiance.
Results: The intraocular irradiances Eintra increases for various increasing applied pressures with the illumination fiber on the eyewall and for various positions inside the eye. With this the irradiances weighted with the photochemical and thermal hazard weighting function, EA-R and EVIR-R, also increases. Differences in Eintra, EA-R and EVIR-R could be found for different pigmented eyes as these values are higher for less pigmented eyes than for strong pigmented ones.
Conclusion: The hazard to the retina during diaphanoscopic illumination of the eye depends on how strong the surgeon presses the illumination fiber on the eyewall. Depending on the applied pressure and the measuring position in the eye, the specified limit for the photochemical hazard to the retina is partly exceeded. The pigmentation of the eye also plays a role. The irradiance in less pigmented eyes appears to be higher than in strongly pigmented eyes. Because of this, the surgeon should be able to adjust the intensity of the light source to the color of the patient's eye.
{"title":"Determination of the intraocular irradiance and potential retinal hazards at various positions in the eye during transscleral equatorial illumination for different applied pressures.","authors":"Nicole Fehler, Christian Lingenfelder, Sebastian Kupferschmid, Martin Hessling","doi":"10.1016/j.zemedi.2022.11.005","DOIUrl":"10.1016/j.zemedi.2022.11.005","url":null,"abstract":"<p><strong>Purpose: </strong>With diaphanoscopic illumination of the eye, the intensity of light entering its interior depends on the transmission properties of the eyewall. Light that passes through the eyewall can cause damage to the retina. Therefore, in this study, the intraocular irradiances are determined at different positions on the retina, directly behind the illuminated eyewall, the opposite eyewall and near the macula of ex-vivo porcine eyes. These irradiances are examined for their dependence on the pressure applied on the eyewall with the illuminating fiber and for the influence of the pigmentation of the eye.</p><p><strong>Methods: </strong>In total 221 ex-vivo porcine eyes were investigated. For transscleral illumination an illumination fiber with a diffusing adapter cap is pressed against the equatorial eyewall. The illumination fiber is pressed onto the eye and the pressure is measured in the anterior chamber. Three different pressures are applied, 23, 78 and 132 mmHg. A detection fiber with diffusing fiber tip is inserted into the eye at the desired position. The eyes were divided in groups with high and less pigmentation to investigate the influence of the pigmentation on the intraocular irradiance.</p><p><strong>Results: </strong>The intraocular irradiances E<sub>intra</sub> increases for various increasing applied pressures with the illumination fiber on the eyewall and for various positions inside the eye. With this the irradiances weighted with the photochemical and thermal hazard weighting function, E<sub>A-R</sub> and E<sub>VIR-R</sub>, also increases. Differences in E<sub>intra</sub>, E<sub>A-R</sub> and E<sub>VIR-R</sub> could be found for different pigmented eyes as these values are higher for less pigmented eyes than for strong pigmented ones.</p><p><strong>Conclusion: </strong>The hazard to the retina during diaphanoscopic illumination of the eye depends on how strong the surgeon presses the illumination fiber on the eyewall. Depending on the applied pressure and the measuring position in the eye, the specified limit for the photochemical hazard to the retina is partly exceeded. The pigmentation of the eye also plays a role. The irradiance in less pigmented eyes appears to be higher than in strongly pigmented eyes. Because of this, the surgeon should be able to adjust the intensity of the light source to the color of the patient's eye.</p>","PeriodicalId":54397,"journal":{"name":"Zeitschrift fur Medizinische Physik","volume":" ","pages":"610-619"},"PeriodicalIF":2.4,"publicationDate":"2024-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"10338500","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-11-01Epub Date: 2022-12-19DOI: 10.1016/j.zemedi.2022.11.007
Achim Langenbucher, Peter Hoffmann, Alan Cayless, Damien Gatinel, Guillaume Debellemanière, Jascha Wendelstein, Nóra Szentmáry
Background: In recent years, some lens manufacturers have committed to providing lens shape data for some of their lens models. The purpose of this study is to present a strategy for prediction of intraocular lens power and residual refraction based on a pseudophakic model eye containing 5 refractive surfaces and to show its applicability using worked examples.
Methods: A pseudophakic model eye with a thin spectacle correction, a thick cornea (radius of curvatures for both surfaces and central thickness) and a thick IOL (either radius of curvatures RLa and RLp for front and back surface or equivalent power PL and Coddington factor CL; and either central thickness LT or edge thickness and optic diameter) was set up. Calculations were performed based on linear Gaussian optics (vergence formulae). Formulae were provided to derive the lens power/shape and the residual equivalent spectacle refraction SEQ. From the lens shape the location of the haptic plane HP, the image sided principal plane of the lens HL, and the ocular magnification OM were extracted.
Results: The calculation of a thick intraocular lens and the prediction of residual refraction is presented with reference to 3 working examples: A) lens varied in PL and shifted with its haptic plane keeping the CL constant, B) lens varied in CL and shifted with its haptic plane keeping PL constant, and C) CL and PL of the lens varied keeping its haptic plane position in the eye constant. For each combination of parameters (PL, CL, or haptic plane shift) the parameters influencing SEQ, OM and HL-HP were analysed.
Conclusion: Some modern optical biometers currently on the market provide the radii of curvature of both corneal surface and all relevant distances in the eye. With additional data on the lens shape, it would be possible to improve lens power calculations by switching from thin to thick lens models for the cornea and for the lens. This would overcome one of the major drawbacks of current lens power calculation methods.
{"title":"Considerations of a thick lens formula for intraocular lens power calculation.","authors":"Achim Langenbucher, Peter Hoffmann, Alan Cayless, Damien Gatinel, Guillaume Debellemanière, Jascha Wendelstein, Nóra Szentmáry","doi":"10.1016/j.zemedi.2022.11.007","DOIUrl":"10.1016/j.zemedi.2022.11.007","url":null,"abstract":"<p><strong>Background: </strong>In recent years, some lens manufacturers have committed to providing lens shape data for some of their lens models. The purpose of this study is to present a strategy for prediction of intraocular lens power and residual refraction based on a pseudophakic model eye containing 5 refractive surfaces and to show its applicability using worked examples.</p><p><strong>Methods: </strong>A pseudophakic model eye with a thin spectacle correction, a thick cornea (radius of curvatures for both surfaces and central thickness) and a thick IOL (either radius of curvatures RLa and RLp for front and back surface or equivalent power PL and Coddington factor CL; and either central thickness LT or edge thickness and optic diameter) was set up. Calculations were performed based on linear Gaussian optics (vergence formulae). Formulae were provided to derive the lens power/shape and the residual equivalent spectacle refraction SEQ. From the lens shape the location of the haptic plane HP, the image sided principal plane of the lens HL, and the ocular magnification OM were extracted.</p><p><strong>Results: </strong>The calculation of a thick intraocular lens and the prediction of residual refraction is presented with reference to 3 working examples: A) lens varied in PL and shifted with its haptic plane keeping the CL constant, B) lens varied in CL and shifted with its haptic plane keeping PL constant, and C) CL and PL of the lens varied keeping its haptic plane position in the eye constant. For each combination of parameters (PL, CL, or haptic plane shift) the parameters influencing SEQ, OM and HL-HP were analysed.</p><p><strong>Conclusion: </strong>Some modern optical biometers currently on the market provide the radii of curvature of both corneal surface and all relevant distances in the eye. With additional data on the lens shape, it would be possible to improve lens power calculations by switching from thin to thick lens models for the cornea and for the lens. This would overcome one of the major drawbacks of current lens power calculation methods.</p>","PeriodicalId":54397,"journal":{"name":"Zeitschrift fur Medizinische Physik","volume":" ","pages":"620-631"},"PeriodicalIF":2.4,"publicationDate":"2024-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"10764181","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-11-01Epub Date: 2023-04-06DOI: 10.1016/j.zemedi.2023.02.006
Clarita Saldarriaga Vargas, Michelle Andersson, Céline Bouvier-Capely, Wei Bo Li, Balázs Madas, Peter Covens, Lara Struelens, Lidia Strigari
Absorbed dose heterogeneity in kidney tissues is an important issue in radiopharmaceutical therapy. The effect of absorbed dose heterogeneity in nephrotoxicity is, however, not fully understood yet, which hampers the implementation of treatment optimization by obscuring the interpretation of clinical response data and the selection of optimal treatment options. Although some dosimetry methods have been developed for kidney dosimetry to the level of microscopic renal substructures, the clinical assessment of the microscopic distribution of radiopharmaceuticals in kidney tissues currently remains a challenge. This restricts the anatomical resolution of clinical dosimetry, which hinders a thorough clinical investigation of the impact of absorbed dose heterogeneity. The potential of absorbed dose-response modelling to support individual treatment optimization in radiopharmaceutical therapy is recognized and gaining attraction. However, biophysical modelling is currently underexplored for the kidney, where particular modelling challenges arise from the convolution of a complex functional organization of renal tissues with the function-mediated dose distribution of radiopharmaceuticals. This article reviews and discusses the heterogeneity of absorbed dose distribution in kidney tissues and the absorbed dose-response modelling of nephrotoxicity in radiopharmaceutical therapy. The review focuses mainly on the peptide receptor radionuclide therapy with beta-particle emitting somatostatin analogues, for which the scientific literature reflects over two decades of clinical experience. Additionally, detailed research perspectives are proposed to address various identified challenges to progress in this field.
{"title":"Heterogeneity of absorbed dose distribution in kidney tissues and dose-response modelling of nephrotoxicity in radiopharmaceutical therapy with beta-particle emitters: A review.","authors":"Clarita Saldarriaga Vargas, Michelle Andersson, Céline Bouvier-Capely, Wei Bo Li, Balázs Madas, Peter Covens, Lara Struelens, Lidia Strigari","doi":"10.1016/j.zemedi.2023.02.006","DOIUrl":"10.1016/j.zemedi.2023.02.006","url":null,"abstract":"<p><p>Absorbed dose heterogeneity in kidney tissues is an important issue in radiopharmaceutical therapy. The effect of absorbed dose heterogeneity in nephrotoxicity is, however, not fully understood yet, which hampers the implementation of treatment optimization by obscuring the interpretation of clinical response data and the selection of optimal treatment options. Although some dosimetry methods have been developed for kidney dosimetry to the level of microscopic renal substructures, the clinical assessment of the microscopic distribution of radiopharmaceuticals in kidney tissues currently remains a challenge. This restricts the anatomical resolution of clinical dosimetry, which hinders a thorough clinical investigation of the impact of absorbed dose heterogeneity. The potential of absorbed dose-response modelling to support individual treatment optimization in radiopharmaceutical therapy is recognized and gaining attraction. However, biophysical modelling is currently underexplored for the kidney, where particular modelling challenges arise from the convolution of a complex functional organization of renal tissues with the function-mediated dose distribution of radiopharmaceuticals. This article reviews and discusses the heterogeneity of absorbed dose distribution in kidney tissues and the absorbed dose-response modelling of nephrotoxicity in radiopharmaceutical therapy. The review focuses mainly on the peptide receptor radionuclide therapy with beta-particle emitting somatostatin analogues, for which the scientific literature reflects over two decades of clinical experience. Additionally, detailed research perspectives are proposed to address various identified challenges to progress in this field.</p>","PeriodicalId":54397,"journal":{"name":"Zeitschrift fur Medizinische Physik","volume":" ","pages":"491-509"},"PeriodicalIF":2.4,"publicationDate":"2024-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"9629564","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-11-01Epub Date: 2023-04-21DOI: 10.1016/j.zemedi.2023.03.003
Adam Ryczkowski, Tomasz Piotrowski, Marcin Staszczak, Marcin Wiktorowicz, Przemysław Adrich
Introduction: The shape of the energy spectrum is an essential component of any electron beam Monte Carlo model. Due to specialized equipment and the long measurement time for the direct methods for determining the energy spectrum, attractive alternatives are backward spectrum reconstructions from the measured data. One such approach is solving the first-degree Fredholm integral equation with appropriate regularization. It makes it possible to calculate the depth distribution as the sum of the distributions from monoenergetic beams. This study aims to determine the optimal value of the regularization parameter for the problem of determining the spectrum of the electron beam produced by a mobile accelerator used during intraoperative radiotherapy.
Material and methods: The Geant4 package was used to generate the distributions of deep doses for monoenergetic beams for two models with different degrees of complexity, i.e. simple (theoretical) and full (for the mobile accelerator). The dose distributions for four different shapes of energy spectrum (for each model) were obtained similarly. They were established as the reference data for further calculations. The Dual Annealing optimization method was used to obtain the reconstructed spectrum. The multiple optimizations that differ by the regularization parameter (ranging from 0 to 1) were performed. For each reconstruction, similarity indicators of the energy spectrum and the dose distribution to the referenced data were calculated to determine the optimal regularization parameters.
Results: Optimal regularization parameters determined by similarity indicators for the spectrum and the dose distribution differ for geometry models considered in the study. The regularization parameter for the simple geometry ranged from 0.03 to 0.05, while for full geometry, they were from 0.05 to 0.06. The results for conventional linear accelerators found in the literature range from 0.5 to 1.1.
Conclusion: The Dual Annealing optimization method can be effectively used to solve the Fredholm equation with Tikhonov regularization to reconstruct an electron beam's energy spectrum. The regularization parameter value depends on the beam-forming system. Its value for the mobile accelerator considered in the study ranges from 0.05 to 0.06, depending on the nominal beam energy value.
{"title":"Optimization of the regularization parameter in the Dual Annealing method used for the reconstruction of energy spectrum of electron beam generated by the AQURE mobile accelerator.","authors":"Adam Ryczkowski, Tomasz Piotrowski, Marcin Staszczak, Marcin Wiktorowicz, Przemysław Adrich","doi":"10.1016/j.zemedi.2023.03.003","DOIUrl":"10.1016/j.zemedi.2023.03.003","url":null,"abstract":"<p><strong>Introduction: </strong>The shape of the energy spectrum is an essential component of any electron beam Monte Carlo model. Due to specialized equipment and the long measurement time for the direct methods for determining the energy spectrum, attractive alternatives are backward spectrum reconstructions from the measured data. One such approach is solving the first-degree Fredholm integral equation with appropriate regularization. It makes it possible to calculate the depth distribution as the sum of the distributions from monoenergetic beams. This study aims to determine the optimal value of the regularization parameter for the problem of determining the spectrum of the electron beam produced by a mobile accelerator used during intraoperative radiotherapy.</p><p><strong>Material and methods: </strong>The Geant4 package was used to generate the distributions of deep doses for monoenergetic beams for two models with different degrees of complexity, i.e. simple (theoretical) and full (for the mobile accelerator). The dose distributions for four different shapes of energy spectrum (for each model) were obtained similarly. They were established as the reference data for further calculations. The Dual Annealing optimization method was used to obtain the reconstructed spectrum. The multiple optimizations that differ by the regularization parameter (ranging from 0 to 1) were performed. For each reconstruction, similarity indicators of the energy spectrum and the dose distribution to the referenced data were calculated to determine the optimal regularization parameters.</p><p><strong>Results: </strong>Optimal regularization parameters determined by similarity indicators for the spectrum and the dose distribution differ for geometry models considered in the study. The regularization parameter for the simple geometry ranged from 0.03 to 0.05, while for full geometry, they were from 0.05 to 0.06. The results for conventional linear accelerators found in the literature range from 0.5 to 1.1.</p><p><strong>Conclusion: </strong>The Dual Annealing optimization method can be effectively used to solve the Fredholm equation with Tikhonov regularization to reconstruct an electron beam's energy spectrum. The regularization parameter value depends on the beam-forming system. Its value for the mobile accelerator considered in the study ranges from 0.05 to 0.06, depending on the nominal beam energy value.</p>","PeriodicalId":54397,"journal":{"name":"Zeitschrift fur Medizinische Physik","volume":" ","pages":"510-520"},"PeriodicalIF":2.4,"publicationDate":"2024-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"9775836","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-11-01Epub Date: 2023-02-20DOI: 10.1016/j.zemedi.2022.11.009
Achim Langenbucher, Nóra Szentmáry, Alan Cayless, Jascha Wendelstein, Peter Hoffmann
Purpose: To implement a fully data driven strategy for identifying outliers in clinical datasets used for formula constant optimisation, in order to achieve proper formula predicted refraction after cataract surgery, and to assess the capabilities of this outlier detection method.
Methods: 2 clinical datasets (DS1/DS2: N = 888/403) of eyes treated with a monofocal aspherical intraocular lens (Hoya XY1/Johnson&Johnson Vision Z9003) containing preoperative biometric data, power of the lens implant and postoperative spherical equivalent (SEQ) were transferred to us for formula constant optimisation. Original datasets were used to generate baseline formula constants. A random forest quantile regression algorithm was set up using bootstrap resampling with replacement. Quantile regression trees were grown and the 25% and 75% quantile, and the interquartile range were extracted from SEQ and formula predicted refraction REF for the SRKT, Haigis and Castrop formulae. Fences were defined from the quantiles and data points outside the fences were marked and removed as outliers before recalculating the formula constants.
Results: NB = 1000 bootstrap samples were derived from both datasets, and random forest quantile regression trees were grown to model SEQ versus REF and to estimate the median and 25% and 75% quantiles. The fence boundaries were defined as being from 25% quantile - 1.5·IQR to 75% quantile + 1.5·IQR, with data points outside the fence being marked as outliers. In total, for DS1 and DS2, 25/27/32 and 4/5/4 data points were identified as outliers for the SRKT/Haigis/Castrop formulae respectively. The respective root mean squared formula prediction errors for the three formulae were slightly reduced from: 0.4370 dpt;0.4449 dpt/0.3625 dpt;0.4056 dpt/and 0.3376 dpt;0.3532 dpt to: 0.4271 dpt;0.4348 dpt/0.3528 dpt;0.3952 dpt/0.3277 dpt;0.3432 dpt for DS1;DS2.
Conclusion: We were able to prove that with random forest quantile regression trees a fully data driven outlier identification strategy acting in the response space is achievable. In a real life scenario this strategy has to be complemented by an outlier identification method acting in the parameter space for a proper qualification of datasets prior to formula constant optimisation.
{"title":"Preconditioning of clinical data for intraocular lens formula constant optimisation using Random Forest Quantile Regression Trees.","authors":"Achim Langenbucher, Nóra Szentmáry, Alan Cayless, Jascha Wendelstein, Peter Hoffmann","doi":"10.1016/j.zemedi.2022.11.009","DOIUrl":"10.1016/j.zemedi.2022.11.009","url":null,"abstract":"<p><strong>Purpose: </strong>To implement a fully data driven strategy for identifying outliers in clinical datasets used for formula constant optimisation, in order to achieve proper formula predicted refraction after cataract surgery, and to assess the capabilities of this outlier detection method.</p><p><strong>Methods: </strong>2 clinical datasets (DS1/DS2: N = 888/403) of eyes treated with a monofocal aspherical intraocular lens (Hoya XY1/Johnson&Johnson Vision Z9003) containing preoperative biometric data, power of the lens implant and postoperative spherical equivalent (SEQ) were transferred to us for formula constant optimisation. Original datasets were used to generate baseline formula constants. A random forest quantile regression algorithm was set up using bootstrap resampling with replacement. Quantile regression trees were grown and the 25% and 75% quantile, and the interquartile range were extracted from SEQ and formula predicted refraction REF for the SRKT, Haigis and Castrop formulae. Fences were defined from the quantiles and data points outside the fences were marked and removed as outliers before recalculating the formula constants.</p><p><strong>Results: </strong>N<sub>B</sub> = 1000 bootstrap samples were derived from both datasets, and random forest quantile regression trees were grown to model SEQ versus REF and to estimate the median and 25% and 75% quantiles. The fence boundaries were defined as being from 25% quantile - 1.5·IQR to 75% quantile + 1.5·IQR, with data points outside the fence being marked as outliers. In total, for DS1 and DS2, 25/27/32 and 4/5/4 data points were identified as outliers for the SRKT/Haigis/Castrop formulae respectively. The respective root mean squared formula prediction errors for the three formulae were slightly reduced from: 0.4370 dpt;0.4449 dpt/0.3625 dpt;0.4056 dpt/and 0.3376 dpt;0.3532 dpt to: 0.4271 dpt;0.4348 dpt/0.3528 dpt;0.3952 dpt/0.3277 dpt;0.3432 dpt for DS1;DS2.</p><p><strong>Conclusion: </strong>We were able to prove that with random forest quantile regression trees a fully data driven outlier identification strategy acting in the response space is achievable. In a real life scenario this strategy has to be complemented by an outlier identification method acting in the parameter space for a proper qualification of datasets prior to formula constant optimisation.</p>","PeriodicalId":54397,"journal":{"name":"Zeitschrift fur Medizinische Physik","volume":" ","pages":"632-640"},"PeriodicalIF":2.4,"publicationDate":"2024-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"10816824","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-11-01Epub Date: 2023-04-22DOI: 10.1016/j.zemedi.2023.03.002
Angelika Altergot, Michaela Schürmann, Tanja Jungert, Hendrik Auerbach, Frank Nüsken, Jan Palm, Christian Rübe, Claudia E Rübe, Yvonne Dzierma
Introduction: Image guided radiotherapy allows for particularly conformal tumour irradiation through precise patient positioning. Becoming the standard for radiotherapy, this increases imaging doses to the patient. The Halcyon 3.0 linear accelerator (Varian Medical Systems, Palo Alto, CA) requires daily imaging due to its geometry. For this reason, the accelerator is equipped with on-line kV and MV imaging. However, daily CBCT images required for irradiation apply additional radiation, which increases the dose to normal tissue and therefore can affect the patient's secondary cancer risk. In this study, actual organ doses were measured for the kV system, and a comparison of normal tissue doses for all available kV CBCT protocols was presented to demonstrate differences in imaging doses across entities and protocols. In addition, effective dose and secondary cancer risk from imaging are evaluated.
Material and methods: Measurements were performed with thermoluminescent dosimeters in an anthropomorphic phantom positioned according to each entity (brain, head and neck, breast, lung, pelvis). CBCT images were obtained, using all available pre-set protocols without further adjustment of the parameters. Measured doses for each position and each protocol were then compared and secondary cancer risk of relevant and specifically radiosensitive organs was calculated.
Results: It was found that imaging doses for protocols such as Pelvis and Head could be reduced by up to half using the corresponding Fast and Low Dose modes, respectively. On the other hand, larger field sizes or the Large mode yielded higher doses than their initial protocols. Image Gently was found to spare normal tissue best, however it is not suitable for certain entities due to low image quality or insufficient projection data.
Discussion: By using appropriate kV-CBCT protocols, it is possible to reduce imaging doses to a significant extent and therefore spare healthy tissue. Combined with studies of image quality, the results of this study could lead to adjustments in workflow regarding the choice of protocols used in daily routine. This could prevent unnecessary radiation exposure and reduce secondary cancer risk.
{"title":"Imaging doses for different CBCT protocols on the Halcyon 3.0 linear accelerator - TLD measurements in an anthropomorphic phantom.","authors":"Angelika Altergot, Michaela Schürmann, Tanja Jungert, Hendrik Auerbach, Frank Nüsken, Jan Palm, Christian Rübe, Claudia E Rübe, Yvonne Dzierma","doi":"10.1016/j.zemedi.2023.03.002","DOIUrl":"10.1016/j.zemedi.2023.03.002","url":null,"abstract":"<p><strong>Introduction: </strong>Image guided radiotherapy allows for particularly conformal tumour irradiation through precise patient positioning. Becoming the standard for radiotherapy, this increases imaging doses to the patient. The Halcyon 3.0 linear accelerator (Varian Medical Systems, Palo Alto, CA) requires daily imaging due to its geometry. For this reason, the accelerator is equipped with on-line kV and MV imaging. However, daily CBCT images required for irradiation apply additional radiation, which increases the dose to normal tissue and therefore can affect the patient's secondary cancer risk. In this study, actual organ doses were measured for the kV system, and a comparison of normal tissue doses for all available kV CBCT protocols was presented to demonstrate differences in imaging doses across entities and protocols. In addition, effective dose and secondary cancer risk from imaging are evaluated.</p><p><strong>Material and methods: </strong>Measurements were performed with thermoluminescent dosimeters in an anthropomorphic phantom positioned according to each entity (brain, head and neck, breast, lung, pelvis). CBCT images were obtained, using all available pre-set protocols without further adjustment of the parameters. Measured doses for each position and each protocol were then compared and secondary cancer risk of relevant and specifically radiosensitive organs was calculated.</p><p><strong>Results: </strong>It was found that imaging doses for protocols such as Pelvis and Head could be reduced by up to half using the corresponding Fast and Low Dose modes, respectively. On the other hand, larger field sizes or the Large mode yielded higher doses than their initial protocols. Image Gently was found to spare normal tissue best, however it is not suitable for certain entities due to low image quality or insufficient projection data.</p><p><strong>Discussion: </strong>By using appropriate kV-CBCT protocols, it is possible to reduce imaging doses to a significant extent and therefore spare healthy tissue. Combined with studies of image quality, the results of this study could lead to adjustments in workflow regarding the choice of protocols used in daily routine. This could prevent unnecessary radiation exposure and reduce secondary cancer risk.</p>","PeriodicalId":54397,"journal":{"name":"Zeitschrift fur Medizinische Physik","volume":" ","pages":"580-595"},"PeriodicalIF":2.4,"publicationDate":"2024-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"9775868","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-11-01Epub Date: 2023-03-25DOI: 10.1016/j.zemedi.2023.02.001
Sepideh Hatamikia, Ander Biguri, Gernot Kronreif, Tom Russ, Joachim Kettenbach, Wolfgang Birkfellner
Precise instrument placement plays a critical role in all interventional procedures, especially percutaneous procedures such as needle biopsies, to achieve successful tumor targeting and increased diagnostic accuracy. C-arm cone beam computed tomography (CBCT) has the potential to precisely visualize the anatomy in direct vicinity of the needle and evaluate the adequacy of needle placement during the intervention, allowing for instantaneous adjustment in case of misplacement. However, even with the most advanced C-arm CBCT devices, it can be difficult to identify the exact needle position on CBCT images due to the strong metal artifacts around the needle. In this study, we proposed a framework for customized trajectory design in CBCT imaging based on Prior Image Constrained Compressed Sensing (PICCS) reconstruction with the goal of reducing metal artifacts in needle-based procedures. We proposed to optimize out-of-plane rotations in three-dimensional (3D) space and minimize projection views while reducing metal artifacts at specific volume of interests (VOIs). An anthropomorphic thorax phantom with a needle inserted inside and two tumor models as the imaging targets were used to validate the proposed approach. The performance of the proposed approach was also evaluated for CBCT imaging under kinematic constraints by simulating some collision areas on the geometry of the C-arm. We compared the result of optimized 3D trajectories using the PICCS algorithm and 20 projections with the result of a circular trajectory with sparse view using PICCS and Feldkamp, Davis, and Kress (FDK), both using 20 projections, and the circular FDK method with 313 projections. For imaging targets 1 and 2, the highest values of structural similarity index measure (SSIM) and universal quality index (UQI) between the reconstructed image from the optimized trajectories and the initial CBCT image at the VOI was calculated 0.7521, 0.7308 and 0.7308, 0.7248 respectively. These results significantly outperformed the FDK method (with 20 and 313 projections) and the PICCS method (20 projections) both using the circular trajectory. Our results showed that the proposed optimized trajectories not only significantly reduce metal artifacts but also suggest a dose reduction for needle-based CBCT interventions, considering the small number of projections used. Furthermore, our results showed that the optimized trajectories are compatible with spatially constrained situations and enable CBCT imaging under kinematic constraints when the standard circular trajectory is not feasible.
{"title":"Source-detector trajectory optimization for CBCT metal artifact reduction based on PICCS reconstruction.","authors":"Sepideh Hatamikia, Ander Biguri, Gernot Kronreif, Tom Russ, Joachim Kettenbach, Wolfgang Birkfellner","doi":"10.1016/j.zemedi.2023.02.001","DOIUrl":"10.1016/j.zemedi.2023.02.001","url":null,"abstract":"<p><p>Precise instrument placement plays a critical role in all interventional procedures, especially percutaneous procedures such as needle biopsies, to achieve successful tumor targeting and increased diagnostic accuracy. C-arm cone beam computed tomography (CBCT) has the potential to precisely visualize the anatomy in direct vicinity of the needle and evaluate the adequacy of needle placement during the intervention, allowing for instantaneous adjustment in case of misplacement. However, even with the most advanced C-arm CBCT devices, it can be difficult to identify the exact needle position on CBCT images due to the strong metal artifacts around the needle. In this study, we proposed a framework for customized trajectory design in CBCT imaging based on Prior Image Constrained Compressed Sensing (PICCS) reconstruction with the goal of reducing metal artifacts in needle-based procedures. We proposed to optimize out-of-plane rotations in three-dimensional (3D) space and minimize projection views while reducing metal artifacts at specific volume of interests (VOIs). An anthropomorphic thorax phantom with a needle inserted inside and two tumor models as the imaging targets were used to validate the proposed approach. The performance of the proposed approach was also evaluated for CBCT imaging under kinematic constraints by simulating some collision areas on the geometry of the C-arm. We compared the result of optimized 3D trajectories using the PICCS algorithm and 20 projections with the result of a circular trajectory with sparse view using PICCS and Feldkamp, Davis, and Kress (FDK), both using 20 projections, and the circular FDK method with 313 projections. For imaging targets 1 and 2, the highest values of structural similarity index measure (SSIM) and universal quality index (UQI) between the reconstructed image from the optimized trajectories and the initial CBCT image at the VOI was calculated 0.7521, 0.7308 and 0.7308, 0.7248 respectively. These results significantly outperformed the FDK method (with 20 and 313 projections) and the PICCS method (20 projections) both using the circular trajectory. Our results showed that the proposed optimized trajectories not only significantly reduce metal artifacts but also suggest a dose reduction for needle-based CBCT interventions, considering the small number of projections used. Furthermore, our results showed that the optimized trajectories are compatible with spatially constrained situations and enable CBCT imaging under kinematic constraints when the standard circular trajectory is not feasible.</p>","PeriodicalId":54397,"journal":{"name":"Zeitschrift fur Medizinische Physik","volume":" ","pages":"565-579"},"PeriodicalIF":2.4,"publicationDate":"2024-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"9199379","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-11-01Epub Date: 2023-01-28DOI: 10.1016/j.zemedi.2023.01.002
Simon Reiss, Kevin Wäscher, Ali Caglar Özen, Thomas Lottner, Timo Heidt, Constantin von Zur Mühlen, Michael Bock
Purpose: To test intra-arterial spin labeling (iASL) using active guiding catheters for myocardial perfusion measurements during magnetic resonance (MR)-guided interventions in a pig study.
Methods: In this work, a single-loop radiofrequency (RF) coil at the tip of a 6F active coronary catheter was used as a transmit coil for local spin labeling. The transmit magnetic RF field (B1) of the coil and the labeling efficiency were determined, and iASL was tested in two pigs after the catheter was engaged in the aortic root, the ostium of the left coronary artery (LCA) under MR-guidance. The iASL effect was assessed by the signal difference between spin-labeling On and control (spin-labeling OFF) images, and in a cross-correlation between ON/Off states of spin-labeling a binary labeling paradigm. In addition, quantitative myocardial perfusion was calculated from the iASL experiments.
Results: The maximum B1 in the vicinity of the catheter coil was 2.1 µT. A strong local labeling effect with a labeling efficiency of 0.45 was achieved with iASL both in vitro and in vivo. In both pigs, the proximal myocardial segments supplied by the LCA showed significant labelling effect up to distances of 60 mm from the aortic root with a relative signal difference of (3.14 ± 2.89)% in the first and (3.50 ± 1.25)% in the second animal. The mean correlation coefficients were R = 0.63 ± 0.22 and 0.42 ± 0.16, respectively. The corresponding computed myocardial perfusion values in this region of the myocardium were similar to those obtained with contrast perfusion methods ((1.2 ± 1.1) mL/min/g and (0.8 ± 0.6) mL/min/g).
Conclusion: The proposed iASL method demonstrates the feasibility of selective myocardial perfusion measurements during MR-guided coronary interventions, which with further technical improvements may provide an alternative to exogenous contrast-based perfusion. Due to the invasive nature of the iASL method, it can potentially be used in concert with MRI-guided coronary angioplasty.
{"title":"Quantifying myocardial perfusion during MR-guided interventions without exogenous contrast agents: intra-arterial spin labeling.","authors":"Simon Reiss, Kevin Wäscher, Ali Caglar Özen, Thomas Lottner, Timo Heidt, Constantin von Zur Mühlen, Michael Bock","doi":"10.1016/j.zemedi.2023.01.002","DOIUrl":"10.1016/j.zemedi.2023.01.002","url":null,"abstract":"<p><strong>Purpose: </strong>To test intra-arterial spin labeling (iASL) using active guiding catheters for myocardial perfusion measurements during magnetic resonance (MR)-guided interventions in a pig study.</p><p><strong>Methods: </strong>In this work, a single-loop radiofrequency (RF) coil at the tip of a 6F active coronary catheter was used as a transmit coil for local spin labeling. The transmit magnetic RF field (B<sub>1</sub>) of the coil and the labeling efficiency were determined, and iASL was tested in two pigs after the catheter was engaged in the aortic root, the ostium of the left coronary artery (LCA) under MR-guidance. The iASL effect was assessed by the signal difference between spin-labeling On and control (spin-labeling OFF) images, and in a cross-correlation between ON/Off states of spin-labeling a binary labeling paradigm. In addition, quantitative myocardial perfusion was calculated from the iASL experiments.</p><p><strong>Results: </strong>The maximum B<sub>1</sub> in the vicinity of the catheter coil was 2.1 µT. A strong local labeling effect with a labeling efficiency of 0.45 was achieved with iASL both in vitro and in vivo. In both pigs, the proximal myocardial segments supplied by the LCA showed significant labelling effect up to distances of 60 mm from the aortic root with a relative signal difference of (3.14 ± 2.89)% in the first and (3.50 ± 1.25)% in the second animal. The mean correlation coefficients were R = 0.63 ± 0.22 and 0.42 ± 0.16, respectively. The corresponding computed myocardial perfusion values in this region of the myocardium were similar to those obtained with contrast perfusion methods ((1.2 ± 1.1) mL/min/g and (0.8 ± 0.6) mL/min/g).</p><p><strong>Conclusion: </strong>The proposed iASL method demonstrates the feasibility of selective myocardial perfusion measurements during MR-guided coronary interventions, which with further technical improvements may provide an alternative to exogenous contrast-based perfusion. Due to the invasive nature of the iASL method, it can potentially be used in concert with MRI-guided coronary angioplasty.</p>","PeriodicalId":54397,"journal":{"name":"Zeitschrift fur Medizinische Physik","volume":" ","pages":"596-609"},"PeriodicalIF":2.4,"publicationDate":"2024-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"10587236","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-11-01Epub Date: 2023-03-29DOI: 10.1016/j.zemedi.2023.02.003
Klara Uher, Stefanie Ehrbar, Stephanie Tanadini-Lang, Riccardo Dal Bello
Purpose: Volumetric modulated arc therapy (VMAT) is a widespread technique for the delivery of normo-fractionated radiation therapy (NFRT) and stereotactic body radiation therapy (SBRT). It is associated with a significant hardware burden requiring dose rate modulation, collimator movement and gantry rotation synchronisation. Patient specific quality assurance (PSQA) guarantees that the linacs can precisely and accurately deliver the planned dose. However, PSQA requires a significant time allocation and class solutions to reduce this while guaranteeing the deliverability of the plans should be investigated.
Methods: In this study, an in-house developed Eclipse Scripting API (ESAPI) script was used to extract five independent plan complexity metrics from N = 667 VMAT treatment fields. The correlation between metrics and portal dosimetry measurements was investigated with Pearson correlation, box plot analysis and receiver operating characteristic curves, which were used to defined the best performing metric and its threshold.
Results: The incidence of fields failing the clinical PSQA criteria of 3%/2mm (NFRT) and 3%/1.5mm (SBRT) was low (N = 1). The mean MLC opening was the metric with the highest correlation with the portal dosimetry data and among the best in discriminating the requirement of PSQA. The thresholds of 16.12 mm (NFRT) and 7.96 mm (SBRT) corresponded to true positive rates higher than 90%.
Conclusions: This work presents a quantitative approach to reduce the time allocation for PSQA by identifying the most complex plans demanding a dedicated measurement. The proposed method requires PSQA for approximately 10% of the plans. The ESAPI script is distributed open-source to ease the investigation and implementation at other institutions.
{"title":"Reduction of patient specific quality assurance through plan complexity metrics for VMAT plans with an open-source TPS script.","authors":"Klara Uher, Stefanie Ehrbar, Stephanie Tanadini-Lang, Riccardo Dal Bello","doi":"10.1016/j.zemedi.2023.02.003","DOIUrl":"10.1016/j.zemedi.2023.02.003","url":null,"abstract":"<p><strong>Purpose: </strong>Volumetric modulated arc therapy (VMAT) is a widespread technique for the delivery of normo-fractionated radiation therapy (NFRT) and stereotactic body radiation therapy (SBRT). It is associated with a significant hardware burden requiring dose rate modulation, collimator movement and gantry rotation synchronisation. Patient specific quality assurance (PSQA) guarantees that the linacs can precisely and accurately deliver the planned dose. However, PSQA requires a significant time allocation and class solutions to reduce this while guaranteeing the deliverability of the plans should be investigated.</p><p><strong>Methods: </strong>In this study, an in-house developed Eclipse Scripting API (ESAPI) script was used to extract five independent plan complexity metrics from N = 667 VMAT treatment fields. The correlation between metrics and portal dosimetry measurements was investigated with Pearson correlation, box plot analysis and receiver operating characteristic curves, which were used to defined the best performing metric and its threshold.</p><p><strong>Results: </strong>The incidence of fields failing the clinical PSQA criteria of 3%/2mm (NFRT) and 3%/1.5mm (SBRT) was low (N = 1). The mean MLC opening was the metric with the highest correlation with the portal dosimetry data and among the best in discriminating the requirement of PSQA. The thresholds of 16.12 mm (NFRT) and 7.96 mm (SBRT) corresponded to true positive rates higher than 90%.</p><p><strong>Conclusions: </strong>This work presents a quantitative approach to reduce the time allocation for PSQA by identifying the most complex plans demanding a dedicated measurement. The proposed method requires PSQA for approximately 10% of the plans. The ESAPI script is distributed open-source to ease the investigation and implementation at other institutions.</p>","PeriodicalId":54397,"journal":{"name":"Zeitschrift fur Medizinische Physik","volume":" ","pages":"555-564"},"PeriodicalIF":2.4,"publicationDate":"2024-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"9279104","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-11-01Epub Date: 2023-05-22DOI: 10.1016/j.zemedi.2023.03.005
Severin Kampfer, Sophie Dobiasch, Stephanie E Combs, Jan J Wilkens
In human radiotherapy a safety margin (PTV margin) is essential for successful irradiation and is usually part of clinical treatment planning. In preclinical radiotherapy research with small animals, most uncertainties and inaccuracies are present as well, but according to the literature a margin is used only scarcely. In addition, there is only little experience about the appropriate size of the margin, which should carefully be investigated and considered, since sparing of organs at risk or normal tissue is affected. Here we estimate the needed margin for preclinical irradiation by adapting a well-known human margin recipe from van Herck et al. to the dimensions and requirements of the specimen on a small animal radiation research platform (SARRP). We adjusted the factors of the described formula to the specific challenges in an orthotopic pancreatic tumor mouse model to establish an appropriate margin concept. The SARRP was used with its image-guidance irradiation possibility for arc irradiation with a field size of 10 × 10 mm2 for 5 fractions. Our goal was to irradiate the clinical target volume (CTV) of at least 90% of our mice with at least 95% of the prescribed dose. By carefully analyzing all relevant factors we gain a CTV to planning target volume (PTV) margin of 1.5 mm for our preclinical setup. The stated safety margin is strongly dependent on the exact setting of the experiment and has to be adjusted for other experimental settings. The few stated values in literature correspond well to our result. Even if using margins in the preclinical setting might be an additional challenge, we think it is crucial to use them to produce reliable results and improve the efficacy of radiotherapy.
{"title":"Development of a PTV margin for preclinical irradiation of orthotopic pancreatic tumors derived from a well-known recipe for humans.","authors":"Severin Kampfer, Sophie Dobiasch, Stephanie E Combs, Jan J Wilkens","doi":"10.1016/j.zemedi.2023.03.005","DOIUrl":"10.1016/j.zemedi.2023.03.005","url":null,"abstract":"<p><p>In human radiotherapy a safety margin (PTV margin) is essential for successful irradiation and is usually part of clinical treatment planning. In preclinical radiotherapy research with small animals, most uncertainties and inaccuracies are present as well, but according to the literature a margin is used only scarcely. In addition, there is only little experience about the appropriate size of the margin, which should carefully be investigated and considered, since sparing of organs at risk or normal tissue is affected. Here we estimate the needed margin for preclinical irradiation by adapting a well-known human margin recipe from van Herck et al. to the dimensions and requirements of the specimen on a small animal radiation research platform (SARRP). We adjusted the factors of the described formula to the specific challenges in an orthotopic pancreatic tumor mouse model to establish an appropriate margin concept. The SARRP was used with its image-guidance irradiation possibility for arc irradiation with a field size of 10 × 10 mm<sup>2</sup> for 5 fractions. Our goal was to irradiate the clinical target volume (CTV) of at least 90% of our mice with at least 95% of the prescribed dose. By carefully analyzing all relevant factors we gain a CTV to planning target volume (PTV) margin of 1.5 mm for our preclinical setup. The stated safety margin is strongly dependent on the exact setting of the experiment and has to be adjusted for other experimental settings. The few stated values in literature correspond well to our result. Even if using margins in the preclinical setting might be an additional challenge, we think it is crucial to use them to produce reliable results and improve the efficacy of radiotherapy.</p>","PeriodicalId":54397,"journal":{"name":"Zeitschrift fur Medizinische Physik","volume":" ","pages":"533-541"},"PeriodicalIF":2.4,"publicationDate":"2024-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"9516757","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}