hyperspy._signals.signal1d module¶

class
hyperspy._signals.signal1d.
LazySignal1D
(*args, **kwargs)¶ Bases:
hyperspy._signals.lazy.LazySignal
,hyperspy._signals.signal1d.Signal1D

class
hyperspy._signals.signal1d.
Signal1D
(*args, **kwargs)¶ Bases:
hyperspy.signal.BaseSignal
,hyperspy._signals.common_signal1d.CommonSignal1D

align1D
(start=None, end=None, reference_indices=None, max_shift=None, interpolate=True, number_of_interpolation_points=5, interpolation_method='linear', crop=True, expand=False, fill_value=nan, also_align=None, mask=None, show_progressbar=None)¶ Estimate the shifts in the signal axis using crosscorrelation and use the estimation to align the data in place. This method can only estimate the shift by comparing unidimensional features that should not change the position.
To decrease memory usage, time of computation and improve accuracy it is convenient to select the feature of interest setting the start and end keywords. By default interpolation is used to obtain subpixel precision.
Parameters:  end (start,) – The limits of the interval. If int they are taken as the axis index. If float they are taken as the axis value.
 reference_indices (tuple of ints or None) – Defines the coordinates of the spectrum that will be used as eference. If None the spectrum at the current coordinates is used for this purpose.
 max_shift (int) – “Saturation limit” for the shift.
 interpolate (bool) – If True, interpolation is used to provide subpixel accuracy.
 number_of_interpolation_points (int) – Number of interpolation points. Warning: making this number too big can saturate the memory
 interpolation_method (str or int) – Specifies the kind of interpolation as a string (‘linear’, ‘nearest’, ‘zero’, ‘slinear’, ‘quadratic, ‘cubic’) or as an integer specifying the order of the spline interpolator to use.
 crop (bool) – If True automatically crop the signal axis at both ends if needed.
 expand (bool) – If True, the data will be expanded to fit all data after alignment. Overrides crop.
 fill_value (float) – If crop is False fill the data outside of the original interval with the given value where needed.
 also_align (list of signals, None) – A list of BaseSignal instances that has exactly the same dimensions as this one and that will be aligned using the shift map estimated using the this signal.
 mask (BaseSignal of bool data type.) – It must have signal_dimension = 0 and navigation_shape equal to the current signal. Where mask is True the shift is not computed and set to nan.
 show_progressbar (None or bool) – If True, display a progress bar. If None the default is set in preferences.
Returns: Return type: An array with the result of the estimation. The shift will be
Raises: SignalDimensionError if the signal dimension is not 1.
See also

calibrate
(display=True, toolkit=None)¶ Calibrate the spectral dimension using a gui. It displays a window where the new calibration can be set by: * Setting the offset, units and scale directly * Selection a range by dragging the mouse on the spectrum figure
andsetting the new values for the given range limits
display: bool
If True, display the user interface widgets. If False, return the widgets container in a dictionary, usually for customisation or testing.
toolkit: str, iterable of strings or NoneIf None (default), all available widgets are displayed or returned. If string, only the widgets of the selected toolkit are displayed if available. If an interable of toolkit strings, the widgets of all listed toolkits are displayed or returned.
For this method to work the output_dimension must be 1. Set the view accordingly Raises —— SignalDimensionError if the signal dimension is not 1.

create_model
(dictionary=None)¶ Create a model for the current data.
Returns: model Return type: Model1D instance.

crop_signal1D
(*args, **kwargs)¶ Crop in place the spectral dimension.
Parameters: righ_value (left_value,) – If int the values are taken as indices. If float they are converted to indices using the spectral axis calibration. If left_value is None crops from the beginning of the axis. If right_value is None crops up to the end of the axis. If both are None the interactive cropping interface is activated enabling cropping the spectrum using a span selector in the signal plot. Raises: SignalDimensionError if the signal dimension is not 1.

estimate_peak_width
(factor=0.5, window=None, return_interval=False, parallel=None, show_progressbar=None)¶ Estimate the width of the highest intensity of peak of the spectra at a given fraction of its maximum.
It can be used with asymmetric peaks. For accurate results any background must be previously substracted. The estimation is performed by interpolation using cubic splines.
Parameters:  factor (0 < float < 1) – The default, 0.5, estimates the FWHM.
 window (None, float) – The size of the window centred at the peak maximum used to perform the estimation. The window size must be chosen with care: if it is narrower than the width of the peak at some positions or if it is so wide that it includes other more intense peaks this method cannot compute the width and a NaN is stored instead.
 return_interval (bool) – If True, returns 2 extra signals with the positions of the desired height fraction at the left and right of the peak.
 parallel ({None, bool}) –
 show_progressbar (None or bool) – If True, display a progress bar. If None the default is set in preferences.
Returns:  width or [width, left, right], depending on the value of
 return_interval.

estimate_shift1D
(start=None, end=None, reference_indices=None, max_shift=None, interpolate=True, number_of_interpolation_points=5, mask=None, parallel=None, show_progressbar=None)¶  Estimate the shifts in the current signal axis using
 crosscorrelation.
This method can only estimate the shift by comparing unidimensional features that should not change the position in the signal axis. To decrease the memory usage, the time of computation and the accuracy of the results it is convenient to select the feature of interest providing sensible values for start and end. By default interpolation is used to obtain subpixel precision. :param start, end: The limits of the interval. If int they are taken as the
axis index. If float they are taken as the axis value.Parameters:  reference_indices (tuple of ints or None) – Defines the coordinates of the spectrum that will be used as eference. If None the spectrum at the current coordinates is used for this purpose.
 max_shift (int) – “Saturation limit” for the shift.
 interpolate (bool) – If True, interpolation is used to provide subpixel accuracy.
 number_of_interpolation_points (int) – Number of interpolation points. Warning: making this number too big can saturate the memory
 mask (BaseSignal of bool data type.) – It must have signal_dimension = 0 and navigation_shape equal to the current signal. Where mask is True the shift is not computed and set to nan.
 parallel ({None, bool}) –
 show_progressbar (None or bool) – If True, display a progress bar. If None the default is set in preferences.
Returns:  An array with the result of the estimation in the axis units. although
 the computation is performed in batches if the signal is lazy, the
 result is computed in memory because it depends on the current state
 of the axes that could change later on in the workflow.
Raises: SignalDimensionError if the signal dimension is not 1.

filter_butterworth
(cutoff_frequency_ratio=None, type='low', order=2, display=True, toolkit=None)¶ Butterworth filter in place.
display: bool
If True, display the user interface widgets. If False, return the widgets container in a dictionary, usually for customisation or testing.
toolkit: str, iterable of strings or NoneIf None (default), all available widgets are displayed or returned. If string, only the widgets of the selected toolkit are displayed if available. If an interable of toolkit strings, the widgets of all listed toolkits are displayed or returned.
SignalDimensionError if the signal dimension is not 1.

find_peaks1D_ohaver
(xdim=None, slope_thresh=0, amp_thresh=None, subchannel=True, medfilt_radius=5, maxpeakn=30000, peakgroup=10, parallel=None)¶ Find peaks along a 1D line (peaks in spectrum/spectra).
Function to locate the positive peaks in a noisy xy data set.
Detects peaks by looking for downward zerocrossings in the first derivative that exceed ‘slope_thresh’.
Returns an array containing position, height, and width of each peak.
‘slope_thresh’ and ‘amp_thresh’, control sensitivity: higher values will neglect smaller features.
peakgroup is the number of points around the top peak to search around
Parameters:  slope_thresh (float (optional)) – 1st derivative threshold to count the peak default is set to 0.5 higher values will neglect smaller features.
 amp_thresh (float (optional)) – intensity threshold above which default is set to 10% of max(y) higher values will neglect smaller features.
 medfilt_radius (int (optional)) – median filter window to apply to smooth the data (see scipy.signal.medfilt) if 0, no filter will be applied. default is set to 5
 peakgroup (int (optional)) – number of points around the “top part” of the peak default is set to 10
 maxpeakn (int (optional)) – number of maximum detectable peaks default is set to 5000
 subpix (bool (optional)) – default is set to True
 parallel ({None, bool}) – Perform the operation in a threaded (parallel) manner.
Returns:  peaks (structured array of shape _navigation_shape_in_array in which)
 each cell contains an array that contains as many structured arrays as
 peaks where found at that location and which fields (position, height,)
 width, contains position, height, and width of each peak.
Raises: SignalDimensionError if the signal dimension is not 1.

gaussian_filter
(FWHM)¶ Applies a Gaussian filter in the spectral dimension in place.
Parameters: FWHM (float) – The Full Width at Half Maximum of the gaussian in the spectral axis units
Raises:  ValueError if FWHM is equal or less than zero.
 SignalDimensionError if the signal dimension is not 1.

hanning_taper
(side='both', channels=None, offset=0)¶ Apply a hanning taper to the data in place.
Parameters:  side ({'left', 'right', 'both'}) –
 channels ({None, int}) – The number of channels to taper. If None 5% of the total number of channels are tapered.
 offset (int) –
Returns: Return type: channels
Raises: SignalDimensionError if the signal dimension is not 1.

integrate_in_range
(signal_range='interactive', display=True, toolkit=None)¶ Sums the spectrum over an energy range, giving the integrated area. The energy range can either be selected through a GUI or the command line. :param signal_range: l and r are the left and right limits of the range. They can be
numbers or None, where None indicates the extremes of the interval. If l and r are floats the signal_range will be in axis units (for example eV). If l and r are integers the signal_range will be in index units. When signal_range is “interactive” (default) the range is selected using a GUI.Returns: integrated_spectrum Return type: BaseSignal subclass See also
integrate_simpson()
Examples
Using the GUI >>> s = hs.signals.Signal1D(range(1000)) >>> s.integrate_in_range() #doctest: +SKIP
Using the CLI >>> s_int = s.integrate_in_range(signal_range=(560,None))
Selecting a range in the axis units, by specifying the signal range with floats. >>> s_int = s.integrate_in_range(signal_range=(560.,590.))
Selecting a range using the index, by specifying the signal range with integers. >>> s_int = s.integrate_in_range(signal_range=(100,120))

interpolate_in_between
(start, end, delta=3, parallel=None, show_progressbar=None, **kwargs)¶ Replace the data in a given range by interpolation. The operation is performed in place. :param start, end: The limits of the interval. If int they are taken as the
axis index. If float they are taken as the axis value.Parameters:  delta ({int  float}) – The windows around the (start, end) to use for interpolation
 show_progressbar (None or bool) – If True, display a progress bar. If None the default is set in preferences.
 parallel ({None, bool}) –
 extra keyword arguments are passed to (All) –
 See the function documentation (scipy.interpolate.interp1d.) –
 details. (for) –
Raises: SignalDimensionError if the signal dimension is not 1.

remove_background
(signal_range='interactive', background_type='Power Law', polynomial_order=2, fast=True, show_progressbar=None, display=True, toolkit=None)¶ Remove the background, either in place using a gui or returned as a new spectrum using the command line. Parameters ——— signal_range : tuple, optional
If this argument is not specified, the signal range has to be selected using a GUI. And the original spectrum will be replaced. If tuple is given, the a spectrum will be returned. background_type : string
 The type of component which should be used to fit the background. Possible components: PowerLaw, Gaussian, Offset, Polynomial If Polynomial is used, the polynomial order can be specified
 polynomial_order : int, default 2
 Specify the polynomial order if a Polynomial background is used.
 fast : bool
 If True, perform an approximative estimation of the parameters. If False, the signal is fitted using nonlinear least squares afterwards.This is slower compared to the estimation but possibly more accurate.
 show_progressbar : None or bool
 If True, display a progress bar. If None the default is set in preferences.
display: bool
If True, display the user interface widgets. If False, return the widgets container in a dictionary, usually for customisation or testing.
toolkit: str, iterable of strings or NoneIf None (default), all available widgets are displayed or returned. If string, only the widgets of the selected toolkit are displayed if available. If an interable of toolkit strings, the widgets of all listed toolkits are displayed or returned.
Using gui, replaces spectrum s >>> s = hs.signals.Signal1D(range(1000)) >>> s.remove_background() #doctest: +SKIP
Using command line, returns a spectrum >>> s1 = s.remove_background(signal_range=(400,450), background_type=’PowerLaw’)
Using a full model to fit the background >>> s1 = s.remove_background(signal_range=(400,450), fast=False)
SignalDimensionError if the signal dimension is not 1.

shift1D
(shift_array, interpolation_method='linear', crop=True, expand=False, fill_value=nan, parallel=None, show_progressbar=None)¶ Shift the data in place over the signal axis by the amount specified by an array.
Parameters:  shift_array (numpy array) – An array containing the shifting amount. It must have axes_manager._navigation_shape_in_array shape.
 interpolation_method (str or int) – Specifies the kind of interpolation as a string (‘linear’, ‘nearest’, ‘zero’, ‘slinear’, ‘quadratic, ‘cubic’) or as an integer specifying the order of the spline interpolator to use.
 crop (bool) – If True automatically crop the signal axis at both ends if needed.
 expand (bool) – If True, the data will be expanded to fit all data after alignment. Overrides crop.
 fill_value (float) – If crop is False fill the data outside of the original interval with the given value where needed.
 parallel ({None, bool}) –
 show_progressbar (None or bool) – If True, display a progress bar. If None the default is set in preferences.
Raises: SignalDimensionError if the signal dimension is not 1.

smooth_lowess
(smoothing_parameter=None, number_of_iterations=None, show_progressbar=None, parallel=None, display=True, toolkit=None)¶ Lowess data smoothing in place. If smoothing_parameter or number_of_iterations are None the method is run in interactive mode. Parameters ——— smoothing_parameter: float or None
Between 0 and 1. The fraction of the data used when estimating each yvalue. number_of_iterations: int or None
 The number of residualbased reweightings to perform.
 show_progressbar : None or bool
 If True, display a progress bar. If None the default is set in preferences.
 parallel : {Bool, None, int}
 Perform the operation parallel
display: bool
If True, display the user interface widgets. If False, return the widgets container in a dictionary, usually for customisation or testing.
toolkit: str, iterable of strings or NoneIf None (default), all available widgets are displayed or returned. If string, only the widgets of the selected toolkit are displayed if available. If an interable of toolkit strings, the widgets of all listed toolkits are displayed or returned.
SignalDimensionError if the signal dimension is not 1. ImportError if statsmodels is not installed. Notes —– This method uses the lowess algorithm from statsmodels. statsmodels is required for this method.

smooth_savitzky_golay
(polynomial_order=None, window_length=None, differential_order=0, parallel=None, display=True, toolkit=None)¶ Apply a SavitzkyGolay filter to the data in place. If polynomial_order or window_length or differential_order are None the method is run in interactive mode. Parameters ——— window_length : int
The length of the filter window (i.e. the number of coefficients). window_length must be a positive odd integer. polynomial_order : int
 The order of the polynomial used to fit the samples. polyorder must be less than window_length.
 differential_order: int, optional
 The order of the derivative to compute. This must be a nonnegative integer. The default is 0, which means to filter the data without differentiating.
 parallel : {bool, None}
 Perform the operation in a threaded manner (parallely).
display: bool
If True, display the user interface widgets. If False, return the widgets container in a dictionary, usually for customisation or testing.
toolkit: str, iterable of strings or NoneIf None (default), all available widgets are displayed or returned. If string, only the widgets of the selected toolkit are displayed if available. If an interable of toolkit strings, the widgets of all listed toolkits are displayed or returned.
More information about the filter in scipy.signal.savgol_filter.

smooth_tv
(smoothing_parameter=None, show_progressbar=None, parallel=None, display=True, toolkit=None)¶ Total variation data smoothing in place. Parameters ——— smoothing_parameter: float or None
Denoising weight relative to L2 minimization. If None the method is run in interactive mode. show_progressbar : None or bool
 If True, display a progress bar. If None the default is set in preferences.
 parallel : {Bool, None, int}
 Perform the operation parallely
display: bool
If True, display the user interface widgets. If False, return the widgets container in a dictionary, usually for customisation or testing.
toolkit: str, iterable of strings or NoneIf None (default), all available widgets are displayed or returned. If string, only the widgets of the selected toolkit are displayed if available. If an interable of toolkit strings, the widgets of all listed toolkits are displayed or returned.
SignalDimensionError if the signal dimension is not 1.

spikes_removal_tool
(signal_mask=None, navigation_mask=None, display=True, toolkit=None)¶ Graphical interface to remove spikes from EELS spectra.
Parameters:  signal_mask (boolean array) – Restricts the operation to the signal locations not marked as True (masked)
 navigation_mask (boolean array) – Restricts the operation to the navigation locations not marked as True (masked)
 display (bool) – If True, display the user interface widgets. If False, return the widgets container in a dictionary, usually for customisation or testing.
 toolkit (str, iterable of strings or None) – If None (default), all available widgets are displayed or returned. If string, only the widgets of the selected toolkit are displayed if available. If an interable of toolkit strings, the widgets of all listed toolkits are displayed or returned.
See also
_spikes_diagnosis()


hyperspy._signals.signal1d.
find_peaks_ohaver
(y, x=None, slope_thresh=0.0, amp_thresh=None, medfilt_radius=5, maxpeakn=30000, peakgroup=10, subchannel=True)¶ Find peaks along a 1D line. Function to locate the positive peaks in a noisy xy data set. Detects peaks by looking for downward zerocrossings in the first derivative that exceed ‘slope_thresh’. Returns an array containing position, height, and width of each peak. Sorted by position. ‘slope_thresh’ and ‘amp_thresh’, control sensitivity: higher values will neglect smaller features. :param y: 1D input array, e.g. a spectrum :type y: array :param x: 1D array describing the calibration of y (must have same shape as y) :type x: array (optional) :param slope_thresh: 1st derivative threshold to count the peak
default is set to 0.5 higher values will neglect smaller features.Parameters:  amp_thresh (float (optional)) – intensity threshold above which default is set to 10% of max(y) higher values will neglect smaller features.
 medfilt_radius (int (optional)) – median filter window to apply to smooth the data (see scipy.signal.medfilt) if 0, no filter will be applied. default is set to 5
 peakgroup (int (optional)) – number of points around the “top part” of the peak default is set to 10
 maxpeakn (int (optional)) – number of maximum detectable peaks default is set to 30000
 subchannel (bool (optional)) – default is set to True
Returns: P – contains position, height, and width of each peak
Return type: structured array of shape (npeaks) and fields: position, width, height
Examples
>>> x = np.arange(0,50,0.01) >>> y = np.cos(x) >>> peaks = find_peaks_ohaver(y, x, 0, 0)
Notes
Original code from T. C. O’Haver, 1995. Version 2 Last revised Oct 27, 2006 Converted to Python by Michael Sarahan, Feb 2011. Revised to handle edges better. MCS, Mar 2011

hyperspy._signals.signal1d.
interpolate1D
(number_of_interpolation_points, data)¶