qgis / QGIS

QGIS is a free, open source, cross platform (lin/win/mac) geographical information system (GIS)
https://qgis.org
GNU General Public License v2.0
10.52k stars 2.99k forks source link

Georeferencer - misalignment #51699

Closed robbitt closed 1 year ago

robbitt commented 1 year ago

What is the bug or the crash?

The Georeferencer output is consistently mis-aligned by about 8,600 meters. The output image appears south of all the control points (see figures below).

This happens every time. I've tried different input images, different transform settings, different control points, different QGIS projects, and different coordinate systems.

My GDAL script from the Georeferencer is:

gdal_translate -of GTiff -gcp 3362.36 -7441.16 -122.402 40.5941 -gcp 6270.55 -996.878 -122.376 40.5425 -gcp 8096.07 -7639.3 -122.354 40.5932 -gcp 1733.51 -5840.23 -122.419 40.5823 -gcp 5345.15 -8702.86 -122.381 40.6028 -gcp 3012.84 -2513.73 -122.408 40.5561 "E:/_DOCS/Projects_Robs_working/Side_Channel_GIS/_Shea_Island/Historic_Imagery/1943-09-15/1AN0000010016.tif" "C:/Users/Owner/AppData/Local/Temp/1AN0000010016.tif"
gdalwarp -r near -order 1 -co COMPRESS=  -t_srs EPSG:4326 "C:/Users/Owner/AppData/Local/Temp/1AN0000010016.tif" "E:/_DOCS/Projects_Robs_working/Side_Channel_GIS/_Shea_Island/Historic_Imagery/1943-09-15/1AN0000010016_modified.tif"

Image_alignment GCP_table

Steps to reproduce the issue

  1. Set both QGIS Project and Georeferencer to EPSG:4326 (WGS 84)

  2. Set Georeferencer transformation Parameters: Transformation type: polynomial 1 (tried many other settings - same result) Resampling method: Nearest Neighbour

  3. Add raster, and set 6 control points relative to a QGIS map

  4. Run

Versions

QGIS version | 3.28.3-Firenze | QGIS code revision | c12bcb2f76c -- | -- | -- | -- Qt version | 5.15.3 Python version | 3.9.5 GDAL/OGR version | 3.6.2 PROJ version | 9.1.1 EPSG Registry database version | v10.076 (2022-08-31) GEOS version | 3.11.1-CAPI-1.17.1 SQLite version | 3.39.4 PDAL version | 2.4.3 PostgreSQL client version | unknown SpatiaLite version | 5.0.1 QWT version | 6.1.6 QScintilla2 version | 2.13.1 OS version | Windows 10 Version 2009   |   |   |   Active Python plugins BulkVectorExport | 1.1 coordinates_converter | 1.0 ee_plugin | 0.0.5 FreehandRasterGeoreferencer | 0.8.3 latlontools | 3.6.6 MagneticDeclination | 2.3.1 pointsamplingtool | 0.5.4 profiletool | 4.2.2 QCopycanvas | 0.5 qfieldsync | v4.2.0 qgis_gee_data_catalog | 0.4.3 qgsAzimuth | 0.9.15 splitmultipart | 1.0.0 db_manager | 0.1.20 grassprovider | 2.12.99 MetaSearch | 0.3.6 processing | 2.12.99 sagaprovider | 2.12.99 QGIS version 3.28.3-Firenze QGIS code revision [c12bcb2f76c](https://github.com/qgis/QGIS/commit/c12bcb2f76c) Qt version 5.15.3 Python version 3.9.5 GDAL/OGR version 3.6.2 PROJ version 9.1.1 EPSG Registry database version v10.076 (2022-08-31) GEOS version 3.11.1-CAPI-1.17.1 SQLite version 3.39.4 PDAL version 2.4.3 PostgreSQL client version unknown SpatiaLite version 5.0.1 QWT version 6.1.6 QScintilla2 version 2.13.1 OS version Windows 10 Version 2009 Active Python plugins BulkVectorExport 1.1 coordinates_converter 1.0 ee_plugin 0.0.5 FreehandRasterGeoreferencer 0.8.3 latlontools 3.6.6 MagneticDeclination 2.3.1 pointsamplingtool 0.5.4 profiletool 4.2.2 QCopycanvas 0.5 qfieldsync v4.2.0 qgis_gee_data_catalog 0.4.3 qgsAzimuth 0.9.15 splitmultipart 1.0.0 db_manager 0.1.20 grassprovider 2.12.99 MetaSearch 0.3.6 processing 2.12.99 sagaprovider 2.12.99 ### Supported QGIS version - [ ] I'm running a supported QGIS version according to the roadmap. ### New profile - [X] I tried with a new QGIS profile ### Additional context _No response_
robbitt commented 1 year ago

I forgot to mention, the input images are un-georeferenced tiffs

robbitt commented 1 year ago

If I run the gdal script copied from the Georeferencer, gdal creates the same mis-aligned output raster.