《遥感图像处理实验论文-基于MODIS影像的西安市PM2.5研究.docx》由会员分享，可在线阅读，更多相关《遥感图像处理实验论文-基于MODIS影像的西安市PM2.5研究.docx（23页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、 遥感技术与数字图像处理实验报告作品名称： 基于MODIS影像的西安市PM2.5研究学 校： 兰州交通大学_ 指导教师： _ _杨树文_ 团队成员： _ _殷成钰_ _ 摘要可入肺颗粒物 PM2.5是影响大气环境质量的重要因子之一，不仅对人体健康造成严重危害，同时也是影响全球气候、区域城市能见度重要因子，对其进行大范围的分布特征研究。本文西安市PM2.5浓度反演实验是建立在气溶胶光学厚度基础之上，利用晴天环境下的MODIS L1B数据，在ENVI环境下，使用云检测以及气溶胶反演扩展工具，将确定的地表反射率反演气溶胶光学厚度，并分析西安市当日的气溶胶光学厚度与PM2.5的空间分布特征.目录摘要1
2、1.MODIS 数据介绍31.1数据概述31.2数据分级31.3产品类型41.4 MODIS L1B数据命名规则42.MODIS数据准备42.1搜索所需数据42.2数据下载73.MODIS数据分析与处理83.1实验技术路线83.2反射率与发射率文件处理93.2.1 发射率文件几何校正93.2.2 反射率文件几何校正113.2.3反射率和发射率数据合成与裁剪113.3角度数据处理123.3.1角度数据合成123.3.2角度数据集的几何校正143.3.3角度数据波段运算与裁剪153.4气溶胶反演163.4.1 云检测163.4.2气溶胶反演173.5反演结果处理184.结果分析与总结194.1西安
3、市PM2.5质量浓度空间特征分析19参考文献211.MODIS 数据介绍1.1数据概述MODIS 数据自 2000 年 4 月正式向全球用户发布以来，NASA 对 MODIS 数据实行全球免费接收的政策，我国自 2001 年 3 月开始接收 MODIS 数据，目前已在国内建立了数个接收站，使得 MODIS 数据的获取更加方便和廉价。以下是MODIS的各个通道、波段的主要用途：通道波段宽度主要用途1-70.459-2.135m陆地、云特征及边界8-190.105-0.965m海洋水色、浮游生物、生物化学及大气水汽20-253.660-4.549m地球表面，云及大气温度26-301.360-9.8
4、80m卷云，水汽及臭氧31-3610.780-14.385m地球表面、云温度及云顶高度1.2数据分级MODIS数据产品分级系统：MODIS标准数据产品分级系统由5级数据构成，它们分别是：0级、1级、2级、3级和4级。 0级数据：卫星地面站直接接收到的、未经处理的、包括全部数据信息在内的原始数据为0级数据。其数据文件扩展名为(.PDS)。 1级数据：对完全分辨率的、没有经过处理的仪器数据进行重建，使用辅助数据注解，进行时间配准，计算和增补后数据产品，包括 MOD01(L1A)、MOD02(L1B)和 MOD03 三类文件，其中 MOD02(L1B)是 MOD0I(L1A)和 MOD03处理之后获
5、取的数据文件，也是用户可直接使用的数据，这些数据文件主要为HDF 格式，文件的扩展名为(.HDF)。2级数据：在1级数据基础上开发出的、具有相同空间分辨率和覆盖相同地理区域的数据为2级数据。 3级数据：3级数据是以统一的时间-空间栅格表达的变量，通常具有一定的完整性和一致性。在3级水平上，将可以集中进行科学研究，如：定点时间序列，来自单一技术的观测方程和通用模型等。4级数据：通过分析模型和综合分析3级以下数据得出的结果数据为4级数据。1.3产品类型MODIS 标准数据产品根据内容总计分解为 44 种标准数据产品类型，其中，2-4级数据产品为三种主要标准数据产品类型，包括：陆地标准数据产品、大气
6、标准数据产品和海洋标准数据产品。1.4 MODIS L1B数据命名规则MOD02QKM 校正过的 250 米分辨率的对地观测数据MOD02HKM 校正过的 500 米分辨率的对地观测数据，包括25米重采样为500 米分辨率数据MOD021KM校正过的1 公里分辨率的对地观测数据，包括250米和500 米重采样为1 公里分辨率的数据。MOD02OBC 星上校正器（OBC）和工程数据。MOD02OBC MYD0 L1B数据文件名的定义规则:例如： MOD021KM.A2016036.0405.005.2016036134251.hdfMOD表示传感器为MODIS，02表示L1B数据，A表示上午星T
7、ERRA，B表示下午星Auqa,1KM表示1公里分辨率，2016036表示数据采集日期为2016年的第36天，134251表示该轨数据是在国际标准时间13时42分51秒入境的。2.MODIS数据准备本实验数据来源于美国LAADS DAAC,网址：https:/ladsweb.nascom.nasa.gov/search/2.1搜索所需数据1）打开网页，选中需要的数据源类型 。在Satellite/Instrument 一栏下拉菜单选 Terra MODIS，Group下拉菜单中选择:Terra Level 1 Products，再选择所需数据类型。注意：本实验选用分辨率为 1000m的 MOD
8、IS L1B 数据，其中，MOD03数据是用于对1KM,QKM,HKM数据进行几何纠正所用。具体如下图。2）确定数据日期。网页中日期类型为：月/日/年 时:分:秒，并且其显示的时间为UTC时间，换算为北京时间为：UTC时间=北京时间-8小时。在此选择2016年2月1日到2016年2月8日范围内的数据。如下图。3）选择在spatial selection 选项中选择“latitude/longitude”,按经纬度形式选择影像范围。 本实验研究范围为西安市，所以选择合适经纬度进行数据搜索查询，具体如下图。4）点击search查到需要的数据： 勾选需要的数据，可以通过数据右侧“+View RGB”
9、查看数据经纬度范围，如图。选择好数据后点击order files now， 输入你接收信息的邮箱，点 order开始订购该数据。5）所订购数据的存放位置信息： 点击Data-Track Orders可以查看所有已订购的数据的状态。如下图。记住所订购数据的Order ID，这是所订购数据的文件夹名。如果 state显示 available即可开始下载。2.2数据下载数据需要通过新建FTP站点下载，根据提示新建FTP站点：FTP:ladsweb.nascom.nasa.gov；username: anonymous；password:（自己申请数据的邮箱）。完成后打开orders文件夹，在里面找到
10、自己订购的数据文件夹拖到左边本地电脑的目录下即可开始下载。如下图。 3.MODIS数据分析与处理3.1实验技术路线本研究气溶胶反演主要用到分辨率为 1KM的 MODIS L1B 数据，其包括一个是发射率Emissive（band20-band36），第二个是辐射率Radiance（band1-band26）,第三个是反射率Reflectance（band1-band26）还有四个角度数据集：卫星天顶角（SensorZ），卫星方位角（SensorA），太阳天顶角（SolarZ），太阳方位角（SolarA）。分析过程具体包括几何校正、波段合成、云检测、气溶胶反演等。采用经典的暗像元法（DDV）也叫
11、浓密植被法。打开角度数据集MODIS L1B 1KM数据发射率几何校正导出GCP控制点反射率几何校正合成裁剪云检测气溶胶反演查找表角度数据合成角度数据波段运算角度数据几何校正反演结果3.2反射率与发射率文件处理3.2.1 发射率文件几何校正(1)ENVI中有专门对MODIS数据进行几何校正的工具，打开【Toolbox】-【Geometric Correction】-【Georeference by Sensor】-【Georeference MODIS】，先选择发射率文件(Emissive)，如下图：(2)在Georeference MODIS Parameters对话框中参数设置如下图：注意
12、:需要将GCP控制点文件保存下来。 (3)单击OK后，在Registration Parameters窗口中设置参数如下图所示：注意：分辨率默认是1000。(4)选择结果输出位置后，单击OK，得到发射率几何校正结果如下图：3.2.2 反射率文件几何校正过程与发射率文件几何校正相似，所以在这里不再过多赘述，参考发射率文件几何校正过程将反射率文件(Reflectance)进行几何校正。此步不需要进行再次保存GCP控制点文件。3.2.3反射率和发射率数据合成与裁剪(1)打开几何校正后的反射率与发射率数据集，选择【Raster Management】-【Layer Stacking】工具，点击【Imp
13、ort File】按钮，选择几何校正后的反射率和发射率文件，如下图：(2)单击Spatial Subset按钮进行裁剪，在弹出的窗口中选择ROI/EVF，然后选择：EVD:西安市.shp进行裁剪，如下图：(3)回到波段合成窗口中单击【Reorder File】调整文件顺序，必须反射率在上，发射率在下，如下图：(4) 选择合成裁剪结果输出位置后，单击OK，得到合成结果如下图：3.3角度数据处理3.3.1角度数据合成(1)打开角度数据集，选择【File】-【Open As】-【Generic Formats】-【HDF4】，选择MOD03(.hdf)数据文件，在列表中选择四个角数据集：卫星天顶角（
14、SensorZ），卫星方位角（SensorA），太阳天顶角（SolarZ），太阳方位角（SolarA），如下图：(2)打开【Raster Management】-【New File Builder】工具，单击【Import File】工具选择打开的四个角度数据，调整其顺序为：卫星天顶角（Sensor Z），卫星方位角（Sensor A），太阳天顶角（Solar Z），太阳方位角（Solar A）。如下图：3.3.2角度数据集的几何校正查看元数据信息，角度数据的行列数是1354*2030，与发射率的相同，因此不需要进行重采样，但需要几何校正。(1)选择【Toolbox】-【Geometric C
15、orrection】-【Registration】-【 Warp from GCPs：Image to Map Registration】工具，打开选择GCP控制点文件的对话框，在此选择之前保存的GCP控制点文件。(2)在Image to Map Registration窗口中，设置投影参数和输出分辨率(与反射率需要相同，这里为1000)，单击OK。(3)选择合成后的角度数据，在Registration Parameters对话框中设置参数，选择几何校正与重采样方法，与发射率校正结果匹配。如下图：3.3.3角度数据波段运算与裁剪(1)我们可以查看角度数据中的值，其扩大了100倍，所以要将角度数
16、据乘以0.01。打开【Tool box】-【Band Ratio】-【Band Math】工具，输入公式：float(b1)*0.01，单击OK。如图：(2)在参数面板中点击Map Variable to Input File，选择角度数据几何校正结果，如图：(3)回到参数设定面板，单击Spatial Subset按钮，在弹出的面板中单击“ROI/EVF”，选择矢量数据进行裁剪。如图：(4)设置好结果输出路径后，单击OK，得到结果如下图：3.4气溶胶反演3.4.1 云检测云检测需要用到云检测扩展工具，该工具实现对反射率和发射率的合成文件进行去云处理。(1)安装的modis_cloud工具，选择
17、反射率与发射率合成裁剪后的结果，并选择结果输出位置，如图：3.4.2气溶胶反演本实验利用气溶胶反演工具(modis_aerosol_inversion)也是扩展工具，在Extensions目录下；(1)打开【Tool box】-【Extensions】-【modis_aerosol_inversion】，先选择云检测结果，单击OK；(2)选择角度最终处理结果数据，单击OK；(3)选择查找表文件，这里使用的查找表是通用的文本文件，适合3-6月份的MODIS数据影像，也可以自己制作查找表；(4)选择结果输出路径和文件名。3.5反演结果处理(1)在Layer Manager中,右键单击反演结果图层，
18、选择Raster Color Slices，在弹出的窗口中选择反演结果波段，单击OK。(2)在Edit Raster Color Slices窗口下编辑分割参数。本实验分为七个等级。(3)最后用裁剪出西安市范围，最终结果如图：4.结果分析与总结4.1西安市PM2.5质量浓度空间特征分析根据2016 年 2 月5日的西安市各环境空气质量监测站点发布的PM2.5 质量浓度数据，计算出西安市 PM2.5 日质量浓度，如表所示：西安市主城区PM2.5质量浓度点位名称PM2.5高压开关厂67兴庆小区48纺织城小寨53市人民体育场72高新西区68经开区93长安区94阎良区129临潼区62曲江文化集团68广
19、运潭84草滩99全市平均77基于以上数据，本文首先利用克里金插值法对西安市主城区PM2.5 质量浓度做插值分析，对当日PM2.5 质量浓度空间分布简单地做出展示，并将预测分析的结果与地统计分析结果作比较。利用 ArcGIS10.2.2 的空间插值工具进行计算得到西安市主城区 PM2.5 质量浓度空间插值分析结果，如图所示：可以看出，西安市主城区PM2.5 质量浓度在空间分布上呈现多峰值的特点，其中阎良区和长安区的 PM2.5 质量浓度达到了全区域的最高值，从整体上来看，西安市主城区南部区域的 PM2.5 质量浓度明显高于北部区域，由北向南呈现先下降后上升的趋势，且在边缘地形起伏地带出现高值连片
20、区域，说明地形的起伏对城市中 PM2.5 浓度空间变化亦有重大影响。通过比较，MODIS 反演的气溶胶光学厚度和颗粒物 PM2.5的浓度具有相关性。用MODIS反演的气溶胶光学厚度可以用于估算区域PM2.5的分布特征，在时间和空间上，相对于地面监测资料具有较大的优势。参考文献1王静，杨复沫，王鼎益，贺克斌.北京市MODIS气溶胶光学厚度和PM2.5质量浓度的特征及其相关性J.中国科学院研究生学院学报，2010，27(1).2闫婧华.MISR和MODIS气溶胶光学厚度产品的数据融合及其在大气细颗粒物监测中的应用研究D.中国海洋大学，2014.3夏志业，刘志红，王永前，陈洪滨，孙明江，焦露.MOD
21、IS气溶胶光学厚度的PM2.5质量浓度遥感反演研究J.高原气象，2015.34(6).4卢亚飞.MIDIS气溶胶光学厚度与PM2. 5浓度相关性分析以北京市为例D.北京建筑大学，2012.5李亚娟.基于MIDIS的兰州市区冬季PM2. 5浓度定量估算分析D.西北师范大学，2013.6张智，白穆，游浩妍.基于MODIS数据的PM2. 5反演在大气污染监测中的应用J.测绘科学，2016，41(9).7范辰乾.基于MODIS数据的气溶胶光学厚度与PM2.5浓度关系研究以济南市为例D.山东师范大学，2014.8陈美娟.基于MODIS数据的沈阳市PM2.5反演研究D.沈阳农业大学，2016.9陈辉，厉青
22、，张玉环，周春艳，王中挺.基于地理加权模型的我国冬季PM2.5遥感估算方法研究J.环境科学学报，2016，36(6).10马宗伟.基于卫星遥感的我国PM2.5时空分布研究J.南京大学，2015.11金洪芳.基于遥感技术的大气颗粒细物PM2.5监测研究进展J.测绘与空间地理信息，2016，39(8).12杨春，石小毛，刘翔飞.江苏省MODIS气溶胶光学厚度和PM2.5质量浓度关系研究J.常州大学学报，2015，27(4).13陈辉，厉青，王中挺，毛慧琴，周春艳，张丽娟，晃雪林.利用MODIS资料监测京津冀地区近地面PM2.5方法研究J.气象与环境学报，2014，30(5).14徐建辉，江洪.长江
23、三角洲PM2.5质量浓度遥感估算与时空分布特征J.环境科学，2015，36(9).15刘欢.重庆市主城区冬季大气细颗粒物（PM2.5）质量浓度的遥感估算模型研究D.重庆师范大学，2016.1. 基于C8051F单片机直流电动机反馈控制系统的设计与研究2. 基于单片机的嵌入式Web服务器的研究 3. MOTOROLA单片机MC68HC（8）05PV8/A内嵌EEPROM的工艺和制程方法及对良率的影响研究 4. 基于模糊控制的电阻钎焊单片机温度控制系统的研制 5. 基于MCS-51系列单片机的通用控制模块的研究 6. 基于单片机实现的供暖系统最佳启停自校正（STR）调节器7. 单片机控制的二级倒立
24、摆系统的研究8. 基于增强型51系列单片机的TCP/IP协议栈的实现 9. 基于单片机的蓄电池自动监测系统 10. 基于32位嵌入式单片机系统的图像采集与处理技术的研究11. 基于单片机的作物营养诊断专家系统的研究 12. 基于单片机的交流伺服电机运动控制系统研究与开发 13. 基于单片机的泵管内壁硬度测试仪的研制 14. 基于单片机的自动找平控制系统研究 15. 基于C8051F040单片机的嵌入式系统开发 16. 基于单片机的液压动力系统状态监测仪开发 17. 模糊Smith智能控制方法的研究及其单片机实现 18. 一种基于单片机的轴快流CO,2激光器的手持控制面板的研制 19. 基于双单
25、片机冲床数控系统的研究 20. 基于CYGNAL单片机的在线间歇式浊度仪的研制 21. 基于单片机的喷油泵试验台控制器的研制 22. 基于单片机的软起动器的研究和设计 23. 基于单片机控制的高速快走丝电火花线切割机床短循环走丝方式研究 24. 基于单片机的机电产品控制系统开发 25. 基于PIC单片机的智能手机充电器 26. 基于单片机的实时内核设计及其应用研究 27. 基于单片机的远程抄表系统的设计与研究 28. 基于单片机的烟气二氧化硫浓度检测仪的研制 29. 基于微型光谱仪的单片机系统 30. 单片机系统软件构件开发的技术研究 31. 基于单片机的液体点滴速度自动检测仪的研制32. 基
26、于单片机系统的多功能温度测量仪的研制 33. 基于PIC单片机的电能采集终端的设计和应用 34. 基于单片机的光纤光栅解调仪的研制 35. 气压式线性摩擦焊机单片机控制系统的研制 36. 基于单片机的数字磁通门传感器 37. 基于单片机的旋转变压器-数字转换器的研究 38. 基于单片机的光纤Bragg光栅解调系统的研究 39. 单片机控制的便携式多功能乳腺治疗仪的研制 40. 基于C8051F020单片机的多生理信号检测仪 41. 基于单片机的电机运动控制系统设计 42. Pico专用单片机核的可测性设计研究 43. 基于MCS-51单片机的热量计 44. 基于双单片机的智能遥测微型气象站 4
27、5. MCS-51单片机构建机器人的实践研究 46. 基于单片机的轮轨力检测 47. 基于单片机的GPS定位仪的研究与实现 48. 基于单片机的电液伺服控制系统 49. 用于单片机系统的MMC卡文件系统研制 50. 基于单片机的时控和计数系统性能优化的研究 51. 基于单片机和CPLD的粗光栅位移测量系统研究 52. 单片机控制的后备式方波UPS 53. 提升高职学生单片机应用能力的探究 54. 基于单片机控制的自动低频减载装置研究 55. 基于单片机控制的水下焊接电源的研究 56. 基于单片机的多通道数据采集系统 57. 基于uPSD3234单片机的氚表面污染测量仪的研制 58. 基于单片机
28、的红外测油仪的研究 59. 96系列单片机仿真器研究与设计 60. 基于单片机的单晶金刚石刀具刃磨设备的数控改造 61. 基于单片机的温度智能控制系统的设计与实现 62. 基于MSP430单片机的电梯门机控制器的研制 63. 基于单片机的气体测漏仪的研究 64. 基于三菱M16C/6N系列单片机的CAN/USB协议转换器 65. 基于单片机和DSP的变压器油色谱在线监测技术研究 66. 基于单片机的膛壁温度报警系统设计 67. 基于AVR单片机的低压无功补偿控制器的设计 68. 基于单片机船舶电力推进电机监测系统 69. 基于单片机网络的振动信号的采集系统 70. 基于单片机的大容量数据存储技
29、术的应用研究 71. 基于单片机的叠图机研究与教学方法实践 72. 基于单片机嵌入式Web服务器技术的研究及实现 73. 基于AT89S52单片机的通用数据采集系统 74. 基于单片机的多道脉冲幅度分析仪研究 75. 机器人旋转电弧传感角焊缝跟踪单片机控制系统 76. 基于单片机的控制系统在PLC虚拟教学实验中的应用研究77. 基于单片机系统的网络通信研究与应用 78. 基于PIC16F877单片机的莫尔斯码自动译码系统设计与研究79. 基于单片机的模糊控制器在工业电阻炉上的应用研究 80. 基于双单片机冲床数控系统的研究与开发 81. 基于Cygnal单片机的C/OS-的研究82. 基于单片
30、机的一体化智能差示扫描量热仪系统研究 83. 基于TCP/IP协议的单片机与Internet互联的研究与实现 84. 变频调速液压电梯单片机控制器的研究 85. 基于单片机-免疫计数器自动换样功能的研究与实现 86. 基于单片机的倒立摆控制系统设计与实现 87. 单片机嵌入式以太网防盗报警系统 88. 基于51单片机的嵌入式Internet系统的设计与实现 89. 单片机监测系统在挤压机上的应用 90. MSP430单片机在智能水表系统上的研究与应用 91. 基于单片机的嵌入式系统中TCP/IP协议栈的实现与应用92. 单片机在高楼恒压供水系统中的应用 93. 基于ATmega16单片机的流量
31、控制器的开发 94. 基于MSP430单片机的远程抄表系统及智能网络水表的设计95. 基于MSP430单片机具有数据存储与回放功能的嵌入式电子血压计的设计 96. 基于单片机的氨分解率检测系统的研究与开发 97. 锅炉的单片机控制系统 98. 基于单片机控制的电磁振动式播种控制系统的设计 99. 基于单片机技术的WDR-01型聚氨酯导热系数测试仪的研制 100. 一种RISC结构8位单片机的设计与实现 101. 基于单片机的公寓用电智能管理系统设计 102. 基于单片机的温度测控系统在温室大棚中的设计与实现103. 基于MSP430单片机的数字化超声电源的研制 104. 基于ADC841单片机
32、的防爆软起动综合控制器的研究105. 基于单片机控制的井下低爆综合保护系统的设计 106. 基于单片机的空调器故障诊断系统的设计研究 107. 单片机实现的寻呼机编码器 108. 单片机实现的鲁棒MRACS及其在液压系统中的应用研究 109. 自适应控制的单片机实现方法及基上隅角瓦斯积聚处理中的应用研究110. 基于单片机的锅炉智能控制器的设计与研究 111. 超精密机床床身隔振的单片机主动控制 112. PIC单片机在空调中的应用 113. 单片机控制力矩加载控制系统的研究 项目论证，项目可行性研究报告，可行性研究报告，项目推广，项目研究报告，项目设计，项目建议书，项目可研报告，本文档支持完整下载，支持任意编辑！选择我们，选择成功！项目论证，项目可行性研究报告，可行性研究报告，项目推广，项目研究报告，项目设计，项目建议书，项目可研报告，本文档支持完整下载，支持任意编辑！选择我们，选择成功！单片机论文，毕业设计，毕业论文，单片机设计，硕士论文，研究生论文，单片机研究论文，单片机设计论文，优秀毕业论文，毕业论文设计，毕业过关论文，毕业设计，毕业设计说明，毕业论文，单片机论文，基于单片机论文，毕业论文终稿，毕业论文初稿，本文档支持完整下载，支持任意编辑！本文档全网独一无二，放心使用，下载这篇文档，定会成功！ 22 / 23