[postgis-users] estimated extent and EPSG:5514
Sandro Santilli
strk at kbt.io
Tue Nov 22 01:08:01 PST 2016
On Tue, Nov 22, 2016 at 09:29:14AM +0100, Martin Landa wrote:
> Hi,
>
> 2016-11-21 17:15 GMT+01:00 Sandro Santilli <strk at kbt.io>:
>
> hm, I tested simplified workflow on personal computer with GDAL 2.2dev
> and PostGIS 2.4:
>
> $ wget http://vdp.cuzk.cz/vymenny_format/soucasna/20161031_ST_UKSH.xml.gz
> $ createdb test_5514
> $ psql test_5514 -c"create extension postgis"
> $ wget https://epsg.io/5514.sql
> $ psql test_5514 -f 5514.sql
> $ ogr2ogr -f PostgreSQL PG:dbname=test_5514 20161031_ST_UKSH.xml.gz kraje
> $ psql test_5514 -c"select st_extent(originalnihranice) from kraje"
> BOX(-904584.86 -1227295.83,-431724.3 -935236.59)
> psql test_5514 -c"vacuum analyze kraje"
> psql test_5514 -c"select st_estimatedextent('kraje', 'originalnihranice')"
> BOX(-906949.1875 -1228756.125,-429359.96875 -933776.25)
>
> Then on the server (where I found this issue) I am able to reproduce the error:
>
> sudo psql test_5514 -U postgres -c"select st_estimatedextent('kraje',
> 'originalnihranice')"
> st_estimatedextent
> -----------------------------------------------------------------
> BOX(-909107.8125 -1233432.375,4522.92431640625 6136.4794921875)
>
> On the server I have GDAL 2.1.2 and PostgreSQL 9.1, PostGIS 2.1.1.
Is "originalnihranice" a raster column ? If not, forget about GDAL
version. Can you reproduce locally by upgrading to PostGIS 2.4 ?
--strk;
More information about the postgis-users
mailing list