pyIntensityFeatures 0.1.0__tar.gz → 0.2.0__tar.gz

pyintensityfeatures-0.2.0/.zenodo.json

@@ -0,0 +1,39 @@
+{
+    "license": {
+      "id": "bsd-license"
+    }, 
+    "notes": "When referencing this package, please cite both the package DOI and the pyIntensityFeatures journal article: Burrell, A. G., Chisham, G., Longden, N., Fritz, B., & Zawdie, K. A. (2025). Automated identification of auroral luminosity boundaries using pyIntensityFeatures. Journal of Geophysical Research: Space Physics, 130, e2025JA034230. doi:10.1029/2025JA034230.",
+    "references": [
+      "Burrell, A. G., Chisham, G., Longden, N., Fritz, B., & Zawdie, K. A. (2025). Automated identification of auroral luminosity boundaries using pyIntensityFeatures. Journal of Geophysical Research: Space Physics, 130, e2025JA034230. doi:10.1029/2025JA034230."
+    ],
+    "keywords": [
+  "aurora",
+  "aurora australis",
+  "aurora borealis",
+  "boundaries",
+  "imager",
+  "ionosphere",
+  "heliophysics",
+  "magnetosphere",
+  "northern lights",
+  "intensity",
+  "southern lights",
+  "features"
+    ],
+    "creators": [
+      {
+        "orcid": "0000-0001-8875-9326", 
+        "affiliation": "Naval Research Laboratory", 
+        "name": "Angeline G. Burrell"
+      },
+      {
+        "affiliation": "KPMG", 
+        "name": "Nicola S. Longden"
+      },
+      {
+        "orcid": "0000-0003-1151-5934", 
+        "affiliation": "British Antarctic Survey", 
+        "name": "Gareth Chisham"
+      }
+    ]
+}

pyintensityfeatures-0.2.0/CHANGELOG.md

@@ -0,0 +1,31 @@
+Changelog
+=========
+
+Summary of all changes made since the first release
+
+0.2.0 (2026-03-20)
+------------------
+* Fixed bug in the CI apexpy installation
+* Cycled Python version support and updated NEP 29 test
+* Expanded unit tests
+* Removed support for older SciPy versions
+* Added journal article citation to the documents and .zenodo.json
+* Added missing badges to the README
+* Added DOI to file headers
+* Updated documentation by fixing incorrect description of keywords,
+  updating hyperlinks, and fixing a typo in the fitting example
+* Updated xarray objects to use `sizes` attribute instead of `dims`
+* Added an error catch to remove RuntimeWarning in `gauss` function
+
+0.1.0 (2025-03-28)
+------------------
+* Added option to specify maximum dayglow
+* Updated `as_datetime` to remove timezone information
+* Added a unique function that rounds to a specified significance level
+* Extracted code that identified the parameters for Gaussian fitting into a
+  separate function
+* Added publication release statements
+
+0.0.1 (2024-03-01)
+------------------
+* Alpha Release

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/CODE_OF_CONDUCT.md

@@ -72,4 +72,4 @@ This Code of Conduct is adapted from the [Contributor Covenant][homepage],
 version 1.4, available at [http://contributor-covenant.org/version/1/4][version]
 
 [homepage]: https://contributor-covenant.org
-[version]: https://contributor-covenant.org/version/1/4/
+[version]: https://www.contributor-covenant.org/version/1/4/code-of-conduct/

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pyIntensityFeatures
-Version: 0.1.0
+Version: 0.2.0
 Summary: Identify features, such as auroral luminosity boundaries, in imager data
 Home-page: https://github.com/aburrell/pyIntensityFeatures
 Author-email: Angeline Burrell <angeline.g.burrell.civ@us.navy.mil>
@@ -49,16 +49,17 @@ Classifier: Programming Language :: Python :: 3.10
 Classifier: Programming Language :: Python :: 3.11
 Classifier: Programming Language :: Python :: 3.12
 Classifier: Programming Language :: Python :: 3.13
+Classifier: Programming Language :: Python :: 3.14
 Classifier: Operating System :: POSIX :: Linux
 Classifier: Operating System :: MacOS :: MacOS X
 Classifier: Operating System :: Microsoft :: Windows
-Requires-Python: >=3.9
+Requires-Python: >=3.10
 Description-Content-Type: text/markdown
 License-File: LICENSE
 Requires-Dist: aacgmv2
 Requires-Dist: numpy
 Requires-Dist: pandas
-Requires-Dist: scipy
+Requires-Dist: scipy>=1.9.0
 Requires-Dist: xarray
 Provides-Extra: apex
 Requires-Dist: apexpy; extra == "apex"
@@ -76,8 +77,8 @@ Requires-Dist: sphinx_rtd_theme; extra == "doc"
 Dynamic: license-file
 
 [![Documentation Status](https://readthedocs.org/projects/pyIntensityFeatures/badge/?version=latest)](http://pyIntensityFeatures.readthedocs.io/en/latest/?badge=latest)
-[DOI]
-[PyPI version]
+[![DOI](https://zenodo.org/badge/954923667.svg)](https://doi.org/10.5281/zenodo.15102100)
+[![PyPI version](https://badge.fury.io/py/pyIntensityFeatures.svg)](https://badge.fury.io/py/pyIntensityFeatures)
 [![Test Status](https://github.com/aburrell/pyIntensityFeatures/actions/workflows/main.yml/badge.svg)](https://github.com/aburrell/pyIntensityFeatures/actions/workflows/main.yml)
 [![Coverage Status](https://coveralls.io/repos/github/aburrell/pyIntensityFeatures/badge.svg?branch=main)](https://coveralls.io/github/aburrell/pyIntensityFeatures?branch=main)
 

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/README.md

@@ -1,6 +1,6 @@
 [![Documentation Status](https://readthedocs.org/projects/pyIntensityFeatures/badge/?version=latest)](http://pyIntensityFeatures.readthedocs.io/en/latest/?badge=latest)
-[DOI]
-[PyPI version]
+[![DOI](https://zenodo.org/badge/954923667.svg)](https://doi.org/10.5281/zenodo.15102100)
+[![PyPI version](https://badge.fury.io/py/pyIntensityFeatures.svg)](https://badge.fury.io/py/pyIntensityFeatures)
 [![Test Status](https://github.com/aburrell/pyIntensityFeatures/actions/workflows/main.yml/badge.svg)](https://github.com/aburrell/pyIntensityFeatures/actions/workflows/main.yml)
 [![Coverage Status](https://coveralls.io/repos/github/aburrell/pyIntensityFeatures/badge.svg?branch=main)](https://coveralls.io/github/aburrell/pyIntensityFeatures?branch=main)
 

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/__init__.py

@@ -1,5 +1,6 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
 # Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is
@@ -8,11 +9,7 @@
 """Package to identify intensity features in imager data."""
 
 import logging
-
-try:
-    from importlib import metadata
-except ImportError:
-    import importlib_metadata as metadata
+from importlib import metadata
 
 # Define a logger object to allow easier log handling
 logging.raiseExceptions = False

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/_main.py

@@ -1,5 +1,7 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
+# Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is
 # unlimited.
@@ -374,7 +376,7 @@ class AuroralBounds(object):
         un_threshold : float
             Maximum acceptable uncertainty value in degrees (default=1.25)
         dayglow_threshold : float
-            Minimum allowable background intensity value in Rayleighs
+            Maximum allowable background intensity value in Rayleighs
             (default=300.0)
         strict_fit : bool
             Enforce positive values for the x-offsets in quadratic-Gaussian fits

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/instruments/__init__.py

@@ -1,5 +1,7 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
+# Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is
 # unlimited.

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/instruments/satellites.py

@@ -1,5 +1,7 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
+# Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is
 # unlimited.
@@ -117,10 +119,6 @@ def get_auroral_slice(time_data, glat_data, glon_data, int_data,
             if itime >= len(time_data):
                 itime = -1
 
-            if len(sweep_times) == 0:
-                # If no times were found, insert the start time
-                sweep_times.append(start_time)
-
             sweep_times.append(time_data[itime])
 
         # The sweep has ended, select the desired time indices

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/proc/__init__.py

@@ -1,5 +1,7 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
+# Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is
 # unlimited.

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/proc/boundaries.py

@@ -1,5 +1,7 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
+# Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is
 # unlimited.

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/proc/fitting.py

@@ -1,5 +1,7 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
+# Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is
 # unlimited.
@@ -315,13 +317,9 @@ def get_gaussian_func_fit(mlat_bins, mlt_bins, mean_intensity, std_intensity,
                             gauss_out[fmask])
                         func_params.append(lsq_result[0])
 
-                        # As of scipy 1.11.0 the output changed
-                        if hasattr(pres, 'statistic'):
-                            fit_pearsonr.append(pres.statistic)
-                            fit_pearsonp.append(pres.pvalue)
-                        else:
-                            fit_pearsonr.append(pres[0])
-                            fit_pearsonp.append(pres[1])
+                        # SciPy versions below 1.9.0 will not work
+                        fit_pearsonr.append(pres.statistic)
+                        fit_pearsonp.append(pres.pvalue)
 
                         fit_cov.append(lsq_result[1])
                         found_fit = True

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/proc/intensity.py

@@ -1,5 +1,7 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
+# Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is
 # unlimited.
@@ -49,7 +51,7 @@ def find_intensity_boundaries(intensity, glat, glon, sweep_times, alt,
     un_threshold : float
         Maximum acceptable uncertainty value in degrees (default=1.25)
     dayglow_threshold : float
-        Minimum allowable background intensity value in Rayleighs (default=300)
+        Maximum allowable background intensity value in Rayleighs (default=300)
     strict_fit : bool
         Enforce positive values for the x-offsets in quadratic-Gaussian fits
         (default=False)

pyintensityfeatures-0.2.0/pyIntensityFeatures/tests/__init__.py

@@ -0,0 +1,6 @@
+#!/usr/bin/env python
+# -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
+# Full license can be found in License.md
+# ----------------------------------------------------------------------------
+"""Initialization file for the unit and integration tests."""

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/tests/test_instruments_satellites.py

@@ -1,5 +1,6 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
 # Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/tests/test_main.py

@@ -1,5 +1,6 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
 # Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is
@@ -215,8 +216,8 @@ class TestAuroralBounds(unittest.TestCase):
             end = pyIntensityFeatures.utils.coords.as_datetime(
                 self.inst_data[self.time_var].values[-1])
         else:
-            start = self.inst_data[self.time_var][0]
-            end = self.inst_data[self.time_var][-1]
+            start = list(self.inst_data[self.time_var])[0]
+            end = list(self.inst_data[self.time_var])[-1]
 
         # Evaluate the times
         self.assertTrue(self.alb.stime == start,
@@ -273,7 +274,7 @@ class TestAuroralBounds(unittest.TestCase):
                              msg="unexpected default attributes")
         self.assertDictEqual({'sweep_start': 1, 'mlt': 48, 'coeff': 6,
                               'lat': 31, 'sweep_end': 1},
-                             dict(self.alb.boundaries.dims),
+                             dict(self.alb.boundaries.sizes),
                              msg="unexpected default dimensions")
         self.assertListEqual(["sweep_start", "sweep_end", "mlt", "hemisphere",
                               "lat"],
@@ -552,7 +553,7 @@ class TestAuroralBounds(unittest.TestCase):
                     self.inst_data = self.inst_data[:-10]
                 else:
                     self.inst_data = xr.Dataset({
-                        var: (self.inst_data[var].dims,
+                        var: (self.inst_data[var].sizes,
                               self.inst_data[var].values[:-10])
                         for var in [self.time_var, self.glat_var, self.glon_var,
                                     self.intensity_var, 'clean_flag']})

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/tests/test_proc_boundaries.py

@@ -1,5 +1,6 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
 # Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/tests/test_proc_fitting.py

@@ -1,5 +1,6 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
 # Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/tests/test_proc_intensity.py

@@ -1,5 +1,6 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
 # Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is
@@ -203,3 +204,53 @@ class TestIntensityFuncs(unittest.TestCase):
                                        intensity.find_intensity_boundaries,
                                        *args)
         return
+
+    def test_find_intensity_bad_eq_boundaries(self):
+        """Test boundary ID in an equatorial boundary failure."""
+        # Change intensity to have a peak where the equatorial boundary will
+        # not be defined
+        self.intensity = mult_gauss_quad(self.glat, *[
+            1.0, 0.01, 0.001, 500.0, 60.0, 3.0])
+
+        for hemi in [-1, 1]:
+            with self.subTest(hemi=hemi):
+                # Get the boundary outputs
+                (self.sweep_end, self.out_data,
+                 self.out_coeff) = intensity.find_intensity_boundaries(
+                     self.intensity, hemi * self.glat, self.glon,
+                     self.sweep_times, self.alt, self.min_mlat_base,
+                     self.max_coeff, mlat_inc=self.mlat_inc,
+                     mlt_inc=self.mlt_inc, strict_fit=self.strict_fit)
+
+                # Evaluate the outputs
+                self.assertEqual(self.sweep_end, self.sweep_times[-1])
+                self.assertLessEqual(self.out_coeff, 12)
+                self.assertGreater(self.out_coeff, 0)
+                self.assertTrue(np.isnan(self.out_data['eq_bounds']).all())
+                self.assertFalse(np.isnan(self.out_data['po_bounds']).all())
+
+        return
+
+    def test_find_intensity_bad_po_boundaries(self):
+        """Test boundary ID in a polar boundary failure."""
+        # Change intensity to have a peak where the polar boundary will not
+        # be defined in the southern hemisphere
+        self.intensity = mult_gauss_quad(self.glat, *[
+            -1.0, -0.1, 0.01, 40.0, 77.0, 5.0])
+
+        # Get the boundary outputs
+        (self.sweep_end, self.out_data,
+         self.out_coeff) = intensity.find_intensity_boundaries(
+             self.intensity, -1 * self.glat, self.glon,
+             self.sweep_times, self.alt, self.min_mlat_base,
+             self.max_coeff, mlat_inc=self.mlat_inc,
+             mlt_inc=self.mlt_inc, strict_fit=self.strict_fit)
+
+        # Evaluate the outputs
+        self.assertEqual(self.sweep_end, self.sweep_times[-1])
+        self.assertLessEqual(self.out_coeff, 12)
+        self.assertGreater(self.out_coeff, 0)
+        self.assertFalse(np.isnan(self.out_data['eq_bounds']).all())
+        self.assertTrue(np.isnan(self.out_data['po_bounds']).all())
+
+        return

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/tests/test_utils_checks.py

@@ -1,5 +1,6 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
 # Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is
@@ -892,13 +893,14 @@ class TestEvalBoundariesMLT(unittest.TestCase):
         """Test raises ValueError with bad MLT coordinates."""
         # Update the MLT coordinates
         self.mlt[0] += 0.01
-        self.bound_data.assign({self.lt_key: self.mlt})
+        self.bound_data = self.bound_data.assign({self.lt_key: self.mlt})
 
         # Run and check the error
         self.assertRaisesRegex(ValueError, "not have a fixed frequency",
                                checks.evaluate_boundary_in_mlt,
                                *[self.bound_data, self.eq_key, self.po_key,
                                  self.lt_key, self.ut_key])
+
         return
 
     def test_eval_success_no_removal(self):

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/tests/test_utils_coords.py

@@ -1,5 +1,6 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
 # Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/tests/test_utils_distributions.py

@@ -1,5 +1,6 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
 # Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/tests/test_utils_grids.py

@@ -1,5 +1,6 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
 # Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/tests/test_utils_output.py

@@ -1,5 +1,6 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
 # Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is
@@ -121,7 +122,7 @@ class TestBoundaryDictFuncs(unittest.TestCase):
         for coord in self.coord_dict.keys():
             self.assertIn(coord, self.dataset.coords)
 
-            if coord in self.dataset.dims:
+            if coord in self.dataset.sizes:
                 if coord.find('sweep') < 0:
                     self.assertTrue(
                         np.all(self.dataset.coords[coord].values
@@ -399,7 +400,7 @@ class TestBoundaryDictFuncs(unittest.TestCase):
                                                     self.data_dict, 0)
 
         # Evalute the empty dataset
-        self.assertEqual(len(self.dataset.dims), 0)
+        self.assertEqual(len(self.dataset.sizes), 0)
         self.assertEqual(len(self.dataset.data_vars), 0)
         return
 

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/utils/__init__.py

@@ -1,5 +1,7 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
+# Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is
 # unlimited.

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/utils/checks.py

@@ -1,5 +1,7 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
+# Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is
 # unlimited.
@@ -194,7 +196,7 @@ def evaluate_dayglow(fit_coeff, locations, thresh=300.0):
     locations : array-like
         Locations at which the dayglow level will be evaluated.
     thresh : float
-        Minimum allowable background intensity value in Rayleighs (default=300)
+        Maximum allowable background intensity value in Rayleighs (default=300)
 
     Returns
     -------

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/utils/coords.py

@@ -1,5 +1,7 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
+# Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is
 # unlimited.

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/utils/distributions.py

@@ -1,5 +1,7 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
+# Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is
 # unlimited.
@@ -57,7 +59,10 @@ def gauss(x, amp, mu, sigma, const):
         Normal value at `x`
 
     """
-    y = amp * np.exp(-np.power((x - mu), 2.0) / (2.0 * sigma**2)) + const
+    if abs(sigma) > 0.0:
+        y = amp * np.exp(-np.power((x - mu), 2.0) / (2.0 * sigma**2)) + const
+    else:
+        y = np.zeros(shape=x.shape, dtype=float)
 
     return y
 

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/utils/grids.py

@@ -1,5 +1,7 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
+# Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is
 # unlimited.

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures/utils/output.py

@@ -1,5 +1,7 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
+# DOI: 10.5281/zenodo.15102100
+# Full license can be found in License.md
 #
 # DISTRIBUTION STATEMENT A: Approved for public release. Distribution is
 # unlimited.

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pyIntensityFeatures
-Version: 0.1.0
+Version: 0.2.0
 Summary: Identify features, such as auroral luminosity boundaries, in imager data
 Home-page: https://github.com/aburrell/pyIntensityFeatures
 Author-email: Angeline Burrell <angeline.g.burrell.civ@us.navy.mil>
@@ -49,16 +49,17 @@ Classifier: Programming Language :: Python :: 3.10
 Classifier: Programming Language :: Python :: 3.11
 Classifier: Programming Language :: Python :: 3.12
 Classifier: Programming Language :: Python :: 3.13
+Classifier: Programming Language :: Python :: 3.14
 Classifier: Operating System :: POSIX :: Linux
 Classifier: Operating System :: MacOS :: MacOS X
 Classifier: Operating System :: Microsoft :: Windows
-Requires-Python: >=3.9
+Requires-Python: >=3.10
 Description-Content-Type: text/markdown
 License-File: LICENSE
 Requires-Dist: aacgmv2
 Requires-Dist: numpy
 Requires-Dist: pandas
-Requires-Dist: scipy
+Requires-Dist: scipy>=1.9.0
 Requires-Dist: xarray
 Provides-Extra: apex
 Requires-Dist: apexpy; extra == "apex"
@@ -76,8 +77,8 @@ Requires-Dist: sphinx_rtd_theme; extra == "doc"
 Dynamic: license-file
 
 [![Documentation Status](https://readthedocs.org/projects/pyIntensityFeatures/badge/?version=latest)](http://pyIntensityFeatures.readthedocs.io/en/latest/?badge=latest)
-[DOI]
-[PyPI version]
+[![DOI](https://zenodo.org/badge/954923667.svg)](https://doi.org/10.5281/zenodo.15102100)
+[![PyPI version](https://badge.fury.io/py/pyIntensityFeatures.svg)](https://badge.fury.io/py/pyIntensityFeatures)
 [![Test Status](https://github.com/aburrell/pyIntensityFeatures/actions/workflows/main.yml/badge.svg)](https://github.com/aburrell/pyIntensityFeatures/actions/workflows/main.yml)
 [![Coverage Status](https://coveralls.io/repos/github/aburrell/pyIntensityFeatures/badge.svg?branch=main)](https://coveralls.io/github/aburrell/pyIntensityFeatures?branch=main)
 

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyIntensityFeatures.egg-info/requires.txt

@@ -1,7 +1,7 @@
 aacgmv2
 numpy
 pandas
-scipy
+scipy>=1.9.0
 xarray
 
 [apex]

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/pyproject.toml

@@ -4,10 +4,10 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "pyIntensityFeatures"
-version = "0.1.0"
+version = "0.2.0"
 description = 'Identify features, such as auroral luminosity boundaries, in imager data'
 readme = "README.md"
-requires-python = ">=3.9"
+requires-python = ">=3.10"
 license = {file = "LICENSE"}
 authors = [
     {name = "Angeline Burrell", email = "angeline.g.burrell.civ@us.navy.mil"}
@@ -24,6 +24,7 @@ classifiers = [
   "Programming Language :: Python :: 3.11",
   "Programming Language :: Python :: 3.12",
   "Programming Language :: Python :: 3.13",
+  "Programming Language :: Python :: 3.14",
   "Operating System :: POSIX :: Linux",
   "Operating System :: MacOS :: MacOS X",
   "Operating System :: Microsoft :: Windows"
@@ -46,7 +47,7 @@ dependencies = [
   "aacgmv2",
   "numpy",
   "pandas",
-  "scipy",
+  "scipy>=1.9.0",
   "xarray"]
 
 [project.optional-dependencies]

{pyintensityfeatures-0.1.0 → pyintensityfeatures-0.2.0}/setup.cfg

@@ -1,6 +1,6 @@
 [metadata]
 name = pyIntensityFeatures
-version = 0.1.0
+version = 0.2.0
 url = https://github.com/aburrell/pyIntensityFeatures
 
 [flake8]

pyintensityfeatures-0.1.0/.zenodo.json

@@ -1,35 +0,0 @@
-{
-    "license": {
-      "id": "bsd-license"
-    },
-    "keywords": [
-  "aurora",
-  "aurora australis",
-  "aurora borealis",
-  "boundaries",
-  "imager",
-  "ionosphere",
-  "heliophysics",
-  "magnetosphere",
-  "northern lights",
-  "intensity",
-  "southern lights",
-  "features"
-    ],
-    "creators": [
-      {
-        "orcid": "0000-0001-8875-9326", 
-        "affiliation": "Naval Research Laboratory", 
-        "name": "Angeline G. Burrell"
-      },
-      {
-        "affiliation": "KPMG", 
-        "name": "Nicola S. Longden"
-      },
-      {
-        "orcid": "0000-0003-1151-5934", 
-        "affiliation": "British Antarctic Survey", 
-        "name": "Gareth Chisham"
-      }
-    ]
-}

pyintensityfeatures-0.1.0/CHANGELOG.md

@@ -1,17 +0,0 @@
-Changelog
-=========
-
-Summary of all changes made since the first release
-
-0.1.0 (2025-03-28)
-------------------
-* Added option to specify maximum dayglow
-* Updated `as_datetime` to remove timezone information
-* Added a unique function that rounds to a specified significance level
-* Extracted code that identified the parameters for Gaussian fitting into a
-  separate function
-* Added publication release statements
-
-0.0.1 (2024-03-01)
-------------------
-* Alpha Release

pyintensityfeatures-0.1.0/pyIntensityFeatures/tests/__init__.py

@@ -1 +0,0 @@
-"""Initialization file for the unit and integration tests."""