fimeval 0.1.54__tar.gz → 0.1.56__tar.gz

{fimeval-0.1.54 → fimeval-0.1.56}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: fimeval
-Version: 0.1.54
+Version: 0.1.56
 Summary: A Framework for Automatic Evaluation of Flood Inundation Mapping Predictions Evaluation
 Author: Surface Dynamics Modeling Lab
 Author-email: Supath Dhital <sdhital@crimson.ua.edu>, Dipshika Devi <ddevi@ua.edu>
@@ -683,6 +683,7 @@ Requires-Dist: notebook<8.0.0,>=7.3.2
 Requires-Dist: boto3<2.0.0,>=1.36.16
 Requires-Dist: geemap
 Requires-Dist: uv
+Requires-Dist: seaborn
 Provides-Extra: dev
 Requires-Dist: pytest; extra == "dev"
 Requires-Dist: black; extra == "dev"

{fimeval-0.1.54 → fimeval-0.1.56}/pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name = "fimeval"
-version = "0.1.54"
+version = "0.1.56"
 description = "A Framework for Automatic Evaluation of Flood Inundation Mapping Predictions Evaluation"
 readme = "README.md"
 requires-python = ">=3.10"
@@ -30,7 +30,8 @@ dependencies = [
   "notebook>=7.3.2,<8.0.0",
   "boto3>=1.36.16,<2.0.0",
   "geemap",
-  "uv"
+  "uv",
+  "seaborn"
 ]
 
 [project.optional-dependencies]

{fimeval-0.1.54 → fimeval-0.1.56}/src/fimeval/BuildingFootprint/evaluationwithBF.py

@@ -6,6 +6,9 @@ import pandas as pd
 from pathlib import Path
 from plotly.subplots import make_subplots
 import plotly.graph_objects as go
+import seaborn as sns
+import matplotlib.pyplot as plt
+import matplotlib.gridspec as gridspec
 
 
 def Changeintogpkg(input_path, output_dir, layer_name):
@@ -20,6 +23,7 @@ def Changeintogpkg(input_path, output_dir, layer_name):
         gdf.to_file(output_gpkg, driver="GPKG")
         return output_gpkg
 
+
 def GetFloodedBuildingCountInfo(
     building_fp_path,
     study_area_path,
@@ -57,7 +61,6 @@ def GetFloodedBuildingCountInfo(
         "True Positive": 0,
     }
 
-    # Count centroids in the contingency map
     def count_centroids_in_contingency(raster_path):
         with rasterio.open(raster_path) as src:
             raster_data = src.read(1)
@@ -74,10 +77,12 @@ def GetFloodedBuildingCountInfo(
 
     count_centroids_in_contingency(contingency_map)
 
-    # Calculate Candidate and Benchmark counts from the contingency map counts
-    centroid_counts["Candidate"] = centroid_counts["True Positive"] + centroid_counts["False Positive"]
-    centroid_counts["Benchmark"] = centroid_counts["True Positive"] + centroid_counts["False Negative"]
-
+    centroid_counts["Candidate"] = (
+        centroid_counts["True Positive"] + centroid_counts["False Positive"]
+    )
+    centroid_counts["Benchmark"] = (
+        centroid_counts["True Positive"] + centroid_counts["False Negative"]
+    )
 
     total_buildings = len(clipped_buildings)
     percentages = {
@@ -92,15 +97,13 @@ def GetFloodedBuildingCountInfo(
     CSI = TP / (TP + FP + FN) if (TP + FP + FN) > 0 else 0
     FAR = FP / (TP + FP) if (TP + FP) > 0 else 0
     POD = TP / (TP + FN) if (TP + FN) > 0 else 0
-
     if centroid_counts["Benchmark"] > 0:
         BDR = (
-            (centroid_counts["Candidate"] - centroid_counts["Benchmark"])
-            / centroid_counts["Benchmark"]
-        )
+            centroid_counts["Candidate"] - centroid_counts["Benchmark"]
+        ) / centroid_counts["Benchmark"]
     else:
-            BDR = 0
-        
+        BDR = 0
+
     counts_data = {
         "Category": [
             "Candidate",
@@ -125,13 +128,14 @@ def GetFloodedBuildingCountInfo(
             f"{BDR:.3f}",
         ],
     }
-
     counts_df = pd.DataFrame(counts_data)
     csv_file_path = os.path.join(
         save_dir, "EvaluationMetrics", f"BuildingCounts_{basename}.csv"
     )
+    os.makedirs(os.path.dirname(csv_file_path), exist_ok=True)
     counts_df.to_csv(csv_file_path, index=False)
 
+    # Plotly interactive visualization only
     third_raster_labels = ["False Positive", "False Negative", "True Positive"]
     third_raster_counts = [
         centroid_counts["False Positive"],
@@ -162,7 +166,6 @@ def GetFloodedBuildingCountInfo(
         row=1,
         col=1,
     )
-
     fig.add_trace(
         go.Bar(
             x=["Benchmark"],
@@ -178,17 +181,17 @@ def GetFloodedBuildingCountInfo(
         col=1,
     )
 
-    for i in range(len(third_raster_labels)):
+    for i, label in enumerate(third_raster_labels):
         fig.add_trace(
             go.Bar(
-                x=[third_raster_labels[i]],
+                x=[label],
                 y=[third_raster_counts[i]],
                 text=[f"{third_raster_counts[i]}"],
                 textposition="auto",
                 marker_color=["#ff5733", "#ffc300", "#28a745"][i],
                 marker_line_color="black",
                 marker_line_width=1,
-                name=f"{third_raster_labels[i]} ({percentages[third_raster_labels[i]]:.2f}%)",
+                name=f"{label} ({percentages[label]:.2f}%)",
             ),
             row=1,
             col=2,
@@ -200,19 +203,94 @@ def GetFloodedBuildingCountInfo(
         yaxis_title="Flooded Building Counts",
         width=1100,
         height=400,
-        plot_bgcolor="rgba(0, 0, 0, 0)",
-        paper_bgcolor="rgba(0, 0, 0, 0)",
+        plot_bgcolor="rgba(0,0,0,0)",
+        paper_bgcolor="rgba(0,0,0,0)",
         showlegend=True,
         font=dict(family="Arial", size=18, color="black"),
     )
+    fig.show()
+
+    # Seaborn for static PNG
+    df_left = pd.DataFrame(
+        {
+            "Category": ["Candidate", "Benchmark"],
+            "Count": [centroid_counts["Candidate"], centroid_counts["Benchmark"]],
+        }
+    )
+    df_right = pd.DataFrame(
+        {
+            "Category": third_raster_labels,
+            "Count": third_raster_counts,
+        }
+    )
 
+    sns.set_theme(style="whitegrid")
+
+    fig_sb = plt.figure(figsize=(10, 3), constrained_layout=True)
+    gs = gridspec.GridSpec(1, 3, figure=fig_sb, width_ratios=[1, 1, 0.4])
+
+    ax0 = fig_sb.add_subplot(gs[0, 0])
+    ax1 = fig_sb.add_subplot(gs[0, 1])
+    ax_leg = fig_sb.add_subplot(gs[0, 2])
+    ax_leg.axis("off")
+
+    def style_axes(ax, title_text, xlab, show_ylabel: bool):
+        ax.set_title(title_text, fontsize=14, pad=15)
+        ax.set_xlabel(xlab, fontsize=13, color="black")
+        if show_ylabel:
+            ax.set_ylabel("Flooded Building Counts", fontsize=13, color="black")
+        else:
+            ax.set_ylabel("")
+
+        for spine in ("left", "bottom"):
+            ax.spines[spine].set_linewidth(1.5)
+            ax.spines[spine].set_color("black")
+
+        sns.despine(ax=ax, right=True, top=True)
+        ax.tick_params(axis="x", labelsize=11, colors="black")
+        ax.tick_params(axis="y", labelsize=11, colors="black")
+
+    # Left panel
+    colors_left = ["#1c83eb", "#a4490e"]
+    sns.barplot(data=df_left, x="Category", y="Count", ax=ax0, palette=colors_left)
+    style_axes(ax0, "Building Counts on Different FIMs", "Inundation Surface", True)
+    for c in ax0.containers:
+        ax0.bar_label(
+            c, fmt="%.0f", label_type="edge", padding=3, fontsize=12, color="black"
+        )
+
+    # Right panel
+    colors_right = ["#ff5733", "#ffc300", "#28a745"]
+    sns.barplot(data=df_right, x="Category", y="Count", ax=ax1, palette=colors_right)
+    style_axes(ax1, "Contingency Flooded Building Counts", "Category", False)
+    for c in ax1.containers:
+        ax1.bar_label(
+            c, fmt="%.0f", label_type="edge", padding=3, fontsize=12, color="black"
+        )
+
+    # Combined legend
+    all_labels = ["Candidate", "Benchmark"] + third_raster_labels
+    all_colors = colors_left + colors_right
+    legend_handles = [
+        plt.Line2D(
+            [0],
+            [0],
+            marker="s",
+            color="w",
+            markerfacecolor=all_colors[i],
+            markersize=12,
+            label=f"{all_labels[i]} ({percentages[all_labels[i]]:.2f}%)",
+        )
+        for i in range(len(all_labels))
+    ]
+    ax_leg.legend(handles=legend_handles, fontsize=12, loc="center left", frameon=True)
     plot_dir = os.path.join(save_dir, "FinalPlots")
-    if not os.path.exists(plot_dir):
-        os.makedirs(plot_dir)
+    os.makedirs(plot_dir, exist_ok=True)
     output_path = os.path.join(plot_dir, f"BuildingCounts_{basename}.png")
-    fig.write_image(output_path, scale=500 / 96, engine="kaleido")
-    print(f"Performance metrics chart is saved as PNG at {output_path}")
-    fig.show()
+    fig_sb.savefig(output_path, dpi=400)
+    plt.close(fig_sb)
+
+    print(f"PNG were saved in : {output_path}")
 
 
 def process_TIFF(

{fimeval-0.1.54 → fimeval-0.1.56}/src/fimeval/ContingencyMap/evaluationFIM.py

@@ -116,7 +116,7 @@ def evaluateFIM(
         bounding_geom = AOI(benchmark_path, shapefile, save_dir)
 
     else:
-        print(f"--- {method.__name__} is processing ---")
+        print(f"**{method.__name__} is processing**")
         bounding_geom = method(smallest_raster_path, save_dir=save_dir)
 
     # Read and process benchmark raster
@@ -277,7 +277,8 @@ def evaluateFIM(
                                 out_transform1,
                             )
                             merged = out_image1 + out_image2_resized
-
+                            merged[merged==7] = 5 
+                            
             # Get Evaluation Metrics
             (
                 unique_values,
@@ -406,7 +407,6 @@ def EvaluateFIM(
         gdf = gpd.read_file(PWB_dir)
 
     # Grant the permission to the main directory
-    print(f"Fixing permissions for {main_dir}...")
     fix_permissions(main_dir)
 
     # runt the process
@@ -425,12 +425,11 @@ def EvaluateFIM(
             for tif_file in tif_files:
                 if "benchmark" in tif_file.name.lower() or "BM" in tif_file.name:
                     benchmark_path = tif_file
-                    print(f"---Benchmark: {tif_file.name}---")
                 else:
                     candidate_path.append(tif_file)
 
         if benchmark_path and candidate_path:
-            print(f"---Flood Inundation Evaluation of {folder_dir.name}---")
+            print(f"**Flood Inundation Evaluation of {folder_dir.name}**")
             Metrics = evaluateFIM(
                 benchmark_path,
                 candidate_path,

{fimeval-0.1.54 → fimeval-0.1.56}/src/fimeval/ContingencyMap/plotevaluationmetrics.py

@@ -2,19 +2,23 @@ import os
 import glob
 import pandas as pd
 import plotly.express as px
+import seaborn as sns
+import matplotlib.pyplot as plt
 
 
 # Function to plot individual metric scores
 def PlotMetrics(csv_path, method_path):
     metrics_df = pd.read_csv(csv_path)
-    # Extract relevant metrics
+
+    # Keep only the desired metrics
     metrics = metrics_df.loc[
         metrics_df["Metrics"].isin(
             ["CSI_values", "POD_values", "Acc_values", "Prec_values", "F1_values"]
         )
     ].copy()
 
-    metrics.loc[:, "Metrics"] = metrics["Metrics"].replace(
+    # Rename for presentation
+    metrics["Metrics"] = metrics["Metrics"].replace(
         {
             "CSI_values": "CSI",
             "POD_values": "POD",
@@ -23,54 +27,87 @@ def PlotMetrics(csv_path, method_path):
             "F1_values": "F1 Score",
         }
     )
+
     value_columns = metrics.select_dtypes(include="number").columns
 
+    # Output directory
+    plot_dir = os.path.join(method_path, "FinalPlots")
+    os.makedirs(plot_dir, exist_ok=True)
+
     for value_column in value_columns:
         metrics[value_column] = metrics[value_column].round(2)
 
-        # Create the bar plot
-        fig = px.bar(
+        # Showing with Plotly
+        fig_plotly = px.bar(
             metrics,
             x=value_column,
             y="Metrics",
-            title=f"Performance Metrics",
-            labels={value_column: "Score"},
             text=value_column,
             color="Metrics",
+            orientation="h",
             color_discrete_sequence=px.colors.qualitative.Set2,
+            title=f"Performance Metrics",
         )
-        fig.update_traces(texttemplate="%{text:.2f}", textposition="outside")
-        fig.update_layout(
-            yaxis_title="Metrics",
-            xaxis_title="Score",
+        fig_plotly.update_traces(texttemplate="%{text:.2f}", textposition="outside")
+        fig_plotly.update_layout(
             showlegend=False,
-            plot_bgcolor="rgba(0, 0, 0, 0)",
-            paper_bgcolor="rgba(0, 0, 0, 0)",
-            margin=dict(l=10, r=10, t=40, b=10),
-            xaxis=dict(showline=True, linewidth=2, linecolor="black"),
-            yaxis=dict(showline=True, linewidth=2, linecolor="black"),
+            plot_bgcolor="rgba(0,0,0,0)",
+            paper_bgcolor="rgba(0,0,0,0)",
             height=350,
             width=900,
+            xaxis=dict(showline=True, linewidth=2, linecolor="black", title="Score"),
+            yaxis=dict(showline=True, linewidth=2, linecolor="black"),
             title_font=dict(family="Arial", size=24, color="black"),
-            xaxis_title_font=dict(family="Arial", size=20, color="black"),
-            yaxis_title_font=dict(family="Arial", size=20, color="black"),
             font=dict(family="Arial", size=18, color="black"),
         )
+        fig_plotly.show()
 
-        # Save each plot as a PNG, using the column name as the filename
-        plot_dir = os.path.join(method_path, "FinalPlots")
-        if not os.path.exists(plot_dir):
-            os.makedirs(plot_dir)
+        # Save with Seaborn
+        sns.set_theme(style="whitegrid")
+        fig, ax = plt.subplots(figsize=(8, 3))
+        sns.barplot(
+            data=metrics,
+            x=value_column,
+            y="Metrics",
+            hue="Metrics",
+            palette="Set2",
+            ax=ax,
+            dodge=False,
+            legend=False,
+        )
 
-        output_filename = f"EvaluationMetrics_{value_column}.png"
-        output_path = os.path.join(plot_dir, output_filename)
+        # Annotate bars
+        for container in ax.containers:
+            ax.bar_label(container, fmt="%.2f", label_type="edge", fontsize=14)
 
-        # Save the plot as PNG
-        fig.write_image(output_path, engine="kaleido", scale=500 / 96)
-        print(
-            f"Performance metrics chart ({value_column}) saved as PNG at {output_path}"
+        # Styling
+        ax.set_title("Performance Metrics", fontsize=16)
+        ax.set_xlabel("Score", fontsize=16, color="black")  # just bigger, not bold
+        ax.set_ylabel("Metrics", fontsize=16, color="black")
+
+        ax.set_xticks([i / 10 for i in range(0, 11, 2)])
+        ax.set_xticklabels(
+            [f"{i/10:.1f}" for i in range(0, 11, 2)], fontsize=14, color="black"
         )
-        fig.show()
+
+        # Increase y-tick label font size
+        ax.tick_params(axis="y", labelsize=12, colors="black")
+        ax.tick_params(axis="x", labelsize=14, colors="black")
+
+        # Force spines black + thicker
+        ax.spines["left"].set_linewidth(1.5)
+        ax.spines["bottom"].set_linewidth(1.5)
+        ax.spines["left"].set_color("black")
+        ax.spines["bottom"].set_color("black")
+
+        sns.despine(right=True, top=True)
+
+        # Save to file
+        save_path = os.path.join(plot_dir, f"EvaluationMetrics_{value_column}.png")
+        plt.tight_layout()
+        fig.savefig(save_path, dpi=400)
+        plt.close(fig)
+        print(f"PNG saved at: {save_path}")
 
 
 def PlotEvaluationMetrics(main_dir, method_name, out_dir):

{fimeval-0.1.54 → fimeval-0.1.56}/src/fimeval/utilis.py

@@ -71,7 +71,6 @@ def reprojectFIMs(src_path, dst_path, target_crs):
                         resampling=Resampling.nearest,
                     )
         else:
-            print(f"Source raster is already in {target_crs}. No reprojection needed.")
             shutil.copy(src_path, dst_path)
         compress_tif_lzw(dst_path)
 
@@ -79,9 +78,7 @@ def reprojectFIMs(src_path, dst_path, target_crs):
 # Resample into the coarser resoution amoung all FIMS within the case
 def resample_to_resolution(src_path, x_resolution, y_resolution):
     src_path = Path(src_path)
-    print(src_path)
     temp_path = src_path.with_name(src_path.stem + "_resampled.tif")
-    print(temp_path)
 
     with rasterio.open(src_path) as src:
         transform = rasterio.transform.from_origin(
@@ -144,14 +141,13 @@ def MakeFIMsUniform(fim_dir, target_crs=None, target_resolution=None):
         if not final_crs:
             if all(is_within_conus(b, c) for b, c in zip(bounds_list, crs_list)):
                 final_crs = "EPSG:5070"
-                print(f"Defaulting to CONUS CRS: {final_crs}")
+                print(f"Defaulting to CONUS CRS: {final_crs}, Reprojecting.")
             else:
                 print(
                     "Mixed or non-CONUS CRS detected. Please provide a valid target CRS."
                 )
                 return
 
-        print(f"Reprojecting all rasters to {final_crs}")
         for src_path in tif_files:
             dst_path = processing_folder / src_path.name
             reprojectFIMs(str(src_path), str(dst_path), final_crs)

{fimeval-0.1.54 → fimeval-0.1.56}/src/fimeval.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: fimeval
-Version: 0.1.54
+Version: 0.1.56
 Summary: A Framework for Automatic Evaluation of Flood Inundation Mapping Predictions Evaluation
 Author: Surface Dynamics Modeling Lab
 Author-email: Supath Dhital <sdhital@crimson.ua.edu>, Dipshika Devi <ddevi@ua.edu>
@@ -683,6 +683,7 @@ Requires-Dist: notebook<8.0.0,>=7.3.2
 Requires-Dist: boto3<2.0.0,>=1.36.16
 Requires-Dist: geemap
 Requires-Dist: uv
+Requires-Dist: seaborn
 Provides-Extra: dev
 Requires-Dist: pytest; extra == "dev"
 Requires-Dist: black; extra == "dev"

{fimeval-0.1.54 → fimeval-0.1.56}/src/fimeval.egg-info/requires.txt

@@ -11,6 +11,7 @@ notebook<8.0.0,>=7.3.2
 boto3<2.0.0,>=1.36.16
 geemap
 uv
+seaborn
 
 [dev]
 pytest

{fimeval-0.1.54 → fimeval-0.1.56}/tests/test_evaluationfim.py

@@ -29,7 +29,7 @@ countryISO = "USA"
 def test_evaluation_framework():
     # Run the evaluation
     # It has the Permanent Water Bodies (PWB) dataset as default for United States
-    fe.EvaluateFIM(Main_dir, method_name, output_dir)
+    # fe.EvaluateFIM(Main_dir, method_name, output_dir)
 
     # OR, If the Evaluation Study Area is outside the US or, user has their own PWB dataset
     # fe.EvaluateFIM(Main_dir, method_name, output_dir)