Dane, I’ve attached a copy of the tiff I’m trying to display. -Amy From: Dane Springmeyer [mailto:[email protected]] Sent: Saturday, April 09, 2011 8:45 PM To: Maples, Amy C. Cc: [email protected] Subject: Re: [Mapnik-users] Displaying raster data with mapnik2
Hi Amy. Would it be possible to post your raster data as well? (or at least a small clip of it). Dane --- \o/ --- Sent from my phone On Apr 7, 2011, at 11:33 AM, "Maples, Amy C." <[email protected]<mailto:[email protected]>> wrote: Mapnik-users, I’m trying to display a raster geotiff image using mapnik2. Based on all the documentation and example code, I’ve written an xml file to create a layer with the appropriate style. In order to test this, I have a simple python script. However, when I call “mapnik2.load_map(m, mapfile)” I get the message “Raster Plugin: exception caught”, which is a little vague. The resulting plot is empty whether I’m using the “us” or “nozoom” levels in the python script (different types of bbox coordinates). I’ve used this python script before for plotting vector data with mapnik, so I believe the problem is with the xml file. Any help on what I might be overlooking would be greatly appreciated. Thanks, Amy Maples Noblis, Inc. Python script: #!/usr/bin/env python # first arguement is xml file name # second argument is zoom level(all lower case)--can be nozoom, us, state, city, locale, street import sys import mapnik2 if not sys.argv[1:]: sys.argv += ['climateTest.xml'] #sys.argv += ["argument1", "argument2", "argument3"] if (len(sys.argv)<=2): zoom='nozoom' else: zoom=sys.argv[2] mapfile = sys.argv[1] map_output = mapfile[:-4]+'_'+zoom+'.png' m = mapnik2.Map(1000,500) mapnik2.load_map(m, mapfile) if (zoom=='nozoom'): bbox = mapnik2.Box2d(mapnik2.Coord(-14911468.431, 2242177.475), mapnik2.Coord(-6158842.809, 8775827.587)) elif (zoom=='us'): bbox = mapnik2.Box2d(mapnik2.Coord(-133.952,219.444), mapnik2.Coord(-55.577,61.647)) elif (zoom=='state'): bbox = mapnik2.Box2d(mapnik2.Coord(-84,37), mapnik2.Coord(-76,40)) elif (zoom=='city'): bbox = mapnik2.Box2d(mapnik2.Coord(-77,39), mapnik2.Coord(-77,39)) elif (zoom=='locale'): bbox = mapnik2.Box2d(mapnik2.Coord(-77.242,38.8399), mapnik2.Coord(-77.195,38.873)) elif (zoom=='street'): bbox = mapnik2.Box2d(mapnik2.Coord(-77.219253,38.858257), mapnik2.Coord(-77.211013,38.866544)) else: bbox = mapnik2.Box2d(mapnik2.Coord(-180.0, -90.0), mapnik2.Coord(180.0, 90.0)) m.zoom_to_box(bbox) mapnik2.render_to_file(m, map_output) climateTest.xml file: <?xml version="1.0" encoding="utf-8"?> <!DOCTYPE Map> <Map srs="+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs"> <Style name="tk"> <Rule> <RasterSymbolizer> <RasterColorizer> <ColorBand value="250" color="rgb(238,130,238)" label="<0"></ColorBand> <ColorBand value="260" color="rgb(75,0,130)" midpoints="199"></ColorBand> <ColorBand value="270" color="rgb(0,0,255)" midpoints="199"></ColorBand> <ColorBand value="280" color="rgb(0,128,0)" midpoints="199"></ColorBand> <ColorBand value="290" color="rgb(255,205,0)" midpoints="199"></ColorBand> <ColorBand value="300" color="rgb(255,127,0)" midpoints="199"></ColorBand> <ColorBand value="310" color="rgb(196,2,51)"></ColorBand> </RasterColorizer> </RasterSymbolizer> </Rule> </Style> <Layer name="RegCM_tk_NCAR_historical_1" srs="+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs"> <StyleName>tk</StyleName> <Datasource> <Parameter name="type">raster</Parameter> <Parameter name="file">/home/amaples/new.tif</Parameter> <Parameter name="lox">-14911468.431</Parameter> <Parameter name="loy">2242177.475</Parameter> <Parameter name="hix">-6158842.809</Parameter> <Parameter name="hiy">8775827.587</Parameter> </Datasource> </Layer> </Map> GDALINFO new.tif: Driver: GTiff/GeoTIFF Files: new.tif Size is 71, 53 Coordinate System is: PROJCS["unnamed", GEOGCS["unnamed ellipse", DATUM["unknown", SPHEROID["unnamed",6378137,0]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433]], PROJECTION["Mercator_1SP"], PARAMETER["central_meridian",0], PARAMETER["scale_factor",1], PARAMETER["false_easting",0], PARAMETER["false_northing",0], UNIT["metre",1, AUTHORITY["EPSG","9001"]]] Origin = (-14911468.430740585550666,8775827.586873944848776) Pixel Size = (123276.417201214062516,-123276.417201214062516) Metadata: AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left (-14911468.431, 8775827.587) (133d57'7.20"W, 61d38'47.76"N) Lower Left (-14911468.431, 2242177.475) (133d57'7.20"W, 19d44'21.62"N) Upper Right (-6158842.809, 8775827.587) ( 55d19'32.97"W, 61d38'47.76"N) Lower Right (-6158842.809, 2242177.475) ( 55d19'32.97"W, 19d44'21.62"N) Center (-10535155.620, 5509002.531) ( 94d38'20.09"W, 44d16'50.94"N) Band 1 Block=71x14 Type=Float64, ColorInterp=Gray _______________________________________________ Mapnik-users mailing list [email protected]<mailto:[email protected]> https://lists.berlios.de/mailman/listinfo/mapnik-users
_______________________________________________ Mapnik-users mailing list [email protected] https://lists.berlios.de/mailman/listinfo/mapnik-users
