On Sun, 2008-11-02 at 16:15 +0200, Rahkonen Jukka wrote:
> Thaks a lot Jon. I feel stupid, I should have been looking closer to the
> BOX3D. By giving the same BOX3D both queries are returning same count of
> features.
>
> What I think now is that perhaps the real reason for the problem is in the
> data. OpenJUMP is suggesting that the feature set returned by the query that
> leads to missing islands is containing invalid polygons. Does somebody know
> how could I verify that from the database or from the shapefiles?
>
> This query should return wrong geometries:
>
> select asbinary(way) as way FROM shoreline_a where way &&
> setSRID('BOX3D(2796623.487632482 8414856.658170387,3648235.036730606
> 8450106.434262134)'::box3d,900913)
>
> I am sorry about my misleading reports. I hope that at least the data error
> is real.
Postgis will identify invalid polygon geometries for you:
=> select osm_id from planet_osm_polygon where not isvalid(way);
These errors may cause the individual polygons to mis-render but should
not make a difference to other data. The positive ID's listed should
match the ways IDs. Negative values are the IDs of multipolygon
relations.
Jon
_______________________________________________
Mapnik-users mailing list
[email protected]
https://lists.berlios.de/mailman/listinfo/mapnik-users
