nmr-processing 9.3.7 → 9.3.8-pre.1681135907
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/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/xyz/Data2D.d.ts +21 -0
- package/lib/xyz/Data2D.js +3 -0
- package/lib/xyz/Data2D.js.map +1 -0
- package/lib/xyz/quadrature.d.ts +8 -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/xyzBidimensionalFFT.d.ts +41 -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/utilities/checkData2DFid.js +6 -0
- package/lib-esm/utilities/checkData2DFid.js.map +1 -0
- package/lib-esm/xyz/Data2D.js +2 -0
- package/lib-esm/xyz/Data2D.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/xyzBidimensionalFFT.js +75 -0
- package/lib-esm/xyz/xyzBidimensionalFFT.js.map +1 -0
- package/package.json +17 -17
- 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/utilities/checkData2DFid.ts +10 -0
- package/src/xyz/Data2D.ts +25 -0
- package/src/xyz/quadrature.ts +69 -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/xyzBidimensionalFFT.ts +123 -0
|
@@ -0,0 +1,41 @@
|
|
|
1
|
+
import { DoubleArray } from 'cheminfo-types';
|
|
2
|
+
import { ApodizationOptions } from '../apodization/apodization';
|
|
3
|
+
import { Data2D } from './Data2D';
|
|
4
|
+
import { DirectDimensionOptions } from './util/fft2d/fftDirectDimension';
|
|
5
|
+
import { IndirectDimensionOptions } from './util/fft2d/fftIndirectDimension';
|
|
6
|
+
import { ZeroFillingOptions } from './util/fft2d/zeroFilling';
|
|
7
|
+
type Info2D = any;
|
|
8
|
+
interface BidimensionalFFTOptions {
|
|
9
|
+
fnMode: string;
|
|
10
|
+
info: Info2D;
|
|
11
|
+
zeroFilling?: ZeroFillingOptions;
|
|
12
|
+
phaseCorrection?: {
|
|
13
|
+
ph0: number;
|
|
14
|
+
ph1: number;
|
|
15
|
+
};
|
|
16
|
+
digitalFilterValue?: number;
|
|
17
|
+
apodization?: ApodizationOptions;
|
|
18
|
+
indirect?: IndirectDimensionOptions;
|
|
19
|
+
direct?: DirectDimensionOptions;
|
|
20
|
+
}
|
|
21
|
+
export declare function xyzBidimensionalFFT(data: Data2D, options: BidimensionalFFTOptions): {
|
|
22
|
+
re: {
|
|
23
|
+
minX: number;
|
|
24
|
+
maxX: number;
|
|
25
|
+
minY: number;
|
|
26
|
+
maxY: number;
|
|
27
|
+
z: DoubleArray[];
|
|
28
|
+
maxZ: number;
|
|
29
|
+
minZ: number;
|
|
30
|
+
};
|
|
31
|
+
im: {
|
|
32
|
+
minX: number;
|
|
33
|
+
maxX: number;
|
|
34
|
+
minY: number;
|
|
35
|
+
maxY: number;
|
|
36
|
+
z: DoubleArray[];
|
|
37
|
+
maxZ: number;
|
|
38
|
+
minZ: number;
|
|
39
|
+
};
|
|
40
|
+
};
|
|
41
|
+
export {};
|
|
@@ -0,0 +1,79 @@
|
|
|
1
|
+
"use strict";
|
|
2
|
+
Object.defineProperty(exports, "__esModule", { value: true });
|
|
3
|
+
exports.xyzBidimensionalFFT = void 0;
|
|
4
|
+
const quadrature_1 = require("./quadrature");
|
|
5
|
+
const fftDirectDimension_1 = require("./util/fft2d/fftDirectDimension");
|
|
6
|
+
const fftIndirectDimension_1 = require("./util/fft2d/fftIndirectDimension");
|
|
7
|
+
function xyzBidimensionalFFT(data, options) {
|
|
8
|
+
ensureData2DFid(data);
|
|
9
|
+
const { info, fnMode, digitalFilterValue = 0, zeroFilling: zeroFillingGeneral, } = options;
|
|
10
|
+
const orderedData = (0, quadrature_1.quadrature)(data, { fnMode });
|
|
11
|
+
const { re: reData } = orderedData;
|
|
12
|
+
const nbRows = reData.length;
|
|
13
|
+
const nbColumns = reData[0].length;
|
|
14
|
+
const { direct = {}, indirect = {} } = options;
|
|
15
|
+
const { zeroFilling: zeroFillingF2 = zeroFillingGeneral } = direct;
|
|
16
|
+
const { factor: f2Factor = 2 } = zeroFillingF2 || {};
|
|
17
|
+
const nbPointsF2 = (zeroFillingF2 === null || zeroFillingF2 === void 0 ? void 0 : zeroFillingF2.nbPoints) || 2 ** Math.round(Math.log2(nbColumns * f2Factor));
|
|
18
|
+
const transformed = (0, fftDirectDimension_1.fftDirectDimension)(orderedData, {
|
|
19
|
+
...direct,
|
|
20
|
+
fnMode,
|
|
21
|
+
digitalFilterValue,
|
|
22
|
+
zeroFilling: { nbPoints: nbPointsF2 },
|
|
23
|
+
});
|
|
24
|
+
const { zeroFilling: zeroFillingF1 } = options;
|
|
25
|
+
const { factor: f1Factor = 2 } = zeroFillingF1 || {};
|
|
26
|
+
const nbPointsF1 = (zeroFillingF1 === null || zeroFillingF1 === void 0 ? void 0 : zeroFillingF1.nbPoints) || 2 ** Math.round(Math.log2(nbRows * f1Factor));
|
|
27
|
+
const { re: reResult, im: imResult } = (0, fftIndirectDimension_1.fftIndirectDimension)(transformed, {
|
|
28
|
+
...indirect,
|
|
29
|
+
fnMode,
|
|
30
|
+
zeroFilling: { nbPoints: nbPointsF1 },
|
|
31
|
+
});
|
|
32
|
+
const minMaxXY = getMinMaxXY(info);
|
|
33
|
+
return {
|
|
34
|
+
re: {
|
|
35
|
+
z: convertToDoubleArray(reResult),
|
|
36
|
+
maxZ: reResult.max(),
|
|
37
|
+
minZ: reResult.min(),
|
|
38
|
+
...minMaxXY,
|
|
39
|
+
},
|
|
40
|
+
im: {
|
|
41
|
+
z: convertToDoubleArray(imResult),
|
|
42
|
+
maxZ: imResult.max(),
|
|
43
|
+
minZ: imResult.min(),
|
|
44
|
+
...minMaxXY,
|
|
45
|
+
},
|
|
46
|
+
};
|
|
47
|
+
}
|
|
48
|
+
exports.xyzBidimensionalFFT = xyzBidimensionalFFT;
|
|
49
|
+
function ensureData2DFid(data) {
|
|
50
|
+
if (!('re' in data)) {
|
|
51
|
+
throw new Error('bidimensional FFT only can be applied to raw data');
|
|
52
|
+
}
|
|
53
|
+
}
|
|
54
|
+
function convertToDoubleArray(data) {
|
|
55
|
+
const result = new Array(data.rows);
|
|
56
|
+
for (let i = 0; i < data.rows; i++) {
|
|
57
|
+
result[i] = new Float64Array(data.getRow(i));
|
|
58
|
+
}
|
|
59
|
+
return result;
|
|
60
|
+
}
|
|
61
|
+
function getMinMaxXY(info) {
|
|
62
|
+
const min = [];
|
|
63
|
+
const max = [];
|
|
64
|
+
for (let i = 0; i < info.baseFrequency.length; i++) {
|
|
65
|
+
const baseFrequency = Number.parseFloat(info.baseFrequency[i]);
|
|
66
|
+
const frequencyOffset = Number.parseFloat(info.frequencyOffset[i]);
|
|
67
|
+
const spectralWidth = Number.parseFloat(info.spectralWidth[i]);
|
|
68
|
+
const offset = frequencyOffset / baseFrequency;
|
|
69
|
+
min.push(offset - spectralWidth * 0.5);
|
|
70
|
+
max.push(offset + spectralWidth * 0.5);
|
|
71
|
+
}
|
|
72
|
+
return {
|
|
73
|
+
minX: min[0],
|
|
74
|
+
maxX: max[0],
|
|
75
|
+
minY: min[1],
|
|
76
|
+
maxY: max[1],
|
|
77
|
+
};
|
|
78
|
+
}
|
|
79
|
+
//# sourceMappingURL=xyzBidimensionalFFT.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"xyzBidimensionalFFT.js","sourceRoot":"","sources":["../../src/xyz/xyzBidimensionalFFT.ts"],"names":[],"mappings":";;;AAMA,6CAA0C;AAC1C,wEAGyC;AACzC,4EAG2C;AAgB3C,SAAgB,mBAAmB,CACjC,IAAY,EACZ,OAAgC;IAEhC,eAAe,CAAC,IAAI,CAAC,CAAC;IAEtB,MAAM,EACJ,IAAI,EACJ,MAAM,EACN,kBAAkB,GAAG,CAAC,EACtB,WAAW,EAAE,kBAAkB,GAChC,GAAG,OAAO,CAAC;IAEZ,MAAM,WAAW,GAAG,IAAA,uBAAU,EAAC,IAAI,EAAE,EAAE,MAAM,EAAE,CAAC,CAAC;IAEjD,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,WAAW,CAAC;IACnC,MAAM,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC;IAC7B,MAAM,SAAS,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;IACnC,MAAM,EAAE,MAAM,GAAG,EAAE,EAAE,QAAQ,GAAG,EAAE,EAAE,GAAG,OAAO,CAAC;IAE/C,MAAM,EAAE,WAAW,EAAE,aAAa,GAAG,kBAAkB,EAAE,GAAG,MAAM,CAAC;IAEnE,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,WAAW,GAAG,IAAA,uCAAkB,EAAC,WAAW,EAAE;QAClD,GAAG,MAAM;QACT,MAAM;QACN,kBAAkB;QAClB,WAAW,EAAE,EAAE,QAAQ,EAAE,UAAU,EAAE;KACtC,CAAC,CAAC;IAEH,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,EAAE,EAAE,QAAQ,EAAE,EAAE,EAAE,QAAQ,EAAE,GAAG,IAAA,2CAAoB,EAAC,WAAW,EAAE;QACvE,GAAG,QAAQ;QACX,MAAM;QACN,WAAW,EAAE,EAAE,QAAQ,EAAE,UAAU,EAAE;KACtC,CAAC,CAAC;IAEH,MAAM,QAAQ,GAAG,WAAW,CAAC,IAAI,CAAC,CAAC;IACnC,OAAO;QACL,EAAE,EAAE;YACF,CAAC,EAAE,oBAAoB,CAAC,QAAQ,CAAC;YACjC,IAAI,EAAE,QAAQ,CAAC,GAAG,EAAE;YACpB,IAAI,EAAE,QAAQ,CAAC,GAAG,EAAE;YACpB,GAAG,QAAQ;SACZ;QACD,EAAE,EAAE;YACF,CAAC,EAAE,oBAAoB,CAAC,QAAQ,CAAC;YACjC,IAAI,EAAE,QAAQ,CAAC,GAAG,EAAE;YACpB,IAAI,EAAE,QAAQ,CAAC,GAAG,EAAE;YACpB,GAAG,QAAQ;SACZ;KACF,CAAC;AACJ,CAAC;AA1DD,kDA0DC;AAED,SAAS,eAAe,CAAC,IAAY;IACnC,IAAI,CAAC,CAAC,IAAI,IAAI,IAAI,CAAC,EAAE;QACnB,MAAM,IAAI,KAAK,CAAC,mDAAmD,CAAC,CAAC;KACtE;AACH,CAAC;AAED,SAAS,oBAAoB,CAAC,IAAY;IACxC,MAAM,MAAM,GAAkB,IAAI,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IACnD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE;QAClC,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;KAC9C;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,SAAS,WAAW,CAAC,IAAY;IAC/B,MAAM,GAAG,GAAa,EAAE,CAAC;IACzB,MAAM,GAAG,GAAa,EAAE,CAAC;IACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QAClD,MAAM,aAAa,GAAG,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC,CAAC;QAC/D,MAAM,eAAe,GAAG,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,eAAe,CAAC,CAAC,CAAC,CAAC,CAAC;QACnE,MAAM,aAAa,GAAG,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC,CAAC;QAC/D,MAAM,MAAM,GAAG,eAAe,GAAG,aAAa,CAAC;QAC/C,GAAG,CAAC,IAAI,CAAC,MAAM,GAAG,aAAa,GAAG,GAAG,CAAC,CAAC;QACvC,GAAG,CAAC,IAAI,CAAC,MAAM,GAAG,aAAa,GAAG,GAAG,CAAC,CAAC;KACxC;IAED,OAAO;QACL,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC;QACZ,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC;QACZ,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC;QACZ,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC;KACb,CAAC;AACJ,CAAC"}
|
|
@@ -1,7 +1,9 @@
|
|
|
1
1
|
import { applyWindow } from './applyWindow';
|
|
2
2
|
import { compose } from './compose';
|
|
3
|
-
export function apodization(data, options) {
|
|
3
|
+
export function apodization(data, options = {}) {
|
|
4
4
|
const { compose: composeOptions, pointsToShift } = options;
|
|
5
|
+
if (!composeOptions)
|
|
6
|
+
return { ...data, windowData: [] };
|
|
5
7
|
const windowData = compose(composeOptions);
|
|
6
8
|
const applyWindowOptions = {
|
|
7
9
|
windowData,
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"apodization.js","sourceRoot":"","sources":["../../src/apodization/apodization.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,WAAW,EAAE,MAAM,eAAe,CAAC;AAC5C,OAAO,EAAE,OAAO,EAAgB,MAAM,WAAW,CAAC;AAiBlD,MAAM,UAAU,WAAW,CAAC,IAAc,EAAE,
|
|
1
|
+
{"version":3,"file":"apodization.js","sourceRoot":"","sources":["../../src/apodization/apodization.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,WAAW,EAAE,MAAM,eAAe,CAAC;AAC5C,OAAO,EAAE,OAAO,EAAgB,MAAM,WAAW,CAAC;AAiBlD,MAAM,UAAU,WAAW,CAAC,IAAc,EAAE,UAA8B,EAAE;IAC1E,MAAM,EAAE,OAAO,EAAE,cAAc,EAAE,aAAa,EAAE,GAAG,OAAO,CAAC;IAE3D,IAAI,CAAC,cAAc;QAAE,OAAO,EAAE,GAAG,IAAI,EAAE,UAAU,EAAE,EAAE,EAAE,CAAC;IAExD,MAAM,UAAU,GAAG,OAAO,CAAC,cAAc,CAAC,CAAC;IAE3C,MAAM,kBAAkB,GAAG;QACzB,UAAU;QACV,aAAa;KACd,CAAC;IACF,MAAM,EAAE,GAAG,WAAW,CAAC,IAAI,CAAC,EAAE,EAAE,kBAAkB,CAAC,CAAC;IACpD,MAAM,EAAE,GAAG,WAAW,CAAC,IAAI,CAAC,EAAE,EAAE,kBAAkB,CAAC,CAAC;IAEpD,OAAO,EAAE,EAAE,EAAE,EAAE,EAAE,UAAU,EAAE,CAAC;AAChC,CAAC"}
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"applyWindow.js","sourceRoot":"","sources":["../../src/apodization/applyWindow.ts"],"names":[],"mappings":"AA0BA;;GAEG;AACH,MAAM,UAAU,WAAW,CAAC,IAAiB,EAAE,OAA2B;IACxE,MAAM,UAAU,GAAG,IAAI,CAAC,MAAM,CAAC;IAE/B,MAAM,EACJ,UAAU,EACV,KAAK,GAAG,CAAC,EACT,MAAM,GAAG,UAAU,EACnB,aAAa,GAAG,CAAC,EACjB,MAAM,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,GAChC,GAAG,OAAO,CAAC;IAEZ,MAAM,aAAa,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,MAAM,EAAE,UAAU,GAAG,aAAa,CAAC,CAAC;IAC3E,KAAK,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,aAAa,EAAE,CAAC,EAAE,EAAE;QACjD,MAAM,CAAC,CAAC,CAAC,IAAI,UAAU,CAAC,CAAC,EAAE,CAAC,CAAC;KAC9B;
|
|
1
|
+
{"version":3,"file":"applyWindow.js","sourceRoot":"","sources":["../../src/apodization/applyWindow.ts"],"names":[],"mappings":"AA0BA;;GAEG;AACH,MAAM,UAAU,WAAW,CAAC,IAAiB,EAAE,OAA2B;IACxE,MAAM,UAAU,GAAG,IAAI,CAAC,MAAM,CAAC;IAE/B,MAAM,EACJ,UAAU,EACV,KAAK,GAAG,CAAC,EACT,MAAM,GAAG,UAAU,EACnB,aAAa,GAAG,CAAC,EACjB,MAAM,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,GAChC,GAAG,OAAO,CAAC;IAEZ,MAAM,aAAa,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,MAAM,EAAE,UAAU,GAAG,aAAa,CAAC,CAAC;IAC3E,KAAK,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,aAAa,EAAE,CAAC,EAAE,EAAE;QACjD,MAAM,CAAC,CAAC,CAAC,IAAI,UAAU,CAAC,CAAC,EAAE,CAAC,CAAC;KAC9B;IACD,KAAK,IAAI,CAAC,GAAG,UAAU,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,GAAG,aAAa,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QAC3E,MAAM,CAAC,CAAC,CAAC,IAAI,UAAU,CAAC,CAAC,EAAE,CAAC,CAAC;KAC9B;IAED,OAAO,MAAM,CAAC;AAChB,CAAC"}
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"compose.js","sourceRoot":"","sources":["../../src/apodization/compose.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,WAAW,EAAE,MAAM,eAAe,CAAC;AAC5C,OAAO,EAAE,WAAW,EAAE,MAAM,eAAe,CAAC;AAE5C,OAAO,EAAE,OAAO,EAAE,MAAM,iBAAiB,CAAC;AAyB1C,MAAM,UAAU,OAAO,CAAC,YAA0B;IAChD,MAAM,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,GAAG,EAAE,EAAE,GAAG,YAAY,CAAC;IACzD,IAAI,IAAI,GAAG,IAAI,YAAY,CAAC,UAAU,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAChD,MAAM,CAAC,OAAO,CAAC,CAAC,OAAO,EAAE,EAAE;QACzB,MAAM,EAAE,KAAK,EAAE,KAAK,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;QACrC,MAAM,EAAE,OAAO,EAAE,YAAY,EAAE,GAAG,KAAK,CAAC;QACxC,MAAM,MAAM,GACV,QAAQ,IAAI,YAAY,CAAC,CAAC,CAAC,YAAY,CAAC,MAAM,CAAC,CAAC,CAAC,UAAU,GAAG,KAAK,CAAC;QAEtE,MAAM,UAAU,GAAG,OAAO,CAAC,WAAW,CAAC,KAAK,CAAC,EAAE,MAAM,CAAC,CAAC;
|
|
1
|
+
{"version":3,"file":"compose.js","sourceRoot":"","sources":["../../src/apodization/compose.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,WAAW,EAAE,MAAM,eAAe,CAAC;AAC5C,OAAO,EAAE,WAAW,EAAE,MAAM,eAAe,CAAC;AAE5C,OAAO,EAAE,OAAO,EAAE,MAAM,iBAAiB,CAAC;AAyB1C,MAAM,UAAU,OAAO,CAAC,YAA0B;IAChD,MAAM,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,GAAG,EAAE,EAAE,GAAG,YAAY,CAAC;IACzD,IAAI,IAAI,GAAG,IAAI,YAAY,CAAC,UAAU,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAChD,MAAM,CAAC,OAAO,CAAC,CAAC,OAAO,EAAE,EAAE;QACzB,MAAM,EAAE,KAAK,EAAE,KAAK,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;QACrC,MAAM,EAAE,OAAO,EAAE,YAAY,EAAE,GAAG,KAAK,CAAC;QACxC,MAAM,MAAM,GACV,QAAQ,IAAI,YAAY,CAAC,CAAC,CAAC,YAAY,CAAC,MAAM,CAAC,CAAC,CAAC,UAAU,GAAG,KAAK,CAAC;QAEtE,MAAM,UAAU,GAAG,OAAO,CAAC,WAAW,CAAC,KAAK,CAAC,EAAE,MAAM,CAAC,CAAC;QAEvD,WAAW,CAAC,IAAI,EAAE;YAChB,UAAU;YACV,MAAM;YACN,KAAK;YACL,MAAM,EAAE,IAAI;SACb,CAAC,CAAC;IACL,CAAC,CAAC,CAAC;IACH,OAAO,IAAI,CAAC;AACd,CAAC"}
|
|
@@ -1,5 +1,6 @@
|
|
|
1
1
|
import { exponential } from './shapes/exponential';
|
|
2
2
|
import { lorentzToGauss } from './shapes/lorentzToGauss';
|
|
3
|
+
import { sineBell } from './shapes/sineBell';
|
|
3
4
|
export function getFunction(shape) {
|
|
4
5
|
const { kind, options } = shape;
|
|
5
6
|
switch (kind) {
|
|
@@ -7,6 +8,8 @@ export function getFunction(shape) {
|
|
|
7
8
|
return exponential(options);
|
|
8
9
|
case 'lorentzToGauss':
|
|
9
10
|
return lorentzToGauss(options);
|
|
11
|
+
case 'sineBell':
|
|
12
|
+
return sineBell(options);
|
|
10
13
|
default:
|
|
11
14
|
throw Error(`Unknown distribution ${kind}`);
|
|
12
15
|
}
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"getFunction.js","sourceRoot":"","sources":["../../src/apodization/getFunction.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,WAAW,EAAE,MAAM,sBAAsB,CAAC;AACnD,OAAO,EAAE,cAAc,EAAE,MAAM,yBAAyB,CAAC;
|
|
1
|
+
{"version":3,"file":"getFunction.js","sourceRoot":"","sources":["../../src/apodization/getFunction.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,WAAW,EAAE,MAAM,sBAAsB,CAAC;AACnD,OAAO,EAAE,cAAc,EAAE,MAAM,yBAAyB,CAAC;AACzD,OAAO,EAAE,QAAQ,EAAE,MAAM,mBAAmB,CAAC;AAE7C,MAAM,UAAU,WAAW,CAAC,KAAsB;IAChD,MAAM,EAAE,IAAI,EAAE,OAAO,EAAE,GAAG,KAAK,CAAC;IAChC,QAAQ,IAAI,EAAE;QACZ,KAAK,aAAa;YAChB,OAAO,WAAW,CAAC,OAAO,CAAC,CAAC;QAC9B,KAAK,gBAAgB;YACnB,OAAO,cAAc,CAAC,OAAO,CAAC,CAAC;QACjC,KAAK,UAAU;YACb,OAAO,QAAQ,CAAC,OAAO,CAAC,CAAC;QAC3B;YACE,MAAM,KAAK,CAAC,wBAAwB,IAAc,EAAE,CAAC,CAAC;KACzD;AACH,CAAC"}
|
|
@@ -0,0 +1,13 @@
|
|
|
1
|
+
export function sineBell(options) {
|
|
2
|
+
const { exponent = 1, offset = 0, end = 1, length } = options;
|
|
3
|
+
if (length === undefined) {
|
|
4
|
+
throw new Error('length options should exists for sineBell shape');
|
|
5
|
+
}
|
|
6
|
+
if (end < 0 || end > 1) {
|
|
7
|
+
throw new Error(`the end parameter should be [0-1]`);
|
|
8
|
+
}
|
|
9
|
+
const c1 = Math.PI * offset;
|
|
10
|
+
const c2 = (Math.PI * (end - offset)) / (length - 1);
|
|
11
|
+
return (i) => Math.pow(Math.sin(c1 + c2 * i), exponent);
|
|
12
|
+
}
|
|
13
|
+
//# sourceMappingURL=sineBell.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"sineBell.js","sourceRoot":"","sources":["../../../src/apodization/shapes/sineBell.ts"],"names":[],"mappings":"AA2BA,MAAM,UAAU,QAAQ,CAAC,OAAwB;IAC/C,MAAM,EAAE,QAAQ,GAAG,CAAC,EAAE,MAAM,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,MAAM,EAAE,GAAG,OAAO,CAAC;IAE9D,IAAI,MAAM,KAAK,SAAS,EAAE;QACxB,MAAM,IAAI,KAAK,CAAC,iDAAiD,CAAC,CAAC;KACpE;IAED,IAAI,GAAG,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,EAAE;QACtB,MAAM,IAAI,KAAK,CAAC,mCAAmC,CAAC,CAAC;KACtD;IACD,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,CAAC;IAC5B,MAAM,EAAE,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,MAAM,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IACrD,OAAO,CAAC,CAAS,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC;AAClE,CAAC"}
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"getData.js","sourceRoot":"","sources":["../../../src/apodization/utils/getData.ts"],"names":[],"mappings":"AAAA;;GAEG;AACH,MAAM,UAAU,OAAO,CAAC,IAA2B,EAAE,MAAc;IACjE,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,CAAC;IACtC,IAAI,GAAG,GAAG,MAAM,CAAC,gBAAgB,CAAC;IAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE;QAC/B,MAAM,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QACtB,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;QAChB,IAAI,KAAK,GAAG,GAAG;YAAE,GAAG,GAAG,KAAK,CAAC;KAC9B;
|
|
1
|
+
{"version":3,"file":"getData.js","sourceRoot":"","sources":["../../../src/apodization/utils/getData.ts"],"names":[],"mappings":"AAAA;;GAEG;AACH,MAAM,UAAU,OAAO,CAAC,IAA2B,EAAE,MAAc;IACjE,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,CAAC;IACtC,IAAI,GAAG,GAAG,MAAM,CAAC,gBAAgB,CAAC;IAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE;QAC/B,MAAM,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QACtB,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;QAChB,IAAI,KAAK,GAAG,GAAG;YAAE,GAAG,GAAG,KAAK,CAAC;KAC9B;IACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE;QAAE,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC;IAChD,OAAO,IAAI,CAAC;AACd,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":"Data2D.js","sourceRoot":"","sources":["../../src/xyz/Data2D.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;AAGhF,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"}
|
|
@@ -0,0 +1,75 @@
|
|
|
1
|
+
import { quadrature } from './quadrature';
|
|
2
|
+
import { fftDirectDimension, } from './util/fft2d/fftDirectDimension';
|
|
3
|
+
import { fftIndirectDimension, } from './util/fft2d/fftIndirectDimension';
|
|
4
|
+
export function xyzBidimensionalFFT(data, options) {
|
|
5
|
+
ensureData2DFid(data);
|
|
6
|
+
const { info, fnMode, digitalFilterValue = 0, zeroFilling: zeroFillingGeneral, } = options;
|
|
7
|
+
const orderedData = quadrature(data, { fnMode });
|
|
8
|
+
const { re: reData } = orderedData;
|
|
9
|
+
const nbRows = reData.length;
|
|
10
|
+
const nbColumns = reData[0].length;
|
|
11
|
+
const { direct = {}, indirect = {} } = options;
|
|
12
|
+
const { zeroFilling: zeroFillingF2 = zeroFillingGeneral } = direct;
|
|
13
|
+
const { factor: f2Factor = 2 } = zeroFillingF2 || {};
|
|
14
|
+
const nbPointsF2 = (zeroFillingF2 === null || zeroFillingF2 === void 0 ? void 0 : zeroFillingF2.nbPoints) || 2 ** Math.round(Math.log2(nbColumns * f2Factor));
|
|
15
|
+
const transformed = fftDirectDimension(orderedData, {
|
|
16
|
+
...direct,
|
|
17
|
+
fnMode,
|
|
18
|
+
digitalFilterValue,
|
|
19
|
+
zeroFilling: { nbPoints: nbPointsF2 },
|
|
20
|
+
});
|
|
21
|
+
const { zeroFilling: zeroFillingF1 } = options;
|
|
22
|
+
const { factor: f1Factor = 2 } = zeroFillingF1 || {};
|
|
23
|
+
const nbPointsF1 = (zeroFillingF1 === null || zeroFillingF1 === void 0 ? void 0 : zeroFillingF1.nbPoints) || 2 ** Math.round(Math.log2(nbRows * f1Factor));
|
|
24
|
+
const { re: reResult, im: imResult } = fftIndirectDimension(transformed, {
|
|
25
|
+
...indirect,
|
|
26
|
+
fnMode,
|
|
27
|
+
zeroFilling: { nbPoints: nbPointsF1 },
|
|
28
|
+
});
|
|
29
|
+
const minMaxXY = getMinMaxXY(info);
|
|
30
|
+
return {
|
|
31
|
+
re: {
|
|
32
|
+
z: convertToDoubleArray(reResult),
|
|
33
|
+
maxZ: reResult.max(),
|
|
34
|
+
minZ: reResult.min(),
|
|
35
|
+
...minMaxXY,
|
|
36
|
+
},
|
|
37
|
+
im: {
|
|
38
|
+
z: convertToDoubleArray(imResult),
|
|
39
|
+
maxZ: imResult.max(),
|
|
40
|
+
minZ: imResult.min(),
|
|
41
|
+
...minMaxXY,
|
|
42
|
+
},
|
|
43
|
+
};
|
|
44
|
+
}
|
|
45
|
+
function ensureData2DFid(data) {
|
|
46
|
+
if (!('re' in data)) {
|
|
47
|
+
throw new Error('bidimensional FFT only can be applied to raw data');
|
|
48
|
+
}
|
|
49
|
+
}
|
|
50
|
+
function convertToDoubleArray(data) {
|
|
51
|
+
const result = new Array(data.rows);
|
|
52
|
+
for (let i = 0; i < data.rows; i++) {
|
|
53
|
+
result[i] = new Float64Array(data.getRow(i));
|
|
54
|
+
}
|
|
55
|
+
return result;
|
|
56
|
+
}
|
|
57
|
+
function getMinMaxXY(info) {
|
|
58
|
+
const min = [];
|
|
59
|
+
const max = [];
|
|
60
|
+
for (let i = 0; i < info.baseFrequency.length; i++) {
|
|
61
|
+
const baseFrequency = Number.parseFloat(info.baseFrequency[i]);
|
|
62
|
+
const frequencyOffset = Number.parseFloat(info.frequencyOffset[i]);
|
|
63
|
+
const spectralWidth = Number.parseFloat(info.spectralWidth[i]);
|
|
64
|
+
const offset = frequencyOffset / baseFrequency;
|
|
65
|
+
min.push(offset - spectralWidth * 0.5);
|
|
66
|
+
max.push(offset + spectralWidth * 0.5);
|
|
67
|
+
}
|
|
68
|
+
return {
|
|
69
|
+
minX: min[0],
|
|
70
|
+
maxX: max[0],
|
|
71
|
+
minY: min[1],
|
|
72
|
+
maxY: max[1],
|
|
73
|
+
};
|
|
74
|
+
}
|
|
75
|
+
//# sourceMappingURL=xyzBidimensionalFFT.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"xyzBidimensionalFFT.js","sourceRoot":"","sources":["../../src/xyz/xyzBidimensionalFFT.ts"],"names":[],"mappings":"AAMA,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAC1C,OAAO,EAEL,kBAAkB,GACnB,MAAM,iCAAiC,CAAC;AACzC,OAAO,EACL,oBAAoB,GAErB,MAAM,mCAAmC,CAAC;AAgB3C,MAAM,UAAU,mBAAmB,CACjC,IAAY,EACZ,OAAgC;IAEhC,eAAe,CAAC,IAAI,CAAC,CAAC;IAEtB,MAAM,EACJ,IAAI,EACJ,MAAM,EACN,kBAAkB,GAAG,CAAC,EACtB,WAAW,EAAE,kBAAkB,GAChC,GAAG,OAAO,CAAC;IAEZ,MAAM,WAAW,GAAG,UAAU,CAAC,IAAI,EAAE,EAAE,MAAM,EAAE,CAAC,CAAC;IAEjD,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,WAAW,CAAC;IACnC,MAAM,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC;IAC7B,MAAM,SAAS,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;IACnC,MAAM,EAAE,MAAM,GAAG,EAAE,EAAE,QAAQ,GAAG,EAAE,EAAE,GAAG,OAAO,CAAC;IAE/C,MAAM,EAAE,WAAW,EAAE,aAAa,GAAG,kBAAkB,EAAE,GAAG,MAAM,CAAC;IAEnE,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,WAAW,GAAG,kBAAkB,CAAC,WAAW,EAAE;QAClD,GAAG,MAAM;QACT,MAAM;QACN,kBAAkB;QAClB,WAAW,EAAE,EAAE,QAAQ,EAAE,UAAU,EAAE;KACtC,CAAC,CAAC;IAEH,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,EAAE,EAAE,QAAQ,EAAE,EAAE,EAAE,QAAQ,EAAE,GAAG,oBAAoB,CAAC,WAAW,EAAE;QACvE,GAAG,QAAQ;QACX,MAAM;QACN,WAAW,EAAE,EAAE,QAAQ,EAAE,UAAU,EAAE;KACtC,CAAC,CAAC;IAEH,MAAM,QAAQ,GAAG,WAAW,CAAC,IAAI,CAAC,CAAC;IACnC,OAAO;QACL,EAAE,EAAE;YACF,CAAC,EAAE,oBAAoB,CAAC,QAAQ,CAAC;YACjC,IAAI,EAAE,QAAQ,CAAC,GAAG,EAAE;YACpB,IAAI,EAAE,QAAQ,CAAC,GAAG,EAAE;YACpB,GAAG,QAAQ;SACZ;QACD,EAAE,EAAE;YACF,CAAC,EAAE,oBAAoB,CAAC,QAAQ,CAAC;YACjC,IAAI,EAAE,QAAQ,CAAC,GAAG,EAAE;YACpB,IAAI,EAAE,QAAQ,CAAC,GAAG,EAAE;YACpB,GAAG,QAAQ;SACZ;KACF,CAAC;AACJ,CAAC;AAED,SAAS,eAAe,CAAC,IAAY;IACnC,IAAI,CAAC,CAAC,IAAI,IAAI,IAAI,CAAC,EAAE;QACnB,MAAM,IAAI,KAAK,CAAC,mDAAmD,CAAC,CAAC;KACtE;AACH,CAAC;AAED,SAAS,oBAAoB,CAAC,IAAY;IACxC,MAAM,MAAM,GAAkB,IAAI,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IACnD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE;QAClC,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;KAC9C;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,SAAS,WAAW,CAAC,IAAY;IAC/B,MAAM,GAAG,GAAa,EAAE,CAAC;IACzB,MAAM,GAAG,GAAa,EAAE,CAAC;IACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,aAAa,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QAClD,MAAM,aAAa,GAAG,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC,CAAC;QAC/D,MAAM,eAAe,GAAG,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,eAAe,CAAC,CAAC,CAAC,CAAC,CAAC;QACnE,MAAM,aAAa,GAAG,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC,CAAC;QAC/D,MAAM,MAAM,GAAG,eAAe,GAAG,aAAa,CAAC;QAC/C,GAAG,CAAC,IAAI,CAAC,MAAM,GAAG,aAAa,GAAG,GAAG,CAAC,CAAC;QACvC,GAAG,CAAC,IAAI,CAAC,MAAM,GAAG,aAAa,GAAG,GAAG,CAAC,CAAC;KACxC;IAED,OAAO;QACL,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC;QACZ,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC;QACZ,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC;QACZ,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC;KACb,CAAC;AACJ,CAAC"}
|