Importing tiles with raster2pgsql

[Regina Obe]

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.

 

If any others are interested in investigating, please don’t wait for me.

 

Thanks,

Regina

 

From: Andreas B <panden at gmail.com> 
Sent: Friday, December 13, 2024 6:48 AM
To: postgis-users at lists.osgeo.org
Subject: Re: Importing tiles with raster2pgsql

 

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 <https://www.dropbox.com/scl/fo/bk8cwktahry0oj0ab2rsz/AHdTD73unb4Efxtgf7exfm0?rlkey=esjtw2iyzysan2hdqklwmj5e3&st=hupv5tsr&dl=0> &st=hupv5tsr&dl=0

 

Best,

Andreas

 

On Thu, Dec 12, 2024 at 3:26 PM Regina Obe <lr at pcorp.us <mailto:lr at pcorp.us> > 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 <panden at gmail.com <mailto:panden at gmail.com> > 
Sent: Thursday, December 12, 2024 6:17 AM
To: postgis-users at lists.osgeo.org <mailto:postgis-users at lists.osgeo.org> 
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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/postgis-users/attachments/20241214/f60c01ab/attachment.htm>