fimeval 0.1.54__py3-none-any.whl → 0.1.55__py3-none-any.whl

fimeval/BuildingFootprint/evaluationwithBF.py

@@ -6,6 +6,8 @@ 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
 
 
 def Changeintogpkg(input_path, output_dir, layer_name):
@@ -19,7 +21,8 @@ def Changeintogpkg(input_path, output_dir, layer_name):
         output_gpkg = os.path.join(output_dir, f"{layer_name}.gpkg")
         gdf.to_file(output_gpkg, driver="GPKG")
         return output_gpkg
-
+    
+    
 def GetFloodedBuildingCountInfo(
     building_fp_path,
     study_area_path,
@@ -57,7 +60,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,11 +76,9 @@ 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"]
 
-
     total_buildings = len(clipped_buildings)
     percentages = {
         key: (count / total_buildings) * 100 if total_buildings > 0 else 0
@@ -92,15 +92,11 @@ 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"]
-        )
+        BDR = (centroid_counts["Candidate"] - centroid_counts["Benchmark"]) / centroid_counts["Benchmark"]
     else:
-            BDR = 0
-        
+        BDR = 0
+
     counts_data = {
         "Category": [
             "Candidate",
@@ -125,13 +121,12 @@ def GetFloodedBuildingCountInfo(
             f"{BDR:.3f}",
         ],
     }
-
     counts_df = pd.DataFrame(counts_data)
-    csv_file_path = os.path.join(
-        save_dir, "EvaluationMetrics", f"BuildingCounts_{basename}.csv"
-    )
+    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"],
@@ -140,80 +135,108 @@ def GetFloodedBuildingCountInfo(
     ]
 
     fig = make_subplots(
-        rows=1,
-        cols=2,
-        subplot_titles=(
-            "Building Counts on Different FIMs",
-            "Contingency Flooded Building Counts",
-        ),
+        rows=1, cols=2,
+        subplot_titles=("Building Counts on Different FIMs", "Contingency Flooded Building Counts"),
     )
 
     fig.add_trace(
         go.Bar(
-            x=["Candidate"],
-            y=[centroid_counts["Candidate"]],
-            text=[f"{centroid_counts['Candidate']}"],
-            textposition="auto",
-            marker_color="#1c83eb",
-            marker_line_color="black",
-            marker_line_width=1,
+            x=["Candidate"], y=[centroid_counts["Candidate"]],
+            text=[f"{centroid_counts['Candidate']}"], textposition="auto",
+            marker_color="#1c83eb", marker_line_color="black", marker_line_width=1,
             name=f"Candidate ({percentages['Candidate']:.2f}%)",
         ),
-        row=1,
-        col=1,
+        row=1, col=1,
     )
-
     fig.add_trace(
         go.Bar(
-            x=["Benchmark"],
-            y=[centroid_counts["Benchmark"]],
-            text=[f"{centroid_counts['Benchmark']}"],
-            textposition="auto",
-            marker_color="#a4490e",
-            marker_line_color="black",
-            marker_line_width=1,
+            x=["Benchmark"], y=[centroid_counts["Benchmark"]],
+            text=[f"{centroid_counts['Benchmark']}"], textposition="auto",
+            marker_color="#a4490e", marker_line_color="black", marker_line_width=1,
             name=f"Benchmark ({percentages['Benchmark']:.2f}%)",
         ),
-        row=1,
-        col=1,
+        row=1, 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]],
-                y=[third_raster_counts[i]],
-                text=[f"{third_raster_counts[i]}"],
-                textposition="auto",
+                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}%)",
+                marker_line_color="black", marker_line_width=1,
+                name=f"{label} ({percentages[label]:.2f}%)",
             ),
-            row=1,
-            col=2,
+            row=1, col=2,
         )
 
     fig.update_layout(
         title="Flooded Building Counts",
         xaxis_title="Inundation Surface",
         yaxis_title="Flooded Building Counts",
-        width=1100,
-        height=400,
-        plot_bgcolor="rgba(0, 0, 0, 0)",
-        paper_bgcolor="rgba(0, 0, 0, 0)",
-        showlegend=True,
-        font=dict(family="Arial", size=18, color="black"),
+        width=1100, height=400,
+        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 saving only
+    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, axes = plt.subplots(1, 2, figsize=(8, 3), constrained_layout=True)
+
+    def style_axes(ax, title_text, xlab, show_ylabel: bool):
+        # Adding a bit of padding so bar labels don’t overlap with the title
+        ax.set_title(title_text, fontsize=16, pad=20)
+        ax.set_xlabel(xlab, fontsize=14, color="black")
+        if show_ylabel:
+            ax.set_ylabel("Flooded Building Counts", fontsize=14, color="black")
+        else:
+            ax.set_ylabel("")
+
+        # Thicker black left/bottom spines
+        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=12, colors="black")
+        ax.tick_params(axis="y", labelsize=12, colors="black")
+
+    # Left panel
+    ax0 = axes[0]
+    sns.barplot(data=df_left, x="Category", y="Count", ax=ax0,
+                palette=["#1c83eb", "#a4490e"])
+    style_axes(ax0, "Building Counts on Different FIMs", "Inundation Surface", show_ylabel=True)
+    for c in ax0.containers:
+        ax0.bar_label(c, fmt="%.0f", label_type="edge", padding=3, fontsize=14, color="black")
+
+    # Right panel
+    ax1 = axes[1]
+    sns.barplot(data=df_right, x="Category", y="Count", ax=ax1,
+                palette=["#ff5733", "#ffc300", "#28a745"])
+    style_axes(ax1, "Contingency Flooded Building Counts", "Category", show_ylabel=False)
+    for c in ax1.containers:
+        ax1.bar_label(c, fmt="%.0f", label_type="edge", padding=3, fontsize=14, color="black")
 
     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(
     tif_files, contingency_files, building_footprint, boundary, method_path

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
@@ -406,7 +406,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 +424,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/ContingencyMap/plotevaluationmetrics.py

@@ -2,75 +2,108 @@ 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"]
-        )
+        metrics_df["Metrics"].isin([
+            "CSI_values", "POD_values", "Acc_values", "Prec_values", "F1_values"
+        ])
     ].copy()
 
-    metrics.loc[:, "Metrics"] = metrics["Metrics"].replace(
-        {
-            "CSI_values": "CSI",
-            "POD_values": "POD",
-            "Acc_values": "Accuracy",
-            "Prec_values": "Precision",
-            "F1_values": "F1 Score",
-        }
-    )
+    # Rename for presentation
+    metrics["Metrics"] = metrics["Metrics"].replace({
+        "CSI_values": "CSI",
+        "POD_values": "POD",
+        "Acc_values": "Accuracy",
+        "Prec_values": "Precision",
+        "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}"
-        )
-        fig.show()
+        # 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")
+
+        # 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/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.dist-info → fimeval-0.1.55.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: fimeval
-Version: 0.1.54
+Version: 0.1.55
 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.dist-info → fimeval-0.1.55.dist-info}/RECORD

@@ -1,17 +1,17 @@
 fimeval/__init__.py,sha256=HZJKq7XEhL6TnwFkhpf8NcEQ5h7zPQ3XJh3z5gF0gQ8,603
-fimeval/utilis.py,sha256=JYpWTJx3WOuC0wRr3nb-cwMPf0SAGIli5dJdKFS23w8,6954
+fimeval/utilis.py,sha256=KdU6pSMS8dPf1zFyS-SOUdDSJr87IQ28uSG3ZuVfq1A,6782
 fimeval/BuildingFootprint/__init__.py,sha256=oP9YWLdo6ANzSQFxYLv7Ku_26AY5NkLNhZLK28ICMLo,109
-fimeval/BuildingFootprint/evaluationwithBF.py,sha256=2QBYbuJg98B5_PKDLnmV7M1dDVuOkJd6_rP0EdmoF-8,15645
+fimeval/BuildingFootprint/evaluationwithBF.py,sha256=84s6y-S3GAuclgAEXQ4DmBAtnnZ3t0LFhIPjeSMtzPw,17224
 fimeval/BuildingFootprint/microsoftBF.py,sha256=73M_e_n_3lsejsnP9eX06hou0sr2-yaWdjuoUwZ8O2Y,4063
 fimeval/ContingencyMap/PWBs3.py,sha256=UFICxO58c2fA9mIffH4ooqphv3ZXe6yX8QzpRjtI6fs,1275
 fimeval/ContingencyMap/__init__.py,sha256=ckps2dyg6aci3TA-3P7oTMcCAcSTz9AA6sndHtZEwdE,259
-fimeval/ContingencyMap/evaluationFIM.py,sha256=JiQfXiiQXZd5T4PBqEyU567PHGaEXr8RUsnGmWQxzjc,17390
+fimeval/ContingencyMap/evaluationFIM.py,sha256=MTh9W5M-g43uzgW8ci0GJ65KKCEsepbqiE8o8Ox0Pj4,17270
 fimeval/ContingencyMap/methods.py,sha256=kbutfo9FUH-yjvnOXxwLpdErUuebMJ8NjCroNWIYCjo,3299
 fimeval/ContingencyMap/metrics.py,sha256=jwOia0Nl7aU7AuGJFAcQ4fVENnp2G_5W6JSJBzo1-_4,1094
-fimeval/ContingencyMap/plotevaluationmetrics.py,sha256=3bKfPKZnMR39dA3teDVpQBeTFKnF9v_2Vku0JNVGggs,3921
+fimeval/ContingencyMap/plotevaluationmetrics.py,sha256=CLw3y3XB3XhGm_X2_U3JSlPq0DYKxATNSW7lKPvMkMA,4838
 fimeval/ContingencyMap/printcontingency.py,sha256=-1H_Ze2TbRSER7vy7bd0HvxnziNzPPOIPOm2YhB7r4A,5422
-fimeval-0.1.54.dist-info/licenses/LICENSE.txt,sha256=hIahDEOTzuHCU5J2nd07LWwkLW7Hko4UFO__ffsvB-8,34523
-fimeval-0.1.54.dist-info/METADATA,sha256=MKaEZodfSX6Ot7pbV0uTqTy-E-dmc7zXC7C8X3YctR8,56090
-fimeval-0.1.54.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-fimeval-0.1.54.dist-info/top_level.txt,sha256=F4QW50msI8sRrX_DK3NQ-s3swQ4-2_5Ty3mfm9ZMc6k,8
-fimeval-0.1.54.dist-info/RECORD,,
+fimeval-0.1.55.dist-info/licenses/LICENSE.txt,sha256=hIahDEOTzuHCU5J2nd07LWwkLW7Hko4UFO__ffsvB-8,34523
+fimeval-0.1.55.dist-info/METADATA,sha256=A_OiFx3NhXA_t9P0mW3qYGFchKw6UlGKuOwXFv7YzYk,56113
+fimeval-0.1.55.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+fimeval-0.1.55.dist-info/top_level.txt,sha256=F4QW50msI8sRrX_DK3NQ-s3swQ4-2_5Ty3mfm9ZMc6k,8
+fimeval-0.1.55.dist-info/RECORD,,