Python爬取全国水系shp数据----文末附数据下载地址

本文分享的是如何用Python从中央台的全国雷达拼图中提取出纯净的水系图层。实现过程中会较多地使用numpy的向量运算技巧，比较适合有扎实的numpy基础的朋友进行综合练习

依赖包

要运行本文的Python程序，需要预先安装以下依赖，使用pip即可安装。

numpy
matplotlib
scipy
webcolors

一些先验信息

在正式开始coding之前，我们需要对图片的一些信息有所了解，首先就是dbz的颜色值，也就是原图右下角图例上显示的13个色块所表示颜色的编码，我们必须要知道这个编码，才能将dbz与底图分离。获取该颜色编码的方式有很多，比如你把图片用Photoshop打开，然后使用取色器直接从色块中取色即可获得它的hex码，或者安装一些小工具（比如MacOS上有一个小软件叫做ColorSlurp）使用取色器取色。由于我已经将色码取出，在这里就可以直接使用：

COLORS = [
    '#AD90F0',
    '#9600B4',
    '#FF00F0',
    '#C00001',
    '#D60100',
    '#FF0200',
    '#FF9000',
    '#E7C000',