integrate_peak(): Handle the situation where there is no intensity

1 files changed, 10 insertions, 3 deletions

diff --git a/src/peaks.c b/src/peaks.c
index 02ad6492..cf0b4ae0 100644
--- a/src/peaks.c
+++ b/src/peaks.c
@@ -160,9 +160,16 @@ static void integrate_peak(struct image *image, int xp, int yp,
 	}
 	}
 
-	*xc = (float)xct / total;
-	*yc = (float)yct / total;
-	*intensity = total;
+	/* The centroid is excitingly undefined if there is no intensity */
+	if ( total != 0 ) {
+		*xc = (float)xct / total;
+		*yc = (float)yct / total;
+		*intensity = total;
+	} else {
+		*xc = (float)xp;
+		*yc = (float)yp;
+		*intensity = 0;
+	}
 }