I have re-fitted the synthetic images with the FFM method, but this time used weighted least squares to perform the fitting. The weight was 1/sigma in each pixel where sigma is the Poisson SD. The previous bias that was seen for single images appears to be gone now. The scatter for single images is also much smaller.