def test_polygon_boundary():
"""Import XYZ polygons from file and test boundary function."""
mypoly = Polygons()
mypoly.from_file(PFILE, fformat="xyz")
boundary = mypoly.get_boundary()
tsetup.assert_almostequal(boundary[0], 460595.6036, 0.0001)
tsetup.assert_almostequal(boundary[4], 2025.952637, 0.0001)
tsetup.assert_almostequal(boundary[5], 2266.996338, 0.0001)
def test_import_export_polygons():
"""Import XYZ polygons from file. Modify, and export."""
mypoly = Polygons()
mypoly.from_file(PFILE, fformat="xyz")
z0 = mypoly.dataframe["Z_TVDSS"].values[0]
tsetup.assert_almostequal(z0, 2266.996338, 0.001)
mypoly.dataframe["Z_TVDSS"] += 100
mypoly.to_file(TMPD + "/polygon_export.xyz", fformat="xyz")
# reimport and check
mypoly2 = Polygons(TMPD + "/polygon_export.xyz")
tsetup.assert_almostequal(z0 + 100, mypoly2.dataframe["Z_TVDSS"].values[0], 0.001)
def test_import_zmap_and_xyz():
"""Import XYZ polygons on ZMAP and XYZ format from file"""
mypol2a = Polygons()
mypol2b = Polygons()
mypol2c = Polygons()
mypol2a.from_file(PFILE1A, fformat="zmap")
mypol2b.from_file(PFILE1B)
mypol2c.from_file(PFILE1C)
assert mypol2a.nrow == mypol2b.nrow
assert mypol2b.nrow == mypol2c.nrow
logger.info(mypol2a.nrow, mypol2b.nrow)
logger.info(mypol2a.dataframe)
logger.info(mypol2b.dataframe)
for col in ["X_UTME", "Y_UTMN", "Z_TVDSS", "POLY_ID"]:
status = np.allclose(
mypol2a.dataframe[col].values, mypol2b.dataframe[col].values
)
assert status is True
def test_import_zmap_and_xyz(testpath):
"""Import XYZ polygons on ZMAP and XYZ format from file"""
mypol2a = Polygons()
mypol2b = Polygons()
mypol2c = Polygons()
mypol2a.from_file(testpath / PFILE1A, fformat="zmap")
mypol2b.from_file(testpath / PFILE1B)
mypol2c.from_file(testpath / PFILE1C)
assert mypol2a.nrow == mypol2b.nrow
assert mypol2b.nrow == mypol2c.nrow
for col in ["X_UTME", "Y_UTMN", "Z_TVDSS", "POLY_ID"]:
assert np.allclose(mypol2a.dataframe[col].values,
mypol2b.dataframe[col].values)