hydroanomaly 0.7.2__tar.gz → 1.0.0__tar.gz

{hydroanomaly-0.7.2 → hydroanomaly-1.0.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: hydroanomaly
-Version: 0.7.2
+Version: 1.0.0
 Summary: A Python package for hydro anomaly detection with simple USGS data retrieval
 Author-email: Ehsan Kahrizi <ehsan.kahrizi@usu.edu>
 License: MIT License

{hydroanomaly-0.7.2 → hydroanomaly-1.0.0}/hydroanomaly/__init__.py

@@ -9,12 +9,12 @@ A simple Python package with just 3 modules:
 That's it - nothing else!
 """
 
-__version__ = "0.7.2"
+__version__ = "1.0.0"
 __author__ = "Ehsan Kahrizi (Ehsan.kahrizi@usu.edu)"
 
 # Import the 3 simple modules
 from .usgs_turbidity import get_turbidity, get_usgs_turbidity
-from .sentinel_bands import get_sentinel_bands, get_satellite_data, get_sentinel, get_sentinel_bands_gee
+from .sentinel_bands import get_sentinel_bands, get_satellite_data, get_sentinel, get_sentinel_bands_gee, show_sentinel_ndwi_map
 from .visualize import plot_timeseries, plot_turbidity, plot_sentinel, plot_comparison, plot, visualize
 
 # Export everything
@@ -28,6 +28,7 @@ __all__ = [
     'get_sentinel_bands',
     'get_satellite_data', 
     'get_sentinel',
+    'show_sentinel_ndwi_map',
 
 
     # Visualization functions

{hydroanomaly-0.7.2 → hydroanomaly-1.0.0}/hydroanomaly/sentinel_bands.py

@@ -4,13 +4,15 @@ Sentinel-2 Satellite Data Retrieval using Google Earth Engine (GEE)
 This module provides a function to retrieve Sentinel-2 satellite band data
 for a specified location and time period, with masking and cloud filtering.
 """
-
+import ee
+import geemap 
 import pandas as pd
 import numpy as np
 from datetime import datetime, timedelta
 import requests
 import warnings
 
+# Retrive data from Google Earth Engine ========================================================
 def get_sentinel_bands_gee(
     latitude: float,
     longitude: float,
@@ -98,6 +100,54 @@ def get_sentinel_bands_gee(
         df = df.set_index('date')
     return df
 
+
+# Showing map with NDWI (Normalized Difference Water Index) ========================================================
+def show_sentinel_ndwi_map(
+    latitude: float,
+    longitude: float,
+    start_date: str,
+    end_date: str,
+    buffer_meters: int = 20,
+    cloudy_pixel_percentage: int = 20,
+    zoom: int = 15
+):
+    """
+    Display an interactive map showing the NDWI, point, and buffer.
+
+    Args:
+        latitude (float): Latitude of the center point.
+        longitude (float): Longitude of the center point.
+        start_date (str): Start date as "YYYY-MM-DD".
+        end_date (str): End date as "YYYY-MM-DD".
+        buffer_meters (int): Buffer radius in meters.
+        cloudy_pixel_percentage (int): Max allowed cloud percentage.
+        zoom (int): Zoom level for map.
+    """
+    point = ee.Geometry.Point([longitude, latitude])
+    buffer = point.buffer(buffer_meters)
+
+    image = (ee.ImageCollection("COPERNICUS/S2_SR_HARMONIZED")
+             .filterBounds(buffer)
+             .filterDate(start_date, end_date)
+             .filter(ee.Filter.lt("CLOUDY_PIXEL_PERCENTAGE", cloudy_pixel_percentage))
+             .median())
+
+    ndwi = image.normalizedDifference(["B3", "B8"]).rename("NDWI")
+
+    ndwi_vis = {
+        'min': 0,
+        'max': 1,
+        'palette': ['white', 'cyan', 'blue']}
+
+    Map = geemap.Map()
+    Map.centerObject(buffer, zoom=zoom)
+    Map.addLayer(ndwi, ndwi_vis, "NDWI (Water)")
+    Map.addLayer(point, {'color': 'yellow'}, 'Point')
+    Map.addLayer(buffer, {'color': 'red'}, 'Buffer')
+    Map.add_colorbar(ndwi_vis, label="NDWI", layer_name="NDWI (Water)")
+    return Map
+
+
 # Aliases for user convenience
 get_sentinel_bands = get_sentinel_bands_gee
 get_satellite_data = get_sentinel_bands_gee

{hydroanomaly-0.7.2 → hydroanomaly-1.0.0}/hydroanomaly/visualize.py

@@ -11,7 +11,7 @@ import pandas as pd
 import numpy as np
 from datetime import datetime
 
-
+# ==============================================================================================
 def plot_timeseries(data: pd.DataFrame, title: str = "Time Series Data", save_file: str = None) -> None:
     """
     Create a simple time series plot.
@@ -26,10 +26,10 @@ def plot_timeseries(data: pd.DataFrame, title: str = "Time Series Data", save_fi
     """
     
     if data.empty:
-        print("❌ No data to plot")
+        print("No data to plot")
         return
     
-    print(f"📊 Creating plot: {title}")
+    print(f"Creating plot: {title}")
     
     # Create figure
     plt.figure(figsize=(12, 6))
@@ -59,12 +59,12 @@ def plot_timeseries(data: pd.DataFrame, title: str = "Time Series Data", save_fi
     # Save if requested
     if save_file:
         plt.savefig(save_file, dpi=300, bbox_inches='tight')
-        print(f"💾 Plot saved as {save_file}")
+        print(f"Plot saved as {save_file}")
     
     plt.show()
-    print("✅ Plot created successfully!")
-
+    print("Plot created successfully!")
 
+# ==============================================================================================
 def plot_turbidity(turbidity_data: pd.DataFrame, save_file: str = None) -> None:
     """
     Create a turbidity-specific plot with appropriate formatting.
@@ -75,10 +75,10 @@ def plot_turbidity(turbidity_data: pd.DataFrame, save_file: str = None) -> None:
     """
     
     if turbidity_data.empty:
-        print("❌ No turbidity data to plot")
+        print("No turbidity data to plot")
         return
     
-    print("🌫️  Creating turbidity plot")
+    print("Creating turbidity plot")
     
     plt.figure(figsize=(12, 6))
     
@@ -92,7 +92,7 @@ def plot_turbidity(turbidity_data: pd.DataFrame, save_file: str = None) -> None:
     plt.axhline(y=25, color='red', linestyle='--', alpha=0.7, label='High (25 NTU)')
     
     # Format plot
-    plt.title('💧 Turbidity Time Series', fontsize=14, fontweight='bold', pad=20)
+    plt.title('Turbidity Time Series', fontsize=14, fontweight='bold', pad=20)
     plt.xlabel('Date', fontsize=12)
     plt.ylabel('Turbidity (NTU)', fontsize=12)
     plt.grid(True, alpha=0.3)
@@ -107,12 +107,12 @@ def plot_turbidity(turbidity_data: pd.DataFrame, save_file: str = None) -> None:
     # Save if requested
     if save_file:
         plt.savefig(save_file, dpi=300, bbox_inches='tight')
-        print(f"💾 Turbidity plot saved as {save_file}")
+        print(f"Turbidity plot saved as {save_file}")
     
     plt.show()
-    print("✅ Turbidity plot created!")
-
+    print("Turbidity plot created!")
 
+# ==============================================================================================
 def plot_sentinel(sentinel_data: pd.DataFrame, save_file: str = None) -> None:
     """
     Create a Sentinel satellite data plot.
@@ -123,10 +123,10 @@ def plot_sentinel(sentinel_data: pd.DataFrame, save_file: str = None) -> None:
     """
     
     if sentinel_data.empty:
-        print("❌ No Sentinel data to plot")
+        print("No Sentinel data to plot")
         return
     
-    print("🛰️  Creating Sentinel bands plot")
+    print("Creating Sentinel bands plot")
     
     plt.figure(figsize=(12, 8))
     
@@ -146,7 +146,7 @@ def plot_sentinel(sentinel_data: pd.DataFrame, save_file: str = None) -> None:
                 label=column, color=color, linewidth=2, marker='o', markersize=4)
     
     # Format plot
-    plt.title('🛰️ Sentinel Satellite Data', fontsize=14, fontweight='bold', pad=20)
+    plt.title('Sentinel Satellite Data', fontsize=14, fontweight='bold', pad=20)
     plt.xlabel('Date', fontsize=12)
     plt.ylabel('Digital Number / Index Value', fontsize=12)
     plt.grid(True, alpha=0.3)
@@ -161,12 +161,12 @@ def plot_sentinel(sentinel_data: pd.DataFrame, save_file: str = None) -> None:
     # Save if requested
     if save_file:
         plt.savefig(save_file, dpi=300, bbox_inches='tight')
-        print(f"💾 Sentinel plot saved as {save_file}")
+        print(f"Sentinel plot saved as {save_file}")
     
     plt.show()
-    print("✅ Sentinel plot created!")
-
+    print("Sentinel plot created!")
 
+# ==============================================================================================
 def plot_comparison(data1: pd.DataFrame, data2: pd.DataFrame, 
                    label1: str = "Dataset 1", label2: str = "Dataset 2",
                    title: str = "Data Comparison", save_file: str = None) -> None:
@@ -183,10 +183,10 @@ def plot_comparison(data1: pd.DataFrame, data2: pd.DataFrame,
     """
     
     if data1.empty and data2.empty:
-        print("❌ No data to plot")
+        print("No data to plot")
         return
     
-    print(f"📊 Creating comparison plot: {title}")
+    print(f"Creating comparison plot: {title}")
     
     fig, (ax1, ax2) = plt.subplots(2, 1, figsize=(12, 10), sharex=True)
     
@@ -215,10 +215,10 @@ def plot_comparison(data1: pd.DataFrame, data2: pd.DataFrame,
     # Save if requested
     if save_file:
         plt.savefig(save_file, dpi=300, bbox_inches='tight')
-        print(f"💾 Comparison plot saved as {save_file}")
+        print(f"Comparison plot saved as {save_file}")
     
     plt.show()
-    print("✅ Comparison plot created!")
+    print("Comparison plot created!")
 
 
 # Simple aliases

{hydroanomaly-0.7.2 → hydroanomaly-1.0.0}/hydroanomaly.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: hydroanomaly
-Version: 0.7.2
+Version: 1.0.0
 Summary: A Python package for hydro anomaly detection with simple USGS data retrieval
 Author-email: Ehsan Kahrizi <ehsan.kahrizi@usu.edu>
 License: MIT License

{hydroanomaly-0.7.2 → hydroanomaly-1.0.0}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "hydroanomaly"
-version = "0.7.2"
+version = "1.0.0"
 authors = [
     {name = "Ehsan Kahrizi", email = "ehsan.kahrizi@usu.edu"},
 ]