[postgis-users] value of spatial index on 300, 000, 000 point geometries

[Luigi Castro Cardeles]

Hi,

here to (on a windows install):
PostgreSQL 8.3.7
PostGIS 1.3.5

teste=# select st_asewkt(st_transform('SRID=3005;POINT(1240227 495643 5)',
4001));
                       st_asewkt
-------------------------------------------------------
 SRID=4001;POINT(-122.689717977579 49.4265525863074 5)
(1 registro)

teste=# select postgis_full_version();
                                postgis_full_version
-------------------------------------------------------------------------------------
 POSTGIS="1.3.5" GEOS="3.0.3-CAPI-1.4.2" PROJ="Rel. 4.6.1, 21 August 2008"
USE_STATS


Luigi Castro Cardeles


2009/5/6 Paul Ramsey <pramsey at opengeo.org>

> Well, our versions differ, but that means nothing :)
>
> medford=# select postgis_proj_version();
>    postgis_proj_version
> ----------------------------
>  Rel. 4.6.1, 21 August 2008
> (1 row)
>
> I'm also running 1.4SVN. And I can't recall if I installed the grid shift
> files.
>
> P.
>
> On Tue, May 5, 2009 at 10:55 PM, Kevin Neufeld <kneufeld at refractions.net>
> wrote:
> > :) So if proj doesn't have z transformation then ... I have a bogus proj
> > install or corrupt spatial ref definition?
> >
> > postgis=# select postgis_proj_version();
> >  postgis_proj_version  -------------------------
> > Rel. 4.4.9, 29 Oct 2004
> > (1 row)
> >
> > Since my answer is different than yours, something is amiss.
> >
> >
> > Paul Ramsey wrote:
> >>
> >> On Tue, May 5, 2009 at 9:58 PM, Kevin Neufeld <kneufeld at refractions.net
> >
> >> wrote:
> >>
> >>>
> >>>  But (and folks can correct me if I'm wrong since the docs don't
> actually
> >>> specify) transforms operate on *all* dimensions.
> >>> postgis=# select st_asewkt(st_transform('SRID=3005;POINT(1240227 495643
> >>> 5)',
> >>> 4001));
> >>>                            st_asewkt
> >>>  ----------------------------------------------------------------------
> >>> SRID=4001;POINT(-122.689717977579 49.4258574870912 533.550271589169)
> >>> (1 row)
> >>>
> >>
> >> That's news to me and my version of PostGIS:
> >>
> >> medford=# select st_asewkt(st_transform('SRID=3005;POINT(1240227
> >> 495643 5)', 4001));
> >>                       st_asewkt
> >> -------------------------------------------------------
> >>  SRID=4001;POINT(-122.689717977579 49.4265525863074 5)
> >> (1 row)
> >>
> >> My intuition was that reprojection would leave Z alone, since proj4
> >> has no z-transform capabilities.
> >>
> >> P.
> >> _______________________________________________
> >> postgis-users mailing list
> >> postgis-users at postgis.refractions.net
> >> http://postgis.refractions.net/mailman/listinfo/postgis-users
> >>
> >
> > _______________________________________________
> > postgis-users mailing list
> > postgis-users at postgis.refractions.net
> > http://postgis.refractions.net/mailman/listinfo/postgis-users
> >
> _______________________________________________
> postgis-users mailing list
> postgis-users at postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/postgis-users/attachments/20090506/4e50012c/attachment.html>