Index: trunk/psphot/src/LocalSky.c
===================================================================
--- trunk/psphot/src/LocalSky.c	(revision 4251)
+++ trunk/psphot/src/LocalSky.c	(revision 4375)
@@ -13,7 +13,8 @@
     srcRegion = psRegionForImage (srcRegion, imdata->image, srcRegion);
 
-//  psImage *subImage = psImageSubset(imdata->image, srcRegion);
-//  psImage *subImageMask = psImageSubset(imdata->mask, srcRegion);
-//  psImage *subImageNoise = psImageSubset(imdata->noise, srcRegion);
+    //  use these when psImageSubset is updated
+    //  psImage *subImage = psImageSubset(imdata->image, srcRegion);
+    //  psImage *subImageMask = psImageSubset(imdata->mask, srcRegion);
+    //  psImage *subImageNoise = psImageSubset(imdata->noise, srcRegion);
 
     psImage *subImage = psImageSubset(imdata->image, srcRegion->x0, srcRegion->y0, srcRegion->x1, srcRegion->y1);
@@ -47,7 +48,9 @@
     psImageMaskRegion (mask, srcRegion, AND, 0x7f);
 
-    source->moments = pmMomentsAlloc();
     psF64 tmpF64;
     p_psGetStatValue(myStats, &tmpF64);
+    if (isnan(tmpF64)) return (false);
+
+    source->moments = pmMomentsAlloc();
     source->moments->Sky = (psF32) tmpF64;