Re: [OSM-talk-nl] Probleem met mapnik rendering voor EPSG:28992
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
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
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
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
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
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
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