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 -
本文参与 腾讯云自媒体分享计划，分享自微信公众号。
原始发表：2020-05-02，如有侵权请联系 cloudcommunity@tencent.com 删除
评论
登录后参与评论
推荐阅读