Hi Ghislain, you could try increasing the "SAMPLE_STEPS" warp option, e.g. gdalwarp -wo "SAMPLE_STEPS=210" which should improve the localization of the source bounding window. Choosing SAMPLE_STEPS as a divisor of the source datasets x-resolution should help avoid this discretization artifact altogether. bye, Manuel