基于Sentinel-2的杞县大蒜提取试验
一、前言
杞县大蒜已有2000多年的栽培历史是全国农产品地理标志之一,2009年9月14日,中华人民共和国农业部正式批准对“杞县大蒜”实施农产品地理标志登记保护。杞县已建成大蒜冷藏保鲜库600多座,年贮藏能力达到25万吨。建成蒜片、蒜米、蒜粒、蒜油等加工企业80多家,年加工能力5万多吨。已在全国20多个大中城市设立了大蒜销售网点,且远销到美国、欧洲、日本、韩国、东南亚、印尼等几十个国家和地区,每年直接间接出口大蒜制品25多万吨。
Sentinel-2(哨兵2)超大的幅宽,双星组合,较短的重访周期以及丰富的波谱信息都为农作物提取提供了宝贵的数据源。今天我们就来通过Sentinel-2来进行杞县大蒜提取试验。
二、数据选择
本次实验选取了三期的Sentinel-2 L1C数据,通过预处理(Sen2Cor大气校正)、SNAP采样输出到ENVI中进行 波段组合、裁剪等工作。
下面是本次实验的数据,分别为2019年3月14日、2018年4月8日、2019年4月28日,三期数据,下图是分别进行 B4B3B2 真彩色组合:
经过反复测试 B11B8B3 的假彩色组合显示最能突出地物间的差别,层次分明:
三、外业采样
本次选用手机APP:LocaSpace Viewer(LSV)软件来辅助外业采集样本,手机端打开GPS+WLAN+移动网络的定位模式,然后采集田块的角点作为坐标点。下面我们举例来说明:例如W1 W2 W3 W4为小麦样本的坐标点,G1 G2 G3 G4为大蒜样本的坐标点:
然后将样本点导入ENVI当中,蓝色点为大蒜地块的样本角点,红色点为小麦样本地块的样本角点:
采集的大蒜样本集:
采集的小麦样本集:
由上面的样本集图可见,大蒜的颜色在B11B8B3组合下呈现蓝紫色,小麦的颜色在B11B8B3组合下呈现亮绿色,下面我们分别来看两者的光谱曲线:
大蒜的光谱曲线:
小麦的光谱曲线:
对比两者的光谱曲线可以看到,在近红外波段小麦的反射率比大蒜要高,在绿波段和红波段大蒜的反射率比小麦要高,在短波红外波段,大蒜的反射率比小麦高。因此在B11 B8 B3的组合下,小麦要比大蒜更绿,大蒜要比小麦更蓝和更红。
四、2019年数据分类
我们先对2019年的数据进行分类,本次分类采用监督分类的方法,采用的是支持向量机的模型进行分类。因为数据时间是2019年3月14日,树木尚未返青,因此我们选取了大蒜、小麦、水域、居民地、空闲地五大类,并对样本进行了可分离度评价,结果如下:
由图可以看到各个样本之间的可分离度良好,基本都在1.9以上:
执行监督分类,下图依次是监督分类的结果