This data set includes the continuous observation data set of soil texture, roughness and surface temperature measured by vehicle borne microwave radiometer from November 17 to 18, 2013 in Wuxing village farmland, Ganzhou District, Zhangye City, Gansu Province. The surface temperature and humidity include four layers of temperature sensor at the soil depth of 1cm, 5cm, 10cm, 20cm, and the observation of soil temperature and soil moisture data at the soil depth of 0-5cm. The time frequency of routine observation of soil temperature and humidity is 5 minutes. Data details: 1. Time: November 17-18, 2013 2. data: Brightness temperature: observed by vehicle mounted multi frequency passive microwave radiometer, including 6.925, 18.7 and 36.5ghz V polarization and H polarization data (10.65ghz band damage) Soil temperature: use sensor installed on dt80 to measure 1cm, 5cm, 10cm, 20cm soil temperature Soil moisture: use h-probe sensor to measure 0-5cm soil moisture, the probe can measure 0-5cm soil temperature at the same time Soil texture: soil samples measured in Beijing Normal University Soil roughness: measured by roughness meter provided by northeast geography 3. Data size: 3.6m 4. Data format:. Xls
ZHAO Shaojie, KOU Xiaokang, YE Qinyu, MA Mingguo
This data set includes the continuous observation data set of soil texture, roughness and surface temperature measured by the vehicle borne microwave radiometer on November 15-16, 2013 in the farmland of jiushe, Kangning, Zhangye City, Gansu Province. The surface temperature includes the soil temperature data observed by the temperature sensor at the soil depth of 0 cm, 1 cm, 3 cm, 5 cm and 10 cm. The time frequency of conventional observation of soil temperature is 5 minutes. Data details: 1. Time: November 15-16, 2013 2. data: Bright temperature: observed by vehicle mounted multi frequency passive microwave radiometer, including 6.925, 18.7 and 36.5ghz v-polarization and H-polarization data (10.65ghz band instrument damaged) Soil temperature: use the sensor installed on dt85 to measure the soil temperature of 0cm, 1cm, 3cm, 5cm and 10cm Soil texture: soil samples measured in Beijing Normal University Soil roughness: measured by roughness meter provided by northeast geography 3. Data size: 4.8m 4. Data format:. Xls
ZHAO Shaojie, KOU Xiaokang, YE Qinyu, MA Mingguo
On July 26, 2012, the airborne ground synchronous observation was carried out in the plmr quadrat in the dense observation area of Daman. Plmr (polarimetric L-band multibeam radiometer) is a dual polarized (H / V) L-band microwave radiometer, with a center frequency of 1.413 GHz, a bandwidth of 24 MHz, a resolution of 1 km (relative altitude of 3 km), six beam simultaneous observations, an incidence angle of ± 7 °, ± 21.5 °, ± 38.5 °, and a sensitivity of < 1K. The flight mainly covers the middle reaches of the artificial oasis eco hydrological experimental area. The local synchronous data set can provide the basic ground data set for the development and verification of passive microwave remote sensing soil moisture inversion algorithm. Quadrat and sampling strategy: The observation area is located in the matrix of the dense observation area of Daman, and the detailed plan with an area of 3.0KM × 2.4km is selected to carry out synchronous observation on the underlying surface of oasis. The selection of the sample is mainly based on the representativeness of the surface coverage, accessibility and observation (road consumption) time, so as to obtain the comparison of brightness and temperature with plmr observation. Considering the resolution of plmr observation, 5 splines (east-west distribution) were collected at an interval of 450 m in the east-west direction. Each line has 31 points (north-south direction) at an interval of 100 m, and 5 hydraprobe data acquisition systems (HDAS, reference 2) were used for simultaneous measurement. Measurement content: About 150 points on the quadrat were obtained, each point was observed twice, that is to say, two times were observed at each sampling point, one time was inside the film (marked as a in the data record) and one time was outside the film (marked as B in the data record). As the HDAS system uses pogo portable soil sensor, the soil temperature, soil moisture (volume moisture content), loss tangent, soil conductivity, real part and imaginary part of soil complex dielectric are observed. Because the vegetation in this area has been sampled and observed once every five days, no special vegetation synchronous sampling has been carried out on that day. Data: This data set consists of two parts: soil moisture observation and vegetation observation. The former saves data in vector file format, and the spatial location is the location of each sampling point (WGS84 + UTM 47N). Soil moisture and other measurement information are recorded in attribute file.
WANG Shuguo, MA Mingguo, LI Xin
Er’ba Reservoir surface temperature of water body can offer in situ calibration data for TASI, WiDAS and L band sensor used in aerospace experiment. Observation Site: This site is 14 KM away from East of ZhangYe city. It’s located in Er’ba village, JianTan town, ZhangYe city. The coordinates of this site: 38°54′57.14" N, 100°36′57.39" E. Observation Instrument: The observation system consists of two SI-111 infrared radiometers (Campbell, USA) and two 109SS temperature probes (Campbell, USA). Two SI-111 sensors, one installed vertically downward to water surface, another face to south of zenith angle 35°. Temperature probes float under water surface at 0 cm. SI-111 sensor installed at 3.0 m height, 3.4 m away from water edge. Observation Time: This site operates from 27 May, 2012 to 27 September, 2012. Observation data laagered by every 5 seconds uninterrupted. Output data contained sample data of every 5 seconds and mean data of 1 minute. Accessory data： Water surface infrared temperature (by SI-111), sky infrared temperature (by SI-111), water surface temperature (by 109ss) can be obtained. Dataset is stored in *.dat file, which can be read by Microsoft excel or other text processing software (UltraEdit, et. al). Table heads meaning: TarT_Atm, Sky infrared temperature (℃) @ facing south of zenith angle 35°; SBT_Atm, body temperature of SI-111 sensor (℃) measured sky; TarT_Sur, water surface infrared temperature @ 3.0 m height; SBT_Sur, body temperature of SI-111 sensor (℃) measured water surface; WaterT_1, WaterT_2, water surface temperature (℃) measured by 109SS temperature probes. Dataset is stored day by day, named as: data format + site name + interval time + date + time. The detailed information about data item showed in data header introduction in dataset.
MA Mingguo
On July 10, 2012, the airborne flight and ground observation was synchronously carried out in the PLMR quadrat of Yingke Oasis and the Huazhaizi Desert. PLMR (Polarimetric L-band Multibeam Radiometer) is a dual-polarized (H/V) L-band microwave radiometer with a center frequency of 1.413 GHz, a bandwidth of 24 MHz, and a resolution of 1 km (relative flight height of 3 km).The radiometer has 6 beams to observe synchronously, and the incident angles are ±7º，±21.5º，±38.5º, and the sensitivity is less than 1K. The flight observation mainly covers the artificial oasis eco-hydrological test area in the middle reaches. This ground-synchronized data set provides a basic ground dataset for developing and validating passive microwave remote sensing inversion soil moisture algorithms. Quadrat and sampling strategy: The observation area is located in the transition zone between the southern margin of Zhangye Oasis and Anyang beach desert, the west side of Zhang (Zhangye)-Da (Daman) highway. It is divided into two parts by the main canal of the Dragon Canal from North to South. The Southwest area is a desert quadrat with the size of 1 km×1 km. The desert is relatively homogeneous, so soil moisture of 5 points in the 1 km quadrat are collected (1 point of each corner and the center point, in the actual measurement process, several extra points can be measured along the road). The four corner points are 600 meters away from each other,except the diagonal direction. And the southwest corner point is Huazhaizi Desert Station, for the convenience of comparison with weather station data. On the northeast side, a large size quadrat of 6 km×1.6 km is selected for simultaneous observation of the oasis underlying surface.In order to obtain the brightness temperature comparison with the PLMR observation, the quadrat was chose based on the following factors :surface coverage representative, avoiding the residential and greenhouses, crossing the oasis farmland and part of the Southern desert, accessibility, and observation time(road consumption). Taking the resolution of PLMR observations into consideration, in the synchronous observation, 11 sampling lines (East-West distribution) were collected with an interval of 160 meters from the East to the West. Each line from the North to the South was separated by 21 points with an interval of 80 meters. And 4 Hydraprobe Data Acquisition System (HDAS, Reference 2) were used to measure simultaneously. Measurement contents: About 230 points of the quadrat were collected, 2 observations were performed on each point, that is, 2 observations were performed on each sampling point of the film mulched corn field, 1 inside the film (marked as a in the data record), 1 outside the film (marked as b in the data record). Since the HDAS system useed the POGO portable soil sensor, the soil temperature, soil moisture (volumetric water content), loss tangent, soil electrical conductivity, soil complex dielectric real part and imaginary part were obtained by observation. No special simultaneous sampling of vegetation was carried out on the same day. Data: The data set includes two parts: soil moisture observation and vegetation observation. The former saves the data as a vector file, the spatial position is the position of each sampling point (WGS84+UTM 47N), and the measurement information of soil moisture is recorded in the attribute file.
WANG Shuguo, LI Xin
On 30 June 2012 (UTC+8), TASI sensor carried by the Harbin Y-12 aircraft was used in a visible near Infrared hyperspectral airborne remote sensing experiment, which is located in the observation experimental area (30×30 km), Linze region and Heihe riverway. The relative flight altitude is 2500 meters. The wavelength of TASI is 8-11.5 μm with a spatial resolution of 3 meters. Through the ground sample points and atmospheric data, the data are recorded in surface radiance processed by geometric correction and atmospheric correction. Land surface temperature (LST) data was retrieved by temperature/emissivity separation algorithm.
XIAO Qing, Wen Jianguang
On 2 August 2012 (UTC+8), a Wide-angle Infrared Dual-mode line/area Array Scanner (WIDAS) carried by the Harbin Y-12 aircraft was used in a visible near Infrared thermal Dual-mode airborne remote sensing experiment, which is located in the artificial oasis eco-hydrology experimental area (30×30 km). WIDAS includes an CCD cameras with spatial resolution 0.26 m, a visible near Infrared multispectral camera with five bands scanner (an maximum observation angle 48° and spatial resolution 1.3 m), and a thermal image camera with spatial resolution 6.3 m. The CCD camera data production are recorded in DN values processed by mosaic and orthorectification. The mutispectral camera data production are recorded in reflectance processed by atmospheric and geometric correction. Thermal image camera data production are recorded in radiation brightness temperature processed by atmospheric and geometric correction.
XIAO Qing, Wen Jianguang
The dataset of airborne Polarimetric L-band Multibeam Radiometers (PLMR) was acquired on 10 July, 2012, located in the middle reaches of the Heihe River Basin. The aircraft took off at 10:30 am (UTC+8) from Zhangye airport and landed at 15:30 pm, with the flight time of 5 hours. The flight was performed in the altitude of about 2500 m and at the speed of about 220-250 km during the observation, corresponding to an expected ground resolution of about 750 m. The PLMR instrument flown on a small aircraft operates at 1.413 GHz (L-band), with both H- and V-polarizations at incidence angles of ±7.5°, ±21.5° and ±38.5°. PLMR ‘warm’ and ‘cold’ calibrations were performed before and after each flight. The processed PLMR data include 2 DAT files (v-pol and h-pol separately) and 1 KMZ file for each flying day. The DAT file contains all the TB values together with their corresponding beam ID, incidence angle, location, time stamp (in UTC) and other flight attitude information as per headings. The KMZ file shows the gridded 1-km TB values corrected to 38.5 degrees together with flight lines. Cautions should be taken when using these data, as the RFI contaminations are often higher than expected at v-polarization.
CHE Tao, Gao Ying, LI Xin
On 29 June 2012 (UTC+8), a CASI/SASI sensor carried by the Harbin Y-12 aircraft was used in a visible near Infrared hyperspectral airborne remote sensing experiment, which is located in the observation experimental area (30×30 km). The relative flight altitude is 3500 meters（an elevation of 3500 meters）, The wavelength of CASI and SASI is 380-1050 nm and 950-2450 nm, respectively. The spatial resolution of CASI and SASI is 1 m and 2.4 m, respectively. Through the ground sample points and atmospheric data, the data are recorded in reflectance processed by geometric correction and atmospheric correction based on 6S model.
XIAO Qing, Wen Jianguang
The dataset of airborne Polarimetric L-band Multibeam Radiometers (PLMR) was acquired on 26 July, 2012, located in the middle reaches of the Heihe River Basin. The aircraft took off at 9:10 am (UTC+8) from Zhangye airport and landed at 13:40 pm, with the flight time of 4.5 hours. The flight was performed in the altitude of about 2300 m and at the speed of about 220-250 km during the observation, corresponding to an expected ground resolution of about 700 m. The PLMR instrument flown on a small aircraft operates at 1.413 GHz (L-band), with both H- and V-polarizations at incidence angles of ±7.5°, ±21.5° and ±38.5°. PLMR ‘warm’ and ‘cold’ calibrations were performed before and after each flight. The processed PLMR data include 2 DAT files (v-pol and h-pol separately) and 1 KMZ file for each flying day. The DAT file contains all the TB values together with their corresponding beam ID, incidence angle, location, time stamp (in UTC) and other flight attitude information as per headings. The KMZ file shows the gridded 1-km TB values corrected to 38.5 degrees together with flight lines. Cautions should be taken when using these data, as the RFI contaminations are often higher than expected at v-polarization.
CHE Tao, Gao Ying, LI Xin
The aim of the simultaneous observation of land surface temperature is obtaining the land surface temperature of different kinds of underlying surface, including greenhouse film, the roof, road, ditch, concrete floor and so on, while the sensor of thermal infrared go into the experimental areas of artificial oases eco-hydrology on the middle stream. All the land surface temperature data will be used for validation of the retrieved land surface temperature from thermal infrared sensor and the analysis of the scale effect of the land surface temperature, and finally serve for the validation of the plausibility checks of the surface temperature product from remote sensing. 1. Observation time and other details On 25 June, 2012, ditch and asphalt road surface temperatures were observed once every five minutes using handheld infrared thermometers recorded. On 26 June, 2012, ditch and asphalt road surface temperatures were observed once every five minutes using handheld infrared thermometers while greenhouse film and concrete floor surface temperatures were observed once every one second using self-recording point thermometer. On 29 June, 2012, concrete floor surface temperatures were observed continuously using handheld infrared thermometers during the sensor of TASI go into the region. At the same time, greenhouse film and concrete floor surface temperatures were observed once every one second using self-recording point thermometer. On 30 June, 2012, asphalt road, ditch, bare soil, melonry and ridge of field surface temperatures were observed continuously using handheld infrared thermometers during the sensor of TASI go into the region. At the same time, greenhouse film and concrete floor surface temperatures were observed once every one second using self-recording point thermometer. On 10 July, 2012, asphalt road, ditch, bare soil, melonry and ridge of field surface temperatures were observed once every one minute using handheld infrared thermometers during the sensor of TASI go into the region. At the same time, concrete floor surface temperatures were observed once every six second using self-recording point thermometer. On 26 July, 2012, asphalt road, concrete floor, bare soil and melonry surface temperatures were observed once every one minute using handheld infrared thermometers during the sensor of WiDAS go into the region. At the same time, greenhouse film and concrete floor surface temperatures were observed once every six second using self-recording point thermometer. On 2 August, 2012, corn field and concrete floor surface temperatures were observed using handheld infrared thermometers. At the same time, greenhouse film and concrete floor surface temperatures were observed once every six second using self-recording point thermometer. For corn field, twelve sites were selected according to the flight strip of the WiDAS sensor, and for each site one plot surface temperatures were recorded continuously during the sensor of WiDAS go into the region. On 3 August, 2012, corn field and concrete floor surface temperatures were observed using handheld infrared thermometers. At the same time, greenhouse film and concrete floor surface temperatures were observed once every six second using self-recording point thermometer. For corn field, fourteen sites were selected according to the flight strip of the WiDAS sensor, and for each site three plots surface temperatures were recorded continuously during the sensor of WiDAS go into the region. 2. Instrument parameters and calibration The field of view of the self-recording point thermometer and the handheld infrared thermometer are 10 and 1 degree, respectively. The emissivity of the latter was assumed to be 0.95. The observation heights of the self-recording point thermometer for the greenhouse film and the concrete floor were 0.5 m and 1 m, respectively. All instruments were calibrated three times (on 6 July, 5 August and 20 September, 2012) using black body during observation. 3. Data storage All the observation data were stored in excel.
GENG Liying, Jia Shuzhen, WANG Haibo, PENG Li, Dong Cunhui
Data content: precipitation data of the Aral Sea basin from 2015 to 2018. Data sources and processing methods: from the new generation of global precipitation measurement (GPM) of NASA (version 06, global precipitation observation program), the daily rainfall can be obtained by adding the three-hour rainfall data, and then the eight day rainfall can be obtained. Data quality: the spatial resolution is 0.1 ° x 0.1 ° and the temporal resolution is 8 days. The value of each pixel is the sum of rainfall in 8 days. Data application results: under the background of climate change, it can be used to analyze the correlation between meteorological elements and vegetation characteristics.
XIAO Qing, Wen Jianguang
On 4 July 2012 (UTC+8), a TASI sensor boarded on the Y-12 aircraft was used to obtain the thermal-infrared hyperspectral image, which is located in the observation experimental area, Linze region and Heihe riverway. The relative flight altitude is 1000 meters. The wavelength of TASI is 8-11.5 μm with a spatial resolution of 3 meters. Through the ground sample points and atmospheric data, the data are recorded in surface radiance processed by geometric correction and atmospheric correction. Land surface temperature (LST) data was retrieved by temperature/emissivity separation algorithm.
XIAO Qing, Wen Jianguang
The dataset of airborne Polarimetric L-band Multibeam Radiometers (PLMR) was acquired on 7 July, 2012, located in the middle reaches of the Heihe River Basin. The aircraft took off at 13:40 pm (UTC+8) from Zhangye airport and landed at 17:40 pm, with the flight time of 4 hours. The flight was performed in the altitude of about 2000 m and at the speed of about 220-250 km during the observation, corresponding to an expected ground resolution of about 600 m. The PLMR instrument flown on a small aircraft operates at 1.413 GHz (L-band), with both H- and V-polarizations at incidence angles of ±7.5°, ±21.5° and ±38.5°. PLMR ‘warm’ and ‘cold’ calibrations were performed before and after each flight. The processed PLMR data include 2 DAT files (v-pol and h-pol separately) and 1 KMZ file for each flying day. The DAT file contains all the TB values together with their corresponding beam ID, incidence angle, location, time stamp (in UTC) and other flight attitude information as per headings. The KMZ file shows the gridded 1-km TB values corrected to 38.5 degrees together with flight lines. Cautions should be taken when using these data, as the RFI contaminations are often higher than expected at v-polarization.
CHE Tao, Gao Ying, LI Xin
On 26 July 2012, a Wide-angle Infrared Dual-mode line/area Array Scanner (WIDAS) carried by the Y-12 aircraft was used in a visible near Infrared thermal Dual-mode airborne remote sensing experiment, which is located in the artificial oasis eco-hydrology experimental area (5×5 km). WIDAS includes an CCD cameras with spatial resolution 0.2 m, a visible near Infrared multispectral camera with five bands scanner (an maximum observation angle 48° and spatial resolution 1 m), and a thermal image camera with spatial resolution 4.8 m. The CCD camera data are recorded in DN values processed by mosaic and orthorectification. The multispectral camera data are recorded in reflectance processed by atmospheric and geometric correction. Thermal image camera data are recorded in radiation brightness temperature processed by atmospheric and geometric correction.
XIAO Qing, Wen Jianguang
This dataset contains the flux measurements from the Zhangye wetland station eddy covariance system (EC) in the flux observation matrix from 25 June to 26 September, 2012. The site (100.44640° E, 38.97514° N) was located in a wetland surface, which is near Zhangye city, Gansu Province. The elevation is 1460.00 m. The EC was installed at a height of 5.2 m; the sampling rate was 10 Hz. The sonic anemometer faced north, and the separation distance between the sonic anemometer and the CO2/H2O gas analyzer (Gill&Li7500A) was 0.25 m. Raw data acquired at 10 Hz were processed using the Eddypro post-processing software (Li-Cor Company, http://www.licor.com/env/products/ eddy_covariance/software.html), including spike detection, lag correction of H2O/CO2 relative to the vertical wind component, sonic virtual temperature correction, angle of attack correction, coordinate rotation (2-D rotation), corrections for density fluctuation (Webb-Pearman-Leuning correction), and frequency response correction. The EC data were subsequently averaged over 30 min periods. Moreover, the observation data quality was divided into three classes according to the quality assessment method of stationarity (Δst) and the integral turbulent characteristics test (ITC), which was proposed by Foken and Wichura [1996]: class 1 (level 0: Δst<30 and ITC<30), class 2 (level 1: Δst<100 and ITC<100), and class 3 (level 2: Δst>100 and ITC>100), representing high-, medium-, and low-quality data, respectively. In addition to the above processing steps, the half-hourly flux data were screened in a four-step procedure: (1) data from periods of sensor malfunction were rejected; (2) data before or after 1 h of precipitation were rejected; (3) incomplete 30 min data were rejected when the missing data constituted more than 3% of the 30 min raw record; and (4) data were rejected at night when the friction velocity (u*) was less than 0.1 m/s. There were 48 records per day; the missing data were replaced with -6999. Moreover, suspicious data were marked in red. The released data contained the following variables: data/time, wind direction (Wdir, °), wind speed (Wnd, m/s), the standard deviation of the lateral wind (Std_Uy, m/s), virtual temperature (Tv, ℃), H2O mass density (H2O, g/m^3), CO2 mass density (CO2, mg/m^3), friction velocity (ustar, m/s), stability (z/L), sensible heat flux (Hs, W/m^2), latent heat flux (LE, W/m^2), carbon dioxide flux (Fc, mg/ (m^2s)), quality assessment of the sensible heat flux (QA_Hs), quality assessment of the latent heat flux (QA_LE), and quality assessment of the carbon flux (QA_Fc). In this dataset, the time of 0:30 corresponds to the average data for the period between 0:00 and 0:30; the data were stored in *.xlsx format. For more information, please refer to Liu et al. (2016) (for multi-scale observation experiment or sites information), Xu et al. (2013) (for data processing) in the Citation section.
LIU Shaomin, LI Xin, XU Ziwei
Trough the select tasking, we obtained the WorldView-2 stereo image data in Dayekou Basin production in mid-May 2012. In the same year from July to August, 27 GPS ground control points (GCP) and checkpoints were measured based on the watershed differential GPS control network. Based on the full-field GCPs, the rational polynomial coefficients (RPC) files of WorldView-2 images were corrected in the digital photogrammetry software system. In the stereo model, 60 high-precision tie points evenly distributed were got through image matching technology, and the 1-m and 2-m resolution digital elevation model (DEM) were rapid extracted. Moreover, the DEM was edited in some key areas, such as the shady forest coverage and Dayekou reservoir. The terrain feature points and line data were added to improve the accuracy of the results in large variation of terrain feature. Check points were composed of GPS points and model confidential points, which used for quantitative validation. And they root mean square errors RMSE were 1.9 meters and 1.2 meters respectively, which achieve the requirements of two degree accuracy of 2.0 m at a scale of 1:2000 in high mountains.
Zhang Yanli, MA Mingguo
Trough the select tasking, we obtained the WorldView-2 stereo image data in Dayekou Watershed production in mid-May 2012. In the same year from July to August, 27 GPS ground control points (GCP) and checkpoints were measured based on the watershed differential GPS control network. Based on the full-field GCPs, the rational polynomial coefficients (RPC) files of WorldView-2 images were corrected in the digital photogrammetry software system. In the stereo model, 60 high-precision tie points evenly distributed were got through image matching technology, and the 1-m and 2-m resolution digital elevation model (DEM) were rapid extracted. Based on collinearity equations, images at nadir were corrected to adjust relief displacements and geometric errors, and the 0.5-m resolution digital orthorectified images DOM were obtained with the principle of digital differential rectification in Dayekou Basin.
Zhang Yanli, MA Mingguo
The dataset include the planting structure and area information of major crops in 11 districts and counties of the Heihe River Basin from 2000 to 2012 (grain, wheat, corn, potato, soybean, cotton, oil, vegetables, etc.)
DENG XiangZheng
We produced surface photosynthetic effective radiation (PAR), solar radiation (SSR) and net radiation (NR) products with 1KM resolution in the heihe basin in 2012.The temporal resolution ranges from instantaneous to hourly and daily.Day-by-day ancillary data were also produced, including aerosol optical thickness, moisture content, NDVI, snow cover, and surface albedo.Among them, PAR and SSR use the method of lookup table to directly invert by combining the stationary weather satellite and polar orbit satellite MODIS product.NR was calculated by analyzing the relationship between net short-wave and net surface radiation.Hourly instantaneous products are weighted by average and integral to obtain hourly and daily cumulative products.
HUANG Guanghui
Contact Support
Northwest Institute of Eco-Environment and Resources, CAS
0931-4967287
poles@itpcas.ac.cn
Links
National Tibetan Plateau Data CenterFollow Us
A Big Earth Data Platform for Three Poles © 2018-2020 No.05000491 | All Rights Reserved
｜ 
No.11010502040845
Tech Support: westdc.cn
