Localizing hot spots in Poisson radiation data matrices: nonnegative tensor factorization and phase congruency

Detecting and delineating hot spots in data from radiation sensors is required in applications ranging from monitoring large geospatial areas to imaging small objects in close proximity. This paper describes a computational method for localizing potential hot spots in matrices of independent Poisson data where, in numerical terms, a hot spot is a cluster of locally higher sample mean values (higher Poisson intensity) embedded in lower sample mean values (lower background intensity). Two numerical algorithms are computed sequentially for a 3D array of 2D matrices of gross Poisson counts: (1) nonnegative tensor factorization of the 3D array to maximize a Poisson likelihood and (2) phase congruency in pertinent matrices. The indicators of potential hot spots are closed contours in phase congruency in these matrices. The method is illustrated for simulated Poisson radiation datasets, including visualization of the phase congruency contours. The method may be useful in other applications in which there are matrices of nonnegative counts, provided that a Poisson distribution fits the dataset.