From 8b2ef6450cd4013e3e957bac912831b476aff6ca Mon Sep 17 00:00:00 2001 From: Aaron LI Date: Mon, 10 Oct 2016 20:55:15 +0800 Subject: utils: zea2healpix(): Implement inpainting missing pixels The missing pixels in the reprojected HEALPix map are filled with the averages of their 8 neighboring pixels (excluding the NaN's if any). Also add the "--inpaint" argument to the executable script. --- bin/zea2healpix | 3 +++ 1 file changed, 3 insertions(+) (limited to 'bin') diff --git a/bin/zea2healpix b/bin/zea2healpix index 1eb8730..33cb694 100755 --- a/bin/zea2healpix +++ b/bin/zea2healpix @@ -37,6 +37,8 @@ def main(): parser.add_argument("-O", "--interp-order", dest="interp_order", type=int, default=1, help="interpolation order (integer, 0-5; default: 1)") + parser.add_argument("-I", "--inpaint", action="store_true", + help="inpaint the missing pixels if present") parser.add_argument("-C", "--clobber", action="store_true", help="overwrite the existing output file") parser.add_argument("-l", "--log", dest="loglevel", default=None, @@ -77,6 +79,7 @@ def main(): hp_data, hp_header, __ = zea2healpix(args.infile1, args.infile2, nside=args.nside, order=args.interp_order, + inpaint=args.inpaint, append_history=history, append_comment=comments) write_fits_healpix(args.outfile, hpmap=hp_data, header=hp_header, -- cgit v1.2.2