Matlab对GeoTiff文件的支持正在增强，不过还没有达到IDL的能力。

函数geotiffread：[A, R] = geotiffread(filename)，A返回图像矩阵，R返回地理参考。A经过一番运算以后，还是按照R地理参考进行保存，若R为Projected Coordinate System投影坐标系，以代码geotiffwrite(filename,A,R)运行返回错误结果：

Error usinggeotiffwrite>validateR (line 838)

The input, R, is a map.rasterref.MapCellsReference object indicating that you are working in a projected coordinate system. If so, then specify a projected coordinate system by setting the appropriate values for the 'CoordRefSysCode' or 'GeoKeyDirectoryTag' optional parameters.

大意：以投影坐标系作为地理参考，在保存之时geotiffwrite函数中应指明

‘CoordRefSysCode’或‘GeoKeyDirectoryTag’，两者二选其一。

示例：mm文件在投影坐标下（UTM zone 50N），经过一番运算得到m*依然试图以原地理参考保存geotiffwrite，这时须指定CoordRefSysCode或GeoKeyDirectoryTag。

查询EPSG（Projected CS Type Codes），CoordRefSysCode=32650；GeoKeyDirectoryTag参数信息来自函数geotiffinfo，代码分别以两种方法保存为m1和m2。

两种方法所得结果的地理参考完全相同，如图 1。

图 1

附测试数据及代码（tt.rar）。

% Method 2

geotiffwrite('m2.tif', image, geo, 'GeoKeyDirectoryTag', info.GeoTIFFTags.GeoKeyDirectoryTag);