IDL生成的遥感图像利用matlab读取其中的数据【讨论】

wudi1201

初学matlab，求围观，求交流
前几日，得到一些由IDL软件生成的tif格式的遥感图像（在windows下没有预览）如：http://bbs.06climate.com/forum.php?mod=attachment&aid=MTI0MDB8YjVlYTA5NjM1YTIzOTg3YjViOWIzN2Q1MTIyODlmNzV8MTcxNDkyOTA2MQ%3D%3D&request=yes&_f=.tif；
参考一片文献http://bbs.06climate.com/forum.php?mod=attachment&aid=MTI0MDF8NzliNDRiOWM2MWEwN2U5MDM2ODQ4MGQxMWFmYTIyNTB8MTcxNDkyOTA2MQ%3D%3D&request=yes&_f=.pdf上得作法，自己写了一下代码；

最初开始想的简单以为这样做：：
%第一步：
[data,map]=imread('SoilMoisture2012001.tif');
%第二步：
[x,y]=size(data);
%范围
%image(data)
%colormap(data)
fi=fopen('E:\data\modis\SoilMoisture\SoilMoisture\2012\SoilMoisture2012001.dat','w');
for i=1:x
    for j=1:y
        fprintf(fi,'%6f',data(i,j));
    end
end

惭愧的是，得到的dat文件fortran不支持，不支持就真心没办法了，想看看里面的内容也就不能了。也曾尝试过转成ascii格式的，但是图像文件不能直接保存为ascii格式


后来细细研读文章发现，原文中的处理遥感图像中提取数据主要是依赖于色标中不同像素对应的数值。而我得到的文件没有色标，由于图像不能直接看，所以只能在matlab中imshow（'SoilMoisture2012001.tif'），得到图像如下：

imshow

没有色标，利用matlab查看图片的info，info = imfinfo('SoilMoisture2012001.tif')，得到信息如下：
>> info = imfinfo('SoilMoisture2012001.tif')
info =
                     Filename: 'D:\MATLAB\data\Modis\SoilMoisture2012001.tif'
                  FileModDate: '02-十一月-2012 14:36:18'
                     FileSize: 501132
                       Format: 'tif'
                FormatVersion: []
                        Width: 411
                       Height: 304
                     BitDepth: 32
                    ColorType: 'grayscale'
              FormatSignature: [73 73 42 0]
                    ByteOrder: 'little-endian'
               NewSubFileType: 0
                BitsPerSample: 32
                  Compression: 'Uncompressed'
    PhotometricInterpretation: 'BlackIsZero'
                 StripOffsets: [76x1 double]
              SamplesPerPixel: 1
                 RowsPerStrip: 4
              StripByteCounts: [76x1 double]
                  XResolution: 100
                  YResolution: 100
               ResolutionUnit: 'Inch'
                     Colormap: []
          PlanarConfiguration: 'Chunky'
                    TileWidth: []
                   TileLength: []
                  TileOffsets: []
               TileByteCounts: []
                  Orientation: 1
                    FillOrder: 1
             GrayResponseUnit: 0.0100
               MaxSampleValue: 4.2950e+009
               MinSampleValue: 0
                 Thresholding: 1
                       Offset: 499784
                 DocumentName: [1x54 char]
             ImageDescription: 'IDL TIFF file'
                     Software: 'IDL 6.1, Research Systems, Inc.'
                     DateTime: '2012:11:02 14:36:17'
                 SampleFormat: 'IEEE floating point'
           ModelPixelScaleTag: [3x1 double]
             ModelTiepointTag: [6x1 double]
           GeoKeyDirectoryTag: [76x1 double]
           GeoDoubleParamsTag: [8x1 double]
            GeoAsciiParamsTag: 'WGS-84|'

上面的信息也看不太懂，其中最主要的几个是：ColorType: 'grayscale'中像素和色标的对应关系;GeoAsciiParamsTag: 'WGS-84|'这个经纬度怎么得到，（原来图片中没有经纬度信息，据说是要利用这个的？？？）


最后：遥感图像资料处理的都有哪些方法？有了解的可以多多交流，私信可以QQ

wudi1201

哎，帖子中添加的附件地址分别是：
1，http://bbs.06climate.com/forum.p ... est=yes&_f=.tif
2，http://bbs.06climate.com/forum.p ... est=yes&_f=.pdf

尉鹏

matlab 有个现成的函数可以读idl结果的吧

wudi1201

尉鹏 发表于 2013-3-8 11:51

                               
登录/注册后可看大图

matlab 有个现成的函数可以读idl结果的吧

你好，有现成的函数？具体是哪个？能给点参考嘛？谢谢了！

尉鹏

X = multibandread(filename, size, precision, offset, interleave, byteorder)
参数到hdr里面找就行了，具体看matlab帮助。
我用过，可以读。

wudi1201

尉鹏 发表于 2013-3-8 13:52

                               
登录/注册后可看大图

X = multibandread(filename, size, precision, offset, interleave, byteorder)
参数到hdr里面找就行了， ...

好的，多谢了。

尉鹏

lpf20092300226

好想学习啊

kongfeng0824

尉鹏 发表于 2013-3-8 11:51

                               
登录/注册后可看大图

matlab 有个现成的函数可以读idl结果的吧

你指的是什么函数呢？

尉鹏

看帖子，不是都说了么