[gdal-dev] Misaligned Geotiff/MrSid

[Brandon Boldenow]

This is a custom software build, so it's my own making. But you're right, I
think I have my projections messed up. Thanks for pointing me in the right
direction. Hopefully it'll help me sort it all out.

On Tue, May 11, 2021 at 1:50 PM Even Rouault <even.rouault at spatialys.com>
wrote:

> 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>
> 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) 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/K1qUkot.png
>> https://pasteboard.co/K1qUsAW.png
>>
>> --
>>
>>
>> _______________________________________________
>> gdal-dev mailing listgdal-dev at lists.osgeo.orghttps://lists.osgeo.org/mailman/listinfo/gdal-dev
>>
>> -- 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/6542d1c7/attachment-0001.html>