留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于MODIS数据的渤海海冰厚度反演算法优化

朱星源 苏洁 宋梅 杨茜 梁韵

朱星源,苏洁,宋梅,等. 基于MODIS数据的渤海海冰厚度反演算法优化[J]. 海洋学报,2022,44(12):70–83 doi: 10.12284/hyxb2022141
引用本文: 朱星源,苏洁,宋梅,等. 基于MODIS数据的渤海海冰厚度反演算法优化[J]. 海洋学报,2022,44(12):70–83 doi: 10.12284/hyxb2022141
Zhu Xingyuan,Su Jie,Song Mei, et al. Optimization of the Bohai Sea ice thickness retrieval algorithm based on MODIS data[J]. Haiyang Xuebao,2022, 44(12):70–83 doi: 10.12284/hyxb2022141
Citation: Zhu Xingyuan,Su Jie,Song Mei, et al. Optimization of the Bohai Sea ice thickness retrieval algorithm based on MODIS data[J]. Haiyang Xuebao,2022, 44(12):70–83 doi: 10.12284/hyxb2022141

基于MODIS数据的渤海海冰厚度反演算法优化

doi: 10.12284/hyxb2022141
基金项目: 国家自然科学基金面上项目 (42176208)。
详细信息
    作者简介:

    朱星源(1997-),男,山东省临沂市人,主要研究方向为卫星海洋学。E-mail:xingxingzhihuo123@foxmail.com

    通讯作者:

    苏洁(1966-),女,教授, 主要从事北极海冰热力学研究。E-mail:sujie@ouc.edu.cn

  • 中图分类号: P731.15

Optimization of the Bohai Sea ice thickness retrieval algorithm based on MODIS data

  • 摘要: 海冰厚度是监测与研究渤海海冰的重要参数。为了获取更加可靠的渤海海冰厚度数据,本研究基于MODIS数据对海冰厚度反演中的冰水分离环节和冰厚计算方法都进行了改进。对于冰水分离环节,本文在Canny边缘检测算子提取海冰基础上,加入了二值化处理、阈值判别等步骤,实现了较高精度的渤海海冰范围自动化提取。通过试验确定了海冰厚度与反照率指数关系模型中的参数,包括海冰衰减系数和海水反照率参数,使其更加符合渤海海区的物理特征。将改进后算法的海冰厚度反演结果与渤海海上石油平台实测数据进行比较,并分析了误差来源。结果表明,经过对算法的改进,海冰厚度与反照率指数关系模型的反演结果与实测数据之间的平均绝对误差由7.05 cm缩小到2.74 cm,相关系数由0.434提高到0.485。
  • 图  1  获取实测数据的石油平台位置

    Fig.  1  Positions of oil platform observing the measured data

    图  2  渤海海冰厚度反演算法流程

    图中红色字部分为本文针对Canny算子提取海冰裂缝[32]新加入的步骤,蓝色字部分为针对海冰厚度指数关系模型[17]参数设置的改进

    Fig.  2  Flow chart of Bohai Sea ice thickness retrieval algorithm

    The red character part in the figure is the new step of the Canny operator to extract sea ice cracks[32], and the blue character part is the improvement of the parameter setting of the sea ice thickness exponential model[17]

    图  3  2010年1月22日渤海MODIS真彩图

    Fig.  3  MODIS true color image of the Bohai Sea on January 22, 2010

    图  4  2010年1月22日渤海冰水分离过程

    Fig.  4  Process diagram of ice-water separation of the Bohai Sea on January 22, 2010

    图  5  频数分布(a)和粗糙度分布(b)

    Fig.  5  Frequency distribution (a) and roughness index distribution (b)

    图  6  渤海海区反照率(a)和海水反照率插值(b)

    Fig.  6  The albedo of Bohai Sea (a) and the albedo of sea water interpolation (b)

    图  7  训练数据集实测冰厚与衰减系数散点图

    绿色虚线为帮助描述收敛过程的辅助线

    Fig.  7  Scatter plot between attenuation coefficient and measured sea ice thickness of training data sets

    The green dotted line is the auxiliary line to help describe the convergence process

    图  8  模型中各参数对反演结果的影响

    Fig.  8  The influence of each parameter in the model on the retrieval result

    图  9  渤海海冰厚度反演结果对比

    Fig.  9  Comparison of retrieval results of Bohai Sea sea ice thickness

    图  10  2020–2021年冬季反演结果与实测数据的时间序列比较(JZ20-2平台和JZ9-3平台)

    日期后用*标注的来自测试数据集,其他为训练数据

    Fig.  10  Time series comparison between measured data and retrieval results in 2020–2021 winter (JZ9-3 and JZ20-2)

    The date with * means the data come from the test dataset, the others come from the trainning dataset

    图  11  2021年1月5日(a, e)、7日(b, f)、8日(c, g)、13日(d, h)辽东湾MODIS真彩图及海冰厚度反演结果

    Fig.  11  MODIS true color image and sea ice thickness retrieval results of Liaodong Bay on January 5 (a, e), 7 (b, f), 8 (c, g), 13 (d, h), 2021

    表  1  MODIS部分波段光谱范围

    Tab.  1  Spectral range of MODIS partial bands

    波段光谱范围/nm分辨率/m
    1620~670250(重采样为1 000)
    2841~876250(重采样为1 000)
    3459~479500(重采样为1 000)
    4545~565500(重采样为1 000)
    51 230~1 250500(重采样为1 000)
    61 628~1 652500(重采样为1 000)
    72 105~2 135500(重采样为1 000)
    3111.770~12.2701 000
    3213.185~13.4851 000
    下载: 导出CSV

    表  2  海冰厚度测试数据集实测数据与反演结果

    Tab.  2  Measured data of test data sets and retrieval results of Bohai Sea sea ice thickness

    日期站位最大实测
    厚度/cm
    平均实测
    厚度/cm
    改进前算法厚
    度(T0)/cm
    改进后算法厚
    度(T1)/cm
    仅改变海水反照
    率厚度(T2)/cm
    仅改变衰减系
    数厚度(T3)/cm
    Yuan等[26]
    法厚度/cm
    2009年12月19日*JZ9-353.511.425.68.067.938.95
    2010年1月6日*JZ9-38616.69.4913.6611.5412.7
    2010年1月11日*JZ20-25314.245.888.479.99.28
    2010年1月24日*JZ9-31512.513.66.399.29.457.99
    2010年2月5日*JZ9-3201424.7614.1620.3717.215.65
    2010年2月16日*JZ20-28611.385.598.057.911.61
    2010年2月16日*JZ9-3644.040.520.742.812.49
    2010年1月29日*JZ20-2435.922.663.834.113.37
    2010年2月12日*JZ9-38610.95.778.37.579.85
    2013年1月26日*JZ20-285.512.236.749.718.511.21
    2013年2月2日*JX1-1107.510.656.859.857.410.87
    2013年2月2日*JZ25-1S3022.517.0911.2916.2511.8714.81
    2013年12月28日*JZ20-2421.710.821.181.191.45
    2014年1月12日JZ20-243.517.8410.7815.5212.413.8
    2014年1月18日JZ9-3104.510.945.578.027.610
    2014年2月11日JZ20-28412.376.849.848.5911.61
    2016年1月7日JZ9-31269.223.975.726.416.85
    2016年1月8日JZ9-313810.444.56.487.257.25
    2016年1月9日JZ20-2197.210.146.088.757.0510.45
    2016年1月12日JZ20-2135.810.295.217.57.158.64
    2016年1月13日JZ20-2135.512.464.66.628.668.01
    2016年1月18日JZ20-210520.5612.0917.4114.2814.42
    2020年12月31日JZ9-314712.687.9911.498.8113.41
    2021年1月1日JZ9-3126.517.298.9512.8912.0114.83
    2021年1月7日JZ20-214714.019.113.19.7414.35
    2021年1月8日JZ9-314822.0213.6219.615.317.68
    2021年1月10日JZ20-212619.599.8214.1313.6114.65
    2021年1月13日JZ9-310718.6310.6915.3912.9515.35
    2021年1月17日JZ9-314915.597.9911.510.8313.78
    平均实测厚度误差平均误差6.660.493.662.574.13
    平均绝对误差7.052.744.723.725.17
    均方根误差8.253.755.84.735.94
    相关系数0.4340.4850.4850.4340.417
    最大实测厚度
    误差
    平均误差2.37−3.810.631.720.16
    平均绝对误差4.734.674.043.873.96
    均方根误差6.026.095.235.254.99
    相关系数0.4350.4800.4800.4350.459
    注:T0为算法改进前的冰厚反演结果;T1为算法改进后的冰厚反演结果;T2为仅改进海水反照率的冰厚反演结果;T3为仅改进衰减系数的冰厚反演结果。日期后标注*的数据来源于参考文献[24, 35]。
    下载: 导出CSV
  • [1] 白珊, 刘钦政, 吴辉碇, 等. 渤海、北黄海海冰与气候变化的关系[J]. 海洋学报, 2001, 23(5): 33−41.

    Bai Shan, Liu Qinzheng, Wu Huiding, et al. Relation of ice conditions with climate change in the Bohai Sea and the northern Huanghai Sea[J]. Haiyang Xuebao, 2001, 23(5): 33−41.
    [2] Gong Daoyi, Wang Shaowu, Zhu Jinhong. East Asian winter monsoon and Arctic Oscillation[J]. Geophysical Research Letters, 2001, 28(10): 2073−2076.
    [3] 唐茂宁, 洪洁莉, 刘煜, 等. 气候因子对渤海冰情影响的统计分析[J]. 海洋通报, 2015, 34(2): 152−157. doi: 10.11840/j.issn.1001-6392.2015.02.005

    Tang Maoning, Hong Jieli, Liu Yu, et al. Statistical analysis of climatic factors impacting on the Bohai Sea ice[J]. Marine Science Bulletin, 2015, 34(2): 152−157. doi: 10.11840/j.issn.1001-6392.2015.02.005
    [4] 药蕾, 苏洁. 渤海海冰与西伯利亚高压之间的关系及与北大西洋涛动之间的可能联系[J]. 中国海洋大学学报(自然科学版), 2018, 48(6): 1−12. doi: 10.16441/j.cnki.hdxb.20170128

    Yao Lei, Su Jie. Relationships between Bohai Sea ice and Siberian high and possible connections between Bohai Sea ice and north Atlantic Oscillation[J]. Periodical of Ocean University of China, 2018, 48(6): 1−12. doi: 10.16441/j.cnki.hdxb.20170128
    [5] 耿淑琴, 王旭. 2001−2002年冬季黄渤海天气气候特征及对渤海海冰的影响[J]. 海洋预报, 2002, 19(4): 38−47. doi: 10.3969/j.issn.1003-0239.2002.04.006

    Geng Shuqin, Wang Xu. Features of weather and climate in winter of 2001−2002 and its effect on sea ice in Bohai[J]. Marine Forecasts, 2002, 19(4): 38−47. doi: 10.3969/j.issn.1003-0239.2002.04.006
    [6] 顾卫, 史培军, 刘杨, 等. 渤海和黄海北部地区负积温资源的时空分布特征[J]. 自然资源学报, 2002, 17(2): 168−173. doi: 10.3321/j.issn:1000-3037.2002.02.007

    Gu Wei, Shi Peijun, Liu Yang, et al. The characteristics of temporal and spatial distribution of negative accumulated temperature in Bohai Sea and North Yellow Sea[J]. Journal of Natural Resources, 2002, 17(2): 168−173. doi: 10.3321/j.issn:1000-3037.2002.02.007
    [7] 李彦青, 苏洁, 汪洋, 等. 渤海海冰外缘线候平均离岸距离的变化及其关键影响因子[J]. 中国海洋大学学报(自然科学版), 2013, 43(7): 7−16. doi: 10.16441/j.cnki.hdxb.2013.07.002

    Li Yanqing, Su Jie, Wang Yang, et al. Variability of the pentadly average distance between the sea ice edge and the coast in the Bohai Sea and its key impact factors[J]. Periodical of Ocean University of China, 2013, 43(7): 7−16. doi: 10.16441/j.cnki.hdxb.2013.07.002
    [8] 刘煜, 吴辉碇. 第1讲: 渤、黄海的海冰[J]. 海洋预报, 2017, 34(3): 94−101. doi: 10.11737/j.issn.1003-0239.2017.03.012

    Liu Yu, Wu Huiding. Sea ice in the Bohai Sea and the northern Yellow Sea[J]. Marine Forecasts, 2017, 34(3): 94−101. doi: 10.11737/j.issn.1003-0239.2017.03.012
    [9] 邓树奇. 渤海海冰灾害及其预防概况[J]. 灾害学, 1986(1): 80.

    Den Shuqi. Sea ice disaster and its prevention in Bohai Sea[J]. Journal of Catastrophology, 1986(1): 80.
    [10] 国家海洋局. 2010年中国海洋灾害公报[EB/OL]. (2011−04−26) [2022−05−06]. http://www.nmdis.org.cn/hygb/zghyzhgb/2010nzghyzhgb/.

    State Oceanic Administration. Bulletin of China Marine Disaster 2010[EB/OL]. (2011−04−26) [2022−05−06]. http://www.nmdis.org.cn/hygb/zghyzhgb/2010nzghyzhgb/.
    [11] 国家海洋局. 2013年中国海洋灾害公报[EB/OL]. (2014−03−24) [2022−05−06]. http://www.nmdis.org.cn/hygb/zghyzhgb/2013nzghyzhgb/.

    State Oceanic Administration. Bulletin of China Marine Disaster 2013[EB/OL]. (2014−03−24) [2022−05−06]. http://www.nmdis.org.cn/hygb/zghyzhgb/2013nzghyzhgb/.
    [12] 郑新江, 邱康睦, 陆风. 定量计算渤海海冰参数的遥感方法[J]. 应用气象学报, 1998, 9(3): 312−316.

    Zheng Xinjiang, Qiu Kangmu, Lu Feng. Quantitative calculation of sea ice over the Bohai Sea using NOAA/AVHRR imagery[J]. Quarterly Journal of Applied Meteorology, 1998, 9(3): 312−316.
    [13] 罗亚威, 张蕴斐, 孙从容, 等. “海洋1号”卫星在海冰监测和预报中的应用[J]. 海洋学报, 2005, 27(1): 7−18.

    Luo Yawei, Zhang Yunfei, Sun Congrong, et al. Application of the “HY-1” satellite to sea ice monitoring and forecasting[J]. Haiyang Xuebao, 2005, 27(1): 7−18.
    [14] 吴龙涛, 吴辉碇, 孙兰涛, 等. MODIS渤海海冰遥感资料反演[J]. 中国海洋大学学报 (自然科学版), 2006, 36(2): 173−179.

    Wu Longtao, Wu Huiding, Sun Lantao, et al. Retrieval of sea ice in the Bohai Sea from MODIS data[J]. Periodical of Ocean University of China, 2006, 36(2): 173−179.
    [15] 刘志强, 苏洁, 时晓旭, 等. 渤海AVHRR多通道海冰密集度反演算法试验研究[J]. 海洋学报, 2014, 36(11): 74−84.

    Liu Zhiqiang, Su Jie, Shi Xiaoxu, et al. Study on the multi-band retrieval algorithm for the Bohai Sea ice concentration using AVHRR data[J]. Haiyang Xuebao, 2014, 36(11): 74−84.
    [16] Zhou Ye, Fan Zhaopeng. Study on Bohai sea ice based on MODIS data[J]. E3S Web of Conferences, 2018, 53: 03038.
    [17] Grenfell T C. A radiative transfer model for sea ice with vertical structure variations[J]. Journal of Geophysical Research, 1991, 96(C9): 16991−17001.
    [18] 谢锋, 顾卫, 袁艺, 等. 辽东湾海冰资源量的遥感估算方法研究[J]. 资源科学, 2003, 25(3): 17−23. doi: 10.3321/j.issn:1007-7588.2003.03.003

    Xie Feng, Gu Wei, Yuan Yi, et al. Estimation of sea ice resources in Liaodong Gulf using remote sensing[J]. Resources Science, 2003, 25(3): 17−23. doi: 10.3321/j.issn:1007-7588.2003.03.003
    [19] Yuan Shuai, Gu Wei, Xu Yingjun, et al. The estimate of sea ice resources quantity in the Bohai Sea based on NOAA/AVHRR data[J]. Acta Oceanologica Sinica, 2012, 31(1): 33−40. doi: 10.1007/s13131-012-0173-4
    [20] Su Hua, Wang Yunpeng. Using MODIS data to estimate sea ice thickness in the Bohai Sea (China) in the 2009−2010 winter[J]. Journal of Geophysical Research: Oceans, 2012, 117(C10): C10018.
    [21] Liu Wensong, Sheng Hui, Zhang Xi. Sea ice thickness estimation in the Bohai Sea using geostationary ocean color imager data[J]. Acta Oceanologica Sinica, 2016, 35(7): 105−112.
    [22] Xu Zhantang, Yang Yuezhong, Wang Guifen, et al. Optical properties of sea ice in Liaodong Bay, China[J]. Journal of Geophysical Research: Oceans, 2012, 117(C3): C03007.
    [23] 陈树果. 黄渤海水体光学性质变化及其影响机制[D]. 青岛: 中国海洋大学, 2015.

    Chen Shuguo. Variations and influencing mechanisms in the optical properties of the waters in the Yellow Sea and Bohai Sea[D]. Qingdao: Ocean University of China, 2015.
    [24] Zeng Tao, Shi Lijian, Marko M, et al. Sea ice thickness analyses for the Bohai Sea using MODIS thermal infrared imagery[J]. Acta Oceanologica Sinica, 2016, 35(7): 96−104.
    [25] Ning Li, Xie Feng, Gu Wei, et al. Using remote sensing to estimate sea ice thickness in the Bohai Sea, China based on ice type[J]. International Journal of Remote Sensing, 2009, 30(17): 4539−4552.
    [26] Yuan Shuai, Gu Wei, Liu Chengyu, et al. Towards a semi-empirical model of the sea ice thickness based on hyperspectral remote sensing in the Bohai Sea[J]. Acta Oceanologica Sinica, 2017, 36(1): 80−89. doi: 10.1007/s13131-017-0996-0
    [27] Yuan Shuai, Liu Chengyu, Liu Xueqin. Practical model of sea ice thickness of Bohai Sea based on MODIS data[J]. Chinese Geographical Science, 2018, 28(5): 863−872.
    [28] 吴奎桥, 徐莹, 郝轶萌. MODIS数据在海冰遥感中的应用[J]. 海洋预报, 2005, 22(S1): 44−49. doi: 10.3969/j.issn.1003-0239.2005.z1.007

    Wu Kuiqiao, Xu Ying, Hao Yimeng. Application in sea ice remote sensing of MODIS data[J]. Marine Forecasts, 2005, 22(S1): 44−49. doi: 10.3969/j.issn.1003-0239.2005.z1.007
    [29] Su Hua, Wang Yunpeng, Xiao Jie, et al. Improving MODIS sea ice detectability using gray level co-occurrence matrix texture analysis method: a case study in the Bohai Sea[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2013, 85: 13−20.
    [30] Zhang Na, Wu Yongshen, Zhang Qinghe. Detection of sea ice in sediment laden water using MODIS in the Bohai Sea: a CART decision tree method[J]. International Journal of Remote Sensing, 2015, 36(6): 1661−1674. doi: 10.1080/01431161.2015.1015658
    [31] Su Hua, Wang Yunpeng, Xiao Jie, et al. Classification of MODIS images combining surface temperature and texture features using the support vector machine method for estimation of the extent of sea ice in the frozen Bohai Bay, China[J]. International Journal of Remote Sensing, 2015, 36(10): 2734−2750. doi: 10.1080/01431161.2015.1041619
    [32] Li Yawen, Yang Daiqin. Extraction of Bohai Sea ice from MODIS data based on multi-constraint endmembers and linear spectral unmixing[J]. International Journal of Remote Sensing, 2020, 41(14): 5525−5548.
    [33] Canny J. A computational approach to edge detection[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 1986, PAMI-8(6): 679−698. doi: 10.1109/TPAMI.1986.4767851
    [34] 刘荣高, 刘洋, 刘纪远. MODIS科学数据处理研究进展[J]. 自然科学进展, 2009, 19(2): 141−147. doi: 10.3321/j.issn:1002-008X.2009.02.003

    Liu Ronggao, Liu Yang, Liu Jiyuan. Research progress of MODIS scientific data processing[J]. Progress in Natural Science, 2009, 19(2): 141−147. doi: 10.3321/j.issn:1002-008X.2009.02.003
    [35] Karvonen J, Shi Lijian, Cheng Bin, et al. Bohai Sea ice parameter estimation based on thermodynamic ice model and earth observation data[J]. Remote Sensing, 2017, 9(3): 234.
    [36] 李微, 方圣辉, 佃袁勇, 等. 基于光谱分析的MODIS云检测算法研究[J]. 武汉大学学报(信息科学版), 2005, 30(5): 435−438, 443.

    Li Wei, Fang Shenghui, Dian Yuanyong, et al. Cloud detection in MODIS data based on spectrum analysis[J]. Geomatics and Information Science of Wuhan University, 2005, 30(5): 435−438, 443.
    [37] 马开玉, 丁裕国, 屠其璞, 等. 气候统计原理与方法[M]. 北京: 气象出版社, 1993.

    Ma Kaiyu, Ding Yuguo, Tu Qipu, et al. Principles and Methods of Climate Statistics[M]. Beijing: China Meteorological Press, 1993.
    [38] 许向阳, 宋恩民, 金良海. Otsu准则的阈值性质分析[J]. 电子学报, 2009, 37(12): 2716−2719. doi: 10.3321/j.issn:0372-2112.2009.12.020

    Xu Xiangyang, Song Enmin, Jin Lianghai. Characteristic analysis of threshold based on Otsu criterion[J]. Acta Electronica Sinica, 2009, 37(12): 2716−2719. doi: 10.3321/j.issn:0372-2112.2009.12.020
    [39] 郭衍游, 焦明连. 利用MODIS数据反演渤海海冰分布[J]. 淮海工学院学报(自然科学版), 2010, 19(1): 84−87.

    Guo Yanyou, Jiao Minglian. Using MODIS data to retrieve distribution of sea ice in Bohai Sea[J]. Journal of Huaihai Institute of Technology (Natural Sciences Edition), 2010, 19(1): 84−87.
    [40] 刘良明, 周军元. MODIS数据的海洋表面温度反演[J]. 地理空间信息, 2006, 4(2): 7−9. doi: 10.3969/j.issn.1672-4623.2006.02.003

    Liu Liangming, Zhou Junyuan. Using MODIS imagery to map sea surface temperature[J]. Geospatial Information, 2006, 4(2): 7−9. doi: 10.3969/j.issn.1672-4623.2006.02.003
    [41] 赵进平, 任敬萍. 从航空数字影像提取北极海冰形态参数的方法研究[J]. 遥感学报, 2000, 21(4): 271−278. doi: 10.11834/jrs.20000406

    Zhao Jinping, Ren Jingping. Study on the method to analyze parameters of Arctic sea ice from airborne digital imagery[J]. Journal of Remote Sensing, 2000, 21(4): 271−278. doi: 10.11834/jrs.20000406
    [42] Allison I, Brandt R E, Warren S G. East Antarctic sea ice: albedo, thickness distribution, and snow cover[J]. Journal of Geophysical Research: Oceans, 1993, 98(C7): 12417−12429.
    [43] Liang Shunlin. Narrowband to broadband conversions of land surface albedo I: Algorithms[J]. Remote Sensing of Environment, 2001, 76(2): 213−238.
    [44] 王维波, 苏洁. 基于形态学的海冰外缘线自动提取[J]. 遥感学报, 2015, 19(6): 983−997.

    Wang Weibo, Su Jie. Sea ice edge automatic retrieval based on morphology[J]. Journal of Remote Sensing, 2015, 19(6): 983−997.
    [45] Kemp K K. Encyclopedia of Geographic Information Science[M]. Los Angeles: SAGE Publications, 2008: 146−147.
    [46] 丁德文. 工程海冰学概论[M]. 北京: 海洋出版社, 1999.

    Ding Dewen. Introduction to Engineering Sea Ice Science[M]. Beijing: China Ocean Press, 1999.
  • 加载中
图(11) / 表(2)
计量
  • 文章访问数:  541
  • HTML全文浏览量:  209
  • PDF下载量:  93
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-01-07
  • 修回日期:  2022-06-12
  • 网络出版日期:  2022-10-24
  • 刊出日期:  2023-01-17

目录

    /

    返回文章
    返回