Multi-modal MRI modalities quantify different, yet complimentary, properties of the brain and its activity. When studied jointly, multi-modal imaging data may improve our understanding of the brain. We aim to study the complex relationships between multiple imaging modalities and map how these relationships vary spatially across different anatomical brain regions. Given a particular location in the brain, we regress an outcome image modality on one or more other modalities using all voxels in a local neighborhood of a target voxel. We apply our method to study how the relationship between local functional connectivity and cerebral blood flow varies spatially.