|
| 1 | + |
| 2 | +from typing import NamedTuple |
| 3 | + |
| 4 | +from qgis.core import ( |
| 5 | + QgsTask, |
| 6 | + QgsVectorLayer |
| 7 | +) |
| 8 | + |
| 9 | +from ... import CanceledError |
| 10 | +from ...utils import ( |
| 11 | + LayerReader, |
| 12 | + SqlAccess, |
| 13 | + tr |
| 14 | +) |
| 15 | +from ._debug import debug_thread |
| 16 | + |
| 17 | + |
| 18 | +class Row(NamedTuple): |
| 19 | + fid: int |
| 20 | + geoid: str |
| 21 | + district: int |
| 22 | + geometry: str = None |
| 23 | + |
| 24 | + |
| 25 | +class AutoAssignUnassignedUnits(SqlAccess, QgsTask): |
| 26 | + def __init__(self, assignments: QgsVectorLayer, distField: str): |
| 27 | + super().__init__(tr("Auto assign unassigned units")) |
| 28 | + self.assignments = assignments |
| 29 | + self.distField = distField |
| 30 | + self.update: dict[int, int] = None |
| 31 | + self.indeterminate: list[Row] = None |
| 32 | + self.retry: list[Row] = None |
| 33 | + self.exception: Exception = None |
| 34 | + |
| 35 | + def run(self): |
| 36 | + debug_thread() |
| 37 | + |
| 38 | + try: |
| 39 | + reader = LayerReader(self.assignments, self) |
| 40 | + missing = reader.read_layer(columns=["fid", "geoid", self.distField], |
| 41 | + filt={self.distField: 0}, read_geometry=True, fid_as_index=True) \ |
| 42 | + .reset_index() \ |
| 43 | + .rename(columns={self.distField: "district"}) |
| 44 | + |
| 45 | + if len(missing) == self.assignments.featureCount(): |
| 46 | + # No units are assigned -- don't waste our time |
| 47 | + raise RuntimeError("Can't infer districts for unassigned units when no units are assigned") |
| 48 | + |
| 49 | + # find polygons that are adjacent by more than a point |
| 50 | + sql = f"""SELECT fid, geoid, {self.distField} AS district FROM assignments |
| 51 | + WHERE ST_relate(geometry, GeomFromText(:geometry), 'F***1****') |
| 52 | + AND fid IN ( |
| 53 | + SELECT id FROM rtree_assignments_geometry r |
| 54 | + WHERE r.minx < st_maxx(GeomFromText(:geometry)) |
| 55 | + AND r.maxx >= st_minx(GeomFromText(:geometry)) |
| 56 | + AND r.miny < st_maxy(GeomFromText(:geometry)) |
| 57 | + AND r.maxy >= st_miny(GeomFromText(:geometry)) |
| 58 | + )""" |
| 59 | + |
| 60 | + with self._connectSqlOgrSqlite(self.assignments.dataProvider()) as db: |
| 61 | + db.row_factory = lambda c, r: Row(*r) |
| 62 | + update: dict[int, int] = {} |
| 63 | + retry: list[Row] = [] |
| 64 | + indeterminate: list[Row] = [] |
| 65 | + count = 0 |
| 66 | + total = len(missing) |
| 67 | + for g in missing.to_wkt().itertuples(index=False, name="Row"): |
| 68 | + neighbors: list[Row] = db.execute(sql, (g.geometry,)).fetchall() |
| 69 | + dists = set(r.district for r in neighbors) |
| 70 | + |
| 71 | + # is the unassigned unit surrounded by units from the same district or |
| 72 | + # unassigned units (but not entirely by unassigned units) |
| 73 | + if len(dists) == 1 and 0 in dists: |
| 74 | + retry.append(g) |
| 75 | + elif len(dists) == 1 or (len(dists) == 2 and 0 in dists): |
| 76 | + newdist = max(dists) |
| 77 | + if g.fid in update: |
| 78 | + print("oops") |
| 79 | + update[g.fid] = newdist |
| 80 | + else: |
| 81 | + # multiple adjacent districts |
| 82 | + indeterminate.append(g) |
| 83 | + |
| 84 | + if self.isCanceled(): |
| 85 | + raise CanceledError() |
| 86 | + |
| 87 | + count += 1 |
| 88 | + self.setProgress(count/total) |
| 89 | + |
| 90 | + # TODO: there's probably a better way to to find the surrounding assigned units than continually looping |
| 91 | + while retry: |
| 92 | + retry_count = len(retry) |
| 93 | + newretry: list[Row] = [] |
| 94 | + for g in retry: |
| 95 | + neighbors: list[Row] = db.execute(sql, (g.geometry,)).fetchall() |
| 96 | + dists = set(update.get(r.fid, r.district) for r in neighbors) |
| 97 | + |
| 98 | + # is the unassigned unit surrounded by units from the same district or |
| 99 | + # unassigned units (but not entirely by unassigned units) |
| 100 | + if len(dists) == 1 and 0 in dists: |
| 101 | + newretry.append(g) |
| 102 | + elif len(dists) == 1 or (len(dists) == 2 and 0 in dists): |
| 103 | + newdist = max(dists) |
| 104 | + if g.fid in update: |
| 105 | + print("oops") |
| 106 | + update[g.fid] = newdist |
| 107 | + else: |
| 108 | + # multiple adjacent districts |
| 109 | + indeterminate.append(g) |
| 110 | + |
| 111 | + if self.isCanceled(): |
| 112 | + raise CanceledError() |
| 113 | + |
| 114 | + retry = newretry |
| 115 | + if len(retry) == retry_count: |
| 116 | + # if we were unable to assign any of the as-yet unassigned units, give up |
| 117 | + break |
| 118 | + |
| 119 | + self.update = update |
| 120 | + self.indeterminate = indeterminate |
| 121 | + self.retry = retry |
| 122 | + except CanceledError: |
| 123 | + return False |
| 124 | + except Exception as e: # pylint: disable=broad-except |
| 125 | + self.exception = e |
| 126 | + return False |
| 127 | + |
| 128 | + return True |
| 129 | + |
| 130 | + def finished(self, result): |
| 131 | + if not result: |
| 132 | + return |
| 133 | + |
| 134 | + if self.update: |
| 135 | + i = self.assignments.fields().indexFromName(self.distField) |
| 136 | + self.assignments.startEditing() |
| 137 | + self.assignments.beginEditCommand("Auto-assign") |
| 138 | + for fid, dist in self.update.items(): |
| 139 | + self.assignments.changeAttributeValue(fid, i, dist, 0) |
| 140 | + self.assignments.endEditCommand() |
0 commit comments