Regina, 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. Note that I'm a PostGIS beginner, so it's possible I'm doing something wrong!
The issues encountered differs a bit from what I described yesterday, but the essence is the same. Link to folder: https://www.dropbox.com/scl/fo/bk8cwktahry0oj0ab2rsz/AHdTD73unb4Efxtgf7exfm0?rlkey=esjtw2iyzysan2hdqklwmj5e3&st=hupv5tsr&dl=0 Best, Andreas On Thu, Dec 12, 2024 at 3:26 PM Regina Obe <[email protected]> wrote: > Andreas, > > > > Yes I would expect them to give the same result. > > > > What version of raster2pgsql are you running? It should tell you if you > run raster2pgsql without any args. > > > > Also what platform are you on? Any chance you have some of those tiles > available so we can check it out? > > > > *From:* Andreas B <[email protected]> > *Sent:* Thursday, December 12, 2024 6:17 AM > *To:* [email protected] > *Subject:* Importing tiles with raster2pgsql > > > > Hi all, > > I have a dataset with several hundred tiles (GeoTIFFs) that make up a > digital terrain model. > Each tile is 0.5d x 0.5d. > > I tried to load these GeoTIFFs to PostGIS with: > > $ raster2pgsql -s 32631 -I -M -F -C -t auto -d -l 2,4,8,16,32 tif/*.tif (1) > > 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. > This didn't look right. > I also noticed that for each tif, the constraints were printed. > I expected this to be done once, after the last tile was loaded. > > I then created a virtual raster with gdalbuildvrt: > > $ gdalbuildvrt merge.vrt tif/*.tif (2) > > And then used raster2pgsql to load in the vrt: > > $ raster2pgsql -s 32631 -I -M -F -C -t auto -d -l 2,4,8,16,32 merge.vrt | > psql (3) > > The data was loaded, and looked good in QGIS. > > Shouldn't commands (1) and (3) give the same results, or am I > misunderstanding? > > Best, Andreas >
