FIT格式图像的读写

在matlab中关于FITS图像只有两个函数，一个是读函数fitread（），一个是读取信息函数fitinfo(),也就是说matlab无法对fit文件进行建立，和写入的工作。
matlab中在显示fit文件的时候，如果文件像素过多，将无法显示，只能显示一部分。我想这对其它格式的图像也有这种情况。
关于FIT图像的读写和操作，c和c++都有成熟的库可以调用。具体参见：http://fits.gsfc.nasa.gov/fits_libraries.html 这里有支持FITS图像操作的各种开发库的介绍。