nmr-processing 9.7.8 → 9.8.0

package/lib/xyz/util/fft2d/fftIndirectDimension.d.ts

@@ -9,6 +9,7 @@ interface PhaseCorrection {
 }
 export interface IndirectDimensionOptions {
     fnMode?: string;
+    reverse?: boolean;
     zeroFilling?: ZeroFillingOptions;
     apodization?: ApodizationFilterOptions;
     phaseCorrection?: PhaseCorrection;
@@ -16,8 +17,5 @@ export interface IndirectDimensionOptions {
 export declare function fftIndirectDimension(data: {
     re: DoubleArray[] | Matrix;
     im: DoubleArray[] | Matrix;
-}, options: IndirectDimensionOptions): {
-    re: Matrix;
-    im: Matrix;
-};
+}, options: IndirectDimensionOptions): Record<string, Matrix>;
 export {};

package/lib/xyz/util/fft2d/fftIndirectDimension.js

@@ -9,63 +9,117 @@ const ml_spectra_processing_1 = require("ml-spectra-processing");
 const apodization_1 = require("../../../apodization/apodization");
 const zeroFilling_1 = require("./zeroFilling");
 function fftIndirectDimension(data, options) {
-    const reData = ml_matrix_1.default.checkMatrix(data.re);
-    const imData = ml_matrix_1.default.checkMatrix(data.im);
-    const nbRows = reData.rows;
-    const nbColumns = reData.columns;
-    const { zeroFilling: zeroFillingF1 } = options;
+    const newData = {
+        re: ml_matrix_1.default.checkMatrix(data.re),
+        im: ml_matrix_1.default.checkMatrix(data.im),
+    };
+    const nbRows = newData.re.rows;
+    const nbColumns = newData.re.columns;
+    const { zeroFilling: zeroFillingF1, reverse } = options;
     const { factor: f1Factor = 2 } = zeroFillingF1 || {};
-    const nbPointsF1 = zeroFillingF1?.nbPoints || 2 ** Math.round(Math.log2(nbRows * f1Factor));
-    const { apodization: apodizationF1 } = options;
-    let { ph0: ph0F1, ph1: ph1F1, mode: pcModeF1, } = options.phaseCorrection || {};
-    const reResult = new ml_matrix_1.default(nbPointsF1, nbColumns);
-    const imResult = new ml_matrix_1.default(nbPointsF1, nbColumns);
-    for (let i = 0; i < nbColumns; i++) {
-        const { re: aRe, im: aIm } = (0, apodization_1.apodization)({
-            re: reData.getColumn(i),
-            im: imData.getColumn(i),
-        }, apodizationF1);
-        const newRe = [];
-        const newIm = [];
-        if (options.fnMode === 'echo') {
-            for (let j = 0; j < aRe.length; j += 2) {
-                newRe.push(aRe[j]);
-                newIm.push(aRe[j + 1]);
-            }
+    const nbPoints = zeroFillingF1?.nbPoints || 2 ** Math.round(Math.log2(nbRows * f1Factor));
+    const { fnMode } = options;
+    const newOptions = {
+        ...options,
+        zeroFilling: { ...zeroFillingF1, nbPoints },
+    };
+    const quadrants = {};
+    if (fnMode === 'QF') {
+        const reResult = new ml_matrix_1.default(nbPoints, nbColumns);
+        const imResult = new ml_matrix_1.default(nbPoints, nbColumns);
+        for (let i = 0; i < nbColumns; i++) {
+            const data = {
+                re: newData.re.getColumn(i),
+                im: reverse
+                    ? newData.im.getColumn(i)
+                    : (0, ml_spectra_processing_1.xMultiply)(newData.im.getColumn(i), -1),
+            };
+            const { re, im } = applyFFT(data, newOptions);
+            reResult.setColumn(i, re);
+            imResult.setColumn(i, im);
         }
-        else {
-            for (let j = 0; j < aRe.length; j++) {
-                newRe.push(aRe[j]);
-                newIm.push(-aIm[j]);
+        quadrants.rr = reResult;
+        quadrants.ir = imResult;
+    }
+    else {
+        // four quadrant processing
+        for (const current in newData) {
+            const reResult = new ml_matrix_1.default(nbPoints, nbColumns);
+            const imResult = new ml_matrix_1.default(nbPoints, nbColumns);
+            for (let i = 0; i < nbColumns; i++) {
+                const data = getColumnFrom(newData[current], { fnMode, index: i });
+                const { re, im } = applyFFT(data, newOptions);
+                reResult.setColumn(i, re);
+                imResult.setColumn(i, im);
             }
-        }
-        let transformed = (0, ml_spectra_processing_1.reimFFT)((0, zeroFilling_1.zeroFilling)({
-            re: newRe,
-            im: newIm,
-        }, { digitalFilterValue: 0, nbPoints: nbPointsF1 }), { applyZeroShift: true });
-        let { re, im } = transformed;
-        if (options.phaseCorrection?.mode === 'mc') {
-            re = (0, ml_spectra_processing_1.reimAbsolute)({ re, im });
-            im = new Float64Array(im.length);
-        }
-        if (pcModeF1 === 'pk') {
-            if (i === 0 && ph0F1 === undefined && ph1F1 === undefined) {
-                //here we can improve it by checking if the experiment is phase sensitive.
-                const phased = (0, ml_spectra_processing_1.reimAutoPhaseCorrection)({ re, im }, { magnitudeMode: false });
-                ph0F1 = phased.ph0;
-                ph1F1 = phased.ph1;
+            if (current === 're') {
+                quadrants.rr = reResult;
+                quadrants.ir = imResult;
             }
-            if (ph0F1 !== undefined && ph1F1 !== undefined) {
-                Object.assign(transformed, (0, ml_spectra_processing_1.reimPhaseCorrection)(transformed, (ph0F1 * Math.PI) / 180, (ph1F1 * Math.PI) / 180));
+            else {
+                quadrants.ri = reResult;
+                quadrants.ii = imResult;
             }
         }
-        reResult.setColumn(i, re);
-        imResult.setColumn(i, im);
     }
-    return {
-        re: reResult,
-        im: imResult,
-    };
+    return quadrants;
 }
 exports.fftIndirectDimension = fftIndirectDimension;
+function applyFFT(data, options) {
+    const { fnMode, apodization: apodizationF1, zeroFilling: zeroFillingF1 = {}, phaseCorrection = {}, } = options;
+    const { re, im } = (0, apodization_1.apodization)(data, apodizationF1);
+    const { nbPoints = 0 } = zeroFillingF1;
+    const transformed = (0, ml_spectra_processing_1.reimFFT)((0, zeroFilling_1.zeroFilling)({ re, im }, { nbPoints }), {
+        applyZeroShift: !['SEQ', 'TPPI'].some((k) => fnMode === k),
+    });
+    if (options.phaseCorrection?.mode === 'mc') {
+        Object.assign(transformed, {
+            re: (0, ml_spectra_processing_1.reimAbsolute)(transformed),
+            im: new Float64Array(nbPoints),
+        });
+    }
+    const { mode: pcModeF1, ph0: ph0F1, ph1: ph1F1 } = phaseCorrection;
+    if (pcModeF1 === 'pk' && ph0F1 !== undefined && ph1F1 !== undefined) {
+        Object.assign(transformed, (0, ml_spectra_processing_1.reimPhaseCorrection)(transformed, (ph0F1 * Math.PI) / 180, (ph1F1 * Math.PI) / 180));
+    }
+    return transformed;
+}
+function getColumnFrom(data, info) {
+    const { fnMode, index } = info;
+    const reverse = info.reverse || fnMode === 'QF';
+    const nbPoints = data.rows / 2;
+    const re = new Float64Array(nbPoints);
+    const im = new Float64Array(nbPoints);
+    switch (fnMode?.toLocaleLowerCase()) {
+        // case 'states':
+        // case 'echo-antiecho':
+        case 'seq':
+        case 'tppi':
+            return {
+                re: data.getColumn(index),
+                im: new Float64Array(data.rows),
+            };
+        default:
+            if (reverse) {
+                for (let i = 0; i < nbPoints; i++) {
+                    re[i] = data.get(2 * i, index);
+                    im[i] = -data.get(2 * i + 1, index);
+                }
+            }
+            else {
+                for (let i = 0; i < nbPoints; i++) {
+                    re[i] = data.get(2 * i, index);
+                    im[i] = data.get(2 * i + 1, index);
+                }
+            }
+    }
+    if (fnMode === 'states-ttpi') {
+        //it is to convert to the equivalent of states
+        for (let i = 0; i < nbPoints; i += 2) {
+            re[i] *= -1;
+            im[i] *= -1;
+        }
+    }
+    return { re, im };
+}
 //# sourceMappingURL=fftIndirectDimension.js.map

package/lib/xyz/util/fft2d/fftIndirectDimension.js.map

@@ -1 +1 @@
-{"version":3,"file":"fftIndirectDimension.js","sourceRoot":"","sources":["../../../../src/xyz/util/fft2d/fftIndirectDimension.ts"],"names":[],"mappings":";;;;;;AACA,0DAA+B;AAC/B,iEAK+B;AAE/B,kEAG0C;AAE1C,+CAAgE;AAehE,SAAgB,oBAAoB,CAClC,IAAgE,EAChE,OAAiC;IAEjC,MAAM,MAAM,GAAG,mBAAM,CAAC,WAAW,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAC3C,MAAM,MAAM,GAAG,mBAAM,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,aAAa,EAAE,QAAQ,IAAI,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,mBAAM,CAAC,UAAU,EAAE,SAAS,CAAC,CAAC;IACnD,MAAM,QAAQ,GAAG,IAAI,mBAAM,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,IAAA,yBAAW,EACtC;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,IAAA,+BAAO,EACvB,IAAA,yBAAW,EACT;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,OAAO,CAAC,eAAe,EAAE,IAAI,KAAK,IAAI,EAAE;YAC1C,EAAE,GAAG,IAAA,oCAAY,EAAC,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,IAAA,+CAAuB,EACpC,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,IAAA,2CAAmB,EACjB,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;AA7FD,oDA6FC"}
+{"version":3,"file":"fftIndirectDimension.js","sourceRoot":"","sources":["../../../../src/xyz/util/fft2d/fftIndirectDimension.ts"],"names":[],"mappings":";;;;;;AACA,0DAA+B;AAC/B,iEAM+B;AAE/B,kEAG0C;AAE1C,+CAAgE;AAgBhE,SAAgB,oBAAoB,CAClC,IAAgE,EAChE,OAAiC;IAEjC,MAAM,OAAO,GAA2B;QACtC,EAAE,EAAE,mBAAM,CAAC,WAAW,CAAC,IAAI,CAAC,EAAE,CAAC;QAC/B,EAAE,EAAE,mBAAM,CAAC,WAAW,CAAC,IAAI,CAAC,EAAE,CAAC;KAChC,CAAC;IAEF,MAAM,MAAM,GAAG,OAAO,CAAC,EAAE,CAAC,IAAI,CAAC;IAC/B,MAAM,SAAS,GAAG,OAAO,CAAC,EAAE,CAAC,OAAO,CAAC;IAErC,MAAM,EAAE,WAAW,EAAE,aAAa,EAAE,OAAO,EAAE,GAAG,OAAO,CAAC;IACxD,MAAM,EAAE,MAAM,EAAE,QAAQ,GAAG,CAAC,EAAE,GAAG,aAAa,IAAI,EAAE,CAAC;IAErD,MAAM,QAAQ,GACZ,aAAa,EAAE,QAAQ,IAAI,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,GAAG,QAAQ,CAAC,CAAC,CAAC;IAC3E,MAAM,EAAE,MAAM,EAAE,GAAG,OAAO,CAAC;IAE3B,MAAM,UAAU,GAA6B;QAC3C,GAAG,OAAO;QACV,WAAW,EAAE,EAAE,GAAG,aAAa,EAAE,QAAQ,EAAE;KAC5C,CAAC;IACF,MAAM,SAAS,GAA2B,EAAE,CAAC;IAC7C,IAAI,MAAM,KAAK,IAAI,EAAE;QACnB,MAAM,QAAQ,GAAG,IAAI,mBAAM,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;QACjD,MAAM,QAAQ,GAAG,IAAI,mBAAM,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;QACjD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,EAAE,CAAC,EAAE,EAAE;YAClC,MAAM,IAAI,GAAG;gBACX,EAAE,EAAE,OAAO,CAAC,EAAE,CAAC,SAAS,CAAC,CAAC,CAAC;gBAC3B,EAAE,EAAE,OAAO;oBACT,CAAC,CAAC,OAAO,CAAC,EAAE,CAAC,SAAS,CAAC,CAAC,CAAC;oBACzB,CAAC,CAAC,IAAA,iCAAS,EAAC,OAAO,CAAC,EAAE,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;aAC3C,CAAC;YACF,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,QAAQ,CAAC,IAAI,EAAE,UAAU,CAAC,CAAC;YAC9C,QAAQ,CAAC,SAAS,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;YAC1B,QAAQ,CAAC,SAAS,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;SAC3B;QACD,SAAS,CAAC,EAAE,GAAG,QAAQ,CAAC;QACxB,SAAS,CAAC,EAAE,GAAG,QAAQ,CAAC;KACzB;SAAM;QACL,2BAA2B;QAC3B,KAAK,MAAM,OAAO,IAAI,OAAO,EAAE;YAC7B,MAAM,QAAQ,GAAG,IAAI,mBAAM,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;YACjD,MAAM,QAAQ,GAAG,IAAI,mBAAM,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;YACjD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,EAAE,CAAC,EAAE,EAAE;gBAClC,MAAM,IAAI,GAAG,aAAa,CAAC,OAAO,CAAC,OAAO,CAAC,EAAE,EAAE,MAAM,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC,CAAC;gBACnE,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,QAAQ,CAAC,IAAI,EAAE,UAAU,CAAC,CAAC;gBAC9C,QAAQ,CAAC,SAAS,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;gBAC1B,QAAQ,CAAC,SAAS,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;aAC3B;YACD,IAAI,OAAO,KAAK,IAAI,EAAE;gBACpB,SAAS,CAAC,EAAE,GAAG,QAAQ,CAAC;gBACxB,SAAS,CAAC,EAAE,GAAG,QAAQ,CAAC;aACzB;iBAAM;gBACL,SAAS,CAAC,EAAE,GAAG,QAAQ,CAAC;gBACxB,SAAS,CAAC,EAAE,GAAG,QAAQ,CAAC;aACzB;SACF;KACF;IACD,OAAO,SAAS,CAAC;AACnB,CAAC;AA7DD,oDA6DC;AAED,SAAS,QAAQ,CAAC,IAAc,EAAE,OAAiC;IACjE,MAAM,EACJ,MAAM,EACN,WAAW,EAAE,aAAa,EAC1B,WAAW,EAAE,aAAa,GAAG,EAAE,EAC/B,eAAe,GAAG,EAAE,GACrB,GAAG,OAAO,CAAC;IAEZ,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,IAAA,yBAAW,EAAC,IAAI,EAAE,aAAa,CAAC,CAAC;IAEpD,MAAM,EAAE,QAAQ,GAAG,CAAC,EAAE,GAAG,aAAa,CAAC;IACvC,MAAM,WAAW,GAAG,IAAA,+BAAO,EAAC,IAAA,yBAAW,EAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,QAAQ,EAAE,CAAC,EAAE;QACjE,cAAc,EAAE,CAAC,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,MAAM,KAAK,CAAC,CAAC;KAC3D,CAAC,CAAC;IAEH,IAAI,OAAO,CAAC,eAAe,EAAE,IAAI,KAAK,IAAI,EAAE;QAC1C,MAAM,CAAC,MAAM,CAAC,WAAW,EAAE;YACzB,EAAE,EAAE,IAAA,oCAAY,EAAC,WAAW,CAAC;YAC7B,EAAE,EAAE,IAAI,YAAY,CAAC,QAAQ,CAAC;SAC/B,CAAC,CAAC;KACJ;IAED,MAAM,EAAE,IAAI,EAAE,QAAQ,EAAE,GAAG,EAAE,KAAK,EAAE,GAAG,EAAE,KAAK,EAAE,GAAG,eAAe,CAAC;IACnE,IAAI,QAAQ,KAAK,IAAI,IAAI,KAAK,KAAK,SAAS,IAAI,KAAK,KAAK,SAAS,EAAE;QACnE,MAAM,CAAC,MAAM,CACX,WAAW,EACX,IAAA,2CAAmB,EACjB,WAAW,EACX,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,EACvB,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CACxB,CACF,CAAC;KACH;IACD,OAAO,WAAW,CAAC;AACrB,CAAC;AAED,SAAS,aAAa,CACpB,IAAY,EACZ,IAA2D;IAE3D,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,GAAG,IAAI,CAAC;IAC/B,MAAM,OAAO,GAAG,IAAI,CAAC,OAAO,IAAI,MAAM,KAAK,IAAI,CAAC;IAChD,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC;IAC/B,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,QAAQ,CAAC,CAAC;IACtC,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,QAAQ,CAAC,CAAC;IAEtC,QAAQ,MAAM,EAAE,iBAAiB,EAAE,EAAE;QACnC,iBAAiB;QACjB,wBAAwB;QACxB,KAAK,KAAK,CAAC;QACX,KAAK,MAAM;YACT,OAAO;gBACL,EAAE,EAAE,IAAI,CAAC,SAAS,CAAC,KAAK,CAAC;gBACzB,EAAE,EAAE,IAAI,YAAY,CAAC,IAAI,CAAC,IAAI,CAAC;aAChC,CAAC;QACJ;YACE,IAAI,OAAO,EAAE;gBACX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE;oBACjC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,KAAK,CAAC,CAAC;oBAC/B,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,KAAK,CAAC,CAAC;iBACrC;aACF;iBAAM;gBACL,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE;oBACjC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,KAAK,CAAC,CAAC;oBAC/B,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,KAAK,CAAC,CAAC;iBACpC;aACF;KACJ;IAED,IAAI,MAAM,KAAK,aAAa,EAAE;QAC5B,8CAA8C;QAC9C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,IAAI,CAAC,EAAE;YACpC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;YACZ,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;SACb;KACF;IAED,OAAO,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;AACpB,CAAC"}

package/lib/xyz/util/fft2d/removeDCOffset.js

@@ -3,7 +3,7 @@ Object.defineProperty(exports, "__esModule", { value: true });
 exports.removeDCOffset = void 0;
 const ml_spectra_processing_1 = require("ml-spectra-processing");
 function removeDCOffset(data, options) {
-    let { digitalFilterValue = 0 } = options;
+    const { digitalFilterValue = 0 } = options;
     const nbPoints = data.re.length;
     const newRe = new Float64Array(data.re);
     const newIm = new Float64Array(data.im);

package/lib/xyz/util/fft2d/removeDCOffset.js.map

@@ -1 +1 @@
-{"version":3,"file":"removeDCOffset.js","sourceRoot":"","sources":["../../../../src/xyz/util/fft2d/removeDCOffset.ts"],"names":[],"mappings":";;;AAAA,iEAAwD;AAExD,SAAgB,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,IAAA,6BAAK,EACrB,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,QAAQ,GAAG,kBAAkB,CAAC,CACrE,CAAC;IACF,MAAM,SAAS,GAAG,IAAA,6BAAK,EACrB,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;AAhBD,wCAgBC"}
+{"version":3,"file":"removeDCOffset.js","sourceRoot":"","sources":["../../../../src/xyz/util/fft2d/removeDCOffset.ts"],"names":[],"mappings":";;;AAAA,iEAAwD;AAExD,SAAgB,cAAc,CAAC,IAAc,EAAE,OAAY;IACzD,MAAM,EAAE,kBAAkB,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;IAC3C,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,IAAA,6BAAK,EACrB,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,QAAQ,GAAG,kBAAkB,CAAC,CACrE,CAAC;IACF,MAAM,SAAS,GAAG,IAAA,6BAAK,EACrB,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;AAhBD,wCAgBC"}

package/lib/xyz/util/fft2d/zeroFilling.js

@@ -2,10 +2,8 @@
 Object.defineProperty(exports, "__esModule", { value: true });
 exports.zeroFilling = void 0;
 function zeroFilling(data, options) {
-    let { nbPoints, factor = 2, digitalFilterValue: grpdly = 0 } = options;
-    if (!nbPoints) {
-        nbPoints = 2 ** Math.round(Math.log2(data.re.length * factor));
-    }
+    const { factor = 2, digitalFilterValue: grpdly = 0 } = options;
+    const nbPoints = options.nbPoints || 2 ** Math.round(Math.log2(data.re.length * factor));
     let pointsToShift;
     if (grpdly > 0) {
         pointsToShift = Math.floor(grpdly);
@@ -14,8 +12,8 @@ function zeroFilling(data, options) {
         pointsToShift = 0;
     }
     const { re, im } = data;
-    let newRE = new Float64Array(nbPoints);
-    let newIM = new Float64Array(nbPoints);
+    const newRE = new Float64Array(nbPoints);
+    const newIM = new Float64Array(nbPoints);
     const length = Math.min(nbPoints, re.length);
     newRE.set(re.slice(0, length - pointsToShift));
     newIM.set(im.slice(0, length - pointsToShift));

package/lib/xyz/util/fft2d/zeroFilling.js.map

@@ -1 +1 @@
-{"version":3,"file":"zeroFilling.js","sourceRoot":"","sources":["../../../../src/xyz/util/fft2d/zeroFilling.ts"],"names":[],"mappings":";;;AAQA,SAAgB,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;AA9BD,kCA8BC"}
+{"version":3,"file":"zeroFilling.js","sourceRoot":"","sources":["../../../../src/xyz/util/fft2d/zeroFilling.ts"],"names":[],"mappings":";;;AAQA,SAAgB,WAAW,CAAC,IAAc,EAAE,OAA2B;IACrE,MAAM,EAAE,MAAM,GAAG,CAAC,EAAE,kBAAkB,EAAE,MAAM,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;IAE/D,MAAM,QAAQ,GACZ,OAAO,CAAC,QAAQ,IAAI,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC;IAE1E,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,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,QAAQ,CAAC,CAAC;IACzC,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,QAAQ,CAAC,CAAC;IAEzC,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;AA7BD,kCA6BC"}

package/lib/xyz/util/formatZone.js

@@ -3,10 +3,10 @@ Object.defineProperty(exports, "__esModule", { value: true });
 exports.formatZones = void 0;
 const uuid_1 = require("@lukeed/uuid");
 function formatZones(signals) {
-    let zones = [];
+    const zones = [];
     for (const signal of signals) {
-        let minMax1 = [Number.MAX_VALUE, 0];
-        let minMax2 = [Number.MAX_VALUE, 0];
+        const minMax1 = [Number.MAX_VALUE, 0];
+        const minMax2 = [Number.MAX_VALUE, 0];
         for (const peak of signal.peaks || []) {
             if (peak.minX < minMax1[0]) {
                 minMax1[0] = peak.minX;

package/lib/xyz/util/formatZone.js.map

@@ -1 +1 @@
-{"version":3,"file":"formatZone.js","sourceRoot":"","sources":["../../../src/xyz/util/formatZone.ts"],"names":[],"mappings":";;;AAAA,uCAAgD;AAKhD,SAAgB,WAAW,CAAC,OAA4B;IACtD,IAAI,KAAK,GAAoB,EAAE,CAAC;IAChC,KAAK,MAAM,MAAM,IAAI,OAAO,EAAE;QAC5B,IAAI,OAAO,GAAG,CAAC,MAAM,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;QACpC,IAAI,OAAO,GAAG,CAAC,MAAM,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;QACpC,KAAK,MAAM,IAAI,IAAI,MAAM,CAAC,KAAK,IAAI,EAAE,EAAE;YACrC,IAAI,IAAI,CAAC,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,EAAE;gBAC1B,OAAO,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;aACxB;YACD,IAAI,IAAI,CAAC,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,EAAE;gBAC1B,OAAO,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;aACxB;YACD,IAAI,IAAI,CAAC,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,EAAE;gBAC1B,OAAO,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;aACxB;YACD,IAAI,IAAI,CAAC,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,EAAE;gBAC1B,OAAO,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;aACxB;SACF;QACD,KAAK,CAAC,IAAI,CAAC;YACT,EAAE,EAAE,IAAA,SAAU,GAAE;YAChB,CAAC,EAAE;gBACD,IAAI,EAAE,OAAO,CAAC,CAAC,CAAC;gBAChB,EAAE,EAAE,OAAO,CAAC,CAAC,CAAC;aACf;YACD,CAAC,EAAE;gBACD,IAAI,EAAE,OAAO,CAAC,CAAC,CAAC;gBAChB,EAAE,EAAE,OAAO,CAAC,CAAC,CAAC;aACf;YACD,OAAO,EAAE,CAAC,MAAM,CAAC;SAClB,CAAC,CAAC;KACJ;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAjCD,kCAiCC"}
+{"version":3,"file":"formatZone.js","sourceRoot":"","sources":["../../../src/xyz/util/formatZone.ts"],"names":[],"mappings":";;;AAAA,uCAAgD;AAKhD,SAAgB,WAAW,CAAC,OAA4B;IACtD,MAAM,KAAK,GAAoB,EAAE,CAAC;IAClC,KAAK,MAAM,MAAM,IAAI,OAAO,EAAE;QAC5B,MAAM,OAAO,GAAG,CAAC,MAAM,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;QACtC,MAAM,OAAO,GAAG,CAAC,MAAM,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;QACtC,KAAK,MAAM,IAAI,IAAI,MAAM,CAAC,KAAK,IAAI,EAAE,EAAE;YACrC,IAAI,IAAI,CAAC,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,EAAE;gBAC1B,OAAO,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;aACxB;YACD,IAAI,IAAI,CAAC,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,EAAE;gBAC1B,OAAO,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;aACxB;YACD,IAAI,IAAI,CAAC,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,EAAE;gBAC1B,OAAO,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;aACxB;YACD,IAAI,IAAI,CAAC,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,EAAE;gBAC1B,OAAO,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;aACxB;SACF;QACD,KAAK,CAAC,IAAI,CAAC;YACT,EAAE,EAAE,IAAA,SAAU,GAAE;YAChB,CAAC,EAAE;gBACD,IAAI,EAAE,OAAO,CAAC,CAAC,CAAC;gBAChB,EAAE,EAAE,OAAO,CAAC,CAAC,CAAC;aACf;YACD,CAAC,EAAE;gBACD,IAAI,EAAE,OAAO,CAAC,CAAC,CAAC;gBAChB,EAAE,EAAE,OAAO,CAAC,CAAC,CAAC;aACf;YACD,OAAO,EAAE,CAAC,MAAM,CAAC;SAClB,CAAC,CAAC;KACJ;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAjCD,kCAiCC"}

package/lib/xyz/util/getMinMaxXY.d.ts

@@ -0,0 +1,6 @@
+export declare function getMinMaxXY(info: Record<string, any>): {
+    minX: number;
+    maxX: number;
+    minY: number;
+    maxY: number;
+};

package/lib/xyz/util/getMinMaxXY.js

@@ -0,0 +1,23 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.getMinMaxXY = void 0;
+function getMinMaxXY(info) {
+    const min = [];
+    const max = [];
+    for (let i = 0; i < info.baseFrequency.length; i++) {
+        const baseFrequency = info.baseFrequency[i];
+        const frequencyOffset = info.frequencyOffset[i];
+        const spectralWidth = info.spectralWidth[i];
+        const offset = frequencyOffset / baseFrequency;
+        min.push(offset - spectralWidth * 0.5);
+        max.push(offset + spectralWidth * 0.5);
+    }
+    return {
+        minX: min[0],
+        maxX: max[0],
+        minY: min[1],
+        maxY: max[1],
+    };
+}
+exports.getMinMaxXY = getMinMaxXY;
+//# sourceMappingURL=getMinMaxXY.js.map

package/lib/xyz/util/getMinMaxXY.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"getMinMaxXY.js","sourceRoot":"","sources":["../../../src/xyz/util/getMinMaxXY.ts"],"names":[],"mappings":";;;AAAA,SAAgB,WAAW,CAAC,IAAyB;IACnD,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,IAAI,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC;QAC5C,MAAM,eAAe,GAAG,IAAI,CAAC,eAAe,CAAC,CAAC,CAAC,CAAC;QAChD,MAAM,aAAa,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC;QAC5C,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;AAlBD,kCAkBC"}

package/lib/xyz/util/padData.js

@@ -2,20 +2,20 @@
 Object.defineProperty(exports, "__esModule", { value: true });
 exports.padData = void 0;
 function padData(spectraData, options) {
-    let { minX, maxX, minY, maxY } = spectraData;
+    const { minX, maxX, minY, maxY } = spectraData;
     const width = options.width;
-    let nbPoints = spectraData.z[0].length;
-    let nbSubSpectra = spectraData.z.length;
-    let yInterval = (maxY - minY) / (nbSubSpectra - 1);
-    let xInterval = (maxX - minX) / (nbPoints - 1);
+    const nbPoints = spectraData.z[0].length;
+    const nbSubSpectra = spectraData.z.length;
+    const yInterval = (maxY - minY) / (nbSubSpectra - 1);
+    const xInterval = (maxX - minX) / (nbPoints - 1);
     let yDiff = width - nbSubSpectra;
     let xDiff = Math.max(width - nbPoints, 0);
     if (xDiff % 2)
         xDiff++;
     if (yDiff % 2)
         yDiff++;
-    let xOffset = xDiff / 2;
-    let yOffset = yDiff / 2;
+    const xOffset = xDiff / 2;
+    const yOffset = yDiff / 2;
     const newMatrix = [];
     for (let i = 0; i < nbSubSpectra + yDiff; i++) {
         newMatrix.push(new Float64Array(nbPoints + xDiff));

package/lib/xyz/util/padData.js.map

@@ -1 +1 @@
-{"version":3,"file":"padData.js","sourceRoot":"","sources":["../../../src/xyz/util/padData.ts"],"names":[],"mappings":";;;AAEA,SAAgB,OAAO,CACrB,WAA6B,EAC7B,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;AAxCD,0BAwCC"}
+{"version":3,"file":"padData.js","sourceRoot":"","sources":["../../../src/xyz/util/padData.ts"],"names":[],"mappings":";;;AAEA,SAAgB,OAAO,CACrB,WAA6B,EAC7B,OAA0B;IAE1B,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,WAAW,CAAC;IAE/C,MAAM,KAAK,GAAG,OAAO,CAAC,KAAK,CAAC;IAE5B,MAAM,QAAQ,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;IACzC,MAAM,YAAY,GAAG,WAAW,CAAC,CAAC,CAAC,MAAM,CAAC;IAE1C,MAAM,SAAS,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC;IACrD,MAAM,SAAS,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,CAAC,CAAC,CAAC;IAEjD,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,MAAM,OAAO,GAAG,KAAK,GAAG,CAAC,CAAC;IAC1B,MAAM,OAAO,GAAG,KAAK,GAAG,CAAC,CAAC;IAE1B,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;AAxCD,0BAwCC"}

package/lib/xyz/util/phaseCorrection/applyOverImag.d.ts

@@ -0,0 +1,24 @@
+import { NmrData2DFt } from 'cheminfo-types';
+import { PhaseCorrectionOptions } from '../../xyzPhaseCorrectionDirectDimension';
+export declare function applyOverImag<T extends NmrData2DFt>(data: T, options: PhaseCorrectionOptions): T & {
+    ri: {
+        z: Float64Array[];
+        maxZ: number;
+        minZ: number;
+        minY: number;
+        maxY: number;
+        minX: number;
+        maxX: number;
+    };
+    ii: {
+        z: Float64Array[];
+        maxZ: number;
+        minZ: number;
+        minY: number;
+        maxY: number;
+        minX: number;
+        maxX: number;
+    };
+    rr: import("cheminfo-types").NmrData2DContent;
+    ir: import("cheminfo-types").NmrData2DContent;
+};

package/lib/xyz/util/phaseCorrection/applyOverImag.js

@@ -0,0 +1,27 @@
+"use strict";
+var __importDefault = (this && this.__importDefault) || function (mod) {
+    return (mod && mod.__esModule) ? mod : { "default": mod };
+};
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.applyOverImag = void 0;
+const getZMinMax_1 = require("./getZMinMax");
+const ml_matrix_1 = __importDefault(require("ml-matrix"));
+const checkNmrData2DFourQuad_1 = require("./checkNmrData2DFourQuad");
+function applyOverImag(data, options) {
+    (0, checkNmrData2DFourQuad_1.checkNmrData2DFourQuad)(data);
+    const { applyPhaseCorrection } = options;
+    const { re, im } = applyPhaseCorrection({ re: ml_matrix_1.default.checkMatrix(data.ri.z), im: ml_matrix_1.default.checkMatrix(data.ii.z) }, options);
+    return {
+        ...data,
+        ri: {
+            ...data.ri,
+            ...(0, getZMinMax_1.getZMinMax)(re),
+        },
+        ii: {
+            ...data.ii,
+            ...(0, getZMinMax_1.getZMinMax)(im),
+        },
+    };
+}
+exports.applyOverImag = applyOverImag;
+//# sourceMappingURL=applyOverImag.js.map

package/lib/xyz/util/phaseCorrection/applyOverImag.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"applyOverImag.js","sourceRoot":"","sources":["../../../../src/xyz/util/phaseCorrection/applyOverImag.ts"],"names":[],"mappings":";;;;;;AAEA,6CAA0C;AAC1C,0DAA+B;AAC/B,qEAAkE;AAElE,SAAgB,aAAa,CAC3B,IAAO,EACP,OAA+B;IAE/B,IAAA,+CAAsB,EAAC,IAAI,CAAC,CAAC;IAC7B,MAAM,EAAE,oBAAoB,EAAE,GAAG,OAAO,CAAC;IACzC,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,oBAAoB,CACrC,EAAE,EAAE,EAAE,mBAAM,CAAC,WAAW,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,mBAAM,CAAC,WAAW,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EACxE,OAAO,CACR,CAAC;IAEF,OAAO;QACL,GAAG,IAAI;QACP,EAAE,EAAE;YACF,GAAG,IAAI,CAAC,EAAE;YACV,GAAG,IAAA,uBAAU,EAAC,EAAE,CAAC;SAClB;QACD,EAAE,EAAE;YACF,GAAG,IAAI,CAAC,EAAE;YACV,GAAG,IAAA,uBAAU,EAAC,EAAE,CAAC;SAClB;KACF,CAAC;AACJ,CAAC;AAtBD,sCAsBC"}

package/lib/xyz/util/phaseCorrection/applyOverReal.d.ts

@@ -0,0 +1,22 @@
+import { NmrData2DFt } from 'cheminfo-types';
+import { PhaseCorrectionOptions } from '../../xyzPhaseCorrectionDirectDimension';
+export declare function applyOverReal<T extends NmrData2DFt>(data: T, options: PhaseCorrectionOptions): T & {
+    rr: {
+        z: Float64Array[];
+        maxZ: number;
+        minZ: number;
+        minY: number;
+        maxY: number;
+        minX: number;
+        maxX: number;
+    };
+    ir: {
+        z: Float64Array[];
+        maxZ: number;
+        minZ: number;
+        minY: number;
+        maxY: number;
+        minX: number;
+        maxX: number;
+    };
+};

package/lib/xyz/util/phaseCorrection/applyOverReal.js

@@ -0,0 +1,27 @@
+"use strict";
+var __importDefault = (this && this.__importDefault) || function (mod) {
+    return (mod && mod.__esModule) ? mod : { "default": mod };
+};
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.applyOverReal = void 0;
+const checkNmrData2DTwoQuad_1 = require("./checkNmrData2DTwoQuad");
+const ml_matrix_1 = __importDefault(require("ml-matrix"));
+const getZMinMax_1 = require("./getZMinMax");
+function applyOverReal(data, options) {
+    const { applyPhaseCorrection } = options;
+    (0, checkNmrData2DTwoQuad_1.checkNmrData2DTwoQuad)(data);
+    const { re, im } = applyPhaseCorrection({ re: ml_matrix_1.default.checkMatrix(data.rr.z), im: ml_matrix_1.default.checkMatrix(data.ir.z) }, options);
+    return {
+        ...data,
+        rr: {
+            ...data.rr,
+            ...(0, getZMinMax_1.getZMinMax)(re),
+        },
+        ir: {
+            ...data.ir,
+            ...(0, getZMinMax_1.getZMinMax)(im),
+        },
+    };
+}
+exports.applyOverReal = applyOverReal;
+//# sourceMappingURL=applyOverReal.js.map

package/lib/xyz/util/phaseCorrection/applyOverReal.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"applyOverReal.js","sourceRoot":"","sources":["../../../../src/xyz/util/phaseCorrection/applyOverReal.ts"],"names":[],"mappings":";;;;;;AACA,mEAAgE;AAChE,0DAA+B;AAC/B,6CAA0C;AAG1C,SAAgB,aAAa,CAC3B,IAAO,EACP,OAA+B;IAE/B,MAAM,EAAE,oBAAoB,EAAE,GAAG,OAAO,CAAC;IACzC,IAAA,6CAAqB,EAAC,IAAI,CAAC,CAAC;IAC5B,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,oBAAoB,CACrC,EAAE,EAAE,EAAE,mBAAM,CAAC,WAAW,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,mBAAM,CAAC,WAAW,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EACxE,OAAO,CACR,CAAC;IAEF,OAAO;QACL,GAAG,IAAI;QACP,EAAE,EAAE;YACF,GAAG,IAAI,CAAC,EAAE;YACV,GAAG,IAAA,uBAAU,EAAC,EAAE,CAAC;SAClB;QACD,EAAE,EAAE;YACF,GAAG,IAAI,CAAC,EAAE;YACV,GAAG,IAAA,uBAAU,EAAC,EAAE,CAAC;SAClB;KACF,CAAC;AACJ,CAAC;AAtBD,sCAsBC"}

package/lib/xyz/util/phaseCorrection/checkNmrData2DFourQuad.d.ts

@@ -0,0 +1,3 @@
+import { NmrData2DFt } from 'cheminfo-types';
+export type NmrData2DFourQuadrant = Required<Pick<NmrData2DFt, 'rr' | 'ir' | 'ri' | 'ii'>>;
+export declare function checkNmrData2DFourQuad(data: NmrData2DFt): asserts data is NmrData2DFourQuadrant;

package/lib/xyz/util/phaseCorrection/checkNmrData2DFourQuad.js

@@ -0,0 +1,10 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.checkNmrData2DFourQuad = void 0;
+function checkNmrData2DFourQuad(data) {
+    if (!('ir' in data)) {
+        throw new Error('four quadrants should has ir, ri and ii quadrants');
+    }
+}
+exports.checkNmrData2DFourQuad = checkNmrData2DFourQuad;
+//# sourceMappingURL=checkNmrData2DFourQuad.js.map

package/lib/xyz/util/phaseCorrection/checkNmrData2DFourQuad.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"checkNmrData2DFourQuad.js","sourceRoot":"","sources":["../../../../src/xyz/util/phaseCorrection/checkNmrData2DFourQuad.ts"],"names":[],"mappings":";;;AAKA,SAAgB,sBAAsB,CACpC,IAAiB;IAEjB,IAAI,CAAC,CAAC,IAAI,IAAI,IAAI,CAAC,EAAE;QACnB,MAAM,IAAI,KAAK,CAAC,mDAAmD,CAAC,CAAC;KACtE;AACH,CAAC;AAND,wDAMC"}

package/lib/xyz/util/phaseCorrection/checkNmrData2DTwoQuad.d.ts

@@ -0,0 +1,3 @@
+import { NmrData2DFt } from 'cheminfo-types';
+export type NmrData2DTwoQuadrant = Required<Pick<NmrData2DFt, 'rr' | 'ir'>>;
+export declare function checkNmrData2DTwoQuad(data: NmrData2DFt): asserts data is NmrData2DTwoQuadrant;

package/lib/xyz/util/phaseCorrection/checkNmrData2DTwoQuad.js

@@ -0,0 +1,10 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.checkNmrData2DTwoQuad = void 0;
+function checkNmrData2DTwoQuad(data) {
+    if (!('ir' in data)) {
+        throw new Error('Two quadrants should has ir quadrant');
+    }
+}
+exports.checkNmrData2DTwoQuad = checkNmrData2DTwoQuad;
+//# sourceMappingURL=checkNmrData2DTwoQuad.js.map

package/lib/xyz/util/phaseCorrection/checkNmrData2DTwoQuad.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"checkNmrData2DTwoQuad.js","sourceRoot":"","sources":["../../../../src/xyz/util/phaseCorrection/checkNmrData2DTwoQuad.ts"],"names":[],"mappings":";;;AAIA,SAAgB,qBAAqB,CACnC,IAAiB;IAEjB,IAAI,CAAC,CAAC,IAAI,IAAI,IAAI,CAAC,EAAE;QACnB,MAAM,IAAI,KAAK,CAAC,sCAAsC,CAAC,CAAC;KACzD;AACH,CAAC;AAND,sDAMC"}

package/lib/xyz/util/phaseCorrection/getZMinMax.d.ts

@@ -0,0 +1,6 @@
+import Matrix from 'ml-matrix';
+export declare function getZMinMax(data: Matrix): {
+    z: Float64Array[];
+    maxZ: number;
+    minZ: number;
+};

package/lib/xyz/util/phaseCorrection/getZMinMax.js

@@ -0,0 +1,13 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.getZMinMax = void 0;
+const ml_spectra_processing_1 = require("ml-spectra-processing");
+function getZMinMax(data) {
+    return {
+        z: data.to2DArray().map(ml_spectra_processing_1.xEnsureFloat64),
+        maxZ: data.max(),
+        minZ: data.min(),
+    };
+}
+exports.getZMinMax = getZMinMax;
+//# sourceMappingURL=getZMinMax.js.map

package/lib/xyz/util/phaseCorrection/getZMinMax.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"getZMinMax.js","sourceRoot":"","sources":["../../../../src/xyz/util/phaseCorrection/getZMinMax.ts"],"names":[],"mappings":";;;AACA,iEAAuD;AAEvD,SAAgB,UAAU,CAAC,IAAY;IACrC,OAAO;QACL,CAAC,EAAE,IAAI,CAAC,SAAS,EAAE,CAAC,GAAG,CAAC,sCAAc,CAAC;QACvC,IAAI,EAAE,IAAI,CAAC,GAAG,EAAE;QAChB,IAAI,EAAE,IAAI,CAAC,GAAG,EAAE;KACjB,CAAC;AACJ,CAAC;AAND,gCAMC"}