I did exactly the same command (and many variations of that) and still there is deviation, when I compare e.g. with OSM over WMS
Last idea: perhaps my gdal installation (I have played with proj.4 config files some time ago) was broken, tried on another machine. And this helped :) /Jaak 2010/1/12 Rahkonen Jukka <[email protected]>: > Hi, > > I tried your command with shapefile output > ogr2ogr -f "ESRI Shapefile" > -s_srs "+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 > +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m > +nadgri...@null +no_defs +over" > -t_srs EPSG:4326 coastline_wgs84.shp processed_p.shp > > Result suits fine with other maps I have. I could check the data up till 70° > North. Perhaps there is something wrong in ogr2ogr handling MapInfo > projection system. Have a try by converting data first to shapefile and read > that in with MapInfo. > > -Jukka Rahkonen- > > Jaak Laineste wrote: >> >> Does it really work? I converted processed_p file with ogr2ogr into >> MapInfo format (as I needed to use it with MapInfo Pro) with following >> command: >> >> ogr2ogr -f "MapInfo File" -s_srs "+proj=merc +a=6378137 +b=6378137 >> +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgri...@null >> +no_defs +over" -t_srs EPSG:4326 processed_mi_wgs processed_p.shp >> >> But the result is about 20 km north of the correct position (in >> Europe, it depends on latitude, in equator error is 0). Looks like an >> ellipsoid error for me. >> >> Actually my final target is to create similar polygons for sea (to >> use as sea mask in top of raster data). Similar tiled way like >> processed_p has for land. Like "invert" the data. Maybe someone >> somewhere has done it already and can share the result? >> >> /Jaak >> >> 2009/12/18 Frederik Ramm <[email protected]>: >> > Hi, >> > >> > Alexrk wrote: >> >>> The projection is spherical mercator (900913), you can >> download a .PRJ >> >>> file from the menu on >> >>> >> >>> http://spatialreference.org/ref/sr-org/6627/ >> >> >> >> I've tried it in ArcMap - doesn't fit with a regular >> Mercator .prj file. >> > >> > Yes; "regular Mercator" is something other than "spherical >> Mercator". >> > >> >> I don't know, what projection this file actually is :( >> > >> > processed_p.shp is spherical mercator (properly EPSG:3857, formerly >> > EPSG:900913, and erroneously EPSG:3857). >> > >> > If you don't like it, you can use ogr2ogr to reproject the shape to >> > something else. >> > >> > Bye >> > Frederik >> > >> > -- >> > Frederik Ramm ## eMail [email protected] ## >> N49°00'09" E008°23'33" >> > >> > _______________________________________________ >> > talk mailing list >> > [email protected] >> > http://lists.openstreetmap.org/listinfo/talk >> > >> >> >> >> -- >> Jaak Laineste >> www.nutiteq.com >> Skype: jaakl3000 >> Mobile: +372 509 2586 >> _______________________________________________ >> Mapnik-users mailing list >> [email protected] >> https://lists.berlios.de/mailman/listinfo/mapnik-users >> > -- Jaak Laineste www.nutiteq.com Skype: jaakl3000 Mobile: +372 509 2586 _______________________________________________ Mapnik-users mailing list [email protected] https://lists.berlios.de/mailman/listinfo/mapnik-users
