1. 农业空间数据认知与QGIS环境搭建

农业空间数据处理的第一步,不是急着打开软件画图,而是先弄清楚"我们手里到底有哪些数据、它们长什么样、从哪来、能回答什么问题"。QGIS只是一个工具,真正决定分析质量的是对数据本身的认知。本节围绕四类最常见的农业空间数据展开,随后完成QGIS的安装、界面熟悉、插件配置与项目创建,为后续所有章节打下可复现的工作环境。

一、四类农业空间数据的定位与特征

农业空间数据可以按"表达对象"和"几何形态"两个维度来理解。地块数据是农业分析的骨架,通常以矢量面(Polygon)表达田块边界,属性表里记录地块编号、面积、承包户、作物类型、种植制度等字段;它的精度直接决定后续面积统计、变量施肥、产量核算的可信度。土壤数据多为采样点(Point)加插值生成的栅格面,或直接使用土壤类型图(面),关键属性包括pH、有机质、速效氮磷钾、质地等,采样点数据必须关注采样密度与空间自相关。气象数据以站点观测(点)或格网再分析产品(栅格/NetCDF)形式存在,涉及降水、积温、辐射等时间序列,使用时要注意时间分辨率与空间插值的匹配。遥感数据是栅格(Raster)的典型代表,如Sentinel-2的10米多光谱影像,用于提取NDVI、作物分类、长势监测,其核心参数是空间分辨率、重访周期、波段组合与云量。

这四类数据在几何形态、坐标系、时间维度上差异明显,混用时最容易出问题的正是坐标参照系(CRS)不一致。一个常见错误是把WGS84经纬度(EPSG:4326)的地块与投影坐标(如CGCS2000 3度带,EPSG:4547)的遥感影像直接叠加,导致位置偏移数百米。因此,在项目开始前统一CRS是必须遵守的约定:本文后续统一以CGCS2000 3度带高斯-克吕格投影作为面积计算与叠加分析的平面坐标系,经纬度数据仅用于展示与交换。

农业空间数据四类对象与QGIS处理链 地块数据 矢量面 Polygon 田块边界 / 作物类型 面积统计基础 土壤数据 采样点 Point / 栅格 pH / 有机质 / 养分 插值成面 气象数据 站点 / 格网 NetCDF 降水 / 积温 / 辐射 时间序列 遥感数据 栅格 Raster Sentinel-2 / NDVI 分辨率 / 云量 统一坐标参照系 CRS CGCS2000 3度带投影(面积计算基准) QGIS 处理链 导入 / 格式转换 重投影与裁剪 叠加与空间分析 制图与导出 成果:地块面积表 / 养分分级图 / 长势监测图 可复现的 .qgz 工程 + 规范命名的数据目录

二、数据来源与格式的实用判断

实际工作中,数据来源往往比格式更值得关注。地块边界可能来自农村土地承包经营权确权登记的矢量成果(常见Shapefile或GeoPackage),也可能来自无人机航测正射影像的人工矢量化;土壤数据多来自测土配方施肥项目的采样点表(Excel坐标点)或区域土壤图;气象数据可从中国气象数据网下载站点日值,也可使用格网产品;遥感数据则通过Copernicus Open Access Hub或国内地理空间数据云获取。格式上,矢量优先使用GeoPackage(.gpkg),它单文件、支持多图层、字段名长度不受限,比Shapefile更稳健;栅格使用GeoTIFF;表格用CSV或Excel,但导入时必须明确X字段、Y字段与CRS。

判断数据是否"可用"的三个快速检查:①打开图层属性看CRS是否已定义(显示"未知"就要警惕);②缩放到全图看范围是否落在合理经纬度或投影带内;③打开属性表看关键字段是否有空值与异常编码(如-9999)。

三、QGIS安装与界面熟悉

建议安装QGIS的LTR(长期支持)版本,稳定性和插件兼容性更好。Windows用户下载独立安装包,安装时勾选所需组件即可;安装完成后首次启动,界面由菜单栏、工具栏、图层面板(左侧)、地图画布(中央)、浏览器面板与处理工具箱(右侧)组成。初学者最该先熟悉三处:图层面板的右键菜单(设置CRS、缩放至图层、属性)、浏览器面板(直接拖拽本地数据入库)、以及"矢量→数据管理工具"和"栅格→投影"等常用菜单入口。

插件管理通过"插件→管理和安装插件"完成。农业空间分析中值得优先启用的包括:QuickMapServices(加载在线底图辅助定位)、QuickOSM(获取道路水系等辅助数据)、Profile Tool(地形剖面)、QGIS Resource Sharing(符号库共享)。插件不是越多越好,启用过多会拖慢启动速度,建议按项目需要逐个开启。

四、项目创建与目录规范

一个可复现的项目,从规范的目录结构开始。建议为每个项目建立如下结构,并在QGIS中通过"项目→属性→常规"设置项目主目录,使相对路径生效:

project_2024/
├── data_raw/        # 原始数据,只读不改
│   ├── parcel/      # 地块矢量
│   ├── soil/        # 土壤采样点与土壤图
│   ├── meteo/       # 气象站点与格网
│   └── rs/          # 遥感影像
├── data_processed/  # 重投影、裁剪后的中间数据
├── output/          # 分析成果与制图
└── project.qgz      # QGIS工程文件

创建项目时,先通过"项目→属性→CRS"将项目CRS设为CGCS2000 3度带对应EPSG代码(示例区域取EPSG:4547),再依次拖入地块、土壤、气象、遥感图层。每加入一个图层,立即检查其CRS并在必要时用"矢量→数据管理工具→重投影图层"或"栅格→投影→变形(重投影)"统一到项目CRS。最后保存为.qgz工程,工程文件只记录图层引用与样式,不打包数据,因此目录结构必须保持稳定。

常见坑:直接双击打开他人发来的.qgz工程时,若数据路径为绝对路径且原目录不存在,图层会显示红色感叹号。解决方法是右键图层→"修复数据源",重新指向本地数据;或在项目属性中勾选"保存路径为相对路径"。