【QGIS从shp文件中筛选目标区域导出为shp】

木叶清风666

已于 2024-03-25 21:22:44 修改

阅读量2k

点赞数 8

分类专栏：地理信息数据处理文章标签： arcgis

于 2024-03-25 19:30:00 首次发布

本文链接：https://blog.csdn.net/qq_38734327/article/details/137003663

版权

地理信息数据处理专栏收录该内容

25 篇文章 1 订阅

订阅专栏

本文介绍了如何在QGIS中通过手动点选或高级过滤功能从shp文件中提取目标区域，并结合Python的regionmask库，生成用于研究的mask数据。

摘要由CSDN通过智能技术生成

1、写在前面

利用shp文件制作研究区域mask，Matlab版本，请点击 Matlab
利用shp文件制作研究区域mask，Python版本，请点击 Python

2、QGIS将shp文件中目标区域输出为shp

2.1、手动点选

将shp文件拖动到QGIS软件中，或者从界面【工程】中打开相应shp文件，结果如下图所示：
点击工具栏【切换编辑模式】进入编辑状态—>点击【框选或单击选择要素】选中目标区域—>在【图层】中右键【导出】—>【选中的要素另存为】，在弹窗中选择ERSI形状文件格式，选择导出路径即可完成，

2.2、高级过滤

除了手动点选之外，可以使用过滤功能进行快速选择，选择图层，右键选择【过滤】。在过滤窗口中选择合适的过滤条件，点击ok，(不知道为啥，我这个点击ok过滤之后，界面显示为空白，但是不影响结果输出)，然后按照上述导出流程进行导入即可。

3、上述shp完成后，配合python的shp文件，即可凸显研究区域了

import regionmask
import numpy as np
import geopandas as gpd

file= "./china2.shp"
countries = gpd.read_file(file)
lon =np.linspace(70,140,7000)
lat =np.linspace(15,60,4500)
mask= regionmask.mask_geopandas(countries, lon, lat).to_numpy()
mask[~np.isnan(mask)]=1