Re: [OSM-talk-nl] Probleem met mapnik rendering voor EPSG:28992

2010-03-06 Berichten over hetzelfde onderwerp Martijn van Exel
Handige info dit, hier ga ik in de toekomst zeker nog een keer wat aan hebben 


Fijn weekend,

Martijn van Exel +++ m...@rtijn.org
Laziness – Impatience – Hubris

http://schaaltreinen.nl
twitter: mvexel
skype: mvexel
flickr: rhodes

On Mar 6, 2010, at 1:07 AM, Just van den Broecke wrote:

> Ja dit was het probleem en dus idd op te lossen door de world_boundaries 
> shapes te clippen/herprojecteren. Bijv. met
> 
> ogr2ogr -f "ESRI Shapefile" -s_srs EPSG:3785 -t_srs EPSG:28992  -spat 
> 311523.765594493 6555476.44574815 822461.515529216 7160903.43417988 
> processed_p_nl.shp processed_p.shp
> 
> waarbij extent verkregen kan worden via
> select ST_Extent(ST_Transform(way,3785)) from planet_osm_line;
> st_extent
> --
>  BOX(311523.765594493 6555476.44574815,822461.515529216 7160903.43417988)
> 
> (Ipv EPSG:900913 EPSG:3785 gebruikt.
> )
> De SRS-en dan naar 28992 zetten in inc/layer-shapefiles.xml.inc en voila.
> 
> nogmaals dank en groet,
> 
> Just van den Broecke
> 
> 
> Just van den Broecke wrote:
>> Ja dit is precies waar ik net mee bezig ging: clippen naar RD extent en 
>> dan herprojecteren. Iets als
>> ogr2ogr -f "ESRI Shapefile" -s_srs EPSG:... -t_srs EPSG:28992 -spat 
>> ${extent}  ${dstShp} ${srcShp}
>> 
>> Bedankt voor de hulp alvast! Ik laat weten hoe dit verder verloopt.
>> 
>> groet,
>> 
>> Just van den Broecke
>> 
>> 
>> Lennard wrote:
>>> Just van den Broecke wrote:
>>> 
 Dit en  zijn de enige plekken waar 28992 wordt gebruikt. Tilecache 
 zal iets soortgelijks doen (gaat via WMS 28992 request).
>>> Na een kort overleg op irc, komen we tot de conclusie dat dit 
>>> waarschijnlijk gerelateerd is aan http://trac.mapnik.org/ticket/308.
>>> 
>>> Met andere woorden: proj4 vindt de shapefile niet 100% goed om te 
>>> herprojecteren, en mapnik dropt hem daardoor.
>>> 
>>> Je kunt proberen om de processed_p/shoreline_300 vooraf al te 
>>> herprojecteren naar 28992. Voor de snelheid is dit sowieso aan te raden. 
>>> Eventueel in ogr2ogr clippen naar de extent van RD, want ver buiten onze 
>>> regio zal het er niet zo best meer uitzien.
>>> 
>>> Extent van planet-benelux in RD:
>>> select ST_Extent(ST_Transform(way,28992)) from planet_osm_line;
>>>  st_extent
>>> ---
>>>  BOX(-68410.3434188713 153865.862287649,295804.476664188 636639.887598609)
>>> 
>>> 
>> 
>> 
>> 
>> 
>> ___
>> Talk-nl mailing list
>> Talk-nl@openstreetmap.org
>> http://lists.openstreetmap.org/listinfo/talk-nl
>> 
> 
> 
> -- 
> Just
> 
> Just van den Broecke  j...@justobjects.nl
> Just Objects B.V. tel +31 65 4268627 Skype: justb4
> The Netherlands   http://www.justobjects.nl
> 
> 
> 
> 
> 
> 
> ___
> Talk-nl mailing list
> Talk-nl@openstreetmap.org
> http://lists.openstreetmap.org/listinfo/talk-nl


___
Talk-nl mailing list
Talk-nl@openstreetmap.org
http://lists.openstreetmap.org/listinfo/talk-nl


Re: [OSM-talk-nl] Probleem met mapnik rendering voor EPSG:28992

2010-03-05 Berichten over hetzelfde onderwerp Just van den Broecke
Ja dit was het probleem en dus idd op te lossen door de world_boundaries 
shapes te clippen/herprojecteren. Bijv. met

ogr2ogr -f "ESRI Shapefile" -s_srs EPSG:3785 -t_srs EPSG:28992  -spat 
311523.765594493 6555476.44574815 822461.515529216 7160903.43417988 
processed_p_nl.shp processed_p.shp

waarbij extent verkregen kan worden via
select ST_Extent(ST_Transform(way,3785)) from planet_osm_line;
 st_extent
--
  BOX(311523.765594493 6555476.44574815,822461.515529216 7160903.43417988)

(Ipv EPSG:900913 EPSG:3785 gebruikt.
)
De SRS-en dan naar 28992 zetten in inc/layer-shapefiles.xml.inc en voila.

nogmaals dank en groet,

Just van den Broecke


Just van den Broecke wrote:
> Ja dit is precies waar ik net mee bezig ging: clippen naar RD extent en 
> dan herprojecteren. Iets als
> ogr2ogr -f "ESRI Shapefile" -s_srs EPSG:... -t_srs EPSG:28992 -spat 
> ${extent}  ${dstShp} ${srcShp}
> 
> Bedankt voor de hulp alvast! Ik laat weten hoe dit verder verloopt.
> 
> groet,
> 
> Just van den Broecke
> 
> 
> Lennard wrote:
>> Just van den Broecke wrote:
>>
>>> Dit en  zijn de enige plekken waar 28992 wordt gebruikt. Tilecache 
>>> zal iets soortgelijks doen (gaat via WMS 28992 request).
>> Na een kort overleg op irc, komen we tot de conclusie dat dit 
>> waarschijnlijk gerelateerd is aan http://trac.mapnik.org/ticket/308.
>>
>> Met andere woorden: proj4 vindt de shapefile niet 100% goed om te 
>> herprojecteren, en mapnik dropt hem daardoor.
>>
>> Je kunt proberen om de processed_p/shoreline_300 vooraf al te 
>> herprojecteren naar 28992. Voor de snelheid is dit sowieso aan te raden. 
>> Eventueel in ogr2ogr clippen naar de extent van RD, want ver buiten onze 
>> regio zal het er niet zo best meer uitzien.
>>
>> Extent van planet-benelux in RD:
>> select ST_Extent(ST_Transform(way,28992)) from planet_osm_line;
>>   st_extent
>> ---
>>   BOX(-68410.3434188713 153865.862287649,295804.476664188 636639.887598609)
>>
>>
> 
> 
> 
> 
> ___
> Talk-nl mailing list
> Talk-nl@openstreetmap.org
> http://lists.openstreetmap.org/listinfo/talk-nl
> 


-- 
Just

Just van den Broecke  j...@justobjects.nl
Just Objects B.V. tel +31 65 4268627 Skype: justb4
The Netherlands   http://www.justobjects.nl






___
Talk-nl mailing list
Talk-nl@openstreetmap.org
http://lists.openstreetmap.org/listinfo/talk-nl


Re: [OSM-talk-nl] Probleem met mapnik rendering voor EPSG:28992

2010-03-05 Berichten over hetzelfde onderwerp Just van den Broecke
Ja dit is precies waar ik net mee bezig ging: clippen naar RD extent en 
dan herprojecteren. Iets als
ogr2ogr -f "ESRI Shapefile" -s_srs EPSG:... -t_srs EPSG:28992 -spat 
${extent}  ${dstShp} ${srcShp}

Bedankt voor de hulp alvast! Ik laat weten hoe dit verder verloopt.

groet,

Just van den Broecke


Lennard wrote:
> Just van den Broecke wrote:
> 
>> Dit en  zijn de enige plekken waar 28992 wordt gebruikt. Tilecache 
>> zal iets soortgelijks doen (gaat via WMS 28992 request).
> 
> Na een kort overleg op irc, komen we tot de conclusie dat dit 
> waarschijnlijk gerelateerd is aan http://trac.mapnik.org/ticket/308.
> 
> Met andere woorden: proj4 vindt de shapefile niet 100% goed om te 
> herprojecteren, en mapnik dropt hem daardoor.
> 
> Je kunt proberen om de processed_p/shoreline_300 vooraf al te 
> herprojecteren naar 28992. Voor de snelheid is dit sowieso aan te raden. 
> Eventueel in ogr2ogr clippen naar de extent van RD, want ver buiten onze 
> regio zal het er niet zo best meer uitzien.
> 
> Extent van planet-benelux in RD:
> select ST_Extent(ST_Transform(way,28992)) from planet_osm_line;
>   st_extent
> ---
>   BOX(-68410.3434188713 153865.862287649,295804.476664188 636639.887598609)
> 
> 




___
Talk-nl mailing list
Talk-nl@openstreetmap.org
http://lists.openstreetmap.org/listinfo/talk-nl


Re: [OSM-talk-nl] Probleem met mapnik rendering voor EPSG:28992

2010-03-05 Berichten over hetzelfde onderwerp Lennard
Just van den Broecke wrote:

> Dit en  zijn de enige plekken waar 28992 wordt gebruikt. Tilecache 
> zal iets soortgelijks doen (gaat via WMS 28992 request).

Na een kort overleg op irc, komen we tot de conclusie dat dit 
waarschijnlijk gerelateerd is aan http://trac.mapnik.org/ticket/308.

Met andere woorden: proj4 vindt de shapefile niet 100% goed om te 
herprojecteren, en mapnik dropt hem daardoor.

Je kunt proberen om de processed_p/shoreline_300 vooraf al te 
herprojecteren naar 28992. Voor de snelheid is dit sowieso aan te raden. 
Eventueel in ogr2ogr clippen naar de extent van RD, want ver buiten onze 
regio zal het er niet zo best meer uitzien.

Extent van planet-benelux in RD:
select ST_Extent(ST_Transform(way,28992)) from planet_osm_line;
  st_extent
---
  BOX(-68410.3434188713 153865.862287649,295804.476664188 636639.887598609)


-- 
Lennard

___
Talk-nl mailing list
Talk-nl@openstreetmap.org
http://lists.openstreetmap.org/listinfo/talk-nl


Re: [OSM-talk-nl] Probleem met mapnik rendering voor EPSG:28992

2010-03-05 Berichten over hetzelfde onderwerp Just van den Broecke
Lennard wrote:
> Just van den Broecke wrote:
> 
>> Ik vermoed ergens een mismatch tussen projecties/extents in osm.xml 
>> (mijn  is srs28992) en datasource extent is 
>> 2.307,50.134,8.752,54.087 met PG data in 4326) of moeten de 
>> world_boundaries een herprojectie naar 28992 krijgen ? Iemand enig idee 
>> ? Bedankt en groet,
> 
> De srs op de shapefile layers moet op &srs900913; blijven staan. Dan 
> herprojecteert mapnik naar de  srs.
> 
Deze heb ik niet gewijzigd:
grep srs inc/layer-shapefiles.xml.inc geeft:
  
  
  
Ook als ik "builtup" op srs900913 zet zie ik geen verschil.

In generate_image.py gebruik ik:
prj = mapnik.Projection("+proj=sterea +lat_0=52.156160 
+lon_0=5.387639 +k=0.079 +x_0=155000 +y_0=463000 
+ellps=bessel 
+towgs84=565.237,50.0087,465.658,-0.406857,0.350733,-1.87035,4.0812 
+units=m +no_defs")

Dit en  zijn de enige plekken waar 28992 wordt gebruikt. Tilecache 
zal iets soortgelijks doen (gaat via WMS 28992 request).

groet,

Just





___
Talk-nl mailing list
Talk-nl@openstreetmap.org
http://lists.openstreetmap.org/listinfo/talk-nl


Re: [OSM-talk-nl] Probleem met mapnik rendering voor EPSG:28992

2010-03-05 Berichten over hetzelfde onderwerp Lennard
Just van den Broecke wrote:

> Ik vermoed ergens een mismatch tussen projecties/extents in osm.xml 
> (mijn  is srs28992) en datasource extent is 
> 2.307,50.134,8.752,54.087 met PG data in 4326) of moeten de 
> world_boundaries een herprojectie naar 28992 krijgen ? Iemand enig idee 
> ? Bedankt en groet,

De srs op de shapefile layers moet op &srs900913; blijven staan. Dan 
herprojecteert mapnik naar de  srs.

-- 
Lennard

___
Talk-nl mailing list
Talk-nl@openstreetmap.org
http://lists.openstreetmap.org/listinfo/talk-nl


[OSM-talk-nl] Probleem met mapnik rendering voor EPSG:28992

2010-03-05 Berichten over hetzelfde onderwerp Just van den Broecke
Hallo,

Ik probeer tiles te renderen met Mapnik in NL projectie 28992 (ja met de 
"goede" proj setting!). Alles lijkt goed te gaan. Via TileCache met 
Mapnik backend worden tiles gegenereerd. Deze kloppen exact met andere 
28992 bronnen qua extent/resolutie opgegeven in 
tilecache.cfg/OpenLayers. Het enige probleem is dat, denk ik, in deze 
28992-projectie de world_boundaries niet worden gerenderd. Dit zie ik 
ook als ik met "generate_image.py" probeer. Zie bijv. 
http://www.geoskating.com/nl.28992.png
Als ik in een andere projectie render bijv. 900913 krijg ik wel een 
beter kaartbeeld: http://www.geoskating.com/nl.90013.png

Ik vermoed ergens een mismatch tussen projecties/extents in osm.xml 
(mijn  is srs28992) en datasource extent is 
2.307,50.134,8.752,54.087 met PG data in 4326) of moeten de 
world_boundaries een herprojectie naar 28992 krijgen ? Iemand enig idee 
? Bedankt en groet,

--Just

Just van den Broecke  j...@justobjects.nl
The Netherlands   http://www.justobjects.nl






___
Talk-nl mailing list
Talk-nl@openstreetmap.org
http://lists.openstreetmap.org/listinfo/talk-nl