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?
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