专栏首页气象杂货铺MATLAB完美白化,完美支持m_map工具箱!

MATLAB完美白化,完美支持m_map工具箱!

最初的matlab完美白化的脚本是3年多之前在气象家园论坛发布的。当时是我们大师姐想要用matlab白化,然后就有了最初的版本,但最初的版本不支持m_map。虽然最后给了一个支持m_map的思路,但是没有给出具体实现的完整代码。

最近又把代码整理了一下,添加了对m_map的支持,几乎完美的支持m_map进行地图投影的白化。对函数代码进行了修改,并按照m_map工具箱函数的命名方式重命名,同时添加了一个 m_mapshow 函数,用于在 m_map 地图投影添加shapefile地图。

关于白化的思路就不多说了,文末可获取支持m_map白化的代码及注意事项

以下是此次测试的代码,仅测试了lambert、mercator、lat-lon投影的白化,lambert和mercator投影基于m_map工具箱,lat-lon仅是简单的常规白化。

clear, clc

z = peaks(1000);
lon = [60 150];
lat = [0 60];
[LON,LAT] = meshgrid(linspace(lon(1),lon(2),1000), linspace(lat(1), lat(2),1000));

figure
m_proj('lambert', 'longitudes', lon, 'latitudes', lat, 'par', [30, 60], 'clo', 105)
m_contourf(LON, LAT, z);
m_maskmap('../data/chinamap/中国行政区_包含沿海岛屿.shp', true, 'lon', lon, 'lat', lat, 'm_map', true);

m_proj('lambert', 'longitudes', lon, 'latitudes', [15, 60], 'par', [30, 60], 'clo', 105)
m_grid('box','on', 'linestyle', 'none', 'tickdir', 'out', 'linewidth', 3);
m_mapshow('../data/chinamap/cnmap/cnhimap.shp')

figure
m_proj('mercator', 'longitudes', lon, 'latitudes', lat)
m_contourf(LON, LAT, z);
m_maskmap('../data/chinamap/中国行政区_包含沿海岛屿.shp', true, 'lon', lon, 'lat', lat, 'm_map', true);
m_mapshow('../data/chinamap/cnmap/cnhimap.shp')
m_grid('box','on');

figure
c = contourf(LON, LAT, z,'linestyle', 'none');
m_maskmap('../data/chinamap/中国行政区_包含沿海岛屿.shp', true, 'lon', lon, 'lat', lat);
mapshow('../data/chinamap/cnmap/cnhimap.shp', 'color', 'k', 'displaytype', 'line')

以下是三个投影的白化测试结果,可以看出,完美白化

从左到右分别为Lambert、Mercator、lat-lon投影的白化图

此脚本在用m_map处理lambert投影的时候存在一点问题,如下图所示。这个暂时没有很好的方法解决。

但是在上述的测试脚本里,给出了简单的解决方案,即扩大坐标轴的范围,然后进行白化,白化完成之后,重新设置投影,最后再设置m_grid类型。这样就能解决上面的小问题了。

lambert投影白化小瑕疵

注意事项:

1)此函数仅支持shapefile文件类型是polygon的文件,如果是line,需要转换为polygon类型;

2)利用lambert投影白化时要注意前后两个投影的par和clo参数要一致。

此次没有设置付费阅读,直接公开源代码。点击阅读原文获取源码和shapefile文件。

好了,收工。

  • End -

本文分享自微信公众号 - 气象杂货铺(meteogs),作者:bugsuse

原文出处及转载信息见文内详细说明,如有侵权,请联系 yunjia_community@tencent.com 删除。

原始发表时间:2020-05-02

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

我来说两句

0 条评论
登录 后参与评论

相关文章

  • 网格气象场插值-NCL版

    通常所说的regridding/remaping/interpolation都是将不同网格的数据映射到新的网格。

    bugsuse
  • 为什么气象学家应该关注生成对抗网络

    近期GRL上发表了一篇利用生成对抗网络(Generative Adversarial Networks,GANs)重构云垂直结构的文章。文中使用MODIS观测资...

    bugsuse
  • 使用 Basemap 和 Cartopy 绘制子图实例

    平时绘制地图时,经常会将多个图放到同一个 figure 中,而这些图的地图范围通常是相同的,所以可以设置共享 x-y 轴。

    bugsuse
  • Xcode不显示历史记录的解决方案

    LinXunFeng
  • 学界 | CMU与谷歌新研究提出文本跳读方法,速度可达标准序贯LSTM的6倍

    选自arXiv 机器之心编译 参与:吴攀 因为很多文本内容都含有一些与主题不相关的东西,所以让机器学会根据上下文进行跳读可以大大节省文本处理的时间和效率。近日...

    机器之心
  • 为什么今天的企业研发中心要设在城市

    本文发表在哈佛商业评论上 亚特兰大中城区已经发生了很大变化。乔治亚科技城周边的校园已经成为美国企业研发中心的圣地。在过去十年中,科技广场,这个位于中城区的区域...

    点滴科技资讯
  • 谨防云计算的隐性成本

    对于企业来说,公共云的总体成本可能具有吸引力,但还有许多其他因素需要包括在内,如果没有仔细的规划,采用云计算的成本很快就会螺旋上升。

    静一
  • 四月小结

    致码DevOps
  • 吴恩达新书《机器学习训练秘籍》六大要点总结 | 资源

    近日,一位名叫Niklas Donges的工程师小哥发表了一篇博客文章,提炼出了《机器学习训练秘籍》中的六大要点。

    量子位
  • 困惑技术人的17个问题,现在都有答案了!

    通过“从软件外包到阿里技术专家再到CTO,他究竟是如何一路晋升?”一文,许多人了解到了环球黑卡CTO崔涛扎实的晋升路。

    养码场

扫码关注云+社区

领取腾讯云代金券