A medical system (100, 300) comprising: - a memory storing machine executable instructions; - a computational system, wherein execution of the machine executable instructions causes the computational system to: - receive a time series of a R2-star map for a brain volume of a subject; - receive a stimulus signal descriptive of an occurrence of a sensory stimulus repeatedly provided to the subject, wherein the stimulus signal is synchronized to the time series of the R2-star map; - receive a selection of one or more seed voxels identified in the time series of the R2-star map; - calculate a denoised time series of the R2-star map by applying a temporal filter algorithm to the time series of the R2-star map; - calculate a correlation map for each voxel of the one or more seed voxels by calculating a pixel wise correlation between each voxel of the one or more seed voxels and the denoised time series of the R2-star map; - determine an activated region of the brain volume by combining voxels identified in the correlation map for each voxel of the one or more seed voxels above a predetermined threshold; - provide a hemodynamic response function for each voxel and each occurrence of the sensory stimulus in the activated region of the brain volume by aligning the time series of the R2-star map with the stimulus signal; and - provide a subject specific hemodynamic response function by averaging the hemodynamic response function for each voxel and each occurrence of the sensory stimulus in the activated region of the brain volume. - a memory storing machine executable instructions; - a computational system, wherein execution of the machine executable instructions causes the computational system to: - receive a time series of a R2-star map for a brain volume of a subject; - receive a stimulus signal descriptive of an occurrence of a sensory stimulus repeatedly provided to the subject, wherein the stimulus signal is synchronized to the time series of the R2-star map; - receive a selection of one or more seed voxels identified in the time series of the R2-star map; - calculate a denoised time series of the R2-star map by applying a temporal filter algorithm to the time series of the R2-star map; - calculate a correlation map for each voxel of the one or more seed voxels by calculating a pixel wise correlation between each voxel of the one or more seed voxels and the denoised time series of the R2-star map; - determine an activated region of the brain volume by combining voxels identified in the correlation map for each voxel of the one or more seed voxels above a predetermined threshold; - provide a hemodynamic response function for each voxel and each occurrence of the sensory stimulus in the activated region of the brain volume by aligning the time series of the R2-star map with the stimulus signal; and - provide a subject specific hemodynamic response function by averaging the hemodynamic response function for each voxel and each occurrence of the sensory stimulus in the activated region of the brain volume. - receive a time series of a R2-star map for a brain volume of a subject; - receive a stimulus signal descriptive of an occurrence of a sensory stimulus repeatedly provided to the subject, wherein the stimulus signal is synchronized to the time series of the R2-star map; - receive a selection of one or more seed voxels identified in the time series of the R2-star map; - calculate a denoised time series of the R2-star map by applying a temporal filter algorithm to the time series of the R2-star map; - calculate a correlation map for each voxel of the one or more seed voxels by calculating a pixel wise correlation between each voxel of the one or more seed voxels and the denoised time series of the R2-star map; - determine an activated region of the brain volume by combining voxels identified in the correlation map for each voxel of the one or more seed voxels above a predetermined threshold; - provide a hemodynamic response function for each voxel and each occurrence of the sensory stimulus in the activated region of the brain volume by aligning the time series of the R2-star map with the stimulus signal; and - provide a subject specific hemodynamic response function by averaging the hemodynamic response function for each voxel and each occurrence of the sensory stimulus in the activated region of the brain volume.