Mosaic prior to any coordinate-system projections. Be sure to use bilinear interpolation when projecting. Warping and subsampling done with Fortran code.