fix(optimizer): fix output and include original indexs

Author: jobo322

src/post/joinBroadPeaks.js

@@ -23,7 +23,7 @@ export function joinBroadPeaks(peakList, options = {}) {
   broadLines.push({ x: Number.MAX_VALUE });
 
   let candidates = [[broadLines[0].x, broadLines[0].y]];
-  let indexes = [0];
+  let indexes = [broadLines[0].index];
 
   for (let i = 1; i < broadLines.length; i++) {
     // console.log(broadLines[i-1].x+" "+broadLines[i].x);
@@ -33,7 +33,7 @@ export function joinBroadPeaks(peakList, options = {}) {
         max = broadLines[i].y;
         maxI = i;
       }
-      indexes.push(i);
+      indexes.push(broadLines[i].index);
       count++;
     } else {
       if (count > 2) {
@@ -49,6 +49,9 @@ export function joinBroadPeaks(peakList, options = {}) {
           x: parameters[0],
           y: parameters[1],
           width: parameters[2],
+          index: Math.floor(
+            indexes.reduce((a, b) => a + b, 0) / indexes.length,
+          ),
           soft: false,
         });
       } else {

src/post/optimizePeaks.js

@@ -48,11 +48,14 @@ export function optimizePeaks(peakList, x, y, options = {}) {
             );
           }
         }
+
         for (let j = 0; j < optPeaks.length; j++) {
+          let { parameters } = optPeaks[j];
           result.push({
-            x: optPeaks[j][0][0],
-            y: optPeaks[j][1][0],
-            width: optPeaks[j][2][0] * factor,
+            x: parameters[0],
+            y: parameters[1],
+            width: parameters[2] * factor,
+            index: peaks[j].index,
           });
         }
       }
@@ -84,11 +87,13 @@ export function optimizePeaks(peakList, x, y, options = {}) {
             );
           }
         }
+
         let { parameters } = fitResult;
         result.push({
           x: parameters[0],
           y: parameters[1],
           width: parameters[2] * factor,
+          index: peaks.index,
         }); // From https://en.wikipedia.org/wiki/Gaussian_function#Properties}
       }
     }