Large-scale distributed analyses of over 30,000 magnetic resonance imaging scans recently detected common genetic variants associated with the volumes of subcortical brain structures. Scaling up these efforts, still greater computational challenges arise in screening the genome for statistical associations at each voxel in the brain, localizing effects using "image-wide genome-wide" testing (voxelwise genome-wide association studies, vGWASs). Here we benefit from distributed computations at multiple sites to metaanalyze genome-wide image-wide data, allowing private genomic data to stay at the site where it was collected. Site-specific tensor-based morphometry is performed with a custom template for each site, using a multichannel registration. A single vGWAS testing 10⁷ variants against 2million voxels can yield hundreds of terabytes (TB) of summary statistics, which would need to be transferred and pooled for metaanalysis. We propose a two-step method, which reduces data transfer for each site to a subset of single-nucleotide polymorphisms and voxels guaranteed to contain all significant hits.