重點是R的設定，地理座標系使用georefcells生成R矩陣，投影座標使用maprefcells。並選取正確的座標系編碼

地理座標：

DTM=vq1;                  %二維矩陣
rasterSize=size(DTM);        %矩陣大小
R = georefcells(latlim,lonlim,rasterSize);    %latlim,lonlim分別為1*2的南北、東西經緯度界限。
CoordRefSysCode=26904;      %PCS_NAD83_UTM_zone_4N =26904    %座標系編碼，查詢地址點選開啟連結
geotiffwrite(fileName, DTM, R,'CoordRefSysCode',CoordRefSysCode);
 

投影座標：

DTM=vq1;                  %二維矩陣
rasterSize=size(DTM);        %矩陣大小
R = maprefcells(latlim,lonlim,rasterSize);    %latlim,lonlim分別為1*2的南北、東西座標界限。
CoordRefSysCode=26904;      %PCS_NAD83_UTM_zone_4N =26904    %座標系編碼，查詢地址點選開啟連結
geotiffwrite(fileName, DTM, R,'CoordRefSysCode',CoordRefSysCode);