210918-放疗gamma通过率matlab代码

最新推荐文章于 2021-12-23 13:59:23 发布

原创

最新推荐文章于 2021-12-23 13:59:23 发布 · 2k 阅读

11 ·

CC 4.0 BY-SA版权

文章标签：

#matlab

本文介绍了一种用于放射治疗的剂量分布分析方法——绝对剂量伽马分析，并提供了详细的MATLAB实现代码。该方法用于比较参考剂量分布图与计算剂量分布之间的差异，以评估放疗计划的准确性。

引文：放疗中常用通过率，一种是使用matrixx软件(比较平面)，ArcCheck软件(比较探头所在平面)，另一种是使用代码，这里仍然来自Niubi高，做此记录，不然下次就找不到了。

gamma2.m

function [Ir,r] = gamma2(img,I,ra,dta,ps)
%绝对剂量伽马分析；img为参考剂量分布图，I为计算剂量分布，Ir为gamma因子分布，r为通过率,ps为pixelspacing
maximg=double(max(img(:))); %小数除以整数会得到错误结果
[sx,sy]=size(img);
img=double(img);I=double(I);ra=double(ra);dta=double(dta);ps=double(ps);
%Ir=zeros(512,512);
Ir=zeros(sx,sy);
detaD=ra*0.01;   %剂量差异标准3%
detad=dta;      %DTA标准3mm 
area=img>maximg*0.05;  %统计有剂量的区域，排除体表外区域  5%
%img=1/max(img(:))*img;  %相对分布
%I=1/max(I(:))*I;
%for i=1:512
    %for j=1:512
for i=1:sx
   for j=1:sy
        
        if area(i,j)==1
                      %对10个像素内的计算,ms,md行的起止处，ns,nd列的起止处
            if i-10<1
                ms=1;  
            else
                ms=i-10;
            end
            if j-10<1
                ns=1;
            else
                ns=j-10;
            end
           % if i+10>512
               % md=