A method and apparatus for irradiation therapy using voxel based function measurements of organs-at-risk (OAR). The method includes determining size and location of each voxel of a plurality of voxels in a reference frame of a radiation device. The method further includes obtaining measurements that relate to utility of tissue type at each voxel. The method further includes determining a subset of the voxels that enclose an organ-at-risk (OAR) volume. The method further includes determining a value of a utility measure fj at each voxel of the subset based on a corresponding value of the measurements. The method further includes determining a series of beam shapes and intensities which minimize a value of an objective function that is based on a computed dose delivered to an OAR voxel multiplied by the utility measure fj for that voxel summed over all voxels.