QgsVectorLayer.select( QgsVectorLayer.extent(), true) does not return all features

[qgib]

Author Name: ersts - (ersts -)
Original Redmine Issue: 2658

Redmine category:vectors
Assignee: nobody -

---

select() does is not returning all features, when using the layer's extent as the input rectangle.

This can be reproduced by grabbing the data provided with ticket #12681. Load the shadefile ( 12 points ) and enter the following in to the python console

lvActiveLayer = qgis.utils.iface.activeLayer()
lvExtent = lvActiveLayer.extent()
lvActiveLayer.select(lvExtent, True)

The western, northern, and southern most points are not selected. I get the same result using the c++ api

However, I have check another point layer ( also a shapefile ) and get similar but inconsistent results. From the python console the southern most point is not selected, but through the c++ api all points are selected.

I am running on stock Ubuntu 10.04, trunk r13316

[qgib]

Author Name: Martin Dobias (@wonder-sk)

---

Peter,

I believe the problem is caused by some subtle differences among the values of floating point numbers of coordinates. For example, extent passed to [[PostGIS]] database is first encoded as a string and then decoded back to double value, so there is obvious space for loosing the precision.

Martin


---

-  resolution  was changed from    to wontfix 
-  status_id  was changed from   Open   to Closed 

[qgib]

Author Name: ersts - (ersts -)

---

Howdy Martin,
I also think it has to do with comparing floats/doubles.

The example I gave happens to be from the same layer. It could just as easily been the extent from layer A being used to select features from layer B. Even if there are other ways to select everything from a layer, it seems that using the layer's extent should always give you all the features.

I am not sure this should just be ignored.

[qgib]

Author Name: Martin Dobias (@wonder-sk)

---

Well, I am not sure what can we do here apart from expanding the extent with some small deltas when doing the select, because the spatial filtering is typically done in third-party libraries.

CGAL library has an interesting article on the topic: http://www.cgal.org/philosophy.html

[qgib]

Author Name: Donkagen2 - (Donkagen2 -)

---

Replying to [comment:3 wonder]:

The original bug #12681 was that eVis Event Browser misses records, sometimes a lot of them.

If a solution isn't found then, for all practical purposes, the Event Browser is terminally broken. Are there not also likely to be many other applications of this function where this bug will cause similar problems because the function doesn't do "what it says on the box"?

[qgib]

Author Name: ersts - (ersts -)

---

no no. There are other ways to select all the features, all is not lost!