nmr-processing 9.5.0 → 9.5.1-pre.1684857934
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- package/lib/apodization/apodization.d.ts +9 -4
- package/lib/apodization/apodization.js +3 -1
- package/lib/apodization/apodization.js.map +1 -1
- package/lib/apodization/applyWindow.js.map +1 -1
- package/lib/apodization/compose.js.map +1 -1
- package/lib/apodization/getFunction.js +3 -0
- package/lib/apodization/getFunction.js.map +1 -1
- package/lib/apodization/shapes/WindowFunctions.d.ts +6 -1
- package/lib/apodization/shapes/sineBell.d.ts +27 -0
- package/lib/apodization/shapes/sineBell.js +17 -0
- package/lib/apodization/shapes/sineBell.js.map +1 -0
- package/lib/apodization/utils/getData.js.map +1 -1
- package/lib/filters/Data1D.d.ts +6 -0
- package/lib/{xyz/Data2D.js → filters/Data1D.js} +1 -1
- package/lib/filters/Data1D.js.map +1 -0
- package/lib/filters/Entry1D.d.ts +10 -0
- package/lib/filters/Entry1D.js +3 -0
- package/lib/filters/Entry1D.js.map +1 -0
- package/lib/filters/ExclusionZone.d.ts +5 -0
- package/lib/filters/ExclusionZone.js +3 -0
- package/lib/filters/ExclusionZone.js.map +1 -0
- package/lib/filters/Filters.d.ts +7 -0
- package/lib/filters/Filters.js +47 -0
- package/lib/filters/Filters.js.map +1 -0
- package/lib/filters/FiltersManager.d.ts +23 -0
- package/lib/filters/FiltersManager.js +130 -0
- package/lib/filters/FiltersManager.js.map +1 -0
- package/lib/filters/MatrixFilter.d.ts +39 -0
- package/lib/filters/MatrixFilter.js +3 -0
- package/lib/filters/MatrixFilter.js.map +1 -0
- package/lib/filters/MatrixOptions.d.ts +28 -0
- package/lib/filters/MatrixOptions.js +3 -0
- package/lib/filters/MatrixOptions.js.map +1 -0
- package/lib/filters/filter1d/apodization.d.ts +26 -0
- package/lib/filters/filter1d/apodization.js +77 -0
- package/lib/filters/filter1d/apodization.js.map +1 -0
- package/lib/filters/filter1d/baselineCorrection.d.ts +32 -0
- package/lib/filters/filter1d/baselineCorrection.js +61 -0
- package/lib/filters/filter1d/baselineCorrection.js.map +1 -0
- package/lib/filters/filter1d/digitalFilter.d.ts +15 -0
- package/lib/filters/filter1d/digitalFilter.js +47 -0
- package/lib/filters/filter1d/digitalFilter.js.map +1 -0
- package/lib/filters/filter1d/equallySpaced.d.ts +29 -0
- package/lib/filters/filter1d/equallySpaced.js +55 -0
- package/lib/filters/filter1d/equallySpaced.js.map +1 -0
- package/lib/filters/filter1d/exclusionZones.d.ts +11 -0
- package/lib/filters/filter1d/exclusionZones.js +35 -0
- package/lib/filters/filter1d/exclusionZones.js.map +1 -0
- package/lib/filters/filter1d/fft.d.ts +14 -0
- package/lib/filters/filter1d/fft.js +94 -0
- package/lib/filters/filter1d/fft.js.map +1 -0
- package/lib/filters/filter1d/phaseCorrection.d.ts +26 -0
- package/lib/filters/filter1d/phaseCorrection.js +81 -0
- package/lib/filters/filter1d/phaseCorrection.js.map +1 -0
- package/lib/filters/filter1d/shiftX.d.ts +16 -0
- package/lib/filters/filter1d/shiftX.js +30 -0
- package/lib/filters/filter1d/shiftX.js.map +1 -0
- package/lib/filters/filter1d/signalProcessing.d.ts +17 -0
- package/lib/filters/filter1d/signalProcessing.js +88 -0
- package/lib/filters/filter1d/signalProcessing.js.map +1 -0
- package/lib/filters/filter1d/utils/padDataToNextPowerOfTwo.d.ts +2 -0
- package/lib/filters/filter1d/utils/padDataToNextPowerOfTwo.js +41 -0
- package/lib/filters/filter1d/utils/padDataToNextPowerOfTwo.js.map +1 -0
- package/lib/filters/filter1d/zeroFilling.d.ts +18 -0
- package/lib/filters/filter1d/zeroFilling.js +68 -0
- package/lib/filters/filter1d/zeroFilling.js.map +1 -0
- package/lib/filters/filter2d/shiftX.d.ts +14 -0
- package/lib/filters/filter2d/shiftX.js +30 -0
- package/lib/filters/filter2d/shiftX.js.map +1 -0
- package/lib/filters/filter2d/shiftY.d.ts +14 -0
- package/lib/filters/filter2d/shiftY.js +30 -0
- package/lib/filters/filter2d/shiftY.js.map +1 -0
- package/lib/index.d.ts +2 -0
- package/lib/index.js +2 -0
- package/lib/index.js.map +1 -1
- package/lib/signals/signalsToFID.d.ts +3 -0
- package/lib/utilities/checkData2DFid.d.ts +3 -0
- package/lib/utilities/checkData2DFid.js +10 -0
- package/lib/utilities/checkData2DFid.js.map +1 -0
- package/lib/utilities/cloneData1D.d.ts +6 -0
- package/lib/utilities/cloneData1D.js +9 -0
- package/lib/utilities/cloneData1D.js.map +1 -0
- package/lib/utilities/cloneData2D.d.ts +2 -0
- package/lib/utilities/cloneData2D.js +28 -0
- package/lib/utilities/cloneData2D.js.map +1 -0
- package/lib/xyz/Entry2D.d.ts +10 -0
- package/lib/xyz/Entry2D.js +3 -0
- package/lib/xyz/Entry2D.js.map +1 -0
- package/lib/xyz/quadrature.d.ts +7 -0
- package/lib/xyz/quadrature.js +54 -0
- package/lib/xyz/quadrature.js.map +1 -0
- package/lib/xyz/util/fft2d/digitalFilter.d.ts +7 -0
- package/lib/xyz/util/fft2d/digitalFilter.js +20 -0
- package/lib/xyz/util/fft2d/digitalFilter.js.map +1 -0
- package/lib/xyz/util/fft2d/fftDirectDimension.d.ts +24 -0
- package/lib/xyz/util/fft2d/fftDirectDimension.js +57 -0
- package/lib/xyz/util/fft2d/fftDirectDimension.js.map +1 -0
- package/lib/xyz/util/fft2d/fftIndirectDimension.d.ts +23 -0
- package/lib/xyz/util/fft2d/fftIndirectDimension.js +72 -0
- package/lib/xyz/util/fft2d/fftIndirectDimension.js.map +1 -0
- package/lib/xyz/util/fft2d/removeDCOffset.d.ts +5 -0
- package/lib/xyz/util/fft2d/removeDCOffset.js +19 -0
- package/lib/xyz/util/fft2d/removeDCOffset.js.map +1 -0
- package/lib/xyz/util/fft2d/zeroFilling.d.ts +10 -0
- package/lib/xyz/util/fft2d/zeroFilling.js +29 -0
- package/lib/xyz/util/fft2d/zeroFilling.js.map +1 -0
- package/lib/xyz/util/padData.d.ts +5 -3
- package/lib/xyz/util/padData.js +7 -7
- package/lib/xyz/util/padData.js.map +1 -1
- package/lib/xyz/xyzAutoSignalsPicking.d.ts +2 -2
- package/lib/xyz/xyzAutoSignalsPicking.js.map +1 -1
- package/lib/xyz/xyzAutoZonesPicking.d.ts +2 -8
- package/lib/xyz/xyzAutoZonesPicking.js.map +1 -1
- package/lib/xyz/xyzBidimensionalFFT.d.ts +40 -0
- package/lib/xyz/xyzBidimensionalFFT.js +79 -0
- package/lib/xyz/xyzBidimensionalFFT.js.map +1 -0
- package/lib-esm/apodization/apodization.js +3 -1
- package/lib-esm/apodization/apodization.js.map +1 -1
- package/lib-esm/apodization/applyWindow.js.map +1 -1
- package/lib-esm/apodization/compose.js.map +1 -1
- package/lib-esm/apodization/getFunction.js +3 -0
- package/lib-esm/apodization/getFunction.js.map +1 -1
- package/lib-esm/apodization/shapes/sineBell.js +13 -0
- package/lib-esm/apodization/shapes/sineBell.js.map +1 -0
- package/lib-esm/apodization/utils/getData.js.map +1 -1
- package/lib-esm/filters/Data1D.js +2 -0
- package/lib-esm/filters/Data1D.js.map +1 -0
- package/lib-esm/filters/Entry1D.js +2 -0
- package/lib-esm/filters/Entry1D.js.map +1 -0
- package/lib-esm/filters/ExclusionZone.js +2 -0
- package/lib-esm/filters/ExclusionZone.js.map +1 -0
- package/lib-esm/filters/Filters.js +21 -0
- package/lib-esm/filters/Filters.js.map +1 -0
- package/lib-esm/filters/FiltersManager.js +124 -0
- package/lib-esm/filters/FiltersManager.js.map +1 -0
- package/lib-esm/filters/MatrixFilter.js +2 -0
- package/lib-esm/filters/MatrixFilter.js.map +1 -0
- package/lib-esm/filters/MatrixOptions.js +2 -0
- package/lib-esm/filters/MatrixOptions.js.map +1 -0
- package/lib-esm/filters/filter1d/apodization.js +70 -0
- package/lib-esm/filters/filter1d/apodization.js.map +1 -0
- package/lib-esm/filters/filter1d/baselineCorrection.js +52 -0
- package/lib-esm/filters/filter1d/baselineCorrection.js.map +1 -0
- package/lib-esm/filters/filter1d/digitalFilter.js +41 -0
- package/lib-esm/filters/filter1d/digitalFilter.js.map +1 -0
- package/lib-esm/filters/filter1d/equallySpaced.js +49 -0
- package/lib-esm/filters/filter1d/equallySpaced.js.map +1 -0
- package/lib-esm/filters/filter1d/exclusionZones.js +29 -0
- package/lib-esm/filters/filter1d/exclusionZones.js.map +1 -0
- package/lib-esm/filters/filter1d/fft.js +88 -0
- package/lib-esm/filters/filter1d/fft.js.map +1 -0
- package/lib-esm/filters/filter1d/phaseCorrection.js +75 -0
- package/lib-esm/filters/filter1d/phaseCorrection.js.map +1 -0
- package/lib-esm/filters/filter1d/shiftX.js +24 -0
- package/lib-esm/filters/filter1d/shiftX.js.map +1 -0
- package/lib-esm/filters/filter1d/signalProcessing.js +58 -0
- package/lib-esm/filters/filter1d/signalProcessing.js.map +1 -0
- package/lib-esm/filters/filter1d/utils/padDataToNextPowerOfTwo.js +37 -0
- package/lib-esm/filters/filter1d/utils/padDataToNextPowerOfTwo.js.map +1 -0
- package/lib-esm/filters/filter1d/zeroFilling.js +62 -0
- package/lib-esm/filters/filter1d/zeroFilling.js.map +1 -0
- package/lib-esm/filters/filter2d/shiftX.js +24 -0
- package/lib-esm/filters/filter2d/shiftX.js.map +1 -0
- package/lib-esm/filters/filter2d/shiftY.js +24 -0
- package/lib-esm/filters/filter2d/shiftY.js.map +1 -0
- package/lib-esm/index.js +2 -0
- package/lib-esm/index.js.map +1 -1
- package/lib-esm/utilities/checkData2DFid.js +6 -0
- package/lib-esm/utilities/checkData2DFid.js.map +1 -0
- package/lib-esm/utilities/cloneData1D.js +5 -0
- package/lib-esm/utilities/cloneData1D.js.map +1 -0
- package/lib-esm/utilities/cloneData2D.js +24 -0
- package/lib-esm/utilities/cloneData2D.js.map +1 -0
- package/lib-esm/xyz/Entry2D.js +2 -0
- package/lib-esm/xyz/Entry2D.js.map +1 -0
- package/lib-esm/xyz/quadrature.js +50 -0
- package/lib-esm/xyz/quadrature.js.map +1 -0
- package/lib-esm/xyz/util/fft2d/digitalFilter.js +16 -0
- package/lib-esm/xyz/util/fft2d/digitalFilter.js.map +1 -0
- package/lib-esm/xyz/util/fft2d/fftDirectDimension.js +50 -0
- package/lib-esm/xyz/util/fft2d/fftDirectDimension.js.map +1 -0
- package/lib-esm/xyz/util/fft2d/fftIndirectDimension.js +65 -0
- package/lib-esm/xyz/util/fft2d/fftIndirectDimension.js.map +1 -0
- package/lib-esm/xyz/util/fft2d/removeDCOffset.js +15 -0
- package/lib-esm/xyz/util/fft2d/removeDCOffset.js.map +1 -0
- package/lib-esm/xyz/util/fft2d/zeroFilling.js +25 -0
- package/lib-esm/xyz/util/fft2d/zeroFilling.js.map +1 -0
- package/lib-esm/xyz/util/padData.js +7 -4
- package/lib-esm/xyz/util/padData.js.map +1 -1
- package/lib-esm/xyz/xyzAutoSignalsPicking.js.map +1 -1
- package/lib-esm/xyz/xyzAutoZonesPicking.js.map +1 -1
- package/lib-esm/xyz/xyzBidimensionalFFT.js +75 -0
- package/lib-esm/xyz/xyzBidimensionalFFT.js.map +1 -0
- package/package.json +7 -2
- package/src/apodization/apodization.ts +6 -4
- package/src/apodization/applyWindow.ts +0 -1
- package/src/apodization/compose.ts +1 -0
- package/src/apodization/getFunction.ts +3 -0
- package/src/apodization/shapes/WindowFunctions.ts +7 -1
- package/src/apodization/shapes/sineBell.ts +41 -0
- package/src/apodization/utils/getData.ts +0 -1
- package/src/filters/Data1D.ts +7 -0
- package/src/filters/Entry1D.ts +11 -0
- package/src/filters/ExclusionZone.ts +5 -0
- package/src/filters/Filters.ts +24 -0
- package/src/filters/FiltersManager.ts +176 -0
- package/src/filters/MatrixFilter.ts +58 -0
- package/src/filters/MatrixOptions.ts +26 -0
- package/src/filters/filter1d/apodization.ts +90 -0
- package/src/filters/filter1d/baselineCorrection.ts +83 -0
- package/src/filters/filter1d/digitalFilter.ts +51 -0
- package/src/filters/filter1d/equallySpaced.ts +61 -0
- package/src/filters/filter1d/exclusionZones.ts +34 -0
- package/src/filters/filter1d/fft.ts +135 -0
- package/src/filters/filter1d/phaseCorrection.ts +118 -0
- package/src/filters/filter1d/shiftX.ts +31 -0
- package/src/filters/filter1d/signalProcessing.ts +79 -0
- package/src/filters/filter1d/utils/padDataToNextPowerOfTwo.ts +48 -0
- package/src/filters/filter1d/zeroFilling.ts +90 -0
- package/src/filters/filter2d/shiftX.ts +31 -0
- package/src/filters/filter2d/shiftY.ts +29 -0
- package/src/index.ts +3 -0
- package/src/utilities/checkData2DFid.ts +10 -0
- package/src/utilities/cloneData1D.ts +6 -0
- package/src/utilities/cloneData2D.ts +27 -0
- package/src/xyz/Entry2D.ts +12 -0
- package/src/xyz/quadrature.ts +68 -0
- package/src/xyz/util/fft2d/digitalFilter.ts +24 -0
- package/src/xyz/util/fft2d/fftDirectDimension.ts +111 -0
- package/src/xyz/util/fft2d/fftIndirectDimension.ts +123 -0
- package/src/xyz/util/fft2d/removeDCOffset.ts +19 -0
- package/src/xyz/util/fft2d/zeroFilling.ts +39 -0
- package/src/xyz/util/padData.ts +13 -7
- package/src/xyz/xyzAutoSignalsPicking.ts +2 -2
- package/src/xyz/xyzAutoZonesPicking.ts +3 -9
- package/src/xyz/xyzBidimensionalFFT.ts +122 -0
- package/lib/xyz/Data2D.d.ts +0 -7
- package/lib/xyz/Data2D.js.map +0 -1
- package/lib-esm/xyz/Data2D.js +0 -2
- package/lib-esm/xyz/Data2D.js.map +0 -1
- package/src/xyz/Data2D.ts +0 -7
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"shiftX.js","sourceRoot":"","sources":["../../../src/filters/filter1d/shiftX.ts"],"names":[],"mappings":"AAGA;;;;GAIG;AAEH,MAAM,CAAC,MAAM,EAAE,GAAG,QAAQ,CAAC;AAC3B,MAAM,CAAC,MAAM,IAAI,GAAG,SAAS,CAAC;AAE9B,MAAM,CAAC,MAAM,mBAAmB,GAAsC;IACpE,aAAa,EAAE,IAAI;IACnB,aAAa,EAAE,KAAK;CACrB,CAAC;AAEF,MAAM,UAAU,KAAK,CAAC,OAAgB,EAAE,UAAU,GAAG,CAAC;IACpD,OAAO,CAAC,IAAI,CAAC,CAAC,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,GAAG,UAAU,CAAC,CAAC;AACjE,CAAC;AAED,MAAM,UAAU,YAAY;IAC1B,OAAO,IAAI,CAAC;AACd,CAAC;AAED,MAAM,UAAU,MAAM,CAAC,aAAqB,EAAE,QAAgB;IAC5D,OAAO;QACL,IAAI,EAAE,IAAI;QACV,MAAM,EAAE,aAAa,GAAG,QAAQ;KACjC,CAAC;AACJ,CAAC"}
|
|
@@ -0,0 +1,58 @@
|
|
|
1
|
+
import * as Filters from 'ml-signal-processing';
|
|
2
|
+
export const id = 'signalProcessing';
|
|
3
|
+
export const name = 'Signal processing';
|
|
4
|
+
export const DOMAIN_UPDATE_RULES = {
|
|
5
|
+
updateXDomain: true,
|
|
6
|
+
updateYDomain: true,
|
|
7
|
+
};
|
|
8
|
+
export function apply(spectrum, options) {
|
|
9
|
+
if (!isApplicable(spectrum)) {
|
|
10
|
+
throw new Error('Signal Processing is not applicable on this data');
|
|
11
|
+
}
|
|
12
|
+
spectrum.data = filterXY(spectrum, options);
|
|
13
|
+
}
|
|
14
|
+
export function isApplicable(spectrum) {
|
|
15
|
+
if (spectrum.info.isFt)
|
|
16
|
+
return true;
|
|
17
|
+
return false;
|
|
18
|
+
}
|
|
19
|
+
export function reduce(previousValue, newValue) {
|
|
20
|
+
return {
|
|
21
|
+
once: true,
|
|
22
|
+
reduce: newValue,
|
|
23
|
+
};
|
|
24
|
+
}
|
|
25
|
+
export function filterXY(spectrum, options) {
|
|
26
|
+
var _a;
|
|
27
|
+
let x = spectrum.data.x.slice(0);
|
|
28
|
+
let cloneX = spectrum.data.x.slice(0);
|
|
29
|
+
let re = spectrum.data.re.slice(0);
|
|
30
|
+
let im = (_a = spectrum.data.im) === null || _a === void 0 ? void 0 : _a.slice(0);
|
|
31
|
+
const filters = options.filters.slice();
|
|
32
|
+
const { range: { from, to }, numberOfPoints, exclusionsZones: exclusions, } = options;
|
|
33
|
+
filters.push({
|
|
34
|
+
name: 'equallySpaced',
|
|
35
|
+
options: { from, to, numberOfPoints, exclusions },
|
|
36
|
+
properties: {},
|
|
37
|
+
});
|
|
38
|
+
for (const filter of filters) {
|
|
39
|
+
// eslint-disable-next-line import/namespace
|
|
40
|
+
const filterFunction = Filters[filter.name];
|
|
41
|
+
if (!filterFunction) {
|
|
42
|
+
throw new Error(`Unknown filter: ${filter.name}`);
|
|
43
|
+
}
|
|
44
|
+
const filterOptions = (filter === null || filter === void 0 ? void 0 : filter.options) || {};
|
|
45
|
+
// @ts-expect-error some method have options and some other ones don't have any options
|
|
46
|
+
const realResult = filterFunction({ x, y: re }, filterOptions);
|
|
47
|
+
x = realResult.data.x;
|
|
48
|
+
re = realResult.data.y;
|
|
49
|
+
if (im) {
|
|
50
|
+
const imaginaryResult = filterFunction({ x: cloneX, y: im },
|
|
51
|
+
// @ts-expect-error some method have options and some other ones don't have any options
|
|
52
|
+
filterOptions);
|
|
53
|
+
im = imaginaryResult.data.y;
|
|
54
|
+
}
|
|
55
|
+
}
|
|
56
|
+
return { x, re, im };
|
|
57
|
+
}
|
|
58
|
+
//# sourceMappingURL=signalProcessing.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"signalProcessing.js","sourceRoot":"","sources":["../../../src/filters/filter1d/signalProcessing.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,OAAO,MAAM,sBAAsB,CAAC;AAMhD,MAAM,CAAC,MAAM,EAAE,GAAG,kBAAkB,CAAC;AACrC,MAAM,CAAC,MAAM,IAAI,GAAG,mBAAmB,CAAC;AAExC,MAAM,CAAC,MAAM,mBAAmB,GAAsC;IACpE,aAAa,EAAE,IAAI;IACnB,aAAa,EAAE,IAAI;CACpB,CAAC;AAEF,MAAM,UAAU,KAAK,CAAC,QAAiB,EAAE,OAAsB;IAC7D,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,EAAE;QAC3B,MAAM,IAAI,KAAK,CAAC,kDAAkD,CAAC,CAAC;KACrE;IAED,QAAQ,CAAC,IAAI,GAAG,QAAQ,CAAC,QAAQ,EAAE,OAAO,CAAC,CAAC;AAC9C,CAAC;AACD,MAAM,UAAU,YAAY,CAAC,QAAiB;IAC5C,IAAI,QAAQ,CAAC,IAAI,CAAC,IAAI;QAAE,OAAO,IAAI,CAAC;IACpC,OAAO,KAAK,CAAC;AACf,CAAC;AAED,MAAM,UAAU,MAAM,CAAC,aAAkB,EAAE,QAAa;IACtD,OAAO;QACL,IAAI,EAAE,IAAI;QACV,MAAM,EAAE,QAAQ;KACjB,CAAC;AACJ,CAAC;AAED,MAAM,UAAU,QAAQ,CAAC,QAAiB,EAAE,OAAsB;;IAChE,IAAI,CAAC,GAAG,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IACjC,IAAI,MAAM,GAAG,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IACtC,IAAI,EAAE,GAAG,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IACnC,IAAI,EAAE,GAAG,MAAA,QAAQ,CAAC,IAAI,CAAC,EAAE,0CAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IAEpC,MAAM,OAAO,GAAG,OAAO,CAAC,OAAO,CAAC,KAAK,EAAE,CAAC;IAExC,MAAM,EACJ,KAAK,EAAE,EAAE,IAAI,EAAE,EAAE,EAAE,EACnB,cAAc,EACd,eAAe,EAAE,UAAU,GAC5B,GAAG,OAAO,CAAC;IAEZ,OAAO,CAAC,IAAI,CAAC;QACX,IAAI,EAAE,eAAe;QACrB,OAAO,EAAE,EAAE,IAAI,EAAE,EAAE,EAAE,cAAc,EAAE,UAAU,EAAE;QACjD,UAAU,EAAE,EAAE;KACf,CAAC,CAAC;IAEH,KAAK,MAAM,MAAM,IAAI,OAAO,EAAE;QAC5B,4CAA4C;QAC5C,MAAM,cAAc,GAAG,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC;QAC5C,IAAI,CAAC,cAAc,EAAE;YACnB,MAAM,IAAI,KAAK,CAAC,mBAAmB,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC;SACnD;QAED,MAAM,aAAa,GAAG,CAAA,MAAM,aAAN,MAAM,uBAAN,MAAM,CAAE,OAAO,KAAI,EAAE,CAAC;QAE5C,uFAAuF;QACvF,MAAM,UAAU,GAAG,cAAc,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,aAAa,CAAC,CAAC;QAC/D,CAAC,GAAG,UAAU,CAAC,IAAI,CAAC,CAAiB,CAAC;QACtC,EAAE,GAAG,UAAU,CAAC,IAAI,CAAC,CAAiB,CAAC;QAEvC,IAAI,EAAE,EAAE;YACN,MAAM,eAAe,GAAG,cAAc,CACpC,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,EAAE;YACpB,uFAAuF;YACvF,aAAa,CACd,CAAC;YACF,EAAE,GAAG,eAAe,CAAC,IAAI,CAAC,CAAiB,CAAC;SAC7C;KACF;IAED,OAAO,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;AACvB,CAAC"}
|
|
@@ -0,0 +1,37 @@
|
|
|
1
|
+
import { xSequentialFill } from 'ml-spectra-processing';
|
|
2
|
+
export function padDataToNextPowerOfTwo(datum1D, digitalFilterApplied) {
|
|
3
|
+
const { x, re, im } = datum1D.data;
|
|
4
|
+
const size = nextPowerOfTwo(x.length);
|
|
5
|
+
let newRE = new Float64Array(size);
|
|
6
|
+
let newIM = new Float64Array(size);
|
|
7
|
+
const pointsToShift = getPointsToShift(datum1D);
|
|
8
|
+
newRE.set(re.slice(0, size - pointsToShift));
|
|
9
|
+
newIM.set(im.slice(0, size - pointsToShift));
|
|
10
|
+
if (pointsToShift > 0 && digitalFilterApplied) {
|
|
11
|
+
newRE.set(re.slice(re.length - pointsToShift), size - pointsToShift);
|
|
12
|
+
newIM.set(im.slice(re.length - pointsToShift), size - pointsToShift);
|
|
13
|
+
}
|
|
14
|
+
const newX = xSequentialFill({
|
|
15
|
+
from: x[0],
|
|
16
|
+
size,
|
|
17
|
+
step: x[1] - x[0],
|
|
18
|
+
});
|
|
19
|
+
datum1D.data = { ...datum1D.data, re: newRE, im: newIM, x: newX };
|
|
20
|
+
}
|
|
21
|
+
function getPointsToShift(datum1D) {
|
|
22
|
+
var _a;
|
|
23
|
+
let grpdly = ((_a = datum1D.info) === null || _a === void 0 ? void 0 : _a.digitalFilter) || 0;
|
|
24
|
+
return grpdly > 0 ? Math.floor(grpdly) : 0;
|
|
25
|
+
}
|
|
26
|
+
function nextPowerOfTwo(n) {
|
|
27
|
+
if (n === 0)
|
|
28
|
+
return 1;
|
|
29
|
+
n--;
|
|
30
|
+
n |= n >> 1;
|
|
31
|
+
n |= n >> 2;
|
|
32
|
+
n |= n >> 4;
|
|
33
|
+
n |= n >> 8;
|
|
34
|
+
n |= n >> 16;
|
|
35
|
+
return n + 1;
|
|
36
|
+
}
|
|
37
|
+
//# sourceMappingURL=padDataToNextPowerOfTwo.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"padDataToNextPowerOfTwo.js","sourceRoot":"","sources":["../../../../src/filters/filter1d/utils/padDataToNextPowerOfTwo.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,eAAe,EAAE,MAAM,uBAAuB,CAAC;AAIxD,MAAM,UAAU,uBAAuB,CACrC,OAAgB,EAChB,oBAA6B;IAE7B,MAAM,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,OAAO,CAAC,IAAI,CAAC;IACnC,MAAM,IAAI,GAAG,cAAc,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;IAEtC,IAAI,KAAK,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,CAAC;IACnC,IAAI,KAAK,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,CAAC;IAEnC,MAAM,aAAa,GAAG,gBAAgB,CAAC,OAAO,CAAC,CAAC;IAEhD,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,GAAG,aAAa,CAAC,CAAC,CAAC;IAC7C,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,GAAG,aAAa,CAAC,CAAC,CAAC;IAE7C,IAAI,aAAa,GAAG,CAAC,IAAI,oBAAoB,EAAE;QAC7C,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,MAAM,GAAG,aAAa,CAAC,EAAE,IAAI,GAAG,aAAa,CAAC,CAAC;QACrE,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,MAAM,GAAG,aAAa,CAAC,EAAE,IAAI,GAAG,aAAa,CAAC,CAAC;KACtE;IAED,MAAM,IAAI,GAAG,eAAe,CAAC;QAC3B,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC;QACV,IAAI;QACJ,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;KAClB,CAAiB,CAAC;IAEnB,OAAO,CAAC,IAAI,GAAG,EAAE,GAAG,OAAO,CAAC,IAAI,EAAE,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC;AACpE,CAAC;AAED,SAAS,gBAAgB,CAAC,OAAgB;;IACxC,IAAI,MAAM,GAAG,CAAA,MAAA,OAAO,CAAC,IAAI,0CAAE,aAAa,KAAI,CAAC,CAAC;IAC9C,OAAO,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAC7C,CAAC;AAED,SAAS,cAAc,CAAC,CAAS;IAC/B,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IACtB,CAAC,EAAE,CAAC;IACJ,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IACZ,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IACZ,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IACZ,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IACZ,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC;IACb,OAAO,CAAC,GAAG,CAAC,CAAC;AACf,CAAC"}
|
|
@@ -0,0 +1,62 @@
|
|
|
1
|
+
export const id = 'zeroFilling';
|
|
2
|
+
export const name = 'Zero Filling';
|
|
3
|
+
export const DOMAIN_UPDATE_RULES = {
|
|
4
|
+
updateXDomain: true,
|
|
5
|
+
updateYDomain: false,
|
|
6
|
+
};
|
|
7
|
+
export function apply(datum1D, options) {
|
|
8
|
+
var _a;
|
|
9
|
+
if (!isApplicable(datum1D)) {
|
|
10
|
+
throw new Error('zeroFilling not applicable on this data');
|
|
11
|
+
}
|
|
12
|
+
let { nbPoints, factor = 2 } = options;
|
|
13
|
+
if (!nbPoints) {
|
|
14
|
+
nbPoints = 2 ** Math.round(Math.log2(datum1D.data.x.length * factor));
|
|
15
|
+
const filter = datum1D.filters.find((filter) => filter.name === id);
|
|
16
|
+
if (filter) {
|
|
17
|
+
filter.value = { nbPoints };
|
|
18
|
+
}
|
|
19
|
+
}
|
|
20
|
+
let digitalFilterApplied = datum1D.filters.some((e) => e.name === 'digitalFilter' && e.flag);
|
|
21
|
+
let grpdly = ((_a = datum1D.info) === null || _a === void 0 ? void 0 : _a.digitalFilter) || 0;
|
|
22
|
+
let pointsToShift;
|
|
23
|
+
if (grpdly > 0 && digitalFilterApplied) {
|
|
24
|
+
pointsToShift = Math.floor(grpdly);
|
|
25
|
+
}
|
|
26
|
+
else {
|
|
27
|
+
pointsToShift = 0;
|
|
28
|
+
}
|
|
29
|
+
const { re, im, x } = datum1D.data;
|
|
30
|
+
let newRE = new Float64Array(nbPoints);
|
|
31
|
+
let newIM = new Float64Array(nbPoints);
|
|
32
|
+
let newX = new Float64Array(nbPoints);
|
|
33
|
+
const length = Math.min(nbPoints, re.length);
|
|
34
|
+
newRE.set(re.slice(0, length - pointsToShift));
|
|
35
|
+
newIM.set(im.slice(0, length - pointsToShift));
|
|
36
|
+
newX.set(x.slice(0, length - pointsToShift));
|
|
37
|
+
let diff = x[1] - x[0];
|
|
38
|
+
let currentX = x[length - pointsToShift - 1];
|
|
39
|
+
for (let i = length - pointsToShift; i < nbPoints; i++) {
|
|
40
|
+
currentX += diff;
|
|
41
|
+
newX[i] = currentX;
|
|
42
|
+
}
|
|
43
|
+
if (pointsToShift > 0 && pointsToShift < nbPoints) {
|
|
44
|
+
newRE.set(re.slice(re.length - pointsToShift), nbPoints - pointsToShift);
|
|
45
|
+
newIM.set(im.slice(re.length - pointsToShift), nbPoints - pointsToShift);
|
|
46
|
+
}
|
|
47
|
+
datum1D.data = { ...datum1D.data, re: newRE, im: newIM, x: newX };
|
|
48
|
+
}
|
|
49
|
+
export function isApplicable(datum1D) {
|
|
50
|
+
if (datum1D.info.isComplex && datum1D.info.isFid) {
|
|
51
|
+
return true;
|
|
52
|
+
}
|
|
53
|
+
return false;
|
|
54
|
+
}
|
|
55
|
+
export function reduce(_previousValue, newValue) {
|
|
56
|
+
delete newValue.factor;
|
|
57
|
+
return {
|
|
58
|
+
once: true,
|
|
59
|
+
reduce: newValue,
|
|
60
|
+
};
|
|
61
|
+
}
|
|
62
|
+
//# sourceMappingURL=zeroFilling.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"zeroFilling.js","sourceRoot":"","sources":["../../../src/filters/filter1d/zeroFilling.ts"],"names":[],"mappings":"AAIA,MAAM,CAAC,MAAM,EAAE,GAAG,aAAa,CAAC;AAChC,MAAM,CAAC,MAAM,IAAI,GAAG,cAAc,CAAC;AACnC,MAAM,CAAC,MAAM,mBAAmB,GAAsC;IACpE,aAAa,EAAE,IAAI;IACnB,aAAa,EAAE,KAAK;CACrB,CAAC;AAOF,MAAM,UAAU,KAAK,CAAC,OAAgB,EAAE,OAA2B;;IACjE,IAAI,CAAC,YAAY,CAAC,OAAO,CAAC,EAAE;QAC1B,MAAM,IAAI,KAAK,CAAC,yCAAyC,CAAC,CAAC;KAC5D;IAED,IAAI,EAAE,QAAQ,EAAE,MAAM,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;IAEvC,IAAI,CAAC,QAAQ,EAAE;QACb,QAAQ,GAAG,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC;QACtE,MAAM,MAAM,GAAG,OAAO,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,MAAM,EAAE,EAAE,CAAC,MAAM,CAAC,IAAI,KAAK,EAAE,CAAC,CAAC;QACpE,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,KAAK,GAAG,EAAE,QAAQ,EAAE,CAAC;SAC7B;KACF;IAED,IAAI,oBAAoB,GAAG,OAAO,CAAC,OAAO,CAAC,IAAI,CAC7C,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,KAAK,eAAe,IAAI,CAAC,CAAC,IAAI,CAC5C,CAAC;IAEF,IAAI,MAAM,GAAG,CAAA,MAAA,OAAO,CAAC,IAAI,0CAAE,aAAa,KAAI,CAAC,CAAC;IAC9C,IAAI,aAAa,CAAC;IAClB,IAAI,MAAM,GAAG,CAAC,IAAI,oBAAoB,EAAE;QACtC,aAAa,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;KACpC;SAAM;QACL,aAAa,GAAG,CAAC,CAAC;KACnB;IAED,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,GAAG,OAAO,CAAC,IAAI,CAAC;IAEnC,IAAI,KAAK,GAAG,IAAI,YAAY,CAAC,QAAQ,CAAC,CAAC;IACvC,IAAI,KAAK,GAAG,IAAI,YAAY,CAAC,QAAQ,CAAC,CAAC;IACvC,IAAI,IAAI,GAAG,IAAI,YAAY,CAAC,QAAQ,CAAC,CAAC;IAEtC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,EAAE,CAAC,MAAM,CAAC,CAAC;IAE7C,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,MAAM,GAAG,aAAa,CAAC,CAAC,CAAC;IAC/C,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,MAAM,GAAG,aAAa,CAAC,CAAC,CAAC;IAC/C,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,MAAM,GAAG,aAAa,CAAC,CAAC,CAAC;IAE7C,IAAI,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACvB,IAAI,QAAQ,GAAG,CAAC,CAAC,MAAM,GAAG,aAAa,GAAG,CAAC,CAAC,CAAC;IAC7C,KAAK,IAAI,CAAC,GAAG,MAAM,GAAG,aAAa,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE;QACtD,QAAQ,IAAI,IAAI,CAAC;QACjB,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC;KACpB;IAED,IAAI,aAAa,GAAG,CAAC,IAAI,aAAa,GAAG,QAAQ,EAAE;QACjD,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,MAAM,GAAG,aAAa,CAAC,EAAE,QAAQ,GAAG,aAAa,CAAC,CAAC;QACzE,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,MAAM,GAAG,aAAa,CAAC,EAAE,QAAQ,GAAG,aAAa,CAAC,CAAC;KAC1E;IAED,OAAO,CAAC,IAAI,GAAG,EAAE,GAAG,OAAO,CAAC,IAAI,EAAE,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC;AACpE,CAAC;AAED,MAAM,UAAU,YAAY,CAC1B,OAAgB;IAEhB,IAAI,OAAO,CAAC,IAAI,CAAC,SAAS,IAAI,OAAO,CAAC,IAAI,CAAC,KAAK,EAAE;QAChD,OAAO,IAAI,CAAC;KACb;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,MAAM,UAAU,MAAM,CACpB,cAAkC,EAClC,QAA4B;IAE5B,OAAO,QAAQ,CAAC,MAAM,CAAC;IAEvB,OAAO;QACL,IAAI,EAAE,IAAI;QACV,MAAM,EAAE,QAAQ;KACjB,CAAC;AACJ,CAAC"}
|
|
@@ -0,0 +1,24 @@
|
|
|
1
|
+
export const id = 'shift2DX';
|
|
2
|
+
export const name = 'Shift 2D X';
|
|
3
|
+
/**
|
|
4
|
+
*
|
|
5
|
+
* @param {Object} datum2d
|
|
6
|
+
* @param {number} [shiftValue=0]
|
|
7
|
+
*/
|
|
8
|
+
export function apply(datum2D, shiftValue = 0) {
|
|
9
|
+
const data = { ...datum2D.data };
|
|
10
|
+
for (const key in data) {
|
|
11
|
+
data[key].minX += shiftValue;
|
|
12
|
+
data[key].maxX += shiftValue;
|
|
13
|
+
}
|
|
14
|
+
}
|
|
15
|
+
export function isApplicable() {
|
|
16
|
+
return true;
|
|
17
|
+
}
|
|
18
|
+
export function reduce(previousValue, newValue) {
|
|
19
|
+
return {
|
|
20
|
+
once: true,
|
|
21
|
+
reduce: previousValue + newValue,
|
|
22
|
+
};
|
|
23
|
+
}
|
|
24
|
+
//# sourceMappingURL=shiftX.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"shiftX.js","sourceRoot":"","sources":["../../../src/filters/filter2d/shiftX.ts"],"names":[],"mappings":"AAIA,MAAM,CAAC,MAAM,EAAE,GAAG,UAAU,CAAC;AAC7B,MAAM,CAAC,MAAM,IAAI,GAAG,YAAY,CAAC;AAEjC;;;;GAIG;AAEH,MAAM,UAAU,KAAK,CAAC,OAAgB,EAAE,UAAU,GAAG,CAAC;IACpD,MAAM,IAAI,GAAkC,EAAE,GAAG,OAAO,CAAC,IAAI,EAAE,CAAC;IAChE,KAAK,MAAM,GAAG,IAAI,IAAI,EAAE;QACtB,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,IAAI,UAAU,CAAC;QAC7B,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,IAAI,UAAU,CAAC;KAC9B;AACH,CAAC;AAED,MAAM,UAAU,YAAY;IAC1B,OAAO,IAAI,CAAC;AACd,CAAC;AAED,MAAM,UAAU,MAAM,CAAC,aAAqB,EAAE,QAAgB;IAC5D,OAAO;QACL,IAAI,EAAE,IAAI;QACV,MAAM,EAAE,aAAa,GAAG,QAAQ;KACjC,CAAC;AACJ,CAAC"}
|
|
@@ -0,0 +1,24 @@
|
|
|
1
|
+
export const id = 'shift2DY';
|
|
2
|
+
export const name = 'Shift 2D Y';
|
|
3
|
+
/**
|
|
4
|
+
*
|
|
5
|
+
* @param {Object} datum2d
|
|
6
|
+
* @param {number} [shiftValue=0]
|
|
7
|
+
*/
|
|
8
|
+
export function apply(datum2D, shiftValue = 0) {
|
|
9
|
+
const data = { ...datum2D.data };
|
|
10
|
+
for (const key in data) {
|
|
11
|
+
data[key].minY += shiftValue;
|
|
12
|
+
data[key].maxY += shiftValue;
|
|
13
|
+
}
|
|
14
|
+
}
|
|
15
|
+
export function isApplicable() {
|
|
16
|
+
return true;
|
|
17
|
+
}
|
|
18
|
+
export function reduce(previousValue, newValue) {
|
|
19
|
+
return {
|
|
20
|
+
once: true,
|
|
21
|
+
reduce: previousValue + newValue,
|
|
22
|
+
};
|
|
23
|
+
}
|
|
24
|
+
//# sourceMappingURL=shiftY.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"shiftY.js","sourceRoot":"","sources":["../../../src/filters/filter2d/shiftY.ts"],"names":[],"mappings":"AAIA,MAAM,CAAC,MAAM,EAAE,GAAG,UAAU,CAAC;AAC7B,MAAM,CAAC,MAAM,IAAI,GAAG,YAAY,CAAC;AACjC;;;;GAIG;AACH,MAAM,UAAU,KAAK,CAAC,OAAgB,EAAE,UAAU,GAAG,CAAC;IACpD,MAAM,IAAI,GAAkC,EAAE,GAAG,OAAO,CAAC,IAAI,EAAE,CAAC;IAChE,KAAK,MAAM,GAAG,IAAI,IAAI,EAAE;QACtB,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,IAAI,UAAU,CAAC;QAC7B,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,IAAI,UAAU,CAAC;KAC9B;AACH,CAAC;AAED,MAAM,UAAU,YAAY;IAC1B,OAAO,IAAI,CAAC;AACd,CAAC;AAED,MAAM,UAAU,MAAM,CAAC,aAAqB,EAAE,QAAgB;IAC5D,OAAO;QACL,IAAI,EAAE,IAAI;QACV,MAAM,EAAE,aAAa,GAAG,QAAQ;KACjC,CAAC;AACJ,CAAC"}
|
package/lib-esm/index.js
CHANGED
|
@@ -32,6 +32,7 @@ export * from './xy/xyPeaksOptimization';
|
|
|
32
32
|
export * from './xyz/xyzAutoZonesPicking';
|
|
33
33
|
export * from './xyz/xyzAutoSignalsPicking';
|
|
34
34
|
export * from './xyz/xyzJResAnalyzer';
|
|
35
|
+
export * from './xyz/xyzBidimensionalFFT';
|
|
35
36
|
export * from './databases/getDatabase';
|
|
36
37
|
export * from './databases/carbonImpurities';
|
|
37
38
|
export * from './databases/protonImpurities';
|
|
@@ -44,4 +45,5 @@ export * from './multiplicity/splitPatterns';
|
|
|
44
45
|
export * from './multiplicity/checkMultiplet';
|
|
45
46
|
export * from './multiplicity/checkMultiplicity';
|
|
46
47
|
export * from './multiplicity/translateMultiplet';
|
|
48
|
+
export * from './filters/FiltersManager';
|
|
47
49
|
//# sourceMappingURL=index.js.map
|
package/lib-esm/index.js.map
CHANGED
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA,cAAc,oBAAoB,CAAC;AACnC,cAAc,wBAAwB,CAAC;AACvC,cAAc,8BAA8B,CAAC;AAE7C,cAAc,+BAA+B,CAAC;AAC9C,cAAc,mBAAmB,CAAC;AAClC,cAAc,kBAAkB,CAAC;AAEjC,cAAc,4BAA4B,CAAC;AAC3C,cAAc,4BAA4B,CAAC;AAC3C,cAAc,0BAA0B,CAAC;AACzC,cAAc,0BAA0B,CAAC;AACzC,cAAc,0BAA0B,CAAC;AACzC,cAAc,yBAAyB,CAAC;AACxC,cAAc,gCAAgC,CAAC;AAE/C,cAAc,+BAA+B,CAAC;AAC9C,cAAc,gCAAgC,CAAC;AAC/C,cAAc,6BAA6B,CAAC;AAE5C,cAAc,sBAAsB,CAAC;AACrC,cAAc,qBAAqB,CAAC;AAEpC,cAAc,2BAA2B,CAAC;AAC1C,cAAc,uBAAuB,CAAC;AACtC,cAAc,wBAAwB,CAAC;AACvC,cAAc,2BAA2B,CAAC;AAC1C,cAAc,wBAAwB,CAAC;AACvC,cAAc,8BAA8B,CAAC;AAE7C,cAAc,uBAAuB,CAAC;AACtC,cAAc,6BAA6B,CAAC;AAC5C,cAAc,0BAA0B,CAAC;AACzC,cAAc,8BAA8B,CAAC;AAE7C,cAAc,yBAAyB,CAAC;AACxC,cAAc,0BAA0B,CAAC;AACzC,cAAc,0BAA0B,CAAC;AAEzC,cAAc,2BAA2B,CAAC;AAC1C,cAAc,6BAA6B,CAAC;AAC5C,cAAc,uBAAuB,CAAC;
|
|
1
|
+
{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA,cAAc,oBAAoB,CAAC;AACnC,cAAc,wBAAwB,CAAC;AACvC,cAAc,8BAA8B,CAAC;AAE7C,cAAc,+BAA+B,CAAC;AAC9C,cAAc,mBAAmB,CAAC;AAClC,cAAc,kBAAkB,CAAC;AAEjC,cAAc,4BAA4B,CAAC;AAC3C,cAAc,4BAA4B,CAAC;AAC3C,cAAc,0BAA0B,CAAC;AACzC,cAAc,0BAA0B,CAAC;AACzC,cAAc,0BAA0B,CAAC;AACzC,cAAc,yBAAyB,CAAC;AACxC,cAAc,gCAAgC,CAAC;AAE/C,cAAc,+BAA+B,CAAC;AAC9C,cAAc,gCAAgC,CAAC;AAC/C,cAAc,6BAA6B,CAAC;AAE5C,cAAc,sBAAsB,CAAC;AACrC,cAAc,qBAAqB,CAAC;AAEpC,cAAc,2BAA2B,CAAC;AAC1C,cAAc,uBAAuB,CAAC;AACtC,cAAc,wBAAwB,CAAC;AACvC,cAAc,2BAA2B,CAAC;AAC1C,cAAc,wBAAwB,CAAC;AACvC,cAAc,8BAA8B,CAAC;AAE7C,cAAc,uBAAuB,CAAC;AACtC,cAAc,6BAA6B,CAAC;AAC5C,cAAc,0BAA0B,CAAC;AACzC,cAAc,8BAA8B,CAAC;AAE7C,cAAc,yBAAyB,CAAC;AACxC,cAAc,0BAA0B,CAAC;AACzC,cAAc,0BAA0B,CAAC;AAEzC,cAAc,2BAA2B,CAAC;AAC1C,cAAc,6BAA6B,CAAC;AAC5C,cAAc,uBAAuB,CAAC;AACtC,cAAc,2BAA2B,CAAC;AAE1C,cAAc,yBAAyB,CAAC;AACxC,cAAc,8BAA8B,CAAC;AAC7C,cAAc,8BAA8B,CAAC;AAE7C,cAAc,4BAA4B,CAAC;AAC3C,cAAc,4BAA4B,CAAC;AAE3C,cAAc,2BAA2B,CAAC;AAE1C,cAAc,kCAAkC,CAAC;AAEjD,cAAc,8BAA8B,CAAC;AAC7C,cAAc,8BAA8B,CAAC;AAC7C,cAAc,+BAA+B,CAAC;AAC9C,cAAc,kCAAkC,CAAC;AACjD,cAAc,mCAAmC,CAAC;AAElD,cAAc,0BAA0B,CAAC"}
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"checkData2DFid.js","sourceRoot":"","sources":["../../src/utilities/checkData2DFid.ts"],"names":[],"mappings":"AAGA,MAAM,UAAU,cAAc,CAC5B,IAAY;IAEZ,IAAI,CAAC,CAAC,IAAI,IAAI,IAAI,IAAI,IAAI,IAAI,IAAI,CAAC,EAAE;QACnC,MAAM,IAAI,KAAK,CAAC,mCAAmC,CAAC,CAAC;KACtD;AACH,CAAC"}
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"cloneData1D.js","sourceRoot":"","sources":["../../src/utilities/cloneData1D.ts"],"names":[],"mappings":"AAEA,MAAM,UAAU,WAAW,CAAC,KAAa;IACvC,MAAM,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,KAAK,IAAI,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;IACzD,OAAO,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,KAAK,EAAE,EAAE,EAAE,EAAE,CAAA,EAAE,aAAF,EAAE,uBAAF,EAAE,CAAE,KAAK,EAAE,KAAI,IAAI,EAAE,CAAC;AACnE,CAAC"}
|
|
@@ -0,0 +1,24 @@
|
|
|
1
|
+
export function cloneData2D(data) {
|
|
2
|
+
const newData = {};
|
|
3
|
+
if (isData2DFt(data)) {
|
|
4
|
+
for (const key in data) {
|
|
5
|
+
addCloneData(newData, key, data[key]);
|
|
6
|
+
}
|
|
7
|
+
}
|
|
8
|
+
else {
|
|
9
|
+
for (const key in data) {
|
|
10
|
+
addCloneData(newData, key, data[key]);
|
|
11
|
+
}
|
|
12
|
+
}
|
|
13
|
+
return newData;
|
|
14
|
+
}
|
|
15
|
+
function isData2DFt(data) {
|
|
16
|
+
return 're' in data;
|
|
17
|
+
}
|
|
18
|
+
function addCloneData(newData, key, data) {
|
|
19
|
+
if (data) {
|
|
20
|
+
const { z, ...other } = data;
|
|
21
|
+
newData[key] = { z: z.slice(), ...other };
|
|
22
|
+
}
|
|
23
|
+
}
|
|
24
|
+
//# sourceMappingURL=cloneData2D.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"cloneData2D.js","sourceRoot":"","sources":["../../src/utilities/cloneData2D.ts"],"names":[],"mappings":"AAEA,MAAM,UAAU,WAAW,CAAC,IAAY;IACtC,MAAM,OAAO,GAAQ,EAAE,CAAC;IAExB,IAAI,UAAU,CAAC,IAAI,CAAC,EAAE;QACpB,KAAK,MAAM,GAAG,IAAI,IAAI,EAAE;YACtB,YAAY,CAAC,OAAO,EAAE,GAAG,EAAE,IAAI,CAAC,GAAqB,CAAC,CAAC,CAAC;SACzD;KACF;SAAM;QACL,KAAK,MAAM,GAAG,IAAI,IAAI,EAAE;YACtB,YAAY,CAAC,OAAO,EAAE,GAAG,EAAE,IAAI,CAAC,GAAsB,CAAC,CAAC,CAAC;SAC1D;KACF;IACD,OAAO,OAAO,CAAC;AACjB,CAAC;AAED,SAAS,UAAU,CAAC,IAAY;IAC9B,OAAO,IAAI,IAAI,IAAI,CAAC;AACtB,CAAC;AAED,SAAS,YAAY,CAAC,OAAY,EAAE,GAAW,EAAE,IAAoB;IACnE,IAAI,IAAI,EAAE;QACR,MAAM,EAAE,CAAC,EAAE,GAAG,KAAK,EAAE,GAAG,IAAI,CAAC;QAC7B,OAAO,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,EAAE,EAAE,GAAG,KAAK,EAAE,CAAC;KAC3C;AACH,CAAC"}
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"Entry2D.js","sourceRoot":"","sources":["../../src/xyz/Entry2D.ts"],"names":[],"mappings":""}
|
|
@@ -0,0 +1,50 @@
|
|
|
1
|
+
import { xAdd, xMultiply, xSubtract } from 'ml-spectra-processing';
|
|
2
|
+
import { checkData2DFid } from '../utilities/checkData2DFid';
|
|
3
|
+
export function quadrature(data, options) {
|
|
4
|
+
checkData2DFid(data);
|
|
5
|
+
const { fnMode } = options;
|
|
6
|
+
switch (fnMode) {
|
|
7
|
+
case 'ttpi':
|
|
8
|
+
return ttpiQuadrature(data);
|
|
9
|
+
case 'states':
|
|
10
|
+
return shrQuadrature(data);
|
|
11
|
+
case 'echo':
|
|
12
|
+
return echoQuadrature(data);
|
|
13
|
+
default:
|
|
14
|
+
return {
|
|
15
|
+
re: data.re.z,
|
|
16
|
+
im: data.im.z,
|
|
17
|
+
};
|
|
18
|
+
}
|
|
19
|
+
}
|
|
20
|
+
function echoQuadrature(data) {
|
|
21
|
+
const { re: { z: reMatrix }, im: { z: imMatrix }, } = data;
|
|
22
|
+
const newRe = [];
|
|
23
|
+
const newIm = [];
|
|
24
|
+
for (let i = 0; i < reMatrix.length; i += 2) {
|
|
25
|
+
let reCosine = reMatrix[i];
|
|
26
|
+
let imCosine = imMatrix[i];
|
|
27
|
+
let reSine = reMatrix[i + 1];
|
|
28
|
+
let imSine = imMatrix[i + 1];
|
|
29
|
+
newRe.push(xSubtract(xMultiply(imSine, -1), imCosine));
|
|
30
|
+
newRe.push(xSubtract(reSine, reCosine));
|
|
31
|
+
newIm.push(xAdd(reSine, reCosine));
|
|
32
|
+
newIm.push(xSubtract(imSine, imCosine));
|
|
33
|
+
}
|
|
34
|
+
return { re: newRe, im: newIm };
|
|
35
|
+
}
|
|
36
|
+
function shrQuadrature(data) {
|
|
37
|
+
const { re, im } = data;
|
|
38
|
+
return { re: re.z, im: im.z };
|
|
39
|
+
}
|
|
40
|
+
function ttpiQuadrature(data) {
|
|
41
|
+
const newRe = [];
|
|
42
|
+
const newIm = [];
|
|
43
|
+
const { re: { z: reMatrix }, im: { z: imMatrix }, } = data;
|
|
44
|
+
for (let i = 0; i < reMatrix.length; i++) {
|
|
45
|
+
newRe.push(i % 2 ? reMatrix[i].slice() : xMultiply(reMatrix[i], -1));
|
|
46
|
+
newIm.push(i % 2 ? imMatrix[i].slice() : xMultiply(imMatrix[i], -1));
|
|
47
|
+
}
|
|
48
|
+
return { re: newRe, im: newIm };
|
|
49
|
+
}
|
|
50
|
+
//# sourceMappingURL=quadrature.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"quadrature.js","sourceRoot":"","sources":["../../src/xyz/quadrature.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,IAAI,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,uBAAuB,CAAC;AAEnE,OAAO,EAAE,cAAc,EAAqB,MAAM,6BAA6B,CAAC;AAEhF,MAAM,UAAU,UAAU,CAAC,IAAY,EAAE,OAA2B;IAClE,cAAc,CAAC,IAAI,CAAC,CAAC;IAErB,MAAM,EAAE,MAAM,EAAE,GAAG,OAAO,CAAC;IAE3B,QAAQ,MAAM,EAAE;QACd,KAAK,MAAM;YACT,OAAO,cAAc,CAAC,IAAI,CAAC,CAAC;QAC9B,KAAK,QAAQ;YACX,OAAO,aAAa,CAAC,IAAI,CAAC,CAAC;QAC7B,KAAK,MAAM;YACT,OAAO,cAAc,CAAC,IAAI,CAAC,CAAC;QAC9B;YACE,OAAO;gBACL,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;gBACb,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;aACd,CAAC;KACL;AACH,CAAC;AAED,SAAS,cAAc,CAAC,IAAuB;IAC7C,MAAM,EACJ,EAAE,EAAE,EAAE,CAAC,EAAE,QAAQ,EAAE,EACnB,EAAE,EAAE,EAAE,CAAC,EAAE,QAAQ,EAAE,GACpB,GAAG,IAAI,CAAC;IAET,MAAM,KAAK,GAAkB,EAAE,CAAC;IAChC,MAAM,KAAK,GAAkB,EAAE,CAAC;IAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC,EAAE;QAC3C,IAAI,QAAQ,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;QAC3B,IAAI,QAAQ,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;QAC3B,IAAI,MAAM,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC7B,IAAI,MAAM,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAE7B,KAAK,CAAC,IAAI,CAAC,SAAS,CAAC,SAAS,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC;QACvD,KAAK,CAAC,IAAI,CAAC,SAAS,CAAC,MAAM,EAAE,QAAQ,CAAC,CAAC,CAAC;QACxC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,EAAE,QAAQ,CAAC,CAAC,CAAC;QACnC,KAAK,CAAC,IAAI,CAAC,SAAS,CAAC,MAAM,EAAE,QAAQ,CAAC,CAAC,CAAC;KACzC;IAED,OAAO,EAAE,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,EAAE,CAAC;AAClC,CAAC;AAED,SAAS,aAAa,CAAC,IAAuB;IAC5C,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,IAAI,CAAC;IACxB,OAAO,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC;AAChC,CAAC;AAED,SAAS,cAAc,CAAC,IAAuB;IAC7C,MAAM,KAAK,GAAkB,EAAE,CAAC;IAChC,MAAM,KAAK,GAAkB,EAAE,CAAC;IAEhC,MAAM,EACJ,EAAE,EAAE,EAAE,CAAC,EAAE,QAAQ,EAAE,EACnB,EAAE,EAAE,EAAE,CAAC,EAAE,QAAQ,EAAE,GACpB,GAAG,IAAI,CAAC;IACT,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACxC,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,SAAS,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QACrE,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,SAAS,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;KACtE;IAED,OAAO,EAAE,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,EAAE,CAAC;AAClC,CAAC"}
|
|
@@ -0,0 +1,16 @@
|
|
|
1
|
+
export function digitalFilter(data, options) {
|
|
2
|
+
let { digitalFilterValue = 0 } = options;
|
|
3
|
+
let re = new Float64Array(data.re);
|
|
4
|
+
let im = new Float64Array(data.im);
|
|
5
|
+
let pointsToShift = Math.floor(digitalFilterValue);
|
|
6
|
+
const skip = 0;
|
|
7
|
+
pointsToShift += 0;
|
|
8
|
+
const newRe = new Float64Array(re.length);
|
|
9
|
+
const newIm = new Float64Array(im.length);
|
|
10
|
+
newRe.set(re.slice(pointsToShift));
|
|
11
|
+
newRe.set(re.slice(skip, pointsToShift), re.length - pointsToShift);
|
|
12
|
+
newIm.set(im.slice(pointsToShift));
|
|
13
|
+
newIm.set(im.slice(skip, pointsToShift), im.length - pointsToShift);
|
|
14
|
+
return { re: newRe, im: newIm };
|
|
15
|
+
}
|
|
16
|
+
//# sourceMappingURL=digitalFilter.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"digitalFilter.js","sourceRoot":"","sources":["../../../../src/xyz/util/fft2d/digitalFilter.ts"],"names":[],"mappings":"AAEA,MAAM,UAAU,aAAa,CAC3B,IAAc,EACd,OAAuC;IAEvC,IAAI,EAAE,kBAAkB,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;IACzC,IAAI,EAAE,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IACnC,IAAI,EAAE,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAEnC,IAAI,aAAa,GAAG,IAAI,CAAC,KAAK,CAAC,kBAAkB,CAAC,CAAC;IAEnD,MAAM,IAAI,GAAG,CAAC,CAAC;IACf,aAAa,IAAI,CAAC,CAAC;IAEnB,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC;IAC1C,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC;IAC1C,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,aAAa,CAAC,CAAC,CAAC;IACnC,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,IAAI,EAAE,aAAa,CAAC,EAAE,EAAE,CAAC,MAAM,GAAG,aAAa,CAAC,CAAC;IACpE,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,aAAa,CAAC,CAAC,CAAC;IACnC,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,IAAI,EAAE,aAAa,CAAC,EAAE,EAAE,CAAC,MAAM,GAAG,aAAa,CAAC,CAAC;IAEpE,OAAO,EAAE,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,EAAE,CAAC;AAClC,CAAC"}
|
|
@@ -0,0 +1,50 @@
|
|
|
1
|
+
import Matrix from 'ml-matrix';
|
|
2
|
+
import { reimAutoPhaseCorrection, reimFFT, reimPhaseCorrection, } from 'ml-spectra-processing';
|
|
3
|
+
import { apodization, } from '../../../apodization/apodization';
|
|
4
|
+
import { digitalFilter } from './digitalFilter';
|
|
5
|
+
import { removeDCOffset } from './removeDCOffset';
|
|
6
|
+
import { zeroFilling } from './zeroFilling';
|
|
7
|
+
export function fftDirectDimension(data, options) {
|
|
8
|
+
const { re: reData, im: imData } = data;
|
|
9
|
+
const nbRows = reData.length;
|
|
10
|
+
const nbColumns = reData[0].length;
|
|
11
|
+
const { zeroFilling: zeroFillingF2, digitalFilterValue = 0 } = options;
|
|
12
|
+
const { factor: f2Factor = 2 } = zeroFillingF2 || {};
|
|
13
|
+
const nbPointsF2 = (zeroFillingF2 === null || zeroFillingF2 === void 0 ? void 0 : zeroFillingF2.nbPoints) || 2 ** Math.round(Math.log2(nbColumns * f2Factor));
|
|
14
|
+
const reFFT = new Matrix(nbRows, nbPointsF2);
|
|
15
|
+
const imFFT = new Matrix(nbRows, nbPointsF2);
|
|
16
|
+
const { apodization: apodizationF2 } = options;
|
|
17
|
+
let { ph0: ph0F2, ph1: ph1F2, mode: pcModeF2 = 'no', } = (options === null || options === void 0 ? void 0 : options.phaseCorrection) || {};
|
|
18
|
+
for (let i = 0; i < nbRows; i++) {
|
|
19
|
+
const newData = digitalFilter({ re: reData[i], im: imData[i] }, { digitalFilterValue });
|
|
20
|
+
Object.assign(newData, removeDCOffset(newData, { digitalFilterValue }));
|
|
21
|
+
if (apodizationF2) {
|
|
22
|
+
const { re: aRe, im: aIm } = apodization(newData, apodizationF2);
|
|
23
|
+
Object.assign(newData, { re: aRe, im: aIm });
|
|
24
|
+
}
|
|
25
|
+
Object.assign(newData, zeroFilling(newData, { digitalFilterValue, nbPoints: nbPointsF2 }));
|
|
26
|
+
Object.assign(newData, reimFFT(newData, { applyZeroShift: true }));
|
|
27
|
+
const dfResidual = (digitalFilterValue - Math.round(digitalFilterValue)) * Math.PI * 2;
|
|
28
|
+
if (dfResidual > 0) {
|
|
29
|
+
Object.assign(newData, reimPhaseCorrection(newData, 0, -dfResidual));
|
|
30
|
+
}
|
|
31
|
+
if (i === 0 &&
|
|
32
|
+
pcModeF2 === 'pk' &&
|
|
33
|
+
ph0F2 === undefined &&
|
|
34
|
+
ph1F2 === undefined) {
|
|
35
|
+
const phased = reimAutoPhaseCorrection(newData, { magnitudeMode: false });
|
|
36
|
+
ph0F2 = phased.ph0;
|
|
37
|
+
ph1F2 = phased.ph1;
|
|
38
|
+
}
|
|
39
|
+
if (pcModeF2 === 'pk' && ph0F2 !== undefined && ph1F2 !== undefined) {
|
|
40
|
+
Object.assign(newData, reimPhaseCorrection(newData, (ph0F2 * Math.PI) / 180, (ph1F2 * Math.PI) / 180));
|
|
41
|
+
}
|
|
42
|
+
reFFT.setRow(i, newData.re);
|
|
43
|
+
imFFT.setRow(i, newData.im);
|
|
44
|
+
}
|
|
45
|
+
return {
|
|
46
|
+
re: reFFT,
|
|
47
|
+
im: imFFT,
|
|
48
|
+
};
|
|
49
|
+
}
|
|
50
|
+
//# sourceMappingURL=fftDirectDimension.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"fftDirectDimension.js","sourceRoot":"","sources":["../../../../src/xyz/util/fft2d/fftDirectDimension.ts"],"names":[],"mappings":"AACA,OAAO,MAAM,MAAM,WAAW,CAAC;AAC/B,OAAO,EACL,uBAAuB,EACvB,OAAO,EACP,mBAAmB,GACpB,MAAM,uBAAuB,CAAC;AAE/B,OAAO,EACL,WAAW,GAEZ,MAAM,kCAAkC,CAAC;AAE1C,OAAO,EAAE,aAAa,EAAE,MAAM,iBAAiB,CAAC;AAChD,OAAO,EAAE,cAAc,EAAE,MAAM,kBAAkB,CAAC;AAClD,OAAO,EAAE,WAAW,EAAsB,MAAM,eAAe,CAAC;AAgBhE,MAAM,UAAU,kBAAkB,CAChC,IAA8C,EAC9C,OAA+B;IAE/B,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,IAAI,CAAC;IACxC,MAAM,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC;IAC7B,MAAM,SAAS,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;IAEnC,MAAM,EAAE,WAAW,EAAE,aAAa,EAAE,kBAAkB,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;IAEvE,MAAM,EAAE,MAAM,EAAE,QAAQ,GAAG,CAAC,EAAE,GAAG,aAAa,IAAI,EAAE,CAAC;IACrD,MAAM,UAAU,GACd,CAAA,aAAa,aAAb,aAAa,uBAAb,aAAa,CAAE,QAAQ,KAAI,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,SAAS,GAAG,QAAQ,CAAC,CAAC,CAAC;IAE9E,MAAM,KAAK,GAAG,IAAI,MAAM,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC;IAC7C,MAAM,KAAK,GAAG,IAAI,MAAM,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC;IAE7C,MAAM,EAAE,WAAW,EAAE,aAAa,EAAE,GAAG,OAAO,CAAC;IAC/C,IAAI,EACF,GAAG,EAAE,KAAK,EACV,GAAG,EAAE,KAAK,EACV,IAAI,EAAE,QAAQ,GAAG,IAAI,GACtB,GAAG,CAAA,OAAO,aAAP,OAAO,uBAAP,OAAO,CAAE,eAAe,KAAI,EAAE,CAAC;IAEnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE;QAC/B,MAAM,OAAO,GAAG,aAAa,CAC3B,EAAE,EAAE,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,EAChC,EAAE,kBAAkB,EAAE,CACvB,CAAC;QAEF,MAAM,CAAC,MAAM,CAAC,OAAO,EAAE,cAAc,CAAC,OAAO,EAAE,EAAE,kBAAkB,EAAE,CAAC,CAAC,CAAC;QAExE,IAAI,aAAa,EAAE;YACjB,MAAM,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,WAAW,CAAC,OAAO,EAAE,aAAa,CAAC,CAAC;YACjE,MAAM,CAAC,MAAM,CAAC,OAAO,EAAE,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC;SAC9C;QAED,MAAM,CAAC,MAAM,CACX,OAAO,EACP,WAAW,CAAC,OAAO,EAAE,EAAE,kBAAkB,EAAE,QAAQ,EAAE,UAAU,EAAE,CAAC,CACnE,CAAC;QAEF,MAAM,CAAC,MAAM,CAAC,OAAO,EAAE,OAAO,CAAC,OAAO,EAAE,EAAE,cAAc,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC;QAEnE,MAAM,UAAU,GACd,CAAC,kBAAkB,GAAG,IAAI,CAAC,KAAK,CAAC,kBAAkB,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;QACtE,IAAI,UAAU,GAAG,CAAC,EAAE;YAClB,MAAM,CAAC,MAAM,CAAC,OAAO,EAAE,mBAAmB,CAAC,OAAO,EAAE,CAAC,EAAE,CAAC,UAAU,CAAC,CAAC,CAAC;SACtE;QAED,IACE,CAAC,KAAK,CAAC;YACP,QAAQ,KAAK,IAAI;YACjB,KAAK,KAAK,SAAS;YACnB,KAAK,KAAK,SAAS,EACnB;YACA,MAAM,MAAM,GAAG,uBAAuB,CAAC,OAAO,EAAE,EAAE,aAAa,EAAE,KAAK,EAAE,CAAC,CAAC;YAC1E,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC;YACnB,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC;SACpB;QACD,IAAI,QAAQ,KAAK,IAAI,IAAI,KAAK,KAAK,SAAS,IAAI,KAAK,KAAK,SAAS,EAAE;YACnE,MAAM,CAAC,MAAM,CACX,OAAO,EACP,mBAAmB,CACjB,OAAO,EACP,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,EACvB,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CACxB,CACF,CAAC;SACH;QAED,KAAK,CAAC,MAAM,CAAC,CAAC,EAAE,OAAO,CAAC,EAAE,CAAC,CAAC;QAC5B,KAAK,CAAC,MAAM,CAAC,CAAC,EAAE,OAAO,CAAC,EAAE,CAAC,CAAC;KAC7B;IAED,OAAO;QACL,EAAE,EAAE,KAAK;QACT,EAAE,EAAE,KAAK;KACV,CAAC;AACJ,CAAC"}
|
|
@@ -0,0 +1,65 @@
|
|
|
1
|
+
import Matrix from 'ml-matrix';
|
|
2
|
+
import { reimAbsolute, reimAutoPhaseCorrection, reimFFT, reimPhaseCorrection, } from 'ml-spectra-processing';
|
|
3
|
+
import { apodization, } from '../../../apodization/apodization';
|
|
4
|
+
import { zeroFilling } from './zeroFilling';
|
|
5
|
+
export function fftIndirectDimension(data, options) {
|
|
6
|
+
var _a;
|
|
7
|
+
const reData = Matrix.checkMatrix(data.re);
|
|
8
|
+
const imData = Matrix.checkMatrix(data.im);
|
|
9
|
+
const nbRows = reData.rows;
|
|
10
|
+
const nbColumns = reData.columns;
|
|
11
|
+
const { zeroFilling: zeroFillingF1 } = options;
|
|
12
|
+
const { factor: f1Factor = 2 } = zeroFillingF1 || {};
|
|
13
|
+
const nbPointsF1 = (zeroFillingF1 === null || zeroFillingF1 === void 0 ? void 0 : zeroFillingF1.nbPoints) || 2 ** Math.round(Math.log2(nbRows * f1Factor));
|
|
14
|
+
const { apodization: apodizationF1 } = options;
|
|
15
|
+
let { ph0: ph0F1, ph1: ph1F1, mode: pcModeF1, } = options.phaseCorrection || {};
|
|
16
|
+
const reResult = new Matrix(nbPointsF1, nbColumns);
|
|
17
|
+
const imResult = new Matrix(nbPointsF1, nbColumns);
|
|
18
|
+
for (let i = 0; i < nbColumns; i++) {
|
|
19
|
+
const { re: aRe, im: aIm } = apodization({
|
|
20
|
+
re: reData.getColumn(i),
|
|
21
|
+
im: imData.getColumn(i),
|
|
22
|
+
}, apodizationF1);
|
|
23
|
+
const newRe = [];
|
|
24
|
+
const newIm = [];
|
|
25
|
+
if (options.fnMode === 'echo') {
|
|
26
|
+
for (let j = 0; j < aRe.length; j += 2) {
|
|
27
|
+
newRe.push(aRe[j]);
|
|
28
|
+
newIm.push(aRe[j + 1]);
|
|
29
|
+
}
|
|
30
|
+
}
|
|
31
|
+
else {
|
|
32
|
+
for (let j = 0; j < aRe.length; j++) {
|
|
33
|
+
newRe.push(aRe[j]);
|
|
34
|
+
newIm.push(-aIm[j]);
|
|
35
|
+
}
|
|
36
|
+
}
|
|
37
|
+
let transformed = reimFFT(zeroFilling({
|
|
38
|
+
re: newRe,
|
|
39
|
+
im: newIm,
|
|
40
|
+
}, { digitalFilterValue: 0, nbPoints: nbPointsF1 }), { applyZeroShift: true });
|
|
41
|
+
let { re, im } = transformed;
|
|
42
|
+
if (((_a = options.phaseCorrection) === null || _a === void 0 ? void 0 : _a.mode) === 'mc') {
|
|
43
|
+
re = reimAbsolute({ re, im });
|
|
44
|
+
im = new Float64Array(im.length);
|
|
45
|
+
}
|
|
46
|
+
if (pcModeF1 === 'pk') {
|
|
47
|
+
if (i === 0 && ph0F1 === undefined && ph1F1 === undefined) {
|
|
48
|
+
//here we can improve it by checking if the experiment is phase sensitive.
|
|
49
|
+
const phased = reimAutoPhaseCorrection({ re, im }, { magnitudeMode: false });
|
|
50
|
+
ph0F1 = phased.ph0;
|
|
51
|
+
ph1F1 = phased.ph1;
|
|
52
|
+
}
|
|
53
|
+
if (ph0F1 !== undefined && ph1F1 !== undefined) {
|
|
54
|
+
Object.assign(transformed, reimPhaseCorrection(transformed, (ph0F1 * Math.PI) / 180, (ph1F1 * Math.PI) / 180));
|
|
55
|
+
}
|
|
56
|
+
}
|
|
57
|
+
reResult.setColumn(i, re);
|
|
58
|
+
imResult.setColumn(i, im);
|
|
59
|
+
}
|
|
60
|
+
return {
|
|
61
|
+
re: reResult,
|
|
62
|
+
im: imResult,
|
|
63
|
+
};
|
|
64
|
+
}
|
|
65
|
+
//# sourceMappingURL=fftIndirectDimension.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"fftIndirectDimension.js","sourceRoot":"","sources":["../../../../src/xyz/util/fft2d/fftIndirectDimension.ts"],"names":[],"mappings":"AACA,OAAO,MAAM,MAAM,WAAW,CAAC;AAC/B,OAAO,EACL,YAAY,EACZ,uBAAuB,EACvB,OAAO,EACP,mBAAmB,GACpB,MAAM,uBAAuB,CAAC;AAE/B,OAAO,EACL,WAAW,GAEZ,MAAM,kCAAkC,CAAC;AAE1C,OAAO,EAAE,WAAW,EAAsB,MAAM,eAAe,CAAC;AAehE,MAAM,UAAU,oBAAoB,CAClC,IAAgE,EAChE,OAAiC;;IAEjC,MAAM,MAAM,GAAG,MAAM,CAAC,WAAW,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAC3C,MAAM,MAAM,GAAG,MAAM,CAAC,WAAW,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAE3C,MAAM,MAAM,GAAG,MAAM,CAAC,IAAI,CAAC;IAC3B,MAAM,SAAS,GAAG,MAAM,CAAC,OAAO,CAAC;IAEjC,MAAM,EAAE,WAAW,EAAE,aAAa,EAAE,GAAG,OAAO,CAAC;IAC/C,MAAM,EAAE,MAAM,EAAE,QAAQ,GAAG,CAAC,EAAE,GAAG,aAAa,IAAI,EAAE,CAAC;IACrD,MAAM,UAAU,GACd,CAAA,aAAa,aAAb,aAAa,uBAAb,aAAa,CAAE,QAAQ,KAAI,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,GAAG,QAAQ,CAAC,CAAC,CAAC;IAC3E,MAAM,EAAE,WAAW,EAAE,aAAa,EAAE,GAAG,OAAO,CAAC;IAC/C,IAAI,EACF,GAAG,EAAE,KAAK,EACV,GAAG,EAAE,KAAK,EACV,IAAI,EAAE,QAAQ,GACf,GAAG,OAAO,CAAC,eAAe,IAAI,EAAE,CAAC;IAElC,MAAM,QAAQ,GAAG,IAAI,MAAM,CAAC,UAAU,EAAE,SAAS,CAAC,CAAC;IACnD,MAAM,QAAQ,GAAG,IAAI,MAAM,CAAC,UAAU,EAAE,SAAS,CAAC,CAAC;IACnD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,EAAE,CAAC,EAAE,EAAE;QAClC,MAAM,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,WAAW,CACtC;YACE,EAAE,EAAE,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC;YACvB,EAAE,EAAE,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC;SACxB,EACD,aAAa,CACd,CAAC;QAEF,MAAM,KAAK,GAAG,EAAE,CAAC;QACjB,MAAM,KAAK,GAAG,EAAE,CAAC;QACjB,IAAI,OAAO,CAAC,MAAM,KAAK,MAAM,EAAE;YAC7B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC,EAAE;gBACtC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBACnB,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;aACxB;SACF;aAAM;YACL,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;gBACnC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBACnB,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;aACrB;SACF;QAED,IAAI,WAAW,GAAG,OAAO,CACvB,WAAW,CACT;YACE,EAAE,EAAE,KAAK;YACT,EAAE,EAAE,KAAK;SACV,EACD,EAAE,kBAAkB,EAAE,CAAC,EAAE,QAAQ,EAAE,UAAU,EAAE,CAChD,EACD,EAAE,cAAc,EAAE,IAAI,EAAE,CACzB,CAAC;QACF,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,WAAW,CAAC;QAE7B,IAAI,CAAA,MAAA,OAAO,CAAC,eAAe,0CAAE,IAAI,MAAK,IAAI,EAAE;YAC1C,EAAE,GAAG,YAAY,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;YAC9B,EAAE,GAAG,IAAI,YAAY,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC;SAClC;QAED,IAAI,QAAQ,KAAK,IAAI,EAAE;YACrB,IAAI,CAAC,KAAK,CAAC,IAAI,KAAK,KAAK,SAAS,IAAI,KAAK,KAAK,SAAS,EAAE;gBACzD,0EAA0E;gBAC1E,MAAM,MAAM,GAAG,uBAAuB,CACpC,EAAE,EAAE,EAAE,EAAE,EAAE,EACV,EAAE,aAAa,EAAE,KAAK,EAAE,CACzB,CAAC;gBACF,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC;gBACnB,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC;aACpB;YACD,IAAI,KAAK,KAAK,SAAS,IAAI,KAAK,KAAK,SAAS,EAAE;gBAC9C,MAAM,CAAC,MAAM,CACX,WAAW,EACX,mBAAmB,CACjB,WAAW,EACX,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,EACvB,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CACxB,CACF,CAAC;aACH;SACF;QAED,QAAQ,CAAC,SAAS,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAC1B,QAAQ,CAAC,SAAS,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;KAC3B;IAED,OAAO;QACL,EAAE,EAAE,QAAQ;QACZ,EAAE,EAAE,QAAQ;KACb,CAAC;AACJ,CAAC"}
|
|
@@ -0,0 +1,15 @@
|
|
|
1
|
+
import { xMean } from 'ml-spectra-processing';
|
|
2
|
+
export function removeDCOffset(data, options) {
|
|
3
|
+
let { digitalFilterValue = 0 } = options;
|
|
4
|
+
const nbPoints = data.re.length;
|
|
5
|
+
const newRe = new Float64Array(data.re);
|
|
6
|
+
const newIm = new Float64Array(data.im);
|
|
7
|
+
const averageRe = xMean(data.re.slice((nbPoints * 0.75) >> 0, nbPoints - digitalFilterValue));
|
|
8
|
+
const averageIm = xMean(data.im.slice((nbPoints * 0.75) >> 0, nbPoints - digitalFilterValue));
|
|
9
|
+
for (let i = digitalFilterValue; i < nbPoints - digitalFilterValue; i++) {
|
|
10
|
+
newRe[i] -= averageRe;
|
|
11
|
+
newIm[i] -= averageIm;
|
|
12
|
+
}
|
|
13
|
+
return { re: newRe, im: newIm };
|
|
14
|
+
}
|
|
15
|
+
//# sourceMappingURL=removeDCOffset.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"removeDCOffset.js","sourceRoot":"","sources":["../../../../src/xyz/util/fft2d/removeDCOffset.ts"],"names":[],"mappings":"AAAA,OAAO,EAAY,KAAK,EAAE,MAAM,uBAAuB,CAAC;AAExD,MAAM,UAAU,cAAc,CAAC,IAAc,EAAE,OAAY;IACzD,IAAI,EAAE,kBAAkB,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;IACzC,MAAM,QAAQ,GAAG,IAAI,CAAC,EAAE,CAAC,MAAM,CAAC;IAChC,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IACxC,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IACxC,MAAM,SAAS,GAAG,KAAK,CACrB,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,QAAQ,GAAG,kBAAkB,CAAC,CACrE,CAAC;IACF,MAAM,SAAS,GAAG,KAAK,CACrB,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,QAAQ,GAAG,kBAAkB,CAAC,CACrE,CAAC;IACF,KAAK,IAAI,CAAC,GAAG,kBAAkB,EAAE,CAAC,GAAG,QAAQ,GAAG,kBAAkB,EAAE,CAAC,EAAE,EAAE;QACvE,KAAK,CAAC,CAAC,CAAC,IAAI,SAAS,CAAC;QACtB,KAAK,CAAC,CAAC,CAAC,IAAI,SAAS,CAAC;KACvB;IACD,OAAO,EAAE,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,EAAE,CAAC;AAClC,CAAC"}
|
|
@@ -0,0 +1,25 @@
|
|
|
1
|
+
export function zeroFilling(data, options) {
|
|
2
|
+
let { nbPoints, factor = 2, digitalFilterValue: grpdly = 0 } = options;
|
|
3
|
+
if (!nbPoints) {
|
|
4
|
+
nbPoints = 2 ** Math.round(Math.log2(data.re.length * factor));
|
|
5
|
+
}
|
|
6
|
+
let pointsToShift;
|
|
7
|
+
if (grpdly > 0) {
|
|
8
|
+
pointsToShift = Math.floor(grpdly);
|
|
9
|
+
}
|
|
10
|
+
else {
|
|
11
|
+
pointsToShift = 0;
|
|
12
|
+
}
|
|
13
|
+
const { re, im } = data;
|
|
14
|
+
let newRE = new Float64Array(nbPoints);
|
|
15
|
+
let newIM = new Float64Array(nbPoints);
|
|
16
|
+
const length = Math.min(nbPoints, re.length);
|
|
17
|
+
newRE.set(re.slice(0, length - pointsToShift));
|
|
18
|
+
newIM.set(im.slice(0, length - pointsToShift));
|
|
19
|
+
if (pointsToShift > 0 && pointsToShift < nbPoints) {
|
|
20
|
+
newRE.set(re.slice(re.length - pointsToShift), nbPoints - pointsToShift);
|
|
21
|
+
newIM.set(im.slice(re.length - pointsToShift), nbPoints - pointsToShift);
|
|
22
|
+
}
|
|
23
|
+
return { re: newRE, im: newIM };
|
|
24
|
+
}
|
|
25
|
+
//# sourceMappingURL=zeroFilling.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"zeroFilling.js","sourceRoot":"","sources":["../../../../src/xyz/util/fft2d/zeroFilling.ts"],"names":[],"mappings":"AAQA,MAAM,UAAU,WAAW,CAAC,IAAc,EAAE,OAA2B;IACrE,IAAI,EAAE,QAAQ,EAAE,MAAM,GAAG,CAAC,EAAE,kBAAkB,EAAE,MAAM,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;IAEvE,IAAI,CAAC,QAAQ,EAAE;QACb,QAAQ,GAAG,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC;KAChE;IAED,IAAI,aAAa,CAAC;IAClB,IAAI,MAAM,GAAG,CAAC,EAAE;QACd,aAAa,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;KACpC;SAAM;QACL,aAAa,GAAG,CAAC,CAAC;KACnB;IAED,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,IAAI,CAAC;IAExB,IAAI,KAAK,GAAG,IAAI,YAAY,CAAC,QAAQ,CAAC,CAAC;IACvC,IAAI,KAAK,GAAG,IAAI,YAAY,CAAC,QAAQ,CAAC,CAAC;IAEvC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,EAAE,CAAC,MAAM,CAAC,CAAC;IAE7C,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,MAAM,GAAG,aAAa,CAAC,CAAC,CAAC;IAC/C,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,MAAM,GAAG,aAAa,CAAC,CAAC,CAAC;IAE/C,IAAI,aAAa,GAAG,CAAC,IAAI,aAAa,GAAG,QAAQ,EAAE;QACjD,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,MAAM,GAAG,aAAa,CAAC,EAAE,QAAQ,GAAG,aAAa,CAAC,CAAC;QACzE,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,MAAM,GAAG,aAAa,CAAC,EAAE,QAAQ,GAAG,aAAa,CAAC,CAAC;KAC1E;IAED,OAAO,EAAE,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,EAAE,CAAC;AAClC,CAAC"}
|
|
@@ -1,4 +1,3 @@
|
|
|
1
|
-
import Matrix from 'ml-matrix';
|
|
2
1
|
export function padData(spectraData, options) {
|
|
3
2
|
let { minX, maxX, minY, maxY } = spectraData;
|
|
4
3
|
const width = options.width;
|
|
@@ -14,14 +13,18 @@ export function padData(spectraData, options) {
|
|
|
14
13
|
yDiff++;
|
|
15
14
|
let xOffset = xDiff / 2;
|
|
16
15
|
let yOffset = yDiff / 2;
|
|
17
|
-
|
|
16
|
+
const newMatrix = [];
|
|
17
|
+
for (let i = 0; i < nbSubSpectra + yDiff; i++) {
|
|
18
|
+
newMatrix.push(new Float64Array(nbPoints + xDiff));
|
|
19
|
+
}
|
|
18
20
|
for (let i = 0; i < nbSubSpectra; i++) {
|
|
19
21
|
for (let j = 0; j < nbPoints; j++) {
|
|
20
|
-
newMatrix
|
|
22
|
+
newMatrix[i + yOffset].set(spectraData.z[i], xOffset);
|
|
21
23
|
}
|
|
22
24
|
}
|
|
23
25
|
return {
|
|
24
|
-
|
|
26
|
+
...spectraData,
|
|
27
|
+
z: newMatrix,
|
|
25
28
|
minX: minX - xOffset * xInterval,
|
|
26
29
|
maxX: maxX + xOffset * xInterval,
|
|
27
30
|
minY: minY - yOffset * yInterval,
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"padData.js","sourceRoot":"","sources":["../../../src/xyz/util/padData.ts"],"names":[],"mappings":"
|
|
1
|
+
{"version":3,"file":"padData.js","sourceRoot":"","sources":["../../../src/xyz/util/padData.ts"],"names":[],"mappings":"AAEA,MAAM,UAAU,OAAO,CACrB,WAA0B,EAC1B,OAA0B;IAE1B,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,WAAW,CAAC;IAE7C,MAAM,KAAK,GAAG,OAAO,CAAC,KAAK,CAAC;IAE5B,IAAI,QAAQ,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;IACvC,IAAI,YAAY,GAAG,WAAW,CAAC,CAAC,CAAC,MAAM,CAAC;IAExC,IAAI,SAAS,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC;IACnD,IAAI,SAAS,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,CAAC,CAAC,CAAC;IAE/C,IAAI,KAAK,GAAG,KAAK,GAAG,YAAY,CAAC;IACjC,IAAI,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,QAAQ,EAAE,CAAC,CAAC,CAAC;IAC1C,IAAI,KAAK,GAAG,CAAC;QAAE,KAAK,EAAE,CAAC;IACvB,IAAI,KAAK,GAAG,CAAC;QAAE,KAAK,EAAE,CAAC;IAEvB,IAAI,OAAO,GAAG,KAAK,GAAG,CAAC,CAAC;IACxB,IAAI,OAAO,GAAG,KAAK,GAAG,CAAC,CAAC;IAExB,MAAM,SAAS,GAAmB,EAAE,CAAC;IACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,YAAY,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE;QAC7C,SAAS,CAAC,IAAI,CAAC,IAAI,YAAY,CAAC,QAAQ,GAAG,KAAK,CAAC,CAAC,CAAC;KACpD;IACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,YAAY,EAAE,CAAC,EAAE,EAAE;QACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE;YACjC,SAAS,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,GAAG,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;SACvD;KACF;IAED,OAAO;QACL,GAAG,WAAW;QACd,CAAC,EAAE,SAAS;QACZ,IAAI,EAAE,IAAI,GAAG,OAAO,GAAG,SAAS;QAChC,IAAI,EAAE,IAAI,GAAG,OAAO,GAAG,SAAS;QAChC,IAAI,EAAE,IAAI,GAAG,OAAO,GAAG,SAAS;QAChC,IAAI,EAAE,IAAI,GAAG,OAAO,GAAG,SAAS;KACjC,CAAC;AACJ,CAAC"}
|