I had missed something essential and simple ... mode=2 must be added to the end of the PG connection string to tell GDAL to treat a PostGIS raster table as one unit. The default mode=1 treats each tile as a subdataset
Sorry for the noise Best regards Andreas Oxenstierna T-Kartor Geospatial AB Olof Mohlins väg 12 Kristianstad mobile: +46 733 206831 mailto: andreas.oxenstie...@t-kartor.com<mailto:andreas.oxenstie...@t-kartor.com> www.t-kartor.com<http://www.t-kartor.com> On 12 May 2021, 08:12 +0200, Andreas Oxenstierna <andreas.oxenstie...@t-kartor.com>, wrote: Tested with: 1 … table=rastertable WHERE=rid>0" dest.tif Returns: WHERE=rid>0 has no raster bands 2 … table=rastertable WHERE=rid=1" dest.tif Returns: a valid but empty tif file 3 … table=<view based on the rastertable" dest.tif Returns: ... has no raster bands gdalwarp seems to treat PostGIS raster tables as one subdataset for each tile. Basically make them useless in my mind - or have I missed something essential/simple ? Hälsningar Andreas Oxenstierna T-Kartor Geospatial AB Olof Mohlins väg 12 Kristianstad mobile: +46 733 206831 mailto: andreas.oxenstie...@t-kartor.com<mailto:andreas.oxenstie...@t-kartor.com> www.t-kartor.com<http://www.t-kartor.com> On 11 May 2021, 08:48 +0200, Andreas Oxenstierna <andreas.oxenstie...@t-kartor.com>, wrote: Thanks but still stumbling (using 2.3.2) gdalwarp -cutline "PG:dbname=gisdb" -csql 'select * from polytest where id=1' -crop_to_cutline -of GTiff -srcnodata -9999 -dstnodata -9999 "PG:<connectinfo> schema=a table=rastertable" dest.tif returns: ERROR 1: Input file <rastertable> has no raster bands. gdalinfo reports that the table is tiled in 360 pieces which potentially is the source for the issue: Subdatasets: SUBDATASET_1_NAME=PG:host= SUBDATASET_360_DESC=PostGIS Raster at Hälsningar Andreas Oxenstierna T-Kartor Geospatial AB Olof Mohlins väg 12 Kristianstad mobile: +46 733 206831 mailto: andreas.oxenstie...@t-kartor.com<mailto:andreas.oxenstie...@t-kartor.com> www.t-kartor.com<http://www.t-kartor.com> On 10 May 2021, 19:31 +0200, Even Rouault <even.roua...@spatialys.com>, wrote: Try with gdalwarp and the -cutline family of arguments Le 10/05/2021 à 19:23, Andreas Oxenstierna a écrit : Hi list Is it possible to clip a PostGIS raster with a PostGIS vector poly using gdal_translate? I cannot find a working syntax to define the clipping table in the WHERE clause. Eg. this returns "Error browsing database for PostGIS Raster properties : ERROR: syntax error at or near "FROM" gdal_translate -of GTiff "PG:<connectinfo> schema=a table=rastertable WHERE='ST_Intersects(rast,c.geom) AND c.id<http://c.id> = 2 FROM a.vectortable c'" a.tif Best regards Andreas Oxenstierna T-Kartor Geospatial AB Olof Mohlins väg 12 Kristianstad mobile: +46 733 206831 mailto: andreas.oxenstie...@t-kartor.com<mailto:andreas.oxenstie...@t-kartor.com> www.t-kartor.com<http://www.t-kartor.com> _______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org<mailto:gdal-dev@lists.osgeo.org> https://lists.osgeo.org/mailman/listinfo/gdal-dev -- http://www.spatialys.com My software is free, but my time generally not.
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev
