Healpy显示网格和中点

Question

使用python healpy，我想要显示(和导出)用于天空细分(在不同级别)的healpix网格。我该怎么做呢？理想情况下，我需要中点，网格的顶点和连接它们的线。

任何帮助都是非常感谢的。

Answer 1

我已经做过了，还有更多的方法：

1)你可以使用IDL库，它有一个函数pix2vec_nest (和pix2vec_ring)，可以有选择地返回向量(n，3，4)中像素的4个顶点。

2)可以使用有HealpixBase.boundaries方法的java版本

3)可以使用具有相同行为的T_Healpix_Base::boundaries的HEALPix_cxx库。

有了边界，您可以显示栅格。请注意，最后一种方法编码得非常好，并且很容易理解和在另一个程序中使用！它只需要像素索引和输出向量：

 template<typename I> void T_Healpix_Base<I>::boundaries(I pix, tsize step,
  vector<vec3> &out) const
  {
  out.resize(4*step);
  int ix, iy, face;
  pix2xyf(pix, ix, iy, face);
  double dc = 0.5 / nside_;
  double xc = (ix + 0.5)/nside_, yc = (iy + 0.5)/nside_;
  double d = 1.0/(step*nside_);
  for (tsize i=0; i<step; ++i)
    {
    double z, phi, sth;
    bool have_sth;
    xyf2loc(xc+dc-i*d, yc+dc, face, z, phi, sth, have_sth);
    out[i] = locToVec3(z, phi, sth, have_sth);
    xyf2loc(xc-dc, yc+dc-i*d, face, z, phi, sth, have_sth);
    out[i+step] = locToVec3(z, phi, sth, have_sth);
    xyf2loc(xc-dc+i*d, yc-dc, face, z, phi, sth, have_sth);
    out[i+2*step] = locToVec3(z, phi, sth, have_sth);
    xyf2loc(xc+dc, yc-dc+i*d, face, z, phi, sth, have_sth);
    out[i+3*step] = locToVec3(z, phi, sth, have_sth);
    }
  }