基于MODIS影像地表温度反演结果的年内时空特征研究
|
摘要
近年来,随着科学技术的发展和环境的不断恶化,人们开始越来越多的关注自己周围环境的温度变化,比如:温室效应、城市热岛等。这些例子所涉及的温度都是指气温,而气温的最直接和最主要热源为地球表面,所以,地温的研究也开始成为必要。地表温度综合了地-气相互作用过程中物质和能量交换的结果,精确地测定地表温度,对于全球气候变化的研究、灾害监测及资源管理都有非常重要的意义。
对本地接收的MODIS HDF格式原始影像数据,进行去条带、几何校正和大气校正等预处理工作。在ENVI软件支持下采用idl语言编写局地地温劈窗算法,计算每幅影像每个像元的地表温度分布。原始影像为全景影像,为了提高程序运行效率,需要进行裁剪,按照研究区的具体范围边界,给定一个研究区的掩膜,按照该掩膜的尺寸进行数据裁切。将裁剪后的结果转换成点阵图,选择出大于200K的值,其它的值为无效值。采用克里金插值法进行插值,将1km分辨率插值重采样成100m×100m,对于已经裁剪好的栅格数据做统计分析。与NASA网站上公布的同期数据产品进行计算结果的对比分析,取得了比较好的结果。计算出每天的地温数据,将每日的地温影像进行旬、月的合成,结合研究区域分析地温的时空演变规律。
根据旬地温分布图,对地温典型区进行分析,探讨城市热岛、冷湖效应、焚风效应以及地温与河流、交通线路等的关系。通过分析地温与其它自然、人文活动的关系,揭示出一些规律,为指导人类的实践活动提供便利条件。
In recent years, as the development of science and technology, and the continuous deterioration of environment, much more people began to concern their ambient temperature changes, such as: greenhouse effect, urban heat island and so on. The temperature involved in these examples is the air temperature, and the most direct and important heat source of the air temperature is the Earth's surface, so the study of the ground temperature became necessary. Land surface temperature combines the matter and energy exchange results of ground-atomoshere in the interaction.The accurate determination of the land surface temperature has a great importance for the study of global climate changes,disaster monitoring and resource management.
For the local received original image data of MODIS HDF format, take the preprocessing of destriping、geometric correction and atmospheric correction. The land surface temperature distribution was calculated for each pixel of image using split window algorithm, which was converted to grid form in ENVI software. In order to improve the running efficiency ,the whole site of original image need to be cut. According to the study area boundary, cut the data according to the size of the study area of the mask. The results were converted into bitmap, then select the value greater than 200K, the other value is invalid. Using Kriging method,resample 1km resolution into 100m×100m for statistical analysis. Compared to the same period of NASA data products published on the website, the results of in this paper achieves good results. Calculate the daily land temperature data, then compose the product of a period of ten days and month, Take the synthesis of the spatial and temporal evolution of land surface temperature combined with the study area.
Base on the temperature distribution, take an analysis of the typical area of land temperature,dicuss the urban heat island,cold lake effect,Foehn wind ,and the relation with the river、transportation lines and so on. By analyzing the land temperature with other natural,human activities, reveal some rules, and facilitate the activities of human’s practice as a guide.
对本地接收的MODIS HDF格式原始影像数据,进行去条带、几何校正和大气校正等预处理工作。在ENVI软件支持下采用idl语言编写局地地温劈窗算法,计算每幅影像每个像元的地表温度分布。原始影像为全景影像,为了提高程序运行效率,需要进行裁剪,按照研究区的具体范围边界,给定一个研究区的掩膜,按照该掩膜的尺寸进行数据裁切。将裁剪后的结果转换成点阵图,选择出大于200K的值,其它的值为无效值。采用克里金插值法进行插值,将1km分辨率插值重采样成100m×100m,对于已经裁剪好的栅格数据做统计分析。与NASA网站上公布的同期数据产品进行计算结果的对比分析,取得了比较好的结果。计算出每天的地温数据,将每日的地温影像进行旬、月的合成,结合研究区域分析地温的时空演变规律。
根据旬地温分布图,对地温典型区进行分析,探讨城市热岛、冷湖效应、焚风效应以及地温与河流、交通线路等的关系。通过分析地温与其它自然、人文活动的关系,揭示出一些规律,为指导人类的实践活动提供便利条件。
In recent years, as the development of science and technology, and the continuous deterioration of environment, much more people began to concern their ambient temperature changes, such as: greenhouse effect, urban heat island and so on. The temperature involved in these examples is the air temperature, and the most direct and important heat source of the air temperature is the Earth's surface, so the study of the ground temperature became necessary. Land surface temperature combines the matter and energy exchange results of ground-atomoshere in the interaction.The accurate determination of the land surface temperature has a great importance for the study of global climate changes,disaster monitoring and resource management.
For the local received original image data of MODIS HDF format, take the preprocessing of destriping、geometric correction and atmospheric correction. The land surface temperature distribution was calculated for each pixel of image using split window algorithm, which was converted to grid form in ENVI software. In order to improve the running efficiency ,the whole site of original image need to be cut. According to the study area boundary, cut the data according to the size of the study area of the mask. The results were converted into bitmap, then select the value greater than 200K, the other value is invalid. Using Kriging method,resample 1km resolution into 100m×100m for statistical analysis. Compared to the same period of NASA data products published on the website, the results of in this paper achieves good results. Calculate the daily land temperature data, then compose the product of a period of ten days and month, Take the synthesis of the spatial and temporal evolution of land surface temperature combined with the study area.
Base on the temperature distribution, take an analysis of the typical area of land temperature,dicuss the urban heat island,cold lake effect,Foehn wind ,and the relation with the river、transportation lines and so on. By analyzing the land temperature with other natural,human activities, reveal some rules, and facilitate the activities of human’s practice as a guide.
引文
[1]覃志豪,高懋芳,秦晓敏等.农业旱灾监测中的地表温度遥感反演方法——以MODIS数据为例[J].自然灾害学报.2005,14(4):64-71.
[2]王建芳.遥感地表温度反演在寒害监测预警中的应用[D].广州:中国科学院研究生院(广州地球化学研究所),2006.
[3]韩秀珍.东亚飞蝗灾害的遥感监测机理与方法研究[D].北京:中国科学院研究生院(遥感应用研究所),2003.
[4]南鹏,秦其明,姚云军。地表温度异常机理及热红外遥感监测——以河南焦作煤层气富集区为例[J],北京大学学报,2008-03-013
[5]陈云浩,王洁,李晓兵.夏季城市热场的卫星遥感分析[ J].国土资源遥感, 2002, (4): 55—59.
[6]Sobrino JA, G mezM, Jiménez-Mu oz J C,et al.A Simple Algorithm to Estimate Evapotranspiration from DAIS Data:Application to the DAISEX Campaigns [J].Journal of Hydrology,2005,315: 117—125.
[7]李建.长白山地区地表温度反演研究[D].长春:吉林大学,2005.
[8]李粉玲.高山高原地区地表温度反演研究[D].兰州:兰州大学,2006.
[9]叶智威覃志豪宫辉力.洪泽湖区的Landsat TM6地表温度遥感反演和空间差异分析[J].首都师范大学学报,2009,01(19):91-98.
[10]李琴,陈曦,包安明等.干旱区/半干旱区MODIS地表温度反演与验证研究[J].遥感技术与应用,2008,23(8):643-647.
[11]宫阿都,陈云浩,李京等.北京城市热岛与土地利用/土地覆被变化研究[J].中国图像图形学报,2007,12(8):1476-1482.
[12]王镇.福州城市扩展与城市热岛效益的遥感动态监测与分析[D].福州市:中国科学院上海冶金研究所,2000.
[13]刘磊.MODIS数据地表温度反演及其在长江三角洲都市圈城市热岛效应研究中的应用[D].南京:南京信息工程大学,2007.
[14]张金区.珠江三角洲地区地表热环境的遥感探测与时空演变研究[D].广州:中国科学院研究生院(广州地球化学研究所),2006.
[15]QinZ, Karnieli A. Progress in the remote sensing of land surface temperature and ground emissivity using NOAA-AVHRR data[J].Int.J.Remote Sens.1999,20:2367 -393.
[16]Qin Z, Dall OlmoG, KarnieliA et al. Derivation of split window algorithm and its sensitivity analysis for retrieving land surface temperature from NOAA-AVHRR data [J].Journal of Geophysical Research,2001,106(D19):22655-22670.
[17] Qin Z, Xu B, Zhang W et al. Comparison of split window algorithms for land surface temperature retrieval from NOAA-AVHRR data [A].IEEE 2004 International Geosciences and Remote Sensing Symposium [C],VI,3740-3743, September20-24, 2004, Anchorage,Alaska,USA.
[18]毛克彪,覃志豪.用MODIS影像反演环渤海地区的大气水汽含量[J].遥感信息. 2004. 4:47-49.
[19]毛克彪,覃志豪,王建明等.针对MODIS数据的大气水汽含量及31和32波段透过率计算[J].土资源遥感,2005,1:26-30.
[20]毛克彪,覃志豪,施建成等.针对MODIS影像的劈窗算法研究[J].武汉大学学报. 2005,30(8):703-707.
[21]丁莉东,覃志豪,毛克彪.基于MODIS影像数据的劈窗算法研究及其参数确定[J].遥感技术应用.2005,20(2):284-289.
[22]Kaufman Y J,Gao B C.Remote Sensing of Water Vapor in the Near IR from EOS/MODIS[J].IEEE Transactions on Geoscience and Remote Sensing,1992,30 (5): 871-884.
[23]高懋芳,覃志豪,徐斌.用MODIS数据反演地表温度的基本参数估计方法[J].干旱区研究.2007.24(1):113-119.
[24]Prabhakara C, Dalu G, Kunde V G. Estimation ofsea temperature from remote sensing in the 11 to 12μm window region [J].J.Geophys.Res. 1974(79):5039-5044.
[25]McMillin L M.Estimation of sea surface temperature from two infrared window measurements with different absorption[J].Journal of Geophysical Research, 1975, 20:5113-5117.
[26]Price J C. Land surface temperature measurements from split window channels of the NOAA-7 advance very high resolution radiometer[J].J Geophys Res,1984,89:7231-7237.
[27]Price J C. Land surface temperature measurements from split window channels of the NOAA-7 AVHRR[J].Geoghys.Res.,1984,79:5039-5044.
[28]Becker F.The impact of spectral emissivity on the measurement of land surface temperature from a satellite[J].Int J Remote Sens,1987,10:1509-1522..
[29]Becker F,Li Z L.Towards a local split window method over land surface[J].Int J Remote Sens.1990,11:369-393.
[30]Sobrino J A,Coll C,Caselles V.Atmospheric corrections for landsurface temperature using AVHRR channel 4 and 5[J].Remote Sens Environ,1991,38:19-34.
[31]Coll C,Caselles V,Sobrino A. Valor E. On the atmospheric dependence of the split-window equation for land surface temperature,Int.J.Remote Sens, 1994, 27:105-122.
[32]Franqa G. B.,Cracknell A. P..Retrieval of land and sea surface temperature using NOAA-11 AVHRR data in northeastern Brazil,Int.J.Remote Sens.1994,15:1695-1712.
[33]Wan Z,Dozier J.A generalized split-window algorithm for retrieving land- surface temperature from space[J].IEEE Trans Geosci Remote Sens,1996,34(4):892-905.
[34]Wan Z M,Li Z L.A Physics-Based Algorithm for Retrieving land-surface emissivisity and temperature from EOS/MODIS data ,IEEE Trans.Geosci.Remote Sens. 1997,35:980-996.
[35]Wan Z M,Zhang Y L,Zhang Q C et al., Validation of the land-surface temperature products retrieved from Terra Moderate Resolution Imaing Spectroradiometer data, Remote Sens.Environ.2002,83:163-180.
[36]Wan Z,Zhang Y,Zhang Q et al., Qulity assessment and validation of the MODIS global land surface temperature.Int.J.Remote Sens.2004,25:261-274.
[37]李召良.地表温度反演算法,先进对地观测技术与应用学术研讨会材料汇编(一),1999.
[38]马伟强,马耀明.西北干旱区地表能量初步分析[J].干旱区研究,2006,23(1):76- 82.
[39]吴锦奎,丁永建,魏智等.干旱区天然低湿牧草地参考作物蒸散量研究——以黑河中游为例[J].干旱区研究,2005,22(4):514-519.
[40]市地方志编纂委员会,石家庄市志第一卷,中国社会出版社[B],1995.12.
[41]袁金国,遥感图像数字处理,中国环境科学出版社[B],2006,2.
[42]Becker F,Li Z L.Surface temperature and emissivity at various scales: Definition,measurements and related problems[J].Remote Sens Reviews,1995,12: 225-253.
[43]张旭,刘新春,肖继东等.EOS/MODIS影像处理在塔里木河下游植被监测中的应用[J].干旱区研究,2005,22(4):532-536.
[44]毛克彪,覃志豪. MODIS影像反演环渤海地区的大气水汽含量[J].遥感信息,2004, 76(4).
[45]Kaufman Y J, and Gao B C. Remote sensing ofwater vapor in the near IR from EOS/MODIS[J].IEEE Transactions on Geosciences and Remote Sensing, 30(5): 871-884.
[46]牛生丽,唐军武,郭茂华等.基于均衡化曲线的条带消除研究[J].遥感信息.2007,3:44-48.
[47]郭广猛.关于MODIS卫星数据的几何校正方法[J].遥感信息应用技术. 2002,3:26-28.
[48]元雪勇,田庆久.光学遥感大气校正进展[J].国土资源遥感.2005,(4):1-6.
[49]宋晓宇,王纪华,刘良云等.基于高光谱遥感影像的大气校正:用AVIRIS数据评价大气校正模块FLAASH[J].遥感技术与应用.2005,20(4):393-398.
[50]Research System Inc. FLAASH User's Guide.ENVI FLAASH Version 1.0,2001:8-40.
[51]肖荣波等:城市热岛时空特征及其影响因素.气象科学[J],2007, 27(2): 230-236.
[52]窦建奇.关于城市“热岛效应”的思考.武汉城市建设学院学报[J], 2001,18(3/4):76-78.
[53]赵英时等.遥感应用分析原理与方法[M].北京∶科学出版社,2003.
[54]连志鸾,高连山,赵彦厂,等.石家庄夏季持续高温灾害的气候特征及形成机理[J].中国农业气象,2008,29(4):387-391.
[55]王红,周大迈,张爱军,张月辰,任士魁.河流故道区砂质土壤水分运移规律及灌溉制度初探[J].中国生态农业学报,2006,14(2):177-179.
[2]王建芳.遥感地表温度反演在寒害监测预警中的应用[D].广州:中国科学院研究生院(广州地球化学研究所),2006.
[3]韩秀珍.东亚飞蝗灾害的遥感监测机理与方法研究[D].北京:中国科学院研究生院(遥感应用研究所),2003.
[4]南鹏,秦其明,姚云军。地表温度异常机理及热红外遥感监测——以河南焦作煤层气富集区为例[J],北京大学学报,2008-03-013
[5]陈云浩,王洁,李晓兵.夏季城市热场的卫星遥感分析[ J].国土资源遥感, 2002, (4): 55—59.
[6]Sobrino JA, G mezM, Jiménez-Mu oz J C,et al.A Simple Algorithm to Estimate Evapotranspiration from DAIS Data:Application to the DAISEX Campaigns [J].Journal of Hydrology,2005,315: 117—125.
[7]李建.长白山地区地表温度反演研究[D].长春:吉林大学,2005.
[8]李粉玲.高山高原地区地表温度反演研究[D].兰州:兰州大学,2006.
[9]叶智威覃志豪宫辉力.洪泽湖区的Landsat TM6地表温度遥感反演和空间差异分析[J].首都师范大学学报,2009,01(19):91-98.
[10]李琴,陈曦,包安明等.干旱区/半干旱区MODIS地表温度反演与验证研究[J].遥感技术与应用,2008,23(8):643-647.
[11]宫阿都,陈云浩,李京等.北京城市热岛与土地利用/土地覆被变化研究[J].中国图像图形学报,2007,12(8):1476-1482.
[12]王镇.福州城市扩展与城市热岛效益的遥感动态监测与分析[D].福州市:中国科学院上海冶金研究所,2000.
[13]刘磊.MODIS数据地表温度反演及其在长江三角洲都市圈城市热岛效应研究中的应用[D].南京:南京信息工程大学,2007.
[14]张金区.珠江三角洲地区地表热环境的遥感探测与时空演变研究[D].广州:中国科学院研究生院(广州地球化学研究所),2006.
[15]QinZ, Karnieli A. Progress in the remote sensing of land surface temperature and ground emissivity using NOAA-AVHRR data[J].Int.J.Remote Sens.1999,20:2367 -393.
[16]Qin Z, Dall OlmoG, KarnieliA et al. Derivation of split window algorithm and its sensitivity analysis for retrieving land surface temperature from NOAA-AVHRR data [J].Journal of Geophysical Research,2001,106(D19):22655-22670.
[17] Qin Z, Xu B, Zhang W et al. Comparison of split window algorithms for land surface temperature retrieval from NOAA-AVHRR data [A].IEEE 2004 International Geosciences and Remote Sensing Symposium [C],VI,3740-3743, September20-24, 2004, Anchorage,Alaska,USA.
[18]毛克彪,覃志豪.用MODIS影像反演环渤海地区的大气水汽含量[J].遥感信息. 2004. 4:47-49.
[19]毛克彪,覃志豪,王建明等.针对MODIS数据的大气水汽含量及31和32波段透过率计算[J].土资源遥感,2005,1:26-30.
[20]毛克彪,覃志豪,施建成等.针对MODIS影像的劈窗算法研究[J].武汉大学学报. 2005,30(8):703-707.
[21]丁莉东,覃志豪,毛克彪.基于MODIS影像数据的劈窗算法研究及其参数确定[J].遥感技术应用.2005,20(2):284-289.
[22]Kaufman Y J,Gao B C.Remote Sensing of Water Vapor in the Near IR from EOS/MODIS[J].IEEE Transactions on Geoscience and Remote Sensing,1992,30 (5): 871-884.
[23]高懋芳,覃志豪,徐斌.用MODIS数据反演地表温度的基本参数估计方法[J].干旱区研究.2007.24(1):113-119.
[24]Prabhakara C, Dalu G, Kunde V G. Estimation ofsea temperature from remote sensing in the 11 to 12μm window region [J].J.Geophys.Res. 1974(79):5039-5044.
[25]McMillin L M.Estimation of sea surface temperature from two infrared window measurements with different absorption[J].Journal of Geophysical Research, 1975, 20:5113-5117.
[26]Price J C. Land surface temperature measurements from split window channels of the NOAA-7 advance very high resolution radiometer[J].J Geophys Res,1984,89:7231-7237.
[27]Price J C. Land surface temperature measurements from split window channels of the NOAA-7 AVHRR[J].Geoghys.Res.,1984,79:5039-5044.
[28]Becker F.The impact of spectral emissivity on the measurement of land surface temperature from a satellite[J].Int J Remote Sens,1987,10:1509-1522..
[29]Becker F,Li Z L.Towards a local split window method over land surface[J].Int J Remote Sens.1990,11:369-393.
[30]Sobrino J A,Coll C,Caselles V.Atmospheric corrections for landsurface temperature using AVHRR channel 4 and 5[J].Remote Sens Environ,1991,38:19-34.
[31]Coll C,Caselles V,Sobrino A. Valor E. On the atmospheric dependence of the split-window equation for land surface temperature,Int.J.Remote Sens, 1994, 27:105-122.
[32]Franqa G. B.,Cracknell A. P..Retrieval of land and sea surface temperature using NOAA-11 AVHRR data in northeastern Brazil,Int.J.Remote Sens.1994,15:1695-1712.
[33]Wan Z,Dozier J.A generalized split-window algorithm for retrieving land- surface temperature from space[J].IEEE Trans Geosci Remote Sens,1996,34(4):892-905.
[34]Wan Z M,Li Z L.A Physics-Based Algorithm for Retrieving land-surface emissivisity and temperature from EOS/MODIS data ,IEEE Trans.Geosci.Remote Sens. 1997,35:980-996.
[35]Wan Z M,Zhang Y L,Zhang Q C et al., Validation of the land-surface temperature products retrieved from Terra Moderate Resolution Imaing Spectroradiometer data, Remote Sens.Environ.2002,83:163-180.
[36]Wan Z,Zhang Y,Zhang Q et al., Qulity assessment and validation of the MODIS global land surface temperature.Int.J.Remote Sens.2004,25:261-274.
[37]李召良.地表温度反演算法,先进对地观测技术与应用学术研讨会材料汇编(一),1999.
[38]马伟强,马耀明.西北干旱区地表能量初步分析[J].干旱区研究,2006,23(1):76- 82.
[39]吴锦奎,丁永建,魏智等.干旱区天然低湿牧草地参考作物蒸散量研究——以黑河中游为例[J].干旱区研究,2005,22(4):514-519.
[40]市地方志编纂委员会,石家庄市志第一卷,中国社会出版社[B],1995.12.
[41]袁金国,遥感图像数字处理,中国环境科学出版社[B],2006,2.
[42]Becker F,Li Z L.Surface temperature and emissivity at various scales: Definition,measurements and related problems[J].Remote Sens Reviews,1995,12: 225-253.
[43]张旭,刘新春,肖继东等.EOS/MODIS影像处理在塔里木河下游植被监测中的应用[J].干旱区研究,2005,22(4):532-536.
[44]毛克彪,覃志豪. MODIS影像反演环渤海地区的大气水汽含量[J].遥感信息,2004, 76(4).
[45]Kaufman Y J, and Gao B C. Remote sensing ofwater vapor in the near IR from EOS/MODIS[J].IEEE Transactions on Geosciences and Remote Sensing, 30(5): 871-884.
[46]牛生丽,唐军武,郭茂华等.基于均衡化曲线的条带消除研究[J].遥感信息.2007,3:44-48.
[47]郭广猛.关于MODIS卫星数据的几何校正方法[J].遥感信息应用技术. 2002,3:26-28.
[48]元雪勇,田庆久.光学遥感大气校正进展[J].国土资源遥感.2005,(4):1-6.
[49]宋晓宇,王纪华,刘良云等.基于高光谱遥感影像的大气校正:用AVIRIS数据评价大气校正模块FLAASH[J].遥感技术与应用.2005,20(4):393-398.
[50]Research System Inc. FLAASH User's Guide.ENVI FLAASH Version 1.0,2001:8-40.
[51]肖荣波等:城市热岛时空特征及其影响因素.气象科学[J],2007, 27(2): 230-236.
[52]窦建奇.关于城市“热岛效应”的思考.武汉城市建设学院学报[J], 2001,18(3/4):76-78.
[53]赵英时等.遥感应用分析原理与方法[M].北京∶科学出版社,2003.
[54]连志鸾,高连山,赵彦厂,等.石家庄夏季持续高温灾害的气候特征及形成机理[J].中国农业气象,2008,29(4):387-391.
[55]王红,周大迈,张爱军,张月辰,任士魁.河流故道区砂质土壤水分运移规律及灌溉制度初探[J].中国生态农业学报,2006,14(2):177-179.