QGIS应用与实用插件使用指南-CN

本指南旨在为中国范围内的QGIS用户提供一份详尽的操作手册，重点介绍在处理行迹追踪、分布图制作、自然地貌与地形分析、流域绘制以及高级文字标注（特别是针对中文环境）等任务时，实用且高效的QGIS功能及插件。指南将尽可能清晰地描述操作过程，确保用户能够复现。

I. QGIS概述与基础

A. QGIS简介

QGIS (曾用名 Quantum GIS) 是一款功能强大、用户友好、自由且开源的地理信息系统 (GIS) 软件 (1)。它支持在多种操作系统上运行，包括Windows、macOS和Linux，并且能够处理和分析各种栅格与矢量地理空间数据格式 (1)。QGIS不仅提供数据查看、编辑、分析和制图打印等核心GIS功能，还因其开源特性而拥有一个充满活力的全球开发者和用户社区。这个社区不断贡献新的功能、修复错误，并开发了海量的插件，极大地扩展了QGIS的应用范围和能力 (2)。对于中国用户而言，QGIS的免费和开源特性使其成为一个极具吸引力的选择，无需高昂的软件采购成本即可获得专业的GIS处理能力。同时，其开放性也为本地化开发和功能定制提供了可能，使其能更好地适应中国用户的特定需求。

B. QGIS版本解析

QGIS拥有规律性的版本发布周期，通常每四个月会发布一个最新版本 (Latest Release)，而每年会指定一个版本作为长期支持版 (Long Term Release, LTR) (3)。截至编写本指南时（参考 3 的版本信息，具体版本号请查阅QGIS官网获取最新情况），用户可以根据需求选择：

• 长期支持版 (LTR): LTR版本更注重稳定性，在发布后会持续获得错误修复，但不会引入新的功能。这使得LTR版本成为生产环境、大型项目或对软件稳定性有较高要求的用户的理想选择。例如，QGIS 3.40.7 'Bratislava' 是一个LTR版本 (3)。
• 最新版本 (Latest Release): 最新版本包含了QGIS开发团队最新的功能和改进。它适合那些希望体验和使用最新GIS技术、对新功能有需求的用户或开发者。例如，QGIS 3.42.3 'Münster' 是一个最新版本 (3)。

QGIS的版本号中，偶数版本号（如3.40, 3.42）代表稳定发布版，而奇数版本号（如3.41, 3.43）则代表开发中的版本 (3)。这种规律性的发布和版本管理策略，体现了QGIS项目管理的成熟度和持续的开发活力。对于依赖特定插件的用户而言，选择版本时需要考虑插件与QGIS核心版本的兼容性。通常，插件开发者会优先确保其插件与LTR版本的兼容性，并逐步适配最新的QGIS版本。

C. 安装指南与系统要求

用户可以从QGIS官方网站 (qgis.org) 下载适用于Windows、Linux、macOS等操作系统的安装程序 (1)。安装过程通常都比较直观，按照安装向导的提示即可完成。

系统要求 (5):

• 最低配置:
  • 处理器 (CPU): Core i3 2.7 GHz
  • 内存 (RAM): 2 GB
  • 硬盘: 500 GB SATA
  • 显卡: 1 GB RAM
  • 操作系统: Windows 7-10, macOS, Linux
• 推荐配置:
  • 处理器 (CPU): Core i7 3.5 GHz
  • 内存 (RAM): 8 GB 或更高
  • 硬盘: SSD 128 GB 或 500 GB
  • 显卡: 2 GB RAM (NVIDIA GeForce)
  • 操作系统: Windows 7-10 (64位), macOS, Linux

值得注意的是，自QGIS 3.20版本起，Windows平台仅提供64位安装包 (4)。虽然QGIS对硬件要求相对适中，但在处理大规模栅格数据、进行复杂空间分析或三维可视化等资源密集型任务时，拥有接近或超过推荐配置的硬件将显著提升软件的运行效率和用户体验 (6 中提及QGIS有时可能对系统资源需求较高)。

D. QGIS用户界面导览

启动QGIS后，用户将看到一个包含多个核心组件的用户界面 (2)：

• 菜单栏 (Menu Bar): 位于界面顶部，包含QGIS所有功能的分类菜单，如项目、编辑、视图、图层、设置、插件、矢量、栅格、数据库、Web、处理和帮助等。
• 工具栏 (Toolbars): 位于菜单栏下方，提供常用功能的快捷图标按钮，用户可以根据需要自定义显示哪些工具栏。
• 浏览器面板 (Browser Panel): 通常位于界面左侧，用于浏览、查找和管理计算机上的地理空间数据文件、数据库连接以及各种在线服务。
• 图层面板 (Layers Panel): 通常位于浏览器面板下方或界面左侧，显示当前项目中加载的所有图层，用户可以在此控制图层的显示顺序、可见性、样式等。
• 地图画布 (Map Canvas): 占据界面中央的主要区域，用于显示和交互地理数据。
• 状态栏 (Status Bar): 位于界面底部，显示当前鼠标指针的坐标、地图比例尺、投影信息以及一些操作提示。
• 处理工具箱 (Processing Toolbox): 通常以面板形式停靠在界面右侧，集成了QGIS内置算法以及来自GDAL, SAGA GIS, GRASS GIS等第三方提供者的分析工具。

熟悉这些核心组件的功能和布局，是高效使用QGIS进行地理空间任务处理的基础。

II. 针对中国地理空间任务的QGIS配置

在中国进行地理空间数据处理与分析，首要任务之一是正确配置QGIS以适应本地的坐标参照系统和常用的在线地图服务。

A. 中国常用坐标参照系统 (CRS)

正确选择和使用坐标参照系统 (CRS) 是保证地理空间数据精度、实现多源数据正确叠加分析以及地图准确发布的基石。在中国，由于历史和技术发展原因，存在多个国家级和地方性的坐标系统。以下是在QGIS中进行中国相关项目时需要重点关注的CRS：

1. 国家大地坐标系2000 (CGCS2000): 这是中国现行法定的国家大地坐标系，推荐在新的测绘项目和数据生产中使用。
   • 地理坐标系: EPSG:4490 (China Geodetic Coordinate System 2000) (7)
   • 大地坐标系 (三维): EPSG:4479 (CGCS2000) (7)
   • 投影坐标系: CGCS2000通常采用高斯-克吕格投影。根据分带不同，EPSG代码也不同，例如：
     • CGCS2000 / Gauss-Krüger zone 20 (适用于东经114°至120°之间的6度带): EPSG:4498 (9)
     • CGCS2000 / 3-degree Gauss-Krüger CM 117E (适用于中央经线为东经117°的3度带): EPSG:4548 (基于 8 中CGCS2000作为基准的3度带命名规律推断)
2. 北京1954坐标系 (Beijing 1954): 这是中国曾长期使用的一个国家大地坐标系，许多历史数据和图件仍采用此系统。
   • 地理坐标系: EPSG:4214 (Beijing 1954) (7)
   • 投影坐标系: 北京1954同样采用高斯-克吕格投影。例如：
     • Beijing 1954 / Gauss-Kruger zone 13 (适用于东经72°至78°之间的6度带): EPSG:21413 (7)
     • Beijing 1954 / 3-degree Gauss-Kruger zone 25 (适用于中央经线为东经75°的3度带): EPSG:2401 (11)
     • Beijing 1954 / 3-degree Gauss-Kruger CM 117E (适用于中央经线为东经117°的3度带): EPSG:2436 (12)
3. 西安1980坐标系 (Xian 1980): 这是在CGCS2000之前的另一个国家大地坐标系，也广泛应用于各种存量数据中。
   • 地理坐标系: EPSG:4610 (Xian 1980) (7)
   • 投影坐标系: 西安1980也采用高斯-克吕格投影。例如：
     • Xian 1980 / Gauss-Kruger zone 18 (适用于东经102°至108°之间的6度带): EPSG:2332 (14)
     • Xian 1980 / 3-degree Gauss-Kruger CM 105E (适用于中央经线为东经105°的3度带): EPSG:2359 (基于西安1980的3度带命名规律)

高斯-克吕格投影带 (Gauss-Krüger Zones):
为了控制投影变形，中国采用高斯-克吕格分带投影。主要有6度带和3度带两种：

• 6度带: 从东经1°30′起，每隔6°为一带，全球共分60带。中国境内主要涉及第13带至第23带。带号N的中央经线计算公式为 L0​=6∘N−3∘。
• 3度带: 从东经1°30′起，每隔3°为一带，全球共分120带。中国境内广泛使用。带号N的中央经线计算公式为 L0​=3∘N。3度带通常用于较大比例尺的测图 (如1:1万至1:50万)。(15 展示了北京1954坐标系下3度带的EPSG示例)。

选择与转换的重要性:
在QGIS中处理来自不同来源、不同时期的中国地理空间数据时，务必首先确认各数据层所采用的坐标参照系统。如果数据没有定义CRS，需要根据元数据或数据来源进行设置；如果数据CRS不一致，则必须进行精确的坐标转换，才能保证空间分析的准确性和地图叠加的正确性。QGIS内置了强大的坐标转换功能，并支持广泛的EPSG代码。用户在项目属性和图层属性中均可设置和查看CRS。对于涉及国家标准和地方标准的复杂转换，可能需要特定的转换参数或插件辅助。
下表总结了中国常用的坐标参照系统及其EPSG代码，供用户快速参考：

表1: 中国常用坐标参照系统及其EPSG代码

B. 添加中国在线地图服务 (Adding Chinese Basemap Services via XYZ Tiles)

在线地图服务为GIS项目提供了重要的地理背景和参考。QGIS通过“XYZ Tiles”功能，可以方便地加载各类在线瓦片地图服务。对于中国用户而言，常用的国内地图服务提供商包括天地图、高德地图和腾讯地图。

手动配置XYZ连接的通用步骤 (17)：

1. 在QGIS的“浏览器”面板中，找到并右键点击“XYZ Tiles”。
2. 选择“新建连接…” (New Connection…)。
3. 在弹出的“XYZ连接”对话框中：
   • 名称 (Name): 为该服务指定一个易于识别的名称 (例如：“天地图矢量图层”)。
   • URL: 输入该地图服务的XYZ瓦片URL。URL中通常包含 {x}, {y}, {z} (或类似 {q}) 作为瓦片坐标和缩放级别的占位符。
   • 最大缩放级别 (Max. Zoom Level): 根据服务提供商的说明设置，较高的值通常意味着更详细的地图。
   • 其他参数（如认证、参照系等）根据具体服务要求配置。
4. 点击“确定”保存连接。之后，该服务将出现在“XYZ Tiles”列表下，双击即可添加到地图画布。

常用中国在线地图服务XYZ URL参考：
由于在线服务的URL可能发生变化，以下提供的URL仅供参考，建议用户在使用前进行验证或查找最新的服务地址。

• 天地图 (TianDiTu) (20)

  • 矢量地图: http://t{0-7}.tianditu.gov.cn/DataServer?T=vec\_w\&x={x}\&y={y}\&l={z}\&tk=您的天地图密钥
  • 矢量注记: http://t{0-7}.tianditu.gov.cn/DataServer?T=cva\_w\&x={x}\&y={y}\&l={z}\&tk=您的天地图密钥
  • 影像地图: http://t{0-7}.tianditu.gov.cn/DataServer?T=img\_w\&x={x}\&y={y}\&l={z}\&tk=您的天地图密钥
  • 影像注记: http://t{0-7}.tianditu.gov.cn/DataServer?T=cia\_w\&x={x}\&y={y}\&l={z}\&tk=您的天地图密钥
  • 地形地图: http://t{0-7}.tianditu.gov.cn/DataServer?T=ter\_w\&x={x}\&y={y}\&l={z}\&tk=您的天地图密钥
  • 地形注记: http://t{0-7}.tianditu.gov.cn/DataServer?T=cta\_w\&x={x}\&y={y}\&l={z}\&tk=您的天地图密钥
  • 注意: 天地图服务通常需要申请并使用个人开发密钥 (tk参数)。请替换URL中的“您的天地图密钥”部分。{0-7} 表示可以使用 t0 到 t7 的任一子域名。

• 高德地图 (Gaode Maps / AutoNavi) (URL参考 18)

  • 矢量地图 (路网): http://wprd0{1-4}.is.autonavi.com/appmaptile?lang=zh\_cn\&size=1\&scale=1\&style=7\&x={x}\&y={y}\&z={z} (其中 {1-4} 表示 wprd01 到 wprd04)
  • 影像地图: http://wprd0{1-4}.is.autonavi.com/appmaptile?lang=zh\_cn\&size=1\&scale=1\&style=6\&x={x}\&y={y}\&z={z}
  • 影像注记 (混合图层): http://wprd0{1-4}.is.autonavi.com/appmaptile?lang=zh\_cn\&size=1\&scale=1\&style=8\&x={x}\&y={y}\&z={z} (此URL可能同时包含影像和路网)

• 腾讯地图 (Tencent Maps) (URL参考 22)

  • 矢量地图: http://rt{0-3}.map.gtimg.com/realtimerender?z={z}\&x={x}\&y={y}\&type=vector\&styleid=0 (其中 {0-3} 表示 rt0 到 rt3)
  • 影像地图: http://p{0-3}.map.gtimg.com/sateTiles/{z}/{x\_16\_dir}/{y\_16\_dir}/{x}\_{y}.jpg (腾讯影像的URL结构较复杂，{x_16_dir} 和 {y_16_dir} 通常是基于x,y计算的目录结构，可能需要专门的转换逻辑或使用特定插件)
  • 地形地图: http://rt{0-3}.map.gtimg.com/tile?z={z}\&x={x}\&y={y}\&type=terrain\&styleid=0

表2: 常用中国在线地图服务XYZ URL (示例)

使用插件加载天地图:
TianDiTu Tools 插件 (24) 为QGIS用户提供了一个便捷的方式来加载天地图的各种服务，包括矢量、影像、地形及其注记图层，并集成了地名搜索、地理编码等功能。安装该插件后，用户可以直接从插件菜单中选择并添加天地图服务，无需手动配置XYZ连接和管理API密钥（插件可能内置了公共密钥或引导用户配置）。
重要提示:

• 坐标偏移: 中国大陆的在线地图服务（如高德、腾讯，以及在中国大陆运营的谷歌地图）通常采用GCJ-02（“火星坐标系”）或BD-09（百度坐标系）等经过加密偏移的坐标系统。直接将WGS-84坐标系下的数据（如GPS轨迹）与这些底图叠加会导致位置不匹配。此时，需要使用坐标转换工具（见下一节）进行校正。
• 服务可用性与政策: 在线地图服务的URL和可用性可能随时发生变化，部分服务可能需要API密钥或有使用限制。用户应查阅服务提供商的最新文档。

C. 中国区实用工具插件

为了更好地处理中国特有的地理空间数据问题，特别是坐标系统转换的挑战，QGIS社区提供了一些非常实用的插件。

1. Chinese Coordinate Converter 插件 (25)
   • 功能: 该插件专注于解决中国大陆广泛使用的几种坐标系统之间的转换问题，包括：
     • WGS-84 (国际通用的GPS坐标系统)
     • GCJ-02 (国测局坐标系，又称“火星坐标系”，为高德地图、腾讯地图等国内主要地图服务商采用)
     • BD-09 (百度坐标系，在GCJ-02基础上进一步加偏)
   • 特点: 支持将矢量数据在上述坐标系之间进行相互转换，并可以将转换结果输出为多种常见格式，如临时图层、Shapefile、GeoJSON、KML和GeoPackage (25)。这对于整合来自不同源头（例如，WGS-84的GPS数据与高德地图的POI数据）的数据至关重要。
   • 安装与使用:
     1. 打开QGIS，通过顶部菜单栏“插件” -> “管理并安装插件…”。
     2. 在插件管理器中搜索“Chinese Coordinate Converter”并安装。
     3. 安装完成后，通常可以在QGIS的“矢量”菜单或“处理工具箱”中找到该插件的工具。
     4. 启动插件后，选择输入图层、源坐标系、目标坐标系以及输出格式和路径。
     5. 执行转换。
2. GeoHey Toolbox 插件 (26)
   • 功能: 这是一个工具集插件，其中包含了针对中国“火星坐标”的转换功能，支持在WGS-84、GCJ-02和百度坐标之间进行转换。
   • 特别说明: 该插件在其描述中明确指出：“因国家法律规定，不提供转WGS84的工具” (26)。这意味着它主要用于从WGS-84转换到GCJ-02或BD-09，或者在GCJ-02和BD-09之间互转，但可能不提供从加密坐标系（GCJ-02, BD-09）直接转换回未经偏移的WGS-84的功能。用户在使用时需注意此限制，并确保操作符合相关法律法规。
   • 安装与使用: 类似于Chinese Coordinate Converter，通过QGIS插件管理器安装和使用。

这些插件的存在极大地便利了在中国环境下使用QGIS的用户。由于GCJ-02和BD-09这类“加偏”坐标系统在中国在线地图和导航服务中的普遍应用，任何需要将国际标准的WGS-84数据（如野外采集的GPS点、国际遥感影像等）与这些服务进行精确叠加，或者在这些不同偏移标准的地图数据之间进行整合时，都离不开这类坐标转换工具。理解这些坐标系统的差异以及合规使用转换工具，是确保空间数据在中国区域内准确应用的关键环节。

III. 行迹追踪与分析

行迹追踪与分析是GIS应用的重要领域之一，广泛应用于交通管理、物流监控、野生动物行为研究、个人户外活动记录等场景。QGIS及其插件为此类任务提供了强大的支持。

A. 轨迹数据简介

轨迹数据记录了物体在一段时间内的空间位置序列。常见的轨迹数据格式包括：

• GPX (GPS Exchange Format): 一种基于XML的开放标准，常用于GPS设备和应用程序之间交换航点、轨迹和路径数据。
• KML (Keyhole Markup Language): 也是一种基于XML的格式，用于在Google Earth等地球浏览器中显示地理数据，可以包含带时间戳的轨迹。
• CSV (Comma-Separated Values) 或其他文本文件: 通常包含经度、纬度、时间戳以及可能的其他属性（如速度、高度等）列。

在QGIS中，轨迹数据通常被表示为点图层（每个点代表一个时空位置）或线图层（连接连续时空点形成的路径），其核心特征是每个要素都关联一个或多个时间属性。

B. 使用QGIS时态控制器可视化时序数据

QGIS从3.14版本开始内置了强大的“时态控制器” (Temporal Controller)，用于动态可视化和分析时间序列数据，它取代了早期版本中广受欢迎的TimeManager插件 (27)。

操作步骤 (27)：

1. 准备数据:
   • 确保您的图层（如导入的GPX轨迹点或带时间戳的CSV文件）包含一个或多个正确格式的时间字段。QGIS时态控制器能够识别多种时间格式，但标准格式如 YYYY-MM-DD HH:MM:SS (例如 2023-10-26 14:35:00) 通常具有最佳兼容性 (29)。
   • 如果轨迹数据是线要素，并且每条线代表一个完整的轨迹段，那么时间属性可能表示该轨迹段的开始和结束时间。
2. 启用图层时态设置:
   • 在“图层面板”中，右键点击需要进行时态可视化的图层，选择“属性…”。
   • 切换到“时态” (Temporal) 选项卡。
   • 勾选“启用动态时态控制” (Dynamic Temporal Control) 复选框。
   • 配置 (Configuration):
     • 如果每个要素只有一个时间戳（例如，一系列独立的GPS点），选择“单一字段含日期/时间” (Single Field with Date/Time)。然后在“字段” (Field) 下拉菜单中选择包含时间戳的属性列。
     • 如果要素有一个开始时间和结束时间（例如，一个事件的持续时间或一个轨迹段），选择“分别指定开始和结束日期/时间的字段” (Separate Fields for Start and End Date/Time)。然后分别指定开始时间和结束时间字段。
   • 事件持续时间 (Event duration) (可选): 如果适用，可以设置事件的持续时间。
   • 点击“应用”和“确定”。启用时态控制后，图层名称旁边会出现一个时钟图标。
3. 使用时态控制面板:
   • 从QGIS主菜单栏选择“视图” -> “面板” -> “时态控制面板”，或点击地图导航工具栏中的时钟图标 (Temporal control panel) 打开该面板。
   • 点击面板上的“动画时态导航” (Animated temporal navigation) 按钮（通常是一个播放图标）以激活时间滑块和动画控制。
   • 设置时间范围 (Range): 点击“设置为完整范围” (Set to full range) 按钮（通常是一个刷新图标），面板会自动根据图层的时间属性计算并设置整个动画的时间跨度。您也可以手动输入开始和结束时间。
   • 设置步长 (Step): 定义动画播放时每个时间步的长度（例如，1小时、1天、1周）。
   • 播放动画: 点击播放按钮开始动画。时间滑块会根据设定的步长移动，地图画布上只显示在当前时间窗内的要素。
   • 调整帧率: 如果动画播放过快或过慢，可以点击时态控制面板右上角的“时态设置” (Temporal settings) 按钮（黄色齿轮图标），调整“每秒帧数” (Frames per second) (27)。
4. 添加时态文本装饰 (可选) (27)：
   • 为了在地图上清晰显示当前动画的时间，可以添加一个文本装饰。
   • 选择主菜单“视图” -> “装饰” -> “标题注记…” (Title Label…)。
   • 勾选“启用标题注记”复选框。
   • 点击“插入表达式…”按钮，输入表达式以显示当前时间。例如，要显示当前年份，可以使用 format_date(@map_start_time, 'yyyy')。@map_start_time 是一个QGIS变量，代表当前时间片开始的时间。
   • 调整字体大小、颜色、背景和放置位置（如右下角）。
5. 导出动画 (27)：
   • 在时态控制面板中，点击“导出动画” (Export Animation) 按钮（通常是一个保存图标）。
   • 输出目录: 选择一个文件夹用于保存动画的图像序列（通常是PNG格式的图片，每帧一张）。
   • 范围: 可以选择导出整个地图画布范围，或根据某个图层的范围进行计算（例如，Calculate from Layer -> ne_10m_land）。
   • 点击“保存”。QGIS将按设定的时间步长逐帧渲染地图并保存为图片。
   • 生成GIF/视频: 导出的图像序列可以使用第三方工具（如在线的ezgif.com，或视频编辑软件）合成为GIF动画或视频文件。

时态控制器是QGIS中进行动态可视化和初步时空模式探索的核心功能。通过它，用户可以直观地观察地理现象随时间的变化，为后续更深入的分析打下基础。其从早期插件演变为核心功能，也体现了QGIS对时空数据处理能力的持续重视和改进。

C. 使用 Trajectools 插件进行高级轨迹处理

当需要对轨迹数据进行更复杂的分析和处理时，例如计算速度、方向、分割轨迹、平滑轨迹或进行更高级的移动模式挖掘，Trajectools 插件 (30) 提供了一套强大的算法，集成在QGIS的“处理工具箱”中。

安装与依赖:

• 核心依赖: Trajectools 插件的核心依赖是 MovingPandas 库，版本要求通常为 0.21.2 或更高（例如，32 提及的 0.22.3）。MovingPandas 是一个专门用于处理和分析移动数据的Python库。
• 可选依赖: 为了使用某些高级功能，Trajectools 还可能需要以下可选库：
  • scikit-mobility: 用于隐私相关的分析。
  • stonesoup: 用于轨迹平滑算法（如卡尔曼滤波）。
  • gtfs_functions: 用于处理GTFS（General Transit Feed Specification，通用公交数据规范）格式的公共交通数据。
• 安装方式:
  1. 推荐: 使用 conda 或 mamba 包管理器来安装这些Python依赖项，因为它们能更好地处理复杂的环境和库之间的兼容性问题 (32)。通常的步骤包括创建一个新的conda环境，激活该环境，然后安装 qgis、movingpandas 及其他可选库。
  2. 备选: 如果无法使用conda，可以尝试在QGIS的Python控制台中使用 pip 安装，但这可能更容易遇到依赖冲突 (32)。
  3. 插件安装: Trajectools 插件本身通过QGIS的“插件管理器”进行搜索和安装。

核心算法与用例 (基于 31 的信息)：
Trajectools 提供了一系列模块化的处理算法，可以灵活组合应用于各种轨迹数据集。主要功能包括：

• 创建轨迹 (Create Trajectory): 从包含时间戳和轨迹ID的点图层生成轨迹线图层。
• 分割轨迹 (Split Trajectory): 根据特定条件分割轨迹，例如：
  • 按观测间断 (Split by observation gap) (31)：当两个连续点之间的时间或空间间隔超过阈值时分割。
  • 按字段值变化 (Split at field value changes) (31)：当轨迹点序列中某个属性字段的值发生变化时分割。
• 轨迹平滑 (Trajectory Smoothing):
  • 卡尔曼滤波 (Kalman Filter Smoothing) (32)：一种常用的轨迹平滑和状态估计算法。
• 轨迹概括与简化 (Trajectory Generalization/Simplification): 减少轨迹中的冗余点，同时保留其主要形状特征。
• 计算轨迹属性:
  • 提取轨迹点或段的统计信息，如速度、加速度、方向、长度、持续时间等。
  • 可以将选定的字段添加到输出的轨迹图层中 (31)。
• GTFS数据处理:
  • 提取GTFS数据中的站点 (Extract stops) (31)。
  • 提取GTFS数据中的路段 (Extract segments) (31)。
  • 为GTFS站点和路段提供默认输出样式 (31)。

通用操作步骤:

1. 准备输入数据: 确保输入数据（通常是点图层）包含明确的轨迹标识符 (Track ID)、每个点的时间戳 (Timestamp)，以及坐标信息。如果是GPX或KML文件，需要先将其导入QGIS转换成标准的矢量图层。
2. 打开处理工具箱: 在QGIS中，通过菜单“处理” -> “工具箱”打开。
3. 查找Trajectools算法: 在处理工具箱中搜索或浏览找到 Trajectools 提供的算法。
4. 配置算法参数: 双击选定的算法，会弹出参数设置对话框。根据算法要求，指定：
   • 输入图层 (Input layer)。
   • 轨迹ID字段 (Track ID field)。
   • 时间戳字段 (Timestamp field)。
   • 其他特定参数 (例如，分割阈值、平滑参数、输出字段等)。
5. 运行分析: 设置完参数后，点击“运行” (Run)。
6. 查看结果: 分析完成后，结果通常会以新图层的形式添加到QGIS地图画布或图层面板中。

Trajectools 插件通过整合MovingPandas等强大的Python库，极大地增强了QGIS在轨迹数据分析方面的专业能力。它不仅能进行可视化，更能深入挖掘移动数据中蕴含的模式和信息。用户在使用时，应特别注意其依赖库的正确安装，并参考插件的最新文档以了解新增功能和具体用法。例如，Trajectools 2.4版本就新增了GTFS算法和改进了用户界面 (31)。

IV. 高效分布图制作

分布图是地理信息可视化中最常用的一类地图，用于直观展示地理现象的空间分布特征，如位置、密度、数量、频率或变化趋势。

A. 分布图制图原理与数据准备

制作一幅优秀的分布图，首先需要理解其核心目的，并对数据进行细致的准备。

• 目的明确: 是要展示离散点的位置，还是区域的统计值，或是某种现象的连续密度？不同的目的对应不同的制图方法。
• 数据质量: 原始数据往往需要清洗，包括处理缺失值、异常值，确保属性数据的准确性和一致性。
• 分类方法 (Classification): 对于需要将连续数据分级显示的地图（如分区统计图），选择合适的分类方法至关重要。常用方法包括自然间断点法 (Natural Breaks/Jenks)、等间距法 (Equal Interval)、分位数法 (Quantile)、标准差法 (Standard Deviation)等。选择不当的分类方法可能会误导地图读者。
• 视觉变量 (Visual Variables): 根据地图主题和数据类型，选择合适的视觉变量来表达信息，如符号的形状、大小、颜色（色相、饱和度、亮度）、方向、纹理等。

B. QGIS符号化技术的运用

QGIS内置了强大且灵活的符号化引擎，支持多种渲染方式来制作不同类型的分布图。这些功能通常在图层属性的“符号化” (Symbology) 选项卡中进行配置。

• 单一符号 (Single Symbol): 所有要素使用相同的符号，适用于仅展示要素位置的地图。
• 分类符号 (Categorized): 根据要素的某个类别属性（如土地利用类型、行政区等级）赋予不同的符号（颜色、形状等）。
• 分级符号 (Graduated): 根据要素的某个数值属性，将其划分为若干等级，并为每个等级赋予不同的符号样式（通常是颜色深浅或符号大小的变化）。这是制作分区统计图 (Choropleth maps) 和比例符号图的常用方法。
  • 值 (Value): 选择用于分级的数值字段。
  • 分类模式 (Mode): 选择上述提到的分类方法（自然间断点、等间距等）。
  • 类 (Classes): 设置分级数量。
  • 色彩渐变 (Color ramp): 选择合适的颜色方案。
• 比例符号 (Proportional Symbols): 符号的大小与要素的某个数值属性值成正比。QGIS中可以通过“分级”渲染器，并选择合适的方法（如“大小”）来实现。
• 点密度图 (Dot Density): 在面要素内部随机或按规则布设若干点，每个点代表一定的数量值，点的密度直观反映了该区域内某现象的分布强度。
• 热力图 (Heatmap): 基于点要素的密度生成连续的栅格表面，颜色深浅表示密度高低。非常适合展示点状事件的集聚区域 (27 中曾用热力图显示海盗事件的分布)。
• 规则为本的符号系统 (Rule-based): 允许用户定义一系列规则，满足特定规则的要素将应用相应的符号样式。这提供了极高的灵活性，可以组合多种条件进行复杂符号化。

操作步骤示例 (以分级符号制作分区统计图为例):

1. 在“图层面板”中右键点击目标矢量图层（通常是面图层），选择“属性…”。
2. 在“图层属性”对话框中，切换到“符号化” (Symbology) 选项卡。
3. 在顶部的下拉菜单中，选择“分级” (Graduated) 渲染器。
4. 在“值” (Value) 下拉菜单中，选择要用于分级的数值型属性字段（例如，人口密度、GDP等）。
5. 选择一个“色彩渐变” (Color ramp)。
6. 在“模式” (Mode) 中选择一个合适的分类方法（如“自然间断点 (Jenks)”）。
7. 设置“类” (Classes) 的数量。
8. 点击“分类” (Classify) 按钮，QGIS会自动计算分级范围并应用颜色。
9. 用户可以进一步调整每个等级的范围、标签和符号颜色。
10. 点击“应用”和“确定”。

QGIS的符号化系统还支持数据定义的覆盖设置 (Data-defined Override)，允许用户使用要素的属性值或表达式来动态控制符号的几乎所有参数（如颜色、大小、旋转角度、透明度等），这为创建高度定制化和信息丰富的地图提供了无限可能。此外，图层混合模式 (Blending Modes) 和绘制效果 (Draw Effects) 也能进一步增强地图的视觉表现力。

C. 使用 Data Plotly 插件制作高级统计图表

除了在地图上直接进行符号化表达，有时还需要结合统计图表来更全面地展示数据的特征和关系。Data Plotly 插件 (28) 使得用户可以在QGIS内部创建交互式的、出版质量的统计图表。

功能与特点:

• 基于Plotly库: 该插件利用了强大的Plotly Python库（其底层可能基于D3.js），能够在QGIS中生成多种类型的动态图表，如散点图、条形图、箱线图、直方图、饼图、三维图等 (34)。
• 交互性: 生成的图表是完全动态的，用户可以进行缩放、平移、悬停显示数据点信息等交互操作 (35)。
• 图表与地图联动: 这是Data Plotly的核心优势之一。在图表上选择一个或多个数据点（或条、块），地图画布上对应的地理要素会自动高亮显示；反之，在地图上选择要素，图表中对应的数据也会被突出。这种双向联动极大地增强了数据的探索性分析能力 (28)。
• 自定义与导出: 用户可以通过图形用户界面或代码对图表的美学特征（如标题、轴标签、颜色、图例等）进行广泛定制。图表可以导出为静态图片（如PNG格式）或包含完整交互性的HTML文件，方便在报告、网页中嵌入或分享 (35)。
• 离线使用: 插件本地保存了Plotly库的JavaScript代码，因此即使在没有互联网连接的情况下也能正常使用 (35)。

安装与使用:

1. 安装:
   • 通过QGIS的“插件” -> “管理并安装插件…”菜单，搜索“DataPlotly”并安装。
   • Data Plotly 可能依赖一些Python库，如 pandas 和 plotly (或 dash)。通常插件管理器会自动处理这些依赖，但如果遇到问题，可能需要手动在QGIS的Python环境中安装它们 (34)。
2. 界面: 安装后，Data Plotly 通常会在QGIS界面中添加一个专用的面板或工具栏按钮。
3. 操作步骤 (通用流程):
   • 打开Data Plotly插件界面。
   • 选择图层: 从当前QGIS项目中选择一个矢量图层作为数据源。
   • 选择图表类型: 根据分析需求选择合适的图表类型（例如，散点图用于展示两个变量的关系，条形图用于比较不同类别的数值，直方图用于显示数据分布）。
   • 配置图表参数:
     • 指定用于X轴、Y轴（以及可能的Z轴、颜色、大小等视觉通道）的属性字段。
     • 设置图表标题、轴标签、图例等。
     • 调整颜色、点型、线型等外观属性。
   • 生成与交互: 点击生成图表。图表会显示在一个新的窗口或面板中。利用其交互功能（如缩放、平移、悬停提示、选择高亮）与地图画布联动，进行数据探索。
   • 导出: 根据需要将图表导出为图片或HTML文件。

Data Plotly 插件将QGIS从一个传统的GIS制图工具，提升为一个更全面的空间数据可视化和探索性分析平台。它特别适用于需要深入理解数据内部结构、变量间关系，或希望制作具有高度交互性的可视化成果的场景。虽然它可以生成多种统计图表，但其核心价值更在于交互式的高级绘图，而非简单地替代静态的分布图符号化 (36)。在处理非常大的数据集时，其性能表现可能需要用户关注 (34)。

V. 自然地貌与地形分析

地形分析是GIS应用的重要组成部分，通过对数字高程模型 (DEM) 的处理，可以提取各种地形参数，理解地表形态特征，服务于水文、生态、灾害、规划等多个领域。

A. 中国区域数字高程模型 (DEM) 的获取与应用

数字高程模型 (DEM) 是地表高程的栅格化表达，是进行一切地形分析的基础。

• 常用DEM数据源:

  • 全球性开源/免费DEM:
  • SRTM (Shuttle Radar Topography Mission): 提供约30米和90米分辨率的全球DEM。可通过USGS EarthExplorer、OpenTopography (37) 等平台下载。
  • ASTER GDEM (Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model): 提供约30米分辨率的全球DEM。
  • Copernicus DEM (COP30, COP90): 由欧洲航天局发布的高质量全球DEM，分辨率分别为30米和90米 (37)。
  • ALOS AW3D30 (ALOS World 3D - 30m): 日本宇宙航空研究开发机构 (JAXA) 提供的全球30米分辨率DEM (37)。
  • 这些数据通常可以通过OpenTopography等门户网站方便地按区域选择和下载 (39)。
  • 商业/专门数据提供商:
  • NextGIS: 提供覆盖中国的30米分辨率DEM数据产品，以及其他相关的基础地理数据 (41)。
  • 中国本土数据源: 用户也可关注中国国家或地方测绘地理信息部门是否发布了可公开获取的DEM产品。
• 选择DEM时的考量:
  • 分辨率 (Resolution): 决定了DEM能表达的地形细节程度。高分辨率DEM能提供更精细的地形特征，但数据量也更大。
  • 精度 (Accuracy): 包括水平精度和垂直精度，影响地形分析结果的可靠性。
  • 覆盖范围 (Coverage): 确保所选DEM能完整覆盖研究区域。
  • 数据格式与获取方式: 常见的DEM格式有GeoTIFF等。
• DEM预处理:
  原始DEM数据可能存在一些问题，如无数据值 (NoData)、异常高程点 (spikes/pits)、条带噪声等。在使用DEM进行分析前，通常需要进行预处理，例如：
  • 填洼 (Sink Filling / Depression Filling): 填平DEM中的洼地，确保水流路径的连续性，这对于水文分析尤为重要。
  • 重采样 (Resampling): 当需要统一多个不同分辨率DEM时，或为特定分析调整分辨率时使用。
  • 拼接与裁剪 (Mosaic & Clip): 将多个DEM瓦片拼接成一个完整区域，或根据研究区边界裁剪DEM。

B. QGIS核心栅格工具进行地形分析

QGIS自身提供了一系列核心的栅格分析工具，可以直接对加载的DEM数据进行操作，提取基本的地形参数。这些工具通常位于主菜单的“栅格” (Raster) -> “分析” (Analysis) 目录下，或者通过“处理工具箱”中的GDAL工具调用。

• 坡度 (Slope): 计算每个像元地表倾斜的陡峭程度，通常以度或百分比表示 (44)。
  • 操作: “栅格” -> “分析” -> “坡度”。输入DEM图层，选择输出格式和路径，设置坡度单位。
• 坡向 (Aspect): 计算每个像元地表坡面的朝向，通常以0-360度表示（0/360度代表正北） (44)。
  • 操作: “栅格” -> “分析” -> “坡向”。输入DEM图层，选择输出格式和路径。
• 山体阴影/晕渲 (Hillshade): 模拟太阳光照在地表产生的阴影效果，增强地形的立体感，使地图更易于判读 (44)。
  • 操作: “栅格” -> “分析” -> “山体阴影”。输入DEM图层，设置太阳方位角 (Azimuth)、太阳高度角 (Altitude)、垂直夸大系数 (Z-factor)等参数。
• 等高线 (Contour): 从DEM生成具有相同高程值的连续曲线 (45)。
  • 操作: “栅格” -> “提取” -> “等高线”。输入DEM图层，设置等高线间距 (Interval)、输出矢量文件名等。
• 地貌崎岖度指数 (Ruggedness Index): 定量描述地形的异质性或起伏变化程度，例如Riley等人提出的方法，通过计算每个像元及其周围3x3邻域内高程变化的总和来度量 (44)。
  • 操作: “栅格” -> “分析” -> “崎岖度指数”。
• 地貌晕渲图 (Relief): 结合高程和阴影信息，创建具有彩色渲染效果的地形图，不同颜色代表不同高程范围，同时叠加阴影以增强立体感 (44)。
  • 操作: “栅格” -> “分析” -> “地貌”。

这些QGIS内置工具操作简便，能够满足许多基本的地形可视化和初步分析需求，是进行更复杂地形建模和水文分析的基础。

C. 通过QGIS Processing工具箱使用SAGA GIS进行高级地形分析

SAGA GIS (System for Automated Geoscientific Analyses) 是一款功能强大的开源GIS软件，尤其在地形分析、水文分析和地貌参数提取方面拥有众多专业算法。QGIS通过“处理工具箱” (Processing Toolbox) 完美集成了SAGA GIS的工具，用户无需离开QGIS环境即可调用这些高级功能 (46)。

DEM预处理 (关键步骤):

• 洼地填充 (Fill Sinks): 这是许多水文分析和地形参数提取前的必要步骤，用于消除DEM中可能存在的伪洼地，保证水流路径的连续性。SAGA提供了多种填洼算法，其中：
  • Fill Sinks (Wang & Liu): 一种常用且效果较好的填洼算法 (48)。
  • Fill Sinks (Planchon & Darboux): 另一种可选算法。
  • 操作: 在处理工具箱中搜索“Fill Sinks”，选择相应算法，输入原始DEM，输出填洼后的DEM。

高级地形参数提取与分析 (部分示例):
SAGA GIS提供了极为丰富的地形分析工具，以下列举一些常用的：

• 地形湿度指数 (Topographic Wetness Index, TWI): 也称地形指数或汇流指数，反映了地形对土壤水分空间分布的影响，值越高表示该区域越容易积水。
  • 工具路径 (示例): SAGA -> Terrain Analysis - Hydrology -> Topographic Wetness Index (TWI)
• 流域盆地 (Basins): 识别和划分地表径流汇集的区域。
  • 工具路径 (示例): SAGA -> Terrain Analysis - Channels -> Drainage Basins
• 河网与河网密度 (Channel Network & Density): 提取地表径流形成的河网，并可计算河网密度。
  • 工具路径 (示例): SAGA -> Terrain Analysis - Channels -> Channel Network
• 地形位置指数 (Topographic Position Index, TPI): 通过比较中心像元高程与其邻域像元的平均高程，来识别山顶、山谷、平坡、陡坡等地形位置 (50)。
  • 工具路径 (示例): SAGA -> Terrain Analysis - Morphometry -> Topographic Position Index (TPI)
• 曲率 (Curvature): 包括平面曲率 (Plan Curvature)、剖面曲率 (Profile Curvature)等，描述地表形态的弯曲程度。
  • 工具路径 (示例): SAGA -> Terrain Analysis - Morphometry -> Slope, Aspect, Curvature
• 相对高程/高程百分位 (Relative Heights and Slope Positions): 分析像元在其局部邻域内的高程位置。
• 山谷深度与山脊高度 (Valley Depth and Ridge Height)。
• 多尺度地形参数 (Multi-scale Terrain Parameters): 在不同分析窗口大小下计算地形参数，以揭示不同尺度的地形特征。

操作步骤 (通用流程):

1. 确保SAGA GIS已配置: QGIS安装时通常会包含SAGA GIS的接口。用户可以在“处理工具箱”的“选项”中检查SAGA GIS的配置情况，确保路径正确。较新版本的QGIS可能使用“SAGA Next Gen provider” (48)。
2. 打开处理工具箱: 通过菜单“处理” -> “工具箱”。
3. 查找SAGA算法: 在工具箱中，展开“SAGA”算法提供者，然后根据地形分析的类别（如Morphometry, Hydrology, Channels）查找所需工具。也可以直接在搜索框中输入工具名称关键词。
4. 选择并运行工具: 双击选中的SAGA工具，会弹出参数对话框。
   • 输入: 指定输入的DEM图层（通常是经过填洼处理的DEM）。
   • 参数: 根据具体算法的要求，设置其他必要的参数（如分析窗口大小、阈值、输出选项等）。
   • 输出: 指定输出栅格或矢量图层的保存路径和文件名，或选择保存为临时文件。
5. 点击“运行”。

重要性与考量:

• 功能扩展: SAGA GIS的集成极大地扩展了QGIS在地形分析方面的专业能力，使得用户可以免费获取并使用许多高级的地学分析算法。
• 算法版本与命名: SAGA GIS自身在不断发展，其工具的名称、参数甚至算法实现都可能随版本更新而变化。QGIS中集成的SAGA工具的名称和行为可能与原生SAGA GIS软件中的略有差异，或者随着QGIS和SAGA版本的更新而调整 (47 曾提及工具名称变化)。用户应关注其使用的QGIS和SAGA版本，并参考相应的文档。
• 工作流依赖: 许多高级地形分析（尤其是水文相关的）对输入DEM的质量非常敏感，正确的预处理（如填洼）是获得可靠结果的前提。

下表总结了一些在QGIS中常用的SAGA GIS地形分析工具，以供参考：

表3: QGIS中SAGA GIS关键地形分析工具 (示例)

通过结合QGIS核心工具和SAGA GIS的强大算法，用户可以对中国区域的自然地貌和地形特征进行全面而深入的分析。

VI. 流域绘制技术

流域绘制，即确定某一出水口点以上所有汇集地表径流的区域边界，是水文学、水资源管理、生态环境保护和灾害风险评估等领域的一项基础而重要的GIS分析任务。基于DEM的流域分析是实现这一目标的主要技术手段。

A. 基于DEM的流域分析基础

• 基本概念:
  • 流域 (Watershed/Catchment/Basin): 指地表径流通过一个共同出水口汇集而成的区域。流域边界通常是分水岭。
  • 汇水区 (Contributing Area/Upslope Area): 流向某一特定点或河段的所有上游区域。
  • 水系网络 (Stream Network/Channel Network): 地表径流汇集形成的河流、溪沟等地表水道系统。
  • 流向 (Flow Direction): DEM中每个像元的水流方向，通常指向其八个邻域中坡度最陡的那个像元。
  • 流量累积 (Flow Accumulation): 每个像元所承接的上游汇流像元的总数量，值越大的像元通常对应于河道位置。
  • 分水岭 (Divide/Water Divide): 分隔相邻流域的地形高地。
• 基本原理:
  基于DEM的流域分析通常遵循一个逻辑流程：首先，对原始DEM进行预处理（主要是填洼），以确保水流路径的正确性。然后，根据每个像元及其邻域的高程计算流向。接下来，基于流向计算每个像元的流量累积值。通过设定一个流量累积阈值，可以从流量累积栅格中提取出河网。最后，基于河网和流向信息，可以划分出各个子流域或整个研究区的流域盆地。如果指定了特定的出水口点，则可以精确 delineation 出该点的上游汇水区。

B. 在QGIS中使用SAGA GIS工具进行流域划分

SAGA GIS提供了一套完整且被广泛应用的工具链，用于基于DEM进行流域和水系网络的自动提取。这些工具同样可以通过QGIS的处理工具箱进行调用。

标准工作流程 (47)：

1. DEM预处理：填洼 (Fill Sinks)
   • 目的：消除DEM中的伪洼地，确保水流能够在地表连续流动。
   • SAGA工具：Fill Sinks (Wang & Liu) 或 Fill Sinks (Planchon & Darboux)。
   • 输入：原始DEM。
   • 输出：填洼后的DEM (Filled DEM)。
2. 流向计算 (Flow Direction)
   • 目的：确定每个像元的水流方向。
   • SAGA工具：通常，许多后续的水文工具（如流量累积）会内部计算流向，或者有专门的流向计算工具（如 Flow Direction D8）。
   • 输入：填洼后的DEM。
   • 输出：流向栅格 (Flow Direction Grid)。
3. 流量累积 (Flow Accumulation / Catchment Area)
   • 目的：计算每个像元汇集的上游像元数量，用于识别河道位置。
   • SAGA工具：Flow Accumulation (Top-Down) (在旧版或QGIS界面中可能仍称为 Catchment Area (47))。
   • 输入：填洼后的DEM (或流向栅格，取决于具体工具实现)。
   • 输出：流量累积栅格 (Flow Accumulation Grid)。
4. 河网提取 (Channel Network Extraction)
   • 目的：根据流量累积值识别和提取河网。
   • SAGA工具：Channel Network and Drainage Basins 或 Strahler Order (47)。
   • 关键参数：通常需要设置一个流量累积阈值 (Initiation Threshold)。只有当像元的流量累积值大于该阈值时，才被认为是河道的一部分。阈值的选择对河网的密度和形态有显著影响，需要根据研究区域的特点和分析目的进行调试。
   • 输入：流量累积栅格，填洼后的DEM。
   • 输出：河网矢量图层 (Channels) 和/或 河流分级栅格 (Strahler Order)，以及可能的初步流域盆地。
5. 流域盆地划分 (Watershed Basins Delineation)
   • 目的：基于河网和流向，将整个区域划分为若干个独立的流域盆地。
   • SAGA工具：Watershed Basins (通常与 Channel Network and Drainage Basins 工具一起输出，或作为独立工具)。
   • 输入：流向栅格，河网连接点或出水口。
   • 输出：流域盆地矢量或栅格图层。
6. （可选）根据特定出水口点提取子汇水区 (Delineate Sub-watershed for a Specific Outlet)
   • 目的：如果研究目标是确定某个特定点（如水文站点、水库坝址）的上游汇水区域。
   • SAGA工具：Upslope Area (49)。
   • 操作：
     1. 首先，需要在地图上精确确定出水口点的位置，并获取其X, Y坐标。可以创建一个点图层来标记出水口。
     2. 运行 Upslope Area 工具。
     3. 输入：填洼后的DEM (作为高程输入)，出水口的X坐标 (Target X Coordinate)，出水口的Y坐标 (Target Y Coordinate)。
     4. 选择计算方法 (Method)，如 Deterministic 8。
     5. 输出：代表该出水口上游汇水区的栅格图层。可以进一步将其矢量化为面要素。

重要考量:

• 工作流的严谨性: 流域分析是一个多步骤、有严格先后顺序的过程。每一步的输出质量直接影响后续步骤的准确性。
• 参数敏感性: 尤其是流量累积阈值的选择，对最终提取的河网密度和流域划分结果有决定性影响。用户需要结合研究区的实际地貌水文特征和先验知识，反复试验以确定最合适的参数值。
• DEM分辨率: DEM的分辨率会影响提取的河网和流域边界的精细程度。高分辨率DEM能提供更详细的结果，但也需要更多的计算资源。

C. 在QGIS中使用GRASS GIS工具进行流域划分

GRASS GIS (Geographic Resources Analysis Support System) 是另一款功能强大的开源GIS软件，同样在水文分析方面拥有成熟的模块和算法。QGIS通过处理工具箱也集成了GRASS GIS的工具 (47)。

GRASS GIS工作流程示例 (主要基于 r.watershed 和 r.water.outlet 工具，参考 47)：

1. 数据准备与导入 (GRASS环境):
   • GRASS GIS有其特定的数据管理结构 (LOCATION, MAPSET)。在QGIS处理工具箱中使用GRASS工具时，QGIS通常会在后台处理与GRASS环境的交互，但用户仍需理解其基本概念。
   • 输入DEM需要是GRASS可识别的栅格格式。如果DEM已在QGIS中加载，处理工具箱通常能直接使用。
2. 运行 r.watershed:
   • 目的：这是GRASS中进行流域分析的核心模块之一，可以同时计算多项水文参数，包括流量累积、流向、河网、流域盆地等。
   • GRASS工具：r.watershed。
   • 输入：
     • Elevation: 输入的DEM栅格图层。
     • Threshold for stream extraction: 定义河网的最小流量累积值（汇流单元数）。这是提取河网的关键参数。
   • 输出 (部分可选)：
     • Accumulation: 流量累积栅格图。
     • Drainage direction: 流向栅格图。
     • Stream segments: 提取的河网栅格/矢量图。
     • Half basins: 初步划分的子流域。
     • 更多其他输出...
   • 操作：在QGIS处理工具箱中找到 r.watershed，配置输入DEM和阈值等参数，指定需要的输出图层。
3. （可选）运行 r.water.outlet:
   • 目的：如果需要根据一个或多个指定的出水口点来精确划分其上游汇水区。
   • GRASS工具：r.water.outlet。
   • 输入：
     • Drainage direction: 由 r.watershed 或其他GRASS工具生成的流向栅格图。
     • Coordinates of outlet point(s): 出水口的东坐标 (Easting) 和北坐标 (Northing)。可以手动输入，或从包含出水口坐标的文本文件导入。
   • 输出：
     • Basin: 代表指定出水口上游汇水区的栅格图层。

重要考量:

• GRASS的工作范式: GRASS GIS的算法和数据处理方式可能与SAGA GIS有所不同，即使是功能相似的工具，其内部实现和参数设置也可能有差异。这可能导致最终结果在细节上有所不同 (52 曾对比不同工具的结果)。
• 工具选择: 用户可以根据个人对GRASS或SAGA的熟悉程度、特定算法的性能表现或项目需求，选择使用任一工具集，或者在某些情况下结合使用。
• 学习曲线: GRASS GIS本身是一个庞大而复杂的系统，尽管QGIS的集成简化了其使用，但要充分发挥其功能，仍需要一定的学习。

D. (可选) 使用 Watershed Generator 3.0 插件

对于希望简化流域划分流程的用户，一些第三方插件提供了更自动化的解决方案。Watershed Generator 3.0 (53) 就是这样一款插件，它旨在通过最少的步骤自动完成流域和河网的划分。

• 核心功能与操作:
  1. 输入: 用户只需提供一个DEM栅格图层。
  2. 定义范围 (Extent): 指定分析的地理范围。
  3. 指定出水口 (Outlet/Pour Point): 在地图上点选或输入出水口的坐标。
  4. 参数调整 (可选): 用户可以调整最小坡度 (minimum slope) 和阈值 (Threshold) 等参数来优化结果。
• 依赖:
  • 该插件的3.0版本明确依赖于QGIS中安装并配置好的 SAGA Next Generation 和 GRASS GIS (53)。这意味着SAGA和GRASS的相应处理模块需要在QGIS中可用。
  • 如果用户使用的是旧版的SAGA Provider (而非SAGA Next Generation)，则需要使用该插件的2.0版本 (53)。
• 安装与使用:
  • 通过QGIS插件管理器搜索并安装 Watershed Generator。
  • 确保满足其对SAGA和GRASS版本的依赖。
  • 启动插件后，按照其界面提示上传DEM、定义范围和出水口即可。

重要考量:

• 简化与黑箱: 此类自动化插件通过封装底层的SAGA或GRASS算法（或自定义算法）来简化用户操作。这对于快速获取结果或不熟悉复杂参数设置的用户非常友好。但同时也可能意味着用户对分析过程的控制力减弱，某些中间步骤和参数被“黑箱化”。
• 依赖管理: 插件对特定版本的外部工具（如SAGA Next Generation）的依赖，可能给用户在QGIS环境配置上带来一些挑战。如果依赖未能正确满足，插件可能无法正常工作。

总结而言，QGIS结合SAGA GIS和GRASS GIS的处理能力，为流域绘制提供了强大而灵活的开源解决方案。用户可以根据自己的技术水平、项目需求和对分析过程控制的期望，选择使用底层的SAGA/GRASS工具链，或尝试使用更自动化的第三方插件。无论选择哪种方法，高质量的DEM数据和对水文分析基本原理的理解都是成功的关键。

VII. 高级文字标注与地图整饰

清晰、美观且信息准确的文字标注是高质量地图不可或缺的组成部分。QGIS提供了强大的核心标注引擎，并辅以插件，能够满足包括复杂中文标注在内的各种高级需求。

A. QGIS核心标注引擎概述

QGIS的核心标注引擎功能非常全面，允许用户对地图要素进行精细化的文字标注控制 (55)。主要功能包括：

• 基于属性的标注内容: 可以直接使用图层属性字段的值作为标注内容，或通过表达式编辑器组合多个字段、添加文本、进行计算和格式化，生成动态的标注文本。
• 丰富的文本格式化选项:
  • 字体 (Font): 选择系统已安装的各种字体。
  • 样式 (Style): 如粗体、斜体、下划线、删除线等。
  • 大小 (Size): 支持多种单位（毫米、磅、像素等）。
  • 颜色 (Color) 与 透明度 (Opacity)。
  • HTML格式支持: 允许使用一部分HTML标签 (如 \<b>, \<i>, \<sup>, \<sub>, \<img>, \<div>, \<p>) 和CSS样式 (如 font-family, color, text-align) 来实现更复杂的文本排版和图文混排 (55)。
• 多样的标注放置规则: 针对不同几何类型（点、线、面）的要素，QGIS提供了多种智能的标注放置策略：
  • 点要素: 可选择笛卡尔式放置（如右上、左下等八方向优先）、环绕点放置、偏移放置等 (55)。
  • 线要素: 可选择平行于线、弯曲沿线、水平放置等，并可设置标注与线的距离、重复标注间距、是否允许越过线段端点等 (55)。
  • 面要素: 可选择在面内部水平放置、自由角度放置、沿周长放置、偏移自质心等多种方式 (55)。
• 碰撞检测与处理: 自动检测标注之间的重叠，并根据设定的规则（如优先级、是否允许叠置）来避免或解决冲突，确保标注的清晰可读 (55)。
• 叠置规则与优先级: 可以为不同图层或图层内的不同要素类别设置标注的优先级，高优先级的标注会优先显示。
• 数据定义的标注属性: 几乎所有的标注属性（如内容、字体、大小、颜色、位置、旋转角度、可见性等）都可以通过要素的属性值或表达式进行数据驱动的动态设置，实现“一标一式”的精细控制 (55)。
• 比例尺相关的可见性: 可以设置标注在特定地图比例尺范围内才显示，以适应不同缩放级别下的地图信息密度。
• 引线 (Callouts): 当标注无法直接放置在要素旁时，可以使用引线将其连接到对应要素。
• 背景、缓冲、阴影: 可以为标注添加背景框、文字缓冲（光晕效果）和阴影，以增强标注在复杂背景下的可读性 (55)。

QGIS核心标注引擎的这些功能，使得用户无需依赖第三方插件即可完成绝大多数常见乃至复杂的地图标注任务。其灵活性和可配置性足以媲美许多商业GIS软件。

B. 中文文字标注最佳实践

在QGIS中进行中文文字标注时，除了利用核心引擎的通用功能外，还需要特别注意以下几点，以确保标注效果既符合中文阅读习惯，又美观清晰：

• 字体选择:
  • 优先选择高质量中文字体: 确保所选字体包含完整的GBK或GB18030字符集，以避免出现缺字（显示为方框“□”）的情况。
  • 推荐字体:
  • 开源字体: 思源黑体 (Source Han Sans / Noto Sans CJK)、思源宋体 (Source Han Serif / Noto Serif CJK) 系列。这些字体由Adobe和Google合作开发，字形优美，字符集完整，且免费可商用，跨平台兼容性好。
  • 系统自带字体: Windows下的微软雅黑、宋体、黑体；macOS下的苹方、黑体-简/繁等。这些字体在对应操作系统上通常表现良好。
  • 版权考量: 在商业项目或公开发布的地图中，务必使用有合法授权的字体。
  • 在QGIS中设置: 在标注属性的“文本”选项卡中选择字体。
• 可读性与易认性:
  • 字号 (Size): 根据地图的最终用途（打印、屏幕显示）和目标比例尺选择合适的字号。中文笔画相对复杂，过小的字号会影响辨识度。通常，地图上的中文标注不宜小于6磅 (pt) (56 提及最小6pt)。
  • 颜色与对比度: 标注颜色应与背景图层形成足够的对比度，以保证清晰可读。WCAG (Web Content Accessibility Guidelines) 建议文本与背景的对比度至少为4.5:1 (56)。可以使用QGIS的颜色选择器辅助判断。
  • 文字缓冲/光晕 (Buffer/Halo): 为中文字符添加一个细小的、颜色与背景对比强烈的缓冲（例如，白色或浅色背景上的深色文字加浅色缓冲），可以显著提高其在复杂背景（如遥感影像、地形图）上的可读性 (55)。
• 放置策略与中文排版习惯:
  • 通用原则: 遵循地图标注的六大基本原则：可读性、优先级、清晰性、关联性、美观性、一致性 (56)。
  • 水平优先: 中文地名通常以水平方向排列最为自然和易读。尽量优先选择水平放置。
  • 避免拥挤: 汉字字符本身所占空间较大，信息密度高。在标注密集区域，要特别注意避免文字相互堆叠或过于靠近，影响整体美感和信息获取。可以适当调整字间距、行间距，或使用更高级的避让策略。
  • 多行标注: 对于较长的中文地名或需要显示多行信息的标注，可以利用QGIS的自动换行功能（设置“换行宽度”）或在文本中使用换行符 (\n)。注意中文换行的自然断句习惯 (55)。
  • 对齐方式: 对于多行标注，选择合适的对齐方式（左对齐、居中对齐、右对齐）。
• 标注优先级与避让:
  • 设置优先级: 对于不同重要程度的地名（如省会城市 vs. 普通县城），应设置不同的标注优先级，确保重要地名优先显示，不易被覆盖 (55)。
  • 要素作为障碍物: 可以将某些图层要素（如行政区边界、河流）设置为标注的障碍物，避免标注文字压盖这些重要地理要素 (55)。
  • 利用QGIS的自动避让引擎: 充分利用QGIS的碰撞检测和候选位置算法，让软件自动优化标注布局。
• 数据定义的样式:
  • 对于需要根据要素属性（如行政级别、人口数量）来区分标注样式（如字体、大小、颜色）的情况，充分利用QGIS的“数据定义覆盖”功能。例如，可以编写表达式，使得省级地名用较大字号和特定颜色，县级地名用较小字号和另一种颜色。

通过综合运用QGIS核心标注引擎的强大功能，并结合中文排版的特点和最佳实践，用户完全可以制作出专业、美观且易读的中文地图。

C. 使用 EasyCustomLabeling 插件增强标注

尽管QGIS的核心标注引擎功能强大，但在某些需要对大量标注进行高度个性化、逐个精细调整的场景下，或者当数据驱动的标注设置逻辑非常复杂时，EasyCustomLabeling 插件 (57) 提供了一种辅助手段，旨在简化这类手动和数据定义标注的工作。

核心功能与工作流程:

1. 复制图层到内存: 该插件的核心操作是将用户选定的原始矢量图层复制为一个新的内存图层 (Memory Layer) (57)。
2. 添加专用标注字段: 在这个新生成的内存图层中，插件会自动添加一系列专门用于控制标注样式的属性字段。这些字段通常对应QGIS标注设置中的各种参数，例如：
   • 标注文本内容 (Label Text)
   • X坐标 (Label X Coordinate)
   • Y坐标 (Label Y Coordinate)
   • 旋转角度 (Rotation)
   • 字体名称 (Font Name)
   • 字体大小 (Font Size)
   • 字体颜色 (Font Color)
   • 是否显示 (Show Label - True/False)
   • 引线属性 (Callout properties)
   • 对齐方式 (Alignment) 等。
3. 激活并准备标注: 生成的内存图层会自动激活，并准备好使用QGIS的标注工具。用户现在可以直接编辑这个内存图层中新增的这些属性字段的值。
4. 数据定义覆盖: 接下来，用户需要进入这个内存图层的标注设置界面，将各个标注属性（如位置、旋转、颜色、字体等）的“数据定义覆盖”选项，分别链接到插件生成的对应属性字段上。
5. 精细控制: 完成上述链接后，用户就可以通过直接修改内存图层中每个要素的那些专用标注属性值，来实现对该要素标注的精确控制。例如，想把某个特定地名的标注向右移动一点，只需修改其“Label X Coordinate”字段的值即可。

重要考量:

• 数据持久化: EasyCustomLabeling 生成的标注控制数据默认保存在内存图层中。这意味着当关闭QGIS项目时，如果未进行特殊处理，这些自定义的标注设置可能会丢失。为了持久化保存这些设置，插件描述中提到需要配合使用 Memory Layer Saver 插件 (版本3.2或更高)。该插件可以将内存图层的内容保存到与QGIS项目文件 (.qgs) 关联的一个.mldata 文件中 (57)。
  • 注意: 插件描述曾提及 Memory Layer Saver 对QGIS 3的兼容性问题 (57)。虽然 EasyCustomLabeling 本身列出的兼容版本覆盖QGIS 3.x (57)，但用户在使用时仍需关注 Memory Layer Saver 的最新兼容情况，以确保工作成果能被正确保存和恢复。
• 工作流: 该插件引入了一种基于“辅助标注图层”的工作流。用户不是直接在原始数据图层上进行所有标注设置，而是通过操作一个专门用于标注的内存副本来实现精细控制。
• 适用场景: EasyCustomLabeling 特别适用于以下情况：
  • 地图上只有少数几个标注需要非常特殊的、手动调整的位置或样式，而其他大部分标注可以使用常规规则。
  • 需要基于非常复杂的逻辑来数据定义标注，直接编写QGIS表达式可能过于繁琐，而通过编辑属性表的方式更为直观。
• 中文支持: 插件本身不直接处理字符编码或字体渲染，它提供的是控制标注属性的字段。因此，只要在QGIS的标注设置中为该内存图层选择了合适的中文字体，并正确设置了文本编码，那么通过该插件控制的中文标注就能正常显示。

总而言之，EasyCustomLabeling 可以作为QGIS核心标注引擎的一个有益补充，尤其是在需要对标注进行大量精细手动调整或实现复杂数据定义样式时，它提供了一种更为便捷的属性驱动控制方式。用户在使用时应注意其内存图层的工作机制和数据持久化依赖。

VIII. 中国关键地理空间数据源摘要

获取准确、适用的地理空间数据是进行任何GIS分析和制图项目的前提。本节将概述一些针对中国区域的关键数据源类型及其获取途径。

A. 行政区划数据 (Administrative Boundaries)

行政区划矢量数据（通常包括省、市、县、乡等级别）是进行社会经济分析、人口统计、区域规划、资源管理等众多应用的基础地理框架。

• CnOpenData (中国开放数据) (60):
  • 数据特点: 提供从2013年至2023年初的中国行政区划矢量数据，格式为Shapefile。该数据集据称以2022年国家基础地理信息数据中的县级区划为矢量基础，并参考了高德地图和天地图的行政区划数据，属性信息则参考了中国民政部历年发布的行政区划代码和名录。其宣称的特点包括时间序列长、属性信息丰富、国界线准确，并且会进行持续更新。
  • 行政级别: 包含省、市、县等级别的矢量数据。截至2023年1月1日，数据涵盖34个省级单位、333个地级单位和2843个县级单位 (60)。
  • 重要说明:
  • 数据更新至2023年2月12日 (60)。
  • 数据不包含行政区划的“托管”或“代管”情况（例如，云南省昆明市托管西双版纳傣族自治州勐腊县磨憨镇），仅反映在民政部正式注册的行政区划调整 (60)。这对于需要精确法定管辖范围的研究非常重要，因为托管关系与正式行政隶属关系在法律和管理上存在差异。
  • 属性字段丰富，例如县级矢量包含地名、拼音、区划代码、所属各级行政单元名称及代码、曾用名、备注、年份等 (60)。
  • 获取方式: 可通过其网站或相关公众号获取。对于需要高质量、详细属性且关注官方区划变更的中国区域研究，CnOpenData是一个值得重点考虑的数据源。
• Simplemaps.com (63):
  • 提供经过简化的中国GIS边界数据，格式为GeoJSON，采用WGS84投影。包含一级行政区（省、自治区、直辖市）的名称和ID。数据简化是为了快速加载，可能损失部分细节。
• OpenDataSoft (world-administrative-boundaries) (64):
  • 这是一个全球行政区划数据集，用户可以筛选出中国的数据，并下载为Shapefile、GeoJSON、KML等多种格式。需要注意的是，免费下载Shapefile格式可能有记录数量的限制（如50,000条）。
• Esri China (Hong Kong) Ltd Open Data (65):
  • 该平台主要提供香港特别行政区的各类开放地理空间数据，如交通网络、环境设施、土地利用等。虽然不直接提供中国大陆的行政区划，但可作为区域性高精度开放数据源的一个参考。

B. 数字高程模型 (Digital Elevation Models - DEMs)

DEM数据是地形分析、水文模拟、三维可视化等应用的核心基础数据。相关数据源已在 V.A 节 中有所介绍，此处简要回顾：

• OpenTopography (37): 这是一个重要的全球及区域高程数据分发门户，用户可以从中获取多种DEM产品，如：
  • SRTM (约30米/90米分辨率)
  • ALOS AW3D30 (约30米分辨率)
  • Copernicus DEM (GLO-30, GLO-90，即30米和90米分辨率)
  • 以及其他区域性高分辨率激光雷达 (LiDAR) 数据（若有覆盖）。 该平台通常提供便捷的按区域选择、数据格式转换和下载服务。
• NextGIS (41):
  • 提供覆盖中国的30米分辨率DEM数据产品，以及基于此生成的等高线和山体阴影图。这些数据通常是商业提供的，但质量和一致性可能较高。

C. 道路网络及其他专题数据 (Road Networks and Other Thematic Data)

• 道路网络数据:
  • MapOG (mapog.com) (66): 该网站声称提供中国的道路线状GIS数据，以及其他多种类型的地理空间数据图层。
  • University of Toronto Libraries - China Highways (67): 提供了一份1992年至2010年间的中国高速公路Shapefile数据。这类历史数据对于研究中国交通基础设施发展变迁非常有价值。
  • NextGIS Basemap Data (42): NextGIS提供的每日更新的矢量基础地图数据中，包含了公路 (Highways) 和铁路 (Railways) 网络图层，以及其他30多种数据类别。
  • OpenStreetMap (OSM): OSM是全球性的众包地图项目，其数据包含丰富的道路信息。可以通过QGIS的QuickOSM插件或直接从Geofabrik等OSM数据分发网站下载中国区域的OSM道路数据。OSM数据的详细程度和准确性在不同地区可能存在差异。
• 其他专题数据:
  • NextGIS Basemap Data (42): 除了交通网络，还提供建筑物（点和面）、水系（湖泊、河流、水道）、土地利用、兴趣点 (POI)、居民地（点和面）、植被、电力线等多种矢量数据图层。
  • Esri China (Hong Kong) Ltd Open Data (65): 如前所述，主要针对香港，但数据类型丰富，包括公共交通站点、电动汽车充电桩、渡轮码头、行人网络、敏感街道、隧道、政府维护的路面类型等。
  • 国家统计局及各部委网站: 中国国家统计局以及自然资源部、生态环境部、水利部等相关部委的官方网站，有时会发布一些与地理空间相关的统计数据或专题图件，可能需要进一步处理才能在GIS中使用。
  • 科研机构与大学: 一些科研机构和大学可能在其研究项目中产生并共享特定区域或主题的地理空间数据。

表4: 中国地理空间数据源概览 (示例)

在选择和使用这些数据源时，用户应仔细阅读其元数据，了解数据的来源、精度、更新时间、坐标系统、使用许可等关键信息，以确保数据能满足项目需求并合规使用。

IX. 总结与进阶学习

本指南详细介绍了QGIS的基础知识、针对中国地理空间任务的特定配置方法，并重点阐述了在行迹追踪、分布图制作、自然地貌与地形分析、流域绘制以及高级文字标注（特别是中文环境）等方面的核心技术和实用插件。通过结合QGIS的内置功能与强大的第三方插件（如Trajectools, Data Plotly, SAGA GIS, GRASS GIS, Chinese Coordinate Converter, EasyCustomLabeling等），用户可以高效地完成各类复杂的地理空间分析和制图任务。

A. 核心插件与技术回顾

• QGIS基础配置: 正确设置中国常用坐标参照系统 (CGCS2000, Beijing 1954, Xian 1980 及其高斯-克吕格投影带) 和添加本地化在线地图服务 (天地图、高德地图、腾讯地图等XYZ Tiles) 是在中国区域高效使用QGIS的前提。Chinese Coordinate Converter 和 GeoHey Toolbox 等插件为处理中国特有的“加偏”坐标系提供了关键支持。
• 行迹追踪: QGIS内置的“时态控制器”为时序数据可视化提供了便捷工具，而Trajectools插件则通过整合MovingPandas库，实现了高级轨迹分析功能，如轨迹创建、分割、平滑和属性计算。
• 分布图制作: QGIS强大的符号化引擎（分类、分级、热力图等）是制作高质量分布图的基础。Data Plotly插件进一步扩展了QGIS的数据可视化能力，允许用户创建交互式的统计图表并与地图联动。
• 地形分析: QGIS核心栅格工具可完成坡度、坡向、山影、等高线等基本地形参数提取。通过处理工具箱集成的SAGA GIS，用户可以进行更高级的地形分析，如填洼、地形湿度指数、地形位置指数等。
• 流域绘制: 结合SAGA GIS或GRASS GIS的工具，QGIS能够实现基于DEM的完整流域分析流程，包括填洼、流向与流量累积计算、河网提取和流域盆地划分。Watershed Generator 3.0等插件尝试简化此流程。
• 文字标注: QGIS核心标注引擎已非常成熟，支持丰富的文本格式化、智能放置、碰撞检测和数据定义样式。针对中文标注，选择合适的字体、遵循排版习惯至关重要。EasyCustomLabeling插件为需要精细手动调整的标注提供了便利。

B. QGIS持续学习资源

GIS技术和QGIS软件本身都在不断发展演进，保持学习的热情和能力是GIS从业者提升技能的关键。以下是一些推荐的QGIS持续学习资源：

• QGIS官方文档: 这是最权威、最全面的学习资料来源。
  • QGIS用户手册 (Desktop User Guide/Manual): 详细介绍了QGIS的各项功能和操作 (46)。
  • QGIS培训手册 (QGIS Training Manual): 提供了一系列循序渐进的练习教程 (46)。
  • PyQGIS开发者手册 (PyQGIS Cookbook): 面向希望通过Python进行QGIS脚本编写和插件开发的开发者 (68)。
  • QGIS服务器指南 (Server Guide/Manual): 针对QGIS Server的配置和使用 (68)。
  • 官方文档通常提供多种语言版本，包括英文和部分中文翻译 (68)。
• QGIS官方网站与社区:
  • QGIS官方博客 (QGIS.org Blog): 发布QGIS的最新动态、版本更新、社区活动等信息。
  • QGIS社区论坛/邮件列表: 用户可以在此提问、交流经验、寻求帮助。
• 在线教程与书籍:
  • QGIS Tutorials and Tips (qgistutorials.com): 提供大量针对不同QGIS版本和功能的免费图文教程 (27)。
  • 推荐书籍: 如 Kurt Menke 的《Discover QGIS 3.x》、Anita Graser 和 Gretchen Peterson 的《QGIS Map Design》(2)，以及 Scott Madry 的《Introduction to QGIS》(46) 等，这些书籍系统地介绍了QGIS的使用方法和制图技巧。
  • YouTube等视频平台: 有大量用户分享的QGIS操作视频教程。
• QGIS插件库 (QGIS Python Plugins Repository):
  • 官方插件库 (plugins.qgis.org) 是查找、安装和了解各种QGIS插件的主要入口 (69)。
  • 鼓励用户主动探索插件库，根据自己的需求发现更多实用的工具。插件库通常提供插件的描述、作者信息、版本历史、评分、下载量以及问题追踪系统的链接。

通过系统学习和不断实践，结合丰富的在线资源和活跃的社区支持，中国的QGIS用户定能充分发挥这款开源GIS软件的强大潜力，高效解决工作和研究中遇到的各类地理空间问题。

引用的著作

1. QGIS - Wikipedia, 访问时间为 五月 25, 2025， https://en.wikipedia.org/wiki/QGIS
2. GIS: Geographic Information Systems: QGIS - Cornell University Research Guides, 访问时间为 五月 25, 2025， https://guides.library.cornell.edu/gis/qgis
3. Road Map · QGIS Web Site, 访问时间为 五月 25, 2025， https://qgis.org/resources/roadmap/
4. Download · QGIS Web Site, 访问时间为 五月 25, 2025， https://qgis.org/download/
5. QGIS - Windows - University of Strathclyde, 访问时间为 五月 25, 2025， https://www.strath.ac.uk/professionalservices/is/software/qgis-windows/
6. qgis dOWNLOAD & INSTALLATION - Publish, 访问时间为 五月 25, 2025， https://publish.illinois.edu/illinoisdrainageguide/files/2022/07/QGIS-DOWNLOAD-INSTALLATION-GUIDE.pdf
7. Coordinate reference systems for "China" - EPSG.io, 访问时间为 五月 25, 2025， https://epsg.io/?q=China
8. China Geodetic Coordinate System 2000 - situx, 访问时间为 五月 25, 2025， https://situx.github.io/proj4rdf/data/def/crs/EPSG/0/4490/index.html
9. EPSG:4498 CGCS2000 / Gauss-Kruger zone 20 - Spatial Reference, 访问时间为 五月 25, 2025， https://spatialreference.org/ref/epsg/4498/
10. EPSG:4214 Beijing 1954 - Spatial Reference, 访问时间为 五月 25, 2025， https://spatialreference.org/ref/epsg/4214/
11. Projected coordinate systems for "China" - EPSG.io, 访问时间为 五月 25, 2025， https://epsg.io/?q=China%20%20kind%3APROJCRS
12. Beijing 1954 / 3-degree Gauss-Kruger CM 117E - EPSG:2436, 访问时间为 五月 25, 2025， https://epsg.io/2436
13. EPSG:4610 Xian 1980 - Spatial Reference, 访问时间为 五月 25, 2025， https://spatialreference.org/ref/epsg/4610/
14. EPSG:2332 Xian 1980 / Gauss-Kruger zone 18 - Spatial Reference, 访问时间为 五月 25, 2025， https://spatialreference.org/ref/epsg/2332/
15. Beijing 1954 / 3-degree Gauss-Kruger zone 29 - situx, 访问时间为 五月 25, 2025， https://situx.github.io/proj4rdf/data/def/crs/EPSG/0/2405/index.html
16. Beijing 1954 / 3-degree Gauss-Kruger zone 35 - situx, 访问时间为 五月 25, 2025， https://situx.github.io/proj4rdf/data/def/crs/EPSG/0/2411/index.html
17. Loading XYZ Tiles in QGIS, 访问时间为 五月 25, 2025， https://cawm.lib.uiowa.edu/xyz.html
18. Add Google Maps to QGIS 3 - SoCalGIS.org, 访问时间为 五月 25, 2025， https://socalgis.org/2019/11/06/add-google-maps-to-qgis-3/
19. How to Add Base Maps in QGIS: A Comprehensive Guide - GIS Courses Online, 访问时间为 五月 25, 2025， https://giscourse.online/base-maps-for-qgis/
20. Which WMS would work within China? - qgis - GIS StackExchange, 访问时间为 五月 25, 2025， https://gis.stackexchange.com/questions/313878/which-wms-would-work-within-china
21. qgis-basemaps/qgis_basemaps.py at main - GitHub, 访问时间为 五月 25, 2025， https://github.com/opengeos/qgis-basemaps/blob/main/qgis_basemaps.py
22. How to Add Base Maps to QGIS - DEV Community, 访问时间为 五月 25, 2025， https://dev.to/pokumars/how-to-add-base-maps-to-qgis-5f4
23. Tencent Maps | Bellingcat's Online Investigation Toolkit - GitBook, 访问时间为 五月 25, 2025， https://bellingcat.gitbook.io/toolkit/more/all-tools/tencent-maps
24. TianDiTu Tools — QGIS Python Plugins Repository, 访问时间为 五月 25, 2025， https://plugins.qgis.org/plugins/tianditu-tools/
25. Chinese Coordinate Converter — QGIS Python Plugins Repository, 访问时间为 五月 25, 2025， https://plugins.qgis.org/plugins/chinese_coordinate_converter/
26. GeoHey Toolbox — QGIS Python Plugins Repository, 访问时间为 五月 25, 2025， https://plugins.qgis.org/plugins/geohey_toolbox/
27. Animating Time Series Data (QGIS3) — QGIS Tutorials and Tips, 访问时间为 五月 25, 2025， https://www.qgistutorials.com/en/docs/3/animating_time_series.html
28. Must-Have QGIS Plugins To Elevate Your Mapping Projects - Geoawesome, 访问时间为 五月 25, 2025， https://geoawesome.com/must-have-qgis-plugins-to-elevate-your-mapping-projects/
29. Creating a Time-Series Animation with QGIS - GitHub, 访问时间为 五月 25, 2025， https://raw.githubusercontent.com/Brown-University-Library/geodata_pdf_tutorials/main/wildlife/bul_time_series_animation.pdf
30. All Plugins — QGIS Python Plugins Repository, 访问时间为 五月 25, 2025， https://plugins.qgis.org/plugins/?page=114\&sort=name\&per_page=449%2F\&dir=desc
31. trajectools - QGIS Plugins planet, 访问时间为 五月 25, 2025， https://plugins.qgis.org/planet/tag/trajectools/
32. movingpandas/qgis-trajectools: Trajectools - trajectory data ... - GitHub, 访问时间为 五月 25, 2025， https://github.com/movingpandas/qgis-processing-trajectory
33. 访问时间为 一月 1, 1970， https://plugins.qgis.org/plugins/trajectools/
34. Power of DataPlotly plugin: Creating graphs in QGIS software - Orbital Africa, 访问时间为 五月 25, 2025， https://orbital.co.ke/power-of-dataplotly-plugin-creating-graphs-in-qgis/
35. Data Plotly — DataPlotly 0 documentation, 访问时间为 五月 25, 2025， https://dataplotly-docs.readthedocs.io/en/latest/intro.html
36. Data Plotly — QGIS Python Plugins Repository, 访问时间为 五月 25, 2025， https://plugins.qgis.org/plugins/DataPlotly/
37. Find Topography Data - OpenTopography, 访问时间为 五月 25, 2025， https://opentopography.org/id/OTLAS.122009.32611.1\&minX=-118.28057687588213\&minY=37.06759770166735\&maxX=-118.2539788299186\&maxY=37.08364255871278
38. Data Catalog - OpenTopography, 访问时间为 五月 25, 2025， https://portal.opentopography.org/dataCatalog
39. How to download DEM in QGIS - SRTM Downloader Plugin # Lesson 8 (Part 2) of 29 # QGIS Tutorial. - YouTube, 访问时间为 五月 25, 2025， https://m.youtube.com/watch?v=XXbLtuOWM5s\&t=9s
40. Download SRTM DEM using QGIS | Open Topography | #part3 | GeoDev - YouTube, 访问时间为 五月 25, 2025， https://www.youtube.com/watch?v=M2aKbH2ysl8
41. China: elevation — download elevation (topography) map, 访问时间为 五月 25, 2025， https://data.nextgis.com/en/region/CN/dem/
42. China: map and data — download vector layers and ready-to-go GIS projects based on OSM, 访问时间为 五月 25, 2025， https://data.nextgis.com/en/region/CN/base/
43. Custom area: elevation — download elevation (topography) map - NextGIS Data, 访问时间为 五月 25, 2025， https://data.nextgis.com/en/region/custom/dem/
44. Raster Terrain Analysis Plugin, 访问时间为 五月 25, 2025， https://api.qgis.org/qgisdata/QGIS-Documentation-2.0/live/html/zh_TW/docs/user_manual/plugins/plugins_raster_terrain.html
45. Application of Qgis in the Generation of Contour Map - ResearchGate, 访问时间为 五月 25, 2025， https://www.researchgate.net/publication/376787182_Application_of_Qgis_in_the_Generation_of_Contour_Map
46. QGIS - GIS & Geospatial Technologies - LibGuides at University of California San Diego, 访问时间为 五月 25, 2025， https://ucsd.libguides.com/gis/tutorials-qgis
47. How to delineate watersheds in QGIS? - Geographic Information Systems Stack Exchange, 访问时间为 五月 25, 2025， https://gis.stackexchange.com/questions/257355/how-to-delineate-watersheds-in-qgis
48. Using SAGA tools in QGIS 3.30 or newer on Windows - YouTube, 访问时间为 五月 25, 2025， https://www.youtube.com/watch?v=dvcXeP0smTE
49. How to Delineate a Watershed in QGIS in 8 Steps - Equator Studios, 访问时间为 五月 25, 2025， https://equatorstudios.com/how-to-delineate-a-watershed-in-qgis/
50. Topographic Position Index in QGIS I topographic position index I TPI I QGIS I terrain analysis - YouTube, 访问时间为 五月 25, 2025， https://www.youtube.com/watch?v=t9Y-o868EWU
51. QGISHydro Webinar 4: Stream & Catchment Delineation - YouTube, 访问时间为 五月 25, 2025， https://www.youtube.com/watch?v=2Ub0c7Ss-T4
52. Stream and Catchment Delineation in QGIS 3 using GRASS tools - YouTube, 访问时间为 五月 25, 2025， https://www.youtube.com/watch?v=aTYHpft1gus\&pp=0gcJCdgAo7VqN5tD
53. Watershed Generator 3.0 - QGIS Model, 访问时间为 五月 25, 2025， https://hub.qgis.org/models/38/
54. Model tagged with: watershed-analysis - QGIS Hub, 访问时间为 五月 25, 2025， https://hub.qgis.org/models/tags/watershed-analysis/
55. 10.3. Setting a label — QGIS Documentation documentation, 访问时间为 五月 25, 2025， https://docs.qgis.org/latest/en/docs/user_manual/style_library/label_settings.html
56. 8 Unique Labeling Approaches for Geographic Features That Enhance Readability, 访问时间为 五月 25, 2025， https://www.maplibrary.org/1403/unique-labeling-approaches-for-geographic-features/
57. EasyCustomLabeling — QGIS Python Plugins Repository, 访问时间为 五月 25, 2025， https://plugins.qgis.org/plugins/EasyCustomLabeling/
58. All Plugins — QGIS Python Plugins Repository, 访问时间为 五月 25, 2025， https://plugins.qgis.org/plugins/?per_page=300
59. All Plugins — QGIS Python Plugins Repository, 访问时间为 五月 25, 2025， https://plugins.qgis.org/plugins/?page=11\&sort=latest_version_date\&per_page=47\&dir=asc
60. China administrative division shp data - CnOpenData, 访问时间为 五月 25, 2025， https://www.cnopendata.com/en/data/chinese-administrative-divisions-shp.html
61. China's administrative divisions data - CnOpenData, 访问时间为 五月 25, 2025， https://www.cnopendata.com/en/data/m/regional/chinese-administrative-divisions.html
62. CnOpenData China administrative division shp data, 访问时间为 五月 25, 2025， https://www.cnopendata.com/en/data/m/public/chinese-administrative-divisions-shp.html
63. Free GIS Maps of China - Resources - Simplemaps.com, 访问时间为 五月 25, 2025， https://simplemaps.com/gis/country/cn
64. World Administrative Boundaries - Countries and Territories — Opendatasoft, 访问时间为 五月 25, 2025， https://public.opendatasoft.com/explore/dataset/world-administrative-boundaries/export/?location=2,14.26438,51.32813\&basemap=jawg.light
65. Open Geospatial Data by Esri China (Hong Kong) Ltd., 访问时间为 五月 25, 2025， https://data-esrihk.opendata.arcgis.com/search?collection=Dataset\&layout=grid\&tags=%2520transport
66. China - GIS-DATA - MAPOG Map Analysis, 访问时间为 五月 25, 2025， https://story.mapog.com/app/gisdata/china/roads_line///
67. China Highways - Map and Data Library - University of Toronto, 访问时间为 五月 25, 2025， https://mdl.library.utoronto.ca/collections/geospatial-data/china-highways
68. QGIS resources · QGIS Web Site, 访问时间为 五月 25, 2025， https://qgis.org/resources/hub/
69. QGIS Plugins, 访问时间为 五月 25, 2025， https://plugins.qgis.org/
70. All Plugins — QGIS Python Plugins Repository, 访问时间为 五月 25, 2025， https://plugins.qgis.org/plugins/