地球科学进展

地球科学进展 ›› 2018, Vol. 33 ›› Issue (6): 641 -652. doi: 10.11867/j.issn.1001-8166.2018.06.0641

研究论文 上一篇    下一篇

GF-1 WFV2传感器数据的缨帽变换系数反演
王帅( ), 徐涵秋 *( ), 施婷婷   
  1. 1.福州大学环境与资源学院,福建 福州 350116
    2.福州大学遥感信息工程研究所,福建 福州 350116
    3.福州大学福建省水土流失遥感监测评估重点实验室,福建 福州 350116
  • 收稿日期:2017-10-27 修回日期:2018-05-16 出版日期:2018-06-20
  • 通讯作者: 徐涵秋 E-mail:wswin8@hotmail.com;hxu@fzu.edu.cn
  • 基金资助:
    *国家自然科学基金项目“星—地协同多源遥感数据不透水面辨识机理与反演模型研究”(编号:41501469);福建省测绘地理信息局项目“国产高分辨率卫星数据的交互比较研究”(编号:2017JX02)资助.

Retrieval of Tasseled Cap Transformation Coefficients for GF-1 WFV2 Sensor Data

Shuai Wang( ), Hanqiu Xu *( ), Tingting Shi   

  1. 1.College of Environment and Resources, Fuzhou University, Fuzhou 350116, China
    2.Institute of Remote Sensing Information Engineering, Fuzhou University, Fuzhou 350116, China
    3.Fujian Provinicial Key Laboratory of Remote Sensing of Soil Erosion, Fuzhou University,Fuzhou 350116, China
  • Received:2017-10-27 Revised:2018-05-16 Online:2018-06-20 Published:2018-07-23
  • Contact: Hanqiu Xu E-mail:wswin8@hotmail.com;hxu@fzu.edu.cn
  • About author:

    First author:Wang Shuai (1992-), male, Fuyang City, Anhui Province, Master student. Research areas include remote sensing of resources and environment. E-mail:wswin8@hotmail.com

  • Supported by:
    Project supported by the National Natural Science Foundation of China “Study on identification mechanism and inversion model of impervious surface for satellite ground collaborative multi-source remote sensing data”(No.41501469);The Science and Technology Development Fund of Fujian Provincial Administration of Surveying, Mapping and Geoinformation “Study on the cross comparison of domestic high-resolution satellite images” (No.2017JX02).
全文 ( 29 ) 下载PDF ( 1251 )

针对国产GF-1 WFV2传感器,提出一种缨帽变换系数的推导方法。利用全国范围内不同地区的18幅GF-1影像和5幅同步的Landsat 8影像,选择了大量覆盖不同地物类型的典型样本并分析光谱变化特征。在此基础上,用3幅Landsat 8同步影像对湿度分量回归,然后再用施密特正交化方法依次推导出两两正交的亮度分量和绿度分量。结果显示,推导的前2个分量能解释缨帽变换总变化量的97%以上,前3个分量能解释总变化量的99%以上。针对不同区域影像的实验结果表明,各分量都具有很好的稳定性,而不同分量之间的散点图也都有相似的分布特征,说明该系数有很好的适用性。与验证用的Landsat 8同步影像对比后发现,推导的缨帽变换结果与Landsat 8相应的分量之间都有很好的相关性和一致性,尽管2种传感器波段设置及光谱响应能力不同,但相关性仍高达0.9。对GF-1 WFV1,WFV3和WFV4传感器影像的实验结果表明,针对GF-1 WFV2传感器提出的缨帽变换系数,也可以很好地应用在WFV3和WFV1传感器上,但不建议在WFV4传感器上使用。

关键词: GF-1, 缨帽变换, 亮度, 绿度, 湿度

Aiming at GF-1 WFV2 sensor, a method for deriving Tasseled Cap Transformation (TCT) coefficients was proposed in this paper. Based on 18 GF-1 images and 5 synchronous Landsat 8 images of different regions throughout the country, a large number of typical samples covering different ground types were selected, and the spectral characteristics of these samples were then analyzed. On this basis, the wetness component of the TCT was retrieved by the regression using three Landsat 8 and WFV2 synchronous images. The orthogonal brightness and greenness components were then derived with the Gram-Schmidt orthogonalization method. The results showed that the derived first two components can explain more than 97% of total variation, and the first three components can explain more than 99% of the total variation. The experimental results from different regions showed that each component had good stability, and the scatter-plots of different components also exhibited similar distribution patterns. This indicates that the derived coefficients have good applicability. Compared with the validation images, the results of TCT are well correlated with corresponding Landsat 8 components. Although these two sensors have different band settings and different spectral responses, their correlation coefficient is still as high as 0.9. The application results for GF-1 WFV1, WFV3 and WFV4 showed that the coefficients derived from WFV2 sensor can also be used for WFV3 and WFV1 sensors, but they are not recommended for WFV4 sensor.

Key words: GF-1, Tasseled Cap Transformation, Brightness, Greenness, Wetness.

中图分类号: 

表1 数据源列表
Table 1 Source lists of data
所属
区域		 地点 成像
日期		 GF-1
产品号		 Landsat 8
列/行
实验
影像
地区		 中部
地区		 驻马店、信阳 20151020 1116843
开封、商丘、亳州 20160504 1562284 123/36
驻马店、信阳、阜阳 20160908 1562295 123/37
驻马店、信阳、阜阳 20160908 1813878
东部
地区		 日照市 20160429 1552922 120/36
南京、滁州 20160429 1552923
南方
地区		 柳州、来宾 20161216 2047782
玉林、贵港 20170401 2277357
柳州、贺州 20170401 2277356
西部
地区		 乌兰乌拉湖,西北 20160721 1712329
乌兰乌拉湖,西北 20161117 1970495
东北地区 营口市 20160507 1569457
验证
影像
地区		 中部
地区		 聊城、濮阳 20160504 1562286 123/35
驻马店、信阳、阜阳 20160516 1586464
东部地区 连云港、淮安 20160429 1552914 120/37
南方地区 贺州、梧州 20161228 2093687
西部地区 乌兰乌拉湖,西南 20161117 1970493
东北地区 营口市 20160919 1839479
表1 数据源列表
Table 1 Source lists of data
所属
区域		 地点 成像
日期		 GF-1
产品号		 Landsat 8
列/行
实验
影像
地区		 中部
地区		 驻马店、信阳 20151020 1116843
开封、商丘、亳州 20160504 1562284 123/36
驻马店、信阳、阜阳 20160908 1562295 123/37
驻马店、信阳、阜阳 20160908 1813878
东部
地区		 日照市 20160429 1552922 120/36
南京、滁州 20160429 1552923
南方
地区		 柳州、来宾 20161216 2047782
玉林、贵港 20170401 2277357
柳州、贺州 20170401 2277356
西部
地区		 乌兰乌拉湖,西北 20160721 1712329
乌兰乌拉湖,西北 20161117 1970495
东北地区 营口市 20160507 1569457
验证
影像
地区		 中部
地区		 聊城、濮阳 20160504 1562286 123/35
驻马店、信阳、阜阳 20160516 1586464
东部地区 连云港、淮安 20160429 1552914 120/37
南方地区 贺州、梧州 20161228 2093687
西部地区 乌兰乌拉湖,西南 20161117 1970493
东北地区 营口市 20160919 1839479
图1 影像成像区分布
Fig.1 Distribution of imaging regions
图1 影像成像区分布
Fig.1 Distribution of imaging regions
图2 不同地物类型的光谱特征
(a)土壤光谱特征1;(b)土壤光谱特征2;(c)土壤光谱特征3;(d)水体光谱特征;(e)生长期植被光谱特征;(f)成熟期植被光谱特征
Fig.2 Spectral signature of different ground types
(a)Soil spectral signature 1;(b)Soil spectral signature 2;(c)Soil spectral signature 3;(d)Water spectral signature;(e)Spectral signature of growth vegetation;(f)Spectral signature of matured vegetation
图2 不同地物类型的光谱特征
(a)土壤光谱特征1;(b)土壤光谱特征2;(c)土壤光谱特征3;(d)水体光谱特征;(e)生长期植被光谱特征;(f)成熟期植被光谱特征
Fig.2 Spectral signature of different ground types
(a)Soil spectral signature 1;(b)Soil spectral signature 2;(c)Soil spectral signature 3;(d)Water spectral signature;(e)Spectral signature of growth vegetation;(f)Spectral signature of matured vegetation
表2 GF-1 WFV2缨帽变换系数
Table 2 TCT coefficient of GF-1 WFV2
缨帽变换分量 蓝光波段 绿光波段 红光波段 近红外波段
1. 亮度分量 0.2428 0.4347 0.4169 0.7604
2. 绿度分量 -0.2637 -0.4480 -0.5584 0.6465
3. 湿度分量 -0.0392 0.7530 -0.6542 -0.0593
4. 第4分量 -0.9327 0.2082 0.2939 0.0177
表2 GF-1 WFV2缨帽变换系数
Table 2 TCT coefficient of GF-1 WFV2
缨帽变换分量 蓝光波段 绿光波段 红光波段 近红外波段
1. 亮度分量 0.2428 0.4347 0.4169 0.7604
2. 绿度分量 -0.2637 -0.4480 -0.5584 0.6465
3. 湿度分量 -0.0392 0.7530 -0.6542 -0.0593
4. 第4分量 -0.9327 0.2082 0.2939 0.0177
图3 原始影像及缨帽变换分量图(实验影像:驻马店;验证影像:营口市、连云港、贺州市)
Fig.3 The original images and their TCT components (Test image: Zhumadian; Validation image: Yingkou, Lianyungang and Hezhou)
图3 原始影像及缨帽变换分量图(实验影像:驻马店;验证影像:营口市、连云港、贺州市)
Fig.3 The original images and their TCT components (Test image: Zhumadian; Validation image: Yingkou, Lianyungang and Hezhou)
图4 地物在不同投影平面上的理论分布特征
(a)植被平面;(b)土壤平面;(c)过渡带平面
Fig.4 Distribution characteristics of ground objects in different projection planes
(a)Plane of vegetation view;(b)Plane of soils view;(c)Plane of transition zone view
图4 地物在不同投影平面上的理论分布特征
(a)植被平面;(b)土壤平面;(c)过渡带平面
Fig.4 Distribution characteristics of ground objects in different projection planes
(a)Plane of vegetation view;(b)Plane of soils view;(c)Plane of transition zone view
表3 GF-1和Landsat 8同步影像验证
Table 3 Validation using Synchronous images between GF-1 and Landsat 8
同步影像 缨帽变换分量 相关系数 均方根误差
亮度 0.8981 0.1079
聊城市 绿度 0.8914 0.0956
湿度 0.8282 0.052
亮度 0.9588 0.098
连云港 绿度 0.9181 0.0671
湿度 0.8758 0.0598
表3 GF-1和Landsat 8同步影像验证
Table 3 Validation using Synchronous images between GF-1 and Landsat 8
同步影像 缨帽变换分量 相关系数 均方根误差
亮度 0.8981 0.1079
聊城市 绿度 0.8914 0.0956
湿度 0.8282 0.052
亮度 0.9588 0.098
连云港 绿度 0.9181 0.0671
湿度 0.8758 0.0598
表4 前2、3个分量能解释的信息变化量
Table 4 Total variance explained by the first 2 or 3 component
卫星传感器 前2个
分量/%		 前3个
分量/%		 第3
分量/%		 参考文献
IKONOS 98.3 99.8 1.5 [20]
QuickBird 2 98.3 99.7 1.4 [21]
CBERS-02B / 98.0 / [25]
HJ-1 A/B 98.0 99.9 1.9 [26]
GF-1 WFV2 97.3 98.9 1.6 本文
表4 前2、3个分量能解释的信息变化量
Table 4 Total variance explained by the first 2 or 3 component
卫星传感器 前2个
分量/%		 前3个
分量/%		 第3
分量/%		 参考文献
IKONOS 98.3 99.8 1.5 [20]
QuickBird 2 98.3 99.7 1.4 [21]
CBERS-02B / 98.0 / [25]
HJ-1 A/B 98.0 99.9 1.9 [26]
GF-1 WFV2 97.3 98.9 1.6 本文
图5 缨帽变换结果在不同投影平面的分布特征 (0~255拉伸)
Fig.5 The tasseled cap distribution in different projection plane (Stretched between 0~255)
图5 缨帽变换结果在不同投影平面的分布特征 (0~255拉伸)
Fig.5 The tasseled cap distribution in different projection plane (Stretched between 0~255)
图6 GF-1其他3个WFV传感器缨帽变换结果
Fig.6 Tasseled cap transformation results of the other three WFV sensors in GF-1
图6 GF-1其他3个WFV传感器缨帽变换结果
Fig.6 Tasseled cap transformation results of the other three WFV sensors in GF-1
表5 与WFV2传感器相近时相的影像对
Table 5 Images pairs between WFV2 and other sensors
地点 WFV1 WFV2 WFV3 WFV4 相差天数/d
淮安 20160503 20160429 / / 4
扬州 20160503 20160429 / / 4
淮安 / 20160429 20160425 / 4
营口 / 20160919 20160920 / 1
淮安 / 20160429 / 20160430 1
聊城 / 20160504 / 20160517 13
表5 与WFV2传感器相近时相的影像对
Table 5 Images pairs between WFV2 and other sensors
地点 WFV1 WFV2 WFV3 WFV4 相差天数/d
淮安 20160503 20160429 / / 4
扬州 20160503 20160429 / / 4
淮安 / 20160429 20160425 / 4
营口 / 20160919 20160920 / 1
淮安 / 20160429 / 20160430 1
聊城 / 20160504 / 20160517 13
表6 不同传感器之间缨帽变换的相关性统计
Table 6 Correlation statistics of tasseled cap transformation from different sensors
对比传感器 影像地区 亮度 绿度 湿度 第4分量
淮安 0.861 0.930 0.886 0.704
WFV1-WFV2 扬州 0.936 0.941 0.865 0.614
平均 0.899 0.936 0.876 0.659
淮安 0.970 0.985 0.955 0.934
WFV3-WFV2 营口 0.908 0.887 0.893 0.791
平均 0.939 0.936 0.924 0.863
淮安 0.810 0.909 0.853 0.706
WFV4-WFV2 聊城 0.751 0.907 0.854 0.665
平均 0.781 0.908 0.854 0.686
表6 不同传感器之间缨帽变换的相关性统计
Table 6 Correlation statistics of tasseled cap transformation from different sensors
对比传感器 影像地区 亮度 绿度 湿度 第4分量
淮安 0.861 0.930 0.886 0.704
WFV1-WFV2 扬州 0.936 0.941 0.865 0.614
平均 0.899 0.936 0.876 0.659
淮安 0.970 0.985 0.955 0.934
WFV3-WFV2 营口 0.908 0.887 0.893 0.791
平均 0.939 0.936 0.924 0.863
淮安 0.810 0.909 0.853 0.706
WFV4-WFV2 聊城 0.751 0.907 0.854 0.665
平均 0.781 0.908 0.854 0.686
图7 GF-1 WFV传感器观测方向
Fig.7 Observation direction of different WFV sensors in GF-1
图7 GF-1 WFV传感器观测方向
Fig.7 Observation direction of different WFV sensors in GF-1
[1] Kauth R J, Thomas G S.The tasselled cap—A graphic description of the spectral-temporal development of agricultural crops as seen by Landsat[C]//LARS Symposia, 1976: 159.
[2] Zanchetta A, Bitelli G, Karnieli A.Monitoring desertification by remote sensing using the Tasselled Cap transform for long-term change detection[J]. Natural Hazards, 2016, 83(1):1-15.
[3] Zanchetta A, Bitelli G, Karnieli A.Tasselled Cap transform for change detection in the drylands: Findings for SPOT and Landsat satellites with FOSS tools[C]//International Conference on Remote Sensing and Geoinformation of the Environment. International Society for Optics and Photonics, 2015:1 287-1 291.
[4] Santra A, Mitra S S.A comparative study of tasselled cap transformation of DMC and ETM+ images and their application in forest classification[J]. Journal of the Indian Society of Remote Sensing, 2014, 42(2):373-381.
doi: 10.1007/s12524-013-0313-0     URL    
[5] Xu Hanqiu.A remote sensing index for assessment of regional ecological changes[J]. China Environmental Science, 2013, 33(5): 889-897.
[徐涵秋. 区域生态环境变化的遥感评价指数[J]. 中国环境科学, 2013, 33(5):889-897.]
doi: 10.3969/j.issn.1000-6923.2013.05.019     URL    
[6] Jong S M D. Derivation of vegetative variables from a landsat TM image for modelling soil erosion[J]. Earth Surface Processes & Landforms, 2010, 19(2):165-178.
doi: 10.1002/esp.3290190207     URL    
[7] Han T, Wulder M A, White J C, et al. An efficient protocol to process landsat images for change detection with tasselled cap transformation[J]. IEEE Geoscience & Remote Sensing Letters, 2007, 4(1):147-151.
[8] Lobser S E, Cohen W B.MODIS tasselled cap: Land cover characteristics expressed through transformed MODIS data[J]. International Journal of Remote Sensing, 2007, 28(22): 5 079-5 101.
doi: 10.1080/01431160701253303     URL    
[9] Song C, Ren H, Huang C.Estimating soil salinity in the Yellow River Delta, Eastern China—An integrated approach using spectral and terrain indices with the generalized additive model[J]. Pedosphere, 2016, 26(5):626-635.
doi: 10.1016/S1002-0160(15)60071-6     URL    
[10] Karl J W, Maurer B A.Multivariate correlations between imagery and field measurements across scales: Comparing pixel aggregation and image segmentation[J]. Landscape Ecology, 2010, 25(4): 591-605.
doi: 10.1007/s10980-009-9439-4     URL    
[11] Chen Chao, Qin Qiming, Wang Jinliang, et al. Comparison of quality evaluation methods for image fusion of farmland remote sensing[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(10):95-100.
[陈超, 秦其明, 王金梁,等. 农地遥感图像融合质量评价方法比较[J]. 农业工程学报, 2011, 27(10):95-100.]
doi: 10.3969/j.issn.1002-6819.2011.10.017     URL    
[12] Crist E P, Cicone R C.A physically-based transformation of thematic mapper data—The TM tasseled cap[J]. IEEE Transactions on Geoscience & Remote Sensing, 1984, 22(3):256-263.
[13] Crist E P, Cicone R C.Application of the Tasseled Cap concept to simulated thematic mapper data[J]. Photogrammetric Engineering & Remote Sensing, 1984, 50(3):343-352.
doi: 10.1016/0031-8663(84)90021-8     URL    
[14] Huang C, Wylie B, Yang L, et al. Derivation of a tasselled cap transformation based on Landsat 7 at-satellite reflectance[J]. International Journal of Remote Sensing, 2002, 23(8): 1 741-1 748.
doi: 10.1080/01431160110106113     URL    
[15] Baig M H A, Zhang L, Tong S,et al.Derivation of a tasselled cap transformation based on Landsat 8 at-satellite reflectance[J]. Remote Sensing Letters, 2014, 5(5):423-431.
doi: 10.1080/2150704X.2014.915434     URL    
[16] Liu Q, Liu G, Huang C, et al. A tasseled cap transformation for Landsat 8 OLI TOA reflectance images[C]//Geoscience and Remote Sensing Symposium. IEEE, 2014:541-544.
[17] Liu Q, Liu G, Huang C, et al. Comparison of tasselled cap transformations based on the selective bands of Landsat 8 OLI TOA reflectance images[J]. International Journal of Remote Sensing, 2015, 36(2):417-441.
doi: 10.1080/01431161.2014.995274     URL    
[18] Li Bolun, Wei Chaopu, Yan Xiaoyuan.Study of derivation of tasseled cap transformation for Landsat 8 OLI images[J]. Science of Surveying and Mapping, 2016, 41(4):102-107.
[李博伦, 违超普, 颜晓元. Landsat8陆地成像仪影像的缨帽变换推导[J]. 测绘科学, 2016, 41(4):102-107.]
[19] Yarbrough L D, Easson G, Kuszmaul J S.Using at-sensor radiance and reflectance tasseled cap transforms applied to change detection for the ASTER sensor[J]. Aster, 2005, 2: 141-145.
[20] Horne J H.A tasseled cap transformation for IKONOS images[C]//ASPRS 2003 Annual Conference Proceedings. 2003.
[21] Yarbrough L D, Easson G, Kuszmaul J S.QuickBird 2 tasseled cap transform coefficients: A comparison of derivation methods[C]//Pecora 16 Global Priorites in Land Remote Sensing, Sioux Falls, South Pakota, 2005, 16: 23-27.
[22] Verdin J, Eckhardt D, Lyford G.Evaluation of SPOT imagery for monitoring irrigated lands[C]//Proceedings International Conference on SPOT 1 Image Utilization, Assessment, Results p 81-91(SEE N 88-28346 22-43), Paris,1987.
[23] Silva D A.Determination of 'tasseled cap' transformation parameters for images obtained by the SPOT satellite[C]//International Symposium on Remote Sensing of Environment, 24th, Rio de Janeiro, Brazil, 1992: 291-300.
[24] Ivits E, Lamb A, Langar F,et al. Orthogonal transformation of segmented SPOT5 images[J]. Photogrammetric Engineering & Remote Sensing, 2008, 74(11): 1 351-1 364.
doi: 10.2352/J.ImagingSci.Technol.(2008)52:6(060508)     URL    
[25] Sheng L, Huang J, Tang X.A tasseled cap transformation for CBERS-02B CCD data[J]. Journal of Zhejiang University (Science B), 2011, 12(9): 780-786.
[26] Chen C, Tang P, Bian Z. Tasseled cap transformation for HJ-1A/B charge coupled device images[J]. Journal of Applied Remote Sensing, 2012, 6(1): 063575-1-063575-11.
doi: 10.1117/1.JRS.6.063575     URL    
[27] Chander G, Markham B L, Helder D L.Summary of current radiometric calibration coefficients for Landsat MSS, TM, ETM+, and EO-1 ALI sensors[J]. Remote Sensing of Environment, 2009, 113(5): 893-903.
doi: 10.1016/j.rse.2009.01.007     URL    
[28] Hanqiu X.Image-based normalization technique used for Landsat TM/ETM+ imagery[J]. Geomatics and Information Science of Wuhan University, 2007, 32(1): 62-66.
doi: 10.3321/j.issn:1671-8860.2007.01.016     URL    
[29] Xu H.Retrieval of reflectance and land surface temperature of the newly-launched Landsat 8 satellite[J]. Chinese Journal of Geophysics, 2015, 58(3):741-747.
[30] Jackson R D.Spectral indices in n-space[J]. Remote Sensing of Environment, 1983, 13(5): 409-421.
doi: 10.1016/0034-4257(83)90010-X     URL    
[1] 梁承弘, 鹿化煜. 风成沉积物叶蜡氢同位素在揭示东亚季风区干湿变化中的原理及应用[J]. 地球科学进展, 2021, 36(1): 45-57.
[2] 念达, 邓琪敏, 付遵涛. 相对湿度及其变化的年循环研究进展[J]. 地球科学进展, 2018, 33(7): 762-774.
[3] 孙丽, 赵姝慧. 探空仪湿度测量误差研究现状及其对云识别的影响[J]. 地球科学进展, 2018, 33(1): 85-92.
[4] 李得勤, 张述文, 文小航, 贺慧. 土壤湿度参数化及对天气和气候模拟影响的研究进展[J]. 地球科学进展, 2016, 31(3): 236-247.
[5] 兰鑫宇, 郭子祺, 田野, 雷霞, 王婕. 土壤湿度遥感估算同化研究综述[J]. 地球科学进展, 2015, 30(6): 668-679.
[6] 毛伏平, 张述文, 叶丹, 杨茜茜. 模式时间关联误差对集合平方根滤波估算土壤湿度的影响[J]. 地球科学进展, 2015, 30(6): 700-708.
[7] 李得勤，段云霞，张述文. 土壤湿度观测、模拟和估算研究[J]. 地球科学进展, 2012, 27(4): 424-434.
[8] 顾凯，施斌，唐朝生，高磊，姜洪涛. 城区与郊区土体湿度监测与对比分析——以南京市为例[J]. 地球科学进展, 2012, 27(2): 185-193.
[9] 苟浩锋,刘彦华,张述文,李得勤. 评估集合卡曼滤波反演土壤湿度廓线的性能[J]. 地球科学进展, 2010, 25(4): 400-407.
[10] 王晓婷,郭维栋,钟中,崔晓燕. 中国东部土壤温度、湿度变化的长期趋势及其与气候背景的联系[J]. 地球科学进展, 2009, 24(2): 181-191.
[11] 关止,赵凯,宋冬生. 利用反射GPS信号遥感土壤湿度[J]. 地球科学进展, 2006, 21(7): 747-750.
[12] 张生雷,谢正辉,田向军,师春香,陈锋. 基于土壤水模型及站点资料的土壤湿度同化方法[J]. 地球科学进展, 2006, 21(12): 1350-1362.
[13] 马柱国，符淙斌，谢力，陈文海，陶树望. 土壤湿度和气候变化关系研究中的某些问题[J]. 地球科学进展, 2001, 16(4): 563-566.
[14] 马柱国,魏和林,符淙斌. 土壤湿度与气候变化关系的研究进展与展望[J]. 地球科学进展, 1999, 14(3): 299-305.
阅读次数
全文


摘要

地址：甘肃省兰州市天水中路8号
QQ：114723517
电话：0931-8762293　E-mail：adearth@lzb.ac.cn
陇ICP备2021001824号-5　 甘公网安备62010202000687