湖北农业遥感估产信息系统

1、湖北省农业遥感估产信息系统详细设计说明书（报名编号K20090241）项目书拟制人_柳扌旨导老师梅新所在学校湖北大学二00九年十月十日详细设计说明书1. 引言1.1编写目的湖北省农业遥感估产系统是基于MapGIS7.0平台进行的二次开发，该软件主要是利用湖北省的卫星影像中农作物生长期的光谱特征信息，借助相关算法，完成对相关农作物的估产监测。本详细设计说明书就是对该项目进行详细设计，进一步明确系统结构，划分各个功能模块，实现农作物产量监测的自动化，效率化。1.2背景a. 待开发系统的名称：湖北省农业遥感估产系统b. 任务提出者：梅新老师开发者：柳菲，董康，黄奎，陈津晶，郑培蓓用户：农业监测相关部

2、门1.3理论基础W avel ength 6m )a:孕穗期 b :抽穗期 c :乳熟期 d :成熟期图水稻穗期不同生育期光谱特征由上图可分析出，在水稻四个生育时期，孕穗期、抽穗期、乳熟期、成熟期，光谱特征 区别最明显的为第一个反射峰附近，在第一个反射峰区，孕穗期与抽穗期光谱特征最为相似， 反射率也最低，乳熟期次之，成熟期最高；在反射谷区，反射率与反射峰区相似。分析原因，孕穗期，水稻剑叶全部露出叶鞘，抽穗期，20%穗子顶端从剑叶叶鞘中露出，这两个时期水稻茎叶发育最茂盛， 故对太阳辐射吸收最多，反射的较少，乳熟期，穗子顶部的子粒达到正常谷粒的大小，颖壳充满乳浆状内物，茎叶由盛转衰，对太阳辐射的吸

3、收开始减少，成熟期， 穗上有80以上的谷粒呈现该品种固有的颜色，叶片开发变黄，对太阳辐射的反射更多，反 射率更高。基于灰色系统理论的水稻产量监测，在水稻移栽时期，水稻种植区的水深一般在 2-15cm，影像中的水稻种植区是表层水体和秧苗的混合光谱特征，在移植的50-60天后，即孕穗期，水稻冠层几乎覆盖整个水稻种植区，影像中的水稻种植区是茂密的植被。选取水稻秧苗移栽时期，对土壤湿度和植被水分含量较敏感的短波红外波段计算得到的陆表水指数（LSWI）,以及水稻孕穗期对植被敏感的归一化植被指数 （NDVI）、增强植被指数（EVI）,从水 田中提取水稻种植面积。利用MODIS影像的时效性强，能很好的监测农

4、作物的长势， 得到农 作物各时期的NDVI和EVI影像，一般选用每年6月到9月的MODIS影像；利用TM影像的精 度高，进行监督分类，配合耕地分布利用图，能很好的提取估产区域农作物的种植面积。2. 系统的结构湖北省农业遥感估产信息系统关于帮助湖北省农业遥感估产系统的结构序号模块名称子模块1子模块21文件打开打开地图打开影像关闭保存另存为版面视图打印输出打印设置打印退出2基本操作Shp文件上载放大缩小移动复位获取AOI动态投影重采用影像格式转换矢量转栅格栅格转矢量3影像显示原图显示反转显示均衡化显示透明显示卷帘显示闪烁显示鹰眼显示漫游显示4影像处理投影变换几何校正太阳高度角校正影像裁剪薄云去除噪

5、声去除5影像分类监督分类监督最小距离分类:监督广义最小距离分类监督平行六面体分类最大似然法BP神经网络分类RBF神经网络分类LVQ2神经网络分类高阶神经网络分类非监督分类ISODATA法分类非监督最小距离分类r非监督广义最小距离分类非监督平行六面体分类ART神经网络分类FuzzyART神经网络分类Konhonen非监督分类分类后处理精度评价r聚类评价去除分析过滤分析分类后比较法作物面积提取6影像光谱信息影像合成波段运算快速计算7遥感估产单产模型一兀线性回归模型拉格朗日插值模型牛顿插值模型曲线拟合模型总产模型8关于帮助关于操作手册3.1所需数据描述影像数据：湖北省 2002年161天到241天M

6、ODIS影像和TM影像矢量数据：湖北省行政区边界，估产采样区域，湖北省行政区耕地利用分布图文本数据：估产样本区域的时间信息和产量信息Excel数据3.2数据采集影像数据可以在网上下载或者购买获取矢量数据需要利用 GPS等仪器实地测量获取文本数据咨询当地相关农业部门，例如统计年鉴等官方文本3.2数据处理1. 影像数据需要做影像校正使相关光谱信息能很好的显示，例如：太阳高度角校正，几何校正，薄云去除，噪声去除等，然后对 MODIS影像进行NDVI, EVI或者LSWI处理，对TM影像 进行分类提取估产区域的处理；2. 对影像数据和矢量数据做投影变换，使之在同一投影系下重叠显示，用湖北省行政区边界和

7、TM影像提取估产区域的影像图，用估产采样区域和 TM影像提取样本区域影像，配合湖北 省耕地利用分布图进行人工修正，获取实际的估产区域；3. 比对MODIS影像和TM影像进行像元值赋值，参照文本数据，利用样本区域像元值，选择模型建模，计算整个估产区域各个像元值对应的产量，最后累加获得总产量。注:本项目采用的估产方法特点， 利用MODIS影像的时效性强，能很好的监测农作物的长势； 利用TM影像的精度高，能很好的提取估产区域的面积。先求出单产，再利用单产*面积累加求得总产，方法主要是利用样本区域的像元值和对应的单产值，选择单产模型来计算其他区域像元值对应的单产值，最后农作物总产量=单产1*种植面积*

8、修正系数+单产n*种植面积*修正系数（不到一个像元的修正系数取 0.5，以后可以修正），其中单产模型有一元线性回 归模型，拉格朗日插值模型，牛顿插值模型，曲线拟合模型。具体方法可参见以下所示的估产流程图：4 .模块设计说明4.1文件模块文件打开关闭保存另存为版面视图打印输岀退岀打开地图打开影像打印设置打印文件模块包括了打开，关闭，保存，另存为，打印输出，退出等功能子模块，该模块实 现了对文件的基本管理操作。4.2基本操作模块基本操作AOI编辑影像融合重采样影像格式转换矢量转栅格栅格转矢量基本操作模块功能对影像的基本管理，可以进行放大，缩小，移动，复位等操作，影像格式转换使影像能转换成本系统可处

9、理的格式，获取AOI功能方面后面对于影像裁剪，监督分类的操作，重采样处理使影像特征显示的更为明显。如图4-2所示。4.3影像显示模块影像显示原图显示卷帘显示反转显示透明显示均衡化显示闪烁显示鹰眼显示漫游显示原图显示：原始图像的显示反转显示：为了突出感兴趣目标所在的灰度区间，相对抑制不感兴趣的灰度区间 均衡化显示：使各灰度级上的像元分布均匀，增强影像上地物间的反差 卷帘显示：使叠加的影像中一幅成可设置值的透明度显示透明显示：将一幅影像设置透明度，与第二幅影像进行对比显示 闪烁显示：将一幅影像进行闪烁显示，与第二幅影像进行对比显示 鹰眼显示：对影像的局部区域进行放大显示漫游显示：将一幅影像以一定的

10、视角和形式进行显示4.4影像处理模块影像处理投影变换几何校正太阳高度角校正影像裁剪薄云去除噪声去除投影变换：对影像进行转投影，使两幅影像能在同一投影系下显示 几何校正：通过计算机对离散结构的数字图像中的每个像元进行校正 太阳高度角校正：消除太阳高度角导致辐射亮度失真的过程 影像裁剪：获取影像中需要的处理的部分 薄云去除：消除影像拍摄时云层带来的影像光谱信息失真过程噪声去除：利用傅立叶变换法、小波变换法和插值法等去除影像上的噪声影响4.5影像分类模块影像分类监督分类分类后处理分监督分类作物面积提取H监督最小距离分类精度评价过滤分析分类后比较法.SODATA分 类最大似然法聚类分析去除分析非监督最

11、小距离法 |该模块功能通过对影像进行特定方式分类，配合耕地利用分布图来提取估产面积区域, 然后通过分类后的分析，评价判断该分类方法的拟合度，并利用重编码功能提取估产信息。4.6影像运算模块影像运算影像合成'波段运算快速计算本模块根据农作物在影像上的特征，对影像进行处理，使该农作物的光谱信息更为明显，便于对相关区域进行信息提取和估产操作。NDVI:检测植被生长状态、植被覆盖度和消除部分辐射误差等。归一化植被指数=（近红外波段-红光波段）/ （近红外波段+红光波段） EVI:提高对高生物量区的敏感度，通过削弱叶冠背景信号和降低大气的影响来改 善对植被的监测。增强植被指数=G*（近红外波段-

12、红光波段）/ （近红外波段+ Ci*红光波段-G*蓝光 波段+L）其中G, C，C, L为相关系数LSWI:利用了对水体敏感的短波红外波段，对土壤湿度和植被水分含量较敏感。 陆表水指数=（近红外波段-短波红外波段）/ （近红外波段+短波红外波段）4.7遥感估产模块遥感估产总产模型单产模型一元线性回归模型拉格朗日插值模型牛顿插值模型曲线拟合模型Ir一元线性回归模型：yt = b0 + bl xt + ut (4.7.1)上式表示变量yt和xt之间的真实关系。其中yt称作被解释变量(或相依变量、因变 量)，xt称作解释变量(或独立变量、自变量)，ut称作随机误差项，b0称作常数项(截距项)，bl称

13、作回归系数。在模型(4.7.1) 中，xt是影响yt变化的重要解释变量。b0和bl也称作回归参数。 这两个量通常是未知的，需要估计。t表示序数。当t表示时间序数时，xt和yt称为时间序列数据。当t表示非时间序数时，xt和yt称为截面数据。ut则包括了除xt以外的影响yt 变化的众多微小因素。ut的变化是不可控的。上述模型可以分为两部分。(1)b0 +b1 xt是非随机部分；(2)ut是随机部分。拉格朗日插值模型：n -1n个节点的拉格朗日插值多项式公式为P(x"w(x)x_xy)w(Xi)n -1n-1其中 yi = f Xi ， w(x) = ： (x - Xi)， w XiX -

14、 Xji =0jT且 j：#牛顿插值模型：牛顿插值是基于下面这些的公式:fxO,x1,.xk=(fx1,.xk-fxO,.xk-1)/(xk-xO)fx=f(x)1)+R n(x)前两个是均差的递推关系式，而后一个就是牛顿插值公式，其中f(x)=fx0+fx0,x1(x-x0)+fx0,x1,x2(x-x0)(x-x1)+.fx0,.x n(x-x0).(x-x n-N(x)=f(x)-R n(x)，即目标 多项式，Rn(x)是n阶 插值余项，我们就是用 N(x)去近似f(x)。曲线拟合模 型：用连续曲线近似 地刻画或比拟平面 上离散点组所表示 的坐标之间的函数关系的一种数据处理方法。通过实验或观测得到量x与y的一组数