matlab line函数_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, clcz = peaks(1000);lon = [60 150];lat = [0 60];[LON,LAT] = meshgrid(linspace(lon(1),lon(2),1000), linspace(lat(1), lat(2),1000));figurem_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')figurem_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');figurec = 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参数要一致。

此次没有设置付费阅读，直接公开源代码。点击

好了，收工。

- End -

今日推荐

推荐｜MATLAB nctoolbox工具箱中文说明

又是一份来自Github的神器，帮你搞定MATLAB