From 7857dbe8a24145f9c3535e4b997d0e5b9c0e83c8 Mon Sep 17 00:00:00 2001 From: Jordan Gilbert Date: Wed, 23 Aug 2023 11:08:53 -0600 Subject: [PATCH] french proj xmls --- python/packages/rsxml/examples/arly_vbet.py | 244 ++++++++++++++++++ .../rsxml/examples/french_rscontext.py | 15 +- python/packages/rsxml/examples/validation.py | 3 +- 3 files changed, 256 insertions(+), 6 deletions(-) create mode 100644 python/packages/rsxml/examples/arly_vbet.py diff --git a/python/packages/rsxml/examples/arly_vbet.py b/python/packages/rsxml/examples/arly_vbet.py new file mode 100644 index 00000000..e829ae04 --- /dev/null +++ b/python/packages/rsxml/examples/arly_vbet.py @@ -0,0 +1,244 @@ +from rsxml.project_xml import ( + Project, + MetaData, + Meta, + ProjectBounds, + Coords, + BoundingBox, + Dataset, + Geopackage, + GeopackageLayer, + GeoPackageDatasetTypes, + Realization, + Log +) +from rsxml import Logger +from datetime import datetime +import sys +import traceback + +if __name__ == '__main__': + + log = Logger('Project') + + # Create a new Riverscapes Project from scratch + project = Project( + name='VBET for Arly', + proj_path='/mnt/c/Users/jordang/Documents/Riverscapes/data/VBET/Arly/project.rs.xml', + project_type='VBET', + description='VBET pour le bassin versant de larly en France', + bounds=ProjectBounds( + centroid=Coords(6.5235, 45.7618), + bounding_box=BoundingBox(6.2, 6.8, 45.6, 46), + filepath='/mnt/c/Users/jordang/Documents/Riverscapes/data/VBET/Arly/project_bounds.json', + ), + meta_data=MetaData(values=[Meta('Date de Création', '2023-07-19'), Meta('Bassin Versant', 'Arly')]), + realizations=[ + Realization( + xml_id='REALIZATION1', + name='VBET for Arly', + product_version='3.0.7', + date_created=datetime(2023, 7, 19), + description='projet personnalisé du VBET', + logs=[ + Log( + xml_id='LOGFILE', + name='Log', + path='vbet.log', + ) + ], + inputs=[ + Dataset( + xml_id='HILLSHADE', + name='DEM Hillshade', + path='inputs/dem_hillshade.tif', + ds_type='Raster' + ), + Dataset( + xml_id='DEM', + name='DEM', + path='inputs/dem.tif', + ds_type='Raster' + ), + Dataset( + xml_id='SLOPE_RASTER', + name='Slope Raster', + path='inputs/slope.tif', + ds_type='Raster' + ), + Geopackage( + xml_id='INPUTS', + name='Inputs', + path='inputs/inputs.gpkg', + layers=[ + GeopackageLayer( + lyr_name = 'flowlines', + name = 'NHD Flowlines', + ds_type = GeoPackageDatasetTypes.VECTOR + ), + GeopackageLayer( + lyr_name = 'channel_area_polygons', + name = 'Channel Area Polygons', + ds_type = GeoPackageDatasetTypes.VECTOR + ) + ] + ) + ], + intermediates=[ + Dataset( + xml_id='PITFILL', + name='TauDEM Pitfill', + path='intermediates/pitfill.tif', + ds_type='Raster' + ), + Dataset( + xml_id='DINFFLOWDIR_ANG', + name='TauDEM D-Inf Flow Directions', + path='intermediates/dinfflowdir_ang.tif', + ds_type='Raster' + ), + Dataset( + xml_id='DINFFLOWDIR_SLP', + name='TauDEM D-Inf Flow Directions Slope', + path='intermediates/dinfflowdir_slp.tif', + ds_type='Raster' + ), + Dataset( + xml_id='VBET_ZONES', + name='VBET LevelPath Zones', + path='intermediates/vbet_level_path_zones.tif', + ds_type='Raster' + ), + Dataset( + xml_id='ACTIVE_FP_ZONES', + name='Active Floodplain LevelPath Zones', + path='intermediates/active_fp_level_path_zones.tif', + ds_type='Raster' + ), + Dataset( + xml_id='INACTIVE_FP_ZONES', + name='Inactive Floodplain LevelPath Zones', + path='intermediates/inactive_fp_level_path_zones.tif', + ds_type='Raster' + ), + Dataset( + xml_id='HAND_RASTER', + name='HAND Raster', + path='intermediates/hand_composite.tif', + ds_type='Raster' + ), + Dataset( + xml_id='TRANSFORMED_HAND', + name='Transformed HAND Evidence', + path='intermediates/hand_transformed.tif', + ds_type='Raster' + ), + Dataset( + xml_id='TRANSFORMED_SLOPE', + name='Transformed Slope', + path='intermediates/slope_transformed.tif', + ds_type='Raster' + ), + Dataset( + xml_id='EVIDENCE_TOPO_INTERIOR', + name='Topo Evidence (Interior)', + path='intermediates/topographic_evidence_interior.tif', + ds_type='Raster' + ), + Dataset( + xml_id='HAND_RASTER_INTERIOR', + name='HAND Raster (Interior)', + path='intermediates/hand_composite_interior.tif', + ds_type='Raster' + ), + Dataset( + xml_id='TRANSFORMED_HAND_INTERIOR', + name='Transformed HAND Evidence (Interior)', + path='intermediates/hand_transformed_interior.tif', + ds_type='Raster' + ), + Dataset( + xml_id='TRANSFORMED_SLOPE_INTERIOR', + name='Transformed Slope (Interior)', + path='intermediates/slope_transformed_interior.tif', + ds_type='Raster' + ), + Geopackage( + xml_id='Intermediates', + name='Intermediates', + path='intermediates/intermediates.gpkg', + layers=[ + GeopackageLayer( + lyr_name = 'vbet_dgos', + name = 'VBET DGO Polygons', + ds_type = GeoPackageDatasetTypes.VECTOR + ) + ] + ) + ], + outputs=[ + Dataset( + xml_id='VBET_EVIDENCE', + name='VBET Evidence Raster', + path='outputs/vbet_evidence.tif', + ds_type='Raster' + ), + Geopackage( + xml_id='VBET_OUTPUTS', + name='VBET', + path='outputs/vbet.gpkg', + layers=[ + GeopackageLayer( + lyr_name = 'vbet_full', + name = 'VBET Full Extent', + ds_type = GeoPackageDatasetTypes.VECTOR + ), + GeopackageLayer( + lyr_name = 'active_valley_bottom', + name = 'VBET Inactive/Active Boundary', + ds_type = GeoPackageDatasetTypes.VECTOR + ), + GeopackageLayer( + lyr_name='vbet_channel_area', + name='VBET Channel Area', + ds_type=GeoPackageDatasetTypes.VECTOR + ), + GeopackageLayer( + lyr_name='active_floodplain', + name='Active Floodplain', + ds_type=GeoPackageDatasetTypes.VECTOR + ), + GeopackageLayer( + lyr_name='inactive_floodplain', + name='Inactive Floodplain', + ds_type=GeoPackageDatasetTypes.VECTOR + ), + GeopackageLayer( + lyr_name='floodplain', + name='Floodplain', + ds_type=GeoPackageDatasetTypes.VECTOR + ), + GeopackageLayer( + lyr_name='vbet_centerlines', + name='VBET Centerline', + ds_type=GeoPackageDatasetTypes.VECTOR + ), + GeopackageLayer( + lyr_name='vbet_igos', + name='Segmentation Points', + ds_type=GeoPackageDatasetTypes.VECTOR + ), + ] + ) + ] + ) + ] + ) + + #Write it to disk + try: + project.write() + except Exception as e: + traceback.print_exc(file=sys.stdout) + + log.info('Finished writing project file') \ No newline at end of file diff --git a/python/packages/rsxml/examples/french_rscontext.py b/python/packages/rsxml/examples/french_rscontext.py index 9885758d..64f3c227 100644 --- a/python/packages/rsxml/examples/french_rscontext.py +++ b/python/packages/rsxml/examples/french_rscontext.py @@ -23,7 +23,9 @@ Realization, ) from rsxml import Logger -from datetime import date +from datetime import datetime +import sys +import traceback if __name__ == '__main__': @@ -32,7 +34,7 @@ # Create a new Riverscapes Project from scratch project = Project( name='Riverscapes Context for Arly', - proj_path='project.rs.xml', + proj_path='/mnt/c/Users/jordang/Documents/Riverscapes/data/France/Isere_River/Arly/RSContext/project.rs.xml', project_type='RSContext', description='Riverscapes Context pour le bassin versant de larly en France', citation='', @@ -46,8 +48,8 @@ Realization( xml_id='REALIZATION1', name='Riverscapes Context for Arly', - product_version='Custom', - date_created=date(2023, 7, 18), + product_version='1.4.2', + date_created=datetime(2023, 7, 18), summary='', description='projet personnalisé du Riverscapes Context', # meta_data=MetaData(values=[Meta('Test', 'Test Value')]), @@ -225,6 +227,9 @@ ) # Write it to disk - project.write() + try: + project.write() + except Exception as e: + traceback.print_exc(file=sys.stdout) log.info('done') diff --git a/python/packages/rsxml/examples/validation.py b/python/packages/rsxml/examples/validation.py index 9457d578..0164f0cd 100644 --- a/python/packages/rsxml/examples/validation.py +++ b/python/packages/rsxml/examples/validation.py @@ -13,7 +13,8 @@ def main(): # ================================== # Open the xml file called data/xml_builder_test.xml in the tests directory # and read it into a string - xml_file = os.path.join(os.path.dirname(__file__), '..', 'tests', 'data', 'xml_builder_test.xml') + #xml_file = os.path.join(os.path.dirname(__file__), '..', 'tests', 'data', 'xml_builder_test.xml') + xml_file = '/mnt/c/Users/jordang/Documents/Riverscapes/data/France/Isere_River/Arly/RSContext/project.rs.xml' result, errors = validate_project_file(xml_file) log.info(f'Passed: {result}')