Hi Flavio,<br><br>even with a depth of between 6 and 10 I get no matches.<br><br>D:\bart\FWTools1.3.6\bin>shptree D:\bart\dtb\dtb_subset\dtb2000_sym.shp 10<br>creating index of new LSB format<br><br>D:\bart\FWTools1.3.6\bin>ogrinfo D:\bart\dtb\dtb_subset\dtb2000_sym.shp dtb2000_sym -spat 60000 400000 90000 410000
<br>INFO: Open of `D:\bart\dtb\dtb_subset\dtb2000_sym.shp'<br> using driver `ESRI Shapefile' successful.<br><br>Layer name: dtb2000_sym<br>Geometry: 3D Point<br>Feature Count: 0<br>Extent: (50000.000000, 390000.000000
) - (99997.706000, 419998.443000)<br>Layer SRS WKT:<br>PROJCS["Rijksdriehoekstelsel_New",<br> GEOGCS["GCS_Amersfoort",<br> DATUM["Amersfoort",<br> SPHEROID["Bessel_1841",
6377397.155,299.1528128]],<br> PRIMEM["Greenwich",0.0],<br> UNIT["Degree",0.017453292519943295]],<br> PROJECTION["Double_Stereographic"],<br> PARAMETER["False_Easting",
155000.0],<br> PARAMETER["False_Northing",463000.0],<br> PARAMETER["Central_Meridian",5.38763888888889],<br> PARAMETER["Scale_Factor",0.9999079],<br> PARAMETER["Latitude_Of_Origin",
52.15616055555555],<br> UNIT["Meter",1.0]]<br>CTE: String (8.0)<br><br>If I delete the qix file all is fine:<br>D:\bart\FWTools1.3.6\bin>del D:\bart\dtb\dtb_subset\dtb2000_sym.qix<br><br>D:\bart\FWTools1.3.6\bin>ogrinfo D:\bart\dtb\dtb_subset\dtb2000_sym.shp dtb2000_
<br>sym -spat 60000 400000 90000 410000<br>INFO: Open of `D:\bart\dtb\dtb_subset\dtb2000_sym.shp'<br> using driver `ESRI Shapefile' successful.<br><br>Layer name: dtb2000_sym<br>Geometry: 3D Point<br>Feature Count: 39100
<br>Extent: (50000.000000, 390000.000000) - (99997.706000, 419998.443000)<br>Layer SRS WKT:<br>PROJCS["Rijksdriehoekstelsel_New",<br> GEOGCS["GCS_Amersfoort",<br> DATUM["Amersfoort",
<br> SPHEROID["Bessel_1841",6377397.155,299.1528128]],<br> PRIMEM["Greenwich",0.0],<br> UNIT["Degree",0.017453292519943295]],<br> PROJECTION["Double_Stereographic"],
<br> PARAMETER["False_Easting",155000.0],<br> PARAMETER["False_Northing",463000.0],<br> PARAMETER["Central_Meridian",5.38763888888889],<br> PARAMETER["Scale_Factor",0.9999079
],<br> PARAMETER["Latitude_Of_Origin",52.15616055555555],<br> UNIT["Meter",1.0]]<br>CTE: String (8.0)<br>OGRFeature(dtb2000_sym):4665<br> CTE (String) = B0307<br> POINT (85063.820000000007 407501.06300000002
8.6)<br><br>OGRFeature(dtb2000_sym):4666<br> CTE (String) = B0307<br> POINT (85056.664000000004 407494.31300000002 8.6)<br><br>OGRFeature(dtb2000_sym):4667<br> CTE (String) = B0307<br> POINT (85060.241999999998 407497.68800000002
8.6)<br><br>I've opened up the following bug:<br><a href="http://trac.osgeo.org/gdal/ticket/1790">http://trac.osgeo.org/gdal/ticket/1790</a><br><br>Best regards,<br>Bart<br><br><div><span class="gmail_quote">On 9/4/07,
<b class="gmail_sendername">Flavio Hendry</b> <<a href="mailto:flavio@tydac.ch">flavio@tydac.ch</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi Bart<br><br>Read this? (I had even MapServer crashing when using default on points):<br><br><a href="http://umn.mapserver.ch/MapServer/en/shapetree.htm">http://umn.mapserver.ch/MapServer/en/shapetree.htm</a>, on the bottom:
<br><br>SYNTAX: shptree <shpfile> [<depth>]<br><br><depth> (optional) is the maximum depth of the index to create, default<br>is 0 meaning that shptree will calculate automatically the default depth<br>(calculated so, that a quadtree cell contains 8 shapes). Do not use the
<br> default depth with point files, here a depth between 6 and 10 works fine.<br><br>Mit freundlichem Gruss / Best Regards<br>Flavio Hendry<br><br>----------------------------------------------------------------<br>TYDAC Web-Site:
<a href="http://www.tydac.ch">http://www.tydac.ch</a><br>TYDAC MapServer: <a href="http://www.mapserver.ch">http://www.mapserver.ch</a><br>----------------------------------------------------------------<br>############ Mit freundlichen Gruessen / Kind Regards
<br>############ mailto:<a href="mailto:flavio@tydac.ch">flavio@tydac.ch</a><br>############ TYDAC AG - <a href="http://www.tydac.ch">http://www.tydac.ch</a><br>#### #### Geographic Information Solutions
<br>#### #### Luternauweg 12 -- CH-3006 Bern<br>############ Tel +41 (0)31 368 0180 - Fax +41 (0)31 368 1860<br>----------------------------------------------------------------<br><br><br>-----Original Message-----
<br>From: Bart van den Eijnden <<a href="mailto:bartvde@GMAIL.COM">bartvde@GMAIL.COM</a>><br>To: <a href="mailto:MAPSERVER-USERS@LISTS.UMN.EDU">MAPSERVER-USERS@LISTS.UMN.EDU</a><br>Date: Mon, 3 Sep 2007 17:21:46 +0200
<br>Subject: [UMN_MAPSERVER-USERS] strange problem with quadtree index<br><br>> Hi list,<br>><br>> I have a 3D POINT shapefile dataset, and if I do a query using<br>> Mapserver (or<br>> ogrinfo for that matter) when the shapefile has a quadtree index file
<br>> associated (level 4), it does not find any features within the box.<br>><br>> If I delete the qix file, it correctly finds all the features in the<br>> box,<br>> but ofcourse performance is poor.<br>>
<br>> Any pointers on this problem? Or is the dataset needed to find out<br>> the cause<br>> of such an issue (it's quite big unfortunately)?<br>><br>> Best regards,<br>> Bart<br><br><br></blockquote>
</div><br>