1. 农业空间数据认知与QGIS环境搭建
农业空间数据处理的第一步,不是急着打开软件画图,而是先弄清楚"我们手里到底有哪些数据、它们长什么样、从哪来、能回答什么问题"。QGIS只是一个工具,真正决定分析质量的是对数据本身的认知。本节围绕四类最常见的农业空间数据展开,随后完成QGIS的安装、界面熟悉、插件配置与项目创建,为后续所有章节打下可复现的工作环境。
一、四类农业空间数据的定位与特征
农业空间数据可以按"表达对象"和"几何形态"两个维度来理解。地块数据是农业分析的骨架,通常以矢量面(Polygon)表达田块边界,属性表里记录地块编号、面积、承包户、作物类型、种植制度等字段;它的精度直接决定后续面积统计、变量施肥、产量核算的可信度。土壤数据多为采样点(Point)加插值生成的栅格面,或直接使用土壤类型图(面),关键属性包括pH、有机质、速效氮磷钾、质地等,采样点数据必须关注采样密度与空间自相关。气象数据以站点观测(点)或格网再分析产品(栅格/NetCDF)形式存在,涉及降水、积温、辐射等时间序列,使用时要注意时间分辨率与空间插值的匹配。遥感数据是栅格(Raster)的典型代表,如Sentinel-2的10米多光谱影像,用于提取NDVI、作物分类、长势监测,其核心参数是空间分辨率、重访周期、波段组合与云量。
这四类数据在几何形态、坐标系、时间维度上差异明显,混用时最容易出问题的正是坐标参照系(CRS)不一致。一个常见错误是把WGS84经纬度(EPSG:4326)的地块与投影坐标(如CGCS2000 3度带,EPSG:4547)的遥感影像直接叠加,导致位置偏移数百米。因此,在项目开始前统一CRS是必须遵守的约定:本文后续统一以CGCS2000 3度带高斯-克吕格投影作为面积计算与叠加分析的平面坐标系,经纬度数据仅用于展示与交换。
二、数据来源与格式的实用判断
实际工作中,数据来源往往比格式更值得关注。地块边界可能来自农村土地承包经营权确权登记的矢量成果(常见Shapefile或GeoPackage),也可能来自无人机航测正射影像的人工矢量化;土壤数据多来自测土配方施肥项目的采样点表(Excel坐标点)或区域土壤图;气象数据可从中国气象数据网下载站点日值,也可使用格网产品;遥感数据则通过Copernicus Open Access Hub或国内地理空间数据云获取。格式上,矢量优先使用GeoPackage(.gpkg),它单文件、支持多图层、字段名长度不受限,比Shapefile更稳健;栅格使用GeoTIFF;表格用CSV或Excel,但导入时必须明确X字段、Y字段与CRS。
三、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工程,工程文件只记录图层引用与样式,不打包数据,因此目录结构必须保持稳定。