def test_luk_punkt(
firedb: FireDb,
punkt: Punkt,
sagsevent: Sagsevent,
observationer: List[Observation],
koordinat: Koordinat,
punktinformationtype: PunktInformationType,
):
# byg et punkt der har tilknyttet geometri, koordinat,
# punktinfo og observationer
geometri = GeometriObjekt(punkt=punkt,
geometri=Point([10, 55]),
sagsevent=sagsevent)
firedb.session.add(geometri)
observationer[0].opstillingspunkt = punkt
observationer[1].sigtepunkt = punkt
koordinat.punkt = punkt
punkt.punktinformationer = [
PunktInformation(infotype=punktinformationtype, sagsevent=sagsevent)
]
firedb.session.commit()
firedb.luk_punkt(punkt, sagsevent)
assert punkt.registreringtil is not None
assert punkt.sagsevent.eventtype == EventType.PUNKT_NEDLAGT
assert punkt.sagseventtilid == sagsevent.id
assert geometri.registreringtil is not None
assert geometri.sagsevent.eventtype == EventType.PUNKT_NEDLAGT
for koordinat in punkt.koordinater:
assert koordinat.registreringtil is not None
assert koordinat.sagsevent.eventtype == EventType.PUNKT_NEDLAGT
assert koordinat.sagseventtilid == sagsevent.id
for punktinfo in punkt.punktinformationer:
assert punktinfo.registreringtil is not None
assert punktinfo.sagsevent.eventtype == EventType.PUNKT_NEDLAGT
assert punktinfo.sagseventtilid == sagsevent.id
for observation in chain(punkt.observationer_fra, punkt.observationer_til):
assert observation.registreringtil is not None
assert observation.sagsevent.eventtype == EventType.PUNKT_NEDLAGT
assert observation.sagseventtilid == sagsevent.id
with pytest.raises(TypeError):
firedb.luk_punkt(999)
def read(self, sags_id):
sag = self.fireDb.hent_sag(sags_id)
namespace = "{http://www.gnu.org/software/gama/gama-local-adjustment}"
tree = ET.parse(self.input_stream)
root = tree.getroot()
beregning = Beregning()
srid = self.fireDb.hent_srid("EPSG:5799")
adjusted_element = root.find(namespace +
"coordinates").find(namespace +
"adjusted")
cov_mat_values = (
root.find(namespace +
"coordinates").find(namespace +
"cov-mat").findall(namespace +
"flt"))
original_index_indicies = (
root.find(namespace +
"coordinates").find(namespace +
"original-index").findall(namespace +
"ind"))
for idx, point in enumerate(adjusted_element.iter(namespace +
"point")):
# Read values from the point
z = point.find(namespace + "z").text
point_id = point.find(namespace + "id").text
# Read the correct entry in cov_mat_values
cov_mat_index = int(original_index_indicies[idx].text) - 1
cov_mat_element = cov_mat_values[cov_mat_index]
# Read value as float
cov_mat_value = float(cov_mat_element.text)
# .. and tale sqrt to find std_dev
std_dev = math.sqrt(cov_mat_value)
koordinat = Koordinat()
koordinat.srid = srid
koordinat.z = z
koordinat.sz = std_dev
koordinat.punkt = self.fireDb.hent_punkt(point_id)
beregning.koordinater.append(koordinat)
observation_id_list = []
observations_element = root.find(namespace + "observations")
for idx, diff in enumerate(
observations_element.iter(namespace + "height-diff")):
observationId = diff.get("extern")
observation_id_list.append(observationId)
observation_list = self.fireDb.hent_observationer(observation_id_list)
beregning.observationer.extend(observation_list)
self.fireDb.indset_sagsevent(
Sagsevent(
sag=sag,
eventtype=EventType.KOORDINAT_BEREGNET,
sagseventinfos=[
SagseventInfo(beskrivelse="Beregning indlæst via GNU Gama")
],
beregninger=[beregning],
koordinater=beregning.koordinater,
))
def read(self, sags_id):
sag = self.fireDb.hent_sag(sags_id)
namespace = "{http://www.gnu.org/software/gama/gama-local-adjustment}"
tree = ET.parse(self.input_stream)
root = tree.getroot()
# In the description
description_element = root.find(namespace + "description")
description = description_element.text
# .. find all obervation ids
# observation_ids_start = description.find("{observation_ids}") + len("{observation_ids} :")
# bservation_ids_end = description.find("{/observation_ids}")
# observation_ids = description[observation_ids_start:observation_ids_end]
# observation_id_list= ast.literal_eval(observation_ids)
# ... and fetch those observations
beregning = Beregning()
srid = self.fireDb.hent_srid("EPSG:5799")
adjusted_element = root.find(namespace +
"coordinates").find(namespace +
"adjusted")
cov_mat_values = (
root.find(namespace +
"coordinates").find(namespace +
"cov-mat").findall(namespace +
"flt"))
original_index_indicies = (
root.find(namespace +
"coordinates").find(namespace +
"original-index").findall(namespace +
"ind"))
for idx, point in enumerate(adjusted_element.iter(namespace +
"point")):
# Read values from the point
z = point.find(namespace + "z").text
point_id = point.find(namespace + "id").text
# Read the correct entry in cov_mat_values
cov_mat_index = int(original_index_indicies[idx].text) - 1
cov_mat_element = cov_mat_values[cov_mat_index]
# Read value as float
cov_mat_value = float(cov_mat_element.text)
# .. and tale sqrt to find std_dev
std_dev = math.sqrt(cov_mat_value)
koordinat = Koordinat()
koordinat.srid = srid
koordinat.z = z
koordinat.sz = std_dev
koordinat.transformeret = "false"
koordinat.punkt = self.fireDb.hent_punkt(point_id)
beregning.koordinater.append(koordinat)
observation_id_list = []
observations_element = root.find(namespace + "observations")
for idx, diff in enumerate(
observations_element.iter(namespace + "height-diff")):
observationId = diff.get("extern")
observation_id_list.append(observationId)
observation_list = self.fireDb.hent_observationer(observation_id_list)
beregning.observationer.extend(observation_list)
self.fireDb.indset_beregning(Sagsevent(sag=sag), beregning)
