全 文 ：第 源愿 卷 第 圆 期
圆 园 员 圆 年 圆 月
林 业 科 学
杂悦陨耘晕栽陨粤 杂陨蕴灾粤耘 杂陨晕陨悦粤耘 灾燥造郾 源愿袁晕燥郾 圆云藻遭援 袁圆 园 员 圆
基于时间序列 酝韵阅陨杂 耘灾陨数据的森林生长异常监测鄢
刘丽娟员袁圆袁猿 摇 庞摇 勇员 摇 在澡葬灶早 载蚤葬燥赠葬灶早源 摇 杂增藻蚤灶 杂燥造遭藻则早缘摇
范文义猿 摇 李增元员 摇 李明泽猿
渊员郾 中国林业科学研究院资源信息研究所摇 北京 员园园园怨员曰 圆郾 杭州师范大学遥感与地球科学研究院摇 杭州 猿员园园猿远曰
猿郾 东北林业大学林学院摇 哈尔滨 员缘园园源园曰 源郾 晕韵粤粤 辕 阅耘杂阅陨杂 辕 杂栽粤砸袁缘圆园园 粤哉贼澡 砸阅袁 悦葬皂责 杂责则蚤灶早泽袁酝阅 圆园苑源远 哉灶蚤贼藻凿 杂贼葬贼藻泽曰
缘郾 晕燥则憎藻早蚤葬灶 云燥则藻泽贼 葬灶凿 蕴葬灶凿泽糟葬责藻 陨灶泽贼蚤贼怎贼藻袁 孕燥泽贼遭燥噪泽 员员缘袁 员源猿员魡杂袁 晕燥则憎葬赠冤
摘摇 要院摇 利用遥感技术袁基于时间序列的 酝韵阅陨杂数据对 圆园园源要圆园园远 年东北三省的林区进行森林生长异常监测遥
首先利用 酝韵阅陨杂数据时间分辨率高的特点袁采用 杂葬增蚤贼扎噪赠鄄郧燥造葬赠滤波函数平滑 愿 天合成的 耘灾陨袁计算生长季面积和
年 耘灾陨曲线熵值袁两指标联合得到 猿 年间变化量大的像素点袁定义为森林生长异常点曰然后抽取异常点的时间序列
曲线进行分析袁并结合森林灾害事件进行比较验证遥 结果表明院 异常点曲线的熵值明显大于正常年袁生长季峰值
低袁并且在生长旺季会出现峰值突然持续下降或双峰等异常现象袁这与该区域森林生长异常发生时的植被反射率
表征一致袁说明用该法对森林生长异常进行监测是基本可行的遥
关键词院摇 时间序列曰 酝韵阅陨杂 耘灾陨曰 生长季面积曰 熵曰 异常监测
中图分类号院 杂苑缘愿援 源摇 摇 摇 文献标识码院粤摇 摇 摇 文章编号院员园园员 原 苑源愿愿渊圆园员圆冤园圆 原 园园缘源 原 园怨
收稿日期院 圆园员园 原 园缘 原 圆源曰 修回日期院 圆园员园 原 园苑 原 圆远遥
基金项目院 国家自然科学基金渊源员园苑员圆苑圆冤曰林业公益性行业科研专项渊圆园园愿园源园园员冤曰栽澡藻 晕燥则憎藻早蚤葬灶 砸藻泽藻葬则糟澡 悦燥怎灶糟蚤造 砸耘酝云韵砸 孕则燥躁藻糟贼曰国防
科工局专项项目渊耘园猿园缘 辕 员员员圆 辕 园员 辕 园员冤遥
鄢庞勇为通讯作者遥
酝燥灶蚤贼燥则蚤灶早 云燥则藻泽贼 郧则燥憎贼澡 阅蚤泽贼怎则遭葬灶糟藻 哉泽蚤灶早 栽蚤皂藻 杂藻则蚤藻泽 酝韵阅陨杂 耘灾陨 阅葬贼葬
蕴蚤怎 蕴蚤躁怎葬灶员袁圆袁猿 摇 孕葬灶早 再燥灶早员 摇 在澡葬灶早 载蚤葬燥赠葬灶早源 摇 杂增藻蚤灶 杂燥造遭藻则早缘 摇 云葬灶 宰藻灶赠蚤猿 摇 蕴蚤 在藻灶早赠怎葬灶员 摇 蕴蚤 酝蚤灶早扎藻猿
渊员郾 陨灶泽贼蚤贼怎贼藻 燥枣 云燥则藻泽贼 砸藻泽燥怎则糟藻泽 陨灶枣燥则皂葬贼蚤燥灶 栽藻糟澡灶蚤择怎藻泽袁悦粤云摇 月藻蚤躁蚤灶早 员园园园怨员曰
圆郾 陨灶泽贼蚤贼怎贼藻 燥枣 砸藻皂燥贼藻 杂藻灶泽蚤灶早 葬灶凿 耘葬则贼澡 杂糟蚤藻灶糟藻泽袁 匀葬灶早扎澡燥怎 晕燥则皂葬造 哉灶蚤增藻则泽蚤贼赠摇 匀葬灶早扎澡燥怎 猿员园园猿远曰
猿郾 悦燥造造藻早藻 燥枣 云燥则藻泽贼则赠袁 晕燥则贼澡藻葬泽贼 云燥则藻泽贼则赠 哉灶蚤增藻则泽蚤贼赠摇 匀葬则遭蚤灶 员缘园园源园曰
源郾 晕韵粤粤 辕 阅耘杂阅陨杂 辕 杂栽粤砸袁缘圆园园 粤哉贼澡 砸阅袁 悦葬皂责 杂责则蚤灶早泽袁酝阅 圆园苑源远 哉灶蚤贼藻凿 杂贼葬贼藻泽曰
缘郾 晕燥则憎藻早蚤葬灶 云燥则藻泽贼 葬灶凿 蕴葬灶凿泽糟葬责藻 陨灶泽贼蚤贼怎贼藻袁 孕燥泽贼遭燥噪泽 员员缘袁 员源猿员魡杂袁 晕燥则憎葬赠冤
粤遭泽贼则葬糟贼院摇 云燥则藻泽贼 早则燥憎贼澡 蚤泽 皂葬蚤灶造赠 糟怎则则藻灶贼造赠 皂燥灶蚤贼燥则藻凿 怎泽蚤灶早 蚤灶鄄泽蚤贼怎 皂藻葬泽怎则藻皂藻灶贼泽 蚤灶 灶燥则贼澡藻葬泽贼 燥枣 悦澡蚤灶葬援 栽燥 藻枣枣藻糟贼蚤增藻造赠
皂燥灶蚤贼燥则 枣燥则藻泽贼 早则燥憎贼澡 凿蚤泽贼怎则遭葬灶糟藻 葬贼 造葬则早藻 泽糟葬造藻袁 憎藻 葬贼贼藻皂责贼藻凿 贼燥 怎泽藻 则藻皂燥贼藻 泽藻灶泽蚤灶早 贼藻糟澡灶蚤择怎藻袁 责葬则贼蚤糟怎造葬则造赠袁 贼蚤皂藻 泽藻则蚤藻泽
酝韵阅陨杂 凿葬贼葬 枣则燥皂 圆园园源 贼燥 圆园园远郾 栽澡藻 葬灶灶怎葬造 贼蚤皂藻 泽藻则蚤藻泽 燥枣 愿鄄凿葬赠 藻灶澡葬灶糟藻凿 增藻早藻贼葬贼蚤燥灶 蚤灶凿藻曾 渊耘灾陨冤 凿葬贼葬泽藻贼 憎葬泽 早藻灶藻则葬贼藻凿
葬灶凿 泽皂燥燥贼澡藻凿 怎泽蚤灶早 葬 杂葬增蚤贼扎噪赠鄄郧燥造葬赠 枣蚤造贼藻则援 栽澡藻 耘灾陨 贼则葬躁藻糟贼燥则赠 凿怎则蚤灶早 早则燥憎贼澡 泽藻葬泽燥灶 憎葬泽 泽蚤皂怎造葬贼藻凿 怎泽蚤灶早 葬 造燥早蚤泽贼蚤糟 皂燥凿藻造援
云则燥皂 贼澡藻 泽蚤皂怎造葬贼藻凿 贼则葬躁藻糟贼燥则赠袁 贼澡藻 耘灾陨 葬则藻葬 燥枣 早则燥憎贼澡 泽藻葬泽燥灶 葬灶凿 葬灶灶怎葬造 耘灾陨 藻灶贼则燥责赠 憎藻则藻 糟葬造糟怎造葬贼藻凿援 栽澡藻泽藻 贼憎燥 枣葬糟贼燥则泽
憎藻则藻 糟燥皂遭蚤灶藻凿 贼燥 皂葬责 贼澡藻 凿蚤泽贼怎则遭葬灶糟藻 则藻早蚤燥灶泽 燥枣 枣燥则藻泽贼 早则燥憎贼澡援 云蚤灶葬造造赠袁 贼澡藻 凿蚤泽贼怎则遭葬灶糟藻 则藻早蚤燥灶泽 憎藻则藻 增藻则蚤枣蚤藻凿 怎泽蚤灶早 葬 泽藻贼 燥枣
则葬灶凿燥皂 泽葬皂责造藻泽援 栽澡藻 则藻泽怎造贼 蚤灶凿蚤糟葬贼藻泽 贼澡葬贼 贼澡藻 凿蚤泽贼怎则遭葬灶糟藻 责燥蚤灶贼泽 澡葬增藻 凿蚤泽贼蚤灶糟贼蚤增藻造赠 澡蚤早澡藻则 藻灶贼则燥责赠 葬灶凿 造燥憎藻则 责藻葬噪援 杂燥皂藻 燥枣
贼澡藻泽藻 责燥蚤灶贼泽 葬造泽燥 泽澡燥憎 葬遭则怎责贼 耘灾陨 凿藻糟造蚤灶藻 凿怎则蚤灶早 贼澡藻 皂蚤凿泽藻葬泽燥灶 燥枣 贼澡藻 责藻葬噪 责澡葬泽藻泽 燥则 凿燥怎遭造藻 责藻葬噪泽援 栽澡蚤泽 葬责责则燥葬糟澡 蚤泽
凿藻皂燥灶泽贼则葬贼藻凿 贼燥 遭藻 枣藻葬泽蚤遭造藻 枣燥则 凿蚤泽贼怎则遭葬灶糟藻 皂燥灶蚤贼燥则蚤灶早 燥枣 枣燥则藻泽贼 早则燥憎贼澡援
运藻赠 憎燥则凿泽院摇 贼蚤皂藻 泽藻则蚤藻泽曰 酝韵阅陨杂 耘灾陨曰 早则燥憎贼澡 泽藻葬泽燥灶 葬则藻葬曰 藻灶贼则燥责赠曰 凿蚤泽贼怎则遭葬灶糟藻 皂燥灶蚤贼燥则蚤灶早
摇 摇 森林的健康状况直接关系到全球生态安全和人
类社会的可持续发展遥 随着自然资源过度开发袁森
林资源锐减袁森林生态系统的服务功能大大减退袁全
球出现了一系列的生态危机袁诸如生物多样性减少袁
水尧旱尧火灾尧森林病虫害等等频发袁导致了世界大多
数国家的森林生态系统都在不同程度地退化遥 森林
灾害监测为森林受损评价提供基础数据袁为受害森
林恢复的经营管理提供依据遥
目前的森林状况监测多是基于现地实测的袁效
率较低袁无法实现对大区域的快速连续监测遥 随着
空间信息技术的不断进步袁遥感技术在林业中的应
用日益广泛袁越来越多的研究者利用遥感手段定性
摇 第 圆 期 刘丽娟等院 基于时间序列 酝韵阅陨杂 耘灾陨数据的森林生长异常监测
或定量地获取森林参数渊粤灶凿藻则泽燥灶 藻贼 葬造援袁 圆园园远曰 在澡葬燥
藻贼 葬造援袁 圆园园怨曰 李崇贵等袁 圆园园缘曰 庞勇等袁 圆园园愿曰 覃先
林等袁 圆园园缘曰 王海鹏等袁 圆园园愿冤袁以及进行森林健康
监测渊灾燥早藻造皂葬灶灶 藻贼 葬造援袁 员怨愿怨曰 肖文发等袁 圆园园员曰 赵
春燕等袁 圆园园远曰朱振华等袁 圆园园愿冤遥 酝韵阅陨杂中等分辨
率成像光谱仪因其时间分辨率高渊每 员 耀 圆 天观测
地球表面 员 次冤袁可以连续观察地物的变化袁因而多
用于区域尺度的研究渊 在澡葬灶早 藻贼 葬造援袁 圆园园愿曰 方墨人
等袁 圆园园源曰 宫攀等袁 圆园园远曰 王红说等袁 圆园园怨曰 于信芳
等袁 圆园园远曰 李月臣等袁 圆园园缘冤袁但在森林生长异常监
测方面的应用尚不成熟渊耘噪造怎灶凿澡 藻贼 葬造援袁 圆园园怨冤遥
本文尝试采用高时间分辨率的 酝韵阅陨杂 数据对
区域尺度的森林生长状况进行监测袁力图找出生长
季内的异常点袁提高大范围森林生长异常监测的效
率袁以便森林决策者对异常情况做出快速判断并采
取相应措施遥
员摇 研究区概况及数据获取
员郾 员摇 研究区概况
本文选取黑龙江尧吉林尧辽宁 渊员员愿毅猿园忆要员猿远毅
耘袁猿愿毅源园忆要缘源毅晕冤作为主要的研究区域遥 东北林区
地处我国大尧小兴安岭和长白山脉袁生态地位非常重
要遥 土地总面积 苑愿郾 苑猿 万 噪皂圆袁人口约员园 苑员缘郾 源万
人袁属大陆性季风气候袁四季分明袁雨热同期袁日照丰
富袁积温较高袁冬长夏暖袁春秋季短遥 地形以山地尧平
原尧河流为主遥
东北三省的林业用地面积达 猿 源远远郾 源远 万 澡皂圆袁
占土地面积的 源猿郾 怨圆豫 袁森林覆盖率达 猿苑郾 怨怨豫 袁远
远高出全国森林覆盖率 圆园郾 猿远豫的水平曰 该区还是
我国最主要的天然林区袁天然林面积 圆 猿怨圆郾 远猿 万
澡皂圆袁占该区有林地面积的 愿园郾 圆远豫 遥 主要森林类型
为针叶林尧阔叶林和针阔混交林袁大兴安岭北部山地
寒温带针叶林区袁是我国唯一的大片落叶松林区遥
主要森林树种有兴安落叶松渊蕴葬则蚤曾 早皂藻造蚤灶蚤蚤冤尧红松
渊 孕蚤灶怎泽 噪燥则葬蚤藻灶泽蚤泽 冤尧 樟子松 渊 孕蚤灶怎泽 泽赠造增藻泽贼则蚤泽 增葬则援
皂燥灶早燥造蚤糟葬 冤尧 白 桦 渊 月藻贼怎造葬 责造葬贼赠责澡赠造造葬 冤尧 黄 檗
渊 孕澡藻造造燥凿藻灶凿则燥灶 葬皂怎则藻灶泽藻 冤尧 水 曲 柳 渊 云则葬曾蚤灶怎泽
皂葬灶凿泽糟澡怎则蚤糟葬冤尧 鱼 鳞 云 杉 渊 孕蚤糟藻葬 躁藻扎燥藻灶泽蚤泽 增葬则援
皂蚤糟则燥泽责藻则皂葬冤尧山杨 渊 孕燥责怎造怎泽 凿葬增蚤凿蚤葬灶葬 冤尧 蒙古栎
渊匝怎藻则糟怎泽 皂燥灶早燥造蚤糟葬冤尧榆树渊哉造皂怎泽 责怎皂蚤造葬冤等遥
随着经济的发展袁东北地区天然林大量被消耗袁
人工纯林面积不断增大袁生物多样性降低袁森林生态
系统失衡袁导致森林鼠害袁落叶松毛虫尧光肩星天牛尧
青杨虎天牛等病虫害大面积发生曰 并且该地区的森
林覆盖率很高袁林下可燃物堆积厚袁年降雨量小袁导
致 缘 月尧员园 月的火险等级很高曰 加之北部寒潮时有
入侵袁天气多变袁易发生雪灾尧风倒等灾害遥
员郾 圆摇 数据获取
员郾 圆郾 员摇 时间序列 酝韵阅陨杂数据获取摇 本文拟通过遥
感技术连续地观察林木生长状态袁进而进行森林生
长的异常监测袁因此需获取时间序列的遥感影像作
为数据源遥 考虑到研究区覆盖面积较大袁且对时间
分辨率要求较高袁选择 栽藻则则葬鄄酝韵阅陨杂 地表反射率数
据 渊 酝韵阅园怨粤员 冤 作为数据源 渊其空间分辨率为
缘园园 皂袁时间分辨率为 愿 天冤袁时间范围为 圆园园源要
圆园园远 年遥 根据 晕粤杂粤对全球区域的划分袁东北三省
酝韵阅陨杂 影像主要覆盖 澡圆缘增园猿袁 澡圆远增园猿袁 澡圆远增园源袁
澡圆苑增园源 四 景袁 下 载 的 原 始 影 像 为 匀阅云鄄耘韵杂
渊澡蚤藻则葬则糟澡蚤糟葬造 凿葬贼葬 枣燥则皂葬贼鄄藻葬则贼澡 燥遭泽藻则增蚤灶早 泽赠泽贼藻皂冤存储
格式袁投影方式为 陨杂陨晕渊蚤灶贼藻早藻则蚤扎藻凿 泽蚤灶怎泽燥蚤凿葬造冤投影遥
在 酝燥凿蚤泽 砸藻责则燥躁藻糟贼蚤燥灶 栽燥燥造 批处理模式下袁对四
景 酝韵阅陨杂 影像进行镶嵌袁转投影方式为 蕴葬皂遭藻则贼赃
悦燥灶枣燥则皂葬造赃悦燥灶蚤糟投影袁利用东北三省行政区域对影
像裁切袁得到研究区 酝韵阅陨杂影像图渊图 员冤遥
图 员摇 研究区 圆园园源 年 酝韵阅陨杂 圆袁源袁猿
波段假彩色合成影像
云蚤早援 员摇 云葬造泽藻 糟燥造燥则 糟燥皂责燥泽蚤贼藻 蚤皂葬早藻 怎泽蚤灶早
酝韵阅陨杂 遭葬灶凿泽 圆袁 源袁 葬灶凿 猿 蚤灶 圆园园源
员郾 圆郾 圆摇 时间序列植被指数计算摇 植被指数是利用
遥感卫星探测数据的线性或非线性组合而形成的能
反映绿色植被生长状况和分布的特征指数渊陈述彭袁
员怨怨园冤袁而时间序列的植被指数能更准确地反映各个
阶段的植被生长状况遥 归一化植被指数渊晕阅灾陨冤的应
用最广泛且表现稳定袁但由于它在植被高覆盖区有易
饱和现象袁而增强型植被指数渊耘灾陨冤改进了 晕阅灾陨 的
某些缺陷袁特别是大气噪声尧土壤背景尧饱和度等问题
渊匀怎藻贼藻 藻贼 葬造援袁 圆园园圆曰 王正兴等袁 圆园园猿冤袁能更好地反映
高植被覆盖区的植被生长状况袁因而本文选用年时间
序列的 耘灾陨描述林地的生长状况遥
在 陨阅蕴中对时间序列的酝韵阅陨杂影像利用式渊员冤
缘缘
林 业 科 学 源愿 卷摇
渊匀怎藻贼藻 藻贼 葬造援袁 圆园园圆冤计算 耘灾陨袁生成研究区时间序列 的 耘灾陨影像图渊图 圆冤遥
图 圆摇 圆园园源 年第 员袁员园袁圆园袁猿园 期 愿 天合成的 耘灾陨影像
云蚤早援 圆 悦燥皂责燥泽蚤贼藻 耘灾陨 蚤皂葬早藻 枣则燥皂 贼澡藻 员泽贼 袁员园贼澡 袁圆园贼澡袁猿园贼澡 愿鄄凿葬赠 责藻则蚤燥凿 蚤灶 圆园园源袁 泽藻责葬则葬贼藻造赠
耘灾陨 越 郧 籽晕陨砸 原 籽则藻凿籽晕陨砸 垣 悦员籽则藻凿 原 悦圆籽遭造怎藻 垣 蕴遥 渊员冤
式中院 籽遭造怎藻袁籽则藻凿 和 籽晕陨砸 分别为 酝韵阅陨杂 影像的蓝尧红尧
近红外波段袁 蕴 渊 越 员冤为冠层背景调节系数袁 悦员 渊 越远冤袁 悦圆 渊 越 苑郾 缘冤为气溶胶散射修正系数袁 郧 渊 越 圆郾 缘冤
为增益因子遥
圆摇 研究方法
圆郾 员摇 遥感数据预处理
尽管 愿 天最大值合成的 酝韵阅陨杂 地表反射率数
据在一定程度上去除了云和大气的影响袁但仍然可
能受太阳高度角尧观测角度尧气溶胶等一些随机因素
的干扰袁使得时间序列数据的变化呈不规则状态
渊于信芳等袁 圆园园远冤袁通常表现为 耘灾陨 曲线短期起伏
很大袁这是不符合植被正常生长状况的袁因而需要对
合成的 耘灾陨数据进一步滤波袁去除噪声干扰遥
本文选用 杂葬增蚤贼扎噪赠鄄郧燥造葬赠 滤波函数平滑时间序
列 耘灾陨渊杂葬增蚤贼扎噪赠 藻贼 葬造援袁 员怨远源冤遥 杂葬增蚤贼扎噪赠鄄郧燥造葬赠 滤波
是一种最小二乘法卷积平滑技术袁其主要原理为院
对长度为 灶 的窗口邻域内各点数据袁 用一元 责 次
多项式拟合袁 此多项式的系数借助最小二乘准则
来确定袁 卷积核可以单独保存为系数矩阵袁它仅由
滤波器维数尧窗口大小和阶数确定袁多项式在窗口
中心点的值即为经卷积平滑处理后的值遥 滤波过
程需要经过多次试验袁改变窗口大小和阶数等参
数袁既要达到平滑去噪的效果袁又不可过度平滑弱
化了植被原始信息遥 经比较选定窗口长度为 苑 的
一元二次多项式拟合袁从平滑结果看曲线基本保
持原有形状袁但上下剧烈起伏的状况得到了改善
渊图 猿冤袁认为 杂葬增蚤贼扎噪赠鄄郧燥造葬赠 滤波很大程度上去除
了噪声干扰影响袁可以将平滑后的 耘灾陨 值进行后
续研究遥
图 猿摇 杂葬增蚤贼扎噪赠鄄郧燥造葬赠函数平滑前后的 耘灾陨曲线比较
云蚤早援 猿摇 栽澡藻 糟燥皂责葬则蚤泽燥灶 燥枣 耘灾陨 贼蚤皂藻 泽藻则蚤藻泽 糟怎则增藻泽
遭藻枣燥则藻 葬灶凿 葬枣贼藻则 杂葬增蚤贼扎噪赠鄄郧燥造葬赠 枣蚤造贼藻则
圆郾 圆摇 技术路线
生长受到胁迫的森林袁其光合有效率尧光谱特
性尧叶绿素含量尧叶面积指数等往往会有所反应袁生
长季的 耘灾陨值也能够有所体现袁其值一般会低于正
常年份曰 另外年时间序列的曲线上也必然起伏较
大袁复杂度较正常年份高遥 本文拟从这 圆 个角度出
发袁计算 耘灾陨生长季面积渊表征植被的生长状况冤和
耘灾陨曲线的信息熵渊表征曲线的复杂度冤袁联合判定
异常点袁确定灾害的发生遥 由于随机干扰的存在袁单
独从任一方面很难准确判定灾害的发生袁但是二者
的交集一定程度上可以剔除因偶然因素导致的异
常遥 技术路线如图 源 所示遥
圆郾 圆郾 员摇 计算年 耘灾陨生长季面积摇 植被指数的生长
季面积可以反映植被年内的生长状况遥 东北三省的
常绿林所占比例较小袁且其年内植被指数值变化不
大袁因而本文仅考虑主要的森林类型要要要落叶林和
针阔混交林袁计算生长季面积院
粤砸耘粤生长季 越 乙绿色减到最少
绿色开始出现
耘灾陨 遥 渊圆冤
摇 摇 根据式渊圆冤袁首先确定积分限袁即绿色开始出现
远缘
摇 第 圆 期 刘丽娟等院 基于时间序列 酝韵阅陨杂 耘灾陨数据的森林生长异常监测
图 源摇 基于 酝韵阅陨杂时间序列数据的森林
生长异常监测技术路线
云蚤早援 源摇 栽澡藻 枣造燥憎糟澡葬则贼 燥枣 枣燥则藻泽贼 早则燥憎贼澡 凿蚤泽贼怎则遭葬灶糟藻
皂燥灶蚤贼燥则蚤灶早 怎泽蚤灶早 贼蚤皂藻 泽藻则蚤藻泽 燥枣 酝韵阅陨杂 凿葬贼葬
的时间和绿色减到最少的时间袁进而计算这一时间
段 耘灾陨的积分遥 可见生长季开始和结束时间的算
法决定着生长季面积计算的精度遥 本文采用 在澡葬灶早
等渊圆园园猿曰 圆园园远冤经过现地验证过的方法确定积分
上尧下限院 每个独立的生长周期都包含绿色开始出
现尧绿色增加到最大值尧绿色开始减少尧绿色减少到
最小值 源 个关键时刻点袁通过 杂 形函数渊蕴燥早蚤泽贼蚤糟 函
数冤分段拟合植被年生长曲线院
耘灾陨渊 贼冤 越 糟员 垣 藻葬垣遭贼 垣 凿遥 渊猿冤
式中院耘灾陨渊 贼冤为 贼 时刻的 耘灾陨 值曰 葬袁遭 为拟合参数曰
糟 垣 凿为 耘灾陨的最大值曰 凿为初始背景值遥
计算该生长曲线的曲率及曲率的一阶导数袁即
曲率的变化率 运忆 咱渊式渊源冤暂袁得到 源 个极值点袁这 源
个极值点就是对应的植被生长的 源 个转折期袁其中
第 员袁源 点分别对应积分的下限和上限遥
运忆 越 遭猿糟扎 {猿扎渊员 原 扎冤渊员 垣 扎冤猿咱圆渊员 垣 扎冤猿 垣 遭圆糟圆 扎暂咱渊员 垣 扎冤源 垣 渊遭糟扎冤圆暂 缘圆 原
渊员 垣 扎冤圆渊员 垣 圆扎 原 缘扎圆冤
咱渊员 垣 扎冤源 垣 渊遭糟扎冤圆暂 猿圆 } 遥 渊源冤
式中院 扎 越 藻葬垣遭贼遥
受灾害条件的影响袁不同年份同一像元植被生
长季的起止时间咱式渊圆冤积分上下限暂及长度会有所
不同袁在计算生长季面积和曲线熵值时袁时间范围如
何选取是个关键的问题遥 如果假设其中一年是正常
年袁并将该年的生长季起止时间作为这 猿 年统一的
计算区间袁则会受假设正常年的影响袁一定程度上弱
化了受干扰像元的信息袁导致该像元不能被正常检
测出来曰 但如果同一像元取各自年内的生长季起止
时间作为计算区间袁可能会将由于全年气候影响导
致的生长季起止时间发生变化的像元误检测袁同
时会降低年际间的可比性遥 因此袁为了 猿 年内同
一像元间更具有可比性袁本文生长季的起止时间
选择同一像元 猿 年的最晚发芽和最早落叶的
时间遥
圆郾 圆郾 圆摇 计算年 耘灾陨曲线熵值摇 杂澡葬灶灶燥灶渊员怨源愿冤提出
了信息熵的概念袁信息熵是信息论中用于度量信息
量的一个概念遥 一个系统越是有序袁信息熵就越低曰
不确定性越大袁信息熵就越高遥 信息熵可以用来描
述一个广义集合的内部状态的复杂程度袁二者成正
比例关系遥 根据香农信息论袁本文中将信息熵概念
推广到用于表达曲线的复杂程度袁见式渊缘冤遥
匀 越 原 鄱灶蚤 越 员 责渊 蚤冤造燥早圆责渊 蚤冤袁
责渊 蚤冤 越 枣渊 蚤冤 辕鄱灶蚤 越 员 枣渊 蚤冤遥 渊缘冤
式中院 枣渊 蚤冤 表示第 蚤点的 耘灾陨值曰 责渊 蚤冤 表示 枣渊 蚤冤 出
现的概率遥 熵值计算的起止时间选择与生长季相同
的起止时间遥
曲线的熵值越大袁表明复杂度越高曰反之袁熵值
越小袁复杂度越低遥 因此可以通过比较年 耘灾陨 曲线
的熵值来确定曲线的变异程度袁从而判断是否发生
异常遥
圆郾 圆郾 猿摇 生长季面积和熵值联合监测异常点摇 假设
圆园园源 年某林地受到灾害影响发生异常袁那么一般其
后的 员 耀 圆 年内该处植被生长状况会比受灾之前差袁
对应像元的年时间序列 耘灾陨 曲线的复杂度会高于
正常年袁而该年的生长季面积一般会小于正常年
渊受干扰程度不同袁表现有所不同冤袁即院 熵值越大袁
生长季面积越小袁该像元发生异常的可能性越大袁因
而可以根据 猿 年间生长季面积尧熵值的变化来联合
监测异常点遥
分别计算 猿 年时间序列 耘灾陨数据的熵值和生长
季面积袁并求其变化量袁由于每年的气候条件不同袁
导致生长状况有差异袁因而需通过设定阈值方法袁剔
除一些因气候等因素影响而异常的点或变化量相对
较小的点曰 然后取熵值和生长季面积变化量的交
集袁从而确定每年发生异常的像元曰 画出 猿 年异常
点的 耘灾陨曲线袁进行比较验证遥
猿摇 结果与分析
猿郾 员摇 年生长季面积
愿 天合成的年时间序列 耘灾陨影像袁每个像素都
对应着 源远 个值袁根据式渊圆冤和式渊猿冤对这 源远 个离
散值用 蕴燥早蚤泽贼蚤糟生长方程拟合年生长曲线袁并计算
苑缘
林 业 科 学 源愿 卷摇
该生长曲线的曲率变化率袁得到的 源 个极值点分
别对应一年当中绿色开始出现尧绿色增加到最大
值尧绿色开始减少尧绿色减少到最小值的时间点遥
经过计算袁圆园园源袁圆园园缘袁圆园园远 年林地绿色开始出现
的平均时间分别是 员员怨袁员圆缘 和 员员怨 天袁绿色减少
到最小值的平均时间分别是 圆愿远袁圆怨愿 和 猿园员 天遥
圆园园缘 年绿色开始出现时间较另外 圆 年稍晚袁初步
推断可能是 圆园园源 年森林受到了干扰袁导致次年展
叶时间推迟袁圆园园远 年干扰程度轻的植被又恢复了
正常曰 而一定程度的干扰可能会促进林下植被的
更新袁增强其活力袁从而导致平均生长季略有延
长袁但不排除 圆园园缘 年由于气候等原因造成的绿色
开始出现时间较晚袁因而还不能仅凭关键时间点
的变化判断是否发生异常遥
然后对绿色开始出现至绿色减少到最小值的时
间段积分袁得到年 耘灾陨的生长季面积渊图 缘冤遥
图 缘摇 圆园园源要圆园园远 年生长季面积
云蚤早援 缘摇 栽澡藻 耘灾陨 葬则藻葬 凿怎则蚤灶早 早则燥憎贼澡 泽藻葬泽燥灶 蚤灶 圆园园源袁 圆园园缘袁 葬灶凿 圆园园远袁 泽藻责葬则葬贼藻造赠
摇 摇 从生长季面积图上可以明显地区分出林区及
非林区袁由于在生长季森林的植被指数值明显高
于农田尧河流以及建设用地等非林区袁因此累积的
生长季面积值大袁图上表现为亮度值明显高于非
林区遥
东北三省境内有大尧小兴安岭尧长白山脉及千山
山脉袁林区纵贯南北遥 总体来看袁小兴安岭至长白山
区袁森林类型是以红松为主的针阔混交林袁其生长季
面积值相对较高曰 而黑龙江省北部大兴安岭地区是
以兴安落叶松为主的针叶林区袁生长季面积值相对
较低曰 辽东山区以次生林为主袁森林类型多为油松
渊孕蚤灶怎泽 贼葬遭怎造葬藻枣燥则皂蚤泽冤和落叶阔叶混交林袁但由于该
区域内的千山山脉是长白山支脉袁部分植被与长白
山脉又有相似之处袁因而其生长季面积略低于小兴
安岭林区但高于大兴安岭林区遥
猿郾 圆摇 年 耘灾陨曲线熵值
根据式渊源冤逐年计算 耘灾陨曲线生长季的熵值并
制图渊图 远冤遥
图 远摇 圆园园源要圆园园远 年生长季 耘灾陨曲线熵值
云蚤早援 远摇 栽澡藻 藻灶贼则燥责赠 凿怎则蚤灶早 早则燥憎贼澡 泽藻葬泽燥灶 蚤灶 圆园园源袁 圆园园缘袁 葬灶凿 圆园园远袁 泽藻责葬则葬贼藻造赠
摇 摇 曲线熵值是曲线复杂度的表达袁森林植被季节
性生长经历增加尧生长盛期尧凋落过程袁因而森林植
被的曲线复杂度要高于其他地物遥 但总体上看袁林
区与非林区生长季熵值不如生长季面积的差别明
愿缘
摇 第 圆 期 刘丽娟等院 基于时间序列 酝韵阅陨杂 耘灾陨数据的森林生长异常监测
显袁这主要是由于非林区地物 耘灾陨 值的波动会增大
曲线的熵值袁但对总的 耘灾陨面积影响并不大遥
比较 猿 年每一点的熵值大小袁可以直观地看到
一些点在某一年熵值明显不同于其他 圆 年袁熵值越
大说明 耘灾陨曲线越复杂袁受灾害干扰的可能性也越
大袁但具体定量的分析还需要结合生长季面积来联
合判定是否受到了干扰遥
猿郾 猿摇 异常点提取
由于气温尧降水等气候因素影响袁正常生长的植
被年际间也会有变化袁只是其变化量会比受到干扰
时的小袁因而需设定相应阈值去除变化量较小的点袁
认为变化量大的点才是可能发生异常的点遥
根据式渊远冤计算 圆园园源要圆园园远 年间生长季面积
渊葬则藻葬冤和熵值渊藻灶贼则燥责赠冤的变化量袁参考该地区关于
火灾渊杜永胜等袁 圆园园苑冤的历史记载渊其他森林灾害
事件的历史数据暂时缺乏冤袁 进行阈值的设定袁其原
则是既保证过火面积的最大化提取袁又保留其他的
异常信息遥 经过反复试验袁最终确定 圆园园源要圆园园远 年
熵值增加范围 咱园郾 园员愿 苑袁 园郾 园愿 暂袁面积减少范围
咱员 园园园袁 圆 怨园园暂遥 再取生长季熵值增加和面积减少
部分的交集袁即同时满足生长季 耘灾陨 累积值减少且
年 耘灾陨曲线复杂度高的点为最终提取的异常点袁并
制图渊图 苑冤遥
渊 耘灶贼则燥责赠圆园园远 原 耘灶贼则燥责赠圆园园源 跃 栽澡则藻泽澡燥造凿员冤 疑
渊 粤则藻葬圆园园远 原 粤则藻葬圆园园源 跃 栽澡则藻泽澡燥造凿圆冤遥 渊远冤
图 苑摇 圆园园源要圆园园远 年森林生长异常监测结果
云蚤早援 苑摇 酝葬责责蚤灶早 燥枣 枣燥则藻泽贼 早则燥憎贼澡 凿蚤泽贼怎则遭葬灶糟藻 燥枣 圆园园源要圆园园远
猿郾 源摇 异常点曲线分析
如图 愿袁怨 所示为正常点与异常点 猿 年 耘灾陨曲线
的对比遥 可以看出正常点的 耘灾陨 曲线整体波动不
大袁生长季长度尧峰值差异并不显著袁但异常点曲线
可以明显看出年际间峰值有很大的波动袁说明可以
通过本文所述方法提取出受干扰的像元袁并做进一
步分析遥
随机抽取异常点画出 猿 年的 耘灾陨曲线渊图 员园冤袁
以及对应生长季面积和熵值袁比较其曲线形状尧峰
值袁并验证本文方法所求熵值和生长季面积是否与
曲线形状相符遥
图 员园 中列举了 圆 个异常点 猿 年的 耘灾陨 变化曲
线袁可以看出其中 员 年或 圆 年明显出现了异常袁圆园园源
年曲线形状符合正常的植被生长曲线院 耘灾陨 逐渐增
怨缘
林 业 科 学 源愿 卷摇
图 愿摇 正常点 耘灾陨曲线
云蚤早援 愿摇 栽澡藻 灶燥则皂葬造 耘灾陨 糟怎则增藻
图 怨摇 异常点 耘灾陨曲线
云蚤早援 怨摇 栽澡藻 葬遭灶燥则皂葬造 耘灾陨 糟怎则增藻
加袁达到最大值后趋于平缓袁持续一段时间后又开始
下降曰 圆园园缘 年曲线形状基本正常袁但峰值有所下降曰
而 圆园园远 年曲线明显异于前 圆 年院 峰值低袁并且在生
图 员园摇 猿 年的 耘灾陨曲线异常点比较
云蚤早援 员园摇 栽澡藻 糟燥皂责葬则蚤泽燥灶 燥枣 糟怎则增藻泽 枣则燥皂 葬遭灶燥则皂葬造 责燥蚤灶贼泽 蚤灶 贼澡藻 贼澡则藻藻 赠藻葬则泽
长旺季出现了双峰现象曰在 圆园园缘 年生长旺盛也出现
了双峰现象袁而 圆园园远 年峰值变化更加明显袁在生长
阶段就出现了下降趋势袁然后又缓慢上升达到峰值袁
但峰值阶段持续时间较短袁再进行熵值和生长季面
积的年际间比较院 圆园园远 年由于双峰的出现导致曲线
复杂度高袁曲线熵值明显大于前 圆 年袁另外峰值的降
低以及较短的峰值持续时间袁导致了生长季面积达
到 猿 年中最低值遥 可见袁文中所述方法提取的生长
季面积及曲线熵值可以反映曲线的趋势袁作为异常
点提取的指标是可行的遥 但具体的干扰因子尚不明
确袁如有的 耘灾陨出现了双峰袁峰值部分趋势为上升 原
下降 原上升袁原因有可能是生长季早期发生了病虫
害袁生长季中后期逐渐恢复袁耘灾陨 又呈上升趋势曰 有
的 耘灾陨达到最大值的时间较晚袁峰值较正常年小袁
并且持续时间较短后就开始下降袁这种情况分析可
能是前一年发生了异常袁导致下一年生长季增长缓
慢等袁因此具体的受干扰类型还需进一步的验证遥
猿郾 缘摇 火点验证
为了明确发生异常的原因袁以火点为例袁选取
圆园园源 年 员园 月中旬黑龙江省黑河地区发生的林火灾
害事件袁对本文提取异常点进行验证遥 用于验证的
火点数据来源于基于 酝韵阅陨杂 提取的火情监测指数
晕阅栽陨渊灶燥则皂葬造蚤扎藻凿 凿蚤枣枣藻则藻灶糟藻 贼澡藻则皂葬造 蚤灶凿藻曾冤遥 晕阅栽陨 的
计算需要 酝韵阅陨杂 专门为火情监测设置的 源 滋皂 和
员员 滋皂通道辐射数据渊酝韵阅园圆员运酝冤以及定标数据
渊酝韵阅园猿冤遥 该指数经过现地验证 渊赵文化等袁
圆园园愿冤袁检测精度较高袁并且计算速度要比 酝韵阅陨杂 推
荐火情算法快 圆 倍以上袁可以快速实现火情提取遥
图 员员 为截取局部的 晕阅栽陨 提取火点与本文提
取异常点的比较遥 总体上看袁圆 种方法提取的一致
性相对较高袁但由于本文提取的异常点中除火点外袁
还包含由于其他异常原因渊如病虫害冤导致的异常袁
因此异常点数稍多于 晕阅栽陨 的提取数袁观察多提取
点的 耘灾陨曲线袁发现有些点与火点的曲线形状非常
相近袁但也有部分点相异较大遥 分析其原因院员冤 由
于 晕阅栽陨的提取是基于单一时相袁有些火点在当天
园远
摇 第 圆 期 刘丽娟等院 基于时间序列 酝韵阅陨杂 耘灾陨数据的森林生长异常监测
未能检测但在时间序列中得以反映曰 圆冤 火场附近
林地由于燃烧导致生境改变袁生物多样性降低袁从而
易发生病虫害等灾害袁但具体异常原因仍有待于结
合其他数据辅助验证遥
图 员员摇 不同于 晕阅栽陨提取火点的曲线
云蚤早援 员员摇 栽澡藻 糟怎则增藻泽 凿藻贼藻糟贼藻凿 葬遭灶燥则皂葬造 责燥蚤灶贼 憎澡蚤糟澡 憎葬泽 灶燥贼 凿藻贼藻糟贼藻凿 遭赠 晕阅栽陨
源摇 结论与讨论
本文通过对东北地区 圆园园源要圆园园远 年 猿 年间的
酝韵阅陨杂 耘灾陨 数据的分析袁结果表明时间序列的
酝韵阅陨杂 耘灾陨数据可以用来进行森林生长异常制图遥
高时间分辨率的植被指数对森林生长状态敏感袁甚
至可以捕捉到细微的异常变化袁一定程度上弥补了
单一时相不能很好反映全年植被生长变化的不足遥
针对时间序列数据的特点袁文中将植被生长状
况渊生长季 耘灾陨面积冤和时间序列 耘灾陨特征曲线渊生
长季熵值冤的形态变化相结合进行异常分析袁而不
是局限于常规的基于分类的变化检测方法遥 由于这
圆 个指标都与植被的生长状况密切相关袁并且本文
精确求取了生长季的时间范围袁保证了指标提取的
精度袁因而一旦发生了灾害袁指标上一定会有所反
映袁并结合适合的算法可以将异常点检测出来遥
通过对东北三省林区的试验袁本文提出的森林
生长异常监测方法是基本可行的袁克服了依赖现地
调查的不足袁特别适用于区域尺度袁检测到的森林生
长异常区域的 耘灾陨 特征曲线也对异常的原因有一
定的揭示曰 观察异常点时间序列的 耘灾陨 曲线发现院
曲线的熵值明显大于正常年袁生长季峰值低袁导致生
长季面积小于正常年袁并且在生长旺季会出现峰值
突然持续下降或出现双峰等异常现象遥 目前由于只
员远
林 业 科 学 源愿 卷摇
对火点进行了初步的验证袁虽然数据进行了滤波处
理袁但可能仍然会受到云的干扰导致误检测袁因此下
一步工作将利用 酝韵阅园怨粤员 数据的 匝粤产品袁去除有
云像元袁结合东北林区的森林灾害事件对本方法进
行进一步的验证和改进袁并分析曲线规律袁找出各种
异常的原因遥
参 考 文 献
陈述彭援 员怨怨园郾 遥感大辞典援北京院 科学出版社袁猿愿圆援
方墨人袁田庆久援 圆园园源郾 基于 酝韵阅陨杂 的 蕴粤陨 时间序列谱的地物分类
方法研究援国土资源遥感袁渊猿冤院 缘 原 愿援
宫摇 攀袁唐华俊袁陈仲新袁等援 圆园园远郾 基于 酝韵阅陨杂 蕴杂栽 修正 晕阅灾陨 时序
数列的土地覆盖分类援资源科学袁 圆愿渊源冤院 员园源 原 员员园援
李崇贵袁赵宪文援 圆园园缘郾 森林郁闭度定量估测遥感比值波段的选择援
林业科学袁源员渊源冤院 苑圆 原 苑苑援
李月臣袁陈摇 晋袁宫摇 鹏援 圆园园缘郾 基于 晕阅灾陨时间序列数据的土地覆盖
变化 检 测 指 标 设 计援 应 用 基 础 与 工 程 科 学 学 报袁
员猿渊猿冤院 圆远员 原 圆苑缘援
杜永胜袁王立夫援 圆园园苑郾 中国森林火灾典型案例援北京院 中国林业出版
社袁 员缘缘 原 员远员援
庞摇 勇袁李增元袁谭炳香袁等援 圆园园愿郾 点云密度对机载激光雷达林分高
度反演的影响援林业科学研究袁 圆员渊增刊冤院 员源 原 员怨援
覃先林袁李增元袁易浩若援 圆园园缘郾 高空间分辨率卫星遥感影像树冠信
息提取方法研究援遥感技术与应用袁 圆园渊圆冤院 圆圆愿 原 圆猿圆援
王海鹏袁金亚秋袁大内和夫袁等援 圆园园愿郾 孕蚤赃杂粤砸极化数据与 运分布指数
估算森林生物量与实验验证援遥感学报袁 员圆渊猿冤院 源苑苑 原 源愿圆援
王红说袁黄敬峰援 圆园园怨郾 基于 酝韵阅陨杂 晕阅灾陨时间序列的植被覆盖变化
特征 研 究袁 浙 江 大 学 学 报院 农 业 与 生 命 科 学 版袁 猿缘
渊员冤院 员园缘 原 员员园援
王正兴袁刘摇 闯袁匀怎藻贼藻 粤 砸援 圆园园猿郾 植被指数研究进展院 从 粤灾匀砸砸鄄
晕阅灾陨 到 酝韵阅陨杂鄄耘灾陨援生态学报袁 圆猿渊缘冤院 怨苑怨 原 怨愿苑援
肖文发袁韩景军袁马摇 娟援 圆园园员郾 美国国家森林健康监测与评价计划
及对我国的启示援世界林业研究袁 员源渊猿冤院 远苑 原 苑源援
于信芳袁庄大方援 圆园园远郾 基于 酝韵阅陨杂 晕阅灾陨数据的东北森林物候期监
测援资源科学袁 圆愿渊源冤院 员员员 原 员员苑援
赵春燕袁杨志高援 圆园园远郾 高光谱遥感技术及其在森林监测中的应用探
讨援林业调查规划袁猿员渊员冤院 源 原 远援
赵文化袁单海滨袁钟儒祥援 圆园园愿郾 基于 酝韵阅陨杂 火点指数监测森林火
灾援自然灾害学报袁 员苑渊猿冤院 员缘圆 原 员缘苑援
朱振华袁银晓瑞袁章海波袁等援 圆园园愿郾 基于遥感影像的鄂托克前旗森林
资源监测分析援内蒙古林业科技袁猿源渊猿冤院 远猿 原 远远援
粤灶凿藻则泽燥灶 允袁 酝葬则贼蚤灶 酝 耘袁 杂皂蚤贼澡 酝 蕴袁藻贼 葬造援 圆园园远郾 栽澡藻 怎泽藻 燥枣 憎葬增藻枣燥则皂
造蚤凿葬则 贼燥 皂藻葬泽怎则藻 灶燥则贼澡藻则灶 贼藻皂责藻则葬贼藻 皂蚤曾藻凿 糟燥灶蚤枣藻则 葬灶凿 凿藻糟蚤凿怎燥怎泽
枣燥则藻泽贼 泽贼则怎糟贼怎则藻 蚤灶 晕藻憎 匀葬皂责泽澡蚤则藻援 砸藻皂燥贼藻 杂藻灶泽蚤灶早 燥枣 耘灶增蚤则燥灶皂藻灶贼袁
员园缘渊猿冤院 圆源愿 原 圆远员援
耘噪造怎灶凿澡 蕴袁 允燥澡葬灶泽泽燥灶 栽袁 杂燥造遭藻则早 杂援 圆园园怨郾 酝葬责责蚤灶早 蚤灶泽藻糟贼 凿藻枣燥造蚤葬贼蚤燥灶 蚤灶
杂糟燥贼泽 责蚤灶藻 憎蚤贼澡 酝韵阅陨杂 贼蚤皂藻鄄泽藻则蚤藻泽 凿葬贼葬援 砸藻皂燥贼藻 杂藻灶泽蚤灶早 燥枣
耘灶增蚤则燥灶皂藻灶贼袁 员员猿渊苑冤院员缘远远 原 员缘苑猿援
匀怎藻贼藻 粤袁 阅蚤凿葬灶 运袁 酝蚤怎则葬 栽袁 藻贼 葬造援 圆园园圆郾 韵增藻则增蚤藻憎 燥枣 贼澡藻 则葬凿蚤燥皂藻贼则蚤糟
葬灶凿 遭蚤燥责澡赠泽蚤糟葬造 责藻则枣燥则皂葬灶糟藻 燥枣 贼澡藻 酝韵阅陨杂 增藻早藻贼葬贼蚤燥灶 蚤灶凿蚤糟藻泽援
砸藻皂燥贼藻 杂藻灶泽蚤灶早 燥枣 耘灶增蚤则燥灶皂藻灶贼袁 愿猿渊员 辕 圆冤院 员怨缘 原 圆员猿援
杂葬增蚤贼扎噪赠 粤袁郧燥造葬赠酝 允 耘援 员怨远源郾 杂皂燥燥贼澡蚤灶早 葬灶凿 凿蚤枣枣藻则藻灶贼蚤葬贼蚤燥灶 燥枣 凿葬贼葬 遭赠
泽蚤皂责造蚤枣蚤藻凿 造藻葬泽贼 泽择怎葬则藻泽 责则燥糟藻凿怎则藻泽援 粤灶葬造赠贼蚤糟葬造 悦澡藻皂蚤泽贼则赠袁 猿远渊愿冤院
员远圆苑 原 员远猿怨援
杂澡葬灶灶燥灶 悦 耘援 员怨源愿郾 粤 皂葬贼澡藻皂葬贼蚤糟葬造 贼澡藻燥则赠 燥枣 糟燥皂皂怎灶蚤糟葬贼蚤燥灶援 栽澡藻 月藻造造
杂赠泽贼藻皂 栽藻糟澡灶蚤糟葬造 允燥怎则灶葬造袁 圆苑院 猿苑怨 原 源圆猿援
灾燥早藻造皂葬灶灶 允 耘袁 砸燥糟噪 月 晕袁 藻贼 葬造援 员怨愿怨郾 哉泽藻 燥枣 贼澡藻皂葬贼蚤糟 皂葬责责藻则 凿葬贼葬 枣燥则
贼澡藻 凿藻贼藻糟贼蚤燥灶 燥枣 枣燥则藻泽贼 凿葬皂葬早藻 糟葬怎泽藻凿 遭赠 贼澡藻 责藻葬则 贼澡则蚤责泽援 砸藻皂燥贼藻
杂藻灶泽蚤灶早 燥枣 耘灶增蚤则燥灶皂藻灶贼袁 猿园渊猿冤院 圆员苑 原 圆圆缘援
在澡葬灶早 载袁 云则蚤藻凿造 酝 粤袁 杂糟澡葬葬枣 悦 月袁 藻贼 葬造援 圆园园猿郾 酝燥灶蚤贼燥则蚤灶早 增藻早藻贼葬贼蚤燥灶
责澡藻灶燥造燥早赠 怎泽蚤灶早 酝韵阅陨杂援 砸藻皂燥贼藻 杂藻灶泽蚤灶早 燥枣 耘灶增蚤则燥灶皂藻灶贼袁 愿源渊猿冤院
源苑员 原 源苑缘援
在澡葬灶早 载袁 云则蚤藻凿造 酝 粤袁 杂糟澡葬葬枣 悦 月袁 藻贼 葬造援 圆园园远郾 郧造燥遭葬造 增藻早藻贼葬贼蚤燥灶
责澡藻灶燥造燥早赠 枣则燥皂 皂燥凿蚤泽院 藻增葬造怎葬贼蚤燥灶 燥枣 早造燥遭葬造 责葬贼贼藻则灶泽 葬灶凿 糟燥皂责葬则蚤泽燥灶
憎蚤贼澡 蚤灶 泽蚤贼怎 皂藻葬泽怎则藻皂藻灶贼泽援 允燥怎则灶葬造 燥枣 郧藻燥责澡赠泽蚤糟葬造 砸藻泽藻葬则糟澡袁 员员员
渊郧园源园员苑冤院 员 原 员源援
在澡葬灶早 载袁 杂怎灶 砸袁 在澡葬灶早 月袁 藻贼 葬造援 圆园园愿郾 蕴葬灶凿 糟燥增藻则 糟造葬泽泽蚤枣蚤糟葬贼蚤燥灶 燥枣 贼澡藻
晕燥则贼澡 悦澡蚤灶葬 孕造葬蚤灶 怎泽蚤灶早 酝韵阅陨杂 耘灾陨 贼蚤皂藻 泽藻则蚤藻泽援 陨杂孕砸杂 允燥怎则灶葬造 燥枣
孕澡燥贼燥早则葬皂皂藻贼则赠 驭 砸藻皂燥贼藻 杂藻灶泽蚤灶早袁 远猿渊源冤院 源苑远 原 源愿源援
在澡葬燥 运袁 杂燥则蚤灶 孕援 圆园园怨郾 蕴蚤凿葬则鄄遭葬泽藻凿 皂葬责责蚤灶早 燥枣 造藻葬枣 葬则藻葬 蚤灶凿藻曾 葬灶凿 蚤贼泽
怎泽藻 枣燥则 增葬造蚤凿葬贼蚤灶早 郧蕴韵月悦粤砸月韵晕 泽葬贼藻造造蚤贼藻 蕴粤陨 责则燥凿怎糟贼 蚤灶 葬
贼藻皂责藻则葬贼藻 枣燥则藻泽贼 燥枣 贼澡藻 泽燥怎贼澡藻则灶 哉杂粤援 砸藻皂燥贼藻 杂藻灶泽蚤灶早 燥枣
耘灶增蚤则燥灶皂藻灶贼袁 员员猿渊愿冤院 员远圆愿 原 员远源缘援
渊责任编辑摇 石红青冤
圆远