<div dir="ltr">No worries.<div><br></div><div>I compiled a recent version of raster2pgsql just to check, but got the same issue:</div><div><br></div><div>$ raster2pgsql<br>RELEASE: 3.5.0 GDAL_VERSION=36 (d2c3ca4)</div><div><br></div><div>Just to reiterate:</div><div>When using raster2pgsql I get 'WARNING'..:</div><div><br></div><div>WARNING: Raster with different alignment found in the set of rasters being converted to PostGIS raster<br></div><div><br></div><div>and 'NOTICE' after the ANALYZE is run:</div><div><br></div><div>NOTICE:  Adding SRID constraint</div>NOTICE:  Adding scale-X constraint<br>NOTICE:  Adding scale-Y constraint<br>NOTICE:  Adding blocksize-X constraint<br>NOTICE:  Adding blocksize-Y constraint<br>NOTICE:  Adding alignment constraint<br>NOTICE:  The rasters (pixel corner coordinates) are not aligned<br>NOTICE:  Unable to add constraint: enforce_same_alignment_rast<br>NOTICE:  SQL used for failed constraint: ALTER TABLE public.merge_raster2pgsql ADD CONSTRAINT enforce_same_alignment_rast CHECK (st_samealignment(rast, '0100000000000000000000494000000000000049C00000000050B3FCC0000000C02246584100000000000000000000000000000000E964000001000100'::raster))<br>NOTICE:  Returned error message: check constraint "enforce_same_alignment_rast" of relation "merge_raster2pgsql" is violated by some row (23514)<br>WARNING:  Unable to add constraint: 'same_alignment'.  Skipping<br>NOTICE:  Adding number of bands constraint<br>NOTICE:  Adding pixel type constraint<br>NOTICE:  Adding nodata value constraint<br>NOTICE:  Adding out-of-database constraint<br>NOTICE:  Adding maximum extent constraint<div><br></div><div>This does not occur when I make a virtual raster (vrt) with gdalbuildvrt, and use that as input to raster2pgsql; no warning or notices, and the data looks good in QGIS.</div><div><br></div><div>Best,</div><div>Andreas<br><div><br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Dec 15, 2024 at 12:35 AM Regina Obe <<a href="mailto:lr@pcorp.us" target="_blank">lr@pcorp.us</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div lang="EN-US"><div><p class="MsoNormal"><span style="font-size:11pt">Sorry haven’t had a chance to look thru this yet.  Going to be tied up for the next two weeks so I might not get a chance until much later.<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt">If any others are interested in investigating, please don’t wait for me.<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt">Thanks,<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt">Regina<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt"><u></u> <u></u></span></p><div style="border-top:none;border-right:none;border-bottom:none;border-left:1.5pt solid blue;padding:0in 0in 0in 4pt"><div><div style="border-right:none;border-bottom:none;border-left:none;border-top:1pt solid rgb(225,225,225);padding:3pt 0in 0in"><p class="MsoNormal"><b><span style="font-size:11pt;font-family:Calibri,sans-serif">From:</span></b><span style="font-size:11pt;font-family:Calibri,sans-serif"> Andreas B <<a href="mailto:panden@gmail.com" target="_blank">panden@gmail.com</a>> <br><b>Sent:</b> Friday, December 13, 2024 6:48 AM<br><b>To:</b> <a href="mailto:postgis-users@lists.osgeo.org" target="_blank">postgis-users@lists.osgeo.org</a><br><b>Subject:</b> Re: Importing tiles with raster2pgsql<u></u><u></u></span></p></div></div><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal">Regina,<u></u><u></u></p><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal" style="margin-bottom:12pt">I'm linking to a folder with a zip containing a subset of tiles that demonstrates my problem, a script to run the commands, and two screenshots.<br>Note that I'm a PostGIS beginner, so it's possible I'm doing something wrong!<u></u><u></u></p></div><div><p class="MsoNormal">The issues encountered differs a bit from what I described yesterday, but the essence is the same.<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">Link to folder: <a href="https://www.dropbox.com/scl/fo/bk8cwktahry0oj0ab2rsz/AHdTD73unb4Efxtgf7exfm0?rlkey=esjtw2iyzysan2hdqklwmj5e3&st=hupv5tsr&dl=0" target="_blank">https://www.dropbox.com/scl/fo/bk8cwktahry0oj0ab2rsz/AHdTD73unb4Efxtgf7exfm0?rlkey=esjtw2iyzysan2hdqklwmj5e3&st=hupv5tsr&dl=0</a><u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">Best,<u></u><u></u></p></div><div><p class="MsoNormal">Andreas<u></u><u></u></p></div></div><p class="MsoNormal"><u></u> <u></u></p><div><div><p class="MsoNormal">On Thu, Dec 12, 2024 at 3:26<span style="font-family:Arial,sans-serif"> </span>PM Regina Obe <<a href="mailto:lr@pcorp.us" target="_blank">lr@pcorp.us</a>> wrote:<u></u><u></u></p></div><blockquote style="border-top:none;border-right:none;border-bottom:none;border-left:1pt solid rgb(204,204,204);padding:0in 0in 0in 6pt;margin-left:4.8pt;margin-right:0in"><div><div><div><p class="MsoNormal"><span style="font-size:11pt">Andreas,</span><u></u><u></u></p><p class="MsoNormal"><span style="font-size:11pt"> </span><u></u><u></u></p><p class="MsoNormal"><span style="font-size:11pt">Yes I would expect them to give the same result.</span><u></u><u></u></p><p class="MsoNormal"><span style="font-size:11pt"> </span><u></u><u></u></p><p class="MsoNormal"><span style="font-size:11pt">What version of raster2pgsql are you running?  It should tell you if you run raster2pgsql without any args.</span><u></u><u></u></p><p class="MsoNormal"><span style="font-size:11pt"> </span><u></u><u></u></p><p class="MsoNormal"><span style="font-size:11pt">Also what platform are you on?  Any chance you have some of those tiles available so we can check it out?</span><u></u><u></u></p><p class="MsoNormal"><span style="font-size:11pt"> </span><u></u><u></u></p><div style="border-top:none;border-right:none;border-bottom:none;border-left:1.5pt solid blue;padding:0in 0in 0in 4pt"><div><div style="border-right:none;border-bottom:none;border-left:none;border-top:1pt solid rgb(225,225,225);padding:3pt 0in 0in"><p class="MsoNormal"><b><span style="font-size:11pt;font-family:Calibri,sans-serif">From:</span></b><span style="font-size:11pt;font-family:Calibri,sans-serif"> Andreas B <<a href="mailto:panden@gmail.com" target="_blank">panden@gmail.com</a>> <br><b>Sent:</b> Thursday, December 12, 2024 6:17 AM<br><b>To:</b> <a href="mailto:postgis-users@lists.osgeo.org" target="_blank">postgis-users@lists.osgeo.org</a><br><b>Subject:</b> Importing tiles with raster2pgsql</span><u></u><u></u></p></div></div><p class="MsoNormal"> <u></u><u></u></p><div><p class="MsoNormal">Hi all,<br><br>I have a dataset with several hundred tiles (GeoTIFFs) that make up a digital terrain model.<br>Each tile is 0.5d x 0.5d.<br><br>I tried to load these GeoTIFFs to PostGIS with:<br><br>$ raster2pgsql -s 32631 -I -M -F -C -t auto -d -l 2,4,8,16,32 tif/*.tif (1)<br><br>The data was loaded, but when viewing in QGIS, it looked like jagged lines with areas of white overlain by pieces of neighboring tiles, etc.<br>This didn't look right.<br>I also noticed that for each tif, the constraints were printed.<br>I expected this to be done once, after the last tile was loaded.<br><br>I then created a virtual raster with gdalbuildvrt:<br><br>$ gdalbuildvrt merge.vrt tif/*.tif (2)<br><br>And then used raster2pgsql to load in the vrt:<br><br>$ raster2pgsql -s 32631 -I -M -F -C -t auto -d -l 2,4,8,16,32 merge.vrt | psql (3)<br><br>The data was loaded, and looked good in QGIS.<br><br>Shouldn't commands (1) and (3) give the same results, or am I misunderstanding?<br><br>Best, Andreas<u></u><u></u></p></div></div></div></div></div></blockquote></div></div></div></div></div></blockquote></div>