<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<div class="moz-cite-prefix">Misschien is dit wel de oorzaak:<br>
<br>
$ gdalsrsinfo epsg:28992<br>
ERROR 1: gdalsrsinfo was compiled against GDAL 1.9 but current
library version is 1.7<br>
<br>
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.<br>
<br>
Groeten,<br>
Peter<br>
<br>
Op 2-3-2013 22:24, Gertjan Idema schreef:<br>
</div>
<blockquote cite="mid:1362259446.6921.4.camel@laptop-ubuntu"
type="cite">
<meta http-equiv="Content-Type" content="text/html;
charset=ISO-8859-1">
<meta name="GENERATOR" content="GtkHTML/4.1.92">
Ik heb de indruk dat gdal zelf ook coordinatensystemen bij houdt.<br>
Op mijn systeem vind ik ze op 3 plaatsen:<br>
/usr/share/gdal/1.8<br>
/usr/share/gdal/1.9<br>
/usr/share/gdal16<br>
Geen idee of daar wat mee gedaan wordt.<br>
<br>
Wel vond ik een handig commando:<br>
gdalsrsinfo epsg:28992<br>
<br>
Bij mij geeft dit:<br>
<br>
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 '<br>
<br>
OGC WKT :<br>
PROJCS["Amersfoort / RD New",<br>
GEOGCS["Amersfoort",<br>
DATUM["Amersfoort",<br>
SPHEROID["Bessel 1841",6377397.155,299.1528128,<br>
AUTHORITY["EPSG","7004"]],<br>
TOWGS84[565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725],<br>
AUTHORITY["EPSG","6289"]],<br>
PRIMEM["Greenwich",0,<br>
AUTHORITY["EPSG","8901"]],<br>
UNIT["degree",0.0174532925199433,<br>
AUTHORITY["EPSG","9122"]],<br>
AUTHORITY["EPSG","4289"]],<br>
PROJECTION["Oblique_Stereographic"],<br>
PARAMETER["latitude_of_origin",52.15616055555555],<br>
PARAMETER["central_meridian",5.38763888888889],<br>
PARAMETER["scale_factor",0.9999079],<br>
PARAMETER["false_easting",155000],<br>
PARAMETER["false_northing",463000],<br>
UNIT["metre",1,<br>
AUTHORITY["EPSG","9001"]],<br>
AXIS["X",EAST],<br>
AXIS["Y",NORTH],<br>
AUTHORITY["EPSG","28992"]]<br>
<br>
Gertjan<br>
<br>
<br>
On Sat, 2013-03-02 at 22:02 +0100, Peter Peterse wrote:
<blockquote type="CITE">
<pre>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
>> <<a moz-do-not-send="true" href="http://osm.peterse-uithuizen.com/%7Epeter/2deMaasvlakte_in_RD_2.png">http://osm.peterse-uithuizen.com/~peter/2deMaasvlakte_in_RD_2.png</a>>
>> 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:
>>>> <<a moz-do-not-send="true" href="http://www.openstreetmap.org/?lat=51.9777&lon=4.001&zoom=14&layers=M">http://www.openstreetmap.org/?lat=51.9777&lon=4.001&zoom=14&layers=M</a>>
>>>>
>>>> Een screenshot van mijn lokale server is te vinden op:
>>>> <<a moz-do-not-send="true" href="http://osm.peterse-uithuizen.com/%7Epeter/2deMaasvlakte_in_RD.png">http://osm.peterse-uithuizen.com/~peter/2deMaasvlakte_in_RD.png</a> <<a moz-do-not-send="true" href="http://osm.peterse-uithuizen.com/%7Epeter/2deMaasvlakte_in_RD.png">http://osm.peterse-uithuizen.com/%7Epeter/2deMaasvlakte_in_RD.png</a>>>
>>>>
>>>> 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
</pre>
</blockquote>
</blockquote>
<br>
</body>
</html>