astrophysics
22 posts
-
healpy 1.20.0b1 beta — harmonic_ud_grade for artifact-safe map downgrading
pip install --pre healpy==1.20.0b1 harmonicudgrade changes map NSIDE through spherical-harmonic transforms with pixel-window and beam transfer corrections following Eq. 1 of Planck 2015 X (arXiv:1502.01588).
-
healpy harmonic_ud_grade vs skytools change_resolution: A Comparison
A comparison of healpy.harmonic_ud_grade and skytools.change_resolution for HEALPix map resolution changes, highlighting API differences, beam/pixel-window support, direct alm input, and agreement levels.
-
harmonic_ud_grade vs ud_grade: A Focused Comparison
harmonic_ud_grade vs ud_grade A Focused Comparison.
-
Why 3 pixels per beam? The mathematics of HEALPix resolution matching
Exploring why the Planck convention uses FWHM = 3 x pixel size for harmonic-space map downgrading, with equations from the literature and simulated demonstrations.
-
Healpy blm_gauss breaking change analysis
Analyzing the breaking change in healpy's blm_gauss function (1.18.1 to 1.19.0) and its impact on Gaussian beam coefficients.
-
Astrophysics Papers Daily Summaries Notebook with the Jetstream LLM Inference Service
Fetch today's astro-ph papers from arXiv, summarize abstracts with llm, and output a Markdown summary. This is an example of using the Jetstream Inference Service, notice that you need first to configure the llm package to access the Jetstream Inference Service via the API
-
Interactive 3D plot of a Planck map with Matplotlib
This is an update of my older tutorial Unfortunately I was not able to make Mayavi properly plot the data, rerunning the same code would give a map that looked like uniform noise.
-
CFITSIO file writing with checkpointing
I wrote a prototype implementation of writing FITS files in C checkpointing using the FITS flush function so that if the code segfaults, the file is not corrupted, see:
-
Simple WCS with astropy modeling and gwcs
Some notes on how to convert a FITS WCS to a WCS from gwcs (Generalized WCS) to be used within the JWST pipeline python import astropy from astropy.io import fits import numpy as np from astropy import units as u
-
Advanced pandas with astrophysics example notebook
Taught a lesson today on advanced python and pandas based on an example application in Astrophysics with simulations of data from the Planck Satellite), features also a Binder button to run it yourself. Jupyter Notebook available at: under CC-BY
-
Accelerate groupby operation on pixels with Numba
Download the original IPython notebook Astrophysics background
-
clviewer, interactive plot of CMB spectra
Today it was HackDay at .Astronomy, so I felt compelled to hack something around myself, creating something I have been thinking for a while after my previous work on Interactive CMB power spectra in the browser
-
Planck CMB map at high resolution
Prompted by a colleague, I created a high-resolution version of the Cosmic Microwave Background map in MollWeide projection released by the Planck collaboration, available on the Planck Data Release Website in FITS format.
-
Interactive figures in the browser - CMB power spectra
For a long time I've been curious about trying out d3.js , the javascript plotting library which is becoming the standard for interactive plotting in the browser.
-
Noise in spectra and map domain
Spectra NET or $\sigma$ is the standard deviation of the noise, measured in mK/sqrt(Hz), typical values for microwave amplifiers are 0.2-5. This is the natural unit of the amplitude spectra (ASD), therefore the high frequency tail of the ASD should get to the expected value of the NET.