Teile eines Arrays integrieren

Also willst du im Grunde die Gesamtintensität jedes einzelnen Objekts aufsummieren, richtig?
Untergrund würde ich auch erstmal mit threshold auf null setzen. Wenn der Untergrund variiert, kann man das auch ganz gut über den ersten Peak beim Intensitäts-Histogramm lösen.
Zusammenhängende Objekte...das wird interessant. Man könnte eine Art "Fressverfahren" machen:

1) Das komplette untergrundbereinigte 2D-Bild-Array nach einem Wert > 0 durchsuchen.
2) Sobald ein Punkt > 0 gefunden wird, wird dieser in ein neues Array als Objektpunktearray eingetragen und der Punkt auf 0 gesetzt.
3) Nachbarpunkte (x+/-1, y+/-1, wenn man ganz auf Nummer sicher gehen will, sollte man auch diagonal springen) nach weiteren Werten > 0 absuchen und dem gleichen Objektpunktearray zuschreiben, Wert wieder auf 0 setzen. Schritt 3 bei den jeweiligen Nachbarpunkten > 0 solange wiederholen, bis alle Nachbarn null sind.
4) Zurück zu Schritt 1) (mit angepasstem Startindex) , das nächste Objektpunktearray vollschreiben und solange wiederholen, bis das 2D-Bild-Array durch ist.

Auswertungsmöglichkeiten:
a) Anzahl des Objektpunktearrays bewerten. Hat das Objekt nur 1 oder 2 Punkte wird's als Cosmic klassifiziert und rausgeschmissen.
b) Mittelwert des Objektpunktearrays bilden
c) Koordinaten des Objektpunkte festhalten (Dann sollte es allerdings eher ein array von Objektpunktclustern werden) und so über min/max die Ausmaße des Objekts bestimmen. Oder auch jeweils das gefundene Einzelobjekt nochmal "plotten".

Hattest du dir sowas in der Art vorgestellt?