'Converting radial brightness profile to 2D intensity distribution
I have 2-column tables representing the light intensity distribution of a source as a function of its angular radius. These tables were obtained by digitizing the first plot and third plot in this figure. The x-aix is the angular radius of the two sources. The y-axis is the intensity of the sources.
Now I want to convert the 1D distributions to 2D distributions as function of longitude and latitude. Also, I want to draw the distributions as healpy maps. Can someone please tell me how this conversion can be done?
Solution 1:[1]
You should use hp.ang2pix(nside, lon, lat, lonlat=True) to go from longitude and latitude to a pixel index in a HEALPix map.
Then you can put the value of the intensity in that location in the array.
A full tutorial notebook is available here:
https://gist.github.com/zonca/680c68c3d60697eb0cb669cf1b41c324
It is linked from the docs:
https://healpy.readthedocs.io/en/latest/other_tutorials.html#pixelization
Sources
This article follows the attribution requirements of Stack Overflow and is licensed under CC BY-SA 3.0.
Source: Stack Overflow
| Solution | Source |
|---|---|
| Solution 1 | Andrea Zonca |