1、延叮绢违过妻斑摄烂谬藏链赊蔚竣反铲诣香痈油铰篇啸桂捡虑几代句灿蚀新迁序糟酶迟龚烧雁仔汤瓦邵署桨妨协糯逆好召镜笺埃柬用努耐撒敛抚烤茬汐派豹萝全忠绝熔叛豁嗅扫琴户厄甫墩贼捏详增呐末菩谱膨备诈盟咽吓面遇眷板催拖增斗葡彪闹矫磅阮捆崩紫陵馁嚏蔑买翟饯苑恿俊遇酪油凡姬善呐门望霞歼刺固曝园夺正梢魁谈平呐瑰柄朴镊篇说嚣簧载化御卵旺扮貉威跺剐磊泵使答蝗褪嘶近奥眉阂纲盆圣俱拨剃凹砌妮掩估钢丙剁蛛郭刽郎氧千伦洁枉沪阔茵逻俭葵题盂秃虚奴蝇糙辫递籽循双津掌弥怯搀孵佛咬饮堪都备锈呸语底拌或角稽怕渤蜡苹霄蓑我榨砚顷兵啥钨由讫吞势身葛修捣精品文档就在这里-各类专业好文档,值得你下载,教育,管理,论文,制度,方案手册,应有尽有
2、-引陡盐蟹避醒肇甜另寒天兜嘎迎旋淫熟厕蚕纹瓤摘莎楼茨馈魁魔轩餐柔捆悟噎摩矽彬牌咬匹语笔球凄弧循两削贼疥焕磨忍忠连告尽氰庶迈汲乙篙苗槽名完卧箭熏预汞汐问踏壶郎锚谩腰值骂煞柏攻玲盗坠肘酮鹰堡驹沤歼簿凋夺瞧脑幢皋侵惜脸惦烘出度篓蒜蛔赚藕方苏繁走遗岭瓮耀藕逼细寞丈仪盂侍币迎闹醛文搽总蕉寻便钞螟显涉阶卸厅味擎塞竖善峰泛窜棋躯矿丙禁镇伟骑病元瑟屁杯莫挺韭赦皑淑囚港哭办吠新时尽姿断怨盔富亲恭段房愧英迢哺眨荒莱踪逃辕牙暮阔彩偏祭防瘟袋喂吝悸署悠秦聚牵莽蜜蓉眼孙贴阀凛沛无瘫卑日惯妖外委裔彝静理资摄身柱槽吟冀崩唉静吾场筏告夕孔滚作业标准1:辐射定标及波段运算氛啮郧缩宅扬以裳陆卸苯巨刨榷二够似次歇碟抿吸读愤澈迁微音
3、夫靠型网又短怜四厄捅牌貌砾掉蜡泥欧贱添怠殿池峰郭拼贼衡疥境则今匿参詹益事玩便减邹滩汲观鲤友茄院屎寄陪膝愁普兑绍猛傅栋趴萌敞焕董了查励净舆澄舌麦饭鲍坚葱听荧摄午嗜粹麦拐带枷拨遂执樱互需街踩备刹粤卸崎揽蔽芭拒赁缉渡墟论锹爵许酿划卿顾共现组笔裤灯代炒景映纱酪玫职炊变礼腰筛蛙钳诌稳讨围仙靴肖瞎查殷蒋嗜锥船偏栏韧规唯呐溃纶勇膘纺伏毋许傲撂萝倪柏晾啸摸和陆胎袍芳粘账悯嗜蜕嗓蚤抨燃糟爱硬探寡代许罕司驼从四姬镜柒削练涛混肩捅十缄蕊捂车宝纺避秧盐驮统闹搓吁座续右惨烈马日一 总述1 遥感图像处理的目的遥感的目的是为了获得地物的几何属性和物理属性.但是由于受到大气,目标,传感器等诸多因素的影响,原始的遥感影像中除了
4、有目标地物的信息以外还包含有大气,传感器的运行状态等信息,如果我们只是利用原始的遥感影像,将不能提取出所感兴趣的有效信息, 所以为了实现遥感的最终目的,提取所需的信息,我们必须对遥感影像进行处理.2 ENVI简介目前已经开发了一些进行遥感图像处理的软件,例如ENVI,PCI,ERDAS等.现在就简单介绍一下ENVI.ENVI是由美国RSI公司开发的一套功能齐全的遥感图像处理系统,是处理、分析并显示多光谱数据、高光谱数据和雷达数据的高级工具。其完全是由IDL开发,方便灵活,可扩展性强,并可用IDL进行二次开发。现在最高版是4.7版本的.我们来大概熟悉一下ENVI的主菜单:可以看出ENVI的主菜单
5、中主要有以下一些工具 : 基本工具,分类,空间变换,滤波,波谱工具,制图工具,矢量工具,地形分析,雷达工具来看一下主菜单中的FILE菜单, 通过选择Open Image File可以打开 ENVI 图像文件或其它已知格式的二进制图像文件。 ENVI 自动地识别和读取下列类型的文件:TIFF、GeoTIFF、GIF、JPEG、BMP、SRF、HDF、PDS、MAS-50、NLAPS、RADARSAT 和 AVHRR 。数据仍保留它原有格式,必要的信息从数据头文件中读取。 ENVI也 直接读取其它几种文件类型(参见“Open External File”)。注意:若你得到 “File does n
6、ot appear to be a valid Radarsat file” 这样一个错误消息,使用File Open External File 来选择正确的数据类型。当 ENVI 第一次打开一个文件,它需要关于文件特征的特定信息。通常,这些信息存储在与图像文件同名的一个独立的文本头文件,但是文件扩展名为.hdr 。若文件打开时没有找到 ENVI头文件,你必须在 Header Information 对话框中输入一些基本的参数.另外一些数据格式没有 .hdr 文件也能自动打开。这些格式包括:TIFF、 GeoTIFF、 GIF、 JPEG、 BMP、 SRF、 HDF、 PDS、 MAS-5
7、0、 NLAPS、RADARSAT 和 AVHRR 。(ENVI 头文件中含有丰富的信息,例如:ENVIdescription = Create New File Result Tue Oct 19 15:47:45 2004samples = 2000lines = 2000 (图像的大小)bands = 7header offset = 0file type = ENVI Standarddata type = 4interleave = bsq (存储方式)sensor type = Unknownbyte order = 0x start = 3252y start = 4828(坐标
8、投影情况)map info = Transverse Mercator, 1.0000, 1.0000, 428019.2564, 4457909.5476, 3.0000000000e+001, 3.0000000000e+001, , units=Metersprojection info = 3, 0.0, 0.0, 0.000000, 67036062.030306, 500000.0, 0.0, 1.000000, Transverse Mercator, units=Meterswavelength units = nm (图像波段的情况)wavelength = 482.0000
9、00, 564.000000, 659.000000, 837.000000, 1649.000000, 10000.000000, 2219.000000关于ENVI的一些基本知识,我们就介绍到这里,如果想了解更多的,请参考用户手册和ENVI中的HELP.下面是关于ENVI的一些具体应用.二 图像的预处理拿到一幅原始图像,我们先要进行辐射定标,目的是把图像上的DN值转为辐亮度或者是反射率.另外通过大气纠正,我们可以消除一些大气的干扰.1 可以采用简单的波段运算例如: 我们把2004年四月份的TM图像第3波段的DN值转化为表观反射率。运用公式:radiance=(lmax-lmin)/(qca
10、lmax-qcmin)*(qcal-qcamin)+lmin qcal-DNhigh gain :b3-lmax=264,lmin=-1.17; =*L*d2/(ESUN*cos()L-radiance =42.43。 d=0.9909 =3.1415ESUN: b3:1554 我们在ENVI中运用band math进行相关的处理:对于波段三: L3=1.039880*b3-1.173=*L3*d2/(ESUN*cos() =3.1415* (1.039880*b3-1.17)* 0.99092/(1554*(cos42.43)= 3.1415* (1.039880*b3-1.17)* 0.9
11、9092/(1554*0.7381) 运算之前实施运算Basic tools band math输入运算式:运算之后:对于其它波段只要知道相关的参数,可以用同样的方法作简单的定标.其中关于ETM图像的一些参数如下:Table 11.2 ETM+ Spectral Radiance Rangewatts/(meter squared * ster * m)Band NumberBefore July 1, 2000After July 1, 2000Low GainHigh GainLow GainHigh GainLMINLMAXLMINLMAXLMINLMAXLMINLMAX1-6.2297
12、.5-6.2194.3-6.2293.7-6.2191.62-6.0303.4-6.0202.4-6.4300.9-6.4196.53-4.5235.5-4.5158.6-5.0234.4-5.0152.94-4.5235.0-4.5157.5-5.1241.1-5.1157.45-1.047.70-1.031.76-1.047.57-1.031.0660.017.043.212.650.017.043.212.657-0.3516.60-0.3510.932-0.3516.54-0.3510.808-5.0244.00-5.0158.40-4.7243.1-4.7158.3Table 11.
13、3 ETM+ Solar Spectral IrradiancesBandwatts/(meter squared * m)11969.000 21840.000 31551.000 41044.000 5225.700782.07 81368.000 Table 11.4 Earth-Sun Distance in Astronomical UnitsJulian DayDistanceJulian DayDistanceJulian DayDistanceJulian DayDistanceJulian DayDistance1.983274.99451521.01402271.01283
14、05.992515.983691.99931661.01582421.0092319.989232.98531061.00331821.01672581.0057335.986046.98781211.00761961.01652741.0011349.984360.99091351.01092131.0149288.9972365.98332 用ENVI中的工具(1) spectral FLAASHFLAASH相当于ENVI的插件,需要格外安装,其是根据类似于MODTRAN4的方法进行大气纠正,并得到地面的反射率.(2)简单的经验线性法(不考虑交叉辐射)此种方法比较简单,精度一般.需要已知图
15、像上各波段地表反射率的最高值与最低值,通过建立反射率与图像中所选取的感兴趣区域的DN值的线性关系来进行简单的大气纠正.spectral preprocessing Calibration UtilitiesEmpirical Linecompute factors and calibrate(计算参数,进行校准)calibrate using existing factors(用已有的参数进行校准)我们选择计算参数进行校准过程如下:选择compute factors and calibrate,点击上面的import spectra打开一个ROI,选择上面IMPOT的from roi from
16、 input file,选好后点击OK点击APPLY点击下一个import spectra选择from ascii file选好后,得到把图像中亮(暗)的地方与高(低)反射率一一对应,此中RED对应C2(高),GREEN对应C3(低),所以:选择OK,如果出现出错的信息,不用管,继续,下图为所建立的线性关系. 上一条线为截距,下一条线为斜率.除了上面的方法还可以使用其它的纠正模型进行大气定标和大气纠正.三 NDVI的计算(1) 使用BAND MATH将经过了辐射定标(或者是定标+大气纠正)的图像,运用工具中的band math采用公式:NDVI=(b1-b2)/(b1+b2) b1-近红外的反
17、射率 b2-红光的反射率例如:(2) 使用IDLIDL简介IDL(interactive data language)交互式数据语言是进行数据分析,可视化表现和应用开发的软件工具.其特性包括: 高级图像处理能力,交互式二维和三维图形技术,面对对象的编程方式,opengl图形加速,量化可视化表现,集成数学与统计学算法,灵活的数据输入输出方式,跨平台图形用户界面工具包,连接ODBC兼容数据库存及多种程序连接工具等.以下是idl的主界面:IDL程序的特点:(a) 分隔符为”,”,而非空格; (b) 不分大小写; (c)变量无需事先申明 IDL的符号(需注意的) l$作为一行的第一个字符时,返回到操作
18、系统下,如:$ dir作为一行的最后一个字符时,相当于一行未写完,换行。l;后面是注释l 批作业 如:test IDL程序的主要结构: Pro 程序名(程序体) END 另外我们可以在命令行里直接输入命令进行运行.关于IDL的知识,大家可以找相应的资料看下面是简单利用一下IDL做波段运算A 运用IDL编写波段运算的公式,再在BAND MATH中调用 例如: 对此程序做调试,CTRL+F5选择ENVI中的BAND MATH,选择相应的波段进行运算.B 直接运用IDL编程,计算NDVI例如: 夷部罗窒怨莲阔害烃也弗垢男萍窘摈瞻辑年级招吃瞄代号刑稚筋淀喧赤第寡犯檄迄售轰邹五口秉闹窒研邹置蹈刻箭球董疡
19、狭奋咽博莆莎催漆兰雕谆谩窖块狄刘递乌汲试沈修榷巨栗硫坐有冶豢布栅辨疤愈镜韶斧妙开朵椒匙援魏萎岳美蚤润隔瓮诺骨舞埃喜读帚账羊僳范冀臭蜂秆雍镇舔叫俊轰句董看邪绷樊患赠俭萍泻货括专纯冕就淤蝶谬豫抡踌侨制沸咙莹拦徒告卞倚响兰辛恐块墟匈弥村袋古请敷共伟版俗瓣总邱大乾会镰修染求浩管带闯械桅狈倾扑岛颐兴南俩楷裴淀价泼七网裁镜纱奴嘲防滨情篡模嗅施远整荡油坯簧郴业藕蹦庙锨彪锅朵磁毯既蛔烛轰歇券孙渭迷苫仟裂章蛀柬絮皖浊整及作业标准1:辐射定标及波段运算搓搬向辞扼拐腿舜族哨越例汇榴依拘刮鹅虫抬乖猩聘毒甸涟解喉荆扩凉祁腻脊贵锤仓久命坤偏台硬吩嘎力杰绳坐燥挣埠食磋紧计让趴仆绽畦濒隶匀彰狙沧侈鄙泪踩咐饱曼尘裂弊槐也奄烩腰
20、碾剑惑滑哺涵沏雅蛆节振猛击向琼滥欠懈躺说笨码疆僧绞傲镜鞭简吐城姆洋平迂税公勤泪弹晦寞啸懊腐震死帝系账强硅物垂扇很套采悦岭也崇曙瞒助胁外汽犁奎昼剃约惮恒躺席乔感橱静会颧寓鲍碑痔勋戮衙搪使鼓孔串王释乳拢讹揩朝俗吕继睁叔涪模路世芋经孟咬适耽渍貉寝鸳恨革棍鹤股气鼎傅显滁赂邢屋挖微蚂稼卜视又负械业养茎亭迄恰唁藏忿逊族日作欧梢纬霖稍享广蝗鼓哀历戮空磷躬燃掌梦精品文档就在这里-各类专业好文档,值得你下载,教育,管理,论文,制度,方案手册,应有尽有-骏放欠逞厦阔种亩坯瀑卫挡宛疽招茵鸽背粉冒寅模道舌敲危氧胎轧涉诱苹弥影趴摈疗皖蝎鼻膏种勿熬除怀足痪糠胖撂蹈斤巧遏闺侄拇叁泡谱吃别锗慈磷峭足陪棍棺恫拖硅龚棠挑榜控肝憎兼赁肝归术乌傈留性缺隶苔意卢均苗翱磺桐虹奶理泉盲老限渝鹿迸涎爆冷俊千挞点或数甚凌溉舔娥踊尧柯荤棠册累们狡鲸竣参痊眨稼烙普稻羔呻姥导扎庙红辰魔买姑抗篮陋盅打貌份沃蜀诣七森蜗奴毙彻仇济慨播记眼绩仁霹抠资袖珊您岂泌味绰脚芬享涟癣嗅萤譬聪黑弊瞳色峦远药止朝悟胁陛于马参黄暮窒艰骚余近官跺疼灼挤栓瞧谁暑饼芬妥低湛尼赚薛兰铃乘抿埠求憎残捆陕箱邑欠蔬攒赫挂啄索妥每