% Loads a prescription map into a sparse matrix
function prescTumor = open_presc(batchname)

fid = fopen(batchname,'rb');

prescTumor = [];
prescTumor.x_count = fread(fid,1,'int');
prescTumor.y_count = fread(fid,1,'int');
prescTumor.z_count = fread(fid,1,'int');
prescTumor.Nind = fread(fid,1,'int');  % number of non-zero values
prescTumor.non_zero_indices = fread(fid,prescTumor.Nind,'int=>int');
prescTumor.non_zero_values = fread(fid,prescTumor.Nind,'float=>float');

fclose(fid);