nmr-processing 9.0.1 → 9.0.2
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/peaks/peaksToXY.d.ts +10 -0
- package/lib/peaks/peaksToXY.js +20 -4
- package/lib/peaks/peaksToXY.js.map +1 -1
- package/lib/peaks/util/convertWidth.d.ts +2 -1
- package/lib/peaks/util/convertWidth.js +4 -4
- package/lib/peaks/util/convertWidth.js.map +1 -1
- package/lib/peaks/util/convertWidthToHz.d.ts +3 -1
- package/lib/peaks/util/convertWidthToHz.js +2 -2
- package/lib/peaks/util/convertWidthToHz.js.map +1 -1
- package/lib/peaks/util/convertWidthToPPM.d.ts +2 -1
- package/lib/peaks/util/convertWidthToPPM.js +2 -2
- package/lib/peaks/util/convertWidthToPPM.js.map +1 -1
- package/lib/xy/xyAutoPeaksPicking.js +4 -5
- package/lib/xy/xyAutoPeaksPicking.js.map +1 -1
- package/lib/xy/xyPeaksOptimization.d.ts +7 -1
- package/lib/xy/xyPeaksOptimization.js +9 -1
- package/lib/xy/xyPeaksOptimization.js.map +1 -1
- package/lib-esm/peaks/peaksToXY.js +20 -4
- package/lib-esm/peaks/peaksToXY.js.map +1 -1
- package/lib-esm/peaks/util/convertWidth.js +4 -4
- package/lib-esm/peaks/util/convertWidth.js.map +1 -1
- package/lib-esm/peaks/util/convertWidthToHz.js +2 -2
- package/lib-esm/peaks/util/convertWidthToHz.js.map +1 -1
- package/lib-esm/peaks/util/convertWidthToPPM.js +2 -2
- package/lib-esm/peaks/util/convertWidthToPPM.js.map +1 -1
- package/lib-esm/xy/xyAutoPeaksPicking.js +4 -5
- package/lib-esm/xy/xyAutoPeaksPicking.js.map +1 -1
- package/lib-esm/xy/xyPeaksOptimization.js +9 -1
- package/lib-esm/xy/xyPeaksOptimization.js.map +1 -1
- package/package.json +1 -1
- package/src/peaks/peaksToXY.ts +37 -4
- package/src/peaks/util/convertWidth.ts +11 -4
- package/src/peaks/util/convertWidthToHz.ts +5 -3
- package/src/peaks/util/convertWidthToPPM.ts +3 -2
- package/src/xy/xyAutoPeaksPicking.ts +10 -9
- package/src/xy/xyPeaksOptimization.ts +17 -2
package/lib/peaks/peaksToXY.d.ts
CHANGED
|
@@ -4,6 +4,16 @@ import { NMRPeak1D } from './NMRPeak1D';
|
|
|
4
4
|
* Convert an array of peaks assuming that width/fwhm is in Hz
|
|
5
5
|
*/
|
|
6
6
|
export interface PeaksToXYOptions extends OptionsSG1D {
|
|
7
|
+
/**
|
|
8
|
+
* First x value (inclusive).
|
|
9
|
+
* @default `firstPeak.x + firstPeak.width * 3`
|
|
10
|
+
*/
|
|
11
|
+
from?: number;
|
|
12
|
+
/**
|
|
13
|
+
* Last x value (inclusive).
|
|
14
|
+
* @default `lastPeak.x + lastPeak.width * 3`
|
|
15
|
+
*/
|
|
16
|
+
to?: number;
|
|
7
17
|
/**
|
|
8
18
|
* frequency of the spectrum
|
|
9
19
|
*/
|
package/lib/peaks/peaksToXY.js
CHANGED
|
@@ -4,15 +4,31 @@ exports.peaksToXY = void 0;
|
|
|
4
4
|
const spectrum_generator_1 = require("spectrum-generator");
|
|
5
5
|
const convertWidthToPPM_1 = require("./util/convertWidthToPPM");
|
|
6
6
|
function peaksToXY(peaks, options) {
|
|
7
|
-
const {
|
|
8
|
-
|
|
7
|
+
const { frequency, nbPoints = 1024, shape } = options;
|
|
8
|
+
const newPeaks = (0, convertWidthToPPM_1.convertWidthToPPM)(peaks, { frequency });
|
|
9
|
+
return (0, spectrum_generator_1.generateSpectrum)(newPeaks, {
|
|
9
10
|
generator: {
|
|
10
|
-
|
|
11
|
-
to,
|
|
11
|
+
...getFromTo(newPeaks, options),
|
|
12
12
|
nbPoints,
|
|
13
13
|
shape,
|
|
14
14
|
},
|
|
15
15
|
});
|
|
16
16
|
}
|
|
17
17
|
exports.peaksToXY = peaksToXY;
|
|
18
|
+
function getFromTo(newPeaks, options) {
|
|
19
|
+
if ('to' in options && 'from' in options) {
|
|
20
|
+
return {
|
|
21
|
+
from: options.from,
|
|
22
|
+
to: options.to,
|
|
23
|
+
};
|
|
24
|
+
}
|
|
25
|
+
newPeaks.sort((a, b) => a.x - b.x);
|
|
26
|
+
const firstPeak = newPeaks[0];
|
|
27
|
+
const lastPeak = newPeaks[newPeaks.length - 1];
|
|
28
|
+
const { from = firstPeak.x - firstPeak.width * 2, to = lastPeak.x + lastPeak.width * 2, } = options;
|
|
29
|
+
return {
|
|
30
|
+
from,
|
|
31
|
+
to,
|
|
32
|
+
};
|
|
33
|
+
}
|
|
18
34
|
//# sourceMappingURL=peaksToXY.js.map
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"peaksToXY.js","sourceRoot":"","sources":["../../src/peaks/peaksToXY.ts"],"names":[],"mappings":";;;AAAA,2DAAmE;AAGnE,gEAA6D;
|
|
1
|
+
{"version":3,"file":"peaksToXY.js","sourceRoot":"","sources":["../../src/peaks/peaksToXY.ts"],"names":[],"mappings":";;;AAAA,2DAAmE;AAGnE,gEAA6D;AAuB7D,SAAgB,SAAS,CAAC,KAAkB,EAAE,OAAyB;IACrE,MAAM,EAAE,SAAS,EAAE,QAAQ,GAAG,IAAI,EAAE,KAAK,EAAE,GAAG,OAAO,CAAC;IACtD,MAAM,QAAQ,GAAG,IAAA,qCAAiB,EAAC,KAAK,EAAE,EAAE,SAAS,EAAE,CAAC,CAAC;IAEzD,OAAO,IAAA,qCAAgB,EAAC,QAAQ,EAAE;QAChC,SAAS,EAAE;YACT,GAAG,SAAS,CAAC,QAAQ,EAAE,OAAO,CAAC;YAC/B,QAAQ;YACR,KAAK;SACN;KACF,CAAC,CAAC;AACL,CAAC;AAXD,8BAWC;AAED,SAAS,SAAS,CAAC,QAAqB,EAAE,OAAyB;IACjE,IAAI,IAAI,IAAI,OAAO,IAAI,MAAM,IAAI,OAAO,EAAE;QACxC,OAAO;YACL,IAAI,EAAE,OAAO,CAAC,IAAI;YAClB,EAAE,EAAE,OAAO,CAAC,EAAE;SACf,CAAC;KACH;IAED,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACnC,MAAM,SAAS,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;IAC9B,MAAM,QAAQ,GAAG,QAAQ,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IAE/C,MAAM,EACJ,IAAI,GAAG,SAAS,CAAC,CAAC,GAAG,SAAS,CAAC,KAAK,GAAG,CAAC,EACxC,EAAE,GAAG,QAAQ,CAAC,CAAC,GAAG,QAAQ,CAAC,KAAK,GAAG,CAAC,GACrC,GAAG,OAAO,CAAC;IAEZ,OAAO;QACL,IAAI;QACJ,EAAE;KACH,CAAC;AACJ,CAAC"}
|
|
@@ -6,8 +6,9 @@ export interface ConvertWidthOptions {
|
|
|
6
6
|
_kind?: 'hz' | 'ppm';
|
|
7
7
|
};
|
|
8
8
|
frequency: number;
|
|
9
|
+
output?: NMRPeak1D[];
|
|
9
10
|
}
|
|
10
11
|
/**
|
|
11
12
|
* convert width and fwhm to ppm
|
|
12
13
|
*/
|
|
13
|
-
export declare function convertWidth(peaks: NMRPeak1D[], options: ConvertWidthOptions):
|
|
14
|
+
export declare function convertWidth(peaks: NMRPeak1D[], options: ConvertWidthOptions): any;
|
|
@@ -5,17 +5,17 @@ exports.convertWidth = void 0;
|
|
|
5
5
|
* convert width and fwhm to ppm
|
|
6
6
|
*/
|
|
7
7
|
function convertWidth(peaks, options) {
|
|
8
|
-
const { frequency, convertTo } = options;
|
|
8
|
+
const { frequency, convertTo, output = JSON.parse(JSON.stringify(peaks)), } = options;
|
|
9
9
|
const convert = getConverter(convertTo, frequency);
|
|
10
|
-
|
|
10
|
+
for (const peak of output) {
|
|
11
11
|
peak.width = convert(peak.width);
|
|
12
12
|
if (hasShape(peak)) {
|
|
13
13
|
if (peak.shape.fwhm) {
|
|
14
14
|
peak.shape.fwhm = convert(peak.shape.fwhm);
|
|
15
15
|
}
|
|
16
16
|
}
|
|
17
|
-
|
|
18
|
-
|
|
17
|
+
}
|
|
18
|
+
return output;
|
|
19
19
|
}
|
|
20
20
|
exports.convertWidth = convertWidth;
|
|
21
21
|
function hasShape(peak) {
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"convertWidth.js","sourceRoot":"","sources":["../../../src/peaks/util/convertWidth.ts"],"names":[],"mappings":";;;
|
|
1
|
+
{"version":3,"file":"convertWidth.js","sourceRoot":"","sources":["../../../src/peaks/util/convertWidth.ts"],"names":[],"mappings":";;;AAUA;;GAEG;AACH,SAAgB,YAAY,CAAC,KAAkB,EAAE,OAA4B;IAC3E,MAAM,EACJ,SAAS,EACT,SAAS,EACT,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC,GAC3C,GAAG,OAAO,CAAC;IAEZ,MAAM,OAAO,GAAG,YAAY,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;IAEnD,KAAK,MAAM,IAAI,IAAI,MAAM,EAAE;QACzB,IAAI,CAAC,KAAK,GAAG,OAAO,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACjC,IAAI,QAAQ,CAAC,IAAI,CAAC,EAAE;YAClB,IAAI,IAAI,CAAC,KAAK,CAAC,IAAI,EAAE;gBACnB,IAAI,CAAC,KAAK,CAAC,IAAI,GAAG,OAAO,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;aAC5C;SACF;KACF;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAlBD,oCAkBC;AAED,SAAS,QAAQ,CAAC,IAAe;IAC/B,OAAO,OAAO,IAAI,IAAI,CAAC;AACzB,CAAC;AAED,SAAS,YAAY,CAAC,SAAiB,EAAE,SAAiB;IACxD,QAAQ,SAAS,EAAE;QACjB,KAAK,KAAK;YACR,OAAO,CAAC,CAAS,EAAE,EAAE,CAAC,CAAC,GAAG,SAAS,CAAC;QACtC,KAAK,IAAI;YACP,OAAO,CAAC,CAAS,EAAE,EAAE,CAAC,CAAC,GAAG,SAAS,CAAC;QACtC;YACE,MAAM,IAAI,KAAK,CAAC,+BAA+B,SAAS,EAAE,CAAC,CAAC;KAC/D;AACH,CAAC"}
|
|
@@ -1,2 +1,4 @@
|
|
|
1
1
|
import { NMRPeak1D } from '../NMRPeak1D';
|
|
2
|
-
|
|
2
|
+
import { ConvertWidthOptions } from './convertWidth';
|
|
3
|
+
export declare type ConvertWidthTo = Omit<ConvertWidthOptions, 'convertTo'>;
|
|
4
|
+
export declare function convertWidthToHz(peaks: NMRPeak1D[], options: ConvertWidthTo): any;
|
|
@@ -2,8 +2,8 @@
|
|
|
2
2
|
Object.defineProperty(exports, "__esModule", { value: true });
|
|
3
3
|
exports.convertWidthToHz = void 0;
|
|
4
4
|
const convertWidth_1 = require("./convertWidth");
|
|
5
|
-
function convertWidthToHz(peaks,
|
|
6
|
-
return (0, convertWidth_1.convertWidth)(peaks, {
|
|
5
|
+
function convertWidthToHz(peaks, options) {
|
|
6
|
+
return (0, convertWidth_1.convertWidth)(peaks, { ...options, convertTo: 'hz' });
|
|
7
7
|
}
|
|
8
8
|
exports.convertWidthToHz = convertWidthToHz;
|
|
9
9
|
//# sourceMappingURL=convertWidthToHz.js.map
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"convertWidthToHz.js","sourceRoot":"","sources":["../../../src/peaks/util/convertWidthToHz.ts"],"names":[],"mappings":";;;AAEA,
|
|
1
|
+
{"version":3,"file":"convertWidthToHz.js","sourceRoot":"","sources":["../../../src/peaks/util/convertWidthToHz.ts"],"names":[],"mappings":";;;AAEA,iDAAmE;AAInE,SAAgB,gBAAgB,CAAC,KAAkB,EAAE,OAAuB;IAC1E,OAAO,IAAA,2BAAY,EAAC,KAAK,EAAE,EAAE,GAAG,OAAO,EAAE,SAAS,EAAE,IAAI,EAAE,CAAC,CAAC;AAC9D,CAAC;AAFD,4CAEC"}
|
|
@@ -1,2 +1,3 @@
|
|
|
1
1
|
import { NMRPeak1D } from '../NMRPeak1D';
|
|
2
|
-
|
|
2
|
+
import { ConvertWidthTo } from './convertWidthToHz';
|
|
3
|
+
export declare function convertWidthToPPM(peaks: NMRPeak1D[], options: ConvertWidthTo): any;
|
|
@@ -2,8 +2,8 @@
|
|
|
2
2
|
Object.defineProperty(exports, "__esModule", { value: true });
|
|
3
3
|
exports.convertWidthToPPM = void 0;
|
|
4
4
|
const convertWidth_1 = require("./convertWidth");
|
|
5
|
-
function convertWidthToPPM(peaks,
|
|
6
|
-
return (0, convertWidth_1.convertWidth)(peaks, {
|
|
5
|
+
function convertWidthToPPM(peaks, options) {
|
|
6
|
+
return (0, convertWidth_1.convertWidth)(peaks, { ...options, convertTo: 'ppm' });
|
|
7
7
|
}
|
|
8
8
|
exports.convertWidthToPPM = convertWidthToPPM;
|
|
9
9
|
//# sourceMappingURL=convertWidthToPPM.js.map
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"convertWidthToPPM.js","sourceRoot":"","sources":["../../../src/peaks/util/convertWidthToPPM.ts"],"names":[],"mappings":";;;AAEA,iDAA8C;
|
|
1
|
+
{"version":3,"file":"convertWidthToPPM.js","sourceRoot":"","sources":["../../../src/peaks/util/convertWidthToPPM.ts"],"names":[],"mappings":";;;AAEA,iDAA8C;AAG9C,SAAgB,iBAAiB,CAAC,KAAkB,EAAE,OAAuB;IAC3E,OAAO,IAAA,2BAAY,EAAC,KAAK,EAAE,EAAE,GAAG,OAAO,EAAE,SAAS,EAAE,KAAK,EAAE,CAAC,CAAC;AAC/D,CAAC;AAFD,8CAEC"}
|
|
@@ -73,14 +73,13 @@ function getPeakList(data, options) {
|
|
|
73
73
|
optimization,
|
|
74
74
|
});
|
|
75
75
|
}
|
|
76
|
-
|
|
77
|
-
|
|
76
|
+
return (0, convertWidthToHz_1.convertWidthToHz)(optimize
|
|
77
|
+
? (0, ml_gsd_1.optimizePeaks)(data, peakList, {
|
|
78
78
|
shape,
|
|
79
79
|
factorLimits,
|
|
80
80
|
optimization,
|
|
81
|
-
})
|
|
82
|
-
|
|
83
|
-
return (0, convertWidthToHz_1.convertWidthToHz)(peakList, frequency);
|
|
81
|
+
})
|
|
82
|
+
: peakList, { frequency });
|
|
84
83
|
}
|
|
85
84
|
function getCutOff(data, options) {
|
|
86
85
|
const { noiseLevel, useSanPlot, thresholdFactor } = options;
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"xyAutoPeaksPicking.js","sourceRoot":"","sources":["../../src/xy/xyAutoPeaksPicking.ts"],"names":[],"mappings":";;;AACA,mCAA4D;AAM5D,iEAI+B;AAG/B,qEAAkE;AA4DlE,SAAgB,kBAAkB,CAChC,IAAY,EACZ,OAAkC;IAElC,MAAM,EACJ,IAAI,EACJ,EAAE,EACF,SAAS,EACT,UAAU,EACV,eAAe,GAAG,CAAC,EACnB,WAAW,GAAG,IAAI,EAClB,UAAU,GAAG,OAAO,EACpB,UAAU,GAAG,KAAK,EAClB,OAAO,GAAG,IAAI,EACd,QAAQ,GAAG,KAAK,EAChB,YAAY,GAAG,CAAC,EAChB,gBAAgB,GAAG,IAAI,EACvB,KAAK,GAAG,EAAE,IAAI,EAAE,UAAU,EAAE,EAC5B,YAAY,GAAG,EAAE,IAAI,EAAE,IAAI,EAAE,EAC7B,UAAU,GAAG,IAAI,EACjB,SAAS,GAAG,EAAE,EACd,SAAS,GAAG,UAAU,GACvB,GAAG,OAAO,CAAC;IAEZ,MAAM,EAAE,UAAU,GAAG,CAAC,EAAE,UAAU,GAAG,CAAC,EAAE,GAAG,SAAS,CAAC;IAErD,IAAI,IAAI,CAAC,CAAC,CAAC,MAAM,GAAG,UAAU;QAAE,OAAO,EAAE,CAAC;IAE1C,IAAI,IAAI,KAAK,SAAS,IAAI,EAAE,KAAK,SAAS,EAAE;QAC1C,IAAI,GAAG,IAAA,iCAAS,EAAC,IAAI,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,IAAI,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;KACnD;IAED,MAAM,MAAM,GAAG,SAAS,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE,UAAU,EAAE,UAAU,EAAE,eAAe,EAAE,CAAC,CAAC;IAE9E,IAAI,cAAc,GAAuB;QACvC,KAAK;QACL,SAAS;QACT,UAAU;QACV,QAAQ;QACR,YAAY;QACZ,WAAW,EAAE,IAAI;QACjB,SAAS,EAAE,EAAE,UAAU,EAAE,UAAU,EAAE;QACrC,WAAW;QACX,UAAU;QACV,UAAU,EAAE,MAAM,CAAC,QAAQ;QAC3B,OAAO;QACP,YAAY;QACZ,gBAAgB;KACjB,CAAC;IAEF,QAAQ,SAAS,EAAE;QACjB,KAAK,UAAU,CAAC,CAAC;YACf,OAAO,WAAW,CAAC,IAAI,EAAE,cAAc,CAAC,CAAC;SAC1C;QACD,KAAK,UAAU,CAAC,CAAC;YACf,OAAO,gBAAgB,CAAC,IAAI,EAAE,MAAM,EAAE,cAAc,CAAC,CAAC;SACvD;QACD,KAAK,MAAM,CAAC,CAAC;YACX,OAAO;gBACL,GAAG,WAAW,CAAC,IAAI,EAAE,cAAc,CAAC;gBACpC,GAAG,gBAAgB,CAAC,IAAI,EAAE,MAAM,EAAE,cAAc,CAAC;aAClD,CAAC;SACH;QACD,OAAO,CAAC,CAAC;YACP,MAAM,SAAS,GAAU,SAAS,CAAC;YACnC,4EAA4E;YAC5E,MAAM,KAAK,CAAC,wBAAwB,SAAS,EAAE,CAAC,CAAC;SAClD;KACF;AACH,CAAC;AArED,gDAqEC;AAED,SAAS,gBAAgB,CACvB,IAAY,EACZ,MAA4B,EAC5B,cAAkC;IAElC,cAAc,CAAC,UAAU,GAAG,MAAM,CAAC,QAAQ,CAAC;IAC5C,cAAc,CAAC,WAAW,GAAG,KAAK,CAAC;IACnC,OAAO,WAAW,CAAC,IAAI,EAAE,cAAc,CAAC,CAAC;AAC3C,CAAC;AAED,SAAS,WAAW,CAAC,IAAY,EAAE,OAA2B;IAC5D,MAAM,EACJ,KAAK,EACL,SAAS,EACT,UAAU,EACV,QAAQ,EACR,WAAW,EACX,YAAY,EACZ,SAAS,EACT,WAAW,EACX,UAAU,EACV,UAAU,EACV,OAAO,EACP,YAAY,EACZ,gBAAgB,GACjB,GAAG,OAAO,CAAC;IAEZ,IAAI,QAAQ,GAAG,IAAA,YAAG,EAAC,IAAI,EAAE;QACvB,KAAK;QACL,SAAS;QACT,WAAW;QACX,WAAW;QACX,UAAU;QACV,OAAO;QACP,gBAAgB;KACjB,CAAC,CAAC;IAEH,IAAI,UAAU,EAAE;QACd,QAAQ,GAAG,IAAA,uBAAc,EAAC,QAAQ,EAAE;YAClC,UAAU;YACV,UAAU;YACV,KAAK;YACL,YAAY;SACb,CAAC,CAAC;KACJ;IAED,
|
|
1
|
+
{"version":3,"file":"xyAutoPeaksPicking.js","sourceRoot":"","sources":["../../src/xy/xyAutoPeaksPicking.ts"],"names":[],"mappings":";;;AACA,mCAA4D;AAM5D,iEAI+B;AAG/B,qEAAkE;AA4DlE,SAAgB,kBAAkB,CAChC,IAAY,EACZ,OAAkC;IAElC,MAAM,EACJ,IAAI,EACJ,EAAE,EACF,SAAS,EACT,UAAU,EACV,eAAe,GAAG,CAAC,EACnB,WAAW,GAAG,IAAI,EAClB,UAAU,GAAG,OAAO,EACpB,UAAU,GAAG,KAAK,EAClB,OAAO,GAAG,IAAI,EACd,QAAQ,GAAG,KAAK,EAChB,YAAY,GAAG,CAAC,EAChB,gBAAgB,GAAG,IAAI,EACvB,KAAK,GAAG,EAAE,IAAI,EAAE,UAAU,EAAE,EAC5B,YAAY,GAAG,EAAE,IAAI,EAAE,IAAI,EAAE,EAC7B,UAAU,GAAG,IAAI,EACjB,SAAS,GAAG,EAAE,EACd,SAAS,GAAG,UAAU,GACvB,GAAG,OAAO,CAAC;IAEZ,MAAM,EAAE,UAAU,GAAG,CAAC,EAAE,UAAU,GAAG,CAAC,EAAE,GAAG,SAAS,CAAC;IAErD,IAAI,IAAI,CAAC,CAAC,CAAC,MAAM,GAAG,UAAU;QAAE,OAAO,EAAE,CAAC;IAE1C,IAAI,IAAI,KAAK,SAAS,IAAI,EAAE,KAAK,SAAS,EAAE;QAC1C,IAAI,GAAG,IAAA,iCAAS,EAAC,IAAI,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,IAAI,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;KACnD;IAED,MAAM,MAAM,GAAG,SAAS,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE,UAAU,EAAE,UAAU,EAAE,eAAe,EAAE,CAAC,CAAC;IAE9E,IAAI,cAAc,GAAuB;QACvC,KAAK;QACL,SAAS;QACT,UAAU;QACV,QAAQ;QACR,YAAY;QACZ,WAAW,EAAE,IAAI;QACjB,SAAS,EAAE,EAAE,UAAU,EAAE,UAAU,EAAE;QACrC,WAAW;QACX,UAAU;QACV,UAAU,EAAE,MAAM,CAAC,QAAQ;QAC3B,OAAO;QACP,YAAY;QACZ,gBAAgB;KACjB,CAAC;IAEF,QAAQ,SAAS,EAAE;QACjB,KAAK,UAAU,CAAC,CAAC;YACf,OAAO,WAAW,CAAC,IAAI,EAAE,cAAc,CAAC,CAAC;SAC1C;QACD,KAAK,UAAU,CAAC,CAAC;YACf,OAAO,gBAAgB,CAAC,IAAI,EAAE,MAAM,EAAE,cAAc,CAAC,CAAC;SACvD;QACD,KAAK,MAAM,CAAC,CAAC;YACX,OAAO;gBACL,GAAG,WAAW,CAAC,IAAI,EAAE,cAAc,CAAC;gBACpC,GAAG,gBAAgB,CAAC,IAAI,EAAE,MAAM,EAAE,cAAc,CAAC;aAClD,CAAC;SACH;QACD,OAAO,CAAC,CAAC;YACP,MAAM,SAAS,GAAU,SAAS,CAAC;YACnC,4EAA4E;YAC5E,MAAM,KAAK,CAAC,wBAAwB,SAAS,EAAE,CAAC,CAAC;SAClD;KACF;AACH,CAAC;AArED,gDAqEC;AAED,SAAS,gBAAgB,CACvB,IAAY,EACZ,MAA4B,EAC5B,cAAkC;IAElC,cAAc,CAAC,UAAU,GAAG,MAAM,CAAC,QAAQ,CAAC;IAC5C,cAAc,CAAC,WAAW,GAAG,KAAK,CAAC;IACnC,OAAO,WAAW,CAAC,IAAI,EAAE,cAAc,CAAC,CAAC;AAC3C,CAAC;AAED,SAAS,WAAW,CAAC,IAAY,EAAE,OAA2B;IAC5D,MAAM,EACJ,KAAK,EACL,SAAS,EACT,UAAU,EACV,QAAQ,EACR,WAAW,EACX,YAAY,EACZ,SAAS,EACT,WAAW,EACX,UAAU,EACV,UAAU,EACV,OAAO,EACP,YAAY,EACZ,gBAAgB,GACjB,GAAG,OAAO,CAAC;IAEZ,IAAI,QAAQ,GAAG,IAAA,YAAG,EAAC,IAAI,EAAE;QACvB,KAAK;QACL,SAAS;QACT,WAAW;QACX,WAAW;QACX,UAAU;QACV,OAAO;QACP,gBAAgB;KACjB,CAAC,CAAC;IAEH,IAAI,UAAU,EAAE;QACd,QAAQ,GAAG,IAAA,uBAAc,EAAC,QAAQ,EAAE;YAClC,UAAU;YACV,UAAU;YACV,KAAK;YACL,YAAY;SACb,CAAC,CAAC;KACJ;IAED,OAAO,IAAA,mCAAgB,EACrB,QAAQ;QACN,CAAC,CAAC,IAAA,sBAAa,EAAC,IAAI,EAAE,QAAQ,EAAE;YAC5B,KAAK;YACL,YAAY;YACZ,YAAY;SACb,CAAC;QACJ,CAAC,CAAC,QAAQ,EACZ,EAAE,SAAS,EAAE,CACd,CAAC;AACJ,CAAC;AAED,SAAS,SAAS,CAAC,IAA6B,EAAE,OAAyB;IACzE,MAAM,EAAE,UAAU,EAAE,UAAU,EAAE,eAAe,EAAE,GAAG,OAAO,CAAC;IAE5D,MAAM,YAAY,GAAG,CAAC,UAAkB,EAAE,EAAE,CAC1C,OAAO,UAAU,KAAK,QAAQ;QAC5B,CAAC,CAAC,EAAE,QAAQ,EAAE,UAAU,EAAE,QAAQ,EAAE,CAAC,UAAU,EAAE;QACjD,CAAC,CAAC,UAAU,CAAC;IAEjB,IAAI,UAAU,EAAE;QACd,OAAO,YAAY,CAAC,UAAU,CAAC,CAAC;KACjC;SAAM;QACL,OAAO,UAAU;YACf,CAAC,CAAC,IAAA,qCAAa,EAAC,IAAI,EAAE,EAAE,SAAS,EAAE,eAAe,EAAE,CAAC;YACrD,CAAC,CAAC,YAAY,CAAC,IAAA,uCAAe,EAAC,IAAI,CAAC,GAAG,eAAe,CAAC,CAAC;KAC3D;AACH,CAAC"}
|
|
@@ -1,3 +1,9 @@
|
|
|
1
1
|
import { DataXY, PeakXYWidth } from 'cheminfo-types';
|
|
2
2
|
import { OptimizePeaksOptions } from 'ml-gsd';
|
|
3
|
-
export
|
|
3
|
+
export interface XYPeaksOptimizationOptions extends OptimizePeaksOptions {
|
|
4
|
+
/**
|
|
5
|
+
* frequency to convert width and fwhm to ppm.
|
|
6
|
+
*/
|
|
7
|
+
frequency: number;
|
|
8
|
+
}
|
|
9
|
+
export declare function xyPeaksOptimization(data: DataXY, peaks: PeakXYWidth[], options: XYPeaksOptimizationOptions): any;
|
|
@@ -2,8 +2,16 @@
|
|
|
2
2
|
Object.defineProperty(exports, "__esModule", { value: true });
|
|
3
3
|
exports.xyPeaksOptimization = void 0;
|
|
4
4
|
const ml_gsd_1 = require("ml-gsd");
|
|
5
|
+
const convertWidthToHz_1 = require("../peaks/util/convertWidthToHz");
|
|
6
|
+
const convertWidthToPPM_1 = require("../peaks/util/convertWidthToPPM");
|
|
5
7
|
function xyPeaksOptimization(data, peaks, options) {
|
|
6
|
-
|
|
8
|
+
const { frequency } = options;
|
|
9
|
+
const newPeaks = (0, convertWidthToPPM_1.convertWidthToPPM)(peaks, { frequency });
|
|
10
|
+
const optimizedPeaks = (0, ml_gsd_1.optimizePeaks)(data, newPeaks, options);
|
|
11
|
+
return (0, convertWidthToHz_1.convertWidthToHz)(optimizedPeaks, {
|
|
12
|
+
frequency,
|
|
13
|
+
output: optimizedPeaks,
|
|
14
|
+
});
|
|
7
15
|
}
|
|
8
16
|
exports.xyPeaksOptimization = xyPeaksOptimization;
|
|
9
17
|
//# sourceMappingURL=xyPeaksOptimization.js.map
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"xyPeaksOptimization.js","sourceRoot":"","sources":["../../src/xy/xyPeaksOptimization.ts"],"names":[],"mappings":";;;AACA,mCAA6D;AAE7D,SAAgB,mBAAmB,CACjC,IAAY,EACZ,KAAoB,EACpB,
|
|
1
|
+
{"version":3,"file":"xyPeaksOptimization.js","sourceRoot":"","sources":["../../src/xy/xyPeaksOptimization.ts"],"names":[],"mappings":";;;AACA,mCAA6D;AAE7D,qEAAkE;AAClE,uEAAoE;AAQpE,SAAgB,mBAAmB,CACjC,IAAY,EACZ,KAAoB,EACpB,OAAmC;IAEnC,MAAM,EAAE,SAAS,EAAE,GAAG,OAAO,CAAC;IAC9B,MAAM,QAAQ,GAAG,IAAA,qCAAiB,EAAC,KAAK,EAAE,EAAE,SAAS,EAAE,CAAC,CAAC;IACzD,MAAM,cAAc,GAAG,IAAA,sBAAa,EAAC,IAAI,EAAE,QAAQ,EAAE,OAAO,CAAC,CAAC;IAC9D,OAAO,IAAA,mCAAgB,EAAC,cAAc,EAAE;QACtC,SAAS;QACT,MAAM,EAAE,cAAc;KACvB,CAAC,CAAC;AACL,CAAC;AAZD,kDAYC"}
|
|
@@ -1,14 +1,30 @@
|
|
|
1
1
|
import { generateSpectrum } from 'spectrum-generator';
|
|
2
2
|
import { convertWidthToPPM } from './util/convertWidthToPPM';
|
|
3
3
|
export function peaksToXY(peaks, options) {
|
|
4
|
-
const {
|
|
5
|
-
|
|
4
|
+
const { frequency, nbPoints = 1024, shape } = options;
|
|
5
|
+
const newPeaks = convertWidthToPPM(peaks, { frequency });
|
|
6
|
+
return generateSpectrum(newPeaks, {
|
|
6
7
|
generator: {
|
|
7
|
-
|
|
8
|
-
to,
|
|
8
|
+
...getFromTo(newPeaks, options),
|
|
9
9
|
nbPoints,
|
|
10
10
|
shape,
|
|
11
11
|
},
|
|
12
12
|
});
|
|
13
13
|
}
|
|
14
|
+
function getFromTo(newPeaks, options) {
|
|
15
|
+
if ('to' in options && 'from' in options) {
|
|
16
|
+
return {
|
|
17
|
+
from: options.from,
|
|
18
|
+
to: options.to,
|
|
19
|
+
};
|
|
20
|
+
}
|
|
21
|
+
newPeaks.sort((a, b) => a.x - b.x);
|
|
22
|
+
const firstPeak = newPeaks[0];
|
|
23
|
+
const lastPeak = newPeaks[newPeaks.length - 1];
|
|
24
|
+
const { from = firstPeak.x - firstPeak.width * 2, to = lastPeak.x + lastPeak.width * 2, } = options;
|
|
25
|
+
return {
|
|
26
|
+
from,
|
|
27
|
+
to,
|
|
28
|
+
};
|
|
29
|
+
}
|
|
14
30
|
//# sourceMappingURL=peaksToXY.js.map
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"peaksToXY.js","sourceRoot":"","sources":["../../src/peaks/peaksToXY.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,gBAAgB,EAAe,MAAM,oBAAoB,CAAC;AAGnE,OAAO,EAAE,iBAAiB,EAAE,MAAM,0BAA0B,CAAC;
|
|
1
|
+
{"version":3,"file":"peaksToXY.js","sourceRoot":"","sources":["../../src/peaks/peaksToXY.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,gBAAgB,EAAe,MAAM,oBAAoB,CAAC;AAGnE,OAAO,EAAE,iBAAiB,EAAE,MAAM,0BAA0B,CAAC;AAuB7D,MAAM,UAAU,SAAS,CAAC,KAAkB,EAAE,OAAyB;IACrE,MAAM,EAAE,SAAS,EAAE,QAAQ,GAAG,IAAI,EAAE,KAAK,EAAE,GAAG,OAAO,CAAC;IACtD,MAAM,QAAQ,GAAG,iBAAiB,CAAC,KAAK,EAAE,EAAE,SAAS,EAAE,CAAC,CAAC;IAEzD,OAAO,gBAAgB,CAAC,QAAQ,EAAE;QAChC,SAAS,EAAE;YACT,GAAG,SAAS,CAAC,QAAQ,EAAE,OAAO,CAAC;YAC/B,QAAQ;YACR,KAAK;SACN;KACF,CAAC,CAAC;AACL,CAAC;AAED,SAAS,SAAS,CAAC,QAAqB,EAAE,OAAyB;IACjE,IAAI,IAAI,IAAI,OAAO,IAAI,MAAM,IAAI,OAAO,EAAE;QACxC,OAAO;YACL,IAAI,EAAE,OAAO,CAAC,IAAI;YAClB,EAAE,EAAE,OAAO,CAAC,EAAE;SACf,CAAC;KACH;IAED,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACnC,MAAM,SAAS,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;IAC9B,MAAM,QAAQ,GAAG,QAAQ,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IAE/C,MAAM,EACJ,IAAI,GAAG,SAAS,CAAC,CAAC,GAAG,SAAS,CAAC,KAAK,GAAG,CAAC,EACxC,EAAE,GAAG,QAAQ,CAAC,CAAC,GAAG,QAAQ,CAAC,KAAK,GAAG,CAAC,GACrC,GAAG,OAAO,CAAC;IAEZ,OAAO;QACL,IAAI;QACJ,EAAE;KACH,CAAC;AACJ,CAAC"}
|
|
@@ -2,17 +2,17 @@
|
|
|
2
2
|
* convert width and fwhm to ppm
|
|
3
3
|
*/
|
|
4
4
|
export function convertWidth(peaks, options) {
|
|
5
|
-
const { frequency, convertTo } = options;
|
|
5
|
+
const { frequency, convertTo, output = JSON.parse(JSON.stringify(peaks)), } = options;
|
|
6
6
|
const convert = getConverter(convertTo, frequency);
|
|
7
|
-
|
|
7
|
+
for (const peak of output) {
|
|
8
8
|
peak.width = convert(peak.width);
|
|
9
9
|
if (hasShape(peak)) {
|
|
10
10
|
if (peak.shape.fwhm) {
|
|
11
11
|
peak.shape.fwhm = convert(peak.shape.fwhm);
|
|
12
12
|
}
|
|
13
13
|
}
|
|
14
|
-
|
|
15
|
-
|
|
14
|
+
}
|
|
15
|
+
return output;
|
|
16
16
|
}
|
|
17
17
|
function hasShape(peak) {
|
|
18
18
|
return 'shape' in peak;
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"convertWidth.js","sourceRoot":"","sources":["../../../src/peaks/util/convertWidth.ts"],"names":[],"mappings":"
|
|
1
|
+
{"version":3,"file":"convertWidth.js","sourceRoot":"","sources":["../../../src/peaks/util/convertWidth.ts"],"names":[],"mappings":"AAUA;;GAEG;AACH,MAAM,UAAU,YAAY,CAAC,KAAkB,EAAE,OAA4B;IAC3E,MAAM,EACJ,SAAS,EACT,SAAS,EACT,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC,GAC3C,GAAG,OAAO,CAAC;IAEZ,MAAM,OAAO,GAAG,YAAY,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;IAEnD,KAAK,MAAM,IAAI,IAAI,MAAM,EAAE;QACzB,IAAI,CAAC,KAAK,GAAG,OAAO,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACjC,IAAI,QAAQ,CAAC,IAAI,CAAC,EAAE;YAClB,IAAI,IAAI,CAAC,KAAK,CAAC,IAAI,EAAE;gBACnB,IAAI,CAAC,KAAK,CAAC,IAAI,GAAG,OAAO,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;aAC5C;SACF;KACF;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,SAAS,QAAQ,CAAC,IAAe;IAC/B,OAAO,OAAO,IAAI,IAAI,CAAC;AACzB,CAAC;AAED,SAAS,YAAY,CAAC,SAAiB,EAAE,SAAiB;IACxD,QAAQ,SAAS,EAAE;QACjB,KAAK,KAAK;YACR,OAAO,CAAC,CAAS,EAAE,EAAE,CAAC,CAAC,GAAG,SAAS,CAAC;QACtC,KAAK,IAAI;YACP,OAAO,CAAC,CAAS,EAAE,EAAE,CAAC,CAAC,GAAG,SAAS,CAAC;QACtC;YACE,MAAM,IAAI,KAAK,CAAC,+BAA+B,SAAS,EAAE,CAAC,CAAC;KAC/D;AACH,CAAC"}
|
|
@@ -1,5 +1,5 @@
|
|
|
1
1
|
import { convertWidth } from './convertWidth';
|
|
2
|
-
export function convertWidthToHz(peaks,
|
|
3
|
-
return convertWidth(peaks, {
|
|
2
|
+
export function convertWidthToHz(peaks, options) {
|
|
3
|
+
return convertWidth(peaks, { ...options, convertTo: 'hz' });
|
|
4
4
|
}
|
|
5
5
|
//# sourceMappingURL=convertWidthToHz.js.map
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"convertWidthToHz.js","sourceRoot":"","sources":["../../../src/peaks/util/convertWidthToHz.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,YAAY,
|
|
1
|
+
{"version":3,"file":"convertWidthToHz.js","sourceRoot":"","sources":["../../../src/peaks/util/convertWidthToHz.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,YAAY,EAAuB,MAAM,gBAAgB,CAAC;AAInE,MAAM,UAAU,gBAAgB,CAAC,KAAkB,EAAE,OAAuB;IAC1E,OAAO,YAAY,CAAC,KAAK,EAAE,EAAE,GAAG,OAAO,EAAE,SAAS,EAAE,IAAI,EAAE,CAAC,CAAC;AAC9D,CAAC"}
|
|
@@ -1,5 +1,5 @@
|
|
|
1
1
|
import { convertWidth } from './convertWidth';
|
|
2
|
-
export function convertWidthToPPM(peaks,
|
|
3
|
-
return convertWidth(peaks, {
|
|
2
|
+
export function convertWidthToPPM(peaks, options) {
|
|
3
|
+
return convertWidth(peaks, { ...options, convertTo: 'ppm' });
|
|
4
4
|
}
|
|
5
5
|
//# sourceMappingURL=convertWidthToPPM.js.map
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"convertWidthToPPM.js","sourceRoot":"","sources":["../../../src/peaks/util/convertWidthToPPM.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,YAAY,EAAE,MAAM,gBAAgB,CAAC;
|
|
1
|
+
{"version":3,"file":"convertWidthToPPM.js","sourceRoot":"","sources":["../../../src/peaks/util/convertWidthToPPM.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,YAAY,EAAE,MAAM,gBAAgB,CAAC;AAG9C,MAAM,UAAU,iBAAiB,CAAC,KAAkB,EAAE,OAAuB;IAC3E,OAAO,YAAY,CAAC,KAAK,EAAE,EAAE,GAAG,OAAO,EAAE,SAAS,EAAE,KAAK,EAAE,CAAC,CAAC;AAC/D,CAAC"}
|
|
@@ -69,14 +69,13 @@ function getPeakList(data, options) {
|
|
|
69
69
|
optimization,
|
|
70
70
|
});
|
|
71
71
|
}
|
|
72
|
-
|
|
73
|
-
|
|
72
|
+
return convertWidthToHz(optimize
|
|
73
|
+
? optimizePeaks(data, peakList, {
|
|
74
74
|
shape,
|
|
75
75
|
factorLimits,
|
|
76
76
|
optimization,
|
|
77
|
-
})
|
|
78
|
-
|
|
79
|
-
return convertWidthToHz(peakList, frequency);
|
|
77
|
+
})
|
|
78
|
+
: peakList, { frequency });
|
|
80
79
|
}
|
|
81
80
|
function getCutOff(data, options) {
|
|
82
81
|
const { noiseLevel, useSanPlot, thresholdFactor } = options;
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"xyAutoPeaksPicking.js","sourceRoot":"","sources":["../../src/xy/xyAutoPeaksPicking.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,GAAG,EAAE,aAAa,EAAE,cAAc,EAAE,MAAM,QAAQ,CAAC;AAM5D,OAAO,EACL,SAAS,EACT,aAAa,EACb,eAAe,GAChB,MAAM,uBAAuB,CAAC;AAG/B,OAAO,EAAE,gBAAgB,EAAE,MAAM,gCAAgC,CAAC;AA4DlE,MAAM,UAAU,kBAAkB,CAChC,IAAY,EACZ,OAAkC;IAElC,MAAM,EACJ,IAAI,EACJ,EAAE,EACF,SAAS,EACT,UAAU,EACV,eAAe,GAAG,CAAC,EACnB,WAAW,GAAG,IAAI,EAClB,UAAU,GAAG,OAAO,EACpB,UAAU,GAAG,KAAK,EAClB,OAAO,GAAG,IAAI,EACd,QAAQ,GAAG,KAAK,EAChB,YAAY,GAAG,CAAC,EAChB,gBAAgB,GAAG,IAAI,EACvB,KAAK,GAAG,EAAE,IAAI,EAAE,UAAU,EAAE,EAC5B,YAAY,GAAG,EAAE,IAAI,EAAE,IAAI,EAAE,EAC7B,UAAU,GAAG,IAAI,EACjB,SAAS,GAAG,EAAE,EACd,SAAS,GAAG,UAAU,GACvB,GAAG,OAAO,CAAC;IAEZ,MAAM,EAAE,UAAU,GAAG,CAAC,EAAE,UAAU,GAAG,CAAC,EAAE,GAAG,SAAS,CAAC;IAErD,IAAI,IAAI,CAAC,CAAC,CAAC,MAAM,GAAG,UAAU;QAAE,OAAO,EAAE,CAAC;IAE1C,IAAI,IAAI,KAAK,SAAS,IAAI,EAAE,KAAK,SAAS,EAAE;QAC1C,IAAI,GAAG,SAAS,CAAC,IAAI,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,IAAI,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;KACnD;IAED,MAAM,MAAM,GAAG,SAAS,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE,UAAU,EAAE,UAAU,EAAE,eAAe,EAAE,CAAC,CAAC;IAE9E,IAAI,cAAc,GAAuB;QACvC,KAAK;QACL,SAAS;QACT,UAAU;QACV,QAAQ;QACR,YAAY;QACZ,WAAW,EAAE,IAAI;QACjB,SAAS,EAAE,EAAE,UAAU,EAAE,UAAU,EAAE;QACrC,WAAW;QACX,UAAU;QACV,UAAU,EAAE,MAAM,CAAC,QAAQ;QAC3B,OAAO;QACP,YAAY;QACZ,gBAAgB;KACjB,CAAC;IAEF,QAAQ,SAAS,EAAE;QACjB,KAAK,UAAU,CAAC,CAAC;YACf,OAAO,WAAW,CAAC,IAAI,EAAE,cAAc,CAAC,CAAC;SAC1C;QACD,KAAK,UAAU,CAAC,CAAC;YACf,OAAO,gBAAgB,CAAC,IAAI,EAAE,MAAM,EAAE,cAAc,CAAC,CAAC;SACvD;QACD,KAAK,MAAM,CAAC,CAAC;YACX,OAAO;gBACL,GAAG,WAAW,CAAC,IAAI,EAAE,cAAc,CAAC;gBACpC,GAAG,gBAAgB,CAAC,IAAI,EAAE,MAAM,EAAE,cAAc,CAAC;aAClD,CAAC;SACH;QACD,OAAO,CAAC,CAAC;YACP,MAAM,SAAS,GAAU,SAAS,CAAC;YACnC,4EAA4E;YAC5E,MAAM,KAAK,CAAC,wBAAwB,SAAS,EAAE,CAAC,CAAC;SAClD;KACF;AACH,CAAC;AAED,SAAS,gBAAgB,CACvB,IAAY,EACZ,MAA4B,EAC5B,cAAkC;IAElC,cAAc,CAAC,UAAU,GAAG,MAAM,CAAC,QAAQ,CAAC;IAC5C,cAAc,CAAC,WAAW,GAAG,KAAK,CAAC;IACnC,OAAO,WAAW,CAAC,IAAI,EAAE,cAAc,CAAC,CAAC;AAC3C,CAAC;AAED,SAAS,WAAW,CAAC,IAAY,EAAE,OAA2B;IAC5D,MAAM,EACJ,KAAK,EACL,SAAS,EACT,UAAU,EACV,QAAQ,EACR,WAAW,EACX,YAAY,EACZ,SAAS,EACT,WAAW,EACX,UAAU,EACV,UAAU,EACV,OAAO,EACP,YAAY,EACZ,gBAAgB,GACjB,GAAG,OAAO,CAAC;IAEZ,IAAI,QAAQ,GAAG,GAAG,CAAC,IAAI,EAAE;QACvB,KAAK;QACL,SAAS;QACT,WAAW;QACX,WAAW;QACX,UAAU;QACV,OAAO;QACP,gBAAgB;KACjB,CAAC,CAAC;IAEH,IAAI,UAAU,EAAE;QACd,QAAQ,GAAG,cAAc,CAAC,QAAQ,EAAE;YAClC,UAAU;YACV,UAAU;YACV,KAAK;YACL,YAAY;SACb,CAAC,CAAC;KACJ;IAED,
|
|
1
|
+
{"version":3,"file":"xyAutoPeaksPicking.js","sourceRoot":"","sources":["../../src/xy/xyAutoPeaksPicking.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,GAAG,EAAE,aAAa,EAAE,cAAc,EAAE,MAAM,QAAQ,CAAC;AAM5D,OAAO,EACL,SAAS,EACT,aAAa,EACb,eAAe,GAChB,MAAM,uBAAuB,CAAC;AAG/B,OAAO,EAAE,gBAAgB,EAAE,MAAM,gCAAgC,CAAC;AA4DlE,MAAM,UAAU,kBAAkB,CAChC,IAAY,EACZ,OAAkC;IAElC,MAAM,EACJ,IAAI,EACJ,EAAE,EACF,SAAS,EACT,UAAU,EACV,eAAe,GAAG,CAAC,EACnB,WAAW,GAAG,IAAI,EAClB,UAAU,GAAG,OAAO,EACpB,UAAU,GAAG,KAAK,EAClB,OAAO,GAAG,IAAI,EACd,QAAQ,GAAG,KAAK,EAChB,YAAY,GAAG,CAAC,EAChB,gBAAgB,GAAG,IAAI,EACvB,KAAK,GAAG,EAAE,IAAI,EAAE,UAAU,EAAE,EAC5B,YAAY,GAAG,EAAE,IAAI,EAAE,IAAI,EAAE,EAC7B,UAAU,GAAG,IAAI,EACjB,SAAS,GAAG,EAAE,EACd,SAAS,GAAG,UAAU,GACvB,GAAG,OAAO,CAAC;IAEZ,MAAM,EAAE,UAAU,GAAG,CAAC,EAAE,UAAU,GAAG,CAAC,EAAE,GAAG,SAAS,CAAC;IAErD,IAAI,IAAI,CAAC,CAAC,CAAC,MAAM,GAAG,UAAU;QAAE,OAAO,EAAE,CAAC;IAE1C,IAAI,IAAI,KAAK,SAAS,IAAI,EAAE,KAAK,SAAS,EAAE;QAC1C,IAAI,GAAG,SAAS,CAAC,IAAI,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,IAAI,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;KACnD;IAED,MAAM,MAAM,GAAG,SAAS,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE,UAAU,EAAE,UAAU,EAAE,eAAe,EAAE,CAAC,CAAC;IAE9E,IAAI,cAAc,GAAuB;QACvC,KAAK;QACL,SAAS;QACT,UAAU;QACV,QAAQ;QACR,YAAY;QACZ,WAAW,EAAE,IAAI;QACjB,SAAS,EAAE,EAAE,UAAU,EAAE,UAAU,EAAE;QACrC,WAAW;QACX,UAAU;QACV,UAAU,EAAE,MAAM,CAAC,QAAQ;QAC3B,OAAO;QACP,YAAY;QACZ,gBAAgB;KACjB,CAAC;IAEF,QAAQ,SAAS,EAAE;QACjB,KAAK,UAAU,CAAC,CAAC;YACf,OAAO,WAAW,CAAC,IAAI,EAAE,cAAc,CAAC,CAAC;SAC1C;QACD,KAAK,UAAU,CAAC,CAAC;YACf,OAAO,gBAAgB,CAAC,IAAI,EAAE,MAAM,EAAE,cAAc,CAAC,CAAC;SACvD;QACD,KAAK,MAAM,CAAC,CAAC;YACX,OAAO;gBACL,GAAG,WAAW,CAAC,IAAI,EAAE,cAAc,CAAC;gBACpC,GAAG,gBAAgB,CAAC,IAAI,EAAE,MAAM,EAAE,cAAc,CAAC;aAClD,CAAC;SACH;QACD,OAAO,CAAC,CAAC;YACP,MAAM,SAAS,GAAU,SAAS,CAAC;YACnC,4EAA4E;YAC5E,MAAM,KAAK,CAAC,wBAAwB,SAAS,EAAE,CAAC,CAAC;SAClD;KACF;AACH,CAAC;AAED,SAAS,gBAAgB,CACvB,IAAY,EACZ,MAA4B,EAC5B,cAAkC;IAElC,cAAc,CAAC,UAAU,GAAG,MAAM,CAAC,QAAQ,CAAC;IAC5C,cAAc,CAAC,WAAW,GAAG,KAAK,CAAC;IACnC,OAAO,WAAW,CAAC,IAAI,EAAE,cAAc,CAAC,CAAC;AAC3C,CAAC;AAED,SAAS,WAAW,CAAC,IAAY,EAAE,OAA2B;IAC5D,MAAM,EACJ,KAAK,EACL,SAAS,EACT,UAAU,EACV,QAAQ,EACR,WAAW,EACX,YAAY,EACZ,SAAS,EACT,WAAW,EACX,UAAU,EACV,UAAU,EACV,OAAO,EACP,YAAY,EACZ,gBAAgB,GACjB,GAAG,OAAO,CAAC;IAEZ,IAAI,QAAQ,GAAG,GAAG,CAAC,IAAI,EAAE;QACvB,KAAK;QACL,SAAS;QACT,WAAW;QACX,WAAW;QACX,UAAU;QACV,OAAO;QACP,gBAAgB;KACjB,CAAC,CAAC;IAEH,IAAI,UAAU,EAAE;QACd,QAAQ,GAAG,cAAc,CAAC,QAAQ,EAAE;YAClC,UAAU;YACV,UAAU;YACV,KAAK;YACL,YAAY;SACb,CAAC,CAAC;KACJ;IAED,OAAO,gBAAgB,CACrB,QAAQ;QACN,CAAC,CAAC,aAAa,CAAC,IAAI,EAAE,QAAQ,EAAE;YAC5B,KAAK;YACL,YAAY;YACZ,YAAY;SACb,CAAC;QACJ,CAAC,CAAC,QAAQ,EACZ,EAAE,SAAS,EAAE,CACd,CAAC;AACJ,CAAC;AAED,SAAS,SAAS,CAAC,IAA6B,EAAE,OAAyB;IACzE,MAAM,EAAE,UAAU,EAAE,UAAU,EAAE,eAAe,EAAE,GAAG,OAAO,CAAC;IAE5D,MAAM,YAAY,GAAG,CAAC,UAAkB,EAAE,EAAE,CAC1C,OAAO,UAAU,KAAK,QAAQ;QAC5B,CAAC,CAAC,EAAE,QAAQ,EAAE,UAAU,EAAE,QAAQ,EAAE,CAAC,UAAU,EAAE;QACjD,CAAC,CAAC,UAAU,CAAC;IAEjB,IAAI,UAAU,EAAE;QACd,OAAO,YAAY,CAAC,UAAU,CAAC,CAAC;KACjC;SAAM;QACL,OAAO,UAAU;YACf,CAAC,CAAC,aAAa,CAAC,IAAI,EAAE,EAAE,SAAS,EAAE,eAAe,EAAE,CAAC;YACrD,CAAC,CAAC,YAAY,CAAC,eAAe,CAAC,IAAI,CAAC,GAAG,eAAe,CAAC,CAAC;KAC3D;AACH,CAAC"}
|
|
@@ -1,5 +1,13 @@
|
|
|
1
1
|
import { optimizePeaks } from 'ml-gsd';
|
|
2
|
+
import { convertWidthToHz } from '../peaks/util/convertWidthToHz';
|
|
3
|
+
import { convertWidthToPPM } from '../peaks/util/convertWidthToPPM';
|
|
2
4
|
export function xyPeaksOptimization(data, peaks, options) {
|
|
3
|
-
|
|
5
|
+
const { frequency } = options;
|
|
6
|
+
const newPeaks = convertWidthToPPM(peaks, { frequency });
|
|
7
|
+
const optimizedPeaks = optimizePeaks(data, newPeaks, options);
|
|
8
|
+
return convertWidthToHz(optimizedPeaks, {
|
|
9
|
+
frequency,
|
|
10
|
+
output: optimizedPeaks,
|
|
11
|
+
});
|
|
4
12
|
}
|
|
5
13
|
//# sourceMappingURL=xyPeaksOptimization.js.map
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"xyPeaksOptimization.js","sourceRoot":"","sources":["../../src/xy/xyPeaksOptimization.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,aAAa,EAAwB,MAAM,QAAQ,CAAC;AAE7D,MAAM,UAAU,mBAAmB,CACjC,IAAY,EACZ,KAAoB,EACpB,
|
|
1
|
+
{"version":3,"file":"xyPeaksOptimization.js","sourceRoot":"","sources":["../../src/xy/xyPeaksOptimization.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,aAAa,EAAwB,MAAM,QAAQ,CAAC;AAE7D,OAAO,EAAE,gBAAgB,EAAE,MAAM,gCAAgC,CAAC;AAClE,OAAO,EAAE,iBAAiB,EAAE,MAAM,iCAAiC,CAAC;AAQpE,MAAM,UAAU,mBAAmB,CACjC,IAAY,EACZ,KAAoB,EACpB,OAAmC;IAEnC,MAAM,EAAE,SAAS,EAAE,GAAG,OAAO,CAAC;IAC9B,MAAM,QAAQ,GAAG,iBAAiB,CAAC,KAAK,EAAE,EAAE,SAAS,EAAE,CAAC,CAAC;IACzD,MAAM,cAAc,GAAG,aAAa,CAAC,IAAI,EAAE,QAAQ,EAAE,OAAO,CAAC,CAAC;IAC9D,OAAO,gBAAgB,CAAC,cAAc,EAAE;QACtC,SAAS;QACT,MAAM,EAAE,cAAc;KACvB,CAAC,CAAC;AACL,CAAC"}
|
package/package.json
CHANGED
package/src/peaks/peaksToXY.ts
CHANGED
|
@@ -8,6 +8,16 @@ import { convertWidthToPPM } from './util/convertWidthToPPM';
|
|
|
8
8
|
*/
|
|
9
9
|
|
|
10
10
|
export interface PeaksToXYOptions extends OptionsSG1D {
|
|
11
|
+
/**
|
|
12
|
+
* First x value (inclusive).
|
|
13
|
+
* @default `firstPeak.x + firstPeak.width * 3`
|
|
14
|
+
*/
|
|
15
|
+
from?: number;
|
|
16
|
+
/**
|
|
17
|
+
* Last x value (inclusive).
|
|
18
|
+
* @default `lastPeak.x + lastPeak.width * 3`
|
|
19
|
+
*/
|
|
20
|
+
to?: number;
|
|
11
21
|
/**
|
|
12
22
|
* frequency of the spectrum
|
|
13
23
|
*/
|
|
@@ -15,14 +25,37 @@ export interface PeaksToXYOptions extends OptionsSG1D {
|
|
|
15
25
|
}
|
|
16
26
|
|
|
17
27
|
export function peaksToXY(peaks: NMRPeak1D[], options: PeaksToXYOptions) {
|
|
18
|
-
const {
|
|
28
|
+
const { frequency, nbPoints = 1024, shape } = options;
|
|
29
|
+
const newPeaks = convertWidthToPPM(peaks, { frequency });
|
|
19
30
|
|
|
20
|
-
return generateSpectrum(
|
|
31
|
+
return generateSpectrum(newPeaks, {
|
|
21
32
|
generator: {
|
|
22
|
-
|
|
23
|
-
to,
|
|
33
|
+
...getFromTo(newPeaks, options),
|
|
24
34
|
nbPoints,
|
|
25
35
|
shape,
|
|
26
36
|
},
|
|
27
37
|
});
|
|
28
38
|
}
|
|
39
|
+
|
|
40
|
+
function getFromTo(newPeaks: NMRPeak1D[], options: PeaksToXYOptions) {
|
|
41
|
+
if ('to' in options && 'from' in options) {
|
|
42
|
+
return {
|
|
43
|
+
from: options.from,
|
|
44
|
+
to: options.to,
|
|
45
|
+
};
|
|
46
|
+
}
|
|
47
|
+
|
|
48
|
+
newPeaks.sort((a, b) => a.x - b.x);
|
|
49
|
+
const firstPeak = newPeaks[0];
|
|
50
|
+
const lastPeak = newPeaks[newPeaks.length - 1];
|
|
51
|
+
|
|
52
|
+
const {
|
|
53
|
+
from = firstPeak.x - firstPeak.width * 2,
|
|
54
|
+
to = lastPeak.x + lastPeak.width * 2,
|
|
55
|
+
} = options;
|
|
56
|
+
|
|
57
|
+
return {
|
|
58
|
+
from,
|
|
59
|
+
to,
|
|
60
|
+
};
|
|
61
|
+
}
|
|
@@ -6,22 +6,29 @@ export type NMRPeak1DWithShape = MakeMandatory<NMRPeak1D, 'shape'>;
|
|
|
6
6
|
export interface ConvertWidthOptions {
|
|
7
7
|
convertTo: string & { _kind?: 'hz' | 'ppm' };
|
|
8
8
|
frequency: number;
|
|
9
|
+
output?: NMRPeak1D[];
|
|
9
10
|
}
|
|
10
11
|
/**
|
|
11
12
|
* convert width and fwhm to ppm
|
|
12
13
|
*/
|
|
13
14
|
export function convertWidth(peaks: NMRPeak1D[], options: ConvertWidthOptions) {
|
|
14
|
-
const {
|
|
15
|
+
const {
|
|
16
|
+
frequency,
|
|
17
|
+
convertTo,
|
|
18
|
+
output = JSON.parse(JSON.stringify(peaks)),
|
|
19
|
+
} = options;
|
|
20
|
+
|
|
15
21
|
const convert = getConverter(convertTo, frequency);
|
|
16
|
-
|
|
22
|
+
|
|
23
|
+
for (const peak of output) {
|
|
17
24
|
peak.width = convert(peak.width);
|
|
18
25
|
if (hasShape(peak)) {
|
|
19
26
|
if (peak.shape.fwhm) {
|
|
20
27
|
peak.shape.fwhm = convert(peak.shape.fwhm);
|
|
21
28
|
}
|
|
22
29
|
}
|
|
23
|
-
|
|
24
|
-
|
|
30
|
+
}
|
|
31
|
+
return output;
|
|
25
32
|
}
|
|
26
33
|
|
|
27
34
|
function hasShape(peak: NMRPeak1D): peak is NMRPeak1DWithShape {
|
|
@@ -1,7 +1,9 @@
|
|
|
1
1
|
import { NMRPeak1D } from '../NMRPeak1D';
|
|
2
2
|
|
|
3
|
-
import { convertWidth } from './convertWidth';
|
|
3
|
+
import { convertWidth, ConvertWidthOptions } from './convertWidth';
|
|
4
4
|
|
|
5
|
-
export
|
|
6
|
-
|
|
5
|
+
export type ConvertWidthTo = Omit<ConvertWidthOptions, 'convertTo'>;
|
|
6
|
+
|
|
7
|
+
export function convertWidthToHz(peaks: NMRPeak1D[], options: ConvertWidthTo) {
|
|
8
|
+
return convertWidth(peaks, { ...options, convertTo: 'hz' });
|
|
7
9
|
}
|
|
@@ -1,7 +1,8 @@
|
|
|
1
1
|
import { NMRPeak1D } from '../NMRPeak1D';
|
|
2
2
|
|
|
3
3
|
import { convertWidth } from './convertWidth';
|
|
4
|
+
import { ConvertWidthTo } from './convertWidthToHz';
|
|
4
5
|
|
|
5
|
-
export function convertWidthToPPM(peaks: NMRPeak1D[],
|
|
6
|
-
return convertWidth(peaks, {
|
|
6
|
+
export function convertWidthToPPM(peaks: NMRPeak1D[], options: ConvertWidthTo) {
|
|
7
|
+
return convertWidth(peaks, { ...options, convertTo: 'ppm' });
|
|
7
8
|
}
|
|
@@ -189,15 +189,16 @@ function getPeakList(data: DataXY, options: GetPeakListOptions) {
|
|
|
189
189
|
});
|
|
190
190
|
}
|
|
191
191
|
|
|
192
|
-
|
|
193
|
-
|
|
194
|
-
|
|
195
|
-
|
|
196
|
-
|
|
197
|
-
|
|
198
|
-
|
|
199
|
-
|
|
200
|
-
|
|
192
|
+
return convertWidthToHz(
|
|
193
|
+
optimize
|
|
194
|
+
? optimizePeaks(data, peakList, {
|
|
195
|
+
shape,
|
|
196
|
+
factorLimits,
|
|
197
|
+
optimization,
|
|
198
|
+
})
|
|
199
|
+
: peakList,
|
|
200
|
+
{ frequency },
|
|
201
|
+
);
|
|
201
202
|
}
|
|
202
203
|
|
|
203
204
|
function getCutOff(data: number[] | Float64Array, options: OptionsGetCutOff) {
|
|
@@ -1,10 +1,25 @@
|
|
|
1
1
|
import { DataXY, PeakXYWidth } from 'cheminfo-types';
|
|
2
2
|
import { optimizePeaks, OptimizePeaksOptions } from 'ml-gsd';
|
|
3
3
|
|
|
4
|
+
import { convertWidthToHz } from '../peaks/util/convertWidthToHz';
|
|
5
|
+
import { convertWidthToPPM } from '../peaks/util/convertWidthToPPM';
|
|
6
|
+
|
|
7
|
+
export interface XYPeaksOptimizationOptions extends OptimizePeaksOptions {
|
|
8
|
+
/**
|
|
9
|
+
* frequency to convert width and fwhm to ppm.
|
|
10
|
+
*/
|
|
11
|
+
frequency: number;
|
|
12
|
+
}
|
|
4
13
|
export function xyPeaksOptimization(
|
|
5
14
|
data: DataXY,
|
|
6
15
|
peaks: PeakXYWidth[],
|
|
7
|
-
options:
|
|
16
|
+
options: XYPeaksOptimizationOptions,
|
|
8
17
|
) {
|
|
9
|
-
|
|
18
|
+
const { frequency } = options;
|
|
19
|
+
const newPeaks = convertWidthToPPM(peaks, { frequency });
|
|
20
|
+
const optimizedPeaks = optimizePeaks(data, newPeaks, options);
|
|
21
|
+
return convertWidthToHz(optimizedPeaks, {
|
|
22
|
+
frequency,
|
|
23
|
+
output: optimizedPeaks,
|
|
24
|
+
});
|
|
10
25
|
}
|