From 65a54acbe7c2fdf5666f80d6145097aa35296379 Mon Sep 17 00:00:00 2001 From: Aaron LI Date: Fri, 25 Aug 2017 11:36:46 +0800 Subject: ps2d.py: Use [arcsec] for pixel size --- astro/ps2d.py | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/astro/ps2d.py b/astro/ps2d.py index d3d988b..3d20720 100755 --- a/astro/ps2d.py +++ b/astro/ps2d.py @@ -60,8 +60,8 @@ class PS2D: window_name=None, window_width="extended"): logger.info("Initializing PS2D instance ...") self.cube = cube - self.pixelsize = pixelsize # [arcmin] - logger.info("Image pixel size: %.2f [arcmin]" % pixelsize) + self.pixelsize = pixelsize # [arcsec] + logger.info("Image pixel size: %.2f [arcsec]" % pixelsize) self.frequencies = np.array(frequencies) # [MHz] self.nfreq = len(self.frequencies) # Central frequency and redshift @@ -167,7 +167,7 @@ class PS2D: @property def k_xy(self): __, ny, nx = self.cube.shape - dxy = self.DMz * np.deg2rad(self.pixelsize / 60.0) # [Mpc] + dxy = self.DMz * np.deg2rad(self.pixelsize / 3600.0) # [Mpc] kx = 2*np.pi * fftpack.fftshift(fftpack.fftfreq(nx, dxy)) ky = 2*np.pi * fftpack.fftshift(fftpack.fftfreq(ny, dxy)) return (kx, ky) # [Mpc^-1] @@ -255,7 +255,7 @@ def main(): help="overwrite existing file") parser.add_argument("-p", "--pixelsize", dest="pixelsize", type=float, required=True, - help="image cube pixel size; unit: [arcmin]") + help="image cube pixel size; unit: [arcsec]") parser.add_argument("--window", dest="window", choices=["nuttall"], help="apply window along frequency axis " + -- cgit v1.2.2