[mapserver-users] Mapserver does not use the spatial index with GeoPackage

Wed Dec 14 06:32:13 PST 2016

Right, using simply DATA "table_name" makes it fast for me. I suppose that when data source is GeoPackage then typical WMS queries with DATA "SELECT..." would usually be much faster if  spatial index is utilized and almost never much slower.

-Jukka-

Even Rouault wrote:
Re: [mapserver-users] Mapserver does not use the spatial index with GeoPackage

On mercredi 14 décembre 2016 12:00:01 CET Rahkonen Jukka (MML) wrote:

> Hi,

>

> I am using Mapserver 7.0 with GDAL 2.1.1, released 2016/07/07

> I have Geopackage as source data:

>

> CONNECTIONTYPE OGR

> CONNECTION "/users/wms/data/vektorit/kiinteistot/palsta.gpkg"

> DATA "select * from palsta"

Using DATA "palsta" should likely make the RTree to be used. The "SELECT ..." case could probably be also improved in that case, but that depends on whether the OGR driver supports efficient spatial filtering on a layer resulting from a SQL SELECT (which should be the case in the GPKG case)

The GPKG case could also potentially benefit from similar optimizations as the ones I did recently for the Spatialite case (would require checking that a recent enough Spatialite lib is available so that the ST_ operators can be used).

>

> Geopackage is created with GDAL 2.2 with default options and it does have a

> spatial index. Unfortunately it seems that it is not utilized:

>

>

> [Wed Dec 14 13:51:12 2016].928976

> msOGRFileOpen(/users/wms/data/vektorit/kiinteistot/palsta.gpkg)... [Wed Dec

> 14 13:51:12 2016].929040

> OGROPen(/users/wms/data/vektorit/kiinteistot/palsta.gpkg) [Wed Dec 14

> 13:51:12 2016].940209

> msConnPoolRegister(palstat,/users/wms/data/vektorit/kiinteistot/palsta.gpkg

> ,0x94a050) [Wed Dec 14 13:51:12 2016].942224 msOGRFileWhichShapes: SQL =

> select * from palsta. [Wed Dec 14 13:51:12 2016].947950 msOGRFileNextShape:

> Returning shape=0, tile=0 [Wed Dec 14 13:51:12 2016].948047

> msOGRFileNextShape: Returning shape=1, tile=0 [Wed Dec 14 13:51:12

> 2016].948098 msOGRFileNextShape: Returning shape=2, tile=0 [Wed Dec 14

> 13:51:12 2016].948132 msOGRFileNextShape: Returning shape=3, tile=0 [Wed

> Dec 14 13:51:12 2016].948173 msOGRFileNextShape: Returning shape=4, tile=0

> [Wed Dec 14 13:51:12 2016].948215 msOGRFileNextShape: Returning shape=5,

> tile=0 [Wed Dec 14 13:51:12 2016].948291 msOGRFileNextShape: Returning

> shape=6, tile=0 [Wed Dec 14 13:51:12 2016].948340 msOGRFileNextShape:

> Returning shape=7, tile=0 [Wed Dec 14 13:51:12 2016].948405

> msOGRFileNextShape: Returning shape=8, tile=0 [Wed Dec 14 13:51:12

> 2016].948464 msOGRFileNextShape: Returning shape=9, tile=0

>

> Is there something that I should do for forcing Mapserver/OGR to use the

> R-tree index?

>

> -Jukka Rahkonen-

--

Spatialys - Geospatial professional services

http://www.spatialys.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20161214/643d524c/attachment-0001.html>