1

我想使用 ccfits 创建一组未平方的图像。我可以在 primaryHDU 中制作一个，如下所示：

long axes[2] = { jmax, imax };   
std::auto_ptr<CCfits::FITS> pFits(0);
pFits.reset ( new CCfits::FITS ( "fitfile.fits", FLOAT_IMG, 2, axes ) );

std::valarray<double> h2a0array ( jmax * imax );
for ( int i = 0 ; i < imax ; i++ 
  for ( int j = 0 ; j < jmax ; j++ )
    h2a0array [ j + jmax * i ] = i + j;  

pFits->pHDU().write  ( fpixel, imax * jmax, h2a0array );

但我不知道如何将其他非平方图像添加到我的 FITS 文件中。我想我必须使用 CCFITS::addImage 函数，但只能使用它获取平方图像：

long fpixel ( 1 );
std::vector<long> extAx ( 2, dim );
CCfits::ExtHDU* imageExt2 = pFits->addImage ( "h2a0array", FLOAT_IMG, extAx );
imageExt2->write ( fpixel, imax * jmax, h2a0array );

extAx 向量仅包含两个值，第一个是要添加到 FITS 文件的图像的维度（1D、2D、3D），第二个是其大小。我不知道将图像添加到 FITS 文件的任何其他方法。如果有人这样做，我们非常欢迎您的帮助！

谢谢，阿诺。

4

1 回答 1

0

addImage 的最后一个参数的向量可以有任意维度和维度中不同的轴长度。不要求轴长度相同（您似乎称之为“方形”）：

vector<long> extAx ;
extAx.push_back(imax) ; 
extAx.push_back(jmax) ;
extAx.push_back(kmax) ;

pFits->addImage("h2a0array", FLOAT_IMG, extAx );
于 2014-06-24T13:16:07.560 回答