From c3e9a6b4d5acdc4236d24ccbcd3a9a690cbca3b8 Mon Sep 17 00:00:00 2001 From: Brian Benjamin Maranville Date: Tue, 6 Aug 2019 17:32:07 -0400 Subject: [PATCH] Bug in logic for 2d Gaussian There is a logic error in the code for quickly fitting a 2d gaussian - the wrong coordinate offset (x and y are swapped) is used in the center-finding along a row and column. --- ipython/FittingData.ipynb | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/ipython/FittingData.ipynb b/ipython/FittingData.ipynb index 65f52d7..a147d95 100644 --- a/ipython/FittingData.ipynb +++ b/ipython/FittingData.ipynb @@ -368,9 +368,9 @@ " x = (X*data).sum()/total\n", " y = (Y*data).sum()/total\n", " col = data[:, int(y)]\n", - " width_x = np.sqrt(np.abs((np.arange(col.size)-y)**2*col).sum()/col.sum())\n", + " width_x = np.sqrt(np.abs((np.arange(col.size)-x)**2*col).sum()/col.sum())\n", " row = data[int(x), :]\n", - " width_y = np.sqrt(np.abs((np.arange(row.size)-x)**2*row).sum()/row.sum())\n", + " width_y = np.sqrt(np.abs((np.arange(row.size)-y)**2*row).sum()/row.sum())\n", " height = data.max()\n", " return height, x, y, width_x, width_y\n", "\n",