森林资源调查中SPOT5遥感图像处理方法探讨王照利、黄生、张敏中、马胜利(国家林业局西北林业规划设计院,遥感计算中心,西安710048)本文发表于<陕西林业科技>2005 摘要: 目前,多光谱、高空间分辨率的SPOT5卫星遥感数据被广泛应用到森林资源调查中。本文结合SPOT5遥感数据的特点,根据森林资源调查的需要,从遥感数据的正射校正、波段组合、融合处理和数据变换处理等方面探讨了SPOT5数据的处理和信息提取。探讨性地提出了适应于森林资源调查的SPOT5遥感数据处理方法。 关键词:SPOT5 遥感数据,森林资源调查、数据处理DISCUSSION ON SPOT5 IMAGE DATA PROCESSING FOR FOREST INVENTORYWang Zhaoli, Huangsheng,Zhangminzhong,Ma Shengli(Northwest Institute for Forest Inventory, Planning &Design, Xi’an China 710048)Abstract: Now days, high spatial resolution and multispectral SPOT5 image data are widely applied in forest inventory in China. Based on the characteristics of SPOT5 image and requirements of forest inventory, this paper discusses the processing procedures of ordering image data, ortho-rectification, image bands composition and image data fusion. The complete steps of image processing for forest inventory are words: SPOT5 image data,forest inventory, data processing 前言 卫星遥感影像具有空间宏观性、视角广、多分辨率(光谱和空间)、多时相、周期性、信息量丰富等特点,所以卫星遥感影像既可以提供森林资源的宏观空间分布信息又能提供局部的详细信息以及随时间、空间变化的信息等[1]。目前在林业领域卫星遥感数据被广泛的应用于不同尺度层次的森林资源调查、资源监测、病虫害、火灾监测等方面。2002年5月法国SPOT地球观测卫星系列之5号卫星(即SPOT5星)发射。SPOT5遥感数据的多光谱波段空间分辨率为10米(短波红外空间分辨率为20米),但全色波段空间分辨率达到米。SPOT5遥感数据的高空间分辨率和多光谱分辨率为森林资源调查提供了丰富的、可靠的、高精度的基础数据源。从性价比分析,在其他高分辨率遥感数据目前比较昂贵的状况下,SPOT5遥感数据比较适宜应用于大面积的森林资源调查,可大幅度的森林调查的减少外业工作量、提高工作效率。在我国SPOT5卫星数据已被大量地应用于森林资源调查工作中,尤其,是在森林资源“二类”调查中被作基本的森林资源信息源提取各类信息。针对于将多光谱分辨率和高空间分辨率的SPOT5遥感数据应用于森林资源调查的数据处理技术和方法鲜有报道。本文总结工作实践,结合SPOT5遥感数据的特点,根据森林资源调查的需要,从遥感数据的订购、正射校正、波段组合、融合处理和数据变换处理等方面探讨了SPOT5数据的基本处理方法。 1.SPOT5卫星遥感数据特点 SPOT卫星系统采用线性阵列传感器和推扫式扫描技术,具有旋转式平面镜可以进行倾斜观察获得倾斜图像和立体像对。采用与太阳同步的近极地的椭圆形轨道,轨道高度约832Km,轨道倾角 ,每天绕地球14圈多,重复覆盖周期26天[2]。由于有倾斜观测功能,使重复覆盖周期减少到2-3天。SPOT5卫星载有2台高分辨率几何成像仪(HRG)、1台高分辨率立体成像装置(HRS)和1台宽视域植被探测仪(VGT)。高分辨率几何成像仪的波段选择是总结了多年的研究成果,认为HRG的波段设置(见表1)足以取得辨别作物和植被类型的最佳效果。本文主要探讨HRG高空间分辨率数据的处理。2.SPOT5数据的处理方法和过程 SPOT5数据处理工作流程: 遥感数据的订购 订购数据时,用户需向数据代理商提供购买区域的四个角的大地坐标或者数据的景号(PATH/ROW)。特别应该注意数据订购时间和用户拿到数据之间有时间差,间隔时间长短因用户的要求、天气、卫星重复覆盖周期而异。相对于其他卫星数据,比较有利的一面是SPOT5卫星装置有旋转式平面镜可以进行倾斜观察,用户可向代理商申请红色编程提前得到调查区域的遥感数据,但要支付编程费。对于遥感数据的时相、云量、入射角、阴影量、是否购买高空间分辨率的全色波段等用户根据自己具体的工作需要向代理商提出限制要求。 根据我们对SPOT5遥感数据的使用,对于森林资源调查,北方9,10月份和11月初的遥感影像比较适宜。代理商向用户提供经过处理的不同级别的影像产品,在森林资源调查中建议购买SPOT1A级产品,用户可根据自己的工作需要进行处理,同时也可减少费用。 基础数据准备 大比例尺地形图和高精度DEM是进行SPOT5遥感数据高精度正射校正必需的基础地理数据。建议购买1:10000地形图和1:25000数字高程模型(DEM)。 将1:1万地形图扫描,扫描分辨率设置为300DPI。将扫描好的地形图进行几何精纠正,纠正精度控制在毫米内。从测绘部门购买的1:1万地形图为北京54坐标系3度分带高斯克吕格投影,而1:万DEM为北京54坐标系6度分带投影。在数据准备时,将校正好的1:1万地形图通过换带转换转成和DEM一致的6度分带投影。 对于没有1:1万地形图的地区,建议使用差分GPS接收机采集地面控制点。 几何正射校正 正射校正过程应用了法国SPOT公司发行的GEOIMAGE软件。GEOIMAGE软件有针对SPOT5卫星数据开发的SPOT5物理模型。模型模块自动读取DEM信息。SPOT 物理模型可读取卫星在获取遥感数据的瞬间状态参数,这些参数存贮在数据的头文件中[3]。卫星状态参数包括:卫星成像瞬间的经纬度、高度、倾角等。卫星状态参数能够帮助提高几何校正的精度。 以校正好的1:1万地形图为基准,在影像图上找出和地形图上地物相匹配的明显地物作为地面控制点。在进行正射校正时,应先进行全色波段数据校正,然后以校正好的全色波段数据为基准进行多光谱数据校正。以全色波段数据为基准校正多光谱波段就比较容易校正,且能提高两者的匹配精度。地面控制点应分布均匀,影像的边缘部分布要有控制点分布,同时在不同的高程范围最好都有控制点。地面控制点的数量因地形地貌的复杂程度而定,根据我们的经验,一景60KmX60Km的SPOT5数据,一般地势平缓的地区20个左右控制点即可达到满意的结果,在高山区25个左右控制点就可使正射校正精度满足要求。重采样方法采用双线性内插法。 辐射校正 用户购买的SPOT5的各级数据,数据提供商已经根据卫星的记录参数对遥感数据做了辐射校正,即消除了传感器自身引起的、大气辐射引起的辐射噪声。若果影像存在薄雾或地形高差较大引起的辐射误差情况,用户应进一步进行辐射校正处理。薄雾的简单消除原理是基于近红外波段不受大气辐射影响,清澈的水体或死阴影区的数值应为零。从各波段数据中减去近红外波段的水体或阴影的不为零值。地形起伏引起的辐射误差校正公式: f (x,y)=g(x,y)/cosa,g(x,y)为坡度为a的倾斜面上的地物影像;f (x,y)为校正后的影像。由于坡度因子参与校正所以需要DEM支持。 波段组合 根据SPOT5数据波谱特征(表1),各波段分别记录反映了植被的不同特征方面:B4(SWIR)短波红外反映植物和土壤的含水量,利于植被水分状况和长势分析;B3(NIR)近红外波段对植被类别、密度、生长力、病虫害等的变化敏感;B2(RED)红光波段对植被的覆盖度、植被的生长状况敏感;B1(VIS)可见光波段对植物的叶绿素和叶绿素浓度敏感。经过比较分析和实际应用发现SPOT5的B3、B4、B2波段组合对植被类型的识别要优于B3、B2和B1的组合。但由于B4波段的空间分辨率为20米,使B342组合对植被空间几何细节表达没有B321组合清晰,例如林缘界线信息表达方面B321要优于B342。 影像数据融合 对于购买有高空间分辨率全色波段数据的用户,进行数据融合是必不可少的。影像数据融合能够综合不同波段、不同空间分辨率数据(层)的特征,融合后的数据具有更丰富、更可靠的信息[4]。 根据影像数据融合的水平阶段,影像融合分为:像元级、特征级和决策级三个层次。为了最大限度的从SPOT5遥感数据中提取森林植被信息,应进行像元级的数据融合,将米的全色波段和10米多光谱数据进行融合。融合得到的新数据既具有全色波段数据的高空间分辨率特征又具有多光谱特征。像元级数据融合的方法多种多样,根据数据融合的目的,即最大限度的突显森林植被信息,应选取B4、B3、B2和PAN波段,根据我们的试验Brovey 融合算法方法比较理想:遥感影像地图 将融合好的数据按Rfused、Gfused、Bfused组合,叠加上行政界线、公里格网、坐标、比例尺等辅助信息,按1:1万地形图分幅生成1:1万纸质图作为外业手图。 3. 结果和讨论 几何精度 利用SPOT5物理模型,采用1:1万地形图和万DEM ,经过正射校正处理,可使影像的几何精度控制在2个像元内(<10米),达到1:1万制图标准要求。为以遥感影像为基础信息源提取林分调查因子、区划林班界线生成大比例尺的林相图、森林分布图提供了几何精度保障。 波段选择 对于没有全色波段的情况,SPOT5数据的B342组合有利于森林植被类型的识别。在应用遥感技术进行森林资源调查区划中,林分类型信息提取是最为重要的环节,所以B342波段组合是小班区划和外业手图的最佳组合。 融合效果 融合数据技术使SPOT5遥感影像既具有全色波段的高空间分辨率又拥有多光谱数据的光谱分辨率,丰富了遥感影像的信息量。采用Brovey算法使SPOT5遥感影像从色彩、纹理等方面增强了影像的可判读性,提高了小班因子正判率和林分小班的区划精度。 参考文献 1.周成虎,杨晓梅,骆剑承等.《遥感影像地学理解与分析》,科学出版社,北京,2001,3-4. 2.赵英时.《遥感应用分析原理与方法》,科学出版社,北京, 3.北京视宝卫星图像有限公司.《专业制图工作室GEOIMAGE用户指南》,2004,68-70. 4.Christine Pohl. Geometric Aspects of Multisensor Image Fusion for Topographic Map Updating in The Humid Tropics, ITC Publication, 1996,世纪遥感与GIS的发展 来源: 李德仁 时间: 2005-08-11-23:09 浏览次数: 79 21世纪遥感与GIS的发展李德仁 (武汉大学测绘遥感信息工程国家重点实验室,武汉市珞瑜路129号,430079)摘要:在20世纪,人类的一大进步是实现了太空对地观测,即可以从空中和太空对人类赖以生存的地球通过非接触传感器的遥感进行观测,并将所得到的数据和信息存储在计算机网络上,为人类社会的可持续发展服务。在短短的30年中,遥感和GIS作为一个边缘交叉学科已发展成为一门科学、技术和经济实体。本文深入地论述了21世纪中遥感的6大发展趋势和GIS的5个发展特征。 关键词:发展趋势;航空航天遥感;地理信息系统;对地观测 中图法分类号:P208; 随着计算机技术、空间技术和信息技术的发展,人类实现了从空中和太空来观测和感知人类赖以生存的地球的理想,并能将所感知到的结果通过计算机网络在全球流通,为人类的生存、繁荣和可持续发展服务。在20世纪后半叶,遥感和地理信息系统作为一门新兴的科学和技术,迅速地成长起来。 1 遥感技术的主要发展趋势 航空航天遥感传感器数据获取技术趋向三多(多平台、多传感器、多角度)和三高(高空间分辨率、高光谱分辨率和高时相分辨率) 从空中和太空观测地球获取影像是20世纪的重大成果之一,短短几十年,遥感数据获取手段迅猛发展。遥感平台有地球同步轨道卫星(35000km)、太阳同步卫星(600—1000km)、太空飞船(200—300km)、航天飞机(240—350km)、探空火箭(200—1000km),并且还有高、中、低空飞机、升空气球、无人飞机等;传感器有框幅式光学相机、缝隙、全景相机、光机扫描仪、光电扫描仪、CCD线阵、面阵扫描仪、微波散射计雷达测高仪、激光扫描仪和合成孔径雷达等,它们几乎覆盖了可透过大气窗口的所有电磁波段。三行CCD阵列可以同时得到3个角度的扫描成像,EOS Terra卫星上的MISR可同时从9个角度对地成像。 卫星遥感的空间分辨率从Ikonos Ⅱ的1m,进一步提高到Quckbird(快鸟)的,高光谱分辨率已达到5—6nm,500—600个波段。在轨的美国EO-1高光谱遥感卫星,具有220个波段,EOS AM-1(Terra)和EOS PM-1(Aqua)卫星上的MODIS具有36个波段的中等分辨率成像光谱仪。时间分辨率的提高主要依赖于小卫星技术的发展,通过发射地球同步轨道卫星和合理分布的小卫星星座,以及传感器的大角度倾斜,可以以1—3d的周期获得感兴趣地区的遥感影像。由于具有全天候、全天时的特点,以及用INSAR和D-INSAR,特别是双天线INSAR进行高精度三位地形及其变化测定的可能性,SAR雷达卫星为全世界各国所普遍关注。例如,美国宇航局的长远计划是要发射一系列太阳同步和地球同步的长波SAR,美国国防部则要发射一系列短波SAR,实现干涉重访问间隔为8d、3d和1d,空间分辨率分别为20m、5m和2m。我国在机载和星载SAR传感器及其应用研究方面正在形成体系。“十五”期间,我国将全方位地推进遥感数据获取的手段,形成自主的高分辨率资源卫星、雷达卫星、测图卫星和对环境与灾害进行实时监测的小卫星群。 航空航天遥感对地定位趋向于不依赖地面控制 确定影像目标的实地位置(三维坐标),解决影像目标在哪儿(Where)是摄影测量与遥感的主要任务之一。在已成功用于生产的全自动化GPS空中三角测量的基础上,利用DGPS和INS惯性导航系统的组合,可形成航空/航天影像传感器的位置与姿态的自动测量和稳定装置(POS),从而可实现定点摄影成像和无地面控制的高精度对地直接定位。在航空摄影条件下的精度可达到dm级,在卫星遥感的条件下,其精度可达到m级。该技术的推广应用,将改变目前摄影测量和遥感的作业流程,从而实现实时测图和实时数据库更新。若与高精度激光扫描仪集成,可实现实时三维测量(LIDAR),自动生成数字表面模型(DSM),并可推算出数字高程模型(DEM)。 美国NASA在1994年和1997年两次将航天激光测高仪(SLA)安装在航天飞机上,企图建立基于SLA的全球控制点数据库,激光点大小为100m,间隔为750m,每秒10个脉冲;随后又提出了地学激光测高系统(GLAS)计划,已于2002年12月19日将该卫星IICESat(cloud and land elevation satellite)发射上天。该卫星装有激光测距系统、GPS接收机和恒星跟踪姿态测定系统。GLAS发射近红外光(1064nm)和可见绿光(532nm)的短脉冲(4ns)。激光脉冲频率为40次/s,激光点大小实地为70m,间隔为170m,其高程精度要明显高于SRTM,可望达到m级。他们的下一步计划是要在2015年之前使星载LIDAR的激光测高精度达到dm和cm级。 法国利用设在全球的54个站点向卫星发射信号,通过测定多普勒频移,以精确解求卫星的空间坐标,具有极高的精度。测定距地球1300km的Topex/Poseidon卫星的高度,精度达到±3cm。用来测定SPOT 4卫星的轨道,3个坐标方向达到±5cm精度,对于SPOT 5和Envisat,可望达到±1m精度。若忽略SPOT 5传感器的角元素,直接进行无地面控制的正射像片制作,精度可达到±15m,完全可以满足国家安全和西部开发的需求。 摄影测量与遥感数据的计算机处理更趋向自动化和智能化 从影像数据中自动提取地物目标,解决它的属性和语义(What)是摄影测量与遥感的另一大任务。在已取得影像匹配成果的基础上,影像目标的自动识别技术主要集中在影像融合技术,基于统计和基于结构的目标识别与分类,处理的对象既包括高分辨率影像,也更加注重高光谱影像。随着遥感数据量的增大,数据融合和信息融合技术逐渐成熟。压缩倍率高、速度快的影像数据压缩方法也已商业化。我国学者在这些方面取得了不少可喜的成果。 利用多时像影像数据自动发现地表覆盖的变化趋向实时化 利用遥感影像自动进行变化监测(What change)关系到我国的经济建设和国防建设。过去人工方法投入大,周期长。随着各类空间数据库的建立和大量新的影像数据源的出现,实时自动化监测已成为研究的一个热点。 自动变化监测研究包括利用新旧影像(DOM)的对比、新影像与旧数字地图(DLS)的对比来自动发现变化和更新数据库。目前的变化监测是先将新影像与旧影像(或数字地图)进行配准,然后再提取变化目标,这在精度、速度与自动化处理方面都有不足之处。笔者提出了把配准与变化监测同步的整体处理[1]。最理想的方法是将影像目标三维重建与变化监测一起进行,实现三维变化监测和自动更新。进一步的发展则是利用智能传感器,将数据处理在轨完成,发送回来的直接为信息,而不一定为影像数据。 摄影测量与遥感在构建“数字地球”、“数字中国”、“数字省市”和“数字文化遗产”中正在发挥愈来愈大的作用 “数字地球”概念是在全球信息化浪潮推进下形成的。1999年12月在北京成功地召开了第一届国际“数字地球”大会后,我国正积极推进“数字中国”和“数字省市”的建设,2001年国家测绘局完成了构建“数字中国”地理空间基础框架的总体战略研究。在已完成1∶100万和1∶25万全国空间数据库的基础上,2001年全国各省市测绘局开始1∶5万空间数据库的建库工作。在这个数据量达11TB的巨型数据库中,摄影测量与遥感将用来建设DOM(数字正射影像)、DEM(数字高程模型)、DLG(数字线划图)和CP(控制点数据库)。如果要建立全国1m分辨率影像数据库,其数据量将达到60TB。如果整个“数字地球”均达到1m分辨率,其数据量之大可想而知。本世纪内可望建成这一分辨率的数字地球。 “数字文化遗产”是目前联合国和许多国家关心的一个问题,涉及到近景成像、计算机视觉和虚拟现实技术。在近景成像和近景三位量测方面,有室内各种三维激光扫描与成像仪器,还可以直接由视频摄像机的系列图像获取目标场三维重建信息。它们所获取的数据经过计算机自动处理后,可以在虚拟现实技术支持下形成文化遗迹的三维仿真,而且可以按照时间序列,将历史文化在时间隧道中再现,对文化遗产保护、复原与研究具有重要意义。 全定量化遥感方法将走向实用 从遥感科学的本质讲,通过对地球表层(包括岩石圈、水圈、大气圈和生物圈4大圈层)的遥感,其目的是为了获得有关地物目标的几何与物理特性,所以需要通过全定量化遥感方法进行反演。几何方程式是有显式表示的数学方程,而物理方程一直是隐式。目前的遥感解译与目标识别并没有通过物理方程反演,而是采用了基于灰度或加上一定知识的统计、结构和纹理的影像分析方法。但随着对成像机理、地物波谱反射特征、大气模型、气溶胶的研究深入和数据积累,多角度、多传感器、高光谱及雷达卫星遥感技术的成熟,相信在21世纪,估计几何与物理方程式的全定量化遥感方法将逐步由理论研究走向实用化,遥感基础理论研究将迈上新的台阶。只有实现了遥感定量化,才可能真正实现自动化和实时化。 2 GIS技术的主要发展趋势 空间数据库趋向图形、影像和DEM三库一体化和面向对象[2] GIS发展曾经历过栅格、矢量两个不同数据结构发展阶段,目前随着高分辨率卫星遥感数据的飞快增长和数字地球、数码城市的需求,形成了面向对象的数据模型和三库(图形矢量库、影像栅格库和DEM格网库)一体化的数据结构。这样的数据库结构使GIS的发展更加趋向自然化、逼真化,更加贴近用户。以面向应用的GIS软件为前台,以大型关系数据库(Oracle 8i,9i等)为后台数据库管理,成为当前GIS技术的主流趋势。 空间数据表达趋向多比例尺、多尺度、动态多位和实时三维可视化 在传统的GIS中,空间数据是以二维形式存储并挂接相应的属性数据。目前,空间数据表达的趋势是基于金字塔和LOD(level of detail)技术的多比例尺空间数据库,在不同尺度表示时可自动显示出相应比例尺或相应分辨率的数据,多比例尺数据集的跨度要比传统地图的比例尺大,在显示不同比例尺数据时,可采用LOD或地图综合技术。真三维GIS的空间数据要存储三维坐标。动态GIS在土地变更调查、土地覆盖变化监测中已有较好的应用,真四维的时空GIS将有望从理论研究转入实用阶段。基于三库一体化的时空3D可视化技术发展势头迅猛,已能再PC机上实现GIS环境下的三维建筑物室外室内漫游、信息查询、空间分析、剖面分析和阴影分析等,基于虚拟现实技术的真三维GIS将使人们在现实空间外,可以同时拥有一个Cyber空间。 空间分析和辅助决策智能化需要利用数据挖掘方法从空间数据库和属性数据库中发现更多的有用知识 GIS是以应用导向的空间信息技术,空间分析与辅助决策支持是GIS的高水平应用,它需要基于知识的智能系统。知识的获取是专家系统中最困难的任务。随着各种类型数据库的建立,从数据库中挖掘知识成为当今计算机界一个非常引人注目的课题。从GIS空间数据库中发现的知识可以有效的支持遥感图像解译,以解决“同物异谱”和“同谱异物”的问题。反过来,从属性数据库中挖掘的知识又具有优化资源配置等一些列空间分析的功能[3]。尽管数据挖掘和知识发现这一命题仍处于理论研究阶段,但随着数据库的快速增大和对数据挖掘工具的深入研究,其应用前景是不可估量的。 通过Web服务器和WAP服务器的互联网和移动GIS将推进联邦数据库和互操作的研究及地学信息服务事业 随着计算机通讯网络(包括有线和无线网)的大容量和高速化,GIS已成为在网络上的分布式异构系统。许多不同单位、不同组织维护管理的既独立又互联互用的联邦数据库,将可提供全社会各行各业的应用需要。因此,联邦数据库和互操作(federal databases & interoperability)问题成为当前国际GIS联合研究的一个热点。互操作意味着数据库中数据的直接共享,GIS规律功能模块的互操作与共享,以及多点之间的相同工作,这方面的研究已显示出明显的成效。未来的GIS用户将可能在网络上缴纳为其需要所选用数据和软件功能模块的使用费,而不必购买这个数据库和整套的GIS软硬件,这些成果产生的直接效果是GIS应用将走向地学信息服务。 目前已兴起的LBS和MLS,即基于位置的服务和移动定位服务,突出地反映了这种变化趋势。它引起的革命性变化使GIS将走出研究院所和政府机关,成为全社会人人具备的信息服务工具。我国目前已有2亿个手机用户,若每人每月为MLS支付10元费用,全国一年的产值将达到240亿。可以预测在不久的将来,地学信息将能随时随地为任何人和任何事情进行4A服务(geo-in-formation for anyone and anything at anywhere and anytime)。 地理信息科学的研究有望在本世纪形成较完整的理论框架体系 笔者曾扼要地叙述了地球空间信息科学的7大理论问题[4]:(1)地球空间信息的基准,包括几何基准、物理基准和时间基准;(2)地球空间信息标准,包括空间数据采集、存储与交换标准、空间数据精度与质量标准、空间信息的分类与代码标准、空间信息的安全、保密及技术服务标准以及元数据标准等;(3)地球空间信息的时空变化理论,包括时空变化发现的方法和对时空变化特征的和规律的研究;(4)地球空间信息的认知,主要通过各目标各要素的位置、结构形态、相互关联等从静态上的形态分析、发生上的成因分析、动态上的过程分析、演化上的力学分析以及时态上的演化分析达到对地球空间的客观认知;(5)地球空间信息的不确定性,包括类型的不确定性、空间位置的不确定性、空间关系的不确定性、逻辑的不一致性和信息的不完备性;(6)地球空间信息的解译与反演,包括定性解译和定量反演,贯穿在信息获取、信息处理和认知过程中;(7)地球空间信息的表达与可视化,涉及到空间数据库多分辨率表示、数字地图自动综合、图形可视化、动态仿真和虚拟现实等。目前,这些方面的研究对建立完备的理论尚嫌不足,需要在今后加强这方面的基础研究。 关于遥感与GIS的集成,涉及到GPS和通信技术的集成,本文未作具体讨论,其具体内容可参见文献[4—6]。 3 结语 遥感与GIS在20世纪出现,在21世纪不仅将形成自身的理论体系和技术体系,而且将形成天地一体化的空间信息服务产业,为国民经济建设、国家安全、社会可持续发展和提高人民生活质量做出愈来愈大的贡献。 参考文献: [1] Li D R, Sui H G. Automatic Change Detection of Geospatial Data from Imagery. The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, 2002,34(II):245—251 [2] 龚健雅. 地理信息系统基础. 北京:科学出版社,2001 [3] 邸凯昌. 空间数据发掘与知识发现(第一版). 武汉:武汉大学出版社,2000. 182 [4] 李德仁,关泽群. 空间信息系统的集成与实现(第一版). 武汉:武汉测绘科技大学出版社,2000. 244 [5] 李德仁,李清泉. 论地球空间信息技术与通信技术的集成. 武汉大学学报(信息科学版),2001,26(1):1—7 [6] 李德
森林资源调查中SPOT5遥感图像处理方法探讨王照利、黄生、张敏中、马胜利(国家林业局西北林业规划设计院,遥感计算中心,西安710048)本文发表于<陕西林业科技>2005 摘要:目前,多光谱、高空间分辨率的SPOT5卫星遥感数据被广泛应用到森林资源调查中。本文结合SPOT5遥感数据的特点,根据森林资源调查的需要,从遥感数据的正射校正、波段组合、融合处理和数据变换处理等方面探讨了SPOT5数据的处理和信息提取。探讨性地提出了适应于森林资源调查的SPOT5遥感数据处理方法。关键词:SPOT5 遥感数据,森林资源调查、数据处理DISCUSSION ON SPOT5 IMAGE DATA PROCESSING FOR FOREST INVENTORYWang Zhaoli, Huangsheng,Zhangminzhong,Ma Shengli(Northwest Institute for Forest Inventory, Planning &Design, Xi’an China 710048)Abstract: Now days, high spatial resolution and multispectral SPOT5 image data are widely applied in forest inventory in China. Based on the characteristics of SPOT5 image and requirements of forest inventory, this paper discusses the processing procedures of ordering image data, ortho-rectification, image bands composition and image data fusion. The complete steps of image processing for forest inventory are words: SPOT5 image data,forest inventory, data processing前言卫星遥感影像具有空间宏观性、视角广、多分辨率(光谱和空间)、多时相、周期性、信息量丰富等特点,所以卫星遥感影像既可以提供森林资源的宏观空间分布信息又能提供局部的详细信息以及随时间、空间变化的信息等[1]。目前在林业领域卫星遥感数据被广泛的应用于不同尺度层次的森林资源调查、资源监测、病虫害、火灾监测等方面。2002年5月法国SPOT地球观测卫星系列之5号卫星(即SPOT5星)发射。SPOT5遥感数据的多光谱波段空间分辨率为10米(短波红外空间分辨率为20米),但全色波段空间分辨率达到米。SPOT5遥感数据的高空间分辨率和多光谱分辨率为森林资源调查提供了丰富的、可靠的、高精度的基础数据源。从性价比分析,在其他高分辨率遥感数据目前比较昂贵的状况下,SPOT5遥感数据比较适宜应用于大面积的森林资源调查,可大幅度的森林调查的减少外业工作量、提高工作效率。在我国SPOT5卫星数据已被大量地应用于森林资源调查工作中,尤其,是在森林资源“二类”调查中被作基本的森林资源信息源提取各类信息。针对于将多光谱分辨率和高空间分辨率的SPOT5遥感数据应用于森林资源调查的数据处理技术和方法鲜有报道。本文总结工作实践,结合SPOT5遥感数据的特点,根据森林资源调查的需要,从遥感数据的订购、正射校正、波段组合、融合处理和数据变换处理等方面探讨了SPOT5数据的基本处理方法。1.SPOT5卫星遥感数据特点SPOT卫星系统采用线性阵列传感器和推扫式扫描技术,具有旋转式平面镜可以进行倾斜观察获得倾斜图像和立体像对。采用与太阳同步的近极地的椭圆形轨道,轨道高度约832Km,轨道倾角 ,每天绕地球14圈多,重复覆盖周期26天[2]。由于有倾斜观测功能,使重复覆盖周期减少到2-3天。SPOT5卫星载有2台高分辨率几何成像仪(HRG)、1台高分辨率立体成像装置(HRS)和1台宽视域植被探测仪(VGT)。高分辨率几何成像仪的波段选择是总结了多年的研究成果,认为HRG的波段设置(见表1)足以取得辨别作物和植被类型的最佳效果。本文主要探讨HRG高空间分辨率数据的处理。2.SPOT5数据的处理方法和过程SPOT5数据处理工作流程: 遥感数据的订购订购数据时,用户需向数据代理商提供购买区域的四个角的大地坐标或者数据的景号(PATH/ROW)。特别应该注意数据订购时间和用户拿到数据之间有时间差,间隔时间长短因用户的要求、天气、卫星重复覆盖周期而异。相对于其他卫星数据,比较有利的一面是SPOT5卫星装置有旋转式平面镜可以进行倾斜观察,用户可向代理商申请红色编程提前得到调查区域的遥感数据,但要支付编程费。对于遥感数据的时相、云量、入射角、阴影量、是否购买高空间分辨率的全色波段等用户根据自己具体的工作需要向代理商提出限制要求。根据我们对SPOT5遥感数据的使用,对于森林资源调查,北方9,10月份和11月初的遥感影像比较适宜。代理商向用户提供经过处理的不同级别的影像产品,在森林资源调查中建议购买SPOT1A级产品,用户可根据自己的工作需要进行处理,同时也可减少费用。 基础数据准备大比例尺地形图和高精度DEM是进行SPOT5遥感数据高精度正射校正必需的基础地理数据。建议购买1:10000地形图和1:25000数字高程模型(DEM)。将1:1万地形图扫描,扫描分辨率设置为300DPI。将扫描好的地形图进行几何精纠正,纠正精度控制在毫米内。从测绘部门购买的1:1万地形图为北京54坐标系3度分带高斯克吕格投影,而1:万DEM为北京54坐标系6度分带投影。在数据准备时,将校正好的1:1万地形图通过换带转换转成和DEM一致的6度分带投影。对于没有1:1万地形图的地区,建议使用差分GPS接收机采集地面控制点。几何正射校正正射校正过程应用了法国SPOT公司发行的GEOIMAGE软件。GEOIMAGE软件有针对SPOT5卫星数据开发的SPOT5物理模型。模型模块自动读取DEM信息。SPOT 物理模型可读取卫星在获取遥感数据的瞬间状态参数,这些参数存贮在数据的头文件中[3]。卫星状态参数包括:卫星成像瞬间的经纬度、高度、倾角等。卫星状态参数能够帮助提高几何校正的精度。以校正好的1:1万地形图为基准,在影像图上找出和地形图上地物相匹配的明显地物作为地面控制点。在进行正射校正时,应先进行全色波段数据校正,然后以校正好的全色波段数据为基准进行多光谱数据校正。以全色波段数据为基准校正多光谱波段就比较容易校正,且能提高两者的匹配精度。地面控制点应分布均匀,影像的边缘部分布要有控制点分布,同时在不同的高程范围最好都有控制点。地面控制点的数量因地形地貌的复杂程度而定,根据我们的经验,一景60KmX60Km的SPOT5数据,一般地势平缓的地区20个左右控制点即可达到满意的结果,在高山区25个左右控制点就可使正射校正精度满足要求。重采样方法采用双线性内插法。 辐射校正用户购买的SPOT5的各级数据,数据提供商已经根据卫星的记录参数对遥感数据做了辐射校正,即消除了传感器自身引起的、大气辐射引起的辐射噪声。若果影像存在薄雾或地形高差较大引起的辐射误差情况,用户应进一步进行辐射校正处理。薄雾的简单消除原理是基于近红外波段不受大气辐射影响,清澈的水体或死阴影区的数值应为零。从各波段数据中减去近红外波段的水体或阴影的不为零值。地形起伏引起的辐射误差校正公式: f (x,y)=g(x,y)/cosa,g(x,y)为坡度为a的倾斜面上的地物影像;f (x,y)为校正后的影像。由于坡度因子参与校正所以需要DEM支持。 波段组合根据SPOT5数据波谱特征(表1),各波段分别记录反映了植被的不同特征方面:B4(SWIR)短波红外反映植物和土壤的含水量,利于植被水分状况和长势分析;B3(NIR)近红外波段对植被类别、密度、生长力、病虫害等的变化敏感;B2(RED)红光波段对植被的覆盖度、植被的生长状况敏感;B1(VIS)可见光波段对植物的叶绿素和叶绿素浓度敏感。经过比较分析和实际应用发现SPOT5的B3、B4、B2波段组合对植被类型的识别要优于B3、B2和B1的组合。但由于B4波段的空间分辨率为20米,使B342组合对植被空间几何细节表达没有B321组合清晰,例如林缘界线信息表达方面B321要优于B342。 影像数据融合对于购买有高空间分辨率全色波段数据的用户,进行数据融合是必不可少的。影像数据融合能够综合不同波段、不同空间分辨率数据(层)的特征,融合后的数据具有更丰富、更可靠的信息[4]。 根据影像数据融合的水平阶段,影像融合分为:像元级、特征级和决策级三个层次。为了最大限度的从SPOT5遥感数据中提取森林植被信息,应进行像元级的数据融合,将米的全色波段和10米多光谱数据进行融合。融合得到的新数据既具有全色波段数据的高空间分辨率特征又具有多光谱特征。像元级数据融合的方法多种多样,根据数据融合的目的,即最大限度的突显森林植被信息,应选取B4、B3、B2和PAN波段,根据我们的试验Brovey 融合算法方法比较理想:遥感影像地图将融合好的数据按Rfused、Gfused、Bfused组合,叠加上行政界线、公里格网、坐标、比例尺等辅助信息,按1:1万地形图分幅生成1:1万纸质图作为外业手图。3. 结果和讨论 几何精度利用SPOT5物理模型,采用1:1万地形图和万DEM ,经过正射校正处理,可使影像的几何精度控制在2个像元内(<10米),达到1:1万制图标准要求。为以遥感影像为基础信息源提取林分调查因子、区划林班界线生成大比例尺的林相图、森林分布图提供了几何精度保障。 波段选择对于没有全色波段的情况,SPOT5数据的B342组合有利于森林植被类型的识别。在应用遥感技术进行森林资源调查区划中,林分类型信息提取是最为重要的环节,所以B342波段组合是小班区划和外业手图的最佳组合。 融合效果融合数据技术使SPOT5遥感影像既具有全色波段的高空间分辨率又拥有多光谱数据的光谱分辨率,丰富了遥感影像的信息量。采用Brovey算法使SPOT5遥感影像从色彩、纹理等方面增强了影像的可判读性,提高了小班因子正判率和林分小班的区划精度。参考文献1.周成虎,杨晓梅,骆剑承等.《遥感影像地学理解与分析》,科学出版社,北京,2001,.赵英时.《遥感应用分析原理与方法》,科学出版社,北京,.北京视宝卫星图像有限公司.《专业制图工作室GEOIMAGE用户指南》,2004,.Christine Pohl. Geometric Aspects of Multisensor Image Fusion for Topographic Map Updating in The Humid Tropics, ITC Publication, 1996,世纪遥感与GIS的发展来源: 李德仁 时间: 2005-08-11-23:09 浏览次数: 7921世纪遥感与GIS的发展李德仁(武汉大学测绘遥感信息工程国家重点实验室,武汉市珞瑜路129号,430079)摘要:在20世纪,人类的一大进步是实现了太空对地观测,即可以从空中和太空对人类赖以生存的地球通过非接触传感器的遥感进行观测,并将所得到的数据和信息存储在计算机网络上,为人类社会的可持续发展服务。在短短的30年中,遥感和GIS作为一个边缘交叉学科已发展成为一门科学、技术和经济实体。本文深入地论述了21世纪中遥感的6大发展趋势和GIS的5个发展特征。关键词:发展趋势;航空航天遥感;地理信息系统;对地观测中图法分类号:P208;随着计算机技术、空间技术和信息技术的发展,人类实现了从空中和太空来观测和感知人类赖以生存的地球的理想,并能将所感知到的结果通过计算机网络在全球流通,为人类的生存、繁荣和可持续发展服务。在20世纪后半叶,遥感和地理信息系统作为一门新兴的科学和技术,迅速地成长起来。1 遥感技术的主要发展趋势 航空航天遥感传感器数据获取技术趋向三多(多平台、多传感器、多角度)和三高(高空间分辨率、高光谱分辨率和高时相分辨率)从空中和太空观测地球获取影像是20世纪的重大成果之一,短短几十年,遥感数据获取手段迅猛发展。遥感平台有地球同步轨道卫星(35000km)、太阳同步卫星(600—1000km)、太空飞船(200—300km)、航天飞机(240—350km)、探空火箭(200—1000km),并且还有高、中、低空飞机、升空气球、无人飞机等;传感器有框幅式光学相机、缝隙、全景相机、光机扫描仪、光电扫描仪、CCD线阵、面阵扫描仪、微波散射计雷达测高仪、激光扫描仪和合成孔径雷达等,它们几乎覆盖了可透过大气窗口的所有电磁波段。三行CCD阵列可以同时得到3个角度的扫描成像,EOS Terra卫星上的MISR可同时从9个角度对地成像。卫星遥感的空间分辨率从Ikonos Ⅱ的1m,进一步提高到Quckbird(快鸟)的,高光谱分辨率已达到5—6nm,500—600个波段。在轨的美国EO-1高光谱遥感卫星,具有220个波段,EOS AM-1(Terra)和EOS PM-1(Aqua)卫星上的MODIS具有36个波段的中等分辨率成像光谱仪。时间分辨率的提高主要依赖于小卫星技术的发展,通过发射地球同步轨道卫星和合理分布的小卫星星座,以及传感器的大角度倾斜,可以以1—3d的周期获得感兴趣地区的遥感影像。由于具有全天候、全天时的特点,以及用INSAR和D-INSAR,特别是双天线INSAR进行高精度三位地形及其变化测定的可能性,SAR雷达卫星为全世界各国所普遍关注。例如,美国宇航局的长远计划是要发射一系列太阳同步和地球同步的长波SAR,美国国防部则要发射一系列短波SAR,实现干涉重访问间隔为8d、3d和1d,空间分辨率分别为20m、5m和2m。我国在机载和星载SAR传感器及其应用研究方面正在形成体系。“十五”期间,我国将全方位地推进遥感数据获取的手段,形成自主的高分辨率资源卫星、雷达卫星、测图卫星和对环境与灾害进行实时监测的小卫星群。 航空航天遥感对地定位趋向于不依赖地面控制确定影像目标的实地位置(三维坐标),解决影像目标在哪儿(Where)是摄影测量与遥感的主要任务之一。在已成功用于生产的全自动化GPS空中三角测量的基础上,利用DGPS和INS惯性导航系统的组合,可形成航空/航天影像传感器的位置与姿态的自动测量和稳定装置(POS),从而可实现定点摄影成像和无地面控制的高精度对地直接定位。在航空摄影条件下的精度可达到dm级,在卫星遥感的条件下,其精度可达到m级。该技术的推广应用,将改变目前摄影测量和遥感的作业流程,从而实现实时测图和实时数据库更新。若与高精度激光扫描仪集成,可实现实时三维测量(LIDAR),自动生成数字表面模型(DSM),并可推算出数字高程模型(DEM)。美国NASA在1994年和1997年两次将航天激光测高仪(SLA)安装在航天飞机上,企图建立基于SLA的全球控制点数据库,激光点大小为100m,间隔为750m,每秒10个脉冲;随后又提出了地学激光测高系统(GLAS)计划,已于2002年12月19日将该卫星IICESat(cloud and land elevation satellite)发射上天。该卫星装有激光测距系统、GPS接收机和恒星跟踪姿态测定系统。GLAS发射近红外光(1064nm)和可见绿光(532nm)的短脉冲(4ns)。激光脉冲频率为40次/s,激光点大小实地为70m,间隔为170m,其高程精度要明显高于SRTM,可望达到m级。他们的下一步计划是要在2015年之前使星载LIDAR的激光测高精度达到dm和cm级。法国利用设在全球的54个站点向卫星发射信号,通过测定多普勒频移,以精确解求卫星的空间坐标,具有极高的精度。测定距地球1300km的Topex/Poseidon卫星的高度,精度达到±3cm。用来测定SPOT 4卫星的轨道,3个坐标方向达到±5cm精度,对于SPOT 5和Envisat,可望达到±1m精度。若忽略SPOT 5传感器的角元素,直接进行无地面控制的正射像片制作,精度可达到±15m,完全可以满足国家安全和西部开发的需求。 摄影测量与遥感数据的计算机处理更趋向自动化和智能化从影像数据中自动提取地物目标,解决它的属性和语义(What)是摄影测量与遥感的另一大任务。在已取得影像匹配成果的基础上,影像目标的自动识别技术主要集中在影像融合技术,基于统计和基于结构的目标识别与分类,处理的对象既包括高分辨率影像,也更加注重高光谱影像。随着遥感数据量的增大,数据融合和信息融合技术逐渐成熟。压缩倍率高、速度快的影像数据压缩方法也已商业化。我国学者在这些方面取得了不少可喜的成果。 利用多时像影像数据自动发现地表覆盖的变化趋向实时化利用遥感影像自动进行变化监测(What change)关系到我国的经济建设和国防建设。过去人工方法投入大,周期长。随着各类空间数据库的建立和大量新的影像数据源的出现,实时自动化监测已成为研究的一个热点。自动变化监测研究包括利用新旧影像(DOM)的对比、新影像与旧数字地图(DLS)的对比来自动发现变化和更新数据库。目前的变化监测是先将新影像与旧影像(或数字地图)进行配准,然后再提取变化目标,这在精度、速度与自动化处理方面都有不足之处。笔者提出了把配准与变化监测同步的整体处理[1]。最理想的方法是将影像目标三维重建与变化监测一起进行,实现三维变化监测和自动更新。进一步的发展则是利用智能传感器,将数据处理在轨完成,发送回来的直接为信息,而不一定为影像数据。 摄影测量与遥感在构建“数字地球”、“数字中国”、“数字省市”和“数字文化遗产”中正在发挥愈来愈大的作用“数字地球”概念是在全球信息化浪潮推进下形成的。1999年12月在北京成功地召开了第一届国际“数字地球”大会后,我国正积极推进“数字中国”和“数字省市”的建设,2001年国家测绘局完成了构建“数字中国”地理空间基础框架的总体战略研究。在已完成1∶100万和1∶25万全国空间数据库的基础上,2001年全国各省市测绘局开始1∶5万空间数据库的建库工作。在这个数据量达11TB的巨型数据库中,摄影测量与遥感将用来建设DOM(数字正射影像)、DEM(数字高程模型)、DLG(数字线划图)和CP(控制点数据库)。如果要建立全国1m分辨率影像数据库,其数据量将达到60TB。如果整个“数字地球”均达到1m分辨率,其数据量之大可想而知。本世纪内可望建成这一分辨率的数字地球。“数字文化遗产”是目前联合国和许多国家关心的一个问题,涉及到近景成像、计算机视觉和虚拟现实技术。在近景成像和近景三位量测方面,有室内各种三维激光扫描与成像仪器,还可以直接由视频摄像机的系列图像获取目标场三维重建信息。它们所获取的数据经过计算机自动处理后,可以在虚拟现实技术支持下形成文化遗迹的三维仿真,而且可以按照时间序列,将历史文化在时间隧道中再现,对文化遗产保护、复原与研究具有重要意义。 全定量化遥感方法将走向实用从遥感科学的本质讲,通过对地球表层(包括岩石圈、水圈、大气圈和生物圈4大圈层)的遥感,其目的是为了获得有关地物目标的几何与物理特性,所以需要通过全定量化遥感方法进行反演。几何方程式是有显式表示的数学方程,而物理方程一直是隐式。目前的遥感解译与目标识别并没有通过物理方程反演,而是采用了基于灰度或加上一定知识的统计、结构和纹理的影像分析方法。但随着对成像机理、地物波谱反射特征、大气模型、气溶胶的研究深入和数据积累,多角度、多传感器、高光谱及雷达卫星遥感技术的成熟,相信在21世纪,估计几何与物理方程式的全定量化遥感方法将逐步由理论研究走向实用化,遥感基础理论研究将迈上新的台阶。只有实现了遥感定量化,才可能真正实现自动化和实时化。2 GIS技术的主要发展趋势 空间数据库趋向图形、影像和DEM三库一体化和面向对象[2]GIS发展曾经历过栅格、矢量两个不同数据结构发展阶段,目前随着高分辨率卫星遥感数据的飞快增长和数字地球、数码城市的需求,形成了面向对象的数据模型和三库(图形矢量库、影像栅格库和DEM格网库)一体化的数据结构。这样的数据库结构使GIS的发展更加趋向自然化、逼真化,更加贴近用户。以面向应用的GIS软件为前台,以大型关系数据库(Oracle 8i,9i等)为后台数据库管理,成为当前GIS技术的主流趋势。 空间数据表达趋向多比例尺、多尺度、动态多位和实时三维可视化在传统的GIS中,空间数据是以二维形式存储并挂接相应的属性数据。目前,空间数据表达的趋势是基于金字塔和LOD(level of detail)技术的多比例尺空间数据库,在不同尺度表示时可自动显示出相应比例尺或相应分辨率的数据,多比例尺数据集的跨度要比传统地图的比例尺大,在显示不同比例尺数据时,可采用LOD或地图综合技术。真三维GIS的空间数据要存储三维坐标。动态GIS在土地变更调查、土地覆盖变化监测中已有较好的应用,真四维的时空GIS将有望从理论研究转入实用阶段。基于三库一体化的时空3D可视化技术发展势头迅猛,已能再PC机上实现GIS环境下的三维建筑物室外室内漫游、信息查询、空间分析、剖面分析和阴影分析等,基于虚拟现实技术的真三维GIS将使人们在现实空间外,可以同时拥有一个Cyber空间。 空间分析和辅助决策智能化需要利用数据挖掘方法从空间数据库和属性数据库中发现更多的有用知识GIS是以应用导向的空间信息技术,空间分析与辅助决策支持是GIS的高水平应用,它需要基于知识的智能系统。知识的获取是专家系统中最困难的任务。随着各种类型数据库的建立,从数据库中挖掘知识成为当今计算机界一个非常引人注目的课题。从GIS空间数据库中发现的知识可以有效的支持遥感图像解译,以解决“同物异谱”和“同谱异物”的问题。反过来,从属性数据库中挖掘的知识又具有优化资源配置等一些列空间分析的功能[3]。尽管数据挖掘和知识发现这一命题仍处于理论研究阶段,但随着数据库的快速增大和对数据挖掘工具的深入研究,其应用前景是不可估量的。 通过Web服务器和WAP服务器的互联网和移动GIS将推进联邦数据库和互操作的研究及地学信息服务事业随着计算机通讯网络(包括有线和无线网)的大容量和高速化,GIS已成为在网络上的分布式异构系统。许多不同单位、不同组织维护管理的既独立又互联互用的联邦数据库,将可提供全社会各行各业的应用需要。因此,联邦数据库和互操作(federal databases & interoperability)问题成为当前国际GIS联合研究的一个热点。互操作意味着数据库中数据的直接共享,GIS规律功能模块的互操作与共享,以及多点之间的相同工作,这方面的研究已显示出明显的成效。未来的GIS用户将可能在网络上缴纳为其需要所选用数据和软件功能模块的使用费,而不必购买这个数据库和整套的GIS软硬件,这些成果产生的直接效果是GIS应用将走向地学信息服务。目前已兴起的LBS和MLS,即基于位置的服务和移动定位服务,突出地反映了这种变化趋势。它引起的革命性变化使GIS将走出研究院所和政府机关,成为全社会人人具备的信息服务工具。我国目前已有2亿个手机用户,若每人每月为MLS支付10元费用,全国一年的产值将达到240亿。可以预测在不久的将来,地学信息将能随时随地为任何人和任何事情进行4A服务(geo-in-formation for anyone and anything at anywhere and anytime)。 地理信息科学的研究有望在本世纪形成较完整的理论框架体系笔者曾扼要地叙述了地球空间信息科学的7大理论问题[4]:(1)地球空间信息的基准,包括几何基准、物理基准和时间基准;(2)地球空间信息标准,包括空间数据采集、存储与交换标准、空间数据精度与质量标准、空间信息的分类与代码标准、空间信息的安全
遥感学报发一篇论文什么水。水平很高的,《遥感学报》已被荷兰斯高帕斯数据库(Scopus)、EBSCO数据库、俄罗斯文摘杂志(AJ)、波兰哥白尼索引(IC)、日本科学技术文献数据库(JST)、美国乌利希期刊指南(UPD)、美国剑桥科学文摘(CSA)等国际知名数据库收录。《遥感学报》的学术水平得到国际认可,受到国际遥感及相关学术领域的更多关注。
森林资源调查中SPOT5遥感图像处理方法探讨王照利、黄生、张敏中、马胜利(国家林业局西北林业规划设计院,遥感计算中心,西安710048)本文发表于<陕西林业科技>2005 摘要: 目前,多光谱、高空间分辨率的SPOT5卫星遥感数据被广泛应用到森林资源调查中。本文结合SPOT5遥感数据的特点,根据森林资源调查的需要,从遥感数据的正射校正、波段组合、融合处理和数据变换处理等方面探讨了SPOT5数据的处理和信息提取。探讨性地提出了适应于森林资源调查的SPOT5遥感数据处理方法。 关键词:SPOT5 遥感数据,森林资源调查、数据处理DISCUSSION ON SPOT5 IMAGE DATA PROCESSING FOR FOREST INVENTORYWang Zhaoli, Huangsheng,Zhangminzhong,Ma Shengli(Northwest Institute for Forest Inventory, Planning &Design, Xi’an China 710048) Abstract: Now days, high spatial resolution and multispectral SPOT5 image data are widely applied in forest inventory in China. Based on the characteristics of SPOT5 image and requirements of forest inventory, this paper discusses the processing procedures of ordering image data, ortho-rectification, image bands composition and image data fusion. The complete steps of image processing for forest inventory are words: SPOT5 image data,forest inventory, data processing 前言 卫星遥感影像具有空间宏观性、视角广、多分辨率(光谱和空间)、多时相、周期性、信息量丰富等特点,所以卫星遥感影像既可以提供森林资源的宏观空间分布信息又能提供局部的详细信息以及随时间、空间变化的信息等[1]。目前在林业领域卫星遥感数据被广泛的应用于不同尺度层次的森林资源调查、资源监测、病虫害、火灾监测等方面。2002年5月法国SPOT地球观测卫星系列之5号卫星(即SPOT5星)发射。SPOT5遥感数据的多光谱波段空间分辨率为10米(短波红外空间分辨率为20米),但全色波段空间分辨率达到米。SPOT5遥感数据的高空间分辨率和多光谱分辨率为森林资源调查提供了丰富的、可靠的、高精度的基础数据源。从性价比分析,在其他高分辨率遥感数据目前比较昂贵的状况下,SPOT5遥感数据比较适宜应用于大面积的森林资源调查,可大幅度的森林调查的减少外业工作量、提高工作效率。在我国SPOT5卫星数据已被大量地应用于森林资源调查工作中,尤其,是在森林资源“二类”调查中被作基本的森林资源信息源提取各类信息。针对于将多光谱分辨率和高空间分辨率的SPOT5遥感数据应用于森林资源调查的数据处理技术和方法鲜有报道。本文总结工作实践,结合SPOT5遥感数据的特点,根据森林资源调查的需要,从遥感数据的订购、正射校正、波段组合、融合处理和数据变换处理等方面探讨了SPOT5数据的基本处理方法。 1.SPOT5卫星遥感数据特点 SPOT卫星系统采用线性阵列传感器和推扫式扫描技术,具有旋转式平面镜可以进行倾斜观察获得倾斜图像和立体像对。采用与太阳同步的近极地的椭圆形轨道,轨道高度约832Km,轨道倾角 ,每天绕地球14圈多,重复覆盖周期26天[2]。由于有倾斜观测功能,使重复覆盖周期减少到2-3天。SPOT5卫星载有2台高分辨率几何成像仪(HRG)、1台高分辨率立体成像装置(HRS)和1台宽视域植被探测仪(VGT)。高分辨率几何成像仪的波段选择是总结了多年的研究成果,认为HRG的波段设置(见表1)足以取得辨别作物和植被类型的最佳效果。本文主要探讨HRG高空间分辨率数据的处理。 2.SPOT5数据的处理方法和过程 SPOT5数据处理工作流程: 遥感数据的订购 订购数据时,用户需向数据代理商提供购买区域的四个角的大地坐标或者数据的景号(PATH/ROW)。特别应该注意数据订购时间和用户拿到数据之间有时间差,间隔时间长短因用户的要求、天气、卫星重复覆盖周期而异。相对于其他卫星数据,比较有利的一面是SPOT5卫星装置有旋转式平面镜可以进行倾斜观察,用户可向代理商申请红色编程提前得到调查区域的遥感数据,但要支付编程费。对于遥感数据的时相、云量、入射角、阴影量、是否购买高空间分辨率的全色波段等用户根据自己具体的工作需要向代理商提出限制要求。 根据我们对SPOT5遥感数据的使用,对于森林资源调查,北方9,10月份和11月初的遥感影像比较适宜。代理商向用户提供经过处理的不同级别的影像产品,在森林资源调查中建议购买SPOT1A级产品,用户可根据自己的工作需要进行处理,同时也可减少费用。 基础数据准备 大比例尺地形图和高精度DEM是进行SPOT5遥感数据高精度正射校正必需的基础地理数据。建议购买1:10000地形图和1:25000数字高程模型(DEM)。 将1:1万地形图扫描,扫描分辨率设置为300DPI。将扫描好的地形图进行几何精纠正,纠正精度控制在毫米内。从测绘部门购买的1:1万地形图为北京54坐标系3度分带高斯克吕格投影,而1:万DEM为北京54坐标系6度分带投影。在数据准备时,将校正好的1:1万地形图通过换带转换转成和DEM一致的6度分带投影。 对于没有1:1万地形图的地区,建议使用差分GPS接收机采集地面控制点。 几何正射校正 正射校正过程应用了法国SPOT公司发行的GEOIMAGE软件。GEOIMAGE软件有针对SPOT5卫星数据开发的SPOT5物理模型。模型模块自动读取DEM信息。SPOT 物理模型可读取卫星在获取遥感数据的瞬间状态参数,这些参数存贮在数据的头文件中[3]。卫星状态参数包括:卫星成像瞬间的经纬度、高度、倾角等。卫星状态参数能够帮助提高几何校正的精度。 以校正好的1:1万地形图为基准,在影像图上找出和地形图上地物相匹配的明显地物作为地面控制点。在进行正射校正时,应先进行全色波段数据校正,然后以校正好的全色波段数据为基准进行多光谱数据校正。以全色波段数据为基准校正多光谱波段就比较容易校正,且能提高两者的匹配精度。地面控制点应分布均匀,影像的边缘部分布要有控制点分布,同时在不同的高程范围最好都有控制点。地面控制点的数量因地形地貌的复杂程度而定,根据我们的经验,一景60KmX60Km的SPOT5数据,一般地势平缓的地区20个左右控制点即可达到满意的结果,在高山区25个左右控制点就可使正射校正精度满足要求。重采样方法采用双线性内插法。 辐射校正 用户购买的SPOT5的各级数据,数据提供商已经根据卫星的记录参数对遥感数据做了辐射校正,即消除了传感器自身引起的、大气辐射引起的辐射噪声。若果影像存在薄雾或地形高差较大引起的辐射误差情况,用户应进一步进行辐射校正处理。薄雾的简单消除原理是基于近红外波段不受大气辐射影响,清澈的水体或死阴影区的数值应为零。从各波段数据中减去近红外波段的水体或阴影的不为零值。地形起伏引起的辐射误差校正公式: f (x,y)=g(x,y)/cosa,g(x,y)为坡度为a的倾斜面上的地物影像;f (x,y)为校正后的影像。由于坡度因子参与校正所以需要DEM支持。 波段组合 根据SPOT5数据波谱特征(表1),各波段分别记录反映了植被的不同特征方面:B4(SWIR)短波红外反映植物和土壤的含水量,利于植被水分状况和长势分析;B3(NIR)近红外波段对植被类别、密度、生长力、病虫害等的变化敏感;B2(RED)红光波段对植被的覆盖度、植被的生长状况敏感;B1(VIS)可见光波段对植物的叶绿素和叶绿素浓度敏感。经过比较分析和实际应用发现SPOT5的B3、B4、B2波段组合对植被类型的识别要优于B3、B2和B1的组合。但由于B4波段的空间分辨率为20米,使B342组合对植被空间几何细节表达没有B321组合清晰,例如林缘界线信息表达方面B321要优于B342。 影像数据融合 对于购买有高空间分辨率全色波段数据的用户,进行数据融合是必不可少的。影像数据融合能够综合不同波段、不同空间分辨率数据(层)的特征,融合后的数据具有更丰富、更可靠的信息[4]。 根据影像数据融合的水平阶段,影像融合分为:像元级、特征级和决策级三个层次。为了最大限度的从SPOT5遥感数据中提取森林植被信息,应进行像元级的数据融合,将米的全色波段和10米多光谱数据进行融合。融合得到的新数据既具有全色波段数据的高空间分辨率特征又具有多光谱特征。像元级数据融合的方法多种多样,根据数据融合的目的,即最大限度的突显森林植被信息,应选取B4、B3、B2和PAN波段,根据我们的试验Brovey 融合算法方法比较理想: 遥感影像地图 将融合好的数据按Rfused、Gfused、Bfused组合,叠加上行政界线、公里格网、坐标、比例尺等辅助信息,按1:1万地形图分幅生成1:1万纸质图作为外业手图。 3. 结果和讨论 几何精度 利用SPOT5物理模型,采用1:1万地形图和万DEM ,经过正射校正处理,可使影像的几何精度控制在2个像元内(<10米),达到1:1万制图标准要求。为以遥感影像为基础信息源提取林分调查因子、区划林班界线生成大比例尺的林相图、森林分布图提供了几何精度保障。 波段选择 对于没有全色波段的情况,SPOT5数据的B342组合有利于森林植被类型的识别。在应用遥感技术进行森林资源调查区划中,林分类型信息提取是最为重要的环节,所以B342波段组合是小班区划和外业手图的最佳组合。 融合效果 融合数据技术使SPOT5遥感影像既具有全色波段的高空间分辨率又拥有多光谱数据的光谱分辨率,丰富了遥感影像的信息量。采用Brovey算法使SPOT5遥感影像从色彩、纹理等方面增强了影像的可判读性,提高了小班因子正判率和林分小班的区划精度。 参考文献 1.周成虎,杨晓梅,骆剑承等.《遥感影像地学理解与分析》,科学出版社,北京,2001,3-4. 2.赵英时.《遥感应用分析原理与方法》,科学出版社,北京, 3.北京视宝卫星图像有限公司.《专业制图工作室GEOIMAGE用户指南》,2004,68-70. 4.Christine Pohl. Geometric Aspects of Multisensor Image Fusion for Topographic Map Updating in The Humid Tropics, ITC Publication, 1996,世纪遥感与GIS的发展 来源: 李德仁 时间: 2005-08-11-23:09 浏览次数: 79 21世纪遥感与GIS的发展李德仁 (武汉大学测绘遥感信息工程国家重点实验室,武汉市珞瑜路129号,430079) 摘要:在20世纪,人类的一大进步是实现了太空对地观测,即可以从空中和太空对人类赖以生存的地球通过非接触传感器的遥感进行观测,并将所得到的数据和信息存储在计算机网络上,为人类社会的可持续发展服务。在短短的30年中,遥感和GIS作为一个边缘交叉学科已发展成为一门科学、技术和经济实体。本文深入地论述了21世纪中遥感的6大发展趋势和GIS的5个发展特征。 关键词:发展趋势;航空航天遥感;地理信息系统;对地观测 中图法分类号:P208; 随着计算机技术、空间技术和信息技术的发展,人类实现了从空中和太空来观测和感知人类赖以生存的地球的理想,并能将所感知到的结果通过计算机网络在全球流通,为人类的生存、繁荣和可持续发展服务。在20世纪后半叶,遥感和地理信息系统作为一门新兴的科学和技术,迅速地成长起来。 1 遥感技术的主要发展趋势 航空航天遥感传感器数据获取技术趋向三多(多平台、多传感器、多角度)和三高(高空间分辨率、高光谱分辨率和高时相分辨率) 从空中和太空观测地球获取影像是20世纪的重大成果之一,短短几十年,遥感数据获取手段迅猛发展。遥感平台有地球同步轨道卫星(35000km)、太阳同步卫星(600—1000km)、太空飞船(200—300km)、航天飞机(240—350km)、探空火箭(200—1000km),并且还有高、中、低空飞机、升空气球、无人飞机等;传感器有框幅式光学相机、缝隙、全景相机、光机扫描仪、光电扫描仪、CCD线阵、面阵扫描仪、微波散射计雷达测高仪、激光扫描仪和合成孔径雷达等,它们几乎覆盖了可透过大气窗口的所有电磁波段。三行CCD阵列可以同时得到3个角度的扫描成像,EOS Terra卫星上的MISR可同时从9个角度对地成像。 卫星遥感的空间分辨率从Ikonos Ⅱ的1m,进一步提高到Quckbird(快鸟)的,高光谱分辨率已达到5—6nm,500—600个波段。在轨的美国EO-1高光谱遥感卫星,具有220个波段,EOS AM-1(Terra)和EOS PM-1(Aqua)卫星上的MODIS具有36个波段的中等分辨率成像光谱仪。时间分辨率的提高主要依赖于小卫星技术的发展,通过发射地球同步轨道卫星和合理分布的小卫星星座,以及传感器的大角度倾斜,可以以1—3d的周期获得感兴趣地区的遥感影像。由于具有全天候、全天时的特点,以及用INSAR和D-INSAR,特别是双天线INSAR进行高精度三位地形及其变化测定的可能性,SAR雷达卫星为全世界各国所普遍关注。例如,美国宇航局的长远计划是要发射一系列太阳同步和地球同步的长波SAR,美国国防部则要发射一系列短波SAR,实现干涉重访问间隔为8d、3d和1d,空间分辨率分别为20m、5m和2m。我国在机载和星载SAR传感器及其应用研究方面正在形成体系。“十五”期间,我国将全方位地推进遥感数据获取的手段,形成自主的高分辨率资源卫星、雷达卫星、测图卫星和对环境与灾害进行实时监测的小卫星群。 航空航天遥感对地定位趋向于不依赖地面控制 确定影像目标的实地位置(三维坐标),解决影像目标在哪儿(Where)是摄影测量与遥感的主要任务之一。在已成功用于生产的全自动化GPS空中三角测量的基础上,利用DGPS和INS惯性导航系统的组合,可形成航空/航天影像传感器的位置与姿态的自动测量和稳定装置(POS),从而可实现定点摄影成像和无地面控制的高精度对地直接定位。在航空摄影条件下的精度可达到dm级,在卫星遥感的条件下,其精度可达到m级。该技术的推广应用,将改变目前摄影测量和遥感的作业流程,从而实现实时测图和实时数据库更新。若与高精度激光扫描仪集成,可实现实时三维测量(LIDAR),自动生成数字表面模型(DSM),并可推算出数字高程模型(DEM)。 美国NASA在1994年和1997年两次将航天激光测高仪(SLA)安装在航天飞机上,企图建立基于SLA的全球控制点数据库,激光点大小为100m,间隔为750m,每秒10个脉冲;随后又提出了地学激光测高系统(GLAS)计划,已于2002年12月19日将该卫星IICESat(cloud and land elevation satellite)发射上天。该卫星装有激光测距系统、GPS接收机和恒星跟踪姿态测定系统。GLAS发射近红外光(1064nm)和可见绿光(532nm)的短脉冲(4ns)。激光脉冲频率为40次/s,激光点大小实地为70m,间隔为170m,其高程精度要明显高于SRTM,可望达到m级。他们的下一步计划是要在2015年之前使星载LIDAR的激光测高精度达到dm和cm级。 法国利用设在全球的54个站点向卫星发射信号,通过测定多普勒频移,以精确解求卫星的空间坐标,具有极高的精度。测定距地球1300km的Topex/Poseidon卫星的高度,精度达到±3cm。用来测定SPOT 4卫星的轨道,3个坐标方向达到±5cm精度,对于SPOT 5和Envisat,可望达到±1m精度。若忽略SPOT 5传感器的角元素,直接进行无地面控制的正射像片制作,精度可达到±15m,完全可以满足国家安全和西部开发的需求。 摄影测量与遥感数据的计算机处理更趋向自动化和智能化 从影像数据中自动提取地物目标,解决它的属性和语义(What)是摄影测量与遥感的另一大任务。在已取得影像匹配成果的基础上,影像目标的自动识别技术主要集中在影像融合技术,基于统计和基于结构的目标识别与分类,处理的对象既包括高分辨率影像,也更加注重高光谱影像。随着遥感数据量的增大,数据融合和信息融合技术逐渐成熟。压缩倍率高、速度快的影像数据压缩方法也已商业化。我国学者在这些方面取得了不少可喜的成果。 利用多时像影像数据自动发现地表覆盖的变化趋向实时化 利用遥感影像自动进行变化监测(What change)关系到我国的经济建设和国防建设。过去人工方法投入大,周期长。随着各类空间数据库的建立和大量新的影像数据源的出现,实时自动化监测已成为研究的一个热点。 自动变化监测研究包括利用新旧影像(DOM)的对比、新影像与旧数字地图(DLS)的对比来自动发现变化和更新数据库。目前的变化监测是先将新影像与旧影像(或数字地图)进行配准,然后再提取变化目标,这在精度、速度与自动化处理方面都有不足之处。笔者提出了把配准与变化监测同步的整体处理[1]。最理想的方法是将影像目标三维重建与变化监测一起进行,实现三维变化监测和自动更新。进一步的发展则是利用智能传感器,将数据处理在轨完成,发送回来的直接为信息,而不一定为影像数据。 摄影测量与遥感在构建“数字地球”、“数字中国”、“数字省市”和“数字文化遗产”中正在发挥愈来愈大的作用 “数字地球”概念是在全球信息化浪潮推进下形成的。1999年12月在北京成功地召开了第一届国际“数字地球”大会后,我国正积极推进“数字中国”和“数字省市”的建设,2001年国家测绘局完成了构建“数字中国”地理空间基础框架的总体战略研究。在已完成1∶100万和1∶25万全国空间数据库的基础上,2001年全国各省市测绘局开始1∶5万空间数据库的建库工作。在这个数据量达11TB的巨型数据库中,摄影测量与遥感将用来建设DOM(数字正射影像)、DEM(数字高程模型)、DLG(数字线划图)和CP(控制点数据库)。如果要建立全国1m分辨率影像数据库,其数据量将达到60TB。如果整个“数字地球”均达到1m分辨率,其数据量之大可想而知。本世纪内可望建成这一分辨率的数字地球。 “数字文化遗产”是目前联合国和许多国家关心的一个问题,涉及到近景成像、计算机视觉和虚拟现实技术。在近景成像和近景三位量测方面,有室内各种三维激光扫描与成像仪器,还可以直接由视频摄像机的系列图像获取目标场三维重建信息。它们所获取的数据经过计算机自动处理后,可以在虚拟现实技术支持下形成文化遗迹的三维仿真,而且可以按照时间序列,将历史文化在时间隧道中再现,对文化遗产保护、复原与研究具有重要意义。 全定量化遥感方法将走向实用 从遥感科学的本质讲,通过对地球表层(包括岩石圈、水圈、大气圈和生物圈4大圈层)的遥感,其目的是为了获得有关地物目标的几何与物理特性,所以需要通过全定量化遥感方法进行反演。几何方程式是有显式表示的数学方程,而物理方程一直是隐式。目前的遥感解译与目标识别并没有通过物理方程反演,而是采用了基于灰度或加上一定知识的统计、结构和纹理的影像分析方法。但随着对成像机理、地物波谱反射特征、大气模型、气溶胶的研究深入和数据积累,多角度、多传感器、高光谱及雷达卫星遥感技术的成熟,相信在21世纪,估计几何与物理方程式的全定量化遥感方法将逐步由理论研究走向实用化,遥感基础理论研究将迈上新的台阶。只有实现了遥感定量化,才可能真正实现自动化和实时化。 2 GIS技术的主要发展趋势 空间数据库趋向图形、影像和DEM三库一体化和面向对象[2] GIS发展曾经历过栅格、矢量两个不同数据结构发展阶段,目前随着高分辨率卫星遥感数据的飞快增长和数字地球、数码城市的需求,形成了面向对象的数据模型和三库(图形矢量库、影像栅格库和DEM格网库)一体化的数据结构。这样的数据库结构使GIS的发展更加趋向自然化、逼真化,更加贴近用户。以面向应用的GIS软件为前台,以大型关系数据库(Oracle 8i,9i等)为后台数据库管理,成为当前GIS技术的主流趋势。 空间数据表达趋向多比例尺、多尺度、动态多位和实时三维可视化 在传统的GIS中,空间数据是以二维形式存储并挂接相应的属性数据。目前,空间数据表达的趋势是基于金字塔和LOD(level of detail)技术的多比例尺空间数据库,在不同尺度表示时可自动显示出相应比例尺或相应分辨率的数据,多比例尺数据集的跨度要比传统地图的比例尺大,在显示不同比例尺数据时,可采用LOD或地图综合技术。真三维GIS的空间数据要存储三维坐标。动态GIS在土地变更调查、土地覆盖变化监测中已有较好的应用,真四维的时空GIS将有望从理论研究转入实用阶段。基于三库一体化的时空3D可视化技术发展势头迅猛,已能再PC机上实现GIS环境下的三维建筑物室外室内漫游、信息查询、空间分析、剖面分析和阴影分析等,基于虚拟现实技术的真三维GIS将使人们在现实空间外,可以同时拥有一个Cyber空间。 空间分析和辅助决策智能化需要利用数据挖掘方法从空间数据库和属性数据库中发现更多的有用知识 GIS是以应用导向的空间信息技术,空间分析与辅助决策支持是GIS的高水平应用,它需要基于知识的智能系统。知识的获取是专家系统中最困难的任务。随着各种类型数据库的建立,从数据库中挖掘知识成为当今计算机界一个非常引人注目的课题。从GIS空间数据库中发现的知识可以有效的支持遥感图像解译,以解决“同物异谱”和“同谱异物”的问题。反过来,从属性数据库中挖掘的知识又具有优化资源配置等一些列空间分析的功能[3]。尽管数据挖掘和知识发现这一命题仍处于理论研究阶段,但随着数据库的快速增大和对数据挖掘工具的深入研究,其应用前景是不可估量的。 通过Web服务器和WAP服务器的互联网和移动GIS将推进联邦数据库和互操作的研究及地学信息服务事业 随着计算机通讯网络(包括有线和无线网)的大容量和高速化,GIS已成为在网络上的分布式异构系统。许多不同单位、不同组织维护管理的既独立又互联互用的联邦数据库,将可提供全社会各行各业的应用需要。因此,联邦数据库和互操作(federal databases & interoperability)问题成为当前国际GIS联合研究的一个热点。互操作意味着数据库中数据的直接共享,GIS规律功能模块的互操作与共享,以及多点之间的相同工作,这方面的研究已显示出明显的成效。未来的GIS用户将可能在网络上缴纳为其需要所选用数据和软件功能模块的使用费,而不必购买这个数据库和整套的GIS软硬件,这些成果产生的直接效果是GIS应用将走向地学信息服务。 目前已兴起的LBS和MLS,即基于位置的服务和移动定位服务,突出地反映了这种变化趋势。它引起的革命性变化使GIS将走出研究院所和政府机关,成为全社会人人具备的信息服务工具。我国目前已有2亿个手机用户,若每人每月为MLS支付10元费用,全国一年的产值将达到240亿。可以预测在不久的将来,地学信息将能随时随地为任何人和任何事情进行4A服务(geo-in-formation for anyone and anything at anywhere and anytime)。 地理信息科学的研究有望在本世纪形成较完整的理论框架体系 笔者曾扼要地叙述了地球空间信息科学的7大理论问题[4]:(1)地球空间信息的基准,包括几何基准、物理基准和时间基准;(2)地球空间信息标准,包括空间数据采集、存储与交换标准、空间数据精度与质量标准、空间信息的分类与代码标准、空间信息的安全
安徽农业大学学报安全与环境学报氨基酸和生物资源半导体光电半导体技术半导体学报包装工程爆破爆炸与冲击北京大学学报. 医学版北京大学学报. 自然科学版北京工业大学学报北京航空航天大学学报北京化工大学学报. 自然科学版北京交通大学学报. 自然科学版北京科技大学学报北京理工大学学报北京林业大学学报北京师范大学学报. 自然科学版北京医学北京邮电大学学报北京中医药大学学报表面技术冰川冻土兵工学报兵器材料科学与工程病毒学报波谱学杂志玻璃钢/复合材料材料保护材料导报材料工程材料科学与工程学报材料科学与工艺材料热处理学报材料研究学报蚕业科学草地学报草业科学草业学报测绘科学测绘通报测绘学报测井技术测控技术茶叶科学长安大学学报. 自然科学版长江科学院院报长江流域资源与环境肠外与肠内营养沉积学报沉积与特提斯地质成都理工大学学报. 自然科学版城市规划学刊城市环境与城市生态传感技术学报传感器与微系统船舶工程纯粹数学与应用数学磁性材料及器件催化学报大地测量与地球动力学大地构造与成矿学大豆科学大连海事大学学报大连理工大学学报大连水产学院学报大气科学弹道学报弹箭与制导学报淡水渔业导弹与航天运载技术低温工程低温物理学报低温与超导地层学杂志地理科学地理科学进展地理学报地理研究地理与地理信息科学地球化学地球科学地球科学进展地球物理学报地球物理学进展地球信息科学地球学报地球与环境地学前缘地震地震地质地震工程与工程振动地震学报地震研究地质科技情报地质科学地质力学学报地质论评地质通报地质学报地质与勘探地质找矿论丛第二军医大学学报第三军医大学学报第四纪研究第四军医大学学报电波科学学报电池电镀与环保电镀与涂饰电工电能新技术电工技术学报电化学电机与控制学报电力电子技术电力系统及其自动化学报电力系统自动化电力自动化设备电路与系统学报电气传动电网技术电源技术电子测量与仪器学报电子技术应用电子科技大学学报电子显微学报电子学报电子与信息学报电子元件与材料东北大学学报. 自然科学版东北林业大学学报东北农业大学学报东北师大学报. 自然科学版东华大学学报. 自然科学版东南大学学报. 自然科学版动力工程动物分类学报动物学报动物学研究动物学杂志毒理学杂志锻压技术发光学报防灾减灾工程学报纺织学报飞行力学非金属矿分析测试学报分析化学分析科学学报分析试验室分析仪器分子催化分子科学学报分子细胞生物学报分子植物育种粉末冶金技术福建林学院学报福建农林大学学报. 自然科学版福建师范大学学报. 自然科学版福州大学学报. 自然科学版辐射防护辐射研究与辐射工艺学报腐蚀科学与防护技术复旦学报. 医学版复旦学报. 自然科学版复合材料学报复杂系统与复杂性科学干旱地区农业研究干旱区地理干旱区研究干旱区资源与环境甘肃农业大学学报钢铁钢铁研究学报高等学校化学学报高等学校计算数学学报高电压技术高分子材料科学与工程高分子通报高分子学报高技术通讯高校地质学报高校化学工程学报高校应用数学学报. A辑高压物理学报高原气象给水排水工程地质学报工程勘察工程力学工程热物理学报工程设计学报工程数学学报工程塑料应用工程图学学报工业工程工业工程与管理工业建筑工业水处理工业微生物工业卫生与职业病公共管理学报公路交通科技功能材料功能材料与器件学报功能高分子学报古地理学报古脊椎动物学报古生物学报固体电子学研究与进展固体火箭技术固体力学学报管理工程学报管理科学管理科学学报管理评论管理世界管理学报灌溉排水学报光电工程光电子·激光光电子技术光谱实验室光谱学与光谱分析光散射学报光通信技术光通信研究光学技术光学精密工程光学学报光子学报广东农业科学广东微量元素科学广东医学广西大学学报. 自然科学版广西植物硅酸盐通报硅酸盐学报贵金属贵州农业科学国防科技大学学报国土资源遥感果树学报过程工程学报哈尔滨工程大学学报哈尔滨工业大学学报哈尔滨医科大学学报海洋地质与第四纪地质海洋工程海洋湖沼通报海洋环境科学海洋科学海洋科学进展海洋通报海洋学报海洋渔业海洋与湖沼含能材料焊接学报航空材料学报航空动力学报航空精密制造技术航空学报航天控制航天医学与医学工程合成化学合成纤维工业合成橡胶工业合肥工业大学学报. 自然科学版河北工业大学学报河北农业大学学报河海大学学报. 自然科学版河南大学学报. 自然科学版河南农业大学学报河南师范大学学报. 自然科学版核电子学与探测技术核动力工程核化学与放射化学核技术核聚变与等离子体物理核科学与工程核农学报黑龙江农业科学红外技术红外与毫米波学报红外与激光工程湖北农业科学湖泊科学湖南大学学报. 自然科学版湖南农业大学学报湖南农业科学华北农学报华东理工大学学报. 自然科学版华东师范大学学报. 自然科学版华南理工大学学报. 自然科学版华南农业大学学报. 自然科学版华西口腔医学杂志华西药学杂志华中科技大学学报. 医学版华中科技大学学报. 自然科学版华中农业大学学报华中师范大学学报. 自然科学版化工环保化工进展化工新型材料化工学报化学反应工程与工艺化学工程化学进展化学世界化学试剂化学通报化学学报化学研究与应用环境工程环境工程学报环境化学环境科学环境科学学报环境科学研究环境科学与技术环境昆虫学报环境污染与防治环境与健康杂志会计研究火工品火力与指挥控制火灾科学火炸药学报机器人机械传动机械工程材料机械工程学报机械科学与技术机械强度机械设计机械设计与研究基础医学与临床基因组学与应用生物学激光技术激光生物学报激光与光电子学进展激光与红外激光杂志吉林大学学报. 地球科学版吉林大学学报. 工学版吉林大学学报. 理学版吉林大学学报. 医学版吉林农业大学学报吉林农业科学极地研究计量学报计算机测量与控制计算机仿真计算机辅助设计与图形学学报计算机工程计算机工程与科学计算机工程与设计计算机工程与应用计算机集成制造系统计算机科学计算机科学与探索计算机系统应用计算机学报计算机研究与发展计算机应用计算机应用研究计算机与应用化学计算力学学报计算数学计算物理暨南大学学报. 自然科学与医学版检验医学建筑材料学报建筑结构建筑结构学报江苏农业科学江苏农业学报江苏医药江西农业大学学报交通运输工程学报交通运输系统工程与信息解放军医学杂志解剖学报解剖学杂志介入放射学杂志金融研究金属热处理金属学报经济数学精细化工精细石油化工军事医学科学院院刊菌物学报科技导报(北京)科学通报科学学研究科学学与科学技术管理科研管理空间结构空间科学学报空军工程大学学报. 自然科学版空气动力学学报控制工程控制理论与应用控制与决策矿床地质矿物学报矿物岩石矿物岩石地球化学通报矿冶工程昆虫分类学报昆虫学报昆虫知识兰州大学学报. 自然科学版离子交换与吸附理化检验. B, 化学分册力学季刊力学进展力学学报力学与实践量子电子学报量子光学学报辽宁农业科学林产工业林产化学与工业林业科学林业科学研究临床儿科杂志临床耳鼻咽喉头颈外科杂志临床放射学杂志临床检验杂志临床麻醉学杂志临床皮肤科杂志临床心血管病杂志临床与实验病理学杂志流体机械麦类作物学报煤炭学报煤炭转化煤田地质与勘探棉纺织技术棉花学报免疫学杂志模糊系统与数学模式识别与人工智能膜科学与技术摩擦学学报内蒙古大学学报. 自然科学版内燃机工程内燃机学报纳米技术与精密工程南方医科大学学报南京大学学报. 自然科学版南京航空航天大学学报南京理工大学学报. 自然科学版南京林业大学学报. 自然科学版南京农业大学学报南京气象学院学报南京师大学报. 自然科学版南京医科大学学报. 自然科学版南京中医药大学学报南开大学学报. 自然科学版南开管理评论南水北调与水利科技泥沙研究农药学学报农业工程学报农业环境科学学报农业机械学报农业经济问题农业生物技术学报农业系统科学与综合研究农业现代化研究气候变化研究进展气候与环境研究气象气象科学气象学报汽车工程强激光与粒子束桥梁建设轻金属清华大学学报. 自然科学版情报科学情报学报燃料化学学报燃烧科学与技术热带海洋学报热带气象学报热带亚热带植物学报热带作物学报热加工工艺热科学与技术热能动力工程人工晶体学报人类工效学人类学学报日用化学工业软件学报色谱山地学报山东大学学报. 工学版山东大学学报. 理学版山东大学学报. 医学版山东农业大学学报. 自然科学版山西大学学报. 自然科学版山西农业科学陕西师范大学学报. 自然科学版上海大学学报. 自然科学版上海海洋大学学报上海环境科学上海交通大学学报上海交通大学学报. 医学版上海农业学报上海天文台年刊上海医学神经解剖学杂志沈阳农业大学学报沈阳药科大学学报肾脏病与透析肾移植杂志生理科学进展生命的化学生命科学生命科学研究生态毒理学报生态环境学报生态科学生态学报生态学杂志生态与农村环境学报生物多样性生物工程学报生物化学与生物物理进展生物技术生物技术通报生物加工过程生物数学学报生物物理学报生物信息学生物学杂志生物医学工程学杂志生殖与避孕声学技术声学学报湿地科学石油地球物理勘探石油化工石油勘探与开发石油炼制与化工石油实验地质石油物探石油学报石油学报. 石油加工石油与天然气地质石油与天然气化工石油钻采工艺时间频率学报时珍国医国药实验力学实验流体力学实用妇产科杂志实用口腔医学杂志实用肿瘤杂志食品工业科技食品科学食品与发酵工业食品与生物技术学报世界地震工程世界地质世界科技研究与发展世界林业研究兽类学报数据采集与处理数理统计与管理数量经济技术经济研究数学的实践与认识数学进展数学年刊. A辑数学物理学报数学学报数学研究数学研究与评论数学杂志数值计算与计算机应用水产学报水处理技术水动力学研究与进展. A辑水科学进展水力发电学报水利水电技术水利水运工程学报水利学报水生生物学报水土保持通报水土保持学报水土保持研究水文水文地质工程地质水资源与水工程学报四川大学学报. 工程科学版四川大学学报. 医学版四川大学学报. 自然科学版四川动物四川农业大学学报四川师范大学学报. 自然科学版苏州大学学报. 医学版塑料塑料工业塑性工程学报台湾海峡太阳能学报炭素技术探测与控制学报特种铸造及有色合金天津大学学报天津医药天然产物研究与开发天然气地球科学天然气工业天然气化工. C1化学与化工天文学报天文学进展天文研究与技术铁道科学与工程学报铁道学报听力学及言语疾病杂志通信学报同济大学学报. 自然科学版涂料工业土木工程学报土木建筑与环境工程土壤土壤通报土壤学报推进技术微波学报微电子学微电子学与计算机微计算机信息微生物学报微生物学通报微生物学杂志微特电机微体古生物学报卫生研究无机材料学报无机化学学报武汉大学学报. 工学版武汉大学学报. 理学版武汉大学学报. 信息科学版武汉理工大学学报武汉植物学研究物理物理化学学报物理学报物理学进展物探化探计算技术物探与化探西安电子科技大学学报西安建筑科技大学学报西安交通大学学报. 医学版西安交通大学学报. 自然科学版西北大学学报. 自然科学版西北地震学报西北工业大学学报西北林学院学报西北农林科技大学学报. 自然科学版西北农业学报西北植物学报西南大学学报.自然科学版西南交通大学学报西南农业学报西南师范大学学报. 自然科学版西南石油大学学报.自然科学版稀土稀有金属稀有金属材料与工程稀有金属与硬质合金系统仿真学报系统工程系统工程理论与实践系统工程学报系统工程与电子技术系统管理学报系统科学与数学细胞生物学杂志细胞与分子免疫学杂志厦门大学学报. 自然科学版纤维素科学与技术现代地质现代化工现代口腔医学杂志现代雷达现代免疫学现代制造工程湘潭大学自然科学学报小型微型计算机系统心理科学心理科学进展心理学报心血管病学进展新疆地质新疆农业科学新疆石油地质新型炭材料信号处理信息工程大学学报信息与控制畜牧兽医学报循证医学压电与声光牙体牙髓牙周病学杂志亚热带植物科学岩矿测试岩石矿物学杂志岩石力学与工程学报岩石学报岩土工程学报岩土力学岩性油气藏研究与发展管理盐湖研究眼科研究扬州大学学报. 农业与生命科学版扬州大学学报. 自然科学版遥感技术与应用遥感信息遥感学报药物分析杂志药物生物技术药学服务与研究药学学报冶金分析液晶与显示医学分子生物学杂志医用生物力学仪表技术与传感器仪器仪表学报遗传印染应用泛函分析学报应用概率统计应用光学应用化工应用化学应用基础与工程科学学报应用激光应用科学学报应用力学学报应用气象学报应用生态学报应用声学应用数学应用数学和力学应用数学学报应用数学与计算数学学报应用与环境生物学报营养学报影像科学与光化学油田化学铀矿地质有机化学有色金属渔业科学进展宇航材料工艺宇航计测技术宇航学报玉米科学预测园艺学报原子核物理评论原子能科学技术原子与分子物理学报云南大学学报. 自然科学版云南农业大学学报云南植物研究运筹学学报运筹与管理杂交水稻灾害学噪声与振动控制浙江大学学报. 工学版浙江大学学报. 理学版浙江大学学报. 农业与生命科学版浙江大学学报. 医学版浙江林学院学报浙江林业科技浙江农业学报针刺研究真空科学与技术学报诊断病理学杂志振动、测试与诊断振动工程学报振动与冲击郑州大学学报. 工学版植物保护植物保护学报植物病理学报植物分类学报植物生理学通讯植物生态学报植物学通报植物研究植物遗传资源学报植物营养与肥料学报植物资源与环境学报制冷学报制造技术与机床制造业自动化质谱学报中草药中成药中风与神经疾病杂志中国癌症杂志中国安全科学学报中国表面工程中国病理生理杂志中国病原生物学杂志中国草地学报中国超声医学杂志中国当代儿科杂志中国地方病学杂志中国地震中国地质中国地质灾害与防治学报中国电机工程学报中国电力中国动脉硬化杂志中国法医学杂志中国腐蚀与防护学报中国给水排水中国工程科学中国工业经济中国公共卫生中国公路学报中国骨质疏松杂志中国管理科学中国惯性技术学报中国海洋大学学报. 自然科学版中国海洋药物中国航海中国环境监测中国环境科学中国机械工程中国激光中国激光医学杂志中国急救医学中国脊柱脊髓杂志中国寄生虫学与寄生虫病杂志中国矫形外科杂志中国介入影像与治疗学中国康复医学杂志中国抗生素杂志中国科技史杂志中国科学. A辑, 数学中国科学. B辑, 化学中国科学. C辑, 生命科学中国科学. D辑, 地球科学中国科学. E辑, 信息科学中国科学. G辑, 物理学,力学,天文学中国科学基金中国科学技术大学学报中国科学院研究生院学报中国空间科学技术中国矿业大学学报. 自然科学版中国老年学杂志中国粮油学报中国临床解剖学杂志中国临床心理学杂志中国临床药理学与治疗学中国临床药理学杂志中国媒介生物学及控制杂志中国免疫学杂志中国内镜杂志中国男科学杂志中国农村经济中国农业大学学报中国农业科学中国农业气象中国普通外科杂志中国全科医学中国人口·资源与环境中国人兽共患病学报中国软科学中国沙漠中国神经精神疾病杂志中国生化药物杂志中国生态农业学报中国生物防治中国生物工程杂志中国生物化学与分子生物学报中国生物医学工程学报中国生物制品学杂志中国石油大学学报. 自然科学版中国实验动物学报中国实用儿科杂志中国实用妇科与产科杂志中国实用内科杂志中国实用外科杂志中国实用眼科杂志中国兽医科学中国兽医学报中国输血杂志中国水产科学中国水稻科学中国塑料中国糖尿病杂志中国陶瓷中国疼痛医学杂志中国铁道科学中国图象图形学报中国土地科学中国土壤与肥料中国危重病急救医学中国微生态学杂志中国卫生统计中国稀土学报中国现代医学杂志中国现代应用药学中国心理卫生杂志中国新药与临床杂志中国新药杂志中国行为医学科学中国修复重建外科杂志中国血吸虫病防治杂志中国循环杂志中国循证儿科杂志中国循证医学杂志中国烟草科学中国烟草学报中国岩溶中国药房中国药科大学学报中国药理学通报中国药理学与毒理学杂志中国药物化学杂志中国药物依赖性杂志中国药学杂志中国野生植物资源中国医科大学学报中国医学计算机成像杂志中国医学科学院学报中国医学物理学杂志中国医学影像技术中国医学影像学杂志中国医药工业杂志中国医药生物技术中国医院药学杂志中国应用生理学杂志中国油料作物学报中国油脂中国有色金属学报中国预防兽医学报中国运动医学杂志中国造船中国造纸学报中国针灸中国中西医结合杂志中国中药杂志中国中医基础医学杂志中国肿瘤临床中国肿瘤生物治疗杂志中国组织工程研究与临床康复中国组织化学与细胞化学杂志中华病理学杂志中华超声影像学杂志中华传染病杂志中华创伤骨科杂志中华创伤杂志中华儿科杂志中华耳鼻咽喉头颈外科杂志中华耳科学杂志中华放射学杂志中华放射医学与防护杂志中华放射肿瘤学杂志中华风湿病学杂志中华妇产科杂志中华肝胆外科杂志中华肝脏病杂志中华高血压杂志中华骨科杂志中华航海医学与高气压医学杂志中华航空航天医学杂志中华核医学杂志中华护理杂志中华急诊医学杂志中华检验医学杂志中华结核和呼吸杂志中华精神科杂志中华口腔医学杂志中华劳动卫生职业病杂志中华老年医学杂志中华流行病学杂志中华麻醉学杂志中华泌尿外科杂志中华内分泌代谢杂志中华内科杂志中华男科学中华皮肤科杂志中华普通外科杂志中华器官移植杂志中华烧伤杂志中华神经科杂志中华神经外科杂志中华肾脏病杂志中华实验和临床病毒学杂志中华实验外科杂志中华手外科杂志中华外科杂志中华微生物学和免疫学杂志中华物理医学与康复杂志中华显微外科杂志中华消化外科杂志中华消化杂志中华小儿外科杂志中华心血管病杂志中华胸心血管外科杂志中华血液学杂志中华眼底病杂志中华眼科杂志中华医学遗传学杂志中华医学杂志中华医院感染学杂志中华预防医学杂志中华整形外科杂志中华中医药杂志中华肿瘤杂志中南大学学报. 医学版中南大学学报. 自然科学版中南林业科技大学学报中山大学学报. 医学科学版中山大学学报. 自然科学版中文信息学报中药材中药新药与临床药理中药药理与临床中医杂志肿瘤肿瘤防治研究种子重庆大学学报. 自然科学版重庆医科大学学报铸造资源科学自动化学报自然科学进展自然科学史研究自然灾害学报自然资源学报作物学报作物杂志
204 浏览 5 回答
128 浏览 5 回答
278 浏览 2 回答
168 浏览 2 回答
89 浏览 5 回答
96 浏览 4 回答
269 浏览 5 回答
162 浏览 4 回答
336 浏览 5 回答
347 浏览 1 回答
314 浏览 6 回答
238 浏览 3 回答
121 浏览 5 回答
236 浏览 5 回答
325 浏览 3 回答