matlab提取多边形区域_[转载]【转】matlab中地图边界与掩膜(去掉边界外区域)的实现（基于shape...

首先是论坛里的一些资源，其实绘制地理类的图形，我一般都推荐用grads、ArcGis(不是很熟悉就不提NCL、surfer和MeteoInfo等等啦，论坛里关于这几个软件的帖子也很多，不一一举例了)~有关这两个软件底图制作和绘制在论坛里面有很详细的教程，大家可以参考(grads中maskout方法，ArcGis底图制作，ArcGis站点绘图参考)，而有关matlab底图方法，我估计m_map工具箱中会有一些解决方法，但是也不是很熟悉，我就用matlab自带的map工具箱中的函数进行处理了(不清楚各版本有没有装，改天看到工具箱下载的时候我共享一下)。

一、预备知识：

1.地理参考(R)：地理参考(具体内容可在help大部分map工具箱中的函数得到)是指定位矩阵用的矩阵，通常为3*2的矩阵，通过公式

[lon lat] =

[row col 1] * R(3*2) (具体实例见第三部分中程序框)

对矩阵的行列坐标与经纬度进行转换(也有1*3的形式，有兴趣的朋友可以自己了解一下)，是map工具箱处理地理数据的基础；

2.矢量数据：一般矢量数据分为多边形、线和点数据，我们这里采用的shape文件的多边形数据是利用线数据定义边界的，在matlab里shape文件以结构数组的形式保存，其中经纬度信息保存在X与Y属性中；

3.2-D矩阵数据：由行标(row)和列标(col)来表示位置的数据，需要通过坐标系来对数据进行对应进行绘图，常用的绘图参考方式有meshgrid函数生成的参考系和map中的地理参考(R)；

4.matlab中NAN(空值)的运算规则：NAN值除