WiscPlan_olderMV/opt/linlsqUtil/pet2dose.m

function tumorPresc = pet2dose(petFileName,m)

% if(map=='linear')
%     m=1;
% elseif(map=='root')
%     m=2;
% elseif(map=='quad')
%     m=3;
% elseif(map=='gomp')
%     m=4;
% end
%     
PET = am2geom(petFileName);

TDP = open_presc('C:\Documents and Settings\Steve\My Documents\Research\optNew\linlsqOpt\input\HN003uniPrescVolNew\PTV70dosePlus.bin');

TDPfull = full3D(TDP);

TDPbinary = single(~~TDPfull);

PETcorrected = PET.data.*TDPbinary;

PETsparse = sparse3D(PETcorrected);

for i=1:PETsparse.Nind;
    if(PETsparse.non_zero_values(i)<0)
        PETsparse.non_zero_values(i)=0;
    end
    if(m==1)
        PETsparse.non_zero_values(i)=70.+30./100*70*(PETsparse.non_zero_values(i)./0.00506392);
    elseif(m==2)
        PETsparse.non_zero_values(i)=70.+30./100*70*(PETsparse.non_zero_values(i)./0.00506392)^(.5);
    elseif(m==3)
        PETsparse.non_zero_values(i)=70.+30./100*70*(PETsparse.non_zero_values(i)./0.00506392)^2;
    elseif(m==4)
        PETsparse.non_zero_values(i)=70.+30./100*70*exp(-15.*exp(-15.*(PETsparse.non_zero_values(i)./0.00506392-0.1)));
    end
end
        
tumorPresc = PETsparse;