从fit文件中提取lick指数的matlab程序

最新推荐文章于 2024-09-21 10:44:24 发布

爱写bug的疼老师

最新推荐文章于 2024-09-21 10:44:24 发布

阅读量1k

点赞数

分类专栏： matlab 文章标签： lick指数光谱数据 matlab

本文链接：https://blog.csdn.net/sinat_32092165/article/details/54375938

版权

本文介绍了一个MATLAB程序，用于从fit文件中提取lick指数。程序首先读取fit文件，然后通过解析头文件获取波长、步长等信息。接着，程序根据预定义的波段边界计算lick指数，涉及25个波段的处理。最终，程序将为每个fit文件生成一个25维lick数组。

摘要由CSDN通过智能技术生成

lick指数的定义及公式，可以参考以下网址：

http://d.wanfangdata.com.cn/periodical/dlxyxb201406010

本程序实现从fit文件中读取lick指数

clear;
clc;
tic;
 readDir = 'E://cas/b6';    %direction of fits








 readPath = [readDir '\*.fit'];
 readList = dir(readPath);
 [m1, n1] = size(readList);  
dim_num=3522;                                                       %total dimensions
x=0; 
 k1=1;%0-5
 k2=1;%5-10
 k3=1;%10-15
 k4=1;%15-20
 k5=1;%20-25
 k6=1;%25-30
 k7=1;%30-40
 k8=1;%40-60
 k9=1;%60-80
 k10=1;
 k11=1;
 k12=1;
 for i1 = 1:m1                                                        %templete
    filename =[readDir '\'  readList(i1, 1).name];    
 %name of file
    %splate1=fitsread(picName);