def query(args):
eal = EAlGIS()
from dataexpr import DataExpression
iters = []
data = {}
for expr in args.equation:
ti = eal.get_table_info(args.geometry)
geometry_source = ti.geometry_source
it = iter(
DataExpression("CLI",
geometry_source,
expr,
"",
3112,
include_geometry=False,
order_by_gid=True).get_query().yield_per(1))
data[expr] = {}
for gid, v in it:
data[expr][gid] = v
all_gids = set()
for expr in data:
all_gids = all_gids.union(set(data[expr].keys()))
w = csv.writer(sys.stdout)
w.writerow(['gid'] + args.equation)
for gid in sorted(all_gids):
row = [gid] + [data[t].get(gid) for t in args.equation]
w.writerow(row)
def georelate(args):
eal = EAlGIS()
left_source = eal.get_table_info(args.geom_left).geometry_source
right_source = eal.get_table_info(args.geom_right).geometry_source
from .georelate import build_relations
build_relations(left_source, right_source)
def georelate(args):
eal = EAlGIS()
left_source = eal.get_table_info(args.geom_left).geometry_source
right_source = eal.get_table_info(args.geom_right).geometry_source
from .georelate import build_relations
build_relations(left_source, right_source)
