前几个专栏我们讨论了几个不规则的投影转换问题,有需要的可以阅读以下文章:
这次遇到的问题是一个墨卡托投影的数据,需要转成标准的WGS经纬度。
使用的数据是HYCOM输出的全球海洋分层温度数据,下载地址:
这个数据有一个问题,就是使用的是墨卡托投影,即在南北极变形特别大。见下图:
在读取数据的过程中,我发现其中的经度明显不正常,即有些超过360。
仔细研究后发现,是因为经度需要进行取模处理,即 modulo = '360 degrees'。
ncdisp显示的变量列表:
因此采用以下的代码进行运算,大于360的,对360度取模,并进行重采样处理,得到规则的经纬度格网。
matlab代码:
结果图
♥欢迎点赞收藏♥