-
摘要: 目的 在地表温度反演中,Jiménez Muoz和Sobrino开发的单通道(SC)算法因其需要的实时大气参数少而被广泛应用。由于SC算法在2003年提出时只提供了针对LandsatTM的大气参数,导致许多后续基于ETM+数据的地表温度计算也都采用TM的大气参数。即使SC算法于2009年提出了改进版,但这一混用现象仍未改观。因此,基于实测地表温度,以Landsat5/7/8号卫星的热红外数据为例比较了2003和2009版的算法,探讨了Planck函数、λ取值和大气参数等常见问题。结果表明,2009版的算法明显提高了ETM+地表温度反演的精度,且以采用Planck函数和USGS提供的λ值时所获得的平均精度最高。在当前,如果要用SC算法来反演Landsat8的地表温度,可在大气水汽含量较低时选用TIRS10波段来单独进行。
-
地形是指地球表面的高低起伏形态,高程是描述地表起伏形态最基本的几何量[1]。地貌演化、滑坡预警、洪水淹没、城市规划等诸多研究和应用领域,都依赖于地表高程作为基础数据支撑。数字高程模型(digital elevation model,DEM)是高程的数字化表达,通过一组有序数值阵列实现对地面地形的数字化模拟,其反映的是不包含自然及人工地物在内的“裸地”高程信息[2]。而数字表面模型(digital surface model,DSM)则是包含了地表建筑物、树木等自然及人工地物高度的模型,DSM在DEM的基础上添加了地表地物的高程信息[2]。DSM常常应用于土地利用覆盖分类[3]、建筑物高度估算[4]以及森林高度估算等研究[5]。
生产DEM产品的方法主要包括使野外地面测量、地形图数字化、立体摄影测量、合成孔径雷达干涉测量以及三维激光扫描等[1,6]。目前国际上已产生了多种免费开放的DEM产品,如美国航空航天局的航天飞机雷达地形任务(shuttle radar topography mission,SRTM)DEM、日本宇宙航空研究开发机构使用先进陆地观测卫星(advance land observing satellite,ALOS)上的全色立体遥感成像仪采集的图像生成的AW3D30、美国航空航天局的先进星载热辐射和反射辐射计(advanced spaceborne thermal emission and reflection radiometer,ASTER)GDEM,以及德国宇航中心发布的TanDEM-X的DEM产品[7]。ASTER和AW3D30产品都是根据光学立体摄影测量技术生产的,而SRTM和TanDEM-X的DEM产品则是应用了合成孔径雷达干涉测量技术[2]。由于测量原理的限制,AW3D30、ASTER、SRTM以及TanDEM-X等产品其实都是包含了树木和建筑物等地表物体的信息,所以这些产品更精确地是表达了地表物体的高程而非地形表面的高程(即DSM而非DEM),但是它们却常常被简化等同于DEM使用。亦有学者对已有的DEM产品进行了误差改正,例如日本Yamazaki等学者提出的去除绝对偏差、条带噪声、斑点噪声以及树高偏差后得到改进的DEM产品(multi-error-removed improved-terrain DEM,MERIT DEM)[8]。以上5种产品都具有全球或者接近全球可用性,被广泛地应用于全球变化[9]、水文模拟[10]、地质与灾害评估[11]、工程勘测[12]等各个领域。
随着开放DEM产品的蓬勃发展,各个产品的精度情况、误差分布以及适用条件越来越受学界关注。由于不同DEM产品包含不同程度的高程误差,因此这些产品作为参数被纳入到后续应用模型时会造成模拟分析的结果误差。例如,文献[13]在评估海平面上升对沿海城市淹没风险的研究中,发现使用SRTM、MERIT、TanDEM-X的产品分别进行淹没模拟时,所得到的越南风险区域面积相差50%以上。文献[14]利用ASTER和SRTM的产品对印度尼西亚进行海啸淹没模拟研究,发现使用ASTER和SRTM开发的海啸淹没模型会低估淹没程度,不宜用于水动力模型评估海啸淹没区。赵冠华等学者在研究气温数据的空间插值时,发现温度插值的误差受到DEM产品高程误差的影响,在中国西南地区,基于SRTM产品的温度插值结果优于基于ASTER产品的温度插值结果[15]。
目前已有的开放DEM产品精度评估研究,主要依靠地面高程测量、星载激光测高以及机载激光扫描3种不同的参考数据。其中,地面高程测量是最为精确的获取参考点高程的方法,常常被用来评价DEM产品的精度。例如文献[16]使用GPS测量作为参考数据对TanDEM-X的DEM产品进行了精度评估。文献[17]使用差分GPS获取高程点的坐标,评估了SRTM产品在美国密西根东南部地区的数据精度,发现在该区域总体标准偏差为3.3 m,精度难以满足裸露地表的误差标准[17]。星载激光测高仪由于具有很高的测高精度,获取的地面高程也常常用来作为DEM产品精度评价的参考数据。例如文献[18]选择中国典型高海拔山地和低海拔沿海平原为实验区,参考冰、云和陆地高程卫星的地球科学激光测高系统(geoscience laser altimeter system,GLAS)测高数据评价了SRTM和ASTER产品的精度,发现两种产品在西南山地的精度远低于在平原区域的精度。文献[19]参考GLAS测高数据,评价了资源三号DSM产品精度,发现其均方根误差在5 m以内。文献[20]基于GLAS测高数据,评价了TanDEM-X 90 m的DEM产品在中国大陆地区的误差空间分布特征,发现该产品在中国区域整体的绝对误差均值为3.89 m,中误差为9.03 m,且误差随着坡度的上升逐渐增大[20]。
参考地面高程测量和星载激光测高数据进行精度评价时,都是点对点的评价,而使用机载激光雷达(light detection and ranging,LiDAR)来评价DEM产品精度时,能够获取面状的数据,有利于分析误差的空间分布结构。LiDAR是近数十年来摄影测量与遥感领域极具革命性的技术之一。它通过主动发射激光脉冲获得地面高密度的三维点云,之后通过多次回波滤波、分类与建模可以同时得到区域的DSM和DEM[1]。机载LiDAR技术突破了单点测量,变为连续自动数据获取,提高了观测速度与效率[21]。文献[22]使用机载LiDAR数据生成的DSM作为参考,研究发现与AW3D30、ASTER和SRTM产品相比,TanDEM-X产品的高程精度更高[22]。文献[23]使用墨西哥局部区域的机载激光扫描数据生成的参考DSM,评价了多种开放DEM产品在表征地表物体高程时的精度,结果表明AW3D30产品的精度优于SRTM和ASTER的产品。文献[24]基于机载LiDAR数据,评价了高寒山区的DEM产品精度,发现SRTM产品的均方根误差(root mean square error,RMSE)为7.62 m,略微优于RMSE为9.90 m的ASTER产品。文献[25]以位于捷克共和国的Krkonose山脉等为研究对象,分别参考机载LiDAR生成的DEM和DSM产品,评价了TanDEM-X的DEM产品精度,研究发现,以DEM产品作为参考时,TanDEM-X呈现出显著高估,在森林和非森林之间的区域平均误差差异很大,而以DSM产品为参考时,则表现为微弱低估。
在已有研究中,分析不同开放DEM产品分别表征地形表面高程以及地表物体高程时的精度优劣、误差空间分布结构的研究还很缺少。因此,本文旨在研究以机载LiDAR获取的高精度DEM和DSM为参考,对ASTER、AW3D30、MERIT、SRTM和TanDEM-X这5种开放全球DEM产品进行高程精度评价与误差分析。
1 研究区与数据
1.1 研究区域
选取美国国家生态观测网络(National Ecological Observatory Network,NEON)的史密森尼学会生物保护研究所(Smithsonian Conservation Biology Institute,SCBI)实验站为研究区(如图1所示)。SCBI实验站位于美国弗吉尼亚州蓝岭山脉山脚下(78.06ºW—78.16ºW,38.08ºN—38.93ºN),气候温热湿润,年平均气温为11.6 °C,年均降水量1 125 mm。区域内地面高程起伏在251~571 m,土地覆盖类型主要包括森林、居民地,以及少量的农田和草坪。森林植被以落叶阔叶林为主,平均冠层高度约35 m,主要树种为橡树、山核桃、白蜡木和郁金香杨树。
1.2 数据介绍
1.2.1 机载LiDAR数据
研究区域的机载 LiDAR数据可以从NEON官网数据中心下载得到。该数据由NEON的航空观测平台搭载RIEGL的LMS Q780激光雷达系统于2019年6月采集,飞行高度约1 000 m,激光发散角0.25 mrad,点云密度约16.5点/m2。本文直接使用NEON数据中心对点云数据进行处理后提供的DEM与DSM产品[26]。根据数据质量控制报告,其空间分辨率为1 m,垂直精度约0.15 m,水平定位精度约0.3 m,坐标系统是WGS_1984_UTM投影坐标系,垂直基准为NAVD88。本文根据DSM减去DEM生成归一化DSM(normalized DSM,nDSM)如图2所示,结合图1实验区的遥感影像可知,区域西南部和东北部主要由高度较高的森林冠层覆盖(约20~40 m),西北角可见高度较矮的居民住宅区(约0~10 m),中部和东南角可见零星分布的住宅、农田和草坪。实验区的坡度起伏如图3所示,可以观察到大部分区域的坡度在20°以下。
1.2.2 开放DEM产品数据
选取ASTER、AW3D30、MERIT、SRTM和TanDEM-X共5种开放DEM产品进行精度比较和误差分析。其中,ASTER的DEM产品是根据ASTER获得的立体像对生成,空间分辨率为30 m,垂直方向精度为13 m,数据产品可从https://lpdaac.usgs.gov/products/astgtmv003/下载。AW3D30是由日本宇宙航空研究开发机构提供的,可以从http://www.eorc.jaxa.jp/ALOS/en/aw3d30/data/index.htm下载得到,空间分辨率为 30 m,高程精度为5 m[27]。MERIT产品去除了绝对偏差、条带噪声、斑点噪声,并进行了树高偏差改正,可从http://hydro.iis.u-tokyo.ac.jp/~yamadai/MERIT_DEM下载,其空间分辨率为90 m,高程精度为5 m。SRTM的DEM产品版本为SRTMGL1v003,可以从美国地质勘探局网站下载(https://lpdaac.usgs.gov/products/srtmgl1v003/),其空间分辨率为30 m。TanDEM-X是德国航宇中心的雷达地球观测任务,其90 m分辨率的DEM产品可从https://download.geoservice.dlr.de/TDM90/网址下载。各产品的参数汇总[28]如表1所示。
表 1 5种全球开放DEM产品具体参数Table 1. Specific Parameters of 5 Global Open DEM ProductsDEM产品 水平基准 垂直基准 空间分辨率/m 高程精度/m 版本 数据获取时间 ASTER GCS_WGS_1984 EGM96 30 13 v003 1999—2008年 AW3D30 GCS_WGS_1984 EGM96 30 5 v3.1 2006—2011年 MERIT GCS_WGS_1984 EGM96 90 5 v1.0.3 1987—2017年 SRTM GCS_WGS_1984 EGM96 30 10 GL1 v003 1999—2000年 TanDEM-X GCS_WGS_1984 WGS84 90 10 v01 2011—2015年 2 研究方法
2.1 坐标转换、重采样及区域裁剪
实验中,研究参考的机载LiDAR生成的DEM和DSM产品的平面基准是WGS_1984_UTM_17N投影坐标系,垂直基准是NAVD88系统,空间分辨率为1 m,而5种开放DEM产品的具体参数如表1所示,因此,后续分析中,统一选择WGS_1984_UTM_17N为水平基准,EGM96为垂直基准,空间分辨率为1 m。首先,使用GDAL命令行工具将机载LiDAR生成的DEM和DSM产品以及TanDEM-X的DEM产品的垂直基准转换为EGM96系统。GDAL是一个开源栅格空间数据转换库,是诸如QGIS等地理信息系统软件的核心组成之一,广泛应用于空间数据处理。然后,使用ArcGIS软件统一将ASTER、AW3D30、MERIT、SRTM以及TanDEM-X 5种产品从地理坐标系GCS_WGS_1984转换到平面投影坐标系WGS_1984_UTM_17N,空间分辨率统一为1 m。由于实验使用的5种开放DEM产品数据所囊括的区域范围都远大于机载LiDAR数据覆盖的范围,所以以机载LiDAR数据的范围为输入,按掩膜提取的方法对DEM产品影像数据进行了区域裁剪。
2.2 误差分析及精度评价
对开放DEM产品的误差分析主要关注其高程精度,通过统计分析、误差直方图和误差空间结构分析3个方面进行。首先,通过逐像素比较ASTER、AW3D30、MERIT、SRTM和TanDEM-X 5种产品和机载参考DEM以及机载参考DSM产品,计算出5种产品分别表征DEM和DSM时的高程误差,计算公式为:
(1) 式中,
是第j个像素的高程误差; 是第j个像素的待评估DEM产品的高程值; 为第j个像素的参考DEM或DSM产品的高程值。之后选取平均误差(mean error,ME)、RMSE,以及90%线概率误差(linear error at 90% probability,LE90)3个统计指标[1]来定量评估各产品精度。ME可以反映5种 DEM产品高程误差的整体趋势[20],ME与0的偏差可以指示开放DEM产品对地形的表达存在高估或低估。RMSE用来衡量观测值同真值之间的偏差,它可以排除误差值正负号的干扰,反映观测值精度的可靠性。此外,考虑到高程误差有时受到异常值的影响,以及可能存在系统偏差而呈现非正态分布,因此本文根据前人研究建议,选择了具有稳健性的统计指标LE90,来定量比较不同DEM产品之间的精度[1,29]。ME、RMSE以及LE90的计算公式为: (2) (3) (4) 式中,M、R、L分别表示ME、RMSE、LE90的值;
表示将所有像素的高程误差绝对值(即 )从小到大排序后,取90%分位数的值。误差直方图是对高程差值的统计结果的直观展示,在计算出逐像素高程误差的基础上,分析误差的频率直方图并辅助拟合曲线,评估5种开放DEM产品的误差是否符合正态分布。此外,由于使用机载激光扫描最大的优势是能提供面状的数据,所以本文制作了误差的空间分布图,结合研究区的坡度图以及归一化数字表面模型,定量分析高程误差随着坡度以及地表物体高度的变化规律。 3 结果分析与讨论
3.1 高程误差直方图及统计分析结果
当以机载LiDAR生成的DEM为参考真值时,5种开放DEM产品的误差频率分布直方图如图4所示。结果显示5种产品中,ASTER、MERIT和SRTM的误差近似呈现正态分布,而AW3D30的误差呈现明显的右偏态分布,TanDEM-X产品则呈现微弱双峰的误差分布。相对参考的地形表面,5种DEM产品都存在明显的正方向偏移,并且以MERIT产品的偏移程度最小(ME为2.82 m),TanDEM-X产品的正向偏移程度最大(ME为10.39 m)。这表明5种开放DEM产品对实验区域的地面高程值都存在明显高估,并且MERIT的高估程度最小,TanDEM-X的高估程度最大。MERIT的曲线分布最为狭窄且陡峭,这是由于MERIT的DEM产品去除了树高,并未出现整体的显著高估现象(图4(c))。5种开放DEM产品中,AW3D30产品相对ASTER和SRTM而言,RMSE和LE90误差都更大,这一发现与文献[28]的结论不一致。文献[28]对比研究AW3D30、SRTM和ASTER这3种DEM产品,得出AW3D30产品的高程精度最高。这种差异很可能是由于文献[28]的精度评价中,使用的参考DEM数据是TanDEM-X的12 m分辨率DEM产品,而本文使用的是机载LiDAR生成的DEM产品。本文与文献[30]发现TanDEM-X 产品的精度低于AW3D30产品的结论相同,SRTM产品相对ASTER产品表征DEM的精度更高,这与文献[31-32]的研究发现一致。值得注意的是,5种开放DEM产品中只有MERIT和SRTM表征DEM时,LE90是小于15 m的。
当以机载LiDAR生成的DSM为参考真值时,5种开放DEM产品的误差频率直方图如图5所示。5种产品与参考机载DSM之间均存在着明显的负方向偏移,表明它们对地表物体的高程都存在低估。并且其中MERIT产品对地表高程的低估程度最高(ME为-11.26 m),这是由于MERIT校正了树高导致,所以其ME绝对值最大、低估现象最为明显,而且由于MERIT产品只校正了树高,没有校正建筑物的高度,所以呈现了双峰分布的现象。此外,5种产品中,TanDEM-X产品的ME绝对值最小,对地表物体表面高程的低估现象最小,这可能是由于其数据采集时间距离实验区机载LiDAR参考数据采集时间2019年最近,所以地表的相对变化概率更小。此外,TanDEM-X产品与AW3D30产品的误差也显著低于其他3种产品,这表明TanDEM-X和AW3D30比其他产品表征DSM的效果更好。SRTM比ASTER和AW3D30的效果更差,可能是由于SRTM的采集时间在2000年,相距数据采集的2019年最为久远,所以误差更大。
所有误差统计指标如表2所示,对比表1可以发现5种产品在本文研究区的实际精度都低于各产品的目标设计精度。参考机载DEM时,MERIT产品的ME最小,SRTM产品的RMSE
表 2 各开放DEM产品的高程误差统计表/mTable 2. Error Statistics of Each Open DEM Product/mDEM产品 参考机载DEM 参考机载DSM ME RMSE LE90 ME RMSE LE90 ASTER 7.91 12.53 21.17 -6.17 11.31 18.49 AW3D30 8.74 12.43 21.30 -5.34 10.10 16.81 MERIT 2.82 7.32 12.19 -11.26 15.61 25.47 SRTM 4.74 7.08 11.75 -9.34 12.91 20.82 TanDEM-X 10.39 13.80 22.83 -3.69 10.18 16.99 以及LE90最小,TanDEM-X产品的3项误差值都最大,这表明5种产品作为DEM使用时,MERIT和SRTM精度更高,而TanDEM-X精度最低。参考机载DSM时,TanDEM-X产品的ME绝对值最小,AW3D30产品的RMSE和LE90最小,MERIT产品的3项误差值最大,这表明5种产品作为DSM使用时,TanDEM-X和AW3D30精度更高,而MERIT精度最低。类似地,文献[22]发现以机载LiDAR生成的DSM为参考时,TanDEM-X产品的精度始终优于SRTM和ASTER产品[22]。
另外,根据表2还可以发现,ASTER和AW3D30两种产品在表示DSM时,比表示DEM时的RMSE更小(ASTER:11.31 m<12.53 m;AW3D30:10.1 m<12.43 m),相应的LE90也更小(ASTER:18.49 m<21.17 m;AW3D30:16.81 m<21.30 m)。这是由于ASTER和AW3D30的测高原理都是三维立体像对测量,所以更适合表达DSM。
3.2 高程误差的空间分布
参考机载DEM时,5种开放DEM产品的高程误差空间分布如图6所示。对比分析可知, TanDEM-X产品中大部分像素呈现红色,表示高估了地形表面的高程。而ASTER产品中,误差大于15 m以及小于或等于-15 m的像素都很多,表明既有显著低估区域也有显著高估区域。对比图2的nDSM和图1研究区遥感影像,可知ASTER、AW3D30以及TanDEM-X的DEM产品都是显著高估了森林覆盖区域的地面高程(>15 m),而MERIT产品由于校正了树高,所以在这一范围内误差较低(≤15 m)。文献[33]也发现在森林覆盖区域,MERIT表征地形表面高程的精度优于SRTM和TanDEM-X产品[33]。SRTM由于采集时间在20年前,很可能由于植被的生长,且雷达具有一定的穿透森林冠层能力,而光学影像不具备,所以其高估现象不如ASTER、AW3D30、TanDEM-X产品显著。此外,在研究区西北角的居民地覆盖区域,MERIT、SRTM和TanDEM-X产品对地面高程的高估现象都不严重(<5 m),这表明误差的空间分布结构受到了土地覆盖类型的影响。
参考机载DSM时,5种开放DEM产品的高程误差空间分布如图7所示。与图6的情况相反,5种产品在表征地表物体高程时,大部分区域都呈现出低估现象(蓝色像素比例高)。此时,MERIT的低估现象最为严重,在森林区域呈现显著低估,在研究区西北角和东南角的居民地区
域,则呈现较为微弱的误差(浅粉色和浅蓝色)。TanDEM-X产品在森林覆盖区域比居民地区域的低估更为严重。SRTM相对ASTER、AW3D30产品的低估现象更严重,应该是由于其数据采集时间最为久远,所以由于植被生长造成了更严重的低估误差。
3.3 高程误差随地形坡度及地表物体高度的变化
本文进一步定量分析了高程误差随地表物体高度和坡度的变化规律。由图8可知,在表征地形表面高程时,随着地表物体高度的增加,5种DEM产品对地面高程的高估程度都在增加(如图8(a)显示ME递增),且RMSE也在增大(如图8(b)显示RMSE递增),其中MERIT产品的ME在每个高度范围内都是最小的,0~20 m范围内,MERIT产品的RSME最小,大于20 m后SRTM产品的RMSE最小。相对而言,坡度对DEM精度的影响不如地表物体高度的影响大。随着坡度增加,5种DEM产品的高估程度都呈现先增加后减少的趋势(如图8(c)显示ME先递增后递减),且RMSE也是先增大后减小(见图8(d)),在坡度约20°处高估程度最高。对比以往研究,文献[30]在分析TanDEM-X产品在中国区域的误差规律时发现,随着坡度从0°到50°增加,TanDEM-X产品的RMSE呈现出明显增大的趋势;文献[34]发
现在山西省中部,SRTM和ASTER的DEM产品的RMSE随着坡度增大一直呈增大趋势;文献[24]在青藏高原东北缘,发现当坡度从0°增长到5°时,SRTM和ASTER的DEM产品的RMSE微弱降低,而随着坡度从5°增长到30°,RMSE逐渐增大。就整体趋势而言,本文和前人研究都发现坡度会制约DEM产品的精度,但是具体的变化规律有所差异,主要原因是本文实验区和前人研究实验区的地形复杂度和地表覆盖类型有所不同。由图3可知,本文研究区内坡度超过20°的区域占比很小,未来实验中可以探索坡度起伏更大的区域下不同DEM产品的精度。
在表征DSM时,据图9可知随着地表物体高度的增加,5种产品对地表物体高程的低估程度都在增强(如图9(a)显示ME递减),而且误差也在增加(如图9(b)显示RMSE递增)。TanDEM-X产品的ME绝对值在0~25 m高度范围内是最小的,表示其低估程度最弱。随着地面坡度的增加,5种产品对地表物体高程的低估程度都呈现先增大后减小的趋势(如图9(c)显示ME先递减后递增),误差也呈现先增大后减小的趋势(如图9(d)显示RMSE先递增后递减)。5种产品中,由于MERIT校正了树高,在每个坡度区间其误差都是最大的。
4 结 语
本文以机载LiDAR数据生成的高精度DEM和DSM作为参考真值,对比分析了5种常用的全球开放DEM产品,包括ASTER(v003)、AW3D30 (v3.1)、MERIT(v1.0.3)、SRTM(GL1 v003)以及TanDEM-X(90 m,v01)分别表征DEM和DSM时的高程误差。以美国史密森尼学会生物保护研究所实验站为研究区,统计5种开放DEM产品与参考高程真值之间的ME、RMSE、LE90,并分析误差频率直方图、误差空间分布,以及误差随地表物体高度和坡度的变化,进行了不同产品的精度评估比较研究。实验结果表明:
1)5种开放DEM产品的高程精度存在明显差异,在本文研究区的实际精度都低于各自产品的目标设计精度。
2)当作为DEM使用时,5种产品都存在显著高估。MERIT(ME为2.82 m,LE90为12.19 m)和SRTM(ME为4.74 m,LE90为11.75 m)精度更高,TanDEM-X精度最低(ME为10.39 m,LE90为22.83 m)。MERIT产品对地面高程的低估程度最低。
3)当作为DSM使用时,5种产品都存在显著低估。TanDEM-X(ME为-3.69 m,LE90为16.99 m)和AW3D30(ME为-5.34 m,LE90为16.81 m)的精度更高,MERIT的精度最低(ME为-11.26 m,LE90为25.47 m)。TanDEM-X的产品对地表物体高程的低估程度最弱。
4)误差的空间分布结构受到土地覆盖类型的影响。作为DEM使用时,ASTER、AW3D30以及TanDEM-X都显著高估了森林覆盖区域的地面高程(>15 m),作为DSM使用时,TanDEM-X产品在森林覆盖区域比在居民地区域的高程低估现象更为严重。
5)地表物体的高度对5种DEM产品的高程误差有强烈影响。不论作为DEM或者DSM使用,5种DEM产品的RMSE和LE90误差都随着地表物体高度的增加,而呈显著增大的趋势。相对而言,坡度因素对5种产品的误差影响不如地表物体的高度显著。
本文充分利用了机载LiDAR数据,分析了开放DEM产品分别表征地形表面高程和地表物体高程时的精度高低与误差规律,根据实验结果推荐表征地形高程时使用MERIT或者SRTM产品,表征地表物体高程时使用TanDEM-X产品。本文的结论可以为今后在类似土地覆盖类型与地形起伏区域的研究中DEM产品的选择提供参考。
-
[1] Georgescu M,Moustaoui M,Mahalov A,etal.Summer timeClimateImpactsofProjected Mega politanExpansioninArizona[J].犖犪狋狌狉犲犆犾犻犿犪狋犲犆犺犪狀犵犲,2013,3(1):37 41[2] NicholJ.RemoteSensingofUrbanHeatIslandsbyDayandNight[J].犘犺狅狋狅犵狉犪犿犿犲狋狉犻犮犈狀犵犻狀犲犲狉犻狀犵犪狀犱犚犲犿狅狋犲犛犲狀狊犻狀犵,2005,71(5):613 621[3] Jiménez MuozJC,SobrinoJA.A GeneralizedSingle channelMethodforRetrievingLandSurfaceTemperaturefromRemoteSensingData[J].犑狅狌狉狀犪犾狅犳 犌犲狅狆犺狔狊犻犮犪犾犚犲狊犲犪狉犮犺,2003,108(D22):4688,doi:10.1029/2003JD003480[4] QinZ,KarnieliA,BerlinerP.A Mono windowAl gorithmforRetrieving LandSurfaceTemperaturefromLandsatTM DataandItsApplicationtotheIsrael EgyptBorderRegion[J].犐狀狋犲狉狀犪狋犻狅狀犪犾犑狅狌狉狀犪犾狅犳犚犲犿狅狋犲犛犲狀狊犻狀犵,2001,22(18):3719 3746[5] Jiménez MuozJC,SobrinoCJ,SoriaJA,etal.RevisionoftheSingle channelAlgorithmforLandSurfaceTemperatureRetrievalfrom LandsatTher mal infraredData[J].犐犈犈犈犜狉犪狀狊犪犮狋犻狅狀狊狅狀犌犲狅狊犮犻犲狀犮犲犪狀犱犚犲犿狅狋犲犛犲狀狊犻狀犵,2009,47(1):339 349[6] SobrinoJA,Jiménez MuozJC,SoriaG,etal.Land Surface Emissivity Retrievalfrom DifferentVNIRandTIRSensors[J].犐犈犈犈犜狉犪狀狊犪犮狋犻狅狀狊狅狀犌犲狅狊犮犻犲狀犮犲犪狀犱犚犲犿狅狋犲犛犲狀狊犻狀犵,2008,46(2):316 327[7] ICESS.MODIS UCSB EmissivityLibrary[OL].http://www.icess.ucsb.edu/modis/EMIS/html/em.html,2013[8] NicholJ.An Emissivity Modulation MethodforSpatialEnhancementofThermalSatelliteImagesinUrbanHeatIslandAnalysis[J].犘犺狅狋狅犵狉犪犿犿犲狋狉犻犮犈狀犵犻狀犲犲狉犻狀犵犪狀犱 犚犲犿狅狋犲犛犲狀狊犻狀犵,2009,75(5):547 556[9] BarsiJA,SchottJR,PalluconiFD,etal.Valida tionofa Web based AtmosphericCorrection ToolforSingleThermalBandInstruments[C].SPIE,SanDiego,CA,2005[10]USGS.Landsat8(L8)OperationalLandImager(OLI)andThermalInfraredSensor(TIRS)Cali brationNotices[OL].http://landsat.usgs.gov/calibration_notices.php,2013[11]USGS.WhataretheBand DesignationsfortheLandsatSatellites? [OL].http://landsat.usgs.gov/band_designations_landsat_satellites.php,2013[12]ChanderG,MarkhamBL,HelderDL.SummaryofCurrentRadiometricCalibrationCoefficientsforLandsatMSS,TM,ETM+,andEO 1ALISen sors[J].犚犲犿狅狋犲犛犲狀狊犻狀犵狅犳犈狀狏犻狉狅狀犿犲狀狋,2009,113(5):893 903[13]CollC,GalveJM,SánchezJM,etal.ValidationofLandsat7/ETM+ ThermalBandCalibrationandAtmosphericCorrection with Ground Based Meas urements[J].犐犈犈犈 犜狉犪狀狊犪犮狋犻狅狀狊狅狀 犌犲狅狊犮犻犲狀犮犲犪狀犱犚犲犿狅狋犲犛犲狀狊犻狀犵,2010,48(1):547 555554 第40卷第4期 徐涵秋等:单通道算法地表温度反演的若干问题讨论———以Landsat系列数据为例犛狅犿犲犐狊狊狌犲狊犻狀犔犪狀犱犛狌狉犳犪犮犲犜犲犿狆犲狉犪狋狌狉犲犚犲狋狉犻犲狏犪犾狅犳犔犪狀犱狊犪狋犜犺犲狉犿犪犾犇犪狋犪狑犻狋犺狋犺犲犛犻狀犵犾犲犮犺犪狀狀犲犾犃犾犵狅狉犻狋犺犿犡犝犎犪狀狇犻狌1 犔犐犖犣犺狅狀犵犾犻1,2 犘犃犖 犠犲犻犺狌犪31 CollegeofEnvironmentandResources,FuzhouUniversity,Fuzhou350108,China2 FujianProvincialKeyLaboratoryofRemoteSensingofSoilErosion,Fuzhou350108,China3 InstituteofMeteorologicalScienceofFujianProvince,Fuzhou350001,China犃犫狊狋狉犪犮狋:Thesingle channelalgorithm(SC)developedbyJiménez MuozandSobrinoin2003hasbeenfrequentlyusedinretrievinglandsurfacetemperature(LST).However,manystudiesemployedtheLandsatTM’satmosphericparametersfortheLSTretrievalofETM+thermaldatabecausetheparametersforETM+ werenotprovidedin2003.AlthoughJiménez MuozandSobrinorevisedtheSCalgorithmandprovidedtheatmosphericparametersspecializedforETM+in2009,suchmisusingstillnotchanges.Therefore,thispapercomparedtheLSTsretrievedfrom LandsatTM/ETM+/TIRSthermaldatawiththeSC’s2003and2009versions,discussedthecorrectusageofthePlanckfunctionandwavelengthvalue,andexploredthepotentialoftheSCforretrievingLSTfromLandsat8TIRSdata.Resultsshowthatthe2009versioncanachievehigheraccuracyaveragelyforETM+thermaldatathanthe2003version,whencalculatedusingwavelengthvalueprovidedbyUSGSandthebrightnesstemperaturecomputedwithPlanckfunction.ThisstudyalsorevealsthatatthisstagetheLSTofLandsat8couldberetrievedusingtheSC withTIRSband10.Thiscouldbeachievedwhentheatmosphericwatervaporcontentisrelativelylow.犓犲狔狑狅狉犱狊:single channelalgorithm;landsurfacetemperature;thermal infraredremotesensing犉犻狉狊狋犪狌狋犺狅狉:XU Hanqiu,PhD,professor,specializesinremotesensingofenvironmentandresources.Email:hxu@fzu.edu.cn犉狅狌狀犱犪狋犻狅狀狊狌狆狆狅狉狋:TheNationalScienceandTechnologySupportingProgram,No.2013BAC08B01 05. -
期刊类型引用(43)
1. 宋婷婷. 浙江山地区村镇景观格局指数与地温相关性分析——以金华市为例. 现代园艺. 2024(03): 41-44 . 百度学术
2. 邢艳春,甯珂,李雪兰. 松辽流域生态环境时空动态评价及驱动因素分析. 中国水利水电科学研究院学报(中英文). 2024(01): 84-96+107 . 百度学术
3. 柳芬芬,徐涵秋. 耦合卫星热红外与气象数据的广州城市热环境演变研究. 遥感技术与应用. 2024(06): 1500-1511 . 百度学术
4. 丁雪,冯婧文,黄园园,施骏骋,王金亮. 2000—2020年滇中城市群生态环境质量动态监测及空间格局演变. 水土保持通报. 2023(03): 96-104+128 . 百度学术
5. 洪婷婷,黄晓辉,邓西鹏,杨义炜,唐翔. 基于MSPA的城市绿色基础设施与热环境关系研究——以福州市中心城区为例. 中国园林. 2023(10): 97-103 . 百度学术
6. 李伟,陈朝霞,周美玲. 耦合夜光遥感的粤港澳不透水面提取及其热岛效应研究. 城市勘测. 2023(06): 31-37 . 百度学术
7. 石希,夏军强,孙健. 基于热红外遥感影像的河流水温反演方法比较——以长江上游流域为例. 湖泊科学. 2022(01): 307-319 . 百度学术
8. 晏红波,杨志高,卢献健,韦晚秋,黎振宝. 基于改进RSEDI的典型喀斯特地区生态环境质量时空变化. 科学技术与工程. 2022(11): 4646-4653 . 百度学术
9. 汪东川,陈星,孙志超,辛燕,王海庆,柴华,王鸿艺. 格尔木长时间序列遥感生态指数变化监测. 生态学报. 2022(14): 5922-5933 . 百度学术
10. 朱菲,苏涛,王建,夏俊,廖晋一. 基于地表参数变化的南京市热岛效应时空演变特征分析. 遥感信息. 2022(04): 106-113 . 百度学术
11. 万红,郭鹏,骆磊,赵颜创,赵燕,王心源. 基于生态多稳态的山地植被垂直自然带定量识别研究——以新疆博格达山为例. 遥感学报. 2022(11): 2234-2247 . 百度学术
12. 胡欣雨,魏凌伟,莫惠萍,陆东芳. 近15年郑州市主城区热环境与下垫面指数关系研究. 河南农业大学学报. 2021(01): 152-163 . 百度学术
13. 郄宇凡,王宁练,吴玉伟,陈安安. 山地冰川表面温度反演算法对比——以祁连山七一冰川为例. 山地学报. 2021(01): 129-142 . 百度学术
14. 王丹,杨超,郑永宏,肖雪莹,赵林,陈正洪. 海绵城市改造对热岛效应减缓效果的估算. 长江流域资源与环境. 2021(04): 968-975 . 百度学术
15. 朱江山,杨丽萍,白宇兴,苏志强,孔金玲. Landsat-8影像地表温度反演算法对比分析. 科学技术与工程. 2021(12): 4785-4792 . 百度学术
16. 段四波,茹晨,李召良,王猛猛,徐涵秋,历华,吴鹏海,占文凤,周纪,赵伟,任华忠,吴骅,唐伯惠,张霞,尚国琲,覃志豪. Landsat卫星热红外数据地表温度遥感反演研究进展. 遥感学报. 2021(08): 1591-1617 . 百度学术
17. 吴蓉蓉. 矿区地表温度反演算法比较. 河南科技. 2020(04): 149-151 . 百度学术
18. 陈梅花. 10年尺度金华市主城区夏季热环境变化及地表覆盖的影响. 浙江师范大学学报(自然科学版). 2020(03): 320-327 . 百度学术
19. 崔文杰,李忠,李家国,朱利,张永红. GF-5卫星模拟单通道SST反演. 遥感信息. 2020(03): 57-62 . 百度学术
20. 周鹏,谢元礼,高志远,蒋广鑫. 西安市不透水面的变化及其驱动力. 水土保持通报. 2020(03): 274-281 . 百度学术
21. 张寒,潘安,罗芳. 基于遥感生态指数的南充市生态环境质量评价. 西华师范大学学报(自然科学版). 2020(04): 395-403 . 百度学术
22. 张娟,文广超,王恩营,谢洪波,刘正疆. 基于遥感生态指数的焦作市生态环境动态监测与评价. 水土保持通报. 2020(06): 107-114 . 百度学术
23. 侯宇初,张冬有. 基于Landsat8遥感影像的地表温度反演方法对比研究. 中国农学通报. 2019(10): 142-147 . 百度学术
24. 陈冰倩,张友水,程璟媛,赵雪. 福州市地表温度热点及时空变化分析. 地球信息科学学报. 2019(05): 710-719 . 百度学术
25. 杨松勇,刘士彪,刘昭华. 中小城市热岛效应及地表参数相关性研究. 江西理工大学学报. 2019(03): 22-29 . 百度学术
26. 林中立,徐涵秋. 近20年来新旧“火炉城市”热岛状况对比研究. 遥感技术与应用. 2019(03): 521-530 . 百度学术
27. 李鸿芝. 基于RSEI的巢湖流域生态环境质量时空变化分析. 安徽农业科学. 2019(16): 82-86 . 百度学术
28. 林中立,徐涵秋. 基于遥感的海岛型城市发展生态效应分析——以厦门岛为例. 福州大学学报(自然科学版). 2019(05): 610-616 . 百度学术
29. 赵珍,李贵仁,陈刚. ArcGIS和RS技术在高植被覆盖率地区地热资源勘查中的应用. 勘察科学技术. 2019(02): 44-48 . 百度学术
30. 侯浩然,丁凤,黎勤生. 近20年来福州城市热环境变化遥感分析. 地球信息科学学报. 2018(03): 385-395 . 百度学术
31. 林中立. 基于新型遥感生态指数的福州市生态状况变化研究. 海峡科学. 2018(06): 33-36 . 百度学术
32. 王植,董斌,陈炟君. 基于Landsat的沈阳城市热岛效应与地表参数变化分析. 测绘与空间地理信息. 2018(10): 4-7+11 . 百度学术
33. 李军,龚围,辛晓洲,高阳华. 重庆地表温度的遥感反演及其空间分异特征. 遥感技术与应用. 2018(05): 820-829 . 百度学术
34. 金点点,宫兆宁. 基于Landsat系列数据地表温度反演算法对比分析——以齐齐哈尔市辖区为例. 遥感技术与应用. 2018(05): 830-841 . 百度学术
35. 王植,陈炟君,董斌. 基于TVDI的2013—2017年沈阳市旱情等级评估. 测绘与空间地理信息. 2018(11): 12-17 . 百度学术
36. 王美雅,徐涵秋,李霞,林中立,张博博,付伟,唐菲. 不透水面时空变化及其对城市热环境影响的定量分析——以福州市建成区为例. 应用基础与工程科学学报. 2018(06): 1316-1326 . 百度学术
37. 林中立,徐涵秋. 基于LCZ的城市热岛强度研究. 地球信息科学学报. 2017(05): 713-722 . 百度学术
38. 张春森,徐肖雷,陈越峰. 基于ETM+数据的煤田火区温度异常信息提取. 国土资源遥感. 2017(02): 201-206 . 百度学术
39. 林中立,徐涵秋. 基于DTR的城市热岛强度指数构建与应用——以福州市为例. 应用基础与工程科学学报. 2017(02): 266-275 . 百度学术
40. 王美雅,徐涵秋,付伟,林中立,李霞,张博博,唐菲. 城市地表水体时空演变及其对热环境的影响. 地理科学. 2016(07): 1099-1105 . 百度学术
41. 陈洁羽,黎勤生,丁凤,吴文丰. 基于IDL的两种地表温度反演算法实现. 福建电脑. 2016(05): 3-5 . 百度学术
42. 王冕,朱明邦,赵之重. 基于Landsat 5热红外波段反演亮度温度的比较研究. 青海大学学报(自然科学版). 2016(02): 40-43+52 . 百度学术
43. 徐涵秋. Landsat 8热红外数据定标参数的变化及其对地表温度反演的影响. 遥感学报. 2016(02): 229-235 . 百度学术
其他类型引用(47)
计量
- 文章访问数: 2856
- HTML全文浏览量: 122
- PDF下载量: 1451
- 被引次数: 90