目的
地表温度反演是一个典型的病态反演问题，深度动态学习神经网络的出现提供了一条新的地表温度反演途径。文章以MODIS中红外和热红外波段作为参照模拟研究对象，利用深度动态学习神经网络和辐射传输模型（MODTRAN）进行地表温度反演研究，选择最适合于MODIS地表温度反演的波段组合，从而为国产卫星风云系列和高分数据红外波段反演地表温度提供参考算法。
方法
根据中红外波段受太阳的影响以及水汽波段的特征，将反演组合波段分成3组。第1组适合白天和晚上同时反演地表温度的组合（MODIS波段29、31、32和33）；第2组适合白天的热红外波段和水汽波段组合（MODIS波段29、31、32、33和水汽波段）；第3组是只适合晚上的中外波段（MODIS 20、22、23）与热红外波段（MODIS 29、31、32和33）的组合。
结果
利用辐射传输模型（MODTRAN）和深度动态神经网络（NN）反演分析表明，深度动态学习神经网络能够被用来精确地从单景MODIS数据中反演地表温度，克服了传统MODIS白天/黑夜产品算法的缺陷。3种类型的组合地表温度的平均反演误差都在1 K以下，最高精度为热红外波段与水汽波段的组合，平均最高精度为0.251 K，标准差是0.255 K，相关系数是1。
结论
利用深度动态学习神经网络和辐射传输模型彻底解决了地表温度和发射率病态反演难题，为风云系列卫星和高分数据地表温度反演算法提供参考算法模式，深度动态学习神经网络与辐射传输模型相结合反演地表温度和发射率在地表温度反演史上具有里程牌意义。
Purpose
Land surface temperature retrieval is a typical ill-posed problem. The deep dynamic learning neural network provides us with a new method to retrieve land surface temperature. In this paper，the mid-infrared and thermal infrared bands of MODIS are used as reference simulation objects. The deep dynamic learning neural network and radiation transmission model （MODTRAN） are used to study the surface temperature retrieval，and select which band combination is most suitable for MODIS surface temperature retrieval.
Methods
According to the influence of the sun for the mid-infrared band and the characteristics of the water vapor band，the retrieval combined band is divided into three types. The first group is suitable for simultaneous retrieval of surface temperature during the day and night （MODIS bands 29，31，32 and 33）；the second group is thermal infrared and water-wave band combinations （MODIS bands 29，31，32，33 and the water vapor band） which is suitable for daytime；the third group is a combination of the mid-infrared band （MODIS 20，22，23） and the thermal infrared band （MODIS 29，31，32 and 33） for the evening.
Results
The retrieval analysis using the radiation transfer model （MODTRAN） and deep dynamic neural network （NN） shows that the deep dynamic learning neural network can be used to accurately retrieve the surface temperature from the MODIS data. The average retrieval error of the three types of combined surface temperatures is below 1K. The highest precision is the combination of the thermal infrared band and the water vapor band. The average highest precision is 0.251 K，the standard deviation is 0.255K，and the correlation coefficient is 1.
Conclusion
It has the significance of mileage in the history of surface temperature retrieval when the deep dynamic learning neural network combined with radiation transmission model is used to retrieve surface temperature and emissivity.
0 引言
温度是表征地球各圈层（岩石圈、水圈、大气圈和生物圈）之间能量传输的重要物理量，地面温度（包括地表和近地表）是研究地表和大气之间物质和能量交换、全球海洋环流、气候变化异常等方面不可或缺的重要参数，涉及众多基础学科和重大应用领域。运用卫星遥感技术快速准确地获取大面积、长时间序列的地表温度，是顺应当前科学技术发展趋势，是应对多种自然灾害等诸多问题迫切需要开展的研究课题，其研究成果既具有重大的科学意义，同时也具有重要的社会经济价值。最近几十年，大尺度上环境研究和管理活动对温度信息的需求已经使得利用遥感手段反演地表温度和发射率的技术有了很大进步。MODIS是对地观测卫星（EOS）上的一个传感器，它拥有36波段。能够覆盖全球、辐射分辨率非常高、动态范围以及精确校正等特点使得MODIS传感器在海面温度、陆地表面温度和大气参数反演方面具有非常大的潜力。
从多个热红外波段的测量中同时反演地表温度和发射率非常困难，因为N个波段至少拥有N+1个未知数（N个波段的发射率和地表温度），这是一个非常典型的病态反演问题。如果不利用任何先验知识，几乎不可能同时从多个热红外数据中反演地表温度和发射率。许多劈窗算法已经被开发来从NOAA/AVHRR和MODIS数据中反演海面温度和地表温度，这些反演方法主要是利用相邻热红外波段对水汽的敏感性不一样来校正大气水汽的影
响［1,2,3,4,5,6,7,8,9,10,11,12,13,14］ 。这些算法形式基本上相同，主要差别在于关键参数的获取及计算。虽然大多数的算法精度都很高，但它们仍然需要做一些假定并利用发射率和大气状态（特别是大气水汽含量）作为已知的先验知识。许多研究人员［15,16,17,18,19,20,21］ 在地表温度和发射率反演方面做了许多工作，特别是Li和Becker利用TISI（Temperature Independent Spectral Index）技术和白天/晚上的影像数据反演地表温度和发射率［22］ ，这个反演方法需要大气剖面信息。Wan和Li提出了一个多波段算法从EOS/MODIS传感器中同时反演地表温度和发射率［23］ ，这个方法受地面的光学性质和大气状态的影响。这两个方法都假定了白天和晚上同一地点的发射率不变，由于云等条件的影响，使得产品的分辨率降低为5 km。毛克彪等提出利用辐射传输方程与人工智能相结合的方法反演地表温度，该方法的优点是只需要使用同一景热红外遥感数据就可以反演地表温度和发射率，克服了白天/晚上算法的缺陷［24,25,26,27,28,29］ 。地球上各种成分不是独立的，而是彼此相互联系的。毛克彪等首次系统地阐述了利用地球物理参数之间的关系构造额外方程消除反演病态问题，并用人工智能神经网络基于辐射传输模型等方法进行优化计
算［29］ 。但毛克彪等只对MODIS白天的热红外波段和水汽波段反演做了评价分析［29］ ，没有对适合于白天和晚上的热红外波段做模拟分析，也没有对近红外波段与热红外波段联合反演晚上的温度做分析。该研究拟利用MODTRAN4和深度动态学习神经网络分别对适合于白天、白天/晚上和晚上的热红外和近红外反演地表温度做单独模拟评价分析，从而为风云系列卫星和高分数据地表温度反演提供参考模式。1 地表温度反演基本原理
地表温度反演是基于辐射传输方程，利用地表热辐射从地表途经大气达到传感器的过程建立热平衡方程。通常讲，地表不是黑体，在计算地表热辐射时需要考虑地表发射率，同时热辐射传输过程中大气对传感器接收到的能量也有贡献。考虑到这些因素，通用的辐射传输方程可以描述成方程（1）。
Bi （Ti ）=τi （θ）［εi （θ）Bi （Ts ）+（1-εi （θ））I↓ i］+I↑ i （1）式（1）中，
Ts 是地表温度，Ti 是通道i的星上亮度温度，τi （θ）是通道i在θ角时的透过率，εi （θ）是通道i在θ角时的地表发射率。Bi （Ts ）是地表辐射强度，I↑ i和I↓ i分别是大气向上和向下的辐射强度。I↑ i和I↓ i可分别用如下公示表示：I↑ i=（1-τi （θ））Bi （Ta ） （2a）I↓ i=（1-τi （θ′））Bi （T↓ a） （2b）Ta 是向上的大气平均作用温度，T↓ a是向下的大气平均作用温度。将I↑ i和I↓ i代入式（1）可得：Bi （Ti ）=εi （θ）τi （θ）Bi （Ts ）+［1-τi （θ′）］［1-εi （θ）］τi （θ）Bi （T↓ a）+［1-τi （θ）］Bi （Ta ） （3）Qin et al
.［11］ 、Mao et al.［13］ 对等式（3）做了一些分析和合理的简化，得出用Ta 代替T↓ a对反演方程不会有太大的影响，所以式（3）可以描述成式（4）。Bi （Ti ）=εi （θ）τi （θ）Bi （Ts ）+［1-τi （θ′）］［1-εi （θ）］τi （θ）Bi （Ta ）+［1-τi （θ）］Bi （Ta ） （4）每个方程中，有3个未知数（1个波段的发射率、地表温度和大气平均作用温度），这是一个典型的病态问题。为了反演地表温度和发射率，需构造其它的方程和限制条件。地球物理参数彼此之间不是独立的，而是相互联系的。在以往的算法中，没有充分利用地球物理参数之间的关系。毛克彪
等［25,26,27,28,29］ 利用波普曲线对不同波段的发射率之间的关系，大气水汽含量与透过率之间的关系，近地表空气温度和大气平均作用温度与地表温度以及星上亮温做了分析，表明这些物理参数之间存在某种内在的联系。2 深度动态学习神经网络和辐射传输模型解病态反演问题
地球物理参数之间的关系大多是非线性的，很难准确地用数学公式描述出来。从遥感数据中反演地球物理参数，由于其中非线性的关系和相互作用的因素很难描述清楚，传统的方法在简化的过程中使得反演精度下降。由于神经网络具备从复杂的和不精确的数据中提取信息，所以神经网络能够被用来提取模式和预
测［30,31］ 。神经网络具备函数近似、分类和优化计算的能力，动态学习神经网络是解病态反演问题的最好方法之一。神经网络与传统的方法不一样，不需要准确知道输入参数和输出参数之间的具体关系。神经网络通过训练数据直接决定了输入数据和输出数据之间的关系。毛克彪等［29］ 分析表明地表温度和发射率能够通过3个热红外波段以及地球物理参数之间的关系来反演。该研究将针对MODIS的中红外和热红外波段进行模拟分析，利用辐射传输模型（MODTRAN4）和深度动态学习神经网络进行反演分析，找到分别适合于MODIS波段20/22/23（3~5 μm）和波段29/31/32/33（8~13.5 μm）数据的白天/晚上、白天、晚上地表温度反演方法。MODIS波段20、22和23在中红外波段的窗口区，波段29~32在热红外波段的窗口区。虽然33波段在热红外波段窗口区的边缘上，但波段30受臭氧影响很大，因此波段30不适合反演地表温度。辐射传输模型（MODTRAN4）和深度动态学习神经网络执行反演分析可以简单地分为以下4步。
（1）利用MODTRAN4模拟训练和测试数据库，同时可以补充高精度的野外测量数据或MODIS温度产品；
（2）建立训练和测试数据库；
（3）计算波段20，22，23，29，30，31，32，33，利用MODIS波段2，5，17，18，19估算大气水汽含量；
（4）利用训练和网络反演地表温度。
利用MODTRAN4M模拟数据最大的优点在于能够很好地保证同一地物在MODIS波段20/22/23（3~5 μm）和波段29/31/32/33（8~13.5 μm）之间的发射率、水汽等不同物理参数之间的对应关系。利用JPL（URL：http：//speclib.jpl.nasa.gov）测量的地物波谱曲线在MODIS波段20/22/23/29/31/32/33的发射率作为MODTRAN4的输入参数，地表温度变化范围为280~320 K和近地表空气温度变化范围为270~320 K。大气水汽含量变化范围0.2~4.0 g/c
m2 ，观测视角0~20°之间。MODIS的观测视角为0~65°之间，可将其每隔20°建立一个训练数据库，根据角度范围选择不同的训练数据库，从而提高精度。为了简练，该文只分析0~20°之间的数据库分析过程情况，将模拟数据随机分成两部分：训练数据（33800组）和测试数据（15840组），然后用深度动态学习神经网络解病态反演问题。首先，用训练数据训练神经网络，然后用测试数据校验神经网络。根据地表类型和反复试验，部分测试数据信息见表1,2,3,4,5,6,7,8,9,10,11。表1 波段29-31-32组合时反演误差
Table 1 The Summery of Retrieval Error for band 29-31-32 combination
29-31-32 隐含节点 隐含层 300 400 500 600 700 M SD R M SD R M SD R M SD R M SD R 3 0.956 0.879 0.995 0.984 0.865 0.995 0.953 0.894 0.995 0.925 0.873 0.995 0.941 0.901 0.995 4 0.967 1.083 0.994 0.93 0.863 0.995 0.917 0.85 0.995 1.034 2.18 0.982 0.926 0.897 0.995 5 0.907 0.85 0.995 0.956 1.226 0.993 0.941 1.032 0.994 0.977 1.265 0.992 0.937 1.366 0.992 6 0.947 1.039 0.994 1.123 4.04 0.95 1.123 6.05 0.9 4.421 81.11 0.149 3.222 34.628 0.354 注：其中：M为平均误差；R为相关系数；SD为标准偏差，以下同
表2 波段29-31-32-33组合时反演误差
Table 2 The Summery of Retrieval Error for band 29-31-32-33 combination
29-31-32-33 隐含节点 隐含层 300 400 500 600 700 M SD R M SD R M SD R M SD R M SD R 3 0.703 0.593 0.997 0.646 0.568 0.998 0.623 0.57 0.998 0.655 0.591 0.998 0.651 0.592 0.998 4 0.705 0.664 0.997 0.857 3.267 0.967 0.713 0.733 0.997 0.722 0.835 0.996 0.702 0.812 0.996 5 0.673 0.636 0.997 0.678 0.665 0.997 0.707 0.684 0.997 0.802 1.429 0.992 0.764 0.721 0.997 6 0.885 3.249 0.966 0.717 0.715 0.997 0.781 1.032 0.995 0.922 2.764 0.975 1.223 7.457 0.855 表3 波段29-31-32-W组合时反演误差
Table 3 The Summery of Retrieval Error for band 29-31-32-W combination
29-31-32-W 隐含节点 隐含层 300 400 500 600 700 M SD R M SD R M SD R M SD R M SD R 3 0.283 0.276 1 0.33 0.271 0.999 0.273 0.274 1 0.289 0.281 0.999 0.271 0.273 1 4 0.294 0.273 1 0.29 0.279 0.999 0.308 0.723 0.998 0.251 0.255 1 0.259 0.265 1 5 0.285 0.293 0.999 0.291 0.369 0.999 0.34 2.769 0.999 0.917 8.877 0.815 0.546 4.605 0.938 6 0.351 0.313 0.999 0.409 2.548 0.98 0.329 1.087 0.996 0.349 0.637 0.998 0.345 0.68 0.998 注：W为大气水汽含量
表4 波段29-31-32-33-W组合时反演误差
Table 4 The Summery of Retrieval Error for band 29-31-32-33-W combination
29-30-31-32-W 隐含节点 隐含层 300 400 500 600 700 M SD R M SD R M SD R M SD R M SD R 3 0.38 0.318 0.999 0.332 0.26 0.999 0.403 0.316 0.999 0.331 0.27 0.999 0.417 2.763 0.999 4 0.426 0.307 0.999 0.341 0.279 0.999 0.302 0.273 0.999 0.337 0.286 0.999 0.338 0.552 0.999 5 0.488 0.901 0.997 0.319 0.289 0.999 0.372 0.347 0.999 0.3 0.28 0.999 0.332 0.284 0.999 6 0.388 0.331 0.999 0.723 9.754 0.796 0.314 0.407 0.999 0.347 0.48 0.999 0.561 2.724 0.977 表5 波段20-31-32组合时反演误差
Table 5 The Summery of Retrieval Error for band 20-31-32 combination
20-31-32 隐含节点 隐含层 300 400 500 600 700 M SD R M SD R M SD R M SD R M SD R 3 1.015 0.928 0.994 0.973 0.921 0.995 0.977 0.928 0.994 1 0.932 0.994 0.964 0.907 0.995 4 0.98 0.932 0.994 0.958 0.912 0.995 0.957 0.906 0.995 0.952 0.906 0.995 0.963 0.906 0.995 5 0.987 0.928 0.994 0.976 0.933 0.994 0.899 0.909 0.995 0.959 1.184 0.993 0.927 0.907 0.995 6 1.129 0.975 0.993 0.929 0.918 0.995 0.936 0.939 0.995 0.903 0.934 0.995 0.925 1.03 0.994 表6 波段20-29-31-32组合时反演误差
Table 6 The Summery of Retrieval Error for band 20-29-31-32 combination
20-29-31-32 隐含节点 隐含层 300 400 500 600 700 M SD R M SD R M SD R M SD R M SD R 3 1.006 0.928 0.994 0.844 0.78 0.996 0.666 0.673 0.997 0.568 0.608 0.998 0.581 0.61 0.998 4 0.642 0.656 0.997 0.607 0.624 0.998 0.638 0.662 0.997 0.58 0.609 0.998 0.611 0.648 0.998 5 1.022 0.866 0.994 0.666 0.668 0.997 1.149 17.11 0.579 0.602 1.039 0.996 0.619 0.828 0.997 6 0.694 1.512 0.991 0.514 0.568 0.998 0.8 2.572 0.978 0.763 2.067 0.985 0.571 1.37 0.993 表7 波段20-22-31-32组合时反演误差
Table 7 The Summery of Retrieval Error for band 20-22-31-32 combination
20-22-31-32 隐含节点 隐含层 300 400 500 600 700 M SD R M SD R M SD R M SD R M SD R 3 0.838 0.821 0.996 1.204 0.954 0.993 0.782 0.809 0.996 0.776 0.813 0.996 0.719 0.794 0.996 4 0.857 0.813 0.996 0.781 0.812 0.996 0.754 0.808 0.996 0.714 0.775 0.997 0.7 0.737 0.997 5 0.914 0.824 0.995 0.785 0.834 0.996 0.779 0.777 0.996 0.714 0.75 0.997 0.966 4.47 0.94 6 0.944 0.876 0.995 0.709 0.789 0.997 0.76 0.89 0.996 0.718 1.091 0.995 0.765 0.968 0.995 表8 波段20-22-29-31-32组合时反演误差
Table 8 The Summery of Retrieval Error for band 20-22-29-31-32 combination
20-22-29-31-32 隐含节点 隐含层 300 400 500 600 700 M SD R M SD R M SD R M SD R M SD R 3 0.803 0.682 0.997 0.644 0.598 0.998 0.638 0.675 0.997 0.539 0.537 0.998 0.734 0.656 0.997 4 0.716 0.733 0.997 0.584 0.568 0.998 0.749 2.463 0.98 0.563 0.586 0.998 0.694 3.075 0.97 5 0.765 2.469 0.98 0.621 0.621 0.998 0.668 1.286 0.994 0.584 0.709 0.997 0.705 2.254 0.983 6 0.957 0.895 0.995 0.928 7.127 0.867 0.509 0.566 0.998 0.607 1.391 0.993 1.26 10.121 0.774 表9 波段20-22-23-29-31-32组合时反演误差
Table 9 The Summery of Retrieval Error for band 20-22-23-29-31-32 combination
20-22-23-29-31-32 隐含节点 隐含层 300 400 500 600 700 M SD R M SD R M SD R M SD R M SD R 3 0.763 1.38 0.992 0.717 0.646 0.997 0.565 0.796 0.997 0.647 0.724 0.997 0.385 0.43 0.999 4 0.531 0.514 0.998 0.938 5.426 0.915 0.716 3.002 0.972 0.61 1.035 0.996 0.563 1.617 0.991 5 1.037 7.441 0.861 0.544 0.793 0.997 0.479 0.558 0.998 0.434 0.642 0.998 0.73 4.54 0.94 6 0.89 0.816 0.995 0.634 0.868 0.996 0.636 4.363 0.945 0.437 0.607 0.998 0.558 3.967 0.953 表10 波段20-22-29-31-32-33组合时反演误差
Table 10 The Summery of Retrieval Error for band 20-22-29-31-32-33 combination
20-22-29-31-32-33 隐含节点 隐含层 300 400 500 600 700 M SD R M SD R M SD R M SD R M SD R 3 0.838 0.672 0.996 0.734 0.619 0.997 0.706 0.587 0.997 0.622 0.527 0.998 0.63 0.527 0.998 4 0.69 0.606 0.997 0.777 0.667 0.997 0.674 0.575 0.998 0.617 0.526 0.998 0.702 0.8 0.997 5 0.751 0.689 0.997 0.626 0.551 0.998 0.956 3.622 0.959 0.888 1.954 0.986 0.847 1.396 0.992 6 0.783 0.854 0.996 0.755 0.735 0.997 0.585 0.563 0.998 0.562 0.564 0.998 0.643 0.783 0.997 表11 波段20-22-23-29-31-32-33组合时反演误差
Table 11 The Summery of Retrieval Error for band 20-22-23-29-31-32-33 combination
20-22-23-29-31-32-33 隐含节点 隐含层 300 400 500 600 700 M SD R M SD R M SD R M SD R M SD R 3 0.804 0.652 0.997 0.736 0.615 0.997 0.717 0.603 0.997 0.596 0.532 0.998 0.524 0.487 0.998 4 0.814 0.65 0.997 0.605 0.511 0.998 0.869 1.715 0.989 0.557 0.534 0.998 0.478 0.444 0.999 5 0.753 0.64 0.997 0.599 0.537 0.998 0.736 0.701 0.997 0.653 1.562 0.991 0.635 0.748 0.997 6 0.757 0.678 0.997 1.47 24.17 0.47 0.704 0.606 0.997 0.642 0.533 0.998 0.699 0.674 0.997 根据热红外波段和中外波段以及水汽波段的特征，划分成3组。第1组，适合白天和晚上的热红外波段组合，波段29/31/32/33（8~13.5 μm）；第2组适合白天的热红外波段29/31/32/33（8~13.5 μm）和水汽波段；第3组是适合晚上的MODIS波段20/22/23（3~5 μm）与热红外波段29/31/32/33（8~13.5 μm）组合。
对于热红外波段29/31/32/33（8~13.5 μm），测试了不同的组合。表1是波段29-31-32组合，从表1可以看出，当隐含层数是5层且隐含节点为300-300-300-300-300时，平均精度最高为0.907 K，标准差为0.85 K，相关系数为0.995。
表2是波段29-31-32-33组合的误差表，从表中可以看出，当隐含层数是3层并且隐含节点为500-500-500时反演精度最高，最高平均精度为0.623 K，标准差为0.57 K和相关系数为0.998。从表1和表2可以看出，热红外波段增加可以提高反演精度。
热红外波段加水汽（W）波段组合，主要适合于白天，因为水汽是从近红外波段估算得到。表3是MODIS波段29-31-32和水汽参数（W）作为深度学习神经网络的输入参数是反演误差部分信息表，从表3中可以看出当隐含层数是4和隐含节点是600-600-600-600时，反演精度最高。平均最高精度是0.251 K，标准差是0.255 K，相关系数是1。
表4是MODIS波段29-31-32-33与W作为输入参数时的部分反演误差信息表，当隐含层数是4层并且隐含节点为500-500-500-500时精度最高，平均精度最高为0.302 K，标准差为0.273 K，相关系数是0.999。从表3和表4可以看出，当有水汽波段时，添加波段33时并没有提高反演精度。可以看出，水汽波段在地表温度反演中起到作用非常大。因此为提高地表温度反演精度，尽量添加大气水汽信息。
中红波段20/22/23（3~5 μm）与热红外波段29/31/32/33（8~13.5 μm）的组合，由于中红外波段白天受太阳的干扰，因此中红外更适合在晚上反演地表温度。表5是波段20-31-32组合的反演误差信息表，当隐含层数是6层和隐含节点是600-600-600-600-600-600精度最高，为0.903 K，标准差为0.934 K，相关系数0.995。
表6是波段20-29-31-32组合时的反演误差，最高精度是在隐含层数为3和隐含节点为600-600-600时，最高平均精度为0.568 K，标准偏差为0.608 K，相关系数为0.998。
表7是波段20-22-31-32组合时的反演误差，最高精度是在隐含层数为4层和隐含节点为700-700-700-700时，最高平均精度为0.7 K，标准偏差为0.737 K，相关系数为0.997。对比表6和表7可知，利用热红外波段29比中红外波段22反演时候精度要高。
表8是波段20-22-29-31-32组合时的反演误差，最高精度是在隐含层数为6层和隐含节点为500-500-500-500-500-500时，最高平均精度为0.509 K，标准偏差为0.566 K，相关系数为0.998。
表9是波段20-22-23-29-31-32组合时的反演误差，最高精度是在隐含层数为5层和隐含节点为500-500-500-500-500时，最高平均精度为0.479 K，标准偏差为0.558 K，相关系数为0.998。
表10是波段20-22-29-31-32-33组合时的反演误差，最高精度是在隐含层数为5层和隐含节点为400-400-400-400-400时，最高平均精度为0.626 K，标准偏差为0.551 K，相关系数为0.998。对比表9和表10可以看出，利用中外波段23比利用红外波段33反演精度要高。
表11是波段20-22-23-29-31-32-33组合时的反演误差，最高精度是在隐含层数为4层和隐含节点为700-700-700-700时，最高平均精度为0.478 K，标准偏差为0.444 K，相关系数为0.999。对比表9和表11可知，利用7个波段（3个中红外波段+4个红外波段）与6个波段（3个中外波段+3个红外波段）的反演精度相当。由于第33波段受二氧化碳的影响比较大，从计算量角度讲，利用20-22-23-29-31-32波段组合更合适。
3 结论
虽然地表温度反演是属于典型的病态问题，但地球物理参数之间存在着相互联系和相互约束的关系。这些关系不能严格地用数学方法来描述，这就决定了大气辐射传输模型和深度动态学习神经网络的集成是解决地球物理参数（地表温度）病态反演问题的最好方法之一。以MODIS中红外波段和热红外波段作为分析对象，利用MODTRAN对MODIS波段20/22/23（3~5 μm）和波段29/31/32/33（8~13.5 μm）数据的白天/晚上、白天、晚上进行了模拟，并利用深度动态学习神经网络进行了反演分析。
对于适合白天和晚上的反演组合，波段29-31-32-33组合比波段29-31-32组合反演精度要高，说明热红外波段增加可以提高反演精度；对于只适合白天反演组合，波段29-31-32与W组合比波段29-31-32-33与W组合精度更高，反演误差都在0.3 K以下，说明水汽参数做温度反演中起的作用非常大，另外当有水汽波段和3个热红外波段时，再增加热红外波段33不会提高反演精度；在利用中红波段20/22/23（3~5 μm）与热红外波段29/31/32/33（8~13.5 μm）反演晚上地表温度的组合中，波段20-22-23-29-31-32组合和波段20-22-23-29-31-32-33组合时的反演误差都在0.5 K以下，7个波段组合精度略高。由于第33波段受二氧化碳的影响比较大，从计算量角度讲，利用20-22-23-29-31-32波段组合更合适。通过分析发现，利用动态学习神经网络和辐射传输模型可以高精度地从单景MODIS数据中反演地表温度和发射率，并且克服美国NASA白天/晚上算法的缺陷，彻底解决了地表温度和发射率反演的难题。深度动态学习神经网络和辐射传输模型的结合在地表温度和发射率反演的研究成果为国产风云和高分热红外传感器地表温度反演提供了参考模式。
参考文献
- 1
Price J. C. Land surface temperature measurements from the split-window channels of the NOAA-7 AVHRR. Journal of Geophysical Research , 1984, 79: 5039~5044.
- 2
Becker F., Li Z-L. Towards a local split window method over land surface. International Journal of Remote Sensing , 1990, 11: 369~393.
- 3
Sobrino J. A., V. Caselles . A methodology for obtaining the crop temperature from NOAA-9 AVHRR data. International Journal of Remote Sensing , 1991, I2: 2461~2475.
- 4
Sobrino J. A., Li Z-L., Stoll M. P., et al. Improvements in the split window technique for land surface temperature determination. IEEE Transaction on Geoscience and Remote sensing , 1994, 32: 243~253.
- 5
Coll C., V. Caselles , A. Sobrino , et al. On the atmospheric dependence of the split-window equation for land surface temperature. International Journal of Remote sensing Environment , 1994, 27: 105~122.
- 6
Vidal A.. Atmosphere and emissivity correction of land surface temperature measured from satellite data. International Journal of Remote Sensing , 1991, 12: 2449~2460.
- 7
Kerr Y. H, J. P. Lagouarde, J. Imbernon. Accurate land surface temperature retrieval from AVHRR data with use of an improved split window algorithm. Remote Sensing of Environment , 1992, 41: 197~209.
- 8
Otlle C., M. Stoll . Effect of atmospheric absorption and surface emissivity on the determination of land temperature from infrared satellite data. International Journal of Remote Sensing , 1993, 14: 2025~2037.
- 9
Prata. A. J.. Land surface temperatures derived from the advanced very high resolution radiometer and the along-track scanning radiometer 2. Experimental results andvalidation of AVHRR algorithms. Journal of Geophysical Research Atmospheres , 1993, 99: 13025~13058.
- 10
Wan Z., J. Dozier . A generalized split-window algorithm for retrieving land surface temperature measurement from space. IEEE Transaction on Geoscience and Remote sensing , 1996, 34: 892~905.
- 11
Qin Zhihao, Giorgio Dall Olmo, Arnon Karnieli . Derivation of split window algorithm and its sensitivity analysis for retrieving land surface temperature from NOAA-advanced very high resolution radiometer data. Geophysical research , 2001, 106(D19): 22655~22670.
- 12
Sobrino J. A., J. El Kharraz , Z. L. Li . Surface temperature and water vapour retrieval from MODIS data. International Journal of Remote Sensing , 2004, 24: 5161~5182.
- 13
Mao K., Z. Qin , J. Shi , et al. A practical split-window algorithm for retriving land surface temperature from MODIS data. International Journal of Remote Sensing , 2005, 15: 3181~3204.
- 14
Lang Xia, Kebiao Mao, Y. Ma , et al. An algorithm for retrieving land surface temperature using VIIRS data in combination with multi-sensors. Sensors , 2014, 14: 21385~21408.
- 15
Kahle A. B., Alley R. E. . Separation of temperature and emittance in remotely sensed radiance measurements. Remote Sensing of Environment , 1992, 42: 1~20.
- 16
Hook S. J., Gabell A. R., Green A. A., et al. A comparison of techniques for extracting emissivity information from thermal infrared data for geologic studies. Remote Sensing of Environment , 1992, 42: 123~135.
- 17
Watson K. . Spectral ratio method for measuring emissivity. Remote Sensing of Environment , 1992, 42: 113~116.
- 18
Kealy P. S., Hook S. . Separating temperature and emissivity in thermal infrared multispectral scanner data: Implication for recovering land surface temperatures. Transaction on Geoscience and Remote Sensing , 1993, 31: 1155~1164.
- 19
Schmugge T., Hook S. J., Coll C. Recovering surface temperature and emissivity from thermal infrared multispectral data. Remote Sensing of Environment , 1998, 65(2): 121~131.
- 20
Schmugge T., French A., Ritchie J. C., et al. Temperature and emissivity separation from multispectral thermal infrared observations, Remote sensing of Environment , 2002, 79(2~3): 189~198.
- 21
Liang Shunlin. An optimization algorithm for separating land surface temperature and emissivity from multispectral thermal infrared imagery. IEEE Transaction on Geoscience and Remote Sensing , 2001, 39: 264~274.
- 22
Li Z-L, Becker F. . Feasibility of land surface temperature and emissivity determination from AVHRR data. Remote Sensing of Environment , 1993, 43: 67~85.
- 23
Wan Zhengming, Zhao-liang Li. A physics-based algorithm for retrieving land-surface emissivity and temperature from EOS/MODIS data. IEEE Transaction on Geoscience and Remote Sensing , 1997, 35: 980~996.
- 24
Kebiao Mao, SanMei Li, DaoLong Wang, et al. Retrieval of land surface temperature and emissivity from ASTER1B data using dynamic learning neural network. International Journal of Remote Sensing , 2011, 32(19): 5413~5423.
- 25
毛克彪, 唐华俊, 陈仲新, 等. 一个用神经网络优化的针对ASTER数据反演地表温度和发射率的多波段算法. 国土资源遥感 , 2007, 73 (3): 18~22.
- 26
毛克彪, 唐华俊, 李丽英, 等. 一个从MODIS数据同时反演地表温度和发射率的神经网络算法. 遥感信息 , 2007, 92(4): 9~15.
- 27
毛克彪, 施建成, 覃志豪, 等. 一个针对ASTER数据同时反演地表温度和比辐射率的四通道算法. 遥感学报 , 2006, 4: 593~599.
- 28
Kebiao Mao, Jiancheng Shi, Huajun Tang, et al. A neural network technique for separating land surface emissivity and temperature from ASTER imagery. IEEE Transation on Geoscience and Remote Sensing , 2008, 46(1): 200~208.
- 29
Kebiao Mao, J. Shi , Z. Li , et al. An RM-NN algorithm for retrieving land surface temperature and emissivity from EOS/MODIS data. Journal of Geophysical Research-atmosphere , 2007, 112(D21): 1~17.
- 30
Hsu Shinyi , T. Masters , M. Olson , et al. Comparavtive analysis of five neural networks models. Remote Sensing Reviews , 1992, 6: 319~329.
- 31
Yu Chang Tzeng, K. S. Chen , Wen-Liang, et al. , A Dynamic learning nerual network for remote sensing applications. IEEE Transation on Geoscience and Remote Sensing, 1994, 32: 1096~1102.
- 1
基本信息
DOI:10.12105/j.issn.1672-0423.20180506
引用信息
基金信息
稿件历史
纸刊出版
: 2018-10-25
收稿日期
: 2018-09-08
角 色：第一作者
Role：First author
作者简介：毛克彪（1977—），男，汉族，湖南沅江人，博士、研究员。研究方向：农业遥感。Email：maokebiao@caas.cn
| 29-31-32 | 隐含节点 | |||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 隐含层 | 300 | 400 | 500 | 600 | 700 | |||||||||||
| M | SD | R | M | SD | R | M | SD | R | M | SD | R | M | SD | R | ||
| 3 | 0.956 | 0.879 | 0.995 | 0.984 | 0.865 | 0.995 | 0.953 | 0.894 | 0.995 | 0.925 | 0.873 | 0.995 | 0.941 | 0.901 | 0.995 | |
| 4 | 0.967 | 1.083 | 0.994 | 0.93 | 0.863 | 0.995 | 0.917 | 0.85 | 0.995 | 1.034 | 2.18 | 0.982 | 0.926 | 0.897 | 0.995 | |
| 5 | 0.907 | 0.85 | 0.995 | 0.956 | 1.226 | 0.993 | 0.941 | 1.032 | 0.994 | 0.977 | 1.265 | 0.992 | 0.937 | 1.366 | 0.992 | |
| 6 | 0.947 | 1.039 | 0.994 | 1.123 | 4.04 | 0.95 | 1.123 | 6.05 | 0.9 | 4.421 | 81.11 | 0.149 | 3.222 | 34.628 | 0.354 | |
| 29-31-32-33 | 隐含节点 | ||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 隐含层 | 300 | 400 | 500 | 600 | 700 | ||||||||||
| M | SD | R | M | SD | R | M | SD | R | M | SD | R | M | SD | R | |
| 3 | 0.703 | 0.593 | 0.997 | 0.646 | 0.568 | 0.998 | 0.623 | 0.57 | 0.998 | 0.655 | 0.591 | 0.998 | 0.651 | 0.592 | 0.998 |
| 4 | 0.705 | 0.664 | 0.997 | 0.857 | 3.267 | 0.967 | 0.713 | 0.733 | 0.997 | 0.722 | 0.835 | 0.996 | 0.702 | 0.812 | 0.996 |
| 5 | 0.673 | 0.636 | 0.997 | 0.678 | 0.665 | 0.997 | 0.707 | 0.684 | 0.997 | 0.802 | 1.429 | 0.992 | 0.764 | 0.721 | 0.997 |
| 6 | 0.885 | 3.249 | 0.966 | 0.717 | 0.715 | 0.997 | 0.781 | 1.032 | 0.995 | 0.922 | 2.764 | 0.975 | 1.223 | 7.457 | 0.855 |
| 29-31-32-W | 隐含节点 | ||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 隐含层 | 300 | 400 | 500 | 600 | 700 | ||||||||||
| M | SD | R | M | SD | R | M | SD | R | M | SD | R | M | SD | R | |
| 3 | 0.283 | 0.276 | 1 | 0.33 | 0.271 | 0.999 | 0.273 | 0.274 | 1 | 0.289 | 0.281 | 0.999 | 0.271 | 0.273 | 1 |
| 4 | 0.294 | 0.273 | 1 | 0.29 | 0.279 | 0.999 | 0.308 | 0.723 | 0.998 | 0.251 | 0.255 | 1 | 0.259 | 0.265 | 1 |
| 5 | 0.285 | 0.293 | 0.999 | 0.291 | 0.369 | 0.999 | 0.34 | 2.769 | 0.999 | 0.917 | 8.877 | 0.815 | 0.546 | 4.605 | 0.938 |
| 6 | 0.351 | 0.313 | 0.999 | 0.409 | 2.548 | 0.98 | 0.329 | 1.087 | 0.996 | 0.349 | 0.637 | 0.998 | 0.345 | 0.68 | 0.998 |
| 29-30-31-32-W | 隐含节点 | |||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 隐含层 | 300 | 400 | 500 | 600 | 700 | |||||||||||
| M | SD | R | M | SD | R | M | SD | R | M | SD | R | M | SD | R | ||
| 3 | 0.38 | 0.318 | 0.999 | 0.332 | 0.26 | 0.999 | 0.403 | 0.316 | 0.999 | 0.331 | 0.27 | 0.999 | 0.417 | 2.763 | 0.999 | |
| 4 | 0.426 | 0.307 | 0.999 | 0.341 | 0.279 | 0.999 | 0.302 | 0.273 | 0.999 | 0.337 | 0.286 | 0.999 | 0.338 | 0.552 | 0.999 | |
| 5 | 0.488 | 0.901 | 0.997 | 0.319 | 0.289 | 0.999 | 0.372 | 0.347 | 0.999 | 0.3 | 0.28 | 0.999 | 0.332 | 0.284 | 0.999 | |
| 6 | 0.388 | 0.331 | 0.999 | 0.723 | 9.754 | 0.796 | 0.314 | 0.407 | 0.999 | 0.347 | 0.48 | 0.999 | 0.561 | 2.724 | 0.977 | |
| 20-31-32 | 隐含节点 | ||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 隐含层 | 300 | 400 | 500 | 600 | 700 | ||||||||||
| M | SD | R | M | SD | R | M | SD | R | M | SD | R | M | SD | R | |
| 3 | 1.015 | 0.928 | 0.994 | 0.973 | 0.921 | 0.995 | 0.977 | 0.928 | 0.994 | 1 | 0.932 | 0.994 | 0.964 | 0.907 | 0.995 |
| 4 | 0.98 | 0.932 | 0.994 | 0.958 | 0.912 | 0.995 | 0.957 | 0.906 | 0.995 | 0.952 | 0.906 | 0.995 | 0.963 | 0.906 | 0.995 |
| 5 | 0.987 | 0.928 | 0.994 | 0.976 | 0.933 | 0.994 | 0.899 | 0.909 | 0.995 | 0.959 | 1.184 | 0.993 | 0.927 | 0.907 | 0.995 |
| 6 | 1.129 | 0.975 | 0.993 | 0.929 | 0.918 | 0.995 | 0.936 | 0.939 | 0.995 | 0.903 | 0.934 | 0.995 | 0.925 | 1.03 | 0.994 |
| 20-29-31-32 | 隐含节点 | ||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 隐含层 | 300 | 400 | 500 | 600 | 700 | ||||||||||
| M | SD | R | M | SD | R | M | SD | R | M | SD | R | M | SD | R | |
| 3 | 1.006 | 0.928 | 0.994 | 0.844 | 0.78 | 0.996 | 0.666 | 0.673 | 0.997 | 0.568 | 0.608 | 0.998 | 0.581 | 0.61 | 0.998 |
| 4 | 0.642 | 0.656 | 0.997 | 0.607 | 0.624 | 0.998 | 0.638 | 0.662 | 0.997 | 0.58 | 0.609 | 0.998 | 0.611 | 0.648 | 0.998 |
| 5 | 1.022 | 0.866 | 0.994 | 0.666 | 0.668 | 0.997 | 1.149 | 17.11 | 0.579 | 0.602 | 1.039 | 0.996 | 0.619 | 0.828 | 0.997 |
| 6 | 0.694 | 1.512 | 0.991 | 0.514 | 0.568 | 0.998 | 0.8 | 2.572 | 0.978 | 0.763 | 2.067 | 0.985 | 0.571 | 1.37 | 0.993 |
| 20-22-31-32 | 隐含节点 | ||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 隐含层 | 300 | 400 | 500 | 600 | 700 | ||||||||||
| M | SD | R | M | SD | R | M | SD | R | M | SD | R | M | SD | R | |
| 3 | 0.838 | 0.821 | 0.996 | 1.204 | 0.954 | 0.993 | 0.782 | 0.809 | 0.996 | 0.776 | 0.813 | 0.996 | 0.719 | 0.794 | 0.996 |
| 4 | 0.857 | 0.813 | 0.996 | 0.781 | 0.812 | 0.996 | 0.754 | 0.808 | 0.996 | 0.714 | 0.775 | 0.997 | 0.7 | 0.737 | 0.997 |
| 5 | 0.914 | 0.824 | 0.995 | 0.785 | 0.834 | 0.996 | 0.779 | 0.777 | 0.996 | 0.714 | 0.75 | 0.997 | 0.966 | 4.47 | 0.94 |
| 6 | 0.944 | 0.876 | 0.995 | 0.709 | 0.789 | 0.997 | 0.76 | 0.89 | 0.996 | 0.718 | 1.091 | 0.995 | 0.765 | 0.968 | 0.995 |
| 20-22-29-31-32 | 隐含节点 | |||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 隐含层 | 300 | 400 | 500 | 600 | 700 | |||||||||||
| M | SD | R | M | SD | R | M | SD | R | M | SD | R | M | SD | R | ||
| 3 | 0.803 | 0.682 | 0.997 | 0.644 | 0.598 | 0.998 | 0.638 | 0.675 | 0.997 | 0.539 | 0.537 | 0.998 | 0.734 | 0.656 | 0.997 | |
| 4 | 0.716 | 0.733 | 0.997 | 0.584 | 0.568 | 0.998 | 0.749 | 2.463 | 0.98 | 0.563 | 0.586 | 0.998 | 0.694 | 3.075 | 0.97 | |
| 5 | 0.765 | 2.469 | 0.98 | 0.621 | 0.621 | 0.998 | 0.668 | 1.286 | 0.994 | 0.584 | 0.709 | 0.997 | 0.705 | 2.254 | 0.983 | |
| 6 | 0.957 | 0.895 | 0.995 | 0.928 | 7.127 | 0.867 | 0.509 | 0.566 | 0.998 | 0.607 | 1.391 | 0.993 | 1.26 | 10.121 | 0.774 | |
| 20-22-23-29-31-32 | 隐含节点 | |||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 隐含层 | 300 | 400 | 500 | 600 | 700 | |||||||||||
| M | SD | R | M | SD | R | M | SD | R | M | SD | R | M | SD | R | ||
| 3 | 0.763 | 1.38 | 0.992 | 0.717 | 0.646 | 0.997 | 0.565 | 0.796 | 0.997 | 0.647 | 0.724 | 0.997 | 0.385 | 0.43 | 0.999 | |
| 4 | 0.531 | 0.514 | 0.998 | 0.938 | 5.426 | 0.915 | 0.716 | 3.002 | 0.972 | 0.61 | 1.035 | 0.996 | 0.563 | 1.617 | 0.991 | |
| 5 | 1.037 | 7.441 | 0.861 | 0.544 | 0.793 | 0.997 | 0.479 | 0.558 | 0.998 | 0.434 | 0.642 | 0.998 | 0.73 | 4.54 | 0.94 | |
| 6 | 0.89 | 0.816 | 0.995 | 0.634 | 0.868 | 0.996 | 0.636 | 4.363 | 0.945 | 0.437 | 0.607 | 0.998 | 0.558 | 3.967 | 0.953 | |
| 20-22-29-31-32-33 | 隐含节点 | ||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 隐含层 | 300 | 400 | 500 | 600 | 700 | ||||||||||
| M | SD | R | M | SD | R | M | SD | R | M | SD | R | M | SD | R | |
| 3 | 0.838 | 0.672 | 0.996 | 0.734 | 0.619 | 0.997 | 0.706 | 0.587 | 0.997 | 0.622 | 0.527 | 0.998 | 0.63 | 0.527 | 0.998 |
| 4 | 0.69 | 0.606 | 0.997 | 0.777 | 0.667 | 0.997 | 0.674 | 0.575 | 0.998 | 0.617 | 0.526 | 0.998 | 0.702 | 0.8 | 0.997 |
| 5 | 0.751 | 0.689 | 0.997 | 0.626 | 0.551 | 0.998 | 0.956 | 3.622 | 0.959 | 0.888 | 1.954 | 0.986 | 0.847 | 1.396 | 0.992 |
| 6 | 0.783 | 0.854 | 0.996 | 0.755 | 0.735 | 0.997 | 0.585 | 0.563 | 0.998 | 0.562 | 0.564 | 0.998 | 0.643 | 0.783 | 0.997 |
| 20-22-23-29-31-32-33 | 隐含节点 | ||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 隐含层 | 300 | 400 | 500 | 600 | 700 | ||||||||||
| M | SD | R | M | SD | R | M | SD | R | M | SD | R | M | SD | R | |
| 3 | 0.804 | 0.652 | 0.997 | 0.736 | 0.615 | 0.997 | 0.717 | 0.603 | 0.997 | 0.596 | 0.532 | 0.998 | 0.524 | 0.487 | 0.998 |
| 4 | 0.814 | 0.65 | 0.997 | 0.605 | 0.511 | 0.998 | 0.869 | 1.715 | 0.989 | 0.557 | 0.534 | 0.998 | 0.478 | 0.444 | 0.999 |
| 5 | 0.753 | 0.64 | 0.997 | 0.599 | 0.537 | 0.998 | 0.736 | 0.701 | 0.997 | 0.653 | 1.562 | 0.991 | 0.635 | 0.748 | 0.997 |
| 6 | 0.757 | 0.678 | 0.997 | 1.47 | 24.17 | 0.47 | 0.704 | 0.606 | 0.997 | 0.642 | 0.533 | 0.998 | 0.699 | 0.674 | 0.997 |
表1 波段29-31-32组合时反演误差
Table 1 The Summery of Retrieval Error for band 29-31-32 combination
表2 波段29-31-32-33组合时反演误差
Table 2 The Summery of Retrieval Error for band 29-31-32-33 combination
表3 波段29-31-32-W组合时反演误差
Table 3 The Summery of Retrieval Error for band 29-31-32-W combination
表4 波段29-31-32-33-W组合时反演误差
Table 4 The Summery of Retrieval Error for band 29-31-32-33-W combination
表5 波段20-31-32组合时反演误差
Table 5 The Summery of Retrieval Error for band 20-31-32 combination
表6 波段20-29-31-32组合时反演误差
Table 6 The Summery of Retrieval Error for band 20-29-31-32 combination
表7 波段20-22-31-32组合时反演误差
Table 7 The Summery of Retrieval Error for band 20-22-31-32 combination
表8 波段20-22-29-31-32组合时反演误差
Table 8 The Summery of Retrieval Error for band 20-22-29-31-32 combination
表9 波段20-22-23-29-31-32组合时反演误差
Table 9 The Summery of Retrieval Error for band 20-22-23-29-31-32 combination
表10 波段20-22-29-31-32-33组合时反演误差
Table 10 The Summery of Retrieval Error for band 20-22-29-31-32-33 combination
表11 波段20-22-23-29-31-32-33组合时反演误差
Table 11 The Summery of Retrieval Error for band 20-22-23-29-31-32-33 combination
其中：M为平均误差；R为相关系数；SD为标准偏差，以下同
无注解
W为大气水汽含量
无注解
无注解
无注解
无注解
无注解
无注解
无注解
无注解
参考文献
- 1
Price J. C. Land surface temperature measurements from the split-window channels of the NOAA-7 AVHRR. Journal of Geophysical Research , 1984, 79: 5039~5044.
- 2
Becker F., Li Z-L. Towards a local split window method over land surface. International Journal of Remote Sensing , 1990, 11: 369~393.
- 3
Sobrino J. A., V. Caselles . A methodology for obtaining the crop temperature from NOAA-9 AVHRR data. International Journal of Remote Sensing , 1991, I2: 2461~2475.
- 4
Sobrino J. A., Li Z-L., Stoll M. P., et al. Improvements in the split window technique for land surface temperature determination. IEEE Transaction on Geoscience and Remote sensing , 1994, 32: 243~253.
- 5
Coll C., V. Caselles , A. Sobrino , et al. On the atmospheric dependence of the split-window equation for land surface temperature. International Journal of Remote sensing Environment , 1994, 27: 105~122.
- 6
Vidal A.. Atmosphere and emissivity correction of land surface temperature measured from satellite data. International Journal of Remote Sensing , 1991, 12: 2449~2460.
- 7
Kerr Y. H, J. P. Lagouarde, J. Imbernon. Accurate land surface temperature retrieval from AVHRR data with use of an improved split window algorithm. Remote Sensing of Environment , 1992, 41: 197~209.
- 8
Otlle C., M. Stoll . Effect of atmospheric absorption and surface emissivity on the determination of land temperature from infrared satellite data. International Journal of Remote Sensing , 1993, 14: 2025~2037.
- 9
Prata. A. J.. Land surface temperatures derived from the advanced very high resolution radiometer and the along-track scanning radiometer 2. Experimental results andvalidation of AVHRR algorithms. Journal of Geophysical Research Atmospheres , 1993, 99: 13025~13058.
- 10
Wan Z., J. Dozier . A generalized split-window algorithm for retrieving land surface temperature measurement from space. IEEE Transaction on Geoscience and Remote sensing , 1996, 34: 892~905.
- 11
Qin Zhihao, Giorgio Dall Olmo, Arnon Karnieli . Derivation of split window algorithm and its sensitivity analysis for retrieving land surface temperature from NOAA-advanced very high resolution radiometer data. Geophysical research , 2001, 106(D19): 22655~22670.
- 12
Sobrino J. A., J. El Kharraz , Z. L. Li . Surface temperature and water vapour retrieval from MODIS data. International Journal of Remote Sensing , 2004, 24: 5161~5182.
- 13
Mao K., Z. Qin , J. Shi , et al. A practical split-window algorithm for retriving land surface temperature from MODIS data. International Journal of Remote Sensing , 2005, 15: 3181~3204.
- 14
Lang Xia, Kebiao Mao, Y. Ma , et al. An algorithm for retrieving land surface temperature using VIIRS data in combination with multi-sensors. Sensors , 2014, 14: 21385~21408.
- 15
Kahle A. B., Alley R. E. . Separation of temperature and emittance in remotely sensed radiance measurements. Remote Sensing of Environment , 1992, 42: 1~20.
- 16
Hook S. J., Gabell A. R., Green A. A., et al. A comparison of techniques for extracting emissivity information from thermal infrared data for geologic studies. Remote Sensing of Environment , 1992, 42: 123~135.
- 17
Watson K. . Spectral ratio method for measuring emissivity. Remote Sensing of Environment , 1992, 42: 113~116.
- 18
Kealy P. S., Hook S. . Separating temperature and emissivity in thermal infrared multispectral scanner data: Implication for recovering land surface temperatures. Transaction on Geoscience and Remote Sensing , 1993, 31: 1155~1164.
- 19
Schmugge T., Hook S. J., Coll C. Recovering surface temperature and emissivity from thermal infrared multispectral data. Remote Sensing of Environment , 1998, 65(2): 121~131.
- 20
Schmugge T., French A., Ritchie J. C., et al. Temperature and emissivity separation from multispectral thermal infrared observations, Remote sensing of Environment , 2002, 79(2~3): 189~198.
- 21
Liang Shunlin. An optimization algorithm for separating land surface temperature and emissivity from multispectral thermal infrared imagery. IEEE Transaction on Geoscience and Remote Sensing , 2001, 39: 264~274.
- 22
Li Z-L, Becker F. . Feasibility of land surface temperature and emissivity determination from AVHRR data. Remote Sensing of Environment , 1993, 43: 67~85.
- 23
Wan Zhengming, Zhao-liang Li. A physics-based algorithm for retrieving land-surface emissivity and temperature from EOS/MODIS data. IEEE Transaction on Geoscience and Remote Sensing , 1997, 35: 980~996.
- 24
Kebiao Mao, SanMei Li, DaoLong Wang, et al. Retrieval of land surface temperature and emissivity from ASTER1B data using dynamic learning neural network. International Journal of Remote Sensing , 2011, 32(19): 5413~5423.
- 25
毛克彪, 唐华俊, 陈仲新, 等. 一个用神经网络优化的针对ASTER数据反演地表温度和发射率的多波段算法. 国土资源遥感 , 2007, 73 (3): 18~22.
- 26
毛克彪, 唐华俊, 李丽英, 等. 一个从MODIS数据同时反演地表温度和发射率的神经网络算法. 遥感信息 , 2007, 92(4): 9~15.
- 27
毛克彪, 施建成, 覃志豪, 等. 一个针对ASTER数据同时反演地表温度和比辐射率的四通道算法. 遥感学报 , 2006, 4: 593~599.
- 28
Kebiao Mao, Jiancheng Shi, Huajun Tang, et al. A neural network technique for separating land surface emissivity and temperature from ASTER imagery. IEEE Transation on Geoscience and Remote Sensing , 2008, 46(1): 200~208.
- 29
Kebiao Mao, J. Shi , Z. Li , et al. An RM-NN algorithm for retrieving land surface temperature and emissivity from EOS/MODIS data. Journal of Geophysical Research-atmosphere , 2007, 112(D21): 1~17.
- 30
Hsu Shinyi , T. Masters , M. Olson , et al. Comparavtive analysis of five neural networks models. Remote Sensing Reviews , 1992, 6: 319~329.
- 31
Yu Chang Tzeng, K. S. Chen , Wen-Liang, et al. , A Dynamic learning nerual network for remote sensing applications. IEEE Transation on Geoscience and Remote Sensing, 1994, 32: 1096~1102.