Skip to content

Commit

Permalink
french proj xmls
Browse files Browse the repository at this point in the history
  • Loading branch information
@jtgilbert
jtgilbert committed Aug 23, 2023
1 parent 251133b commit 7857dbe
Show file tree
Hide file tree
Showing 3 changed files with 256 additions and 6 deletions.
244 changes: 244 additions & 0 deletions python/packages/rsxml/examples/arly_vbet.py
View file
Original file line number Diff line number Diff line change
@@ -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')
15 changes: 10 additions & 5 deletions python/packages/rsxml/examples/french_rscontext.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -23,7 +23,9 @@
Realization,
)
from rsxml import Logger
from datetime import date
from datetime import datetime
import sys
import traceback

if __name__ == '__main__':

Expand All @@ -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='',
Expand All @@ -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')]),
Expand Down Expand Up @@ -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')
3 changes: 2 additions & 1 deletion python/packages/rsxml/examples/validation.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -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}')
Expand Down

0 comments on commit 7857dbe

Please sign in to comment.