WiscPlan_olderMV/opt/linlsqOpt/pet2dose.asv

function tumorPresc = pet2dose(petFileName,map)

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\Desktop\optNew\linlsqOpt\input\uniformPresc\TumorDosePlus.bin');

TDPfull = full3D(TDP);

TDPbinary = single(~~TDPfull);

PETcorrected = PET.data.*TDPbinary;

PETsparse = sparse3D(PETcorrected);

for i=1:PETsparse.Nind;
    if(m==1)
        PETsparse.non_zero_values(i)=100.+30.*(PETsparse.non_zero_values(i)./0.0372814);
    elseif(m==2)
        PETsparse.non_zero_values(i)=100.+30.*(PETsparse.non_zero_values(i)./0.0372814)^(.5);
    elseif(m==3)
        PETsparse.non_zero_values(i)=100.+30.*(PETsparse.non_zero_values(i)./0.0372814)^2;
    elseif(m==4)
        PETsparse.non_zero_values(i)=100.+30.*exp(-30.*exp(-(PETsparse.non_zero_values(i)./0.0372814-0.25)));
    end
end
        
tumorPresc = PETsparse;

        end