ras-commander 0.46.0__py3-none-any.whl → 0.48.0__py3-none-any.whl

ras_commander/HdfBndry.py

@@ -70,9 +70,9 @@ class HdfBndry:
                 return gpd.GeoDataFrame(
                     {
                         "bc_line_id": bc_line_ids,
-                        "name": names,
+                        "Name": names,
                         "mesh_name": mesh_names,
-                        "type": types,
+                        "Type": types,
                         "geometry": geoms,
                     },
                     geometry="geometry",
@@ -110,7 +110,7 @@ class HdfBndry:
                 )
                 geoms = HdfBndry._get_polylines(hdf_file, breaklines_path)
                 return gpd.GeoDataFrame(
-                    {"bl_id": bl_line_ids, "name": names, "geometry": geoms},
+                    {"bl_id": bl_line_ids, "Name": names, "geometry": geoms},
                     geometry="geometry",
                     crs=HdfUtils.projection(hdf_file),
                 )
@@ -158,7 +158,7 @@ class HdfBndry:
                             )
                         )
                 return gpd.GeoDataFrame(
-                    {"rr_id": rr_ids, "name": names, "geometry": geoms},
+                    {"rr_id": rr_ids, "Name": names, "geometry": geoms},
                     geometry="geometry",
                     crs=HdfUtils.projection(hdf_file),
                 )
@@ -242,9 +242,9 @@ class HdfBndry:
                 return gpd.GeoDataFrame(
                     {
                         "refln_id": refline_ids,
-                        "refln_name": names,
-                        "mesh_name": mesh_names,
-                        "type": types,
+                        "Name": names,
+                        "mesh-name": mesh_names,
+                        "Type": types,
                         "geometry": geoms,
                     },
                     geometry="geometry",
@@ -286,9 +286,9 @@ class HdfBndry:
                 return gpd.GeoDataFrame(
                     {
                         "refpt_id": range(attributes.shape[0]),
-                        "refpt_name": names,
+                        "Name": names,
                         "mesh_name": mesh_names,
-                        "cell_id": cell_id,
+                        "Cell Index": cell_id,
                         "geometry": list(map(Point, points)),
                     },
                     geometry="geometry",

ras_commander/HdfFluvialPluvial.py

@@ -3,7 +3,6 @@ import pandas as pd
 import geopandas as gpd
 import matplotlib.pyplot as plt
 from collections import defaultdict
-from rtree import index
 from shapely.geometry import LineString, MultiLineString
 from tqdm import tqdm
 
@@ -236,6 +235,11 @@ class HdfFluvialPluvial:
 
     @staticmethod
     def _process_cell_adjacencies(cell_polygons_gdf: gpd.GeoDataFrame) -> Tuple[Dict[int, List[int]], Dict[int, Dict[int, LineString]]]:
+        """
+        Process cell adjacencies and common edges using R-tree indexing.
+        """
+        # Install rtree using pip install rtree if not already installed
+        from rtree import index
         cell_adjacency = defaultdict(list)
         common_edges = defaultdict(dict)
         idx = index.Index()
@@ -314,4 +318,9 @@ class HdfFluvialPluvial:
         - GeoDataFrame with the specified projection.
         """
         gdf = gpd.GeoDataFrame(df, geometry='geometry', crs=projection)
+        
         return gdf
+
+
+
+

ras_commander/HdfMesh.py

@@ -64,7 +64,7 @@ class HdfMesh:
                     return list()
                 return list(
                     [
-                        HdfUtils.convert_ras_hdf_string(n)
+                        HdfUtils.convert_ras_hdf_string(n.decode('utf-8'))  # Decode as UTF-8
                         for n in hdf_file["Geometry/2D Flow Areas/Attributes"][()]["Name"]
                     ]
                 )
@@ -174,12 +174,11 @@ class HdfMesh:
         except Exception as e:
             logger.error(f"Error reading mesh cell polygons from {hdf_path}: {str(e)}")
             return GeoDataFrame()
-
     @staticmethod
     @standardize_input(file_type='plan_hdf')
     def mesh_cell_points(hdf_path: Path) -> GeoDataFrame:
         """
-        Return 2D flow mesh cell points.
+        Return 2D flow mesh cell center points.
 
         Parameters
         ----------
@@ -189,29 +188,29 @@ class HdfMesh:
         Returns
         -------
         GeoDataFrame
-            A GeoDataFrame containing the 2D flow mesh cell points.
+            A GeoDataFrame containing the 2D flow mesh cell center points.
         """
         try:
             with h5py.File(hdf_path, 'r') as hdf_file:
                 mesh_area_names = HdfMesh.mesh_area_names(hdf_path)
                 if not mesh_area_names:
                     return GeoDataFrame()
+                
                 pnt_dict = {"mesh_name": [], "cell_id": [], "geometry": []}
-                for i, mesh_name in enumerate(mesh_area_names):
-                    starting_row, count = hdf_file["Geometry/2D Flow Areas/Cell Info"][()][i]
-                    cell_pnt_coords = hdf_file["Geometry/2D Flow Areas/Cell Points"][()][
-                        starting_row : starting_row + count
-                    ]
-                    pnt_dict["mesh_name"] += [mesh_name] * cell_pnt_coords.shape[0]
-                    pnt_dict["cell_id"] += range(count)
+                for mesh_name in mesh_area_names:
+                    cell_center_coords = hdf_file[f"Geometry/2D Flow Areas/{mesh_name}/Cells Center Coordinate"][()]
+                    cell_count = len(cell_center_coords)
+                    
+                    pnt_dict["mesh_name"] += [mesh_name] * cell_count
+                    pnt_dict["cell_id"] += range(cell_count)
                     pnt_dict["geometry"] += list(
                         np.vectorize(lambda coords: Point(coords), signature="(n)->()")(
-                            cell_pnt_coords
+                            cell_center_coords
                         )
                     )
                 return GeoDataFrame(pnt_dict, geometry="geometry", crs=HdfUtils.projection(hdf_path))
         except Exception as e:
-            self.logger.error(f"Error reading mesh cell points from {hdf_path}: {str(e)}")
+            logger.error(f"Error reading mesh cell points from {hdf_path}: {str(e)}")
             return GeoDataFrame()
 
     @staticmethod
@@ -293,7 +292,7 @@ class HdfMesh:
                         try:
                             value = d2_flow_area[name][()]
                             if isinstance(value, bytes):
-                                value = value.decode('utf-8')
+                                value = value.decode('utf-8')  # Decode as UTF-8
                             result[name] = value
                         except Exception as e:
                             logger.warning(f"Error converting attribute '{name}': {str(e)}")
@@ -305,7 +304,6 @@ class HdfMesh:
             logger.error(f"Error reading 2D flow area attributes from {hdf_path}: {str(e)}")
             return {}
 
-
     @staticmethod
     @standardize_input(file_type='geom_hdf')
     def get_face_property_tables(hdf_path: Path) -> Dict[str, pd.DataFrame]:
@@ -340,10 +338,10 @@ class HdfMesh:
                         for z, area, wetted_perimeter, mannings_n in face_values:
                             face_data.append({
                                 'Face ID': face_id,
-                                'Z': z,
-                                'Area': area,
-                                'Wetted Perimeter': wetted_perimeter,
-                                "Manning's n": mannings_n
+                                'Z': str(z),
+                                'Area': str(area), 
+                                'Wetted Perimeter': str(wetted_perimeter),
+                                "Manning's n": str(mannings_n)
                             })
                     
                     result[mesh_name] = pd.DataFrame(face_data)
@@ -352,4 +350,4 @@ class HdfMesh:
 
         except Exception as e:
             logger.error(f"Error extracting face property tables from {hdf_path}: {str(e)}")
-            return {}
+            return {}

ras_commander/RasPrj.py

@@ -107,7 +107,7 @@ class RasPrj:
         # Add Geom_File to plan_df
         self.plan_df['Geom_File'] = self.plan_df.apply(lambda row: self._get_geom_file_for_plan(row['plan_number']), axis=1)
 
-    @log_call
+    
     def _get_geom_file_for_plan(self, plan_number):
         """
         Get the geometry file path for a given plan number.
@@ -133,7 +133,7 @@ class RasPrj:
             logger.error(f"Error reading plan file for geometry: {e}")
         return None
 
-    @log_call
+
     def _parse_plan_file(self, plan_file_path):
         """
         Parse a plan file and extract critical information.
@@ -193,9 +193,6 @@ class RasPrj:
         
         return plan_info
     
-
-    
-    @log_call
     def _get_prj_entries(self, entry_type):
         """
         Extract entries of a specific type from the HEC-RAS project file.
@@ -242,7 +239,6 @@ class RasPrj:
 
         return pd.DataFrame(entries)
 
-    @log_call
     def _parse_unsteady_file(self, unsteady_file_path):
         """
         Parse an unsteady flow file and extract critical information.
@@ -617,7 +613,6 @@ class RasPrj:
         
         return merged_df
 
-    @log_call
     def _parse_boundary_condition(self, block: str, unsteady_number: str, bc_number: int) -> Tuple[Dict, str]:
         lines = block.split('\n')
         bc_info = {

{ras_commander-0.46.0.dist-info → ras_commander-0.48.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: ras-commander
-Version: 0.46.0
+Version: 0.48.0
 Summary: A Python library for automating HEC-RAS operations
 Home-page: https://github.com/billk-FM/ras-commander
 Author: William M. Katzenmeyer
@@ -61,7 +61,7 @@ Create a virtual environment with conda or venv (ask ChatGPT if you need help)
 
 In your virtual environment, install ras-commander using pip:
 ```
-pip install h5py numpy pandas requests tqdm scipy
+pip install h5py numpy pandas requests tqdm scipy xarray geopandas matplotlib ras-commander ipython psutil shapely fiona pathlib rtree
 pip install --upgrade ras-commander
 ```
 

{ras_commander-0.46.0.dist-info → ras_commander-0.48.0.dist-info}/RECORD

@@ -1,8 +1,8 @@
 ras_commander/Decorators.py,sha256=i5AEQbe7JeI8Y3O_dQ5OO4Ab0KO5SiZTiysFBGxqTRU,4978
 ras_commander/HdfBase.py,sha256=HV5ccV9QH2lz4ZRYqK2d7_S833cTSUcostzxxSPb4O4,7129
-ras_commander/HdfBndry.py,sha256=LWaDMHeo0V_VOpx7D9q7W_0WvWsD6NQThiSLYPY8ApE,20761
-ras_commander/HdfFluvialPluvial.py,sha256=Ca0Yd8Vk7Em2xJz9sM5IXkvCiH5Qm0DKrJcookGYmgw,13765
-ras_commander/HdfMesh.py,sha256=1P_780k6CaTOcOTdzMnG8qzs6U7-_JyNDOw0sgQaY10,14648
+ras_commander/HdfBndry.py,sha256=29SAcn-LGf3QiNCeRmlT7JsGI94Dwl9NMm-zHNBuXlA,20752
+ras_commander/HdfFluvialPluvial.py,sha256=5OwYYEoLRJ4gfzi4kOBPO9_LVI_OyqW-Z8Jr7POOr34,13965
+ras_commander/HdfMesh.py,sha256=lNzEN7npA4KXUB9ZLyCMRk-EuL0oPD6J5MPZYamV7xU,14657
 ras_commander/HdfPipe.py,sha256=StL7D5hgdknzwS8VRSd9Rg_U9VOstepDvYvquLYCLOs,34348
 ras_commander/HdfPlan.py,sha256=DDnkMZ4Px4_4Netb7rLyxpDG1N-C742Yce4fDmNZAVo,6977
 ras_commander/HdfPump.py,sha256=LPQPjOL6dWpJBW2F9pQVzNyqQL1XvNvUmjtf_hDygV4,10486
@@ -18,13 +18,13 @@ ras_commander/RasExamples.py,sha256=FDit6FWzTQbv42Xx0gtSJB17HDeol_eykxV0qVpnwIk,
 ras_commander/RasGeo.py,sha256=x7eFvUQ_dPZxAVq2Gw_GIRzn-Cf7DLKZDTwk1fjrkf8,5395
 ras_commander/RasGpt.py,sha256=1Nitlf5QyeBQo8Q8I8pJKaTtQ0GdsBQA71LIKUAOycA,839
 ras_commander/RasPlan.py,sha256=fVuWQL6bnRCoHm-WiWv2mWZ8vTsDXn1BxcwUf03ma8s,51276
-ras_commander/RasPrj.py,sha256=ePLcPNcKVMgijvJKkJIWwQslQi2D6zq2UTW-8XCFmto,36159
+ras_commander/RasPrj.py,sha256=Ppmv3Yx_iN1ARKWC_dfSGVhdi4tHIuapDSdLbuAAGLU,36084
 ras_commander/RasToGo.py,sha256=UiUbvpjcAxAkWdSc-9Yywy7eOkGJaTHs24p4cjiL7co,588
 ras_commander/RasUnsteady.py,sha256=4j5Ga5hP8LepgR93uywjUjOEDhznLF7UnNtRmBS5d48,30779
 ras_commander/RasUtils.py,sha256=NBMxTHWHoTH2MJzqJ0y1_00fgKSS1GnNuEikwZ3Pqzs,34153
 ras_commander/__init__.py,sha256=-8mrN9qNbKiK0RfgCJasORkzuIKPHmY6RfYLzZqm2ok,1804
-ras_commander-0.46.0.dist-info/LICENSE,sha256=_pbd6qHnlsz1iQ-ozDW_49r86BZT6CRwO2iBtw0iN6M,457
-ras_commander-0.46.0.dist-info/METADATA,sha256=zmHXWeiuu3WRcnmcJJXiBLCvkTACli-Dm4uUpwJ9F9k,15935
-ras_commander-0.46.0.dist-info/WHEEL,sha256=OVMc5UfuAQiSplgO0_WdW7vXVGAt9Hdd6qtN4HotdyA,91
-ras_commander-0.46.0.dist-info/top_level.txt,sha256=i76S7eKLFC8doKcXDl3aiOr9RwT06G8adI6YuKbQDaA,14
-ras_commander-0.46.0.dist-info/RECORD,,
+ras_commander-0.48.0.dist-info/LICENSE,sha256=_pbd6qHnlsz1iQ-ozDW_49r86BZT6CRwO2iBtw0iN6M,457
+ras_commander-0.48.0.dist-info/METADATA,sha256=sEeBUJL0yGQez4hOPSqMr39IWLOzqMK90hLiaVFT45k,16020
+ras_commander-0.48.0.dist-info/WHEEL,sha256=P9jw-gEje8ByB7_hXoICnHtVCrEwMQh-630tKvQWehc,91
+ras_commander-0.48.0.dist-info/top_level.txt,sha256=i76S7eKLFC8doKcXDl3aiOr9RwT06G8adI6YuKbQDaA,14
+ras_commander-0.48.0.dist-info/RECORD,,

{ras_commander-0.46.0.dist-info → ras_commander-0.48.0.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (75.2.0)
+Generator: setuptools (75.3.0)
 Root-Is-Purelib: true
 Tag: py3-none-any