函数geotiffread：[A, R] = geotiffread(filename)，A返回图像矩阵，R返回地理参考；

    函数geotiffwrite：geotiffwrite(filename,A,R)，filename：输入文件的名称和位置，文件扩展名必须是 “.tif” 或者 “.TIF”；A：要导出的影像数据；R：地理参考信息；

    重点是R的设置，地理坐标系使用georefcells生成R矩阵，投影坐标使用maprefcells。并选取正确的坐标系编码！

示例（将重采样后的数据DInSAR形变数据写成goetiff格式的数据）：

    将形变数据重采样到1500*2200的格网之中；

    在matlab中将其画出，如图所示。

%运行以下代码%
DTM=A_dinsar;
rasterSize=size(DTM);
latlim=[35.2812232971191,36.1139984130859];
lonlim=[-118.138000488281,-116.916336059570];
R=georefcells(latlim,lonlim,rasterSize,'ColumnsStartFrom','north');
geotiffwrite('a_dinsar.tif',A_dinsar,R);

    就得到了a_dinsar.tif文件，其中R为下图所示。

    用ARCGIS等软件打开生成的tif文件，如图所示。