本数据集包括祁连山区域1985年至2017年每5年一期的30m耕地和建筑用地分布产品。该产品基于Landsat系列数据生产。生产耕地产品时,利用夏季合成的Landsat系列数据,使用归一化植被指数NDVI,并辅助以作物的物候期、作物种植类别等先验知识设定产品生产规则;生产建筑用地产品时,利用夏季合成的Landsat系列数据,使用归一化植被指数NDVI、NDBI、MNDWI等指数,辅助以DEM、灯光数据等设定产品生产规则。以2015年的产品为例,由Google Earth高清影像和实地调研数据进行精度评价,得出2015年耕地和建筑用地产品的精度分别为89.43%和91.89%。
吴俊君, 仲波, 角坤升
本数据集包括祁连山区域2018年的30m土地覆盖分类产品。该产品首先利用Landsat-8/OLI构造2015年时间序列数据,针对各类地物随时间变化呈现的NDVI时间序列曲线不同,对不同地物特征进行知识归纳,设定提取规则不同地物信息,得到2015年的土地覆盖分类图。分类系统参考了IGBP分类系统和FROM_LC分类系统,共分为耕地、林地、草地、灌丛、湿地、水体、不透水面、裸地、冰川和积雪共10大类。由Google Earth高清影像和实地调研数据进行精度评价,得出2015年土地覆盖分类产品的总体精度高达92.19%。以2015年的土地覆盖分类产品为底图,按各类别的比例选取大量样本,基于Google Earth Engine平台的Landsat系列数据和强大地数据处理能力,利用深度学习的思想,选取随机森林分类器,对波段信息和NDVI、MNDWI、NDBI等指数进行样本训练,生产出2018年的土地覆盖分类产品。对分类产品进行比较,得出基于Google Earth Engine平台生产的土地覆盖分类产品与基于时间序列方法得到的分类产品具有很好的一致性。总之,祁连山核心区的土地覆盖数据集具有较高的总体精度,且基于Google Earth Engine平台样本训练的方法能够在时间和空间上对现有的分类产品进行扩展,能够在长时间序列上反映更多的土地覆盖类型变化信息。
仲波
该数据集包含了2018年8月31日至2018年12月24日青海湖流域地表过程综合观测网高寒草甸草原混合草原超级站气象要素梯度观测系统数据。站点位于青海省天峻县苏里路旁侧,下垫面是高寒草甸和高寒草原的混合。观测点经纬度为:东经 98°35′41.62″E,北纬 37°42′11.47″N,海拔3718m。风速/风向、空气温度、相对湿度传感器分别架设在3m、5m、10m、15m、20m、30m、40m处,共7层,朝向正北;气压计安装在3m处;翻斗式雨量计安装在塔西偏北侧10m平台上;四分量辐射仪安装在6m处,朝向正南;两个红外温度计安装在6m处,朝向正南,探头朝向是垂直向下;土壤热流板(自校正式)(3块)依次埋设在地下6cm处,朝向正南距离塔体2m处;土壤温度探头埋设在地下5cm、10cm、20cm、40cm、80cm、120cm、200cm、300cm和400cm处,在距离气象塔2m的正东方;土壤水分传感器分别埋设在地下5cm、10cm、20cm、40cm、80cm、120cm、200cm、300cm和400cm处,在距离气象塔2m的正东方;光合有效辐射仪安装在6m处,探头垂直向上和向下方向各一个,朝向正南。 观测项目有:风速(WS_3m、WS_5m、WS_10m、WS_15m、WS_20m、WS_30m、WS_40m)(单位:米/秒)、风向(WD_3m、WD_5m、WD_10m、WD_15m、WD_20m、WD_30m、WD_40m)(单位:度)、空气温湿度(Ta_3m、Ta_5m、Ta_10m、Ta_15m、Ta_20m、Ta_30m、Ta_40m和RH_3m、RH_5m、RH_10m、RH_15m、RH_20m、RH_30m、RH_40m)(单位:摄氏度、百分比)、气压(Press)(单位:百帕)、降水量(Rain)(单位:毫米)、四分量辐射(DR、UR、DLR_Cor、ULR_Cor、Rn)(单位:瓦/平方米)、地表辐射温度(IRT_1、IRT_2)(单位:摄氏度)、土壤热通量(Gs_1、Gs_2、Gs_3)(单位:瓦/平方米)、土壤水分(Ms_5cm、Ms_10cm、Ms_20cm、Ms_40cm、Ms_80cm、Ms_120cm、Ms_200cm、Ms_300cm、Ms_400cm)(单位:百分比)、土壤温度(Ts_5cm、Ts_10cm、Ts_20cm、Ts_40cm、Ts_80cm、Ts_120cm、Ts_200cm、Ts_300cm、Ts_400cm)(单位:摄氏度) 、向上与向下光合有效辐射(PAR_U_up、PAR_U_down)(单位:微摩尔/平方米秒)。 观测数据的处理与质量控制:(1)确保每天144个数据(每10min),若出现数据的缺失,则由-6999标示;(2)剔除有重复记录的时刻;(3)删除了明显超出物理意义或超出仪器量程的数据;(4)数据中以红字标示的部分为有疑问的数据;(5)日期和时间的格式统一,并且日期、时间在同一列。如,时间为:2018/8/31 10:30。
李小雁
该数据集包含了2018年1月1日至2018年12月31日的黑河流域地表过程综合观测网下游混合林站涡动相关仪观测数据。站点位于内蒙古额济纳旗四道桥,下垫面是胡杨与柽柳。观测点的经纬度是101.1335E, 41.9903N,海拔874 m。涡动相关仪的架高22m,采样频率是10Hz,超声朝向是正北向,超声风速温度仪(CSAT3)与CO2/H2O分析仪(Li7500)之间的距离是17cm。 涡动相关仪的原始观测数据为10Hz,发布的数据是采用Eddypro软件处理的30分钟数据,其处理的主要步骤包括:野点值剔除,延迟时间校正,坐标旋转(二次坐标旋转),频率响应修正,超声虚温修正和密度(WPL)修正等。同时对各通量值进行质量评价,主要是大气平稳性(Δst)和湍流相似性特征(ITC)的检验。对Eddypro软件输出的30min通量值也进行了筛选:(1)剔除仪器出错时的数据;(2)剔除降水前后1h的数据;(3)剔除10Hz原始数据中每30min内缺失率大于10%的数据。观测数据的平均周期为30分钟,一天48个数据,缺失数据标记为-6999。2月7日-11日由于供电问题,导致数据缺失。 发布的观测数据包括:日期/时间Date/Time,风向Wdir(°),水平风速Wnd(m/s),侧向风速标准差Std_Uy(m/s),超声虚温Tv(℃),水汽密度H2O(g/m3),二氧化碳浓度CO2(mg/m3),摩擦速度Ustar(m/s),奥布霍夫长度L(m),感热通量Hs(W/m2),潜热通量LE(W/m2),二氧化碳通量Fc(mg/(m2s)),感热通量的质量标识QA_Hs,潜热通量的质量标识QA_LE,二氧化碳通量的质量标识QA_Fc。感热、潜热、二氧化碳通量的质量标识分为九级(质量标识1-3数据质量好,4-6数据质量较好,7-8数据质量较差(较插补数据好);9数据质量差))。数据时间的含义,如0:30代表0:00-0:30的平均;数据以*.xls格式存储。 黑河流域地表过程综合观测网或站点信息请参考Liu et al. (2018),观测数据处理请参考Liu et al. (2011)。
刘绍民, 李新, 车涛, 徐自为, 任志国, 谭俊磊
本数据集包括祁连山区域1990年至2017年每5年一期的30m土地覆盖分类产品。该产品首先利用Landsat-8/OLI构造2015年时间序列数据,针对各类地物随时间变化呈现的NDVI时间序列曲线不同,对不同地物特征进行知识归纳,设定提取规则不同地物信息,得到2015年的土地覆盖分类图。分类系统参考了IGBP分类系统和FROM_LC分类系统,共分为耕地、林地、草地、灌丛、湿地、水体、不透水面、裸地、冰川和积雪共10大类。由Google Earth高清影像和实地调研数据进行精度评价,得出2015年土地覆盖分类产品的总体精度高达92.19%。以2015年的土地覆盖分类产品为底图,按各类别的比例选取大量样本,基于Google Earth Engine平台的Landsat系列数据和强大地数据处理能力,利用深度学习的思想,选取随机森林分类器,对波段信息和NDVI、MNDWI、NDBI等指数进行样本训练,生产出1985-2017年每5年一期的土地覆盖分类产品。对2套2015年的分类产品进行比较,得出基于Google Earth Engine平台生产的土地覆盖分类产品与基于时间序列方法得到的分类产品具有很好的一致性。总之,祁连山核心区的土地覆盖数据集具有较高的总体精度,且基于Google Earth Engine平台样本训练的方法能够在时间和空间上对现有的分类产品进行扩展,每5年一期的频次能够在长时间序列上反映更多的土地覆盖类型变化信息。
吴俊君, 仲波, 角坤升
植被调查数据是研究生态系统结构与功能必不可少的数据。藏北地区蕴含广袤的草地生态系统,主要包括高寒草甸、高寒草地、以及高寒荒漠化的草地。由于独特的地理位置以及高海拔缺氧的环境条件,在藏北高原的群落调查数据较为稀少。本课题组基于前期工作的积累,在2017年生长季对整个藏北高原15个县域开展了较为全面的植被调查。本数据集包括藏北样带上从那曲到日土县23个采样点的围栏内外的生物量数据。本数据集可用于生产力的空间分析与模型的校准工作。
张宪洲, 牛犇
本数据集根据NASA EOSDIS LP DAAC 和美国地质调查 USGS EROS共同发布的第六版MODIS均一化植被指数产品(2001-2014)。NDVI的时间分辨率是16天,空间分辨率0.05度,我们先将NDVI数据产品从0.05度空间分辨率重采样到0.5度,然后对每年的时间序列采用double-logistic方法进行平滑,并计算平滑后的曲率,选取春季曲率最大值作为植被的春季返青期,该数据可分析泛北极植被春季物候的时空特征。
美国国家航空航天局, 徐希燕
降水强烈的时空变化常使得常规地基台站的降水观测不能准确把握降水的空间分布和强度变化。而卫星微波遥感可以克服此局限,实现全球尺度降水和云的观测,而且相对于红外/可见光只能反映云厚、云高等信息而言,微波能够穿透云体,利用云内降水粒子和云粒子与微波的相互作用对云、雨进行更为直接的探测。 本数据以GPM搭载的DPR双波段降水雷达获取的地表降水量为真值,以NDVI、DEM、ERA5中的土壤温/湿度为参考数据,利用GMI的多波段被动亮温数据反演青藏高原地区暖季(5月-9月)瞬时降水强度,将结果重采样至0.1°空间分辨率后累加到日。
许时光
该数据集为全球生态系统呼吸数据,包含生态系统自养呼吸(Ra)和异养呼吸(Rh)两部分,由耦合模式比较计划第6阶段(CMIP6)中CNRM-CM6-1模式在Historical情景下模拟得到。数据时间范围为1850-2014年,时间分辨率为月,空间分辨率约为1.406°×1.389°。模拟数据详细说明可见链接http://www.umr-cnrm.fr/cmip6/spip.php?article11。
美国气候模式诊断和对比计划委员会
本数据集来源于论文:Zhang, J. F., Xu, B., Turner, F., Zhou, L., Gao, P., Lü, X., & Nesje, A. (2017). Long-term glacier melt fluctuations over the past 2500 yr in monsoonal High Asia revealed by radiocarbon-dated lacustrine pollen concentrates. Geology, 45(4), 359-362. 在本文中,中国科学院青藏高原研究所、地球科学卓越创新中心徐柏青研究员及其博士后张继峰与来自北京大学等单位的合作者,对高原南部枪勇冰川冰前湖沉积物进行了多方法(植物残体、孢粉浓缩物、全有机质)放射性碳测年,提出了一个重建古冰川融化强度的新指标(“老孢粉效应”,即沉积物孢粉年龄与沉积物真实年龄的差值)。该研究发现北半球温度及西风环流活动可能是高原季风区冰川百年尺度波动的主控因素,高原近代的冰川融化强度达到过去2500年以来最强,超过了历史上的中世纪暖期和罗马暖期。 数据由论文作者提供,数据包含了基于老孢粉效应(ΔAgepollen)重建的过去2500年枪勇冰川融化强度变化数据。 研究人员从枪勇错冰前湖获得了一根3.06米长的湖芯(QYL09-4)和一根1.06m长的平行重力钻湖芯(QY-3),使用新的复合提取及纯化程序,从沉积物中获得了相对纯的孢粉浓缩物和植物残体浓缩物(PRC;> 125μm)。对全有机质,PRC和孢粉浓缩物分别进行了14C年代测定。所有14C年龄都使用IntCal13(Reimer et al., 2013)进行了校准。年龄深度模型基于210Pb、137Cs年龄及五个PRC的14C年龄。使用Oxcal 4.2(Bronk Ramsey,2008)中的P_Sequence算法构建岩芯的年龄深度模型。将校准的孢粉年龄中减去根据沉积模型得出的真实沉积物年龄,从而得出老孢粉效应值(ΔAgepollen)。 数据为湖芯(QYL09-4)的放射性碳测年与老孢粉效应数据。 数据包含字段如下: Lab No.:样本编号 Dating Material:测年材料 Depth (cm):深度(厘米) 14C age (yr BP):碳14年龄(年 距今) ΔAgepollen (≥95.4 % yrs):孢粉年龄与估算的沉积物年龄间的差值(≥95.4 % 年) Sediment Age (CE):沉积物年龄(公元) 数据详细信息参见附件:ZhangJF et al. 2017 GEOLOGY_Long-term glacier melt fluctuations over the past 2500 yr on the Tibetan Plateau.pdf。
张继峰
本数据集包括祁连山区域1986、1990、1995、2000、2005、2010、2015和2017年月度合成30m×30m地表植被指数产品。采用最大值合成 (Max value composition, MVC) 方法,利用 Landsat5, Landsat8和sentinel 2红光和近红外两个通道的反射率数据,实现对地表月度NDVI产品的合成。数据通过Google Earth Engine云平台对反射率进行月度合成,通过模型计算指数,对于缺失像素进行插补,质量较好,可用于环境变化监测等领域。
仲波, 吴俊君
该数据集包含了2018年1月1日至2018年12月31日黑河流域地表过程综合观测网上游大沙龙站气象要素观测数据。站点位于青海省祁连县西侧沙龙滩地区,下垫面是沼泽化高寒草甸。观测点的经纬度是98.9406E, 38.8399N,海拔3739m。空气温度、相对湿度传感器架设在5m处,朝向正北;气压计安装在地面上的防撬箱内;翻斗式雨量计安装在10m处;风速与风向传感器架设在10m,朝向正北;四分量辐射仪安装在6m处,朝向正南;两个红外温度计安装在6m处,朝向正南,探头朝向是垂直向下;土壤温度探头埋设在地表0cm和地下4cm、10cm、20cm、40cm、80cm、120cm、160cm处,在距离气象塔2m的正南方;土壤水分探头埋设在地下4cm、10cm、20cm、40cm、80cm、120cm、160cm处,在距离气象塔2m的正南方;土壤热流板(3块)依次埋设在地下6cm处,并距离气象塔2m的正南方。 观测项目有:空气温湿度(Ta_5m、RH_5m)(单位:摄氏度、百分比)、气压(Press)(单位:百帕)、降水量(Rain)(单位:毫米)、风速(WS_10m)(单位:米/秒)、风向(WD_10m)(单位:度)、四分量辐射(DR、UR、DLR_Cor、ULR_Cor、Rn)(单位:瓦/平方米)、地表辐射温度(IRT_1、IRT_2)(单位:摄氏度)、土壤热通量(Gs_1、Gs_2、Gs_3)(单位:瓦/平方米)、土壤温度(Ts_0cm、Ts_4cm、Ts_10cm、Ts_20cm、Ts_40cm、Ts_80cm、Ts_120cm、Ts_160cm)(单位:摄氏度)、土壤水分(Ms_4cm、Ms_10cm、Ms_20cm、Ms_40cm、Ms_80cm、Ms_120cm、Ms_160cm)(单位:体积含水量,百分比)。 观测数据的处理与质量控制:(1)确保每天144个数据(每10min),若出现数据的缺失,则由-6999标示;(2)剔除有重复记录的时刻;(3)删除了明显超出物理意义或超出仪器量程的数据;(4)数据中以红字标示的部分为有疑问的数据;(5)日期和时间的格式统一,并且日期、时间在同一列。如,时间为:2018-9-10 10:30;(6)命名规则为:AWS+站点名称。 黑河综合观测网或站点信息请参考Liu et al. (2018),观测数据处理请参考Liu et al. (2011)。
刘绍民, 李新, 车涛, 谭俊磊, 任志国, 张阳, 徐自为
该数据集包含了2018年1月1日至2018年10月12日青海湖流域水文气象观测网青海湖鱼雷发射基地站气象要素梯度观测系统数据。站点位于青海省青海湖二郎剑景区鱼雷发射基地,下垫面是青海湖水面。观测点经纬度为:东经 100° 29' 59.726'' E,北纬 36° 35' 27.337'' N,海拔3209m。风速/风向架设在距湖面14m处,共1层,朝向正北;空气温度、相对湿度传感器分别架设在距湖面12m、12.5m处,共2层,朝向正北;翻斗式雨量计安装在距湖面10m处;四分量辐射仪安装在距湖面10m处,朝向正南;一个红外温度计安装在距湖面10m处,朝向正南,探头朝向是垂直向下;湖水温度探头设在水下0.2, 0.5, 1.0, 2.0, and 3.0 m处;光合有效辐射仪安装在距湖面10m处,探头朝向是垂直向下,朝向正南。 观测项目有:风速(WS_14m)(单位:米/秒)、风向(WD_14m)(单位:度)、空气温湿度(Ta_12m、Ta_12.5m和RH_12m、RH_12.5m)(单位:摄氏度、百分比)、降水量(Rain)(单位:毫米)、四分量辐射(DR、UR、DLR_Cor、ULR_Cor、Rn)(单位:瓦/平方米)、湖表辐射温度(IRT_1)(单位:摄氏度)、湖水温度(Tw_20cm、Tw_50cm、Tw_100cm、Tw_200cm、Tw_300cm)(单位:摄氏度) 、光合有效辐射(PAR)(单位:微摩尔/平方米秒)。 观测数据的处理与质量控制:(1)确保每天144个数据(每10min),若出现数据的缺失,则由-6999标示;2018.1.1-10.12由于由于采集器的问题,除四分量外的气象数据均无记录;(2)剔除有重复记录的时刻;(3)删除了明显超出物理意义或超出仪器量程的数据;(4)数据中以红字标示的部分为有疑问的数据;(5)日期和时间的格式统一,并且日期、时间在同一列。如,时间为:2018-1-1 10:30。
李小雁
本数据集包括祁连山区域1986、1990、1995、2000、2005、2010、2015和2017年月度合成30m×30m地表植被覆盖度产品。采用最大值合成 (Max value composition, MVC) 方法,利用 Landsat5, Landsat8和sentinel 2红光和近红外两个通道的反射率数据,实现对地表月度NDVI产品的合成,进而计算FVC。数据通过Google Earth Engine云平台对反射率进行月度合成,通过模型计算指数,对于缺失像素进行插补,质量较好,可用于环境变化监测等领域。
仲波, 吴俊君
本数据集包括祁连山区域2018年月度合成30m×30m地表植被指数产品。采用最大值合成 (Max value composition, MVC) 方法,利用Landsat 8和sentinel 2红光和近红外两个通道的反射率数据,实现对地表月度NDVI产品的合成。数据通过Google Earth Engine云平台对反射率进行月度合成,通过模型计算指数,对于缺失像素进行插补,质量较好,可用于环境变化监测等领域。
吴金华, 仲波, 吴俊君
本数据集包括祁连山区域2018年的30m耕地和建筑用地分布产品。该产品基于Landsat-8/OLI数据生产。生产耕地产品时,利用夏季合成的OLI数据,使用归一化植被指数NDVI,并辅助以作物的物候期、作物种植类别等先验知识设定产品生产规则;生产建筑用地产品时,利用夏季合成的Landsat系列数据,使用归一化植被指数NDVI、NDBI、MNDWI等指数,辅助以DEM、灯光数据等设定产品生产规则。由Google Earth高清影像和实地调研数据进行精度评价,得出2018年耕地和建筑用地产品的精度分别为90.05%和90.97%。
吴俊君, 仲波, 角坤升
数据集包括:北极地区人口及GDP数据(1990-2015)、第三极(甘肃、青海、西藏)地区县级人口及GDP数据(1970-2016)。 社会经济统计属性包括:人口(万人)、GDP(万元)、工农业生产总值(万元)、农业总产值(万元)、工业生产总值(万元) 北极人口数据主要来自经济社会局《世界人口展望:2017年修订版》按照 区域和国家划分的人口总数。 第三极数据主要参考甘肃省统计年鉴、青海省统计年鉴、西藏自治区统计年鉴;甘肃省、青海省、西藏自治区各县县志。
经济和社会事务部, 国家统计局, 青海省统计局
该数据集是三极地区9个通量站点的30分钟涡度相关通量观测数据,包括生态系统净碳交换量(NEE)、总初级生产力(GPP)和生态系统呼吸(ER)数据,时间范围覆盖2000-2016年。该数据预处理主要步骤包括野点去除(±3σ)、坐标轴旋转(三维风旋转)、Webb-Pearman-Leuning校正、异常值剔除、碳通量插补与分解等,缺失数据通过CO2通量值(Fc)与环境因子之间的非线性经验公式进行插补。
张扬建, 牛犇
该数据集包含了2018年1月1日至2018年12月31日黑河流域地表过程综合观测网中游花寨子荒漠站气象要素观测数据。站点位于甘肃省张掖市花寨子,下垫面是盐爪爪山前荒漠。观测点的经纬度是100.3201E, 38.7659N,海拔1731m。空气温度、相对湿度传感器架设在5m、10m处,朝向正北;气压计安装防水箱内;翻斗式雨量计安装在10m处;风速风向传感器架设在5m、10m处,朝向正北;四分量辐射仪安装在6m处,朝向正南;两个红外温度计安装在6m处,朝向正南,探头朝向是垂直向下;土壤温度探头埋设在地表0cm和地下2cm、4cm、10cm、20cm、40cm、60cm和100cm处,在距离气象塔2m的正南方;土壤水分传感器分别埋设在地下2cm、4cm、10cm、20cm、40cm、60cm和100cm处,在距离气象塔2m的正南方;土壤热流板(3块)依次埋设在地下6cm处。 观测项目有:空气温湿度(Ta_5m、RH_5m、Ta_10m、RH_10m)(单位:摄氏度、百分比)、气压(Press)(单位:百帕)、降水量(Rain)(单位:毫米)、风速(WS_5m、WS_10m)(单位:米/秒)、风向(WD_5m、WD_10m)(单位:度)、四分量辐射(DR、UR、DLR_Cor、ULR_Cor、Rn)(单位:瓦/平方米)、地表辐射温度(IRT_1、IRT_2)(单位:摄氏度)、土壤热通量(Gs_1、Gs_2、Gs_3)(单位:瓦/平方米)、土壤水分(Ms_2cm、Ms_4cm、Ms_10cm、Ms_20cm、Ms_40cm、Ms_60cm、Ms_100cm)(单位:体积含水量,百分比)和土壤温度(Ts_0cm、Ts_2cm、Ts_4cm、Ts_10cm、Ts_20cm、Ts_40cm、Ts_60cm、Ts_100cm)(单位:摄氏度)。 观测数据的处理与质量控制:(1)确保每天144个数据(每10min),若出现数据的缺失,则由-6999标示;1.1-1.7及8.22-8.31,9.4-9.12间由于传感器的问题,4cm土壤水分出错;(2)剔除有重复记录的时刻;(3)删除了明显超出物理意义或超出仪器量程的数据;(4)数据中以红字标示的部分为有疑问的数据;(5)日期和时间的格式统一,并且日期、时间在同一列。如,时间为:2018-6-10 10:30;(6)命名规则为:AWS+站点名称 黑河综合观测网或站点信息请参考Liu et al. (2018),观测数据处理请参考Liu et al. (2011)。
刘绍民, 李新, 车涛, 徐自为, 任志国, 谭俊磊
本数据集包括祁连山区域1986、1990、1995、2000、2005、2010、2015和2017年月度合成30m×30m地表NPP产品。采用最大值合成 (Max value composition, MVC) 方法,利用 Landsat5, Landsat8和sentinel 2红光和近红外两个通道的反射率数据,实现对地表月度NDVI产品的合成,进而计算NPP。 数据通过Google Earth Engine云平台对反射率进行月度合成,通过模型计算指数,对于缺失像素进行插补,质量较好。
吴金华, 仲波
联系方式
关注我们
时空三极环境大数据平台 © 2018-2020 陇ICP备05000491号 | All Rights Reserved
| 京公网安备11010502040845号
数据中心技术支持: 数云软件