search

RadioAstronomySoftwareGroup / pyradiosky

Python objects and interfaces for representing diffuse, extended and compact astrophysical radio sources
https://pyradiosky.readthedocs.io/en/latest/
BSD 2-Clause "Simplified" License
14 stars 5 forks source link

write healpix FITS files #78

Open bhazelton opened 4 years ago

bhazelton commented 4 years ago

(from notes in the code that I removed)

mkolopanis commented 4 years ago

Do we have copies of the specifications we want? Fits and healpix go together fairly well, the CMB community has been doing this for a while and should not require a ton of work.

bhazelton commented 4 years ago

We need to talk about this. Standard HEALPix FITS files do not support a frequency axis. We need to decide whether to make a new non-standard file type or make different HDUs or files for different frequencies.

mkolopanis commented 4 years ago

Yeah we could end up with a large number of HDUs if all the frequencies are put into one file. I've seen things similar to this before, but usually integrated bandpowers which result in 5ish HDUs (e.g. WMAP cmb maps).

jpober commented 3 years ago

This may be particularly useful for Ruby's catalog. We'll probably need a read method also.

jpober commented 3 years ago

@rlbyrne Have your plans for releasing your catalog/map become any clearer? Do you still plan on a skyh5 and a more standard HEALPix release? For the latter, how do you plan to create that file?

rlbyrne commented 3 years ago

@jpober Yes, I'm planning to release the data in both formats. For the Healpix fits format, I have code separate from pyradiosky. That could get integrated into pyradiosky at some point, but Bryna's concerns about the frequency axis still stand.

jpober commented 3 years ago

Ok, I think we will want to consider that code for integration into pyradiosky (at some point). I can still see at being nice functionality for single frequency maps.