Iterative polynomial fitting along image rows and columns has recently been used to remove curvature bias in multispectral image sets of the human forearm and phantoms. However, this method is only applicable if foreground and background features satisfy strong separation conditions. In this method, we verify that the iterative polynomial approach converges toward bivariate polynomial fitting, and, hence, the resulting fit corresponds to low-pass filtering the image. In contrast to the iterative fitting, the bivariate polynomial fit can be performed on images with missing or excluded parts. Indeed, our observation enables us to modify the scheme and significantly weaken the required assumptions on foreground/background separation allowing a wider range of application.