使用shp文件对数据进行mask

iovedys · 发表于 2021-5-15 13:47:22

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 iovedys 于 2022-3-18 10:21 编辑

小白画图，研究了很久如何抠出一个省的数据，网上很多教程是先画图再对地图白化而不是处理数据，笔者研究了很久找到一个可以用shp来直接对数组进行mask的方法。

主要用到的库：使用geopandas来读取shp文件，再使用salem库来mask数据。

这两个库的安装比较麻烦，笔者参考以下链接后成功安装geopandas和salem，总结如下：

首先使用 pip uninstall (包名）清理Shapely、GDAL、Fiona、pyproj、rasterio、geopandas、regionmask、salem。

从这里Python Extension Packages for Windows - Christoph Gohlke (uci.edu)下载自己python版本对应的上述包的whl文件，如64位Windows和python3.8需要对应下载Shapely-1.7.1-cp38-cp38-win_amd64.whl等(放在文末附件中）。

下载好后使用 pip install D:\xxx\xxx.whl 依次安装Shapely、GDAL、Fiona、pyproj、rasterio。

最后使用 pip install geopandas ， pip install regionmask， pip install salem。

注：第一次 import salem 会自动下载一个约40M的salem-example文件，需要耐心等待。

程序中以某月降水为例，可以看到江苏省省外的数据都被mask成了nan值：
QQ截图20210515132844.png

绘图效果：

MaskoutProvince.rar (45.94 MB, 下载次数: 1414)

灭火器 · 发表于 2021-5-16 08:34:05

学习一下。不过我更好奇有没有什么原生的、不用安装太多库的方法。

VV09 · 发表于 2021-5-15 14:20:24

不错，学习了~~

kingjames · 发表于 2021-5-15 16:22:47

优秀应该和maskout是异曲同工之妙

khaled · 发表于 2021-5-15 16:24:59

学习学习

{:eb512:}

随缘 · 发表于 2021-5-15 17:39:12

不错，学习了~~

tulalang · 发表于 2021-5-16 07:14:17

来看看学习一下

墨家大宝 · 发表于 2021-5-16 09:05:59

灭火器发表于 2021-5-16 08:34
学习一下。不过我更好奇有没有什么原生的、不用安装太多库的方法。

beiruoyiqiao · 发表于 2021-5-16 11:42:38

感谢分享~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

zhuyidong · 发表于 2021-5-16 17:38:39

不错，学习了。。。。。

		自动登录	找回密码
密码			立即注册

[经验总结] 使用shp文件对数据进行mask