MATLAB中自定义栅格数据地理坐标R,利用geotifwrite写入tif

场景描述:

有时候将nc格式的数据转成tiff,或者是将一个矩阵输出成带有地理坐标信息tiff数据时,常常涉及到空间参考的定义和geotiffwrite()函数。

问题描述:

以全球数据为例,今天发现在matlab中对矩阵进行显示后,该数据是正常的,但是添加了空间参考输出后,数据的南北维度颠倒了。如图1是正常显示的结果,图2是输出后的南北纬被颠倒的结果。图下方是一段示例代码。

 正常显示的结果应该如下所示:

上述倒置的图像的代码是:

A=zeros(180,360);  //创建一个180行,360列的矩阵
A=function;  //对A进行复制计算,这里就用function来指代
Ref=georasterref('RasterSize',[180,360],'Latlim',[-90,90],'Lonlim',[-180,180]);  //定义空间参考信息
geotiffwrite(outfilepath, A, Ref); //将A输出到tiff文件中

解决方法:

在完成了矩阵A的赋值后,应该在后面添加一个rot90()函数将原来的数据倒置一下,所以纠正后的代码应该如下。

A=zeros(180,360);  //创建一个180行,360列的矩阵
A=function;  //对A进行复制计算,这里就用function来指代
A=rot90(A');  //将矩阵A倒置一下,即南北维度调换
Ref=georasterref('RasterSize',[180,360],'Latlim',[-90,90],'Lonlim',[-180,180]);  //定义空间参考信息
geotiffwrite(outfilepath, A, Ref); //将A输出到tiff文件中

最近更新

  1. docker php8.1+nginx base 镜像 dockerfile 配置

    2024-05-04 07:36:01       94 阅读
  2. Could not load dynamic library ‘cudart64_100.dll‘

    2024-05-04 07:36:01       101 阅读
  3. 在Django里面运行非项目文件

    2024-05-04 07:36:01       82 阅读
  4. Python语言-面向对象

    2024-05-04 07:36:01       91 阅读

热门阅读

  1. pytorch笔记:ReplicationPad1d

    2024-05-04 07:36:01       34 阅读
  2. json转excel

    2024-05-04 07:36:01       34 阅读
  3. DRF路由组件分析

    2024-05-04 07:36:01       30 阅读
  4. docker搭建redis集群三主三从

    2024-05-04 07:36:01       31 阅读
  5. 【QT】QStackedWidget一利器,常在项目中的使用

    2024-05-04 07:36:01       33 阅读
  6. AWTK 和 QT 资源占用不完全对比

    2024-05-04 07:36:01       34 阅读
  7. 设计模式-02 设计模式-工厂模式factory

    2024-05-04 07:36:01       29 阅读
  8. 2024牛客五一集训派对day2

    2024-05-04 07:36:01       26 阅读