Use nanmin/nanmax instead of min/max

Author: claudiodsf

sourcespec/ssp_build_spectra.py

@@ -495,7 +495,7 @@ def _build_weight_from_ratio(spec, specnoise, smooth_width_decades):
     weight.data = np.log10(weight.data)
     # Weight spectrum is smoothed once more
     _smooth_spectrum(weight, smooth_width_decades)
-    weight.data /= np.max(weight.data)
+    weight.data /= np.nanmax(weight.data)
     # slightly taper weight at low frequencies, to avoid overestimating
     # weight at low frequencies, in cases where noise is underestimated
     cosine_taper(
@@ -522,7 +522,7 @@ def _build_weight_from_noise(config, spec, specnoise):
             spec, specnoise, config.spectral_smooth_width_decades)
     # interpolate to log-frequencies
     weight.make_logspaced_from_linear()
-    weight.data_logspaced /= np.max(weight.data_logspaced)
+    weight.data_logspaced /= np.nanmax(weight.data_logspaced)
     # Make sure weight is positive
     weight.data_logspaced[weight.data_logspaced <= 0] = 0.001
     return weight

sourcespec/ssp_util.py

@@ -25,8 +25,8 @@
 # MISC ------------------------------------------------------------------------
 def spec_minmax(amp, freq, amp_minmax=None, freq_minmax=None):
     """Get minimum and maximum values of spectral amplitude and frequency."""
-    amp_min = amp.min()
-    amp_max = amp.max()
+    amp_min = np.nanmin(amp)
+    amp_max = np.nanmax(amp)
     if amp_minmax is None:
         amp_minmax = [amp_min, amp_max]
     else: