Skip to content

Instantly share code, notes, and snippets.

View Sunmish's full-sized avatar

Stefan Duchesne Sunmish

View GitHub Profile
@Sunmish
Sunmish / basic_axes.py
Created March 28, 2024 08:00
Basic WCS axes plotting.
#! /usr/bin/env python
import os
import numpy as np
from astropy.io import fits
from astropy.wcs import WCS
from astropy.visualization import ZScaleInterval, AsymmetricPercentileInterval, simple_norm
from astropy.wcs.utils import proj_plane_pixel_scales
@Sunmish
Sunmish / multi.py
Created February 5, 2024 04:22
Use miriad to regrid/reproject a set of images to make a nice stack.
#! /usr/bin/env python
from __future__ import print_function, division
import numpy as np
from subprocess import Popen
import os
import shutil
@Sunmish
Sunmish / get_mwa_pb_lobes.py
Created May 5, 2023 02:27
Get the lobes of the MWA FEE primary beam.
#! /usr/bin/env python
# extra comment
import argparse
import numpy as np
from scipy import ndimage # For lobe finding in the primary beam images
from scipy.spatial import distance
#! /usr/bin/env python
import os
from argparse import ArgumentParser
import numpy as np
from astropy.coordinates import SkyCoord, EarthLocation, AltAz, FK5
from astropy import units as u
from astropy.io import fits
from astropy.time import Time
from subprocess import Popen
#! /usr/bin/env python
"""
Python implementation of the diffuse image filtering technique described by Rudnick (2002):
https://ui.adsabs.harvard.edu/abs/2002NewAR..46..101R/abstract
"""
from argparse import ArgumentParser
@Sunmish
Sunmish / make_stokes_cube.py
Created November 7, 2022 06:34
Create full Stokes cube from individual Stokes images
#! /usr/bin/env python
from astropy.io import fits
import numpy as np
from argparse import ArgumentParser
def get_args():
ps = ArgumentParser("Create a full Stokes cube from individual Stokes images.")
ps.add_argument("i_image", help="Stokes I image. Must be 2-dimensional (i.e. not existing Stokes or frequency axis.")
@Sunmish
Sunmish / trim_tile.py
Created May 11, 2022 06:08
Trim a tiled ASKAP image.
#! /usr/bin/env python
from argparse import ArgumentParser
import numpy as np
from astropy.io import fits
from scipy import ndimage
import logging
logging.basicConfig(format="%(levelname)s (%(module)s): %(message)s")
@Sunmish
Sunmish / merge_sourcelists.py
Last active June 14, 2022 07:27
Merge source-lists to create a catalogue.
#! /usr/bin/env python
from argparse import ArgumentParser
import os
import numpy as np
from astropy.io import fits
from astropy.table import Table, vstack
from astropy.coordinates import SkyCoord
@Sunmish
Sunmish / clip_image.py
Created June 14, 2021 09:23
Clip a FITS image using a different FITS image (e.g. of RMS noise or sensitivity).
#! /usr/bin/env python
from astropy.io import fits
import numpy as np
from argparse import ArgumentParser
def clip(image, clip_image, clip_level, outname,
fill_value=0):
"""
"""
@Sunmish
Sunmish / ms_uvcoverage.py
Created August 31, 2020 14:05
UV coverage and plot. Takes a while. Needs a lot of RAM.
#! /usr/bin/env python
from __future__ import print_function, division
import numpy as np
import sys
from astropy.constants import c; c = c.value
import matplotlib as mpl