注册 登录  
 加关注
查看详情
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

面朝大海 春暖花开

一点心得,转载本博客文章,请注明转帖,谢谢!

 
 
 

日志

 
 
关于我

中国科学院博士 主要从事遥感机理、定量反演、数据处理以及GIS应用研究。ArcGIS、Envi 、ERDAS、Ecognition软件、IDL语言、6S、SAIL

多光谱遥感对全球陆地冰空间监测的贡献(转)  

2012-02-15 15:09:37|  分类: 遥感 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
http://hi.baidu.com/joeyzone/blog/item/096b4caf68a4c5fbfaed5071.html

2011-09-15 19:43:15|  分类: 冰川 |  标签: |字号大中小 订阅

English name: Multispectral imaging contributions to global land ice measurements from space

Author: Jeffrey S. Kargel, Michael J. Abrams, Michael P. Bishop, Andrew Bush,

Gordon Hamilton, Hester Jiskoot, Andreas Kaab, Hugh H. Kieffer, Ella M. Lee,

Frank Paul, Frank Rau, Bruce Raup, John F. Shroder, Deborah Soltesz,

David Stainforth, Leigh Stearns, Rick Wessels

Magazine: Remote Sensing of Environment

Date: Remote Sensing of Environment 99(2005) 187-219

Received 8 March 2005; received in revised form 30 June 2005; accepted 5 July 2005

《多光谱遥感对全球陆地冰空间监测的贡献》阅读报告

一、摘要翻译

全球陆地冰空间监测组织(Global Land Ice Measurements from space:GLIMS)是一个以获取全球冰川卫星遥感影像,对获得影像加以分析以获得冰川范围和变化信息,并对变化的数据对人类和环境的影响与意义加以评估为目的而建立起来的国际组织。尽管GLIMS利用了多种不同的遥感系统,ASTER(Advanced Spaceborne Thermal Emission and reflection Radiometer:先进热发射与反射辐射计)系统是对于许多的观测任务而言最优的系统,这些任务包括冰川边界制图、物质表面监测、冰川表面动态监测(例如流动矩形区域和冰湖演化)。GLIMS研发的软件适合于多种冰川监测工作,包括:对净雪和对岩屑覆盖冰川的制图;地形分类(强调雪、冰、水、冰与岩屑混合物的区分);多时态变化分析;影像与派生数据的可视化;派生数据的解译和存档。一个全球冰川数据库已由国家雪冰数据中心(National Snow and Ice Data Center: NSIDC, Boulder, Colorado)设计;相关的参数由全球冰川目录(World Glacier Inventory: WGI)扩展而来并与之相兼容。相应的技术工作将在本文中得到概略介绍,细节则讲在其他地方加以介绍。在本文中我们所展示的内容从属于一个大的问题:ASTER与其他卫星多光谱遥感数据如何被用于制图、监测和发现冰川状态和动态的特征,如何通过这些特征理解冰川变化对20世纪及21世纪全球气候变化的反映?我们所列举的示例结果在冰川学和气候学上并不是有充分代表性的。我们的初步的结果,在展现问题的复杂性的同时,在整体上与冰川学学术界的结论(气候变化引起了全球冰川的反应,主要表现为冰川的退缩)是相一致的。无论单个的冰川个体是前进抑或退缩,整体的综合的冰川变化的结果必然是由气候变化所引起的,而其他非气候因素的偏离是可以综合消除的。我们在研究中发现了冰川(对气候变化)响应行为的区域空间模式;这些模式可以被归因于气候变化差异和不同地区冰川在大小和响应时间上的差异。在许多情况下,冰川长度变化从较低的程度上反映了冰川消融,因为冰川变薄(有时并不伴随有即时的长度变化)也是相当重要的。扩展的系统化且统一的针对更多冰川的分析对于今后的工作是必须的,这些工作包括:分离冰川响应中对气候变化响应和对非气候的偏离的响应组分;对冰川变化区域间差异的定量分析;对未来区域冰川变化(特别是与水资源,冰川灾害和全球海平面变化)的预测。我们期望这些完整的监测与评估(将逐步完成)会帮助我们区分成功的气候变化模型(能够解释近期的冰川变化和变化模式,并且能够对未来的变化加以描述)和不成功的模型。

二、引言
1、问题的提出

现有的研究表明,全球范围而言,冰川处于退缩中。研究表明,许多小冰川已经消失,即使是两极的冰盖也出现了大范围的消退。通过对冰川变化状态的监测和研究,可以加强对冰川灾害的监控,为干旱地区的水资源的认识利用和管理提供基础数据,促进对全球海平面变化的研究,并且对于全球气候变化和气候模型的构建有着重要的意义。

目前对与冰川研究的方法主要是野外实地观测数据和一些遥感卫星数据的分析。但是就目前为止,还只有较少的冰川得到比较完整的监测、研究与分析。现有的数据并不足以很好的支持对全球冰川变化状态的研究。到目前为止的利用多光谱卫星影像数据,特别ASTER数据对冰川的变化进行研究,并且构建全球冰川数据库的工作及其进展是本篇论文将要介绍的。

2、研究现状
(1)冰川变化的复杂性

通过对全球冰川变化综合数据的统计学分析,我们可以得出结论目前全球冰川物质平衡处在负平衡状态上。但是对单独一个冰川而言,其变化状态可能是前进抑或是退缩的。由于气候变化的非均一性,以及不同冰川间环境、变化过程和其他冰川属性的差异,冰川变化的状态是非常复杂的。在没有对全球冰川进行完整的统一的监测,并对其变化进行评估前,凭借目前的冰川数据并不足以在冰川演化与气候变化间构建一个特定的联系。

(2)实地考察工作

实地的冰川考察能够提供最为完整的、最可信赖的冰川信息。但是由于冰川数量巨大,且都位于偏远的难于到达的地区,实地考察工作开展相当困难并且耗资巨大,因而不可能对所有的冰川都进行实地考察。结合实地考察结果与遥感数据是对全球陆地冰变化进行长期的完整的监测的主要手段。

(3)现有的基于遥感数据的冰川变化分析研究

A、现有的遥感冰川监测数据源简介

现有的遥感冰川监测中大量使用了星载多光谱可见光/近红外多光谱影像,例如ASTER、ETM+、和SPOT影像。同时机载或星载的雷达和激光雷达数据也在冰川遥感监测中发挥重要作用。

B、ASTER数据介绍

ASTER项目的主要目的是深入了解地球表面或近地面以及较低大气层发生的各种局部和区域尺度过程,其中包括地表和大气的相 互作用。在ASTER主页上明确指出ASTER影像在冰川、水文、城市扩展、火山预报、蒸散/地表温度、地质等六个方面有着广阔的应用前景。

ASTER影像主要技术参数

波段序号 波长范围(μm) 地面分辨率(m)

1  0.52~0.60  15

2  0.63~0.69  15

3  0.76~0.86  15

4  1.60~1.70  30

5  2.145~2.185  30

6  2.185~2.225  30

7  2.235~2.285  30

8  2.295~2.365  30

9  2.36~2.43  30

10  8.125~8.475  90

11  8.475~8.825  90

12  8.925~9.275  90

13  10.25~10.95  90

14  10.95~11.65 90

比较ASTER与 TM不难发现:无论在可见光/近红外、短波红外还是在热红外部分,ASTER的光谱分辨率都高于后者,在近红外部分,ASTER的空间分辨率与光谱分辨率 都高于TM、SPOT影像。ASTER的时间分辨率为1次/15天。与TM相比,ASTER数据的价格亦低。

ASTER影像数据具有同轨道立体观测能力,这一能力是其他的多光谱影像数据所缺乏的。利用ASTER数据中的 3N和3B波段,可以构建地面DEM。生成的地面DEM可以用于遥感影像的正射校正,对冰川湖、冰川终点和雪线进行监测与制图。生成的DEM也可以用于与已有的地形图对比,以对快速的冰川变薄现象进行评估。

ASTER的3个15米分辨率的可见光/近红外波段和六个30米分辨率的短波红外波段已经被用于对冰、水的监测中;而其6个90米的热红外波段也开始被尝试用于冰川遥感监测工作。

ASTER影像在冰川监测上也存在一些缺点。其垂直分辨率较低,因而不能使用多时相ASTER数据生成的DEM对比以进行冰川表面变化变化监测。ASTER数据也缺乏穿透云雨的能力。

C、目前遥感冰川监测中存在的问题

当前遥感冰川监测主要存在两个问题。

一是对于岩屑覆盖的冰川的分类。这类冰川具有非线性的光谱特征,传统的基于光谱统计的多光谱影像分类方法对此无能为力。例如,在冰川地区影像分类中,对于污化雪、浊水、覆盖有岩屑的冰川和冰碛,目前的多光谱影像分类方法无法进行自动分类。在当前的工作中,对于这些类型地物通常采用人工目视解译分类的方法。

二是对于冰川变薄的监测。就目前应用而言,对于冰川的长度变化,利用多光谱影像能够进行较好的观察与监测。但是对于冰川物质平衡变化而言同样重要的冰川厚度变化,利用多光谱影像难以进行直接的监测。目前只能对于某些特定地区,在特定条件下可以使用多光谱影像,通过对冰川湖泊或冰川表面等的观测间接监测冰川厚度的变化。

(4)GLIMS概述

全球陆地冰空间监测组织(Global Land Ice Measurements from space:GLIMS)是由世界冰川监测组织(World Glacier Monitoring Service)、世界冰川目录(World Glacier Inventory)和美国国家地质调查局世界冰川卫星影像图集项目(USGS series Satellite Image Atlas of the Glacier s of the world)建立的。其目的在于帮助世界范围内的冰川研究人员对冰川变化进行记录。GLIMS现已有超过92名研究人员,62所研究机构,分布于25个国家。

GLIMS正在建立全球冰川数据库,数据库中存储冰川影像及一系列冰川属性信息。GLIMS研发提供了GLIMS view 等软件以辅助对冰川影像的人工或半自动的解译。对数据库的测试正在进行中,以便对数据的系统性、单一冰川数据内或多冰川数据间误差(地理位置误差、解译误差、由算法问题带来的误差或差异)进行分析与评估。

3、文章概况

该论文是一篇综述性的文章,主要从现有的技术进展概况,科学研究案例介绍,研究结果的实际作用与对气候变化研究的重要作用这些方面对于多光谱遥感数据特别是ASTER数据在全球陆地冰监测(GLIMS)的应用进行了介绍。

三、技术进展与实际研究案例
1、地形分类
(1)概述

DEM数据在可以用于多光谱影像的大气纠正,正射校正、获取冰川三维属性信息、对岩屑覆盖的冰川进行多维数据分类和其他的地貌学分析工作。

(2)ASTER数据在地形分析中的应用

ASTER数据的3N和3B波段可以做同轨道立体像对观测,生成相应的DEM。在特定良好情况下,利用生成的DEM可以对冰川的垂直厚度变化进行监测,并由此对冰川物质平衡变化进行评价与估计。但是在绝大多数情况下,由于ASTER影像生成的DEM数据在垂直分辨率上的不足,这样的分析并不能总是被进行。一般的说来ASTER数据更适用于对冰川的变化的长度、面积和表面状况进行监测。常用的DEM数据可以提取自于地形图、SRTM或ASTER数据。

引入DEM参与冰川地区分类也是提高冰川区域分类的重要手段。目前的研究中,研究者普遍认识到多光谱影像数据对于光谱特性呈非线性的岩屑覆盖冰川的无能为力。结合DEM数据及其派生的地形因子,利用邻域分析方法并引入热红外数据,在某些情况下可以提高冰川分类的精度。

ASTER影像RGB假彩色合成,叠加于由ASTER生成DTM之上

地点:LlewellynGlacier,BritishColumbia.Scaleis

(3)实例

以利用ASTER数据的极地冰川监测(原文3.3节)为例。ASTER数据在该研究中展现了其优点:

l 对南北纬86度内区域的覆盖并连续性的获得高分辨率的可见光波段的数据

l ASTER同轨道立体像对观测能力

l ASTER数据的观测计划意味这其数据是按需获取的。这样使得在极地研究中可以方便的设置年际或更频繁的数据获取计划。

在研究中,对于ASTER影像的正射校正工作利用了ASTER自身的立体观测所得到的数据,简化了工作。

2、冰川表面运动速度监测
(1)概述

冰川表面运动速度的监测通常是利用经过正射校正的高光谱数据(例如ASTER、ETM+/TM、MSS、SPOT)的相关分析而完成的。

(2)实例

以利用ASTER数据的极地冰川监测(原文3.3节)为例。在研究中,利用了ASTER数据按需获取的优势,得到了自2001年到2002年的影像数据。再对影像数据叠加分析,对相同地物自动追踪,以得到研究时间段内冰川表面运动速度的数据。

1978年经纬仪测得表面流动速度与2001ASTER影像得到的表面流动速度比较

地点:Byrd Glacier, Antarctica

以山地地貌与阿尔卑斯山地冰川分类研究(原文3.5节)为例。在研究中,利用ASTER数据对阿尔卑斯山地干流冰川与支流冰川的运动进行了监测,根据监测得到的支流冰川与干流冰川活动性之比(Ra),将冰川分为五类:干流冰川与支流并行(Ra=1);干流与之流处于分离中(Ra<<1);支流膨胀(Ra>1);支流越过干流(Ra>>1);干流与支流已分离。

3、综合历史地图与ASTER及其他影像数据
(1)概述

结合历史地图与遥感数据并进行定性与定量的比较,对变化进行分析对误差进行评估,这对于冰川遥感界是关键的也是经常进行的工作。

历史地图数据可以与遥感影像数据相结合进行分析。今年来的地图数据可以弥补遥感影像覆盖上的空缺;早期的地图数据可以提供在卫星遥感出现以前的历史数据。

如何有效的利用历史地图数据,特别是早期的纸质地图数据是研究中面临的重要问题。在特定情况下,纸质地图存储中产生的形变导致的误差使得数据难以甚至不能被用于定量的分析。另外在某些历史地形图中,地图数据本身在测绘学上是准确的,但其标注的冰川却并不完全正确,而是与积雪、海冰有所混淆。

除了历史地形图外,早期的航片和卫片也是遥感冰川变化分析的重要数据来源。

(2)实例

以对阿富汗和巴基斯坦的冰川变化研究为例(原文3.1节)。在研究中使用了历史地形图、航片、MSS、TM影像和ASTER影像对二十世纪70年代以来该地区的冰川变化进行了分析。

4、温度制图以发现冰湖变化
(1)概述

冰川的温度图像对于判断冰川表面状况非常有效。基于温度的物质区分方法包括:

a) 热惯量分析(对物体温度在一天之内进行多次观测,通常情况下不能由ASTER数据提供)

b) 不同日期同一白昼时刻温度的差异

c) 对水蒸发所造成的热焓的观测

利用冰川湖泊与周围温度的差异,可以将其与岩屑覆盖的冰川等其他类型地物区分出来。同一地区不同的冰川湖泊之间也存在的温度的差异,这些差异反映了不同湖泊冰水之间热交换的差异。

(2)实例

以对阿拉斯加的冰川研究(原文2.4节)为例。对于被岩屑覆盖了的冰川和位于冰川尾的堆满了冰山的冰湖,利用温度制图可以清楚的将二者区分出来,因为这些冰水湖的温度接近0摄氏度,而覆盖有岩屑的冰川其温度则相对较高。同时某些冰川湖泊相对于冰点有1-2k的温度的异常,这些异常也是可以被观测到的。另外,对于Martin River Glacier停滞的冰川末端的诸多冰川湖泊,观察各个湖泊的温度的不同以考察冰水间热交换的差异程度。

ASTER VNIRTIR波段对比图

地点:east of the CopperRiver, Chugach Range (Alaska).

以利用ASTER影像对秘鲁的冰川湖泊动态变化研究(原文3.6节)为例。由于季节融雪水进入冰川湖泊,而融雪水中所携带的物质的对于冰川湖泊的光谱特性的影响,使冰川湖泊的颜色和温度在不同季节有所不同。通过这些差异的观测可以对冰川表面状态、季节融雪情况进行观测;在实际应用上,对于防范冰湖溃决等冰川灾害也有很大的作用。

5、GLIMSView 与 GLIMS Database
(1)GLIMSView

GLIMSView是开放源代码的跨平台的软件,其开发目的在于辅助GLIMS项目的冰川数字化工作并为之提供相应的标准。在GLIMSView软件帮助下,用户可以将冰川数字化数据及获得的相关信息打包传送到国家雪冰数据中心,加入GLIMS数据库中以供今后的分析和进一步的信息提取。GLIMS允许用户浏览多种不同类型的卫星影像数据,数字化冰川轮廓与其他影像中的地物,在轮廓线段上添加GLIMS特定的属性信息。GLIMSView处理后的数据可以被保存为shape格式文件。

(2)GLIMS Database

GLIMS 冰川数据库是被设计在World Glacier Monitoring Service和World Glacier Inventory的基础上的。数据库中存储了多时态的类树状的关于一系列地物的信息。生成相关的信息的分析工作包括了对多种不同数据源自动的或人工的解译。数据库中存储了这些信息生成时所采用的分析方法和相关的分析人员的元数据。

国家雪冰数据中心(National Snow and Ice Data Center: NSIDC, Boulder, Colorado)设计并维护GLIMS Database。数据库的实现是基于开放源代码的地理信息数据库工具,包括PostgreSQL,PostGis和MapServer。数据库对外提供了符合OGC组织WMS和WFS标准的接口。

四、实际应用及对全球气候变化的反映
1、冰川灾害
       冰川与冻土灾害越来越威胁到高海拔地区人类的生活。相关的灾害将会造成大量的人员伤亡和财产损失。而随着冰冻圈对气候变化的动态反应和人类在高海拔的危险区的定居与其他活动的增加,冰川与冻土灾害的危险超过了以往任何时候。

       由于冰川区域大多难于接近,且由于冰川灾害评估与预警对于数据的快速获取的需求,遥感数据在高山地区灾害管理与灾害制图上有着极大的重要性。特别的,ASTER数据在冰川灾害管理与受灾情况评估上可以发挥重大的作用。一般情况下,ASTER的时间分辨率是15天,但是在紧急情况下,ASTER的VNIR传感器可以进行跨轨观测,其时间分辨率将缩短为2天内。这对于快速灾害情况评估和对危险状况监测与分析的工作有着特别的意义。

       下面列出一些最为重要和典型的冰川灾害以及ASTER数据对于这些冰川灾害可以做出的监测与评估。

l  冰川及冰川流动将会导致湖泊的形成,而这些湖泊的崩溃将产生洪水和岩屑流。利用ASTER的多光谱数据和DEM数据可以对冰川湖泊的形成进行监测并对未来的冰川湖泊的形成进行预测。

l  悬崖冰川的并崩裂将导致灾难性的冰崩,并且伴随有雪崩和岩石崩裂。ASTER的多光谱数据和DEM数据可以支持对这类灾难的监测与管理,,并对相关的可能导致这类灾难的因素进行建模。

l  冰川滑动可能导致下部基础的崩溃,冰川下水的喷出或是雪崩。利用ASTER数据重返率相对较高特点,使用时间上连续的ASTER影像,可以对冰川表面位移、冰川表面运动速度以及冰川表面岩屑分布模式的变化进行监测,从而对冰川滑动及可能伴随而来的灾害予以观测。

ASTER数据也可以用于对冰川湖泊进行监测(例原文3.6节)。利用ASTER数据可以对湖泊的边界和湖泊水色进行观测。冰川湖泊水色的观测,标志着季节性裂冰与融水流入湖中;湖泊边界的变化也标志着进入湖泊中的融水量的变化。但是ASTER数据的时间分辨率使得它还不足以被用于构建一个冰川湖泊灾害预警系统。MODIS数据在时间分辨率上足够,但却其空间分辨率使其无法被用于冰川湖泊灾害的预警。未来的可能的ASTER/ETM+的后继传感器或许能够被用于构建这样的预警系统。

2、全球气候变化与冰川变化
       现有研究已经揭示冰川对于全球气候变化有所响应,且主要表现为全球冰川退缩。但是全球冰川对气候变化的响应是非均一的,不同的冰川由于地区的差异、环境的不同、冰川自身属性的差异而有不同的响应,有的表现为冰川退缩,也有的冰川前进,抑或有冰川虽然其长度没有明显变化但其厚度变薄。

       由于冰川对气候变化响应模式的复杂性,仅仅凭现有的冰川观测数据难以对全球气候变化与冰川变化间的机理予以揭示也难以在二者间构建特定的联系。利用ASTER数据和其他遥感数据,对全球陆地冰变化进行监测与评估是GLIMS的目的。通过GLIMS获取的全球陆地冰川变化的完整的统一的监测数据,研究人员将能够对现有的GCM进行评价。

五、结论与讨论
1、ASTER数据在冰川变化监测中的应用
       由于其自身的特性(高分辨率,多光谱,高时间分辨率,具备同轨道立体像对观测能力),ASTER数据在冰川变化监测中已经有着广泛的应用,而且新的应用方法也正在被发展中。本篇论文中只是举出了部分实例。

       利用ASTER的同轨道立体像对观测生成的DEM可以方便数据的正射校正等工作。结合DEM和ASTER的多光谱波段数据,可以对冰川表面状态进行监测,从而间接的得到一些冰川的属性信息。在特定的条件下,可以对冰川的厚度等三维信息进行监测。

2、冰川对气候变化的响应
       从目前研究上看,全球冰川对气候变化有所响应,且整体上呈现退缩的态势。但是冰川对气候变化的响应是非均一的。冰川变化的影响因素非常多,机理也非常复杂并未为人们完全认识。凭借现有的对少数冰川的观测结果不能构建冰川变化与全球气候变化间的特定联系。GLIMS致力于对全球陆地冰川的变化进行监测和记录,这些数据将对气候变化和冰川响应的研究提供基础。

3、当前研究中存在的问题
       当前对于冰川变化监测与分析的遥感研究主要采用了多光谱数据。由于岩屑覆盖冰川、污化雪、混浊水之间光谱关系呈非线性,传统的基于统计的多光谱影像分类方法对之无能为力。这导致了冰川遥感难以自动化。另外,多光谱遥感数据对于冰川表面状态以及冰川厚度上也缺乏监测的手段。ASTER影像以及其他类似的遥感影像系统的引入将带来一定的改善。

六、本人观点
1.         当前最大的问题在于现有的遥感数据对于某些冰川分析工作并不适用。例如,在冰川遥感中,确定冰川轮廓和冰川边界是最重要的工作。这就不可避免的涉及到了对岩屑覆盖冰川、污化雪和冰碛的分类。目前的多光谱遥感数据在传统的基于光谱统计的分类方法下对这几类地物难以区分。这直接导致了现有的多光谱冰川遥感无法自动化。在文中作者反复提及这一点,并引述了诸多文献,而我本人在一些相关的工作中也有感受。这一点使得冰川遥感的工作效率无法提高。今后的研究中,是否有可能引入一些其他类型的,对于岩屑覆盖冰川、污化雪和冰碛有较大差异的数据?例如ASTER的热红外波段,或是主动微波遥感数据。或是将人的地学知识引入冰川遥感分类工作中,实现对冰川的自动分类。

2.         综合DEM数据参与冰川遥感分析。随着ASTER和SRTM数据使用逐渐普及,DEM数据在冰川遥感分析工作中也将得到更广泛的利用。结合DEM和多光谱数据可以对冰川表面情况予以分析。同时,通过高精度的DEM数据也能够得到一些原先只能通过实地考察得到的冰川三维数据。GLIMSView软件已经集成提供了一些功能,方便用户在冰川遥感解译工作中集成DEM和多光谱数据,为冰川遥感半自动解译提供了辅助。

3.         结合历史数据是冰川遥感分析中的常见工作。依靠遥感数据与历史数据的对比或时序分析,可以获得冰川进退变化信息。如何评判衡量历史地形图数据(特别是较早的纸质地形图数据)的误差是一个重要的问题。另外,现有大量存档的早期航空胶片和卫星胶片如何利用也是一个重要的问题。

4.         空间数据库与空间数据分析在冰川遥感工作中的应用。对于全球范围的陆地冰川进行监测将生产大量的数据。GLIMS已经采用了开源的软件集构建空间数据以存储这些数据。同时GLIMS在数据库中引入了数据标准,使得未来对数据的空间分析成为可能;而GLIMS中的记录了数据生产方法与分析人员的元数据也将方便未来的分析工作。这些对于现在的冰川遥感工作和冰川目录编目工作都是一个很大的进步。未来如何利用好空间数据库存储的数据,加以空间分析也是重要的问题。

  评论这张
 
阅读(903)| 评论(0)
推荐 转载

历史上的今天

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2018