def image_map_gen(cls, records): """Render the image map""" for record in records: # Récupère l'étendu de la zone de garden parcelle, _envelope, _area = get_as_epsg4326([record.parcelle.geom]) # Léger dézoom pour afficher correctement les points qui touchent la bbox envelope = [_envelope[0] - 0.001, _envelope[1] + 0.001, _envelope[2] - 0.001, _envelope[3] + 0.001] if envelope is None: continue # Map title title = u"Plan local du jardin\n" title += date.today().strftime("%02d/%02m/%Y") m = MapRender(1024, 768, envelope, True) # Ajoute le fond de carte m.add_bg() # Ajoute la pracelle m.plot_geom(parcelle[0], record.name, u"Parcelle", color=cls.COLOR, bgcolor=cls.BGCOLOR) data_nl = m.render() m.plot_legend() m.plot_compass() m.plot_scaling() m.plot_title(title) data = m.render() cls.write([record], {"image_map": buffer(data)})
def generate(cls, records): for record in records: if record.address is None: continue town, envelope, area = get_as_epsg4326([record.address.my_city.contour]) # Calcule de la bbox contenant tout les points _envelope = None for points in cls.search([]): _points, envelope, _area = get_as_epsg4326([points.geom]) if envelope is None: continue _envelope = envelope_union(envelope, _envelope) # Léger dézoom pour afficher correctement les points qui touchent la bbox envelope = [ _envelope[0] - 0.001, _envelope[1] + 0.001, _envelope[2] - 0.001, _envelope[3] + 0.001, ] m = MapRender(640, 480, envelope, True) m.add_bg() for entry in cls.search([]): points, _envelope, _area = get_as_epsg4326([entry.geom]) if len(points) == 0: continue if record == entry: m.plot_geom(points[0], None, None, color=cls.COLOR, bgcolor=cls.BGCOLOR) else: m.plot_geom(points[0], None, None, color=(0, 0, 1, 1), bgcolor=cls.BGCOLOR) data = m.render() cls.write([record], {'image_map': buffer(data)})
def generate(cls, records): for record in records: if record.num is None: continue # Récupère les placettes de mesure du dispositif EmpObj = Pool().get(record.__name__) objs = EmpObj.search([('dispositif', '=', record.dispositif.id)]) pts, envelope, area = get_as_epsg4326([obj.geom for obj in objs]) # Placette en cours points, _envelope, _area = get_as_epsg4326([record.geom]) # Léger dézoom pour afficher correctement les points qui touchent la bbox envelope = [ _envelope[0] - 0.001, _envelope[1] + 0.001, _envelope[2] - 0.001, _envelope[3] + 0.001, ] m = MapRender(640, 480, envelope, True) m.add_bg() # Ajoute les placettes du dispositif for entry in pts: if len(pts) == 0: continue if entry == get_as_epsg4326([record.geom])[0][0]: m.plot_geom(entry, None, None, color=(0, 0, 1, 1), bgcolor=record.BGCOLOR) else: m.plot_geom(entry, None, None, color=(0, 0, 1, 0.5), bgcolor=record.BGCOLOR) m.plot_geom(points[0], str(record.num), None, color=(1, 1, 1, 1), bgcolor=record.BGCOLOR) data = m.render() cls.write([record], {'image_map': buffer(data)})
def generate(cls, records): for record in records: if record.code is None: continue lines, envelope, _line = get_as_epsg4326([record.equipement.geom]) aires, _envelope, _aire = get_as_epsg4326([record.portee.geom]) areas, _envelope, _area = get_as_epsg4326([record.geom]) points = [point.geom for point in record.arbre] points, _points_bbox, _points_area = get_as_epsg4326(points) # Léger dézoom pour afficher correctement les zones qui touchent la bbox envelope = [ _envelope[0] - 0.0001, _envelope[1] + 0.0001, _envelope[2] - 0.0001, _envelope[3] + 0.0001, ] m = MapRender(640, 480, envelope, True) m.add_bg() m.plot_geom(lines[0], None, None, color=(0, 0, 1, 0.3), bgcolor=(0, 0, 1, 0.1)) m.plot_geom(aires[0], None, None, color=(0, 0, 1, 0.3), bgcolor=(0, 0, 1, 0.1)) # Ajoute les points for point, rec in zip(points, record.arbre): m.plot_geom(point, rec.code , None, color=(1, 1, 1, 1), bgcolor=(1, 1, 1, 1)) m.plot_geom(areas[0], record.code, None, color=cls.COLOR, bgcolor=cls.BGCOLOR) data = m.render() cls.write([record], {'image_map': buffer(data)})
def situation_map_gen(cls, records): """Render the situation map""" CadPlot = Pool().get('cadastre.parcelle') Plot = Pool().get('forest.plot') for record in records: cad_plots = [plot.geom for plot in record.cad_plots] cad_plots, envelope, cad_area = get_as_epsg4326(cad_plots) plots = [plot.geom for plot in record.plots] plots, plot_bbox, _plots_area = get_as_epsg4326(plots) if envelope is None: continue # Compute the envelope if plot_bbox is not None: envelope = envelope_union(envelope, plot_bbox) # Include the geometry of the town in the bbox of the map if record.address is not None and record.address.my_city is not None: # Include the town from the address in the bbox town_geo = osr_geo_from_field(record.address.my_city.contour) dst = osr.SpatialReference() dst.SetWellKnownGeogCS("EPSG:4326") town_geo.TransformTo(dst) envelope = envelope_union(envelope, town_geo.GetEnvelope()) # Map title title = u'Plan de situation\n' title += u'Propriétaire: %s\n' % record.owner.name if record.address is not None \ and record.address.city is not None \ and record.address.my_city is not None: city = record.address.city dep = record.address.my_city.subdivision.parent.code.split('-')[1] title += u'Commune: %s (%s)\n' % (city, dep) title += u'Surface: %02i ha %02i a %02i ca\n\nLe ' % cls._area_to_a(cad_area) title += date.today().strftime('%02d/%02m/%Y') m = MapRender(1024, 768, envelope, True) m.add_bg() for plot in cad_plots: m.plot_geom(plot, None, u'Bois de La Forêt', color=CadPlot.COLOR, bgcolor=CadPlot.BGCOLOR) for plot in plots: m.plot_geom(plot, None, u'Parcelle forestière', linestyle='--', color=Plot.COLOR, bgcolor=Plot.BGCOLOR) data_nl = m.render() m.plot_legend() m.plot_compass() m.plot_scaling() cls._plot_logo(m) m.plot_title(title) data = m.render() cls.write([record], { 'situation_map': buffer(data), 'situation_map_nl': buffer(data_nl), })
def situation_map_gen(cls, records): """Render the situation map""" for record in records: # Récupère l'étendu de la zone de travaux areas, _envelope, _area = get_as_epsg4326([record.geom]) # Léger dézoom pour afficher correctement les points qui touchent la bbox envelope = [ _envelope[0] - 0.001, _envelope[1] + 0.001, _envelope[2] - 0.001, _envelope[3] + 0.001, ] if envelope is None: continue # Include the geometry of the town in the bbox of the map if record.commune is not None and record.commune.name is not None: # Include the town from the commune in the bbox town_geo = osr_geo_from_field(record.commune.contour) dst = osr.SpatialReference() dst.SetWellKnownGeogCS("EPSG:4326") town_geo.TransformTo(dst) envelope = envelope_union(envelope, town_geo.GetEnvelope()) # Map title title = u'Plan de situation communal\n' title += u'Propriétaire: %s\n' % record.owner.name if record.commune is not None \ and record.commune.name is not None \ and record.commune.name is not None: city = record.commune.name dep = record.commune.subdivision.parent.code.split('-')[1] title += u'Commune: %s (%s)\n' % (city, dep) title += u'Surface: %02i ha %02i a %02i ca\n\nLe ' % cls._area_to_a(_area) title += date.today().strftime('%02d/%02m/%Y') m = MapRender(1024, 768, envelope, True) # Ajoute le fond de carte m.add_bg() # Ajoute la zone de chantier m.plot_geom(areas[0], None, u'Site', color=cls.COLOR, bgcolor=cls.BGCOLOR) data_nl = m.render() m.plot_legend() m.plot_compass() m.plot_scaling() cls._plot_logo(m) m.plot_title(title) data = m.render() cls.write([record], { 'situation_map': buffer(data), })
def situation_closeup_map_gen(cls, records): """Render the situation map closeup""" CadPlot = Pool().get('cadastre.parcelle') Plot = Pool().get('forest.plot') for record in records: cad_plots = [plot.geom for plot in record.cad_plots] cad_plots, envelope, cad_area = get_as_epsg4326(cad_plots) plots = [plot.geom for plot in record.plots] plots, plot_bbox, _plots_area = get_as_epsg4326(plots) if envelope is None: continue # Compute the envelope if plot_bbox is not None: envelope = envelope_union(envelope, plot_bbox) # Map title title = u'Plan de situation\n' title += u'Propriétaire: %s\n' % record.owner.name if record.address is not None \ and record.address.city is not None \ and record.address.my_city is not None: city = record.address.city dep = record.address.my_city.subdivision.parent.code.split('-')[1] title += u'Commune: %s (%s)\n' % (city, dep) title += u'Surface: %02i ha %02i a %02i ca\n\nLe ' % cls._area_to_a(cad_area) title += date.today().strftime('%02d/%02m/%Y') m = MapRender(1024, 768, envelope, True) m.add_bg() cls._plot_misc_areas(m, record, False) for plot in cad_plots: m.plot_geom(plot, None, u'Bois de La Forêt', color=CadPlot.COLOR, bgcolor=CadPlot.BGCOLOR) for plot, rec in zip(plots, record.plots): m.plot_geom(plot, rec.short_name, None, linestyle='--', color=Plot.COLOR, bgcolor=Plot.BGCOLOR) cls._plot_misc_points(m, record, False) data_nl = m.render() m.plot_legend() m.plot_compass() m.plot_scaling() cls._plot_logo(m) m.plot_title(title) data = m.render() cls.write([record], { 'situation_closeup_map': buffer(data), 'situation_closeup_map_nl': buffer(data_nl), })
def situation_map_gen(cls, records): """Render the situation map""" for record in records: town, envelope_town, area_town = get_as_epsg4326([record.address.my_city.contour]) # Récupère l'étendu de la zone de garden section, envelope_section, area_section = get_as_epsg4326([record.section.geom]) lieudit, envelope_lieudit, area_lieudit = get_as_epsg4326([record.lieudit.geom]) parcelle, envelope_parcelle, area_parcelle = get_as_epsg4326([record.parcelle.geom]) # Léger dézoom pour afficher correctement les points qui touchent la bbox envelope = bbox_aspect(envelope_section, 640, 480) if envelope is None: continue # Map title title = u"Plan de situation du jardin\n" title += date.today().strftime("%02d/%02m/%Y") m = MapRender(1024, 768, envelope, True) # Ajoute le fond de carte m.add_bg() # Ajoute le contour de la ville m.plot_geom(town[0], None, u"Commune", color=(0, 0, 1, 1), bgcolor=(0, 0, 0, 0)) # Ajoute la section m.plot_geom(section[0], None, u"Section", color=(0, 0, 1, 0.3), bgcolor=(0, 0, 1, 0.3)) # Ajoute le lieud dit m.plot_geom(lieudit[0], None, u"Lieu-dit", color=(0, 1, 1, 0.3), bgcolor=(0, 1, 1, 0.3)) # Ajoute la pracelle m.plot_geom(parcelle[0], record.name, u"Parcelle", color=cls.COLOR, bgcolor=cls.BGCOLOR) data_nl = m.render() m.plot_legend() m.plot_compass() m.plot_scaling() m.plot_title(title) data = m.render() cls.write([record], {"situation_map": buffer(data)})
def situation_map_gen(cls, records): """Render the situation map""" for record in records: # Récupère l'étendu de la zone de secteur areas, _envelope, _area = get_as_epsg4326([record.geom]) # Léger dézoom pour afficher correctement les points qui touchent la bbox envelope = [ _envelope[0] - 0.001, _envelope[1] + 0.001, _envelope[2] - 0.001, _envelope[3] + 0.001, ] if envelope is None: continue # Map title title = u'Plan de situation de la zone urbanisée\n' title += u'Surface: %02i ha %02i a %02i ca\n\nLe ' % cls._area_to_a(_area) title += date.today().strftime('%02d/%02m/%Y') m = MapRender(1024, 768, envelope, True) # Ajoute le fond de carte m.add_bg() # Ajoute la zone urbanisée m.plot_geom(areas[0], None, u'Zone urbanisé', color=cls.COLOR, bgcolor=cls.BGCOLOR) data_nl = m.render() m.plot_legend() m.plot_compass() m.plot_scaling() cls._plot_logo(m) m.plot_title(title) data = m.render() cls.write([record], { 'situation_map': buffer(data), })
def generate(cls, records): for record in records: if record.name is None: continue areas, _envelope, _area = get_as_epsg4326([record.geom]) # Léger dézoom pour afficher correctement les points qui touchent la bbox envelope = [ _envelope[0] - 0.001, _envelope[1] + 0.001, _envelope[2] - 0.001, _envelope[3] + 0.001, ] m = MapRender(640, 480, envelope, True) m.add_bg() m.plot_geom(areas[0], None, None, color=cls.COLOR, bgcolor=cls.BGCOLOR) data = m.render() cls.write([record], {'image_map': buffer(data)})
def generate(cls, records): for record in records: if record.code is None: continue lines, envelope, _line = get_as_epsg4326([record.equipement.geom]) points, _envelope, _point = get_as_epsg4326([record.geom]) # Léger dézoom pour afficher correctement les zones qui touchent la bbox envelope = [ _envelope[0] - 0.001, _envelope[1] + 0.001, _envelope[2] - 0.001, _envelope[3] + 0.001, ] m = MapRender(640, 480, envelope, True) m.add_bg() m.plot_geom(lines[0], None, None, color=(0, 0, 1, 0.3), bgcolor=(0, 0, 1, 0.3)) m.plot_geom(points[0], record.code, None, color=cls.COLOR, bgcolor=cls.BGCOLOR) data = m.render() cls.write([record], {'image_map': buffer(data)})
def image_map_gen(cls, records): """Render the image map""" for record in records: # Récupère l'étendu de la zone de travaux areas, _envelope, _area = get_as_epsg4326([record.geom]) aires = [aire.geom for aire in record.misc_obj_poly] aires, _aires_bbox, _aires_area = get_as_epsg4326(aires) lignes = [ligne.geom for ligne in record.misc_obj_line] lignes, _lignes_bbox, _lignes_area = get_as_epsg4326(lignes) points = [point.geom for point in record.misc_obj_point] points, _points_bbox, _points_area = get_as_epsg4326(points) # Léger dézoom pour afficher correctement les points qui touchent la bbox envelope = [ _envelope[0] - 0.001, _envelope[1] + 0.001, _envelope[2] - 0.001, _envelope[3] + 0.001, ] if envelope is None: continue # Map title title = u'Plan de situation chantier\n' title += u'Propriétaire: %s\n' % record.owner.name if record.commune is not None \ and record.commune.name is not None \ and record.commune.name is not None: city = record.commune.name dep = record.commune.subdivision.parent.code.split('-')[1] title += u'Commune: %s (%s)\n' % (city, dep) title += u'Surface: %02i ha %02i a %02i ca\n\nLe ' % cls._area_to_a(_area) title += date.today().strftime('%02d/%02m/%Y') m = MapRender(1024, 768, envelope, True) # Ajoute le fond de carte m.add_bg() # Ajoute les polygones for aire, rec in zip(aires, record.misc_obj_poly): m.plot_geom(aire, rec.name, u'Zones', color=(1, 1, 1, 1), bgcolor=(0, 0, 1, 1)) # Ajoute les polylignes for ligne, rec in zip(lignes, record.misc_obj_line): m.plot_geom(ligne, rec.name, None, color=(1, 1, 1, 1), bgcolor=(1, 1, 1, 1)) # Ajoute les points for point, rec in zip(points, record.misc_obj_point): m.plot_geom(point, rec.name, None, color=(1, 1, 1, 1), bgcolor=(1, 1, 1, 1)) # Ajoute la zone de chantier m.plot_geom(areas[0], None, u'Chantier', color=cls.COLOR, bgcolor=cls.BGCOLOR) data_nl = m.render() m.plot_legend() m.plot_compass() m.plot_scaling() cls._plot_logo(m) m.plot_title(title) data = m.render() cls.write([record], { 'image_map': buffer(data), })
def situation_map_gen(cls, records): """Render the situation map""" for record in records: _envelope = None # Récupère les géométries de la zone d'information aires = [aire.geom for aire in record.pres_obj_poly] aires, _aires_bbox, _aires_area = get_as_epsg4326(aires) if _aires_bbox is None: _envelope = _envelope else: _envelope = envelope_union(_aires_bbox, _envelope) lignes = [ligne.geom for ligne in record.pres_obj_line] lignes, _lignes_bbox, _lignes_area = get_as_epsg4326(lignes) if _lignes_bbox is None: _envelope = _envelope else: _envelope = envelope_union(_lignes_bbox, _envelope) points = [point.geom for point in record.pres_obj_point] points, _points_bbox, _points_area = get_as_epsg4326(points) if _points_bbox is None: _envelope = _envelope else: _envelope = envelope_union(_points_bbox, _envelope) # Léger dézoom pour afficher correctement les points qui touchent la bbox envelope = [ _envelope[0] - 0.001, _envelope[1] + 0.001, _envelope[2] - 0.001, _envelope[3] + 0.001, ] # Map title title = u'Plan de situation\n' title += u'Surface: %02i ha %02i a %02i ca\n\nLe ' % cls._area_to_a(_aires_area) title += date.today().strftime('%02d/%02m/%Y') m = MapRender(1024, 768, envelope, True) # Ajoute le fond de carte m.add_bg() # Ajoute les polygones for aire, rec in zip(aires, record.pres_obj_poly): m.plot_geom(aire, rec.name, u'Zones', color=(1, 1, 1, 1), bgcolor=(0, 0, 1, 0.2)) # Ajoute les polylignes for ligne, rec in zip(lignes, record.pres_obj_line): m.plot_geom(ligne, rec.name, None, color=(1, 1, 1, 1), bgcolor=(1, 1, 1, 1)) # Ajoute les points for point, rec in zip(points, record.pres_obj_point): m.plot_geom(point, rec.name, None, color=(1, 1, 1, 1), bgcolor=(1, 1, 1, 1)) data_nl = m.render() m.plot_legend() m.plot_compass() m.plot_scaling() cls._plot_logo(m) m.plot_title(title) data = m.render() cls.write([record], { 'situation_map': buffer(data), })