高光谱数据Mat文件转换为TIFF文件

布满荆棘的人生 2022-10-10 13:57 232阅读 0赞

> 由于近期需要使用Hyperspectral Remote Sensing Scenes中的高光谱数据，所以特地记录一下数据处理的过程。

> 由于Mat文件时MATLAB中特有的文件，所以需要敲一些代码，但是整个过程是比较简单的，如下所示。

clc
    clear
    close all
    %获取Mat文件
    [fileName,pathName]=uigetfile('*.mat','Input Data-File');
    
    if isempty(fileName) || length(fileName) == 1
        fprintf("未选择mat文件！\n");
        return;
    end
    fullPath = [pathName,fileName];
    % load(fullPath);
    image = importdata(fullPath);
    shape = size(image);
    dirName = fileName(1:length(fileName)-4);
    dirPath = [pathName,dirName];
    if exist(dirPath,'dir')==0
    	mkdir(dirPath);
    end
    
    if length(shape) < 3
        fprintf("该文件不是多光谱数据\n");
        image = uint16(image);         %将灰度值转换为16bit
        str = [dirPath,'\',dirName,'.tif'];
        imwrite(image,str,'tif');
        return;
    end
    fprintf("输入数据的波段数量%d\n",shape(3));
    
    for i=1:shape(3)
        band = image(:,:,i);
        band = uint16(band);         %将灰度值转换为16bit
        str=sprintf('B%d.tif',i);
        str = [dirPath,'\',dirName,str];
        imwrite(band,str,'tif');
        fprintf("输出波段%d\n",i);
    end

> 处理之后的效果：

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70]  
![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 1]

> 特别要注意的是，灰度值要转换成16bit深度，否则有的影像呈现出来会全是黑色。此外就是这种转换出的tif文件是没有地理坐标的，所以如果要使用ENVI进行波段融合，还需要手动的进行地理配准\*^\*。

[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70]: /images/20221005/7e26f39c869f4d8387d0dd453ccaeade.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70 1]: /images/20221005/dbff5010b08e43c29a528774fd8a6188.png