Paul,
Ah, yes. You're right. I forgot that 3dNetCorr wonderfully gives me correlation matrices based on every region of interest's average time series. Thank you for patiently reawakening those circuits in my brain.
Yes, my main result will be the plot generated by 3dLME from the region-to-region correlation matrices that will be fed into 3dLME from 3dNetCorr via the fetal fat_lme_prep.py program. Hmm, that does make spatial smoothing moot in a way. I will need to ponder that a bit as well.
Sincerely,
Dante
P.S.: Although the region-to-region correlations will be my main result, it does seem to me at present that the region-to-voxelwise correlations that I am getting via the -ts_wb_Z option must also be used to show whole-brain maps in the publication for a region or two, so it's not totally moot.