Merge pull request #34 from quishqa/quishqa/prep_pypi

Update setup.py install_requires

Author: quishqa

requirements.txt

@@ -1,6 +1,7 @@
 numpy==1.26.4
 bottleneck==1.3.8
 xarray==2024.2.0
+geopandas==0.14.3
 netcdf4==1.6.5
 dask==2024.2.1
 hdbscan==0.8.33

setup.py

@@ -10,8 +10,7 @@
     author="Mario Gavidia-Calderón",
     author_email="mario.calderon@iag.usp.br",
     packages=find_packages(),
-    setup_requires=["numpy"],
-    install_requires=["numpy"],
+    install_requires=["xarray", "pyproj", "osmnx", "pseudonetcdf", "geopandas"],
     long_description=long_description,
     long_description_content_type="text/markdown",
     license="MIT",

siem/point.py

@@ -16,7 +16,6 @@
 import pandas as pd
 import geopandas as gpd
 import xarray as xr
-from netCDF4 import Dataset
 from siem.proxy import create_wrf_grid, configure_grid_spatial
 import pyproj
 
@@ -101,22 +100,35 @@ def retrive_proj_from(geogrid_path: str):
     a = 6370000.0
     b = 6370000.0
 
-    lcc = pyproj.Proj(proj='lcc', lat_1=geo_ds.TRUELAT1, lat_2=geo_ds.TRUELAT2, 
-                      lat_0=geo_ds.MOAD_CEN_LAT, lon_0=geo_ds.STAND_LON, 
-                      a=a, b=b)
-    merc = pyproj.Proj(proj='merc', lon_0=geo_ds.STAND_LON, lat_ts=geo_ds.TRUELAT1,
-                       a=a, b=b)
-    stere = pyproj.Proj(proj='stere', lat_0=geo_ds.TRUELAT1, lon_0=geo_ds.STAND_LON,
-                        lat_ts=geo_ds.TRUELAT1, a=a, b=b)
-    latlon = pyproj.Proj(proj='longlat', lon_0=geo_ds.STAND_LON, a=a, b=b)
+    lcc = pyproj.Proj(
+        proj="lcc",
+        lat_1=geo_ds.TRUELAT1,
+        lat_2=geo_ds.TRUELAT2,
+        lat_0=geo_ds.MOAD_CEN_LAT,
+        lon_0=geo_ds.STAND_LON,
+        a=a,
+        b=b,
+    )
+    merc = pyproj.Proj(
+        proj="merc", lon_0=geo_ds.STAND_LON, lat_ts=geo_ds.TRUELAT1, a=a, b=b
+    )
+    stere = pyproj.Proj(
+        proj="stere",
+        lat_0=geo_ds.TRUELAT1,
+        lon_0=geo_ds.STAND_LON,
+        lat_ts=geo_ds.TRUELAT1,
+        a=a,
+        b=b,
+    )
+    latlon = pyproj.Proj(proj="longlat", lon_0=geo_ds.STAND_LON, a=a, b=b)
 
     proj_codes = {
-        1: lcc, # lambert
-        2: stere, # polar
-        3: merc, # merc 
-        6: latlon # latlon
+        1: lcc,  # lambert
+        2: stere,  # polar
+        3: merc,  # merc
+        6: latlon,  # latlon
     }
-    
+
     wrf_proj = proj_codes[geo_ds.MAP_PROJ]
     wrf_crs = pyproj.CRS.from_proj4(str(wrf_proj))
     return wrf_crs