matlab读取并且处理遥感影像数据（tif等栅格）

最新推荐文章于 2025-02-21 09:34:19 发布
最新推荐文章于 2025-02-21 09:34:19 发布
阅读量3.1w 收藏 69
点赞数 6
原文链接：https://blog.csdn.net/themingyi/article/details/84842473
版权
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接： https://blog.csdn.net/themingyi/article/details/84842473 

   
   
  1. 
     
     
    
      
      
    
     
     
    
     
     
    
      
      
    
       
       Start_Latitude = 0.6;
      
      
    
     
     
  2. 
     
     
    
      
      
    
     
     
    
     
     
    
      
      
     
      
      
    
     
     
  3. 
     
     
    
      
      
    
     
     
    
     
     
    
      
      
    
       
       End_Latitude = 0.8;
      
      
    
     
     
  4. 
     
     
    
      
      
    
     
     
    
     
     
    
      
      
     
      
      
    
     
     
  5. 
     
     
    
      
      
    
     
     
    
     
     
    
      
      
    
       
       Start_Lontitude = 0.1;
      
      
    
     
     
  6. 
     
     
    
      
      
    
     
     
    
     
     
    
      
      
     
      
      
    
     
     
  7. 
     
     
    
      
      
    
     
     
    
     
     
    
      
      
    
       
       End_Lontitude = 0.3;

Matlab读取遥感影像的函数是：geotiffread（针对的是tif数据文件）

一、简介

（1）数据读取 

[A,R] = geotiffread(filename) ;

Reads a georeferenced grayscale, RGB, or multispectral image or data grid from the GeoTIFF file specified by filename into A and creates a spatial referencing object, R.（官方解释）。

描述一下就是 A 代表的数据文件（高程之类的数据），R 为空间坐标信息（包括数据的大小，经纬度范围与起始经纬度，数据单位、方向等）。

读取部分到此为止！

（2）数据处理

我这里主要介绍的数据的裁剪（比如要从一块数中选择其中的某一部分画图等）。

我们应该很清楚在做数据裁剪时会用到Matlab的冒号（：）函数，这个也是很简单。但是我们在对含有地理坐标的数据进行裁剪时，经纬度信息也要对应数据裁剪的办法进行裁剪。


 
 
  1. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Start_Latitude = 0.6;
    
    
    
   
   
  2. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  3. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     End_Latitude = 0.8;
    
    
    
   
   
  4. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  5. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Start_Lontitude = 0.1;
    
    
    
   
   
  6. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  7. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     End_Lontitude = 0.3;

以上四个部分是裁剪的核心部分，只要修改对应的这四个值就可以选择裁剪对应的数据区域。比如在一个10*10的矩阵中，上面的四个数字（0.6，0.8 ；0.1，0.3）对应的区域如图1中的绿色区域。

       

                                                                                     图1 示例图

（3）数据保存

geotiffwrite(filename,A,R) writes a georeferenced image or data grid, A, spatially referenced by R, into an output file, filename.

操作与读取类似，不再赘述！

二、例子

                                                                                        图2 原始数据

 

                                                   

                                                                             图3 原始数据地理坐标信息

 

                                               

                                                                            图4 截取之后的地理坐标信息

（1）Code：


 
 
  1. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     function Geographic_RPW()
    
    
    
   
   
  2. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     tic %Record the running time
    
    
    
   
   
  3. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     %% Read file
    
    
    
   
   
  4. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     [Data,R] = geotiffread(
     
     'D:\20181102\CSDN-博文\Matlab处理地理数据\ndvi2.tif');
    
    
    
   
   
  5. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     %% 
     
     Select interception position
    
    
    
   
   
  6. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Start_Latitude = 
     
     0.6;
    
    
    
   
   
  7. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     End_Latitude = 
     
     0.8;
    
    
    
   
   
  8. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Start_Lontitude = 
     
     0.1;
    
    
    
   
   
  9. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     End_Lontitude = 
     
     0.3;
    
    
    
   
   
  10. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     %需要补充一下此四个数的意义
    
    
    
   
   
  11. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  12. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     %% cuttiing 
     
     into certain shape
    
    
    
   
   
  13. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     Posion = [Start_Latitude,End_Latitude;Start_Lontitude,End_Lontitude];
    
    
    
   
   
  14. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     [row, col] = size(Data);
    
    
    
   
   
  15. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     first_row = floor(Posion(
     
     1,
     
     1)*row)+
     
     1;%起始维度(加
     
     1为了避免first_row=
     
     0，使用ceil函数可以不用加
     
     1)
    
    
    
   
   
  16. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     last_row = floor(Posion(
     
     1,
     
     2)*row);%结束维度
    
    
    
   
   
  17. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     first_col = floor(Posion(
     
     2,
     
     1)*col)+
     
     1;%起始经度(加
     
     1为了避免first_row=
     
     0，使用ceil函数可以不用加
     
     1)
    
    
    
   
   
  18. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     last_col = floor(Posion(
     
     2,
     
     2)*col);%结束经度
    
    
    
   
   
  19. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subImage = Data(first_row:last_row, first_col:last_col);%数据截取
    
    
    
   
   
  20. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  21. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     %% 
     
     Set GeographicCells Reference parameters
    
    
    
   
   
  22. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     %文件的地理坐标信息还有样本大小信息同样也是需要截取的和更新的，截取的基本理论类似与数据的截取办法。
    
    
    
   
   
  23. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     % mainly 
     
     set the range 
     
     of latitude 
     
     and lontitude
    
    
    
   
   
  24. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subR = R;%原始数据的地理坐标信息和数据大小等信息。
    
    
    
   
   
  25. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subR.RasterSize = size(subImage);%跟新截取后的数据大小信息
    
    
    
   
   
  26. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     del_La = subR.LatitudeLimits(
     
     2) - subR.LatitudeLimits(
     
     1);% the scale 
     
     of latitude（纬度的最大值减去最小值）
    
    
    
   
   
  27. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     del_Lon = subR.LongitudeLimits(
     
     2) - subR.LongitudeLimits(
     
     1);% the scale 
     
     of lontitude（经度的最大值减去最小值）
    
    
    
   
   
  28. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subR.LatitudeLimits = [subR.LatitudeLimits(
     
     1)+Posion(
     
     1,
     
     1)*del_La,...
    
    
    
   
   
  29. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
         subR.LatitudeLimits(
     
     2) - del_La*(
     
     1 - Posion(
     
     1,
     
     2))];%截取位置的维度范围
    
    
    
   
   
  30. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     subR.LongitudeLimits = [subR.LongitudeLimits(
     
     1)+Posion(
     
     2,
     
     1)*del_Lon,...
    
    
    
   
   
  31. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
         subR.LongitudeLimits(
     
     2) - del_Lon*(
     
     1 - Posion(
     
     2,
     
     2))];%截取位置的经度范围
    
    
    
   
   
  32. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  33. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     %% Data processing
    
    
    
   
   
  34. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     %一些最基本的数据处理，在此要根据数据的性质（大小与分布规律看一下）
    
    
    
   
   
  35. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     %在matlab调试中看一下截取后的数据的值，下面这个处理过程只是针对我
    
    
    
   
   
  36. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     %所要解决的一个问题。
    
    
    
   
   
  37. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     [row,col] = size(subImage);
    
    
    
   
   
  38. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     for i = 
     
     1:row
    
    
    
   
   
  39. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
       
     
     for j = 
     
     1:col
    
    
    
   
   
  40. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
          
     
     if subImage(i,j)<
     
     0
    
    
    
   
   
  41. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
               subImage(i,j) = 
     
     0;
    
    
    
   
   
  42. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
          
     
     elseif subImage(i,j)>
     
     0
    
    
    
   
   
  43. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
               subImage(i,j) = 
     
     1;
    
    
    
   
   
  44. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
          
     
     end
    
    
    
   
   
  45. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
       
     
     end
    
    
    
   
   
  46. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     end
    
    
    
   
   
  47. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  48. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     for i = floor(
     
     0.45*row):floor(
     
     0.8*row)
    
    
    
   
   
  49. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
       
     
     for j = floor(
     
     0.14*col):floor(
     
     0.5*col)
    
    
    
   
   
  50. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
          
     
     if subImage(i,j)<=
     
     0
    
    
    
   
   
  51. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
               subImage(i,j) = 
     
     0;
    
    
    
   
   
  52. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
          
     
     elseif subImage(i,j)>
     
     0
    
    
    
   
   
  53. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
               subImage(i,j) = 
     
     2;
    
    
    
   
   
  54. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
          
     
     end
    
    
    
   
   
  55. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
       
     
     end
    
    
    
   
   
  56. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     end
    
    
    
   
   
  57. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     %% Preview the interception posion
    
    
    
   
   
  58. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     figure(
     
     1);
    
    
    
   
   
  59. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imshow(subImage);%黑白图
    
    
    
   
   
  60. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     figure(
     
     2);
    
    
    
   
   
  61. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     imagesc(subImage);%彩色图
    
    
    
   
   
  62. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
     
    
    
    
   
   
  63. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     %% Writiing the data 
     
     to a 
     
     new GEOTIFF file
    
    
    
   
   
  64. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     filename = 
     
     'C:\Users\Administrator\Desktop\ndvi-new.tif';
    
    
    
   
   
  65. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     geotiffwrite(filename,subImage,subR);%保存文件
    
    
    
   
   
  66. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     toc%Recording the time 
     
     when the program 
     
     end running
    
    
    
   
   
  67. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     %% Writed 
     
     By themingyi
    
    
    
   
   
  68. 
   
   
    
    
    
    
   
   
    
   
   
    
    
    
    
     
     %
     
     2018/
     
     12/
     
     05

（2）处理结果

                   

 

                                                                                             图5 imshow展示结果

 

                       

                                                                                             图6 imagesc 展示结果

 

                                

                                                                                          图7 AI修改后的最终结果

 

PS：本程序以及所有使用的数据已经放到百度云供下载：

链接：https://pan.baidu.com/s/1o2vC_1QKyWmwvW8-g-E09Q 
提取码：s97g 
欢迎讨论指正，本人e-mail：themingyi123@gmail.com

 

参考：

1、https://www.mathworks.com/help/map/ref/geotiffread.html

2、https://blog.csdn.net/liyanzhong/article/details/52850581

Matlab读取带有地理信息的图像

3、http://www.ilovematlab.cn/forum.php?mod=viewthread&tid=213403

MATLAB 遥感影像的读写操作，带地理坐标信息

4、https://blog.csdn.net/qq_26906835/article/details/80841823

Matlab（一）：二维矩阵转换为geoTiff

5、http://www.mathworks.com/help/map/ref/geotiffwrite.html

 

 

高弋斌
关注
matlab遥感图像处理
04-27
包括 遥感图像读取NDVI 主成份 kt变换 IHS变换 聚类分离 傅立叶变换
matlab读取遥感图像的程序用于图像处理-readenvi.m
08-13
matlab读取遥感图像的程序用于图像处理-readenvi.m 读取遥感图片，不使用envi软件，用于图像处理。
MATLAB | 设置滑动窗口计算栅格数据的CV变异系数
最新发布
qq_48539301的博客
02-21 344
变异系数（CV）是衡量数据稳定性的重要指标，表示数据的波动程度，计算方式是标准差与均值的比值。在栅格数据分析中，较低的变异系数意味着数据变化较小、稳定性较高，而较高的变异系数则表明数据波动较大、稳定性较低。本次我们采用滑动窗口方法来计算每个时间段内栅格数据的变异系数。通过这种方法，我们可以识别出在不同时间段内变化较大或稳定的区域。同时，此代码还能提取每个滑动窗口栅格数据变异性的最小值、最大值和平均值。
Matlab 读取—处理—保存遥感影像（带有地理信息的图像）
热门推荐
themingyi的博客
12-05 4万+
Start_Latitude = 0.6; End_Latitude = 0.8; Start_Lontitude = 0.1; End_Lontitude = 0.3; Matlab读取遥感影像的函数是：geotiffread（针对的是tif数据文件） 一、简介 （1）数据读取  [A,R] = geotiffread(filename) ; Reads a georefer...
matlab读取并显示遥感img图像，包含单波段显示与多波段显示
10-03
matlab读取并显示envi标注格式img图像，支持显示单波段图像、多波段图像，用户可以根据需要选择不同波段组合显示，提供三种不同图像显示方法可供学习，另附详细注释说明
MatLab基于给定的行列号范围读取遥感影像
arcangle_123的博客
09-20 1210
问题背景 最近基于Matlab处理遥感数据时发现，遥感数据往往很大，直接基于matlab的geotiffread函数读取时，会将整景影像全部读取进内存中，消耗的时间较长，如果遥感数据过大，还会存在内存不足而读取失败的情况。因此，我打算写一个基于给的的行列号范围切片读取遥感影像的函数以解决这个问题。 function [inputRaster,R,info] = readTiff(inputRasterPath,rows,cols) % readTiff 基于给定的行列号范围读取遥感影像 % input
matlab读取并显示遥感dat图像
10-17
主要使用multibandread函数读取dat文件，显示单波段图像，多波段可以按波段显示，详细方法前面有解释，multibandread函数的参数可以根据自己的hdr文件输入
遥感图像处理.zip_matlab tif_tif图像_遥感_遥感图像 matlab_遥感图像处理
07-14
matlab代码实现遥感图像处理，图像为栅格tif类型
基于MATLAB批量读取栅格影像及处理思路
weixin_45909963的博客
05-26 5576
我个人本是要处理60景遥感影像，需逐像元处理，但是由于刚开始我对matlab掌握不深，所以不会编程，就全部用最低级的方法也就是ArcGIS来处理，two thousand years later…我用arcgis处理完毕。 后来又研究了一番，写出了完整的代码，想想都是泪/(ㄒoㄒ)/~~ 以下代码是批量读取数据后，求了一下均值，其他针对多景影像逐像元栅格计算等处理的都是同样的思路。 代码思路： （1）第一部分：批量读取文件夹下的所有文件，通用； （2）第二部分（核心）：利用for循环创建三维矩阵。一景栅格影
利用shp矢量数据裁剪遥感影像的Python实现
02-25
# 读取遥感影像 with rasterio.open('path_to_your_raster_file.tif') as src: # 获取影像的基本信息 meta = src.meta # 获取影像的band数据 data = src.read() ``` 为了裁剪遥感影像，我们需要将Shp文件中的...
matlab .nc数据格式转tif格式
03-06
本代码实现.nc数据格式转tif 、本数据使用的是地面温度和降水地面观测数据的再分析数据。数据时间分辨率是6h。
在matlab中读取遥感图像BSQ数据代码
10-17
此代码描述了如何在matlab中读取遥感中以BSQ类型存储的数据，代码描述详细，应该比较好懂
遥感数字图像处理（内含Matlab)代码
01-25
1 读取BIP 、BIL、 BSQ文件 2 均值滤波 中值滤波 3 边缘信息提取 4 DFT FFT 5 主成分变换 6 缨帽变换 7 图像分类（K—均值分类、最小距离分类、最大似然分类） 8 大气校正 反射率 地表温度的反演 9 Habib教授课程总结
matlab读取envi影像（亲测可用）
06-01
matlab读取envi标准格式影像
Matlab读取 tif文件并绘成三维图
09-14
读取 tif文件并绘成三维图，用Matlab实现，是一个很简单有效的代码哟
高光谱遥感图像MATLAB工具箱
04-11
高光谱遥感图像MATLAB工具箱 程序齐全 非常好用！！！
批量读取遥感影像数据的MATLAB技巧
5. 数据处理函数：在读取遥感影像数据后，可能需要进行数据格式转换、重采样、裁剪、归一化等预处理操作。MATLAB提供了大量的图像处理工具箱函数，例如`imresize`、`imcrop`等。 6. 脚本文件：文件SIFtiqu.m、...
MATLAB读取遥感影像
abeiceidei的博客
12-30 9640
%read envi data without head clear all fid = fopen('l7_pan','rb'); bandnum = 1; row = 400; column = 400; for i=1:bandnum    XX(1:column, 1:row) = fread(fid,[column row], 'uint8');    X(1:row,
MATLAB读取TIFF卫星遥感图像和地理信息，并转化为点云
LrFramework的博客
09-28 896
本文将介绍如何使用MATLAB读取TIFF格式的卫星遥感图像，并提取出地理信息，最后将其转化为点云数据。综上所述，本文介绍了如何使用MATLAB读取TIFF格式的卫星遥感图像，并提取出地理信息，最后将其转化为点云数据。通过上述代码，我们可以将图像的像素值从原始的无符号整型转换为双精度浮点数，并对其进行归一化处理。通过上述代码，我们可以将卫星遥感图像读入MATLAB，并获取到与该图像相关的地理信息。然后，我们将根据图像中每个像素的灰度值和对应的地理坐标，生成点云数据。最后，我们可以可视化生成的点云数据。
评论 9
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值