id summary reporter owner description type status priority milestone component version severity resolution keywords cc 2615 TPS interpolation doesn't work with negative coordinates Even Rouault Even Rouault "The following patch must be applied in thinplatespline.cpp. The initialization of the xmax and ymax is currently wrong when all GCP have negative x or y. {{{ --- gdal.ori/alg/thinplatespline.cpp 2008-07-21 16:27:40.000000000 +0200 +++ gdal.wrk/alg/thinplatespline.cpp 2008-10-17 08:51:46.000000000 +0200 @@ -201,7 +201,7 @@ // More than 2 points - first we have to check if it is 1D or 2D case - double xmax = FLT_MIN, xmin = FLT_MAX, ymax = FLT_MIN, ymin = FLT_MAX; + double xmax = x[0], xmin = x[0], ymax = y[0], ymin = y[0]; double delx, dely; double xx, yy; double sumx = 0.0f, sumy= 0.0f, sumx2 = 0.0f, sumy2 = 0.0f, sumxy = 0.0f; }}} " defect closed normal 1.5.4 GDAL_Raster 1.4.0 normal fixed tps warmerdam