ax1.plot(row['geom_origin_pt3d_%s' % projstr].x, row['geom_origin_pt3d_%s' % projstr].y,
            '*', color='lightgray', markersize=12, label='origin pt')