上回说到要是能根据关键词确定数组指标就更好了（http://blog.sciencenet.cn/home.php?mod=space& uid=117333&do=blog&id=641804）。经过搜索发现，Position函数可以实现这一点。
        对于
header=Import["xxx.fits","Metadata"]
得到的header列表，要得到NAXIS3所在的行数，尝试
c=Position[header,"NAXIS3"]
得到
{{6,1}}
再尝试
d=c[[1]][[1]]
就得到行数6了。这样就可以结合之前用Mathematica读FITS文件数据的操作进行些简单的数据处理了。比如说得到(x,y)=(10,9)处的频谱
a=Import["xxx.fits","RawData"];
header=Import["xxx.fits","Metadata"];
c=Position[header,"NAXIS3"];
d=c[[1]][[1]];
nz=header[[d]][[2]][[1]];
spec=a[[Table[n,{n,1,nz}],9,10]];

转载本文请联系原作者获取授权，同时请注明本文来自钱磊科学网博客。
链接地址：https://blog.sciencenet.cn/blog-117333-641854.html

上一篇：Mathematica处理分子云数据的可能性（一）读FITS文件
下一篇：Mathematica处理分子云数据的可能性（三）二维图像显示