pgGetRast() error

[lucas-johnson]

Context

• I have added several rasters to a pgsql DB using pgWriteRast(..., constraints = FALSE). After all my rasters have been added to the db, I then manually run:

SELECT AddRasterConstraints('my_rasters'::name, 'rast'::name);

• I add constraints after loading all rasters because it takes too much time to add constraints while loading each individual raster.

• Everything is added successfully:

• Retrieving any of the rasters using pgGetRast does not work...

Error

Error in pgSRID(conn, rout@crs) : 
  No SRID matches found. Re-run with 'create.srid = TRUE' to create new SRID entry in spatial_ref_sys.

Question

• Do I really need to re-run with create.srid = TRUE? Or is this a bug?
• Shouldn't running the AddRasterConstraints manually be enough?

[lucas-johnson]

UPDATE:

This error message is printed, but the raster is returned correctly. Should the message be a warning instead?

[dnbucklin]

Hi @lucas-johnson . To answer your questions directly:

• No you shouldn't need to do that. That error message should not be displayed, that is a bug.
• I wouldn't think that your constraints workflow is affecting this. Have you tested pgWriteRast with constraints = TRUE and reloaded without seeing that error?

Since you used pgWriteRast to write to PostGIS, and it looks like SRID was correctly(?) applied as 5070, it's not obvious to me why that error is showing. Could you check the SRID(s) returned from both the original raster, and the pgGetRast-loaded raster?

r1 <- raster("original_raster_file")
pgSRID(conn, r1@crs) # Does this return `5070` ?
rgdal::showEPSG(as.character(r1@crs)) # what does this return?

r2 <- pgGetRast(conn, 'my_rasters')
pgSRID(conn, r2@crs)  # This should produce the error you were seeing...

[dnbucklin]

Error printing fix applied in 74a6ac7 , closing...