Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Clip data at 85S/85N when importing to EPSG:3857 #3

Open
olt opened this issue Sep 17, 2013 · 3 comments
Open

Clip data at 85S/85N when importing to EPSG:3857 #3

olt opened this issue Sep 17, 2013 · 3 comments

Comments

@olt
Copy link
Member

olt commented Sep 17, 2013

OSM data at the poles for example will be inserted with invalid EPSG:3857 values, which will result in NaN values when transformed back to EPSG:4326 (see #2).

@tds4u
Copy link

tds4u commented Jan 23, 2014

Are there any news about this? Planet import not even possible with this and if we run "-optimize" it says e.g.

pq: relation "osm_landusages_geom_geohash" already exists

Maybe an option to force optimization with igonoring errors will solve this?

@olt
Copy link
Member Author

olt commented Jan 23, 2014

On 23.01.2014, at 07:16, tds4u wrote:

Are there any news about this? Planet import not even possible with this and if we run "-optimize" it says e.g.

No, this is not implemented. But a planet import should work without clustering (-optimize).

pq: relation "osm_landusages_geom_geohash" already exists

Maybe an option to force optimization with igonoring errors will solve this?

You can't run -optimize a second time.
Implementing a force option would be more work than implementing the clipping.

@tds4u
Copy link

tds4u commented Jan 23, 2014

We can make indizes and clustering on tables by our own. For deleting mistaken polygons we can use:

DELETE FROM osm_transport_points WHERE NOT ST_Intersects(geometry, ST_GeomFromText('POLYGON((-20037508.34 20037508.34, -20037508.34 20037508.34, 20037508.34 20037508.34, 20037508.34 -20037508.34))', '3857'))

and

DELETE FROM osm_places WHERE NOT ST_Intersects(geometry, ST_GeomFromText('POLYGON((-20037508.34 20037508.34, -20037508.34 20037508.34, 20037508.34 20037508.34, 20037508.34 -20037508.34))', '3857'))

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

2 participants