MATLAB中自定义栅格数据地理坐标并写出数据到tiff文件

MATLAB中自定义栅格数据地理坐标并写出数据到tiff文件

MATLAB中自定义栅格数据地理坐标并写出数据到tiff文件

场景描述：

问题描述：

原因分析：

解决方法：

场景描述：

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

问题描述：

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

图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文件中

原因分析：

在调试完代码后，发现问题出现在没有添加rot90()这个函数。在matlab中用imshow()显示出的结果正确，那说明在该矩阵左上角的坐标是(90°，-180°)，右下角的坐标是(-90°，180°)。而空间坐标是从-90到90维度，-180到180经度，我猜测georasterref()函数内部是按照矩阵A的规格构建了一个坐标矩阵，如图3：

图3. 坐标矩阵

如果按照原始代码运行，那意味着，矩阵第一行维度是-90°，但是实际上应该是90°，所以输出的数据南北维度出现了颠倒的现象。

解决方法：

在完成了矩阵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文件中

MATLAB中自定义栅格数据地理坐标并写出数据到tiff文件相关教程

37.自定义导航栏开发 uni-app

37.自定义导航栏开发 uni-app 参考：https://uniapp.dcloud.io/collocation/pages?id=app-plus 禁止原生导航 1.新建发布页面add-input,在pages里写 2.在index.html写原生导航按钮的监听事件函数 onNavigationBarButtonTap(e){ console.log(JSON.stringify(e)

自定义属性操作/案例:自定义tab栏切换/H5自定义属性

自定义属性操作/案例:自定义tab栏切换/H5自定义属性 自定义属性操作 div id = 'demo' /div script var div = document.querySelector(div); //1.获取元素的属性值 //(1)element.属性 console.log(div.id); //(2)element.getAttribute('属性') console.log(d

MATLAB绘制纹理柱状图

MATLAB绘制纹理柱状图 https://ww2.mathworks.cn/matlabcentral/fileexchange/24021-hatch-fill-patterns-plus-color-invert 需要下载matlab包 applyhatch_pluscolor.m 是用来绘制的 makehatch_plus.m 当中包含了填充纹理的样式 testforapplyhatch.m 是用来

计算机概论之计算

计算机概论之计算 计算机概论计算篇 这里写自定义目录标题 进制运算的基础 进制概述 二进制运算的基础 二进制转换十进制：按权展开法 十进制转换二进制 有符号数与无符号数 原码表示法 二进制的补码表示法 二进制的反码表示法 定点数与浮点数 定点数的表示方

直纹面的简要介绍 以及matlab程序实现

直纹面的简要介绍 以及matlab程序实现 直纹面 是一类特殊的曲面，它可以由一族直线“织成”，即过曲线上每一点都存在过该点的直线落在该曲面上。 直纹面在 几何造型 中的应用非常广泛。直纹面是直线段在空间中沿某一定曲线运动所形成的轨迹. 定曲线称为准线，

mapreduce练习8过滤日志自定义输出

mapreduce练习8过滤日志自定义输出 1.FilterMapper import java.io.IOException;import org.apache.hadoop.io.LongWritable;import org.apache.hadoop.io.NullWritable;import org.apache.hadoop.io.Text;import org.apache.hadoop.mapreduce.Mapper; public

Vue 自定义组件directive使用总结

Vue 自定义组件directive使用总结 !DOCTYPE htmlhtml lang=enhead meta charset=UTF-8 meta name=viewport content=width=device-width, initial-scale=1.0 script src=https://cdn.jsdelivr.net/npm/vue/dist/vue.js/script title自定义组件directive/title/

MATLAB求解峰值、均值、平均幅值、方差、均方根、方根幅值、最值

MATLAB求解峰值、均值、平均幅值、方差、均方根、方根幅值、最值、峭度、斜度、裕度指标、波形指标、脉冲指标、峰值指标、峭度指标 MATLAB求解峰值、均值、平均幅值、方差、均方根、方根幅值、最值、峭度、斜度、裕度指标、波形指标、脉冲指标、峰值指标、峭