Hi,
at what resolution are you resampling the images? Have you tried to subset the images before writting them? If you are not going to use some bands, if you remove them using the subset operator, it should be faster (for example, there are two angle bands for every spectral band and perhaps you are not going to use them)