我有两列表格,表示光源的光强度分布与其角半径的函数关系。这些表格是通过对第一幅图和第三幅图进行数字化而获得的在此图中。 x-aix 是两个源的角半径。 y 轴是源的强度。
现在我想将一维分布转换为二维分布作为经度和纬度的函数。另外,我想将分布绘制为healpy 地图。有人可以告诉我如何完成这种转换吗?
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?
发布评论
评论(1)
您应该使用
hp.ang2pix(nside,lon,lat,lonlat = true)
从healpix地图中的经度和纬度转到像素索引。然后,您可以将强度的值放在数组中的该位置。
完整的教程笔记本可在此处提供:
它是从文档链接的:
https:> https:// en/lestest/other_tutorials.html#像素化
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