QGIS提取区域边界

最新推荐文章于 2023-10-14 16:43:57 发布

彭博锐

最新推荐文章于 2023-10-14 16:43:57 发布

阅读量2k

点赞数 7

文章标签：开源软件学习

本文链接：https://blog.csdn.net/pengbr/article/details/132611468

版权

QGIS提取区域边界——以华北平原为例

之前在我的博客中已经实现了华北平原的地理配准QGIS地理配准，这次我们进一步提取出华北平原的边界。

首先，我们要创建点shp文件，具体如下图：

接着我们需要对点shp文件进行编辑，具体如下图：
在这里插入图片描述
最后点击ok，就可以得到一个点shp文件。接着对point图层进行编辑。

点了一个点后，发现需要我们输入要素属性，先暂时不用管它，后面会一起调整，直接点击ok。这里就点5个点作为示例。

然后我们点击打开属性表，对之前的要素属性做一些修改。
在这里插入图片描述
打开属性表后会出现如下界面，我们需要添加经纬度字段并更新id。添加纬度字段的表达式是 $y ，添加经度字段的表达式是$ x。

结果如下：

接下来我们可以导出要素属性：

下面是我之前点的所有点的结果，也按照上述方式获取了所有点的经纬度：
在这里插入图片描述

利用GMT绘制华北平原边界：

#!/usr/bin/env -S bash -e
set -e
# GMT modern mode bash template
# Date:    2023-08-31T22:01:33
# User:    br
# Purpose: Purpose of this script
export GMT_SESSION_NAME=$$	# Set a unique session name
gmt begin figurename
	# Place modern session commands here
	gmt set MAP_GRID_PEN_PRIMARY 0.25p,gray,2_2
    # 绘制区域地图
    gmt coast -JM5c -R112/120/34/42 -Bafg -G244/243/239 -S167/194/223
    gmt plot CN-border-La.gmt -W0.1p

    # 绘制华北平原，我将经纬度数据存储为了txt文件，经度列在前，纬度列在后
    gmt plot huabeipinyuan.txt -W0.5p,red
gmt end show