[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