FITS基本文件格式

您所在的位置:网站首页 fz格式 FITS基本文件格式

FITS基本文件格式

2024-07-14 17:31| 来源: 网络整理| 查看: 265

FITS (Flexible Image Transport System) 是国际天文学会(IAU)1982 年确定的世界各天文台之间用于数据传输、交换的统一标准格式。它描述了数据的定义和数据编码的一般方法。它是与机器无关的,用磁带作为标准传输介质的独立方法。它提供了图象的单值转换,精度包括符号在内可以达到 32 位。对一维、二维、三维、甚至多维的数据类型都提供了合适的转换。它不仅适用于天文数据, 对其他学科的数据也是可用的。 下图是FITS的基本格式: FITS文件格式

下面的实例程序是要求用户输入待分析的fits文件,然后将fits文件分成2880byte的块,然后写入filenam.out的文件,这样就可以全局分析这个fits的文件了。

#!/usr/bin/python #author yxy import os from os.path import getsize filename= raw_input('Please input you fits file name: ') filesize = int(getsize(filename)) blockNum = filesize/2880 try: fin = open(filename, 'r') fout = open(filename + '.out', 'w') except IOError, e: print 'could not open file:', e for i in range(blockNum): line = fin.read(2880) fout.write('\n' * 5) fout.write(str(i) * 10) fout.write('\n' * 5) j=0 for i in range(36): fout.write(line[j:j+80]) fout.write('\n') j +=80 fout.close()

基本fits的文件格式大体上如上图所示,但是还有一个关键点是文件头中36行标题的格式。文件头中必须存在的五个关键词是: SIMPLE,BITPIX,NAXIS,NAXISn 和 END。它们的含义是: SIMPLE 逻辑变量,指明文件是否符合基本 FITS 标准。 BITPIX 整数变量,指明用于表示每一个象元值的位数。 NAXIS 整数变量,指明图象里坐标轴数。 NAXIS1 整数变量,指明在数组内沿最快变化轴的象元数。 NAXIS2 整数变量,指明在数组内沿次快变化轴的象元数。

下面是一个例子:

SIMPLE = T BITPIX = 16 NAXIS = 2 NAXIS1 = 320 NAXIS2 = 512 END

SIMPLE=T 其中逻辑值 T(True)表明该磁带文件是简单的 FITS 数据文件。BITPIX=16 表明每一象元值为16 位带符号的补码形式,每一象元值占用两个字节, 每一记录可记 1440 个象元值(23040÷16= 1440)。NAXIS=2 表明图象为两维阵列NAXIS1=320 表明第一轴 (列变化方向)的长度为 320 个象元。NAXIS2=512 表明第二轴(行变化方向)的长度为 512行。END表明标题结束。记录后面的所有零字节都用 ASCII 码空格填满。 这五个关键次的顺序必须如上所示。不能颠倒。

除了这五个必须的关键次还有一些其他的关键词: ——BSCALE(浮点数)和 BZERO(浮点数):用于将磁带上的象元值转换为真实值的参数:真实值= ( 磁带值×BSCALE)+BZERO ——BUNIT(字符):单位 ——BLANK (整数), 如果某象元未定义值, 则磁带上该象元对应位置上所写入的就是由 BLANK 定义的数值。 ——OBJECT (字符): 图象名 ——DATE (字符): 写 FITS 文件的日期(‘日/月/年’)。 ——DATE_OBS (字符): 数据的观测日期(‘日/月/年’)。 ——ORIGIN (字符): 写磁带的研究单位。 ——INSTRUME (字符): 数据的探测仪器。 ——TELESCOP (字符): 所使用的望远镜。 ——OBSERVER (字符):观测者姓名。 ——CRVALn (浮点数), CRPIXn(浮点数), CDELTn(浮点数), CTYPEn(字符)和 CROTAn(浮点数):这些关键词是用于定义第 n个轴的象元坐标与真实的物理坐标之间的关系。 ——DATAMAX(浮点数)和 DATAMIN(浮点数): 给出文件中数据的最大值和最小值。 ——EPOCH(浮点数): 坐标系统的历元(以年为单位)。 ——空白(无)、COMMENT(无)和 HISTORY(无): 用于说明其它一些信息表明第 9-80 列是注释用的 ASCII 文本。

注意, 图象阵中象元的二进制数据是高位在前, 低位在后。

参考文献: http://www.cnblogs.com/tyjsjl/p/3952823.html



【本文地址】


今日新闻


推荐新闻


CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3