|
| 1 | +import geopandas as gpd |
| 2 | +import pandas as pd |
| 3 | +from shapely.geometry import Point |
| 4 | + |
| 5 | +gdf_tata_zout = gpd.GeoDataFrame( |
| 6 | + pd.DataFrame( |
| 7 | + data=[ |
| 8 | + ["A1", -180.0, -120.0, 1.0, 1.0, Point(0, 1)], |
| 9 | + ["A2", -180.0, -120.0, 1.0, 1.0, Point(0, 1)], |
| 10 | + ], |
| 11 | + columns=["Name", "botm", "top", "Q_Mm3/a", "Q_m3/d", "geometry"], |
| 12 | + ), |
| 13 | + crs="EPSG:28992", |
| 14 | +) |
| 15 | + |
| 16 | +# retreives top and botm from the model |
| 17 | +gdf_tata_zoet = gpd.GeoDataFrame( |
| 18 | + pd.DataFrame( |
| 19 | + data=[ |
| 20 | + ["B1", 1.0, 1.0, Point(0, 1)], |
| 21 | + ["B2", 1.0, 1.0, Point(0, 1)], |
| 22 | + ], |
| 23 | + columns=["Name", "Q_Mm3/a", "Q_m3/d", "geometry"], |
| 24 | + ), |
| 25 | + crs="EPSG:28992", |
| 26 | +) |
| 27 | + |
| 28 | +gdf_tata_zout.to_file("tata_zoutwaterbronnen.geojson", driver="GeoJSON") |
| 29 | +gdf_tata_zoet.to_file("tata_zoetwaterbronnen.geojson", driver="GeoJSON") |
| 30 | + |
| 31 | +zoet_fp = "/workspaces/NHFLO/models/modelscripts/09pwnmodel2/data/zoetwaterbronnen_tata.shp" |
| 32 | +zout_fp = "/workspaces/NHFLO/models/modelscripts/09pwnmodel2/data/zoutwaterbronnen_tata.shp" |
| 33 | + |
| 34 | +zout = gpd.read_file(zout_fp)[["Name", "botm", "top", "Q_Mm3/a", "Q_m3/d", "geometry"]] |
| 35 | +zoet = gpd.read_file(zoet_fp)[["Name", "Q_Mm3/a", "Q_m3/d", "geometry"]] |
| 36 | + |
| 37 | +zoetzout = pd.concat((zoet, zout)) |
| 38 | +print(f"Extent: {zoetzout.total_bounds}") |
| 39 | + |
| 40 | +zout.to_file("tata_zoutwaterbronnen_unaltered.geojson", driver="GeoJSON") |
| 41 | +zoet.to_file("tata_zoetwaterbronnen_unaltered.geojson", driver="GeoJSON") |
0 commit comments