hyperspy._signals package¶
Submodules¶
hyperspy._signals.dielectric_function module¶

class
hyperspy._signals.dielectric_function.
DielectricFunction
(*args, **kwards)¶ Bases:
hyperspy._signals.signal1d.Signal1D

get_electron_energy_loss_spectrum
(zlp, t)¶

get_number_of_effective_electrons
(nat, cumulative=False)¶  Compute the number of effective electrons using the Bethe fsum
rule.
The Bethe fsum rule gives rise to two definitions of the effective number (see [Egerton2011]): $n_{mathrm{eff}}left(Imleft(epsilon^{1}
ight) ight)$ that
we’ll call neff1 and $n_{mathrm{eff}}left(epsilon_{2} ight)$ that we’ll call neff2. This
method computes both.
 nat: float
 Number of atoms (or molecules) per unit volume of the sample.
 cumulative : bool
 If False calculate the number of effective electrons up to the higher energyloss of the spectrum. If True, calculate the number of effective electrons as a function of the energyloss up to the higher energyloss of the spectrum. True is only supported by SciPy newer than 0.13.2.
 neff1, neff2: Signal
 Signal instances containing neff1 and neff2. The signal and navigation dimensions are the same as the current signal if cumulative is True, otherwise the signal dimension is 0 and the navigation dimension is the same as the current signal.
[Egerton2011] Ray Egerton, “Electron EnergyLoss Spectroscopy in the Electron Microscope”, SpringerVerlag, 2011.

hyperspy._signals.eds module¶

class
hyperspy._signals.eds.
EDSSpectrum
(*args, **kwards)¶ Bases:
hyperspy.signals.Spectrum

add_elements
(elements)¶ Add elements and the corresponding Xray lines.
The list of elements is stored in metadata.Sample.elements
Parameters: elements (list of strings) – The symbol of the elements. Examples
>>> s = hs.datasets.example_signals.EDS_SEM_Spectrum() >>> print s.metadata.Sample.elements >>> s.add_elements(['Ar']) >>> print s.metadata.Sample.elements ['Al' 'C' 'Cu' 'Mn' 'Zr'] ['Al', 'Ar', 'C', 'Cu', 'Mn', 'Zr']
See also

add_lines
(lines=(), only_one=True, only_lines=('a', ))¶ Add Xrays lines to the internal list.
Although most functions do not require an internal list of Xray lines because they can be calculated from the internal list of elements, ocassionally it might be useful to customize the Xray lines to be use by all functions by default using this method. The list of Xray lines is stored in metadata.Sample.xray_lines
Parameters:  lines (list of strings) – A list of valid element Xray lines to add e.g. Fe_Kb. Additionally, if metadata.Sample.elements is defined, add the lines of those elements that where not given in this list. If the list is empty (default), and metadata.Sample.elements is defined, add the lines of all those elements.
 only_one (bool) – If False, add all the lines of each element in metadata.Sample.elements that has not line defined in lines. If True (default), only add the line at the highest energy above an overvoltage of 2 (< beam energy / 2).
 only_lines ({None, list of strings}) – If not None, only the given lines will be added.
Examples
>>> s = hs.datasets.example_signals.EDS_SEM_Spectrum() >>> s.add_lines() >>> print s.metadata.Sample.xray_lines ['Al_Ka', 'C_Ka', 'Cu_La', 'Mn_La', 'Zr_La']
>>> s = hs.datasets.example_signals.EDS_SEM_Spectrum() >>> s.set_microscope_parameters(beam_energy=30) >>> s.add_lines() >>> print s.metadata.Sample.xray_lines ['Al_Ka', 'C_Ka', 'Cu_Ka', 'Mn_Ka', 'Zr_La']
>>> s = hs.datasets.example_signals.EDS_SEM_Spectrum() >>> s.add_lines() >>> print s.metadata.Sample.xray_lines >>> s.add_lines(['Cu_Ka']) >>> print s.metadata.Sample.xray_lines ['Al_Ka', 'C_Ka', 'Cu_La', 'Mn_La', 'Zr_La'] ['Al_Ka', 'C_Ka', 'Cu_Ka', 'Cu_La', 'Mn_La', 'Zr_La']
See also

estimate_background_windows
(line_width=[2, 2], windows_width=1, xray_lines=None)¶ Estimate two windows around each Xray line containing only the background.
Parameters:  line_width (list of two floats) – The position of the two windows around the Xray line is given by the line_width (left and right) times the calculated FWHM of the line.
 windows_width (float) – The width of the windows is is the windows_width times the calculated FWHM of the line.
 xray_lines (None or list of string) – If None, use metadata.Sample.elements.xray_lines. Else, provide an iterable containing a list of valid Xray lines symbols.
Returns: windows_position – The position of the windows in energy. Each line corresponds to a Xray line. In a line, the two first values correspond to the limits of the left window and the two last values correspond to the limits of the right window.
Return type: 2D array of float
Examples
>>> s = hs.datasets.example_signals.EDS_TEM_Spectrum() >>> s.add_lines() >>> bw = s.estimate_background_windows(line_width=[5.0, 2.0]) >>> s.plot(background_windows=bw) >>> s.get_lines_intensity(background_windows=bw, plot_result=True) Fe_Ka at 6.4039 keV : Intensity = 2754.00 Pt_La at 9.4421 keV : Intensity = 15090.00
See also

estimate_integration_windows
(windows_width=2.0, xray_lines=None)¶ Estimate a window of integration for each Xray line.
Parameters:  windows_width (float) – The width of the integration windows is the ‘windows_width’ times the calculated FWHM of the line.
 xray_lines (None or list of string) – If None, use ‘metadata.Sample.elements.xray_lines’. Else, provide an iterable containing a list of valid Xray lines symbols.
Returns: integration_windows – The positions of the windows in energy. Each row corresponds to a Xray line. Each row contains the left and right value of the window.
Return type: 2D array of float
Examples
>>> s = hs.datasets.example_signals.EDS_TEM_Spectrum() >>> s.add_lines() >>> iw = s.estimate_integration_windows() >>> s.plot(integration_windows=iw) >>> s.get_lines_intensity(integration_windows=iw, plot_result=True) Fe_Ka at 6.4039 keV : Intensity = 3710.00 Pt_La at 9.4421 keV : Intensity = 15872.00
See also

get_lines_intensity
(xray_lines=None, integration_windows=2.0, background_windows=None, plot_result=False, only_one=True, only_lines=('a', ), **kwargs)¶ Return the intensity map of selected Xray lines.
The intensities, the number of Xray counts, are computed by suming the spectrum over the different Xray lines. The sum window width is calculated from the energy resolution of the detector as defined in ‘energy_resolution_MnKa’ of the metadata. Backgrounds average in provided windows can be subtracted from the intensities.
Parameters:  xray_lines ({None, "best", list of string}) – If None, if metadata.Sample.elements.xray_lines contains a list of lines use those. If metadata.Sample.elements.xray_lines is undefined or empty but metadata.Sample.elements is defined, use the same syntax as add_line to select a subset of lines for the operation. Alternatively, provide an iterable containing a list of valid Xray lines symbols.
 integration_windows (Float or array) – If float, the width of the integration windows is the ‘integration_windows_width’ times the calculated FWHM of the line. Else provide an array for which each row corresponds to a Xray line. Each row contains the left and right value of the window.
 background_windows (None or 2D array of float) – If None, no background subtraction. Else, the backgrounds average in the windows are subtracted from the return intensities. ‘background_windows’ provides the position of the windows in energy. Each line corresponds to a Xray line. In a line, the two first values correspond to the limits of the left window and the two last values correspond to the limits of the right window.
 plot_result (bool) – If True, plot the calculated line intensities. If the current object is a single spectrum it prints the result instead.
 only_one (bool) – If False, use all the lines of each element in the data spectral range. If True use only the line at the highest energy above an overvoltage of 2 (< beam energy / 2).
 only_lines ({None, list of strings}) – If not None, use only the given lines.
 kwargs – The extra keyword arguments for plotting. See utils.plot.plot_signals
Returns: intensities – A list containing the intensities as Signal subclasses.
Return type: list
Examples
>>> s = hs.datasets.example_signals.EDS_SEM_Spectrum() >>> s.get_lines_intensity(['Mn_Ka'], plot_result=True) Mn_La at 0.63316 keV : Intensity = 96700.00
>>> s = hs.datasets.example_signals.EDS_SEM_Spectrum() >>> s.plot(['Mn_Ka'], integration_windows=2.1) >>> s.get_lines_intensity(['Mn_Ka'], >>> integration_windows=2.1, plot_result=True) Mn_Ka at 5.8987 keV : Intensity = 53597.00
>>> s = hs.datasets.example_signals.EDS_SEM_Spectrum() >>> s.set_elements(['Mn']) >>> s.set_lines(['Mn_Ka']) >>> bw = s.estimate_background_windows() >>> s.plot(background_windows=bw) >>> s.get_lines_intensity(background_windows=bw, plot_result=True) Mn_Ka at 5.8987 keV : Intensity = 46716.00

get_take_off_angle
()¶ Calculate the takeoffangle (TOA).
TOA is the angle with which the Xrays leave the surface towards the detector. Parameters are read in ‘SEM.tilt_stage’, ‘Acquisition_instrument.SEM.Detector.EDS.azimuth_angle’ and ‘SEM.Detector.EDS.elevation_angle’ in ‘metadata’.
Returns: take_off_angle – in Degree Return type: float Examples
>>> s = hs.datasets.example_signals.EDS_SEM_Spectrum() >>> s.get_take_off_angle() 37.0 >>> s.set_microscope_parameters(tilt_stage=20.) >>> s.get_take_off_angle() 57.0
See also
hs.eds.take_off_angle()
Notes
Defined by M. Schaffer et al., Ultramicroscopy 107(8), pp 587597 (2007)

plot
(xray_lines=False, only_lines=('a', 'b'), only_one=False, background_windows=None, integration_windows=None, **kwargs)¶ Plot the EDS spectrum. The following markers can be added
 The position of the Xray lines and their names.
 The background windows associated with each Xray lines. A black line
links the left and right window with the average value in each window.
Parameters:  xray_lines ({False, True, 'from_elements', list of string}) – If not False, indicate the position and the name of the Xray lines. If True, if metadata.Sample.elements.xray_lines contains a list of lines use those. If metadata.Sample.elements.xray_lines is undefined or empty or if xray_lines equals ‘from_elements’ and metadata.Sample.elements is defined, use the same syntax as add_line to select a subset of lines for the operation. Alternatively, provide an iterable containing a list of valid Xray lines symbols.
 only_lines (None or list of strings) – If not None, use only the given lines (eg. (‘a’,’Kb’)). If None, use all lines.
 only_one (bool) – If False, use all the lines of each element in the data spectral range. If True use only the line at the highest energy above an overvoltage of 2 (< beam energy / 2).
 background_windows (None or 2D array of float) – If not None, add markers at the position of the windows in energy. Each line corresponds to a Xray lines. In a line, the two first value corresponds to the limit of the left window and the two last values corresponds to the limit of the right window.
 integration_windows (None or 'auto' or float or 2D array of float) – If not None, add markers at the position of the integration windows. If ‘auto’ (or float), the width of the integration windows is 2.0 (or float) times the calculated FWHM of the line. see ‘estimate_integration_windows’. Else provide an array for which each row corresponds to a Xray line. Each row contains the left and right value of the window.
 kwargs – The extra keyword arguments for plot()
Examples
>>> s = hs.datasets.example_signals.EDS_SEM_Spectrum() >>> s.plot()
>>> s = hs.datasets.example_signals.EDS_SEM_Spectrum() >>> s.plot(True)
>>> s = hs.datasets.example_signals.EDS_TEM_Spectrum() >>> s.add_lines() >>> bw = s.estimate_background_windows() >>> s.plot(background_windows=bw)
>>> s = hs.datasets.example_signals.EDS_SEM_Spectrum() >>> s.plot(['Mn_Ka'], integration_windows='auto')
>>> s = hs.datasets.example_signals.EDS_TEM_Spectrum() >>> s.add_lines() >>> bw = s.estimate_background_windows() >>> s.plot(background_windows=bw, integration_windows=2.1)

rebin
(new_shape)¶ Rebins the data to the new shape
Parameters: new_shape (tuple of ints) – The new shape must be a divisor of the original shape Examples
>>> s = hs.datasets.example_signals.EDS_SEM_Spectrum() >>> print s >>> print s.rebin([512]) <EDSSEMSpectrum, title: EDS SEM Spectrum, dimensions: (1024)> <EDSSEMSpectrum, title: EDS SEM Spectrum, dimensions: (512)>

set_elements
(elements)¶ Erase all elements and set them.
Parameters: elements (list of strings) – A list of chemical element symbols. See also
Examples
>>> s = hs.datasets.example_signals.EDS_SEM_Spectrum() >>> print s.metadata.Sample.elements >>> s.set_elements(['Al']) >>> print s.metadata.Sample.elements ['Al' 'C' 'Cu' 'Mn' 'Zr'] ['Al']

set_lines
(lines, only_one=True, only_lines=('a', ))¶ Erase all Xrays lines and set them.
See add_lines for details.
Parameters:  lines (list of strings) – A list of valid element Xray lines to add e.g. Fe_Kb. Additionally, if metadata.Sample.elements is defined, add the lines of those elements that where not given in this list.
 only_one (bool) – If False, add all the lines of each element in metadata.Sample.elements that has not line defined in lines. If True (default), only add the line at the highest energy above an overvoltage of 2 (< beam energy / 2).
 only_lines ({None, list of strings}) – If not None, only the given lines will be added.
Examples
>>> s = hs.datasets.example_signals.EDS_SEM_Spectrum() >>> s.add_lines() >>> print s.metadata.Sample.xray_lines >>> s.set_lines(['Cu_Ka']) >>> print s.metadata.Sample.xray_lines ['Al_Ka', 'C_Ka', 'Cu_La', 'Mn_La', 'Zr_La'] ['Al_Ka', 'C_Ka', 'Cu_Ka', 'Mn_La', 'Zr_La']
See also

sum
(axis)¶ Sum the data over the given axis.
Parameters: axis ({int, string}) – The axis can be specified using the index of the axis in axes_manager or the axis name. Returns: s Return type: Signal See also
sum_in_mask()
,mean()
Examples
>>> s = hs.datasets.example_signals.EDS_SEM_Spectrum() >>> s.sum(0).data array(1000279)

hyperspy._signals.eds_sem module¶

class
hyperspy._signals.eds_sem.
EDSSEMSpectrum
(*args, **kwards)¶ Bases:
hyperspy._signals.eds.EDSSpectrum

get_calibration_from
(ref, nb_pix=1)¶ Copy the calibration and all metadata of a reference.
Primary use: To add a calibration to ripple file from INCA software
Parameters:  ref (signal) – The reference contains the calibration in its metadata
 nb_pix (int) – The live time (real time corrected from the “dead time”) is divided by the number of pixel (spectrums), giving an average live time.
Examples
>>> ref = hs.datasets.example_signals.EDS_SEM_Spectrum() >>> s = hs.signals.EDSSEMSpectrum( >>> hs.datasets.example_signals.EDS_SEM_Spectrum().data) >>> print s.axes_manager[0].scale >>> s.get_calibration_from(ref) >>> print s.axes_manager[0].scale 1.0 0.01

set_microscope_parameters
(beam_energy=None, live_time=None, tilt_stage=None, azimuth_angle=None, elevation_angle=None, energy_resolution_MnKa=None)¶ Set the microscope parameters.
If no arguments are given, raises an interactive mode to fill the values.
Parameters:  beam_energy (float) – The energy of the electron beam in keV
 live_time (float) – In second
 tilt_stage (float) – In degree
 azimuth_angle (float) – In degree
 elevation_angle (float) – In degree
 energy_resolution_MnKa (float) – In eV
Examples
>>> s = hs.datasets.example_signals.EDS_SEM_Spectrum() >>> print('Default value %s eV' % >>> s.metadata.Acquisition_instrument. >>> SEM.Detector.EDS.energy_resolution_MnKa) >>> s.set_microscope_parameters(energy_resolution_MnKa=135.) >>> print('Now set to %s eV' % >>> s.metadata.Acquisition_instrument. >>> SEM.Detector.EDS.energy_resolution_MnKa) Default value 130.0 eV Now set to 135.0 eV

hyperspy._signals.eds_tem module¶

class
hyperspy._signals.eds_tem.
EDSTEMSpectrum
(*args, **kwards)¶ Bases:
hyperspy._signals.eds.EDSSpectrum

decomposition
(normalize_poissonian_noise=True, navigation_mask=1.0, closing=True, *args, **kwargs)¶ Decomposition with a choice of algorithms
The results are stored in self.learning_results
Parameters:  normalize_poissonian_noise (bool) – If True, scale the SI to normalize Poissonian noise
 navigation_mask (None or float or boolean numpy array) – The navigation locations marked as True are not used in the decompostion. If float is given the vacuum_mask method is used to generate a mask with the float value as threshold.
 closing (bool) – If true, applied a morphologic closing to the maks obtained by vacuum_mask.
 algorithm ('svd'  'fast_svd'  'mlpca'  'fast_mlpca'  'nmf' ) – ‘sparse_pca’  ‘mini_batch_sparse_pca’
 output_dimension (None or int) – number of components to keep/calculate
 centre (None  'variables'  'trials') – If None no centring is applied. If ‘variable’ the centring will be performed in the variable axis. If ‘trials’, the centring will be performed in the ‘trials’ axis. It only has effect when using the svd or fast_svd algorithms
 auto_transpose (bool) – If True, automatically transposes the data to boost performance. Only has effect when using the svd of fast_svd algorithms.
 signal_mask (boolean numpy array) – The signal locations marked as True are not used in the decomposition.
 var_array (numpy array) – Array of variance for the maximum likelihood PCA algorithm
 var_func (function or numpy array) – If function, it will apply it to the dataset to obtain the var_array. Alternatively, it can a an array with the coefficients of a polynomial.
 polyfit –
 reproject (None  signal  navigation  both) – If not None, the results of the decomposition will be projected in the selected masked area.
Examples
>>> s = hs.datasets.example_signals.EDS_TEM_Spectrum() >>> si = hs.stack([s]*3) >>> si.change_dtype(float) >>> si.decomposition()
See also

get_calibration_from
(ref, nb_pix=1)¶ Copy the calibration and all metadata of a reference.
Primary use: To add a calibration to ripple file from INCA software
Parameters:  ref (signal) – The reference contains the calibration in its metadata
 nb_pix (int) – The live time (real time corrected from the “dead time”) is divided by the number of pixel (spectrums), giving an average live time.
Examples
>>> ref = hs.datasets.example_signals.EDS_TEM_Spectrum() >>> s = hs.signals.EDSTEMSpectrum( >>> hs.datasets.example_signals.EDS_TEM_Spectrum().data) >>> print s.axes_manager[0].scale >>> s.get_calibration_from(ref) >>> print s.axes_manager[0].scale 1.0 0.020028

quantification
(intensities, kfactors, composition_units='weight', navigation_mask=1.0, closing=True, plot_result=False, **kwargs)¶ Quantification of intensities to return elemental composition
Method: CliffLorimer
Parameters:  intensities (list of signal) – the intensitiy for each Xray lines.
 kfactors (list of float) – The list of kfactor in same order as intensities. Note that intensities provided by hyperspy are sorted by the aplhabetical order of the Xray lines. eg. kfactors =[0.982, 1.32, 1.60] for [‘Al_Ka’,’Cr_Ka’, ‘Ni_Ka’].
 composition_units ('weight' or 'atomic') – Quantification returns weight percent. By choosing ‘atomic’, the return composition is in atomic percent.
 navigation_mask (None or float or signal) – The navigation locations marked as True are not used in the quantification. If int is given the vacuum_mask method is used to generate a mask with the int value as threhsold. Else provides a signal with the navigation shape.
 closing (bool) – If true, applied a morphologic closing to the mask obtained by vacuum_mask.
 plot_result (bool) – If True, plot the calculated composition. If the current object is a single spectrum it prints the result instead.
 kwargs – The extra keyword arguments are passed to plot.
Returns:  A list of quantified elemental maps (signal) giving the composition of
 the sample in weight or atomic percent.
Examples
>>> s = hs.datasets.example_signals.EDS_TEM_Spectrum() >>> s.add_lines() >>> kfactors = [1.450226, 5.075602] #For Fe Ka and Pt La >>> bw = s.estimate_background_windows(line_width=[5.0, 2.0]) >>> s.plot(background_windows=bw) >>> intensities = s.get_lines_intensity(background_windows=bw) >>> res = s.quantification(intensities, kfactors, plot_result=True, >>> composition_units='atomic') Fe (Fe_Ka): Composition = 15.41 atomic percent Pt (Pt_La): Composition = 84.59 atomic percent
See also

set_microscope_parameters
(beam_energy=None, live_time=None, tilt_stage=None, azimuth_angle=None, elevation_angle=None, energy_resolution_MnKa=None)¶ Set the microscope parameters.
If no arguments are given, raises an interactive mode to fill the values.
Parameters:  beam_energy (float) – The energy of the electron beam in keV
 live_time (float) – In second
 tilt_stage (float) – In degree
 azimuth_angle (float) – In degree
 elevation_angle (float) – In degree
 energy_resolution_MnKa (float) – In eV
Examples
>>> s = hs.datasets.example_signals.EDS_TEM_Spectrum() >>> print(s.metadata.Acquisition_instrument. >>> TEM.Detector.EDS.energy_resolution_MnKa) >>> s.set_microscope_parameters(energy_resolution_MnKa=135.) >>> print(s.metadata.Acquisition_instrument. >>> TEM.Detector.EDS.energy_resolution_MnKa) 133.312296 135.0

vacuum_mask
(threshold=1.0, closing=True, opening=False)¶ Generate mask of the vacuum region
Parameters:  threshold (float) – For a given pixel, maximum value in the energy axis below which the pixel is considered as vacuum.
 closing (bool) – If true, applied a morphologic closing to the mask
 opnening (bool) – If true, applied a morphologic opening to the mask
Examples
>>> # Simulate a spectrum image with vacuum region >>> s = hs.datasets.example_signals.EDS_TEM_Spectrum() >>> s_vac = hs.signals.Simulation(np.ones_like(s.data, dtype=float))*0.005 >>> s_vac.add_poissonian_noise() >>> si = hs.stack([s]*3 + [s_vac]) >>> si.vacuum_mask().data array([False, False, False, True], dtype=bool)
Returns: mask – The mask of the region Return type: signal

hyperspy._signals.eels module¶

class
hyperspy._signals.eels.
EELSSpectrum
(*args, **kwards)¶ Bases:
hyperspy.signals.Spectrum

add_elements
(elements, include_pre_edges=False)¶ Declare the elemental composition of the sample.
The ionisation edges of the elements present in the current energy range will be added automatically.
Parameters:  elements (tuple of strings) – The symbol of the elements. Note this input must always be in the form of a tuple. Meaning: add_elements((‘C’,)) will work, while add_elements((‘C’)) will NOT work.
 include_pre_edges (bool) – If True, the ionization edges with an onset below the lower energy limit of the SI will be incluided
Examples
>>> s = hs.signals.EELSSpectrum(np.arange(1024)) >>> s.add_elements(('C', 'O')) Adding C_K subshell Adding O_K subshell
Raises: ValueError

align_zero_loss_peak
(calibrate=True, also_align=[], print_stats=True, subpixel=True, mask=None, show_progressbar=None, **kwargs)¶ Align the zeroloss peak.
This function first aligns the spectra using the result of estimate_zero_loss_peak_centre and afterward, if subpixel is True, proceeds to align with subpixel accuracy using align1D. The offset is automatically correct if calibrate is True.
Parameters:  calibrate (bool) – If True, set the offset of the spectral axis so that the zeroloss peak is at position zero.
 also_align (list of signals) – A list containing other spectra of identical dimensions to align using the shifts applied to the current spectrum. If calibrate is True, the calibration is also applied to the spectra in the list.
 print_stats (bool) – If True, print summary statistics the ZLP maximum before the aligment.
 subpixel (bool) – If True, perform the alignment with subpixel accuracy using crosscorrelation.
 mask (Signal 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.
See also
estimate_zero_loss_peak_centre()
,align1D()
,estimate_shift1D.()
Notes
Any extra keyword arguments are passed to align1D. For more information read its docstring.

create_model
(ll=None, auto_background=True, auto_add_edges=True, GOS=None)¶ Create a model for the current EELS data.
Parameters:  ll (EELSSpectrum, optional) – If an EELSSpectrum is provided, it will be assumed that it is a lowloss EELS spectrum, and it will be used to simulate the effect of multiple scattering by convolving it with the EELS spectrum.
 auto_background (boolean, default True) – If True, and if spectrum is an EELS instance adds automatically a powerlaw to the model and estimate the parameters by the twoarea method.
 auto_add_edges (boolean, default True) – If True, and if spectrum is an EELS instance, it will automatically add the ionization edges as defined in the Spectrum instance. Adding a new element to the spectrum using the components.EELSSpectrum.add_elements method automatically add the corresponding ionisation edges to the model.
 GOS ({'hydrogenic'  'HartreeSlater'}, optional) – The generalized oscillation strenght calculations to use for the coreloss EELS edges. If None the HartreeSlater GOS are used if available, otherwise it uses the hydrogenic GOS.
Returns: model
Return type: EELSModel instance.

estimate_elastic_scattering_intensity
(threshold, show_progressbar=None)¶ Rough estimation of the elastic scattering intensity by truncation of a EELS lowloss spectrum.
Parameters:  threshold ({Signal, float, int}) – Truncation energy to estimate the intensity of the elastic scattering. The threshold can be provided as a signal of the same dimension as the input spectrum navigation space containing the threshold value in the energy units. Alternatively a constant threshold can be specified in energy/index units by passing float/int.
 show_progressbar (None or bool) – If True, display a progress bar. If None the default is set in preferences.
Returns: I0 – The elastic scattering intensity.
Return type:

estimate_elastic_scattering_threshold
(window=10.0, tol=None, window_length=5, polynomial_order=3, start=1.0)¶ Calculate the first inflexion point of the spectrum derivative within a window.
This method assumes that the zeroloss peak is located at position zero in all the spectra. Currently it looks for an inflexion point, that can be a local maximum or minimum. Therefore, to estimate the elastic scattering threshold start + window must be less than the first maximum for all spectra (often the bulk plasmon maximum). If there is more than one inflexion point in energy the window it selects the smoother one what, often, but not always, is a good choice in this case.
Parameters:  window ({None, float}) – If None, the search for the local inflexion point is performed using the full energy range. A positive float will restrict the search to the (0,window] energy window, where window is given in the axis units. If no inflexion point is found in this spectral range the window value is returned instead.
 tol ({None, float}) – The threshold tolerance for the derivative. If “auto” it is automatically calculated as the minimum value that guarantees finding an inflexion point in all the spectra in given energy range.
 window_length (int) – If non zero performs order three SavitzkyGolay smoothing to the data to avoid falling in local minima caused by the noise. It must be an odd interger.
 polynomial_order (int) – SavitzkyGolay filter polynomial order.
 start (float) – Position from the zeroloss peak centre from where to start looking for the inflexion point.
Returns: threshold – A Signal of the same dimension as the input spectrum navigation space containing the estimated threshold. Where the threshold couldn’t be estimated the value is set to nan.
Return type: See also
estimate_elastic_scattering_intensity()
,align_zero_loss_peak()
,find_peaks1D_ohaver()
,fourier_ratio_deconvolution.()
Notes
The main purpose of this method is to be used as input for estimate_elastic_scattering_intensity. Indeed, for currently achievable energy resolutions, there is not such a thing as a elastic scattering threshold. Therefore, please be aware of the limitations of this method when using it.

estimate_thickness
(threshold, zlp=None)¶ Estimates the thickness (relative to the mean free path) of a sample using the logratio method.
The current EELS spectrum must be a lowloss spectrum containing the zeroloss peak. The hyperspectrum must be well calibrated and aligned.
Parameters:  threshold ({Signal, float, int}) – Truncation energy to estimate the intensity of the elastic scattering. The threshold can be provided as a signal of the same dimension as the input spectrum navigation space containing the threshold value in the energy units. Alternatively a constant threshold can be specified in energy/index units by passing float/int.
 zlp ({None, EELSSpectrum}) – If not None the zeroloss peak intensity is calculated from the ZLP spectrum supplied by integration using Simpson’s rule. If None estimates the zeroloss peak intensity using estimate_elastic_scattering_intensity by truncation.
Returns: s – The thickness relative to the MFP. It returns a Spectrum, Image or a Signal, depending on the currenct spectrum navigation dimensions.
Return type: Notes
For details see: Egerton, R. Electron EnergyLoss Spectroscopy in the Electron Microscope. SpringerVerlag, 2011.

estimate_zero_loss_peak_centre
(mask=None)¶ Estimate the posision of the zeroloss peak.
This function provides just a coarse estimation of the position of the zeroloss peak centre by computing the position of the maximum of the spectra. For subpixel accuracy use estimate_shift1D.
Parameters: mask (Signal 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. Returns: zlpc – The estimated position of the maximum of the ZLP peak. Return type: Signal subclass Notes
This function only works when the zeroloss peak is the most intense feature in the spectrum. If it is not in most cases the spectrum can be cropped to meet this criterium. Alternatively use estimate_shift1D.
See also
estimate_shift1D()
,align_zero_loss_peak()

fourier_log_deconvolution
(zlp, add_zlp=False, crop=False)¶ Performs fourierlog deconvolution.
Parameters:  zlp (EELSSpectrum) – The corresponding zeroloss peak.
 add_zlp (bool) – If True, adds the ZLP to the deconvolved spectrum
 crop (bool) – If True crop the spectrum to leave out the channels that have been modified to decay smoothly to zero at the sides of the spectrum.
Returns: Return type: An EELSSpectrum containing the current data deconvolved.
Notes
For details see: Egerton, R. Electron EnergyLoss Spectroscopy in the Electron Microscope. SpringerVerlag, 2011.

fourier_ratio_deconvolution
(ll, fwhm=None, threshold=None, extrapolate_lowloss=True, extrapolate_coreloss=True)¶ Performs Fourierratio deconvolution.
 The coreloss should have the background removed. To reduce
 the noise amplication the result is convolved with a
Gaussian function.
Parameters:  ll (EELSSpectrum) – The corresponding lowloss (ll) EELSSpectrum.
 fwhm (float or None) – Fullwidth halfmaximum of the Gaussian function by which the result of the deconvolution is convolved. It can be used to select the final SNR and spectral resolution. If None, the FWHM of the zeroloss peak of the lowloss is estimated and used.
 threshold ({None, float}) –
Truncation energy to estimate the intensity of the elastic scattering. If None the threshold is taken as the
first minimum after the ZLP centre.  extrapolate_coreloss (extrapolate_lowloss,) – If True the signals are extrapolated using a power law,
Notes
For details see: Egerton, R. Electron EnergyLoss Spectroscopy in the Electron Microscope. SpringerVerlag, 2011.

generate_subshells
(include_pre_edges=False)¶ Calculate the subshells for the current energy range for the elements present in self.elements
Parameters: include_pre_edges (bool) – If True, the ionization edges with an onset below the lower energy limit of the SI will be incluided

kramers_kronig_analysis
(zlp=None, iterations=1, n=None, t=None, delta=0.5, full_output=False)¶ Calculate the complex dielectric function from a single scattering distribution (SSD) using the KramersKronig relations.
It uses the FFT method as in [Egerton2011]. The SSD is an EELSSpectrum instance containing SSD lowloss EELS with no zeroloss peak. The internal loop is devised to approximately subtract the surface plasmon contribution supposing an unoxidized planar surface and neglecting coupling between the surfaces. This method does not account for retardation effects, instrumental broading and surface plasmon excitation in particles.
Note that either refractive index or thickness are required. If both are None or if both are provided an exception is raised.
Parameters:  zlp ({None, number, Signal}) – ZLP intensity. It is optional (can be None) if t is None and n is not None and the thickness estimation is not required. If t is not None, the ZLP is required to perform the normalization and if t is not None, the ZLP is required to calculate the thickness. If the ZLP is the same for all spectra, the integral of the ZLP can be provided as a number. Otherwise, if the ZLP intensity is not the same for all spectra, it can be provided as i) a Signal of the same dimensions as the current signal containing the ZLP spectra for each location ii) a Signal of signal dimension 0 and navigation_dimension equal to the current signal containing the integrated ZLP intensity.
 iterations (int) – Number of the iterations for the internal loop to remove the surface plasmon contribution. If 1 the surface plasmon contribution is not estimated and subtracted (the default is 1).
 n ({None, float}) – The medium refractive index. Used for normalization of the SSD to obtain the energy loss function. If given the thickness is estimated and returned. It is only required when t is None.
 t ({None, number, Signal}) – The sample thickness in nm. Used for normalization of the SSD to obtain the energy loss function. It is only required when n is None. If the thickness is the same for all spectra it can be given by a number. Otherwise, it can be provided as a Signal with signal dimension 0 and navigation_dimension equal to the current signal.
 delta (float) – A small number (0.10.5 eV) added to the energy axis in specific steps of the calculation the surface loss correction to improve stability.
 full_output (bool) – If True, return a dictionary that contains the estimated thickness if t is None and the estimated surface plasmon excitation and the spectrum corrected from surface plasmon excitations if iterations > 1.
Returns: eps (DielectricFunction instance) –
 The complex dielectric function results,
$epsilon = epsilon_1 + i*epsilon_2$,
contained in an DielectricFunction instance.
output (Dictionary (optional)) – A dictionary of optional outputs with the following keys:
thickness
The estimated thickness in nm calculated by normalization of the SSD (only when t is None)
surface plasmon estimation
The estimated surface plasmon excitation (only if iterations > 1.)
Raises: ValuerError
– If both n and t are undefined (None).AttribureError
– If the beam_energy or the collection semiangle are not defined in metadata.
Notes
This method is based in Egerton’s Matlab code [Egerton2011] with some minor differences:
 The integrals are performed using the simpsom rule instead of using a summation.
 The wraparound problem when computing the ffts is workarounded by padding the signal instead of substracting the reflected tail.
[Egerton2011] Ray Egerton, “Electron EnergyLoss Spectroscopy in the Electron Microscope”, SpringerVerlag, 2011.

power_law_extrapolation
(window_size=20, extrapolation_size=1024, add_noise=False, fix_neg_r=False)¶ Extrapolate the spectrum to the right using a powerlaw
Parameters:  window_size (int) – The number of channels from the right side of the spectrum that are used to estimate the power law parameters.
 extrapolation_size (int) – Size of the extrapolation in number of channels
 add_noise (bool) – If True, add poissonian noise to the extrapolated spectrum.
 fix_neg_r (bool) – If True, the negative values for the “components.PowerLaw” parameter r will be flagged and the extrapolation will be done with a constant zerovalue.
Returns: Return type: A new spectrum, with the extrapolation.

richardson_lucy_deconvolution
(psf, iterations=15, mask=None, show_progressbar=None)¶ 1D RichardsonLucy Poissonian deconvolution of the spectrum by the given kernel.
Parameters:  iterations (int) – Number of iterations of the deconvolution. Note that increasing the value will increase the noise amplification.
 psf (EELSSpectrum) – It must have the same signal dimension as the current spectrum and a spatial dimension of 0 or the same as the current spectrum.
 show_progressbar (None or bool) – If True, display a progress bar. If None the default is set in preferences.
 Notes –
  –
 details on the algorithm see Gloter, A., A. Douiri, (For) –
 Tence, and C. Colliex. “Improving Energy Resolution of (M.) –
 Spectra (EELS) –
 96, no. 3–4 (September 2003) (Ultramicroscopy) –

set_microscope_parameters
(beam_energy=None, convergence_angle=None, collection_angle=None)¶ Set the microscope parameters that are necessary to calculate the GOS.
If not all of them are defined, raises in interactive mode raises an UI item to fill the values
 beam_energy: float
 The energy of the electron beam in keV
 convengence_angle : float
 The microscope convergence semiangle in mrad.
 collection_angle : float
 The collection semiangle in mrad.

hyperspy._signals.eels_spectrum_simulation module¶

class
hyperspy._signals.eels_spectrum_simulation.
EELSSpectrumSimulation
(*args, **kwargs)¶ Bases:
hyperspy._signals.spectrum_simulation.SpectrumSimulation
,hyperspy._signals.eels.EELSSpectrum
hyperspy._signals.image_simulation module¶

class
hyperspy._signals.image_simulation.
ImageSimulation
(*args, **kwargs)¶ Bases:
hyperspy._signals.simulation.Simulation
,hyperspy._signals.signal2d.Signal2D
hyperspy._signals.signal1d module¶

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

create_model
()¶ Create a model for the current data.
Returns: model Return type: Model instance.

spikes_removal_tool
(signal_mask=None, navigation_mask=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)
See also
_spikes_diagnosis()

to_signal2D
()¶ Returns the one dimensional signal as a two dimensional signal.
See also
as_signal2D()
 a method for the same purpose with more options.
signals.Signal1D.to_signal2D()
 performs the inverse operation on images.
Raises: DataDimensionError: when data.ndim < 2


class
hyperspy._signals.signal1d.
Signal1DTools
¶ Bases:
object

align1D
(start=None, end=None, reference_indices=None, max_shift=None, interpolate=True, number_of_interpolation_points=5, interpolation_method='linear', crop=True, fill_value=nan, also_align=[], 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.
 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) – A list of Signal 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 (Signal 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.
Raises: SignalDimensionError if the signal dimension is not 1.
See also

calibrate
(*args, **kwargs)¶ 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
Notes
For this method to work the output_dimension must be 1. Set the view accordingly
Raises: SignalDimensionError if the signal dimension is not 1.

crop_spectrum
(*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, 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.
 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, 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.
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
 mask (Signal 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:  An array with the result of the estimation in the axis units.
 Raises
 SignalDimensionError if the signal dimension is not 1.

filter_butterworth
(cutoff_frequency_ratio=None, type='low', order=2)¶ Butterworth filter in place.
Raises: 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)¶ 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
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
(*args, **kwargs)¶ 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
(*args, **kwargs)¶ 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')¶ 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.
Parameters: signal_range ({a tuple of this form (l, r), "interactive"}) – 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: Signal subclass See also
integrate_simpson()
Examples
Using the GUI
>>> s.integrate_in_range()
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, show_progressbar=None, **kwargs)¶ Replace the data in a given range by interpolation.
The operation is performed in place.
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.
 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.
 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='PowerLaw', polynomial_order=2, show_progressbar=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.
 show_progressbar (None or bool) – If True, display a progress bar. If None the default is set in preferences.
Examples
>>>> s.remove_background() # Using gui, replaces spectrum s >>>> s2 = s.remove_background(
signal_range=(400,450), background_type=’PowerLaw’) #Using cli, returns a spectrumRaises: SignalDimensionError if the signal dimension is not 1.

shift1D
(shift_array, interpolation_method='linear', crop=True, fill_value=nan, 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.
 fill_value (float) – If crop is False fill the data outside of the original interval with the given value where needed.
 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)¶ 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.
Raises:  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)¶ 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.
Notes
More information about the filter in scipy.signal.savgol_filter.

smooth_tv
(smoothing_parameter=None, show_progressbar=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.
Raises: SignalDimensionError if the signal dimension is not 1.


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.
Parameters:  y (array) – 1D input array, e.g. a spectrum
 x (array (optional)) – 1D array describing the calibration of y (must have same shape as y)
 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 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 = arange(0,50,0.01) >>> y = cos(x) >>> one_dim_findpeaks(y, x,0,0) array([[ 1.68144859e05, 9.99999943e01, 3.57487961e+00], [ 6.28318614e+00, 1.00000003e+00, 3.57589018e+00], [ 1.25663708e+01, 1.00000002e+00, 3.57600673e+00], [ 1.88495565e+01, 1.00000002e+00, 3.57597295e+00], [ 2.51327421e+01, 1.00000003e+00, 3.57590284e+00], [ 3.14159267e+01, 1.00000002e+00, 3.57600856e+00], [ 3.76991124e+01, 1.00000002e+00, 3.57597984e+00], [ 4.39822980e+01, 1.00000002e+00, 3.57591479e+00]])
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)¶
hyperspy._signals.signal2d module¶

class
hyperspy._signals.signal2d.
Signal2D
(*args, **kw)¶ Bases:
hyperspy.signal.BaseSignal
,hyperspy._signals.signal2d.Signal2DTools

plot
(colorbar=True, scalebar=True, scalebar_color='white', axes_ticks=None, saturated_pixels=0, vmin=None, vmax=None, no_nans=False, centre_colormap='auto', **kwargs)¶ Plot image.
For multidimensional datasets an optional figure, the “navigator”, with a cursor to navigate that data is raised. In any case it is possible to navigate the data using the sliders. Currently only signals with signal_dimension equal to 0, 1 and 2 can be plotted.
Parameters:  navigator ({"auto", None, "slider", "spectrum", Signal}) – If “auto”, if navigation_dimension > 0, a navigator is provided to explore the data. If navigation_dimension is 1 and the signal is an image the navigator is a spectrum obtained by integrating over the signal axes (the image). If navigation_dimension is 1 and the signal is a spectrum the navigator is an image obtained by stacking horizontally all the spectra in the dataset. If navigation_dimension is > 1, the navigator is an image obtained by integrating the data over the signal axes. Additionaly, if navigation_dimension > 2 a window with one slider per axis is raised to navigate the data. For example, if the dataset consists of 3 navigation axes X, Y, Z and one signal axis, E, the default navigator will be an image obtained by integrating the data over E at the current Z index and a window with sliders for the X, Y and Z axes will be raised. Notice that changing the Zaxis index changes the navigator in this case. If “slider” and the navigation dimension > 0 a window with one slider per axis is raised to navigate the data. If “spectrum” and navigation_dimension > 0 the navigator is always a spectrum obtained by integrating the data over all other axes. If None, no navigator will be provided. Alternatively a Signal instance can be provided. The signal dimension must be 1 (for a spectrum navigator) or 2 (for a image navigator) and navigation_shape must be 0 (for a static navigator) or navigation_shape + signal_shape must be equal to the navigator_shape of the current object (for a dynamic navigator). If the signal dtype is RGB or RGBA this parameters has no effect and is always “slider”.
 axes_manager ({None, axes_manager}) – If None axes_manager is used.
 colorbar (bool, optional) – If true, a colorbar is plotted for nonRGB images.
 scalebar (bool, optional) – If True and the units and scale of the x and y axes are the same a scale bar is plotted.
 scalebar_color (str, optional) – A valid MPL color string; will be used as the scalebar color.
 axes_ticks ({None, bool}, optional) – If True, plot the axes ticks. If None axes_ticks are only plotted when the scale bar is not plotted. If False the axes ticks are never plotted.
 saturated_pixels (scalar) – The percentage of pixels that are left out of the bounds. For example, the low and high bounds of a value of 1 are the 0.5% and 99.5% percentiles. It must be in the [0, 100] range.
 vmax (vmin,) – vmin and vmax are used to normalize luminance data.
 no_nans (bool, optional) – If True, set nans to zero for plotting.
 centre_colormap ({"auto", True, False}) – If True the centre of the color scheme is set to zero. This is specially useful when using diverging color schemes. If “auto” (default), diverging color schemes are automatically centred.
 optional (**kwargs,) – Additional key word arguments passed to matplotlib.imshow()

to_signal1D
()¶ Returns the image as a spectrum.
See also
as_spectrum()
 a method for the same purpose with more options.
signals.Signal1D.to_signal1D()
 performs the inverse operation on one
dimensional()


class
hyperspy._signals.signal2d.
Signal2DTools
¶ Bases:
object

align2D
(crop=True, fill_value=nan, shifts=None, roi=None, sobel=True, medfilter=True, hanning=True, plot=False, normalize_corr=False, reference='current', dtype='float', correlation_threshold=None, chunk_size=30, interpolation_order=1)¶ Align the images in place using user provided shifts or by estimating the shifts.
Please, see estimate_shift2D docstring for details on the rest of the parameters not documented in the following section
Parameters:  crop (bool) – If True, the data will be cropped not to include regions with missing data
 fill_value (int, float, nan) – The areas with missing data are filled with the given value. Default is nan.
 shifts (None or list of tuples) – If None the shifts are estimated using estimate_shift2D.
 interpolation_order (int, default 1.) – The order of the spline interpolation. Default is 1, linear interpolation.
Returns: shifts – The shifts are returned only if shifts is None
Return type: np.array
Notes
The statistical analysis approach to the translation estimation when using `reference`=’stat’ roughly follows [1]_ . If you use it please cite their article.
References
[1] Schaffer, Bernhard, Werner Grogger, and Gerald Kothleitner. “Automated Spatial Drift Correction for EFTEM Image Series.” Ultramicroscopy 102, no. 1 (December 2004): 27–36.

crop_image
(top=None, bottom=None, left=None, right=None)¶ Crops an image in place.
top, bottom, left, right : int or float
If int the values are taken as indices. If float the values are converted to indices.crop

estimate_shift2D
(reference='current', correlation_threshold=None, chunk_size=30, roi=None, normalize_corr=False, sobel=True, medfilter=True, hanning=True, plot=False, dtype='float', show_progressbar=None)¶ Estimate the shifts in a image using phase correlation
This method can only estimate the shift by comparing bidimensional features that should not change position between frames. To decrease the memory usage, the time of computation and the accuracy of the results it is convenient to select a region of interest by setting the roi keyword.
Parameters:  reference ({'current', 'cascade' ,'stat'}) – If ‘current’ (default) the image at the current coordinates is taken as reference. If ‘cascade’ each image is aligned with the previous one. If ‘stat’ the translation of every image with all the rest is estimated and by performing statistical analysis on the result the translation is estimated.
 correlation_threshold ({None, 'auto', float}) – This parameter is only relevant when reference is ‘stat’. If float, the shift estimations with a maximum correlation value lower than the given value are not used to compute the estimated shifts. If ‘auto’ the threshold is calculated automatically as the minimum maximum correlation value of the automatically selected reference image.
 chunk_size ({None, int}) – If int and `reference`==’stat’ the number of images used as reference are limited to the given value.
 roi (tuple of ints or floats (left, right, top bottom)) – Define the region of interest. If int(float) the position is given axis index(value).
 sobel (bool) – apply a sobel filter for edge enhancement
 medfilter (bool) – apply a median filter for noise reduction
 hanning (bool) – Apply a 2d hanning filter
 plot (bool) – If True plots the images after applying the filters and the phase correlation
 dtype (str or dtype) – Typecode or datatype in which the calculations must be performed.
 show_progressbar (None or bool) – If True, display a progress bar. If None the default is set in preferences.
Returns: Return type: list of applied shifts
Notes
The statistical analysis approach to the translation estimation when using `reference`=’stat’ roughly follows [1]_ . If you use it please cite their article.
References
[1] Schaffer, Bernhard, Werner Grogger, and Gerald Kothleitner. “Automated Spatial Drift Correction for EFTEM Image Series.” Ultramicroscopy 102, no. 1 (December 2004): 27–36.


hyperspy._signals.signal2d.
estimate_image_shift
(ref, image, roi=None, sobel=True, medfilter=True, hanning=True, plot=False, dtype='float', normalize_corr=False)¶ Estimate the shift in a image using phase correlation
This method can only estimate the shift by comparing bidimensional features that should not change the position in the given axis. To decrease the memory usage, the time of computation and the accuracy of the results it is convenient to select a region of interest by setting the roi keyword.
Parameters:  roi (tuple of ints (top, bottom, left, right)) – Define the region of interest
 sobel (bool) – apply a sobel filter for edge enhancement
 medfilter (bool) – apply a median filter for noise reduction
 hanning (bool) – Apply a 2d hanning filter
 plot (bool) – If True plots the images after applying the filters and the phase correlation
 reference ('current'  'cascade') – If ‘current’ (default) the image at the current coordinates is taken as reference. If ‘cascade’ each image is aligned with the previous one.
 dtype (str or dtype) – Typecode or datatype in which the calculations must be performed.
 normalize_corr (bool) – If True use phase correlation instead of standard correlation
Returns:  shifts (np.array) – containing the estimate shifts
 max_value (float) – The maximum value of the correlation

hyperspy._signals.signal2d.
fft_correlation
(in1, in2, normalize=False)¶ Correlation of two Ndimensional arrays using FFT.
Adapted from scipy’s fftconvolve.
Parameters:  in2 (in1,) –
 normalize (bool) – If True performs phase correlation

hyperspy._signals.signal2d.
hanning2d
(M, N)¶ A 2D hanning window created by outer product.

hyperspy._signals.signal2d.
shift_image
(im, shift, interpolation_order=1, fill_value=nan)¶

hyperspy._signals.signal2d.
sobel_filter
(im)¶

hyperspy._signals.signal2d.
triu_indices_minus_diag
(n)¶ Returns the indices for the uppertriangle of an (n, n) array excluding its diagonal
Parameters: n (int) – The length of the square array
hyperspy._signals.simulation module¶

class
hyperspy._signals.simulation.
Simulation
(*args, **kwargs)¶ Bases:
hyperspy.signal.BaseSignal

add_gaussian_noise
(*args, **kwargs)¶ Add Gaussian noise to the data :param std: :type std: float

add_poissonian_noise
(*args, **kwargs)¶ Add Poissonian noise to the data

hyperspy._signals.spectrum_simulation module¶

class
hyperspy._signals.spectrum_simulation.
SpectrumSimulation
(*args, **kwargs)¶ Bases:
hyperspy._signals.simulation.Simulation
,hyperspy._signals.signal1d.Signal1D