------ Forwarded Message From: Reinhold Klatt <[email protected]> Date: Thu, 21 Jul 2011 14:41:22 +0200 To: Artem Pavlenko <[email protected]> Conversation: [Mapnik-users] Mapnik and conic projectionsP Subject: Re: [Mapnik-users] Mapnik and conic projectionsP
Hi Artem, I'm using the version 0.7.1 from mapnik.org Reinhold On 21.07.2011 14:23 Uhr, "Artem Pavlenko" <[email protected]> wrote: > Which version of Mapnik are you using? > Artem > > On 21 July 2011 11:42, Reinhold Klatt <[email protected]> wrote: >> Hi all, >> I'm trying to render maps in a conic projection. The problem is, that i get >> a blank image as result. >> If I generate in mercator or latlon there is no problem. I get what I want. >> But I can't generate anything in a conical projection. >> I modified a script I found, to test different projections. Does anybody >> know why it's not working? What's my mistake? Are there parameters that have >> to be set? >> >> Thankyou in advance >> >> >> #!/usr/bin/python >> import mapnik >> >> mapfile = "style9_nodem.xml" >> >> projections = { >> "latlon": "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs", #EPSG:4326 >> "google": "+proj=merc +lon_0=0 +lat_ts=0 +x_0=0 +y_0=0 +a=6378137 >> +b=6378137 +units=m +no_defs", # SR-ORG:95 >> "mercator_world": "+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 >> +datum=WGS84 +units=m +no_defs", # SR-ORG:16 >> "lambert": "+proj=lcc +lat_1=49.83333333333334 +lat_2=51.16666666666666 >> +lat_0=90 +lon_0=0 +x_0=150000 +y_0=5400000 +ellps=intl +pm=brussels >> +units=m +no_defs", >> "albert": "+proj=aea +lat_1=29.83333333333334 +lat_2=45.83333333333334 >> +lat_0=37.5 +lon_0=-96 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m >> +no_defs", >> "pytest": "+proj=lcc +lat_1=45.89893890000052 +lat_2=47.69601440000037 >> +lat_0=46.8 +lon_0=2.33722917 +x_0=600000 +y_0=200000 +a=6378249.145 >> +b=6356514.96582849 +pm=2.337229167 +units=m +no_defs", >> } >> >> ll = (8.487, 49.9857, 8.6237, 50.0622) >> >> imgx = 1000 >> imgy = 1000 >> >> # Render an image for every projection >> for projName, projDef in sorted(projections.items()): >> >> print "Projection: " + projName >> >> # Load the map data >> m = mapnik.Map(imgx, imgy) >> mapnik.load_map(m, mapfile) >> >> # Override projection defined in osm.xml >> m.srs = projDef >> >> # Calculate projected boundaries >> prj = mapnik.Projection(projDef) >> c0 = prj.forward(mapnik.Coord(ll[0], ll[1])) >> c1 = prj.forward(mapnik.Coord(ll[2], ll[3])) >> >> # Apply bounding box >> bbox = mapnik.Envelope(c0.x, c0.y, c1.x, c1.y) >> m.zoom_to_box(bbox) >> >> # Render image >> im = mapnik.Image(imgx, imgy) >> mapnik.render(m, im) >> view = im.view(0, 0, imgx, imgy) >> view.save("img/" + projName + ".png", "png256") >> >> >> _______________________________________________ >> Mapnik-users mailing list >> [email protected] >> https://lists.berlios.de/mailman/listinfo/mapnik-users >> > ------ End of Forwarded Message _______________________________________________ Mapnik-users mailing list [email protected] https://lists.berlios.de/mailman/listinfo/mapnik-users
