diff options
author | Aaron LI <aly@aaronly.me> | 2018-05-26 17:09:25 +0800 |
---|---|---|
committer | Aaron LI <aly@aaronly.me> | 2018-05-26 17:09:25 +0800 |
commit | 2455aff9436831768535c0ecb672aa377867185e (patch) | |
tree | 2b4f2463f8206dd6db136f264a04c41e5a96e40e | |
parent | 1168808b44bc43a94183566485bedafa37d6c1e0 (diff) | |
download | atoolbox-2455aff9436831768535c0ecb672aa377867185e.tar.bz2 |
astro/fitsimage.py: ignore nan
-rwxr-xr-x | astro/fits/fitsimage.py | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/astro/fits/fitsimage.py b/astro/fits/fitsimage.py index db7030e..5f867d3 100755 --- a/astro/fits/fitsimage.py +++ b/astro/fits/fitsimage.py @@ -180,12 +180,12 @@ def show_info(filename, abs_=None, center=None): rc, cc = rows//2, cols//2 cs1, cs2 = center//2, (center+1)//2 data = data[(rc-cs1):(rc+cs2), (cc-cs1):(cc+cs2)] - mean = np.mean(data) - median = np.median(data) - std = np.std(data) - iqr = np.diff(np.percentile(data, q=(25, 75))) - mad = np.median(np.abs(data - median)) - rms = np.sqrt(np.mean(data**2)) + mean = np.nanmean(data) + median = np.nanmedian(data) + std = np.nanstd(data) + iqr = np.diff(np.nanpercentile(data, q=(25, 75))) + mad = np.nanmedian(np.abs(data - median)) + rms = np.sqrt(np.nanmean(data**2)) print("mean: %13.6e" % mean) print("median: %13.6e" % median) print("std: %13.6e (standard deviation)" % std) |