matlab遥感图像读取与tiff输出

1、数据读取与部分矩阵变换处理

[A,R] = geotiffread('E:gree.tif') ;
k=A(1,1);
A(A == k)=nan; % 替换空值
[m,n]=size(A);
GRE=A(:);%二维矩阵转列向量

2、二维矩阵输出为tiff格式遥感图像

DTM=result_net; %二维矩阵
EPSG=32650; %坐标系编码
geotiffwrite('E:pre_net_result.tif',DTM ,R,'CoordRefSysCode', EPSG);

注：若R在前面数据处理中没有被读入，需在输出时加以设置

        地理坐标系对应编码查询网址geotiff.maptools.org/spec/geotiff6.html#6.3.2.1

3、二维矩阵在matlab中画图

imagesc(result_rf);
colorbar;
dlmwrite('result_rf.txt',result_rf);