The detection of potential biomarkers of Alzheimers disease (AD) is crucial for its early prediction, diagnosis, and treatment. Voxelwise genome-wide association study (VGWAS) is a commonly used method in imaging genomics and usually applied to detect AD biomarkers in imaging and genetic data. However, existing VGWAS methods entail large computational cost and disregard spatial correlations within imaging data. A novel method is proposed to solve these issues.