RS2: add half-pixel shift to reported GCP line and column numbers (#1668

) (fixes #1666)
OSGeo · Jun 24, 2019 · ddcf131 · ddcf131
1 parent f1f1682
commit ddcf131
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/gdal/frmts/rs2/rs2dataset.cpp b/gdal/frmts/rs2/rs2dataset.cpp
@@ -1313,9 +1313,9 @@ GDALDataset *RS2Dataset::Open( GDALOpenInfo * poOpenInfo )
             psGCP->pszId = CPLStrdup( szID );
             psGCP->pszInfo = CPLStrdup("");
             psGCP->dfGCPPixel =
-                CPLAtof(CPLGetXMLValue(psNode,"imageCoordinate.pixel","0"));
+                CPLAtof(CPLGetXMLValue(psNode,"imageCoordinate.pixel","0")) + 0.5;
             psGCP->dfGCPLine =
-                CPLAtof(CPLGetXMLValue(psNode,"imageCoordinate.line","0"));
+                CPLAtof(CPLGetXMLValue(psNode,"imageCoordinate.line","0")) + 0.5;
             psGCP->dfGCPX =
                 CPLAtof(CPLGetXMLValue(psNode,"geodeticCoordinate.longitude",""));
             psGCP->dfGCPY =