Misschien is dit wel de oorzaak: $ gdalsrsinfo epsg:28992 ERROR 1: gdalsrsinfo was compiled against GDAL 1.9 but current library version is 1.7
Ik ben laatst van GDAL 1.9 terug gegaan naar GDAL 1.7 omdat Mapnik niet meer functioneerde. Tijd om weer eens van A tot Z op te schonen en naar de laatste versies over te stappen. Groeten, Peter Op 2-3-2013 22:24, Gertjan Idema schreef: > Ik heb de indruk dat gdal zelf ook coordinatensystemen bij houdt. > Op mijn systeem vind ik ze op 3 plaatsen: > /usr/share/gdal/1.8 > /usr/share/gdal/1.9 > /usr/share/gdal16 > Geen idee of daar wat mee gedaan wordt. > > Wel vond ik een handig commando: > gdalsrsinfo epsg:28992 > > Bij mij geeft dit: > > PROJ.4 : '+proj=sterea +lat_0=52.15616055555555 > +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 > +ellps=bessel > +towgs84=565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725 > +units=m +no_defs ' > > OGC WKT : > PROJCS["Amersfoort / RD New", > GEOGCS["Amersfoort", > DATUM["Amersfoort", > SPHEROID["Bessel 1841",6377397.155,299.1528128, > AUTHORITY["EPSG","7004"]], > > TOWGS84[565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725], > AUTHORITY["EPSG","6289"]], > PRIMEM["Greenwich",0, > AUTHORITY["EPSG","8901"]], > UNIT["degree",0.0174532925199433, > AUTHORITY["EPSG","9122"]], > AUTHORITY["EPSG","4289"]], > PROJECTION["Oblique_Stereographic"], > PARAMETER["latitude_of_origin",52.15616055555555], > PARAMETER["central_meridian",5.38763888888889], > PARAMETER["scale_factor",0.9999079], > PARAMETER["false_easting",155000], > PARAMETER["false_northing",463000], > UNIT["metre",1, > AUTHORITY["EPSG","9001"]], > AXIS["X",EAST], > AXIS["Y",NORTH], > AUTHORITY["EPSG","28992"]] > > Gertjan > > > On Sat, 2013-03-02 at 22:02 +0100, Peter Peterse wrote: >> Hallo, >> >> het is een proj 4.7.1 zelf gebouwd. Een andere versie staat er niet op. >> Zoals te zien in mijn eerdere mail bevat deze wel de optie +towgs84= >> >> locate epsg vond enkel dit bestand dat relevant was voor dit issue. >> >> Groeten, >> Peter >> >> >> Op 2-3-2013 21:47, Cartinus schreef: >> > Hallo, >> > >> > Staat er niet ook nog een bestand /usr/share/proj/epsg op je computer? >> > Dat is waar dacht ik de meeste linux distributies het laten. De >> > standaard packages bevatten meestal versie 4.7.0 (of ouder) van proj en >> > daarin miste de +towgs84 parameter. >> > >> > m.v.g., >> > Martijn >> > >> > On 03/02/2013 09:31 PM, Peter Peterse wrote: >> >> Hallo Gertjan, >> >> >> >> bedankt voor het antwoord. Het werkt. Zie resultaat >> >> <http://osm.peterse-uithuizen.com/~peter/2deMaasvlakte_in_RD_2.png >> >> <http://osm.peterse-uithuizen.com/%7Epeter/2deMaasvlakte_in_RD_2.png>> >> >> Echter begrijp ik het niet. In het bestand /usr/local/share/proj/epsg >> >> van proj zie ik: >> >> <28992> +proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 >> >> +k=0.999908 +x_0=155000 +y_0=463000 +ellps=bessel +units=m >> >> +towgs84=565.2369,50.0087,465.658,-0.406857330322398,0.350732676542563,-1.8703473836068,4.0812 >> >> +no_defs no_defs >> >> Deze bevat gelijke waarden als jij in het forum noemt. >> >> Een andere schuldige kan ik zo niet vinden. >> >> >> >> Nogmaals bedankt, >> >> Peter >> >> >> >> >> >> Op 2-3-2013 12:52, Gertjan Idema schreef: >> >>> Hallo Peter, >> >>> >> >>> Het lijkt er op de projectie parameters voor RD niet compleet zijn. >> >>> Dat geeft een afwijking van zo'n 100m naar het zuid-zuidwesten. In het >> >>> volgende forum artikel onder #7 vind je een mogelijke work-around: >> >>> forum.openstreetmap.org/viewtopic.php?id=12204 >> >>> >> >>> Laat even weten of het werkt. >> >>> >> >>> Groeten, Gertjan >> >>> >> >>> On Fri, 2013-03-01 at 16:45 +0100, Peter Peterse wrote: >> >>>> Hallo allemaal, >> >>>> >> >>>> ik probeer de 2de maasvlakte in RD tiles weer te geven. Echter zie ik >> >>>> een >> >>>> verschil t.o.v. de officiele site: >> >>>> <http://www.openstreetmap.org/?lat=51.9777&lon=4.001&zoom=14&layers=M> >> >>>> >> >>>> Een screenshot van mijn lokale server is te vinden op: >> >>>> <http://osm.peterse-uithuizen.com/~peter/2deMaasvlakte_in_RD.png >> >>>> <http://osm.peterse-uithuizen.com/%7Epeter/2deMaasvlakte_in_RD.png> >> >>>> <http://osm.peterse-uithuizen.com/%7Epeter/2deMaasvlakte_in_RD.png>> >> >>>> >> >>>> Tussen de weg en het strand bevind zich bij mij water. Ik heb de >> >>>> world_boundaries bijgewerkt en geherprojecteerd met de commando's: >> >>>> >> >>>> extent="225743.58199723 6343505.39917919 816703.80429861 >> >>>> 7125169.31386459" >> >>>> ogr2ogr -f "ESRI Shapefile" -s_srs EPSG:900913 -t_srs EPSG:${srs} \ >> >>>> -spat ${extent} builtup_area_28992.shp builtup_area.shp >> >>>> ogr2ogr -f "ESRI Shapefile" -s_srs EPSG:900913 -t_srs EPSG:${srs} \ >> >>>> -spat ${extent} processed_p_28992.shp processed_p.shp >> >>>> ogr2ogr -f "ESRI Shapefile" -s_srs EPSG:900913 -t_srs EPSG:${srs} \ >> >>>> -spat ${extent} shoreline_300_28992.shp >> >>>> shoreline_300.shp >> >>>> >> >>>> Heeft iemand een idee wat er mis kan zijn? >> >>>> >> >>>> Alvast bedankt, >> >>>> Peter >>
_______________________________________________ Talk-nl mailing list Talk-nl@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-nl
