GIS数据文件的抽稀利器mapshaper20140902

今天真有点小激动，因为发现了这个GIS数据文件---目前针对shape格式的文件，但这已经足够---的抽稀利器mapshaper。

在早几年的博客当中我多次提到对shape文件的抽稀问题的迷惑，这迷惑既包括对抽稀理论和算法上的迷惑，也包括对没有趁手的抽稀工具的无奈。开始请教过一位GIS的博士，但可惜术业有专攻的博士对这个问题没有研究，然后自己写程序对shape文件抽稀，但缺乏理论的支持结果并不理想，再然后使用ArcGIS中的抽稀工具，效果堪用，但仍然觉得不够理想。

今天偶然的机会遇见mapshaper，我想问题可以就此终结。颇有“踏破铁鞋无觅处，得来全不费功夫”的感触，又觉得也有“暮然回首那人却在灯火阑珊处”的意境。

先解释一下什么是抽稀，百度百科是这么解释的：通过某种规则，在保证矢量曲线形状不变的情况下，最大限度地减少GIS数据中点、线、多边形等元素的个数，这个过程称为抽稀。

我们还是直观地看一下抽稀是怎么回事，下面的截图是在GIS软件中打开的从国家基础地理信息中心下载的中国轮廓图boul_4p.shp文件：

这幅地图从GIS的角度来说没有问题，但从地理教学的角度来说，就很不美观，细节过多，如果再放大一点就更难看了。

我们可以使用ArcGIS把这个shape文件进行抽稀操作，但是前面说了，我觉得结果不太理想。

下面我们使用mapshaper来进行抽稀操作。mapshaper是一个网页应用程序，因此我们需要在浏览器中输入该程序的地址：mapshaper.org，界面如下：

一个简单明了的界面，也正是我喜欢的类型，用户只需要点击“select”按钮选择需要抽稀的shape文件就来到下个界面：

这个操作界面同样简洁明了，用户只需要移动抽稀百分比的滑竿就可以了，这相当于支持对shape文件的无级抽稀，也是mapshaper强大之处。当选择好抽稀百分比后，点击“Export”后面的文件类型，就可以得到经过抽稀后的GIS数据文件，可以仍然是shape格式，也可以是新近流行的GeoJSON格式或者TopoJSON格式。这点也值得欢呼，我正在寻找一个shape文件到GeoJSON文件的转换工具，天上就掉下来了馅饼！

下面是我再次在GIS软件中打开抽稀后的中国轮廓地图：

对比前面的图，是不是整洁干净了许多？

mapshaper是一个开源软件，采用流行的Javascripts语言开发，需要Node.js的支持，其官方开发主页是：

https://github.com/mbloch/mapshaper

另外，这篇文章对mapshaper作了详细介绍，无论是理论还是实践应用都有很高的参考价值：

http://www.cartogis.org/docs/proceedings/2006/bloch_harrower.pdf