首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 390 毫秒
1.
In ground-based GPS meteorology, weighted mean temperature is the key parameter to calculate the conversion factor which will be used to map zenith wet delay to precipitable water vapor. In practical applications, we can hardly obtain the vertical profiles of meteorological parameters over the site, thus cannot use the integration method to calculate weighted mean temperature. In order to exactly calculate weighted mean temperature from a few meteorological parameters, this paper studied the relation between weighted mean temperature and surface temperature, surface water vapor pressure and surface pressure, and determined the relationship between, on the one hand, the weighted mean temperature, and, on the other hand, the surface temperature and surface water vapor pressure. Considering the seasonal and geographic variations in the relationship, we employed the trigonometry functions with an annual cycle and a semi-annual cycle to fit the residuals (seasonal and geographic variations are reflected in the residuals). Through the above work, we finally established the GTm-I model and the PTm-I model with a $2^{\circ }\times 2.5^{\circ }(\mathrm{lat}\times \mathrm{lon})$ resolution. Test results show that the two models both show a consistent high accuracy around the globe, which is about 1.0 K superior to the widely used Bevis weighted mean temperature–surface temperature relationship in terms of root mean square error.  相似文献   

2.
ABSTRACT

Digitizing the land surface temperature (Ts) and surface soil moisture (mv) is essential for developing the intelligent Digital Earth. Here, we developed a two parameter physical-based passive microwave remote sensing model for jointly retrieving Ts and mv using the dual-polarized Tb of Aqua satellite advanced microwave scanning radiometer (AMSR-E) C-band (6.9 GHz) based on the simplified radiative transfer equation. Validation using in situ Ts and mv in southern China showed the average root mean square errors (RMSE) of Ts and mv retrievals reach 2.42 K (R2 = 0.61, n = 351) and 0.025 g cm?3 (R2 = 0.68, n = 663), respectively. The results were also validated using global in situ Ts (n = 2362) and mv (n = 1657) of International Soil Moisture Network. The corresponding RMSE are 3.44 k (R2 = 0.86) and 0.039 g cm?3 (R2 = 0.83), respectively. The monthly variations of model-derived Ts and mv are highly consistent with those of the Moderate Resolution Imaging Spectroradiometer Ts (R2 = 0.57; RMSE = 2.91 k) and ECV_SM mv (R2 = 0.51; RMSE = 0.045 g cm?3), respectively. Overall, this paper indicates an effective way to jointly modeling Ts and mv using passive microwave remote sensing.  相似文献   

3.
Determination of Geopotential of Local Vertical Datum Surface   总被引:1,自引:0,他引:1  
1 IntroductionEachcountryoreachgroupofcountriesselectsmeansealev elatadefinedtidegaugeoratagroupofgaugesforitsverti caldatumsurface .Itisrealized ,however,thatthelocalmeansealevelisusuallydepartedfromthegeoid ,whichshouldbetheidealdatumsurfaceforheight,ow…  相似文献   

4.
In geopotential space, the fundamental geodetic parameter W 0 defines the Gauss–Listing geoid which can be used to best represent the Earth’s mean sea level (MSL) and hence specifies a conventional zero height level to unify vertical datums employed by mapping agencies throughout the world. Further, W 0 cannot be considered invariant as the parameter varies temporally as a direct response to sea level change and mass redistributions. This study determines W 0 and its rate, dW 0/dt, by utilizing altimetric MSL models and an independent mean dynamic topography (MDT) model to define points on the geoid. W 0 and dW 0/dt are estimated by two approaches: (i) by means of a global gravity field model (GGM) and (ii) within normal gravity field space as the geopotential value of the best fitting reference ellipsoid. The study shows that uncertainty in W 0 is mainly influenced by MDT while the choice of methodology, GGM and MSL data coverage are not significant within reason. Our estimate W 0 =?62636854.2 ± 0.2 m2?s?2 at epoch 2005.0 differs by 1.8?m2s?2 from the International Astronomical Union reference value. This study shows that, at a sub-decadal time scale, the time variation dW 0/dt stems mainly from sea level change with negligible effect from gravity field variations. dW 0/dt =?(?2.70 ± 0.03)?×?10?2?m2?s?2?year?1, corresponding to a MSL rise of 2.9?mm?year?1, is evaluated from sea level change based on 16?years of TOPEX and Jason-1 data.  相似文献   

5.
Soil respiration (Rs) is of great importance to the global carbon balance. Remote sensing of Rs is challenging because of (1) the lack of long-term Rs data for model development and (2) limited knowledge of using satellite-based products to estimate Rs. Using 8-years (2002–2009) of continuous Rs measurements with nonsteady-state automated chamber systems at a Canadian boreal black spruce stand (SK-OBS), we found that Rs was strongly correlated with the product of the normalized difference vegetation index (NDVI) and the nighttime land surface temperature (LSTn) derived from Moderate Resolution Imaging Spectroradiometer (MODIS) imagery. The coefficients of the linear regression equation of this correlation between Rs and NDVI × LSTn could be further calibrated using the MODIS leaf area index (LAI) product, resulting in an algorithm that is driven solely by remote sensing observations. Modeled Rs closely tracked the seasonal patterns of measured Rs and explained 74–92% of the variance in Rs with a root mean square error (RMSE) less than 1.0 g C/m2/d. Further validation of the model from SK-OBS site at another two independent sites (SK-OA and SK-OJP, old aspen and old jack pine, respectively) showed that the algorithm can produce good estimates of Rs with an overall R2 of 0.78 (p < 0.001) for data of these two sites. Consequently, we mapped Rs of forest landscapes of Saskatchewan using entirely MODIS observations for 2003 and spatial and temporal patterns of Rs were well modeled. These results point to a strong relationship between the soil respiratory process and canopy photosynthesis as indicated from the greenness index (i.e., NDVI), thereby implying the potential of remote sensing data for detecting variations in Rs. A combination of both biological and environmental variables estimated from remote sensing in this analysis may be valuable in future investigations of spatial and temporal characteristics of Rs.  相似文献   

6.
Stability, which is significantly related to the loop parameters, is an important factor in the traditional GPS tracking loop design. Through the analysis of phase margin values in the discrete GPS PLL tacking loop, we are able to theoretically reveal the relationship between loop stability, equivalent noise bandwidth B n , predetection integration time T, and loop parameters. We calculate the theoretical limitations for B n T, that is, the product of equivalent noise bandwidth multiplied by predetection integration time, for second- and third-order phase-locked loop, respectively. The results are verified by actual data from GPS receivers.  相似文献   

7.
对流层顶与加权平均温度和可降水量(precipitable water vapor,PWV)之间存在很强的相关性,然而目前几乎没有学者讨论它对加权平均温度和PWV的影响。针对对流层顶对GNSS(global navigation satellite system)气象模型影响研究所存在的空白,基于已有学者提出的联系对流层顶和加权平均温度的公式,首次讨论中国区域对流层顶变化对加权平均温度的影响。在不损失公式精度的前提下,将该公式整理后得到对流层顶与加权平均温度的二次函数关系,分析了中国不同纬度区域的对流层顶对加权平均温度的影响,利用探空站观测数据得到了影响分布图。该图可以预测对流层顶对加权平均温度和PWV的影响。  相似文献   

8.
A 2×2 arc-minute resolution geoid model, CARIB97, has been computed covering the Caribbean Sea. The geoid undulations refer to the GRS-80 ellipsoid, centered at the ITRF94 (1996.0) origin. The geoid level is defined by adopting the gravity potential on the geoid as W 0=62 636 856.88 m2/s2 and a gravity-mass constant of GM=3.986 004 418×1014 m3/s2. The geoid model was computed by applying high-frequency corrections to the Earth Gravity Model 1996 global geopotential model in a remove-compute-restore procedure. The permanent tide system of CARIB97 is non-tidal. Comparison of CARIB97 geoid heights to 31 GPS/tidal (ITRF94/local) benchmarks shows an average offset (hHN) of 51 cm, with an Root Mean Square (RMS) of 62 cm about the average. This represents an improvement over the use of a global geoid model for the region. However, because the measured orthometric heights (H) refer to many differing tidal datums, these comparisons are biased by localized permanent ocean dynamic topography (PODT). Therefore, we interpret the 51 cm as partially an estimate of the average PODT in the vicinity of the 31 island benchmarks. On an island-by-island basis, CARIB97 now offers the ability to analyze local datum problems which were previously unrecognized due to a lack of high-resolution geoid information in the area. Received: 2 January 1998 / Accepted: 18 August 1998  相似文献   

9.
This paper presents a technique for ingesting ground- and space-based dual-frequency GPS observations into a semi-empirical global electron density model. The NeQuick-2 model is used as the basis for describing the global electron density distribution. This model is mainly driven by the F2 ionosphere layer parameters (i.e. the electron density, N m F2, and the height, h m F2 of the F2 peak), which, in the absence of directly measured values, are computed from the ITU-R database (ITU-R 1997). This database was established using observations collected from 1954 to 1958 by a network of around 150 ionospheric sounders with uneven global coverage. It allows computing monthly median values of N m F2 and h m F2 (intra-month variations are averaged), for low and high solar activity. For intermediate solar activity a linear interpolation must be performed. Ground-based GNSS observations from a global network of ~350 receivers are pre-processed in order to retrieve slant total electron content (sTEC) information, and space-based GPS observations (radio occultation data from the FORMOSAT-3/COSMIC constellation) are pre-processed to retrieve electron density (ED) information. Both, sTEC and ED are ingested into the NeQuick-2 model in order to adapt N m F2 and h m F2, and reduce simultaneously both, the observed minus computed sTEC and ED differences. The first experimental results presented in this paper suggest that the data ingestion technique is self consistent and able to reduce the observed minus computed sTEC and ED differences to ~25–30% of the values computed from the ITU-R database. Although sTEC and ED are both derived from GPS observations, independent algorithm and models are used to compute their values from ground-based GPS observations and space-based FORMOSAT-3/COSMIC radio occultations. This fact encourages us to pursue this research with the aim to improve the results presented here and assess their accuracy in a reliable way.  相似文献   

10.
An approach for estimating soil moisture is presented and tested by using surface-temperature-based soil evaporation transfer coefficient (ha), a coefficient recently proposed through the equation ha = (Ts − Ta)/(Tsd − Ta), where Ts, Tsd, and Ta are land surface temperature (LST), reference soil (dry soil without evaporation) surface temperature, and air temperature respectively. Our analysis and controllable experiment indicated that ha closely related to soil moisture, and therefore, a relationship between field soil moisture and ha could be developed for soil moisture estimation. Field experiments were carried out to test the relationship between ha and soil moisture. Time series Aqua-MODIS images were acquired between 11 Sep. 2006 and 1 Nov. 2007. Then, MODIS derived ha and simultaneous measured soil moisture for different soil depths were used to establish the relations between the two variables. Results showed that there was a logarithmic relationship between soil moisture and ha (P < 0.01). These logarithmic models were further validated by introducing another ground-truth data gathered from 46 meteorological stations in Hebei Province. Good agreement was observed between the measured and estimated soil moisture with RMSE of 0.0374 cm3/cm3 and 0.0503 cm3/cm3 for surface energy balance method at two soil depths (10 cm and 20 cm), with RMSE of 0.0467 cm3/cm3 and 0.0581 cm3/cm3 for maximum temperature method at two soil depths. For vegetated surfaces, the ratio of ha and NDVI suggested to be considered. The proposed approach has a great potential for soil moisture and drought evaluation by remote sensing.  相似文献   

11.
The gross primary production (GPP) at individual CO2 eddy covariance flux tower sites (GPPTower) in Dali (DL), Wenjiang (WJ) and Linzhi (LZ) around the southeastern Tibetan Plateau were determined by the net ecosystem exchange of CO2 (NEE) and ecosystem respiration (Re). The satellite remote sensing-VPM model estimates of GPP values (GPPMODIS) used the satellite-derived 8-day surface reflectance product (MOD09A1), including satellite-derived enhanced vegetation index (EVI) and land surface water index (LSWI). In this paper, we assembled a subset of flux tower data at these three sites to calibrate and test satellite-VPM model estimated GPPMODIS, and introduced the satellite data and site-level environmental factors to develop four new assimilation models. The new assimilation models’ estimates of GPP values were compared with GPPMODIS and GPPTower, and the final optimum model among the four assimilation models was determined and used to calibrate GPPMODIS. The results showed that GPPMODIS had similar temporal variations to the GPPTower, but GPPMODlS were commonly higher in absolute magnitude than GPPTower with relative error (RE) about 58.85%. While, the assimilation models’ estimates of GPP values (GPPMODEL) were much more closer to GPPTower with RE approximately 6.98%, indicating that the capacity of the simulation in the new assimilation model was greatly improved, the R2 and root mean square error (RMSE) of the new assimilation model were 0.57–4.90% higher and 0.74–2.47 g C m−2 s−1 lower than those of the GPPMODIS, respectively. The assimilation model was used to predicted GPP dynamics around the Tibetan Plateau and showed a reliable result compared with other researches. This study demonstrated the potential of the new assimilation model for estimating GPP around the Tibetan Plateau and the performances of site-level biophysical parameters in related to satellite-VPM model GPP.  相似文献   

12.
 A methodology for precise determination of the fundamental geodetic parameter w 0, the potential value of the Gauss–Listing geoid, as well as its time derivative 0, is presented. The method is based on: (1) ellipsoidal harmonic expansion of the external gravitational field of the Earth to degree/order 360/360 (130 321 coefficients; http://www.uni-stuttgard.de/gi/research/ index.html projects) with respect to the International Reference Ellipsoid WGD2000, at the GPS positioned stations; and (2) ellipsoidal free-air gravity reduction of degree/order 360/360, based on orthometric heights of the GPS-positioned stations. The method has been numerically tested for the data of three GPS campaigns of the Baltic Sea Level project (epochs 1990.8,1993.4 and 1997.4). New w 0 and 0 values (w 0=62 636 855.75 ± 0.21 m2/s2, 0=−0.0099±0.00079 m2/s2 per year, w 0/&γmacr;=6 379 781.502 m,0/&γmacr;=1.0 mm/year, and &γmacr;= −9.81802523 m2/s2) for the test region (Baltic Sea) were obtained. As by-products of the main study, the following were also determined: (1) the high-resolution sea surface topography map for the Baltic Sea; (2) the most accurate regional geoid amongst four different regional Gauss–Listing geoids currently proposed for the Baltic Sea; and (3) the difference between the national height datums of countries around the Baltic Sea. Received: 14 August 2000 / Accepted: 19 June 2001  相似文献   

13.
地基GPS气象学的关键技术是在于通过垂直方向上GPS信号的湿分量延迟值来确定大气可降水量-PWV,而这两个物理量间进行转换时需要用到一个关键性的参数——大气加权平均温度。本文首先讨论了估算加权平均温度的几种方法,然后利用成都地区2005年全年的大气探空数据,采用回归分析方法,建立了适合成都地区的大气加权平均温度模型,精度为±2.21K。最后将该公式应用到成都地区PWV的计算中,得到了非常理想的效果。  相似文献   

14.
Based upon a data set of 25 points of the Baltic Sea Level Project, second campaign 1993.4, which are close to mareographic stations, described by (1) GPS derived Cartesian coordinates in the World Geodetic Reference System 1984 and (2) orthometric heights in the Finnish Height Datum N60, epoch 1993.4, we have computed the primary geodetic parameter W 0(1993.4) for the epoch 1993.4 according to the following model. The Cartesian coordinates of the GPS stations have been converted into spheroidal coordinates. The gravity potential as the additive decomposition of the gravitational potential and the centrifugal potential has been computed for any GPS station in spheroidal coordinates, namely for a global spheroidal model of the gravitational potential field. For a global set of spheroidal harmonic coefficients a transformation of spherical harmonic coefficients into spheroidal harmonic coefficients has been implemented and applied to the global spherical model OSU 91A up to degree/order 360/360. The gravity potential with respect to a global spheroidal model of degree/order 360/360 has been finally transformed by means of the orthometric heights of the GPS stations with respect to the Finnish Height Datum N60, epoch 1993.4, in terms of the spheroidal “free-air” potential reduction in order to produce the spheroidal W 0(1993.4) value. As a mean of those 25 W 0(1993.4) data as well as a root mean square error estimation we computed W 0(1993.4)=(6 263 685.58 ± 0.36) kgal × m. Finally a comparison of different W 0 data with respect to a spherical harmonic global model and spheroidal harmonic global model of Somigliana-Pizetti type (level ellipsoid as a reference, degree/order 2/0) according to The Geodesist's Handbook 1992 has been made. Received: 7 November 1996 / Accepted: 27 March 1997  相似文献   

15.
In this study, sensible heat (H) calculation using remote sensing data over an alpine grass landscape is conducted from May to September 2010, and the calculation is validated using LAS (large aperture scintillometers) measurements. Data from two remote sensing sensors (FY3A-VIRR and TERRA-MODIS) are analysed. Remote sensing data, combined with the ground meteorological observations (pressure, temperature, wind speed, humidity) are fed into the SEBS (Surface Energy Balance System) model. Then the VIRR-derived sensible heat (VIRR_SEBS_H) and MODIS-derived sensible heat (MODIS_SEBS_H) are compared with the LAS-estimated H, which are obtained at the respective satellite overpass time. Furthermore, the similarities and differences between the VIRR_SEBS_H and MODIS_SEBS_H values are investigated. The results indicate that VIRR data quality is as good as MODIS data for the purpose of H estimation. The root mean square errors (rmse) of the VIRR_SEBS_H and MODIS_SEBS_H values are 45.1098 W/m2 (n = 64) and 58.4654 W/m2 (n = 71), respectively. The monthly means of the MODIS_SEBS_H are marginally higher than those of VIRR_SEBS_H because the satellite overpass time of the TERRA satellite lags by 25 min to that of the FT3A satellite. Relative evaporation (EFr), which is more time-independent, shows a higher agreement between MODIS and VIRR. Many common features are shared by the VIRR_SEBS_H and the MODIS_SEBS_H, which can be attributed to the SEBS model performance. In May–June, H is over-estimated with more fluctuations and larger rmse, whereas in July–September, H is under-estimated with fewer fluctuations and smaller rmse. Sensitivity analysis shows that potential temperature gradient (delta_T) plays a dominant role in determining the magnitude and fluctuation of H. The largest rmse and over-estimation in H occur in June, which could most likely be attributed to high delta_T, high wind speed, and the complicated thermodynamic state during the transitional period when bare land transforms to dense vegetation cover.  相似文献   

16.
大气加权平均温度的准确获取对高精度的GPS水汽反演至关重要。文中基于线性回归理论,在分析加权平均温度与地面温度间相关性的基础上,采用一元线性拟合的方法,建立大气加权平均温度经验模型。最后,采用香港地区2006-2015年无线电探空资料对经验模型进行验证。实验结果表明,文中模型计算加权平均温度的整体均方根误差为2.356 K,较Bevis模型精度提高了41.94%,且季节变化对加权平均温度计算的影响并不明显;对于GPS水汽反演,采用本文经验模型反演水汽的均方根误差为1.807 mm,平均偏差为1.362 mm,能够满足GPS可降水量反演的精度,且优于Bevis模型。   相似文献   

17.
Surface soil moisture (SSM) is a critical variable for understanding the energy and water exchange between the land and atmosphere. A multi-linear model was recently developed to determine SSM using ellipse variables, namely, the center horizontal coordinate (x0), center vertical coordinate (y0), semi-major axis (a) and rotation angle (θ), derived from the elliptical relationship between diurnal cycles of land surface temperature (LST) and net surface shortwave radiation (NSSR). However, the multi-linear model has a major disadvantage. The model coefficients are calculated based on simulated data produced by a land surface model simulation that requires sufficient meteorological measurements. This study aims to determine the model coefficients directly using limited meteorological parameters rather than via the complicated simulation process, decreasing the dependence of the model coefficients on meteorological measurements. With the simulated data, a practical algorithm was developed to estimate SSM based on combined optical and thermal infrared data. The results suggest that the proposed approach can be used to determine the coefficients associated with all ellipse variables based on historical meteorological records, whereas the constant term varies daily and can only be determined using the daily maximum solar radiation in a prediction model. Simulated results from three FLUXNET sites over 30 cloud-free days revealed an average root mean square error (RMSE) of 0.042 m3/m3 when historical meteorological records were used to synchronously determine the model coefficients. In addition, estimated SSM values exhibited generally moderate accuracies (coefficient of determination R2 = 0.395, RMSE = 0.061 m3/m3) compared to SSM measurements at the Yucheng Comprehensive Experimental Station.  相似文献   

18.
The estimates of total zenith delay are derived using Bernese GPS Software V4. 2 based on GPS data every 30 s from the first measurement experiment of a ground-based GPS network in Chengdu Plain of Southwest China during the period from July to September 2004. Then the estimates of 0.5 hourly precipitable water vapor (PWV) derived from global positioning system (GPS) are obtained using meteorological data from automatic weather stations (AWS). The comparison of PWV derived from GPS and those from radiosonde observations is given for the Chengdu station, with RMS (root mean square) differences of 3.09m. The consis- tency of precipitable water vapor derived from GPS to those from radiosonde is good. It is concluded that Bevis’ empirical formula for estimating the weighted atmospheric mean temperature can be applicable in Chengdu area because the relationship of GPS PWV with Bevis’ formula and GPS PWV with radiosonde method shows a high correlation. The result of this GPS measurement experiment is helpful both for accumu- lating the study of precipitable water vapor derived from GPS in Chengdu areas located at the eastern side of the Tibetan Plateau and for studying spatial-temporal variations of regional atmospheric water vapor through many disciplines cooperatively.  相似文献   

19.
The estimates of total zenith delay are derived using Bernese GPS Software V4. 2 based on GPS data every 30 s from the first measurement experiment of a ground-based GPS network in Chengdu Plain of Southwest China during the period from July to September 2004. Then the estimates of 0.5 hourly precipitable water vapor (PWV) derived from global positioning system (GPS) are obtained using meteorological data from automatic weather stations (AWS). The comparison of PWV derived from GPS and those from radiosonde observations is given for the Chengdu station, with RMS (root mean square) differences of 3.09m. The consistency of precipitable water vapor derived from GPS to those from radiosonde is good. It is concluded that Bevis’ empirical formula for estimating the weighted atmospheric mean temperature can be applicable in Chengdu area because the relationship of GPS PWV with Bevis’ formula and GPS PWV with radiosonde method shows a high correlation. The result of this GPS measurement experiment is helpful both for accumulating the study of precipitable water vapor derived from GPS in Chengdu areas located at the eastern side of the Tibetan Plateau and for studying spatial-temporal variations of regional atmospheric water vapor through many disciplines cooperatively.  相似文献   

20.
SPOT satellites have been imaging Earth's surface since SPOT 1 was launched in 1986. It is argued that absolute atmospheric correction is a prerequisite for quantitative remote sensing. Areas where land cover changes are occurring rapidly are also often areas most lacking in situ data which would allow full use of radiative transfer models for reflectance factor retrieval (RFR). Consequently, this study details the proposed historical empirical line method (HELM) for RFR from multi-temporal SPOT imagery. HELM is designed for use in landscape level studies in circumstances where no detailed overpass concurrent atmospheric or meteorological data are available, but where there is field access to the research site(s) and a goniometer or spectrometer is available. SPOT data are complicated by the ±27° off-nadir cross track viewing. Calibration to nadir only surface reflectance factor (ρs) is denoted as HELM-1, whilst calibration to ρs modelling imagery illumination and view geometries is termed HELM-2. Comparisons of field measured ρs with those derived from HELM corrected SPOT imagery, covering Helsinki, Finland, and Taita Hills, Kenya, indicated HELM-1 RFR absolute accuracy was ±0.02ρs in the visible and near infrared (VIS/NIR) bands and ±0.03ρs in the shortwave infrared (SWIR), whilst HELM-2 performance was ±0.03ρs in the VIS/NIR and ±0.04ρs in the SWIR. This represented band specific relative errors of 10–15%. HELM-1 and HELM-2 RFR were significantly better than at-satellite reflectance (ρSAT), indicating HELM was effective in reducing atmospheric effects. However, neither HELM approach reduced variability in mean ρs between multi-temporal images, compared to ρSAT. HELM-1 calibration error is dependent on surface characteristics and scene illumination and view geometry. Based on multiangular ρs measurements of vegetation-free ground targets, calibration error was negligible in the forward scattering direction, even at maximum off-nadir view. However, error exceeds 0.02ρs where off-nadir viewing was ≥20° in the backscattering direction within ±55° azimuth of the principal plane. Overall, HELM-1 results were commensurate with an identified VIS/NIR 0.02ρs accuracy benchmark. HELM thus increases applicability of SPOT data to quantitative remote sensing studies.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号