matlab：二维矩阵转换为geoTiff

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

本文链接： https://blog.csdn.net/qq_26906835/article/details/80841823

重点是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);