grdwindinversion 0.2.6__py3-none-any.whl → 0.3.1__py3-none-any.whl

grdwindinversion/utils.py

@@ -1,7 +1,22 @@
+import os
+import time
 import logging
 import xsarsea
 
 
+logging.basicConfig(level=logging.INFO,
+                    format='%(asctime)s - %(name)s - %(levelname)s - %(message)s')
+logger = logging.getLogger('grdwindinversion')
+
+
+mem_monitor = True
+try:
+    from psutil import Process
+except ImportError:
+    logger.warning("psutil module not found. Disabling memory monitor")
+    mem_monitor = False
+
+
 def check_incidence_range(incidence, models, **kwargs):
     """
     Check if the incidence range of the dataset is within the range of the LUT of the model.
@@ -34,7 +49,7 @@ def check_incidence_range(incidence, models, **kwargs):
                 f"GMF {model_name} inc_range will be changed by kwargs to {kwargs['inc_range']}")
             lut_range = kwargs['inc_range']
 
-        inc_range = [incidence.values.min(), incidence.values.max()]
+        inc_range = [incidence.min(), incidence.max()]
         if (inc_range[0] >= lut_range[0] and inc_range[1] <= lut_range[1]):
             rets.append(True)
         else:
@@ -81,3 +96,28 @@ def get_pol_ratio_name(model_co):
             return "not_written_in_lut"
     else:
         return '/'
+
+
+def timing(logger=logger.debug):
+    """provide a @timing decorator() for functions, that log time spent in it"""
+
+    def decorator(f):
+        # @wraps(f)
+        def wrapper(*args, **kwargs):
+            mem_str = ''
+            process = None
+            if mem_monitor:
+                process = Process(os.getpid())
+                startrss = process.memory_info().rss
+            starttime = time.time()
+            result = f(*args, **kwargs)
+            endtime = time.time()
+            if mem_monitor:
+                endrss = process.memory_info().rss
+                mem_str = 'mem: %+.1fMb' % ((endrss - startrss) / (1024 ** 2))
+            logger(
+                'timing %s : %.2fs. %s' % (f.__name__, endtime - starttime, mem_str))
+            return result
+        wrapper.__doc__ = f.__doc__
+        return wrapper
+    return decorator

grdwindinversion-0.3.1.dist-info/METADATA

@@ -0,0 +1,67 @@
+Metadata-Version: 2.1
+Name: grdwindinversion
+Version: 0.3.1
+Summary: Package to perform Wind inversion from GRD Level-1 SAR images
+Author-email: Antoine Grouazel <antoine.grouazel@ifremer.fr>
+License: MIT
+Keywords: xarray,earth-observation,remote-sensing,satellite-imagery,Sentinel-1,RCM,RadarSat2,sar,synthetic-aperture-radar
+Classifier: Development Status :: 4 - Beta
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Operating System :: OS Independent
+Classifier: Programming Language :: Python
+Classifier: Programming Language :: Python :: 3 :: Only
+Classifier: Programming Language :: Python :: 3.11
+Classifier: Programming Language :: Python :: 3.12
+Classifier: Intended Audience :: Science/Research
+Classifier: Topic :: Scientific/Engineering
+Requires-Python: >=3.9
+Description-Content-Type: text/markdown
+License-File: LICENSE
+License-File: AUTHORS.rst
+Requires-Dist: xsar
+Requires-Dist: xsarsea
+Requires-Dist: xarray
+Requires-Dist: xarray-datatree
+Requires-Dist: pyyaml
+Requires-Dist: numpy
+Requires-Dist: scipy
+Requires-Dist: fsspec
+Requires-Dist: aiohttp
+
+
+[![Python Version](https://img.shields.io/pypi/pyversions/grdwindinversion.svg)](https://pypi.org/project/grdwindinversion/)
+[![Dependencies Status](https://img.shields.io/badge/dependencies-up%20to%20date-brightgreen.svg)](https://github.com/umr-lops/grdwindinversion/pulls?utf8=%E2%9C%93&q=is%3Apr%20author%3Aapp%2Fdependabot)
+
+Package to perform Wind inversion from GRD Level-1 SAR images
+
+-   Free software: MIT license
+-   Documentation: https://grdwindinversion.readthedocs.io.
+
+## Usage
+
+```python
+
+   SAR_L1-to-L2_wind_processor -h
+   usage: SAR_L1-to-L2_wind_processor [-h] --input_file INPUT_FILE [--config_file CONFIG_FILE] --outputdir OUTPUTDIR [--verbose] [--overwrite]
+
+   Perform inversion from S1(L1-GRD) SAFE, L1-RCM, L1-RS2 ; using xsar/xsarsea tools
+
+   options:
+     -h, --help            show this help message and exit
+     --input_file INPUT_FILE
+                           input file path
+     --config_file CONFIG_FILE
+                           config file path [if not provided will take config file based on input file]
+     --outputdir OUTPUTDIR
+     --verbose
+     --overwrite           overwrite existing .nc files [default is False]
+```
+
+## Features
+
+This Python library (based on `xarray`) allows to perform wind inversion from level-1 GRD (projected magnitude image).
+Mission supported:
+
+-   Sentinel-1
+-   RCM
+-   RadarSat-2

grdwindinversion-0.3.1.dist-info/RECORD

@@ -0,0 +1,22 @@
+grdwindinversion/.gitignore,sha256=vmDRHGeESYckUdUztsPr7u6ZNfBYMCVR3GE3AJCBnsY,1204
+grdwindinversion/__init__.py,sha256=VQc2V_j124NX2Gp06Go1oSZDSF4lBIVPYT-W9ZkNBXE,415
+grdwindinversion/config_prod.yaml,sha256=VEAgtJt894v3u4AhYOlcCWjb379vPDMfwBb4PKeDmAw,1314
+grdwindinversion/config_prod_recal.yaml,sha256=74n2J4tzSx1I5vaat3ydwig1jyfDgxGTGDxxw8ALC-U,1222
+grdwindinversion/config_prod_recal_streaks_nrcsmod.yaml,sha256=ob0YmzHrwkI4GLxaR9HzidYtXBrTOAlEYwPYdCK7AqA,1092
+grdwindinversion/config_prod_streaks.yaml,sha256=lJMl4qH6XYqO8ich7yW0WDQsQDrjmWThlLyx1xZy1DA,1313
+grdwindinversion/config_prod_streaks_nrcsmod.yaml,sha256=4yC4KcsC7rUumYHwhpKxM8_LDMe_Je-RIu4-9e1ayp8,1312
+grdwindinversion/data_config.yaml,sha256=FnglUHbAtGwHWg8w86hwZt3-vo-dY8uIjxqyI1jZpv8,567
+grdwindinversion/gradientFeatures.py,sha256=pBmsk6L8V9OrfSaAlFnxWztowSMEq2uoNdK93-drCuQ,17234
+grdwindinversion/inversion.py,sha256=Tup2TbjKzxHA4aF-0BOXFrUX_ILHbrKtMi2SVZ4JMgo,53972
+grdwindinversion/load_config.py,sha256=o0zVRpw3cIk-0fAsgOAb5JYzbHd3wdKmFskT8B3FAok,875
+grdwindinversion/main.py,sha256=0DLDW0i14CWhDKAGegzOKlWnCkag_cGn3vzMIAKcv4g,2658
+grdwindinversion/utils.py,sha256=bowCr2wKqXE7j0Uo5xkFMjWkr7NhnmjCaPFTO4MYNV4,3803
+grdwindinversion/utils_memory.py,sha256=1N3Kh4qVZPELPU6I4onbmkur1CZd7EHZqfbSAa6eaVc,1480
+grdwindinversion/.github/ISSUE_TEMPLATE.md,sha256=qiM_a7CCUz3fSrz3Q20Se1nwPNFS8QCc8tkwK_0DSCo,327
+grdwindinversion-0.3.1.dist-info/AUTHORS.rst,sha256=KmhW_5LBKGTIGwWEVkoTm1qx_bvdDR3yYL-1cwbDOFQ,218
+grdwindinversion-0.3.1.dist-info/LICENSE,sha256=-B8mBiTeY3J7OLuayiV1myqmc7yeijBc7s34kc8RTmg,1075
+grdwindinversion-0.3.1.dist-info/METADATA,sha256=-7HAO9FEm6AXR0aFdVtNptBOzVXbRQxzFHMuhBCysmw,2441
+grdwindinversion-0.3.1.dist-info/WHEEL,sha256=PZUExdf71Ui_so67QXpySuHtCi3-J3wvF4ORK6k_S8U,91
+grdwindinversion-0.3.1.dist-info/entry_points.txt,sha256=2rjvlVCy0iasRXjOz3kOIGuy2OCGQ-VTNuwuViQ6cMM,95
+grdwindinversion-0.3.1.dist-info/top_level.txt,sha256=z6lPix3QPEYOo37qq8plA2hY7S3C8MQZY81agRlksMI,17
+grdwindinversion-0.3.1.dist-info/RECORD,,

{grdwindinversion-0.2.6.dist-info → grdwindinversion-0.3.1.dist-info}/WHEEL

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

grdwindinversion/streaks.py

@@ -1,79 +0,0 @@
-import xarray as xr
-import xsarsea.gradients
-import xarray as xr
-from scipy.ndimage import binary_dilation
-import numpy as np
-
-
-def get_streaks(xr_dataset, xr_dataset_100):
-    """
-    Get the streaks from the wind field.
-
-    Parameters
-    ----------
-    xr_dataset : xarray.Dataset
-        dataset at user resolution.
-    xr_dataset_100 : xarray.Dataset
-        dataset at 100m resolution.
-
-    Returns
-    -------
-    xarray.Dataset
-        Extract wind direction from Koch Method using xsarsea tools.
-    """
-
-    # return empy dataArray, waiting for solution
-    return xr.DataArray(data=np.nan * np.ones([len(xr_dataset.coords[dim]) for dim in ['line','sample']]),
-                        dims=['line','sample'],
-                        coords=[xr_dataset.coords[dim] for dim in ['line','sample']])
-    #
-
-    """
-    gradients = xsarsea.gradients.Gradients(xr_dataset_100['sigma0_detrend'], windows_sizes=[
-                                            1600, 3200], downscales_factors=[1, 2], window_step=1)
-
-    # get gradients histograms as an xarray dataset
-    hist = gradients.histogram
-
-    # get orthogonals gradients
-    hist['angles'] = hist['angles'] + np.pi/2
-
-    # mean
-    hist_mean = hist.mean(['downscale_factor', 'window_size', 'pol'])
-
-    # smooth
-    hist_mean_smooth = hist_mean.copy()
-    hist_mean_smooth['weight'] = xsarsea.gradients.circ_smooth(
-        hist_mean['weight'])
-
-    # smooth only
-    # hist_smooth = hist.copy()
-    # hist_smooth['weight'] = xsarsea.gradients.circ_smooth(hist_smooth['weight'])
-
-    # select histogram peak
-    iangle = hist_mean_smooth['weight'].fillna(0).argmax(dim='angles')
-    streaks_dir = hist_mean_smooth.angles.isel(angles=iangle)
-    streaks_weight = hist_mean_smooth['weight'].isel(angles=iangle)
-    streaks = xr.merge(
-        [dict(angle=streaks_dir, weight=streaks_weight)]).drop('angles')
-
-    # streaks are [0, pi]. Remove ambiguity with anciallary wind
-    ancillary_wind = xr_dataset_100['ancillary_wind'].sel(line=streaks.line,
-                                                          sample=streaks.sample,
-                                                          method='nearest').compute()
-    streaks_c = streaks['weight'] * np.exp(1j * streaks['angle'])
-    diff_angle = xr.apply_ufunc(np.angle, ancillary_wind / streaks_c)
-    streaks_c = xr.where(np.abs(diff_angle) > np.pi/2, -streaks_c, streaks_c)
-    streaks['weight'] = np.abs(streaks_c)
-    streaks['angle'] = xr.apply_ufunc(np.angle, streaks_c)
-
-    streaks_dir = xr.apply_ufunc(
-        np.angle, streaks_c.interp(line=xr_dataset.line, sample=xr_dataset.sample))
-    streaks_dir = xr.where(
-        xr_dataset['land_mask'], np.nan, streaks_dir)
-    streaks_dir.attrs['comment'] = 'angle in radians, anticlockwise, 0=line'
-    streaks_dir.attrs['description'] = 'wind direction estimated from local gradient, and direction ambiguity removed with ancillary wind'
-    
-    return streaks_dir
-    
-    """

grdwindinversion-0.2.6.dist-info/METADATA

@@ -1,83 +0,0 @@
-Metadata-Version: 2.1
-Name: grdwindinversion
-Version: 0.2.6
-Summary: Package to perform Wind inversion from GRD Level-1 SAR images
-License: MIT
-Classifier: Development Status :: 2 - Pre-Alpha
-Classifier: Intended Audience :: Developers
-Classifier: License :: OSI Approved :: MIT License
-Classifier: Natural Language :: English
-Classifier: Programming Language :: Python :: 3
-Classifier: Programming Language :: Python :: 3.6
-Classifier: Programming Language :: Python :: 3.7
-Classifier: Programming Language :: Python :: 3.8
-Classifier: Programming Language :: Python :: 3.9
-Classifier: Programming Language :: Python :: 3.10
-Classifier: Programming Language :: Python :: 3.11
-Requires-Python: >=3.9
-Description-Content-Type: text/markdown
-License-File: LICENSE
-License-File: AUTHORS.rst
-Requires-Dist: xsar
-Requires-Dist: xsarsea
-Requires-Dist: xarray
-Requires-Dist: xarray-datatree
-Requires-Dist: pyyaml
-Requires-Dist: numpy
-Requires-Dist: scipy
-Requires-Dist: fsspec
-Requires-Dist: aiohttp
-Requires-Dist: xarray-safe-s1
-Requires-Dist: xradarsat2
-Requires-Dist: xarray-safe-rcm
-
-
-# grdwindinversion
-
-
-
-[![Python Version](https://img.shields.io/pypi/pyversions/grdwindinversion.svg)](https://pypi.org/project/grdwindinversion/)
-[![Dependencies Status](https://img.shields.io/badge/dependencies-up%20to%20date-brightgreen.svg)](https://github.com/umr-lops/grdwindinversion/pulls?utf8=%E2%9C%93&q=is%3Apr%20author%3Aapp%2Fdependabot)
-
-
-
-
-
-Package to perform Wind inversion from GRD Level-1 SAR images
-
-
-* Free software: MIT license
-* Documentation: https://grdwindinversion.readthedocs.io.
-
-
-## Usage
-
-
- ```python
-
-    SAR_L1-to-L2_wind_processor -h
-    usage: SAR_L1-to-L2_wind_processor [-h] --input_file INPUT_FILE [--config_file CONFIG_FILE] --outputdir OUTPUTDIR [--verbose] [--overwrite]
-
-    Perform inversion from S1(L1-GRD) SAFE, L1-RCM, L1-RS2 ; using xsar/xsarsea tools
-
-    options:
-      -h, --help            show this help message and exit
-      --input_file INPUT_FILE
-                            input file path
-      --config_file CONFIG_FILE
-                            config file path [if not provided will take config file based on input file]
-      --outputdir OUTPUTDIR
-      --verbose
-      --overwrite           overwrite existing .nc files [default is False]
-```
-
-
-## Features
-
-This Python library (based on `xarray`) allows to perform wind inversion from level-1 GRD (projected magnitude image).
-Mission supported:
- * Sentinel-1
- * RCM
- * RadarSat-2
-
-

grdwindinversion-0.2.6.dist-info/RECORD

@@ -1,19 +0,0 @@
-grdwindinversion/.gitignore,sha256=vmDRHGeESYckUdUztsPr7u6ZNfBYMCVR3GE3AJCBnsY,1204
-grdwindinversion/__init__.py,sha256=VQc2V_j124NX2Gp06Go1oSZDSF4lBIVPYT-W9ZkNBXE,415
-grdwindinversion/config_prod.yaml,sha256=sYF1TAYM5ifFqlGdVwz8ygydo3_6KnEqYb7eBWYeDiw,1047
-grdwindinversion/config_prod_recal.yaml,sha256=WzucOsKVCK-lw5kRn-4llONwBxl8TogYrReeLycXC6k,1006
-grdwindinversion/data_config.yaml,sha256=YbbgxqpgasR5RINO29Ep_1cELdZotoylLzn6Qh7f6LM,473
-grdwindinversion/inversion.py,sha256=CgaOH4uC-iY2pM0bhaLWduKaGarE7ZTVHfSLM0bdRqM,47304
-grdwindinversion/load_config.py,sha256=bt7EXTmXfbuK0oNR-qwtijv_iZBQzi5BtY92RGKEo4Y,651
-grdwindinversion/main.py,sha256=0DLDW0i14CWhDKAGegzOKlWnCkag_cGn3vzMIAKcv4g,2658
-grdwindinversion/streaks.py,sha256=SucALdajwAr-tcqwJs6KOj8CQ7qQPp3JUKSh2dHAJTY,2950
-grdwindinversion/utils.py,sha256=WmHL-NFv3B0PDbCauhBMeqW2ue_340iALmPH1kSd_Pc,2593
-grdwindinversion/utils_memory.py,sha256=1N3Kh4qVZPELPU6I4onbmkur1CZd7EHZqfbSAa6eaVc,1480
-grdwindinversion/.github/ISSUE_TEMPLATE.md,sha256=qiM_a7CCUz3fSrz3Q20Se1nwPNFS8QCc8tkwK_0DSCo,327
-grdwindinversion-0.2.6.dist-info/AUTHORS.rst,sha256=KmhW_5LBKGTIGwWEVkoTm1qx_bvdDR3yYL-1cwbDOFQ,218
-grdwindinversion-0.2.6.dist-info/LICENSE,sha256=-B8mBiTeY3J7OLuayiV1myqmc7yeijBc7s34kc8RTmg,1075
-grdwindinversion-0.2.6.dist-info/METADATA,sha256=_M5G_r0Dc6bu2K-9kmMxbApsmXuHF9Of76IBn-p6CHI,2479
-grdwindinversion-0.2.6.dist-info/WHEEL,sha256=OVMc5UfuAQiSplgO0_WdW7vXVGAt9Hdd6qtN4HotdyA,91
-grdwindinversion-0.2.6.dist-info/entry_points.txt,sha256=2rjvlVCy0iasRXjOz3kOIGuy2OCGQ-VTNuwuViQ6cMM,95
-grdwindinversion-0.2.6.dist-info/top_level.txt,sha256=z6lPix3QPEYOo37qq8plA2hY7S3C8MQZY81agRlksMI,17
-grdwindinversion-0.2.6.dist-info/RECORD,,