[gdal-dev] Misaligned Geotiff/MrSid

[Even Rouault]

ok, so the 2 layers are not in the same projection, nor in the same datum.

You didn't mention which software does the overlay between the 2 layers, 
and what was the amount of shift you observed.

Software based on recent PROJ versions will typically apply HPGN grids 
to do NAD83 to WGS 84 coordinate transformations, when they are 
available, but the amount of correction due to those grids seem to be at 
most of 1 metre.

Even

Le 11/05/2021 à 21:36, Brandon Boldenow a écrit :
> Sorry for not including this earlier, let me know if you need anything 
> else:
>
> Info Geotiff:
>
> Driver: GTiff/GeoTIFF
> Size is 16855, 12232
> PROJCS["Lambert Conformal Conic",
> GEOGCS["NAD83",
> DATUM["North_American_Datum_1983",
> SPHEROID["GRS 1980",6378137,298.257222101004,
> AUTHORITY["EPSG","7019"]],
> AUTHORITY["EPSG","6269"]],
> PRIMEM["Greenwich",0],
> UNIT["degree",0.0174532925199433,
> AUTHORITY["EPSG","9122"]],
> AUTHORITY["EPSG","4269"]],
> PROJECTION["Lambert_Conformal_Conic_2SP"],
> PARAMETER["latitude_of_origin",30.125],
> PARAMETER["central_meridian",-100],
> PARAMETER["standard_parallel_1",30.6666666666667],
> PARAMETER["standard_parallel_2",25.3333333333333],
> PARAMETER["false_easting",0],
> PARAMETER["false_northing",0],
> UNIT["metre",1,
> AUTHORITY["EPSG","9001"]],
> AXIS["Easting",EAST],
> AXIS["Northing",NORTH]]
> Origin = (-368445.13723378151, 238489.02048349666)
> Pixel Size = (42.3354999959360455, -42.336150730462649)
>
> Info MrSid:
>
> Driver: MrSid
> Size is: 28000, 28000
> PROJCS["WGS 84 / UTM zone 13N",
> GEOGCS["WGS 84",
> DATUM["WGS_1984",
> SPHEROID["WGS 84",6378137,298.257223563,
> AUTHORITY["EPSG","7030"]],
> AUTHORITY["EPSG","6326"]],
> PRIMEM["Greenwich",0],
> UNIT["degree",0.0174532925199433],
> AUTHORITY["EPSG","4326"]],
> PROJECTION["Transverse_Mercator"],
> PARAMETER["latitude_of_origin",0],
> PARAMETER["central_meridian",-105],
> PARAMETER["scale_factor",0.9996],
> PARAMETER["false_easting",500000],
> PARAMETER["false_northing",0],
> UNIT["metre",1,
> AUTHORITY["EPSG","9001"]],
> AXIS["Easting",EAST],
> AXIS["Northing",NORTH],
> AUTHORITY["EPSG","32613"]]
> Origin = (767738.00007868814, 3336219.99999975129)
> Pixel Size = (0.5, -0.5)
>
> On Tue, May 11, 2021 at 12:37 PM Even Rouault 
> <even.rouault at spatialys.com <mailto:even.rouault at spatialys.com>> wrote:
>
>     Brandon,
>
>     The output of gdalinfo on 2 misaligned rasters would be more telling.
>
>     Even
>
>     Le 11/05/2021 à 20:22, Brandon Boldenow a écrit :
>>     Hi All,
>>
>>     I've been having an issue where VFR Raster Sectionals (from
>>     faa.gov <http://faa.gov>) and MrSid imagery are misaligned when
>>     compared to my other raster layers.In other software these appear
>>     correctly, so I figure there must be something I'm doing wrong,
>>     but I haven't been able to figure it out.
>>
>>     Do you know what could cause this? I've attached some screenshots
>>     with a reference point to show the misalignment - the black and
>>     white satellite imagery is representative of what I see for most
>>     raster formats.
>>
>>     Thanks in advance!
>>
>>     https://pasteboard.co/K1qU2X2.png
>>     <https://pasteboard.co/K1qU2X2.png>
>>     https://pasteboard.co/K1qUkot.png <https://pasteboard.co/K1qUkot.png>
>>     https://pasteboard.co/K1qUsAW.png <https://pasteboard.co/K1qUsAW.png>
>>
>>     -- 
>>
>>
>>     _______________________________________________
>>     gdal-dev mailing list
>>     gdal-dev at lists.osgeo.org  <mailto:gdal-dev at lists.osgeo.org>
>>     https://lists.osgeo.org/mailman/listinfo/gdal-dev  <https://lists.osgeo.org/mailman/listinfo/gdal-dev>
>
>     -- 
>     http://www.spatialys.com  <http://www.spatialys.com>
>     My software is free, but my time generally not.
>
>
>
> -- 
>
-- 
http://www.spatialys.com
My software is free, but my time generally not.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20210511/e2ae465a/attachment.html>