Opened 12 years ago
Closed 11 years ago
#2244 closed defect (fixed)
[raster]: ST_SetGeoReference screws up out db
Reported by: | robe | Owned by: | dustymugs |
---|---|---|---|
Priority: | critical | Milestone: | PostGIS 2.1.0 |
Component: | raster | Version: | master |
Keywords: | history | Cc: |
Description
Not sure what causes this issue whether a doc bug or something not right in out db support. I thought we had this probelm before and it was fixed, but I just updated to the latest.
PostgreSQL 9.2.2, compiled by Visual C++ build 1600, 64-bit POSTGIS="2.1.0SVN r11197" GEOS="3.4.0dev-CAPI-1.8.0 r0" PROJ="Rel. 4.8.0, 6 March 2012" GDAL="GDAL 1.9.2, released 2012/10/08" LIBXML="2.7.8" LIBJSON="UNKNOWN" RASTER
I imported this image twice — one as in db and one as outdb http://download.osgeo.org/postgis/logo_suite/adbadge_tall/adbadge_tall.png
SET PGPORT=5442 SET PGHOST=localhost SET PGUSER=postgres SET PGDATABASE=testpostgis21 raster2pgsql C:/pics/adbadge_tall.png test_pele | psql raster2pgsql -e -R -a C:/pics/adbadge_tall.png test_pele | psql
-- After georeferencing SELECT rid, (h).* FROM (SELECT rid, ST_Histogram(ST_SetGeoReference(rast, '1 0 0 -1 445139 5415000'),1,4) As h from test_pele) As a; rid | min | max | count | percent -----+-----+-----+--------+-------------------- 1 | 15 | 75 | 13456 | 0.025542900531511 1 | 75 | 135 | 23572 | 0.0447456340167046 1 | 135 | 195 | 50711 | 0.0962623386484434 1 | 195 | 255 | 439061 | 0.833449126803341 2 | 0 | 0 | 526800 | -1 (5 rows)
But before that they return same answer
--before georeferencing -- SELECT rid, (h).* FROM (SELECT rid, ST_Histogram(rast,1,4) As h from test_pele) As a; rid | min | max | count | percent -----+-----+-----+--------+-------------------- 1 | 15 | 75 | 13456 | 0.025542900531511 1 | 75 | 135 | 23572 | 0.0447456340167046 1 | 135 | 195 | 50711 | 0.0962623386484434 1 | 195 | 255 | 439061 | 0.833449126803341 2 | 15 | 75 | 13456 | 0.025542900531511 2 | 75 | 135 | 23572 | 0.0447456340167046 2 | 135 | 195 | 50711 | 0.0962623386484434 2 | 195 | 255 | 439061 | 0.833449126803341 (8 rows)
haven't tested on 2.0 to see if similar issue.
Change History (6)
comment:1 by , 12 years ago
Component: | postgis → raster |
---|---|
Owner: | changed from | to
comment:2 by , 12 years ago
comment:3 by , 12 years ago
I think we discussed we should just throw an error for outdb then instead of causing mass confusion and having people have to refer to docs when weirdness ensues.
comment:6 by , 11 years ago
Keywords: | history added |
---|---|
Resolution: | → fixed |
Status: | new → closed |
Fixed in r11298. Feel free to change/expand the note I left in the docs.
That is expected behavior. There is an expectation that the metadata stored in database and the actual raster file's metadata match (specifically the geotransform matrix). If the metadata doesn't match, the output from a function (such as ST_Histogram) is unknown. Non-georeferenced rasters are always set the default geotransform (scaleX: 1 scaleY: -1 skewX/Y: 0 upperleftX/Y: 0 srid: 0).