#!/usr/bin/python3 """ Convert PostGIS geometries to an image. To scale. Accepts a few geometry fine-tuning parameters. """ import argparse import geopandas import psycopg2 from matplotlib import rc import matplotlib.pyplot as plt CMAP = 'tab20c' # 'Set3' # is nice too PSQL_CREDS = "host=127.0.0.1 dbname=osm user=osm password=osm" COLORS = { 'black': '#000000', 'green': '#1b9e77', 'orange': '#d95f02', 'purple': '#7570b3', } # see `NOTICE` in the LaTeX document; this is the width of the main text block. TEXTWIDTH_CM = 12.12364 QUADRANTS = {'tl':1, 'tr':2, 'br':3, 'bl':4} def color(string): return COLORS[string if string else 'black'] def inch(cm): return cm / 2.54 def parse_args(): kwcolor = {'type': color, 'default': 'black'} parser = argparse.ArgumentParser(description=__doc__) parser.add_argument('--g1-select') parser.add_argument('--g1-linestyle') parser.add_argument('--g1-label') parser.add_argument('--g1-color', **kwcolor) parser.add_argument('--g2-select') parser.add_argument('--g2-linestyle') parser.add_argument('--g2-label') parser.add_argument('--g2-color', **kwcolor) parser.add_argument('--g3-select') parser.add_argument('--g3-linestyle') parser.add_argument('--g3-label') parser.add_argument('--g3-color', **kwcolor) parser.add_argument('--legend', help="Legend location, following matplotlib rules", default='best') parser.add_argument('--widthdiv', default=1, type=float, help="Divide the width by this number " "(useful when two images are laid horizontally " "in the resulting file") parser.add_argument('--quadrant', choices=QUADRANTS.keys(), help="Image is comprised of 4 quadrants. This variable, " "when non-empty, will clip and return the requested quadrant") parser.add_argument('--outfile', metavar='', help="If unset, displayed on the screen") return parser.parse_args() def read_layer(select, width, maybe_quadrant): if not select: return way = "way" if maybe_quadrant: way = "wm_quadrant(way, {})".format(QUADRANTS[maybe_quadrant]) conn = psycopg2.connect(PSQL_CREDS) sql = "SELECT {way} as way1 FROM {select}".format(way=way, select=select) return geopandas.read_postgis(sql, con=conn, geom_col='way1') def plot_args(geom, color, maybe_linestyle, maybe_label): if geom is None: return # polygons either have fillings or lines if geom.geom_type[0] == 'Polygon': if maybe_linestyle: return { 'edgecolor': 'black', 'linestyle': maybe_linestyle, 'color': (0, 0, 0, 0), } else: return {'cmap': CMAP, 'alpha': .25} r = {'color': color} if maybe_linestyle == 'invisible': r['color'] = (0, 0, 0, 0) elif maybe_linestyle: r['linestyle'] = maybe_linestyle if maybe_label: r['label'] = '\\normalfont %s' % maybe_label return r def main(): args = parse_args() width = TEXTWIDTH_CM / args.widthdiv g1 = read_layer(args.g1_select, width, args.quadrant) g2 = read_layer(args.g2_select, width, args.quadrant) g3 = read_layer(args.g3_select, width, args.quadrant) c1 = plot_args(g1, args.g1_color, args.g1_linestyle, args.g1_label) c2 = plot_args(g2, args.g2_color, args.g2_linestyle, args.g2_label) c3 = plot_args(g3, args.g3_color, args.g3_linestyle, args.g3_label) rc('text', usetex=True) rc('text.latex', preamble='\\usepackage{numprint}\n') fig, ax = plt.subplots(constrained_layout=True) fig.set_figwidth(inch(width)) g1 is not None and g1.plot(ax=ax, linewidth=.75, **c1) g2 is not None and g2.plot(ax=ax, linewidth=.75, **c2) g3 is not None and g3.plot(ax=ax, linewidth=.75, **c3) ax.legend(loc=args.legend, frameon=False) ax.axis('off') ax.margins(0, 0) if args.outfile: fig.savefig(args.outfile, bbox_inches='tight', dpi=600) else: plt.show() if __name__ == '__main__': main()