rdtools 2.2.0b1__tar.gz → 2.2.0b2__tar.gz

{rdtools-2.2.0b1/rdtools.egg-info → rdtools-2.2.0b2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: rdtools
-Version: 2.2.0b1
+Version: 2.2.0b2
 Summary: Functions for reproducible timeseries analysis of photovoltaic systems.
 Home-page: https://github.com/NREL/rdtools
 Author: Rdtools Python Developers

{rdtools-2.2.0b1 → rdtools-2.2.0b2}/rdtools/_version.py

@@ -8,11 +8,11 @@ import json
 
 version_json = '''
 {
- "date": "2022-12-07T11:38:49-0700",
+ "date": "2023-12-01T15:24:38-0700",
  "dirty": false,
  "error": null,
- "full-revisionid": "adf9f44c9a65b8955ffb95b3d38b0c29e1760ee3",
- "version": "2.2.0-beta.1"
+ "full-revisionid": "250e412bda8199491d8dc45673752374913b9c65",
+ "version": "2.2.0-beta.2"
 }
 '''  # END VERSION_JSON
 

{rdtools-2.2.0b1 → rdtools-2.2.0b2}/rdtools/clearsky_temperature.py

@@ -57,7 +57,9 @@ def get_clearsky_tamb(times, latitude, longitude, window_size=40,
         freq_actual = times.freq
 
     dt_daily = pd.date_range(times.date[0] - buffer, times.date[-1] + buffer,
-                             freq='D', tz=times.tz)
+                             freq='D')
+    dt_daily = dt_daily.tz_localize(times.tz, ambiguous='infer',
+                                    nonexistent='shift_forward')
 
     f = h5py.File(filepath, "r")
 

{rdtools-2.2.0b1 → rdtools-2.2.0b2}/rdtools/degradation.py

@@ -38,7 +38,7 @@ def degradation_ols(energy_normalized, confidence_level=68.2):
 
     # calculate a years column as x value for regression, ignoring leap years
     day_diffs = (df.index - df.index[0])
-    df['days'] = day_diffs.astype('timedelta64[s]') / (60 * 60 * 24)
+    df['days'] = day_diffs / pd.Timedelta('1d')
     df['years'] = df.days / 365.0
 
     # add intercept-constant to the exogeneous variable
@@ -123,22 +123,14 @@ def degradation_classical_decomposition(energy_normalized,
 
     # calculate a years column as x value for regression, ignoring leap years
     day_diffs = (df.index - df.index[0])
-    df['days'] = day_diffs.astype('timedelta64[s]') / (60 * 60 * 24)
+    df['days'] = day_diffs / pd.Timedelta('1d')
     df['years'] = df.days / 365.0
 
     # Compute yearly rolling mean to isolate trend component using
     # moving average
-    it = df.iterrows()
-    energy_ma = []
-    for i, row in it:
-        if row.years - 0.5 >= min(df.years) and \
-           row.years + 0.5 <= max(df.years):
-            roll = df[(df.years <= row.years + 0.5) &
-                      (df.years >= row.years - 0.5)]
-            energy_ma.append(roll.energy_normalized.mean())
-        else:
-            energy_ma.append(np.nan)
-
+    energy_ma = df['energy_normalized'].rolling('365d', center=True).mean()
+    has_full_year = (df['years'] >= df['years'][0] + 0.5) & (df['years'] <= df['years'][-1] - 0.5)
+    energy_ma[~has_full_year] = np.nan
     df['energy_ma'] = energy_ma
 
     # add intercept-constant to the exogeneous variable
@@ -290,7 +282,7 @@ def degradation_year_on_year(energy_normalized, recenter=True,
                        tolerance=pd.Timedelta('8D')
                        )
 
-    df['time_diff_years'] = (df.dt - df.dt_right).astype('timedelta64[h]') / 8760.0
+    df['time_diff_years'] = (df.dt - df.dt_right) / pd.Timedelta('365d')
     df['yoy'] = 100.0 * (df.energy - df.energy_right) / (df.time_diff_years)
     df.index = df.dt
 
@@ -395,10 +387,8 @@ def _mk_test(x, alpha=0.05):
     n = len(x)
 
     # calculate S
-    s = 0
-    for k in range(n - 1):
-        for j in range(k + 1, n):
-            s += np.sign(x[j] - x[k])
+    x = np.array(x)
+    s = np.sum(np.triu(np.sign(-np.subtract.outer(x, x)), 1))
 
     # calculate the unique data
     unique_x = np.unique(x)

{rdtools-2.2.0b1 → rdtools-2.2.0b2}/rdtools/filtering.py

@@ -424,8 +424,9 @@ def logic_clip_filter(power_ac,
     # series sampling frequency is less than 95% consistent.
     _check_data_sampling_frequency(power_ac)
     # Get the sampling frequency of the time series
-    time_series_sampling_frequency = power_ac.index.to_series().diff()\
-        .astype('timedelta64[m]').mode()[0]
+    time_series_sampling_frequency = (
+        power_ac.index.to_series().diff() / pd.Timedelta('60s')
+    ).mode()[0]
     # Make copies of the original inputs for the cases that the data is
     # changes for clipping evaluation
     original_time_series_sampling_frequency = time_series_sampling_frequency
@@ -651,8 +652,7 @@ def xgboost_clip_filter(power_ac,
     # series sampling frequency is less than 95% consistent.
     _check_data_sampling_frequency(power_ac)
     # Get the most common sampling frequency
-    sampling_frequency = int(power_ac.index.to_series().diff()
-                             .astype('timedelta64[m]').mode()[0])
+    sampling_frequency = int((power_ac.index.to_series().diff() / pd.Timedelta('60s')).mode()[0])
     freq_string = str(sampling_frequency) + "T"
     # Min-max normalize
     # Resample the series based on the most common sampling frequency

{rdtools-2.2.0b1 → rdtools-2.2.0b2}/rdtools/normalization.py

@@ -380,7 +380,7 @@ def irradiance_rescale(irrad, irrad_sim, max_iterations=100,
     '''
 
     if method == 'iterative':
-        def _rmse(fact):
+        def _rmse(fact, filt):
             """
             Calculates RMSE with a given rescale fact(or) according to global
             filt(er)
@@ -392,10 +392,9 @@ def irradiance_rescale(irrad, irrad_sim, max_iterations=100,
 
         def _single_rescale(irrad, irrad_sim, guess):
             "Optimizes rescale factor once"
-            global filt
             csi = irrad / (guess * irrad_sim)  # clear sky index
             filt = (csi >= 0.8) & (csi <= 1.2) & (irrad > 200)
-            min_result = minimize(_rmse, guess, method='Nelder-Mead')
+            min_result = minimize(_rmse, guess, (filt), method='Nelder-Mead')
 
             factor = min_result['x'][0]
             return factor

{rdtools-2.2.0b1 → rdtools-2.2.0b2}/rdtools/soiling.py

@@ -1762,11 +1762,11 @@ class CODSAnalysis():
             self.soiling_loss = [0, 0, (1 - result_df.soiling_ratio).mean()]
             self.small_soiling_signal = True
             self.errors = (
-                'Soiling signal is small relative to the noise.'
-                'Iterative decomposition not possible.\n'
-                'Degradation found by RdTools YoY')
-            print(self.errors)
-            return
+                'Soiling signal is small relative to the noise. '
+                'Iterative decomposition not possible. '
+                'Degradation found by RdTools YoY.')
+            warnings.warn(self.errors)
+            return self.result_df, self.degradation, self.soiling_loss
         self.small_soiling_signal = False
 
         # Aggregate all bootstrap samples
@@ -2507,7 +2507,8 @@ def _make_seasonal_samples(list_of_SCs, sample_nr=10, min_multiplier=0.5,
     ''' Generate seasonal samples by perturbing the amplitude and the phase of
         a seasonal components found with the fitted CODS model '''
     samples = pd.DataFrame(index=list_of_SCs[0].index,
-                           columns=range(int(sample_nr*len(list_of_SCs))))
+                           columns=range(int(sample_nr*len(list_of_SCs))),
+                           dtype=float)
     # From each fitted signal, we will generate new seaonal components
     for i, signal in enumerate(list_of_SCs):
         # Remove beginning and end of signal

{rdtools-2.2.0b1 → rdtools-2.2.0b2/rdtools.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: rdtools
-Version: 2.2.0b1
+Version: 2.2.0b2
 Summary: Functions for reproducible timeseries analysis of photovoltaic systems.
 Home-page: https://github.com/NREL/rdtools
 Author: Rdtools Python Developers

{rdtools-2.2.0b1 → rdtools-2.2.0b2}/rdtools.egg-info/requires.txt

@@ -1,12 +1,12 @@
 matplotlib>=3.0.0
-numpy>=1.16
-pandas!=0.24.0,!=1.0.0,!=1.0.1,>=0.23.2
+numpy>=1.17.3
+pandas<2.1,>=1.3.0
 statsmodels>=0.11.1
-scipy>=1.1.0
+scipy>=1.2.0
 h5py>=2.8.0
 plotly>=4.0.0
 xgboost>=1.3.3
-pvlib<0.10.0,>=0.7.0
+pvlib<0.11.0,>=0.7.0
 scikit-learn>=0.22.0
 arch>=4.11
 filterpy>=1.4.2
@@ -17,7 +17,7 @@ flake8
 ipython
 nbsphinx-link==1.3.0
 nbsphinx==0.8.8
-nbval
+nbval==0.9.6
 pytest>=3.6.3
 pytest-mock
 sphinx-gallery==0.8.1
@@ -36,5 +36,5 @@ sphinx-gallery==0.8.1
 pytest>=3.6.3
 coverage
 flake8
-nbval
+nbval==0.9.6
 pytest-mock

{rdtools-2.2.0b1 → rdtools-2.2.0b2}/setup.py

@@ -35,21 +35,23 @@ TESTS_REQUIRE = [
     'pytest >= 3.6.3',
     'coverage',
     'flake8',
-    'nbval',
+    'nbval==0.9.6',  # https://github.com/computationalmodelling/nbval/issues/194
     'pytest-mock',
 ]
 
 INSTALL_REQUIRES = [
     'matplotlib >= 3.0.0',
-    'numpy >= 1.16',
-    # exclude pandas==1.0.0 & 1.0.1 for GH142, and 0.24.0 for GH114
-    'pandas >= 0.23.2,!=0.24.0,!=1.0.0,!=1.0.1',
+    'numpy >= 1.17.3',
+    # pandas restricted to <2.1 until
+    # https://github.com/pandas-dev/pandas/issues/55794
+    # is resolved
+    'pandas >= 1.3.0, <2.1',
     'statsmodels >= 0.11.1',
-    'scipy >= 1.1.0',
+    'scipy >= 1.2.0',
     'h5py >= 2.8.0',
     'plotly>=4.0.0',
     'xgboost >= 1.3.3',
-    'pvlib >= 0.7.0, <0.10.0',
+    'pvlib >= 0.7.0, <0.11.0',
     'scikit-learn >= 0.22.0',
     'arch >= 4.11',
     'filterpy >= 1.4.2'