def newGetMap(self, params):
# HACK: check if the image should be strechted
bbox_ratio = float(params['bbox'][2] - params['bbox'][0]) / float(params['bbox'][3] - params['bbox'][1])
image_ratio = float(params['width']) / float(params['height'])
img_height = params['height']
resize = False
if int(bbox_ratio * 100) != int(image_ratio * 100):
params['height'] = int(params['height'] / bbox_ratio)
resize = True
m = self._buildMap(params)
im = Image(params['width'], params['height'])
render(m, im)
format = PIL_TYPE_MAPPING[params['format']]
if resize:
import Image as PILImage
size = params['width'], params['height']
im = PILImage.open(StringIO(im.tostring(format)))
size = params['width'], img_height
im = im.resize(size)
output = StringIO()
im.save(output, format=format)
return Response(params['format'].replace('8',''), output.getvalue())
return Response(params['format'].replace('8',''), im.tostring(format))
def GetMap(self, params):
m = self._buildMap(params)
im = Image(params['width'], params['height'])
render(m, im)
format = PIL_TYPE_MAPPING[params['format']]
if mapnik_version() >= 200300:
# Mapnik 2.3 uses png8 as default, use png32 for backwards compatibility
if format == 'png':
format = 'png32'
return Response(params['format'].replace('8', ''), im.tostring(format))
def GetMap(self, params):
m = self._buildMap(params)
im = Image(params['width'], params['height'])
render(m, im)
im = fromstring('RGBA', (params['width'], params['height']),
rawdata(im))
fh = StringIO()
im.save(fh, PIL_TYPE_MAPPING[params['format']], quality=100)
fh.seek(0)
return Response(params['format'], fh.read())
def cached_image(self):
"""
If caching is enabled, ensure that only one image object is
generated, otherwise generate a fresh one.
"""
if self._image is not None:
return self._image
i = Image(self.image_width, self.image_height)
if self.cache:
self._image = i
return i
def GetMap(self, params):
if params['bbox'][0] >= params['bbox'][2]:
raise OGCException(
"BBOX values don't make sense. minx is greater than maxx.")
if params['bbox'][1] >= params['bbox'][3]:
raise OGCException(
"BBOX values don't make sense. miny is greater than maxy.")
if params.has_key('styles') and len(params['styles']) != len(
params['layers']):
raise OGCException('STYLES length does not match LAYERS length.')
m = Map(params['width'], params['height'])
if params.has_key('transparent') and params['transparent'] == 'FALSE':
m.background = params['bgcolor']
else:
m.background = Color(0, 0, 0, 0)
maplayers = self.mapfactory.layers
mapstyles = self.mapfactory.styles
for layername in params['layers']:
try:
layer = maplayers[layername]
except KeyError:
raise OGCException('Layer "%s" not defined.' % layername,
'LayerNotDefined')
for stylename in layer.styles:
if stylename in mapstyles.keys():
m.append_style(stylename, mapstyles[stylename])
else:
raise ServerConfigurationError(
'Layer "%s" refers to non-existent style "%s".' %
(layername, stylename))
m.layers.append(layer)
m.zoom_to_box(
Envelope(params['bbox'][0], params['bbox'][1], params['bbox'][2],
params['bbox'][3]))
im = Image(params['width'], params['height'])
render(m, im)
im = fromstring('RGBA', (params['width'], params['height']),
rawdata(im))
fh = StringIO()
im.save(fh, PIL_TYPE_MAPPING[params['format']], quality=100)
fh.seek(0)
return Response(params['format'], fh.read())
def __call__(self, layers, system):
request = system['request']
# get image width and height
width = 256
height = 256
# get image bbox
z = int(request.matchdict['z'])
x = int(request.matchdict['x'])
y = int(request.matchdict['y'])
step = max / (2**(int(z) - 1))
xmin = x * step - max
ymin = max - y * step
xmax = (x + 1) * step - max
ymax = max - (y + 1) * step
bbox = Box2d(xmin, ymax, xmax, ymin)
m = Map(width, height)
load_map(m, abspath_from_asset_spec('osmtm:views/map.xml'))
for l in layers:
m.layers.append(l)
m.zoom_to_box(bbox)
format = request.matchdict['format']
if format == 'png':
im = Image(width, height)
render(m, im, 1, 1)
request.response_content_type = 'image/png'
return im.tostring('png')
elif format == 'json':
grid = Grid(width, height)
render_layer(m, grid, layer=0, fields=['id'])
utfgrid = grid.encode('utf', resolution=4)
return json.dumps(utfgrid)
def _render(value, system):
request = system['request']
if not isinstance(value, tuple):
value = (None, value)
layer_name, collection = value
if not hasattr(collection, 'features'):
raise ValueError('renderer is not passed a feature collection')
# get image width and height
try:
img_width = int(request.params.get('img_width', 256))
except:
request.response_status = 400
return 'incorrect img_width'
try:
img_height = int(request.params.get('img_height', 256))
except:
request.response_status = 400
return 'incorrect img_height'
# get image format
try:
img_format = request.params.get(
'img_format', request.matchdict.get('format', 'png'))
img_format = str(img_format)
except:
request.response_status = 400
return 'incorrect img_format'
# get image bbox
img_bbox = request.params.get('img_bbox',
request.params.get('bbox'))
if img_bbox:
try:
img_bbox = map(float, img_bbox.split(','))
except ValueError:
request.response_status = 400
return 'incorrect img_bbox'
img_bbox = Box2d(*img_bbox)
m = Map(img_width, img_height)
load_map(m, mapfile)
if len(m.layers) == 0:
raise ValueError('no layer in the mapnik map')
# if no layer_name is provided then, by convention, use
# the first layer in the mapnik map
if layer_name is None:
layer_name = m.layers[0].name
layer = self._set_layer_in_map(m, layer_name)
layer.datasource = self._create_datasource(collection)
m.zoom_to_box(img_bbox or layer.envelope())
im = Image(img_width, img_height)
render(m, im, 1, 1)
# get the image format from the request
request.response.content_type = 'image/%s' % img_format
return im.tostring(img_format)
if __name__ == "__main__":
mapfile = "my_styles/mapnik2normal.xml"
map_uri = "map.png"
lat = 49.25
lon = 7.0
if len(sys.argv) == 2:
zoom = int(sys.argv[1])
else:
zoom = 13
imgx = 500
imgy = 400
ll = center_to_bbox(lat, lon, zoom, imgx, imgy)
m = Map(imgx, imgy)
load_map(m, mapfile)
prj = Projection(
"+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 +no_defs +over")
c0 = prj.forward(Coord(ll[0], ll[1]))
c1 = prj.forward(Coord(ll[2], ll[3]))
bbox = Box2d(c0.x, c0.y, c1.x, c1.y)
m.zoom_to_box(bbox)
im = Image(imgx, imgy)
render(m, im)
view = im.view(0, 0, imgx, imgy) # x,y,width,height
view.save(map_uri, 'png')
