首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 215 毫秒
1.
Methods to derive wind speed and sea state by simple empirical models from synthetic aperture radar (SAR) data are presented and applied for use in high-resolution numerical modelling for coastal applications. The new radar satellite, TerraSAR-X (TS-X), images the surface of the sea with a high resolution up to 1 m. Therefore, not only wind information and integrated sea state parameters but also individual ocean waves with wavelengths down to 30 m are detectable. Two-dimensional information on the ocean surface retrieved using TS-X data is validated for different oceanographic applications: derivation of finely resolved wind fields (XMOD algorithm) and integrated sea state parameters (XWAVE algorithm). Both algorithms are capable of taking into account fine-scale effects in coastal areas. Wind and sea state information retrieved from SAR data are applied as the input for a wave numerical spectral model (wind forcing and boundary condition) running at a fine spatial horizontal resolution of 100 m. Results are compared to collocated buoy measurements. Studies are carried out for varying wind speeds and comparisons against wave height, simulated using original TS-X-derived wind data, showing the sensitivity of waves to local wind variation and thus the importance of local wind effects on wave behaviour in coastal areas. Examples for the German Bight (North Sea) are shown. The TS-X satellite scenes render well-developed ocean wave patterns of developed swell at the sea surface. Refraction of individual long swell waves at a water depth shallower than about 70 m, caused by the influence of underwater topography in coastal areas, is imaged on the radar scenes. A technique is developed for tracking wave rays depending on changes in swell wavelength and direction. We estimate the wave energy flux along wave tracks from deep water to the coastline based on SAR information: wave height and wavelength are derived from TS-X data.  相似文献   

2.
The wave pattern generated by a moving ship is formed by two dominant features: the turbulent wake and a 'V'-shaped pattern trailing the ship, consisting of the two Kelvin arms. In this paper we investigate the radar imaging mechanism of Kelvin arms, which are formed by the cusp waves. A composite surface model for the radar backscattering at the ocean surface is used. The radar signatures of Kelvin arms can be attributed to tilt and hydrodynamic modulation of Bragg waves by the cusp waves. The proposed model allows the computation of the normalized radar backscattering cross-section (NRCS) as a function of radar frequency, polarization, incidence angle, wind speed and direction, and wavelength, direction, and slope of the cusp waves. By using this imaging model, radar signatures of cusp waves are calculated for several spaceborne Synthetic Aperture Radars (SARs): (1) the SEASAT L-band HH-polarized SAR, (2) the ERS-1/-2 VV-polarized SAR, (3) the RADARSAT C-band HH-polarized SAR, and (4) the X-, C- and L-band multipolarization SARs of the Space Radar Laboratory flown on the space shuttle during the SIRC/X-SAR mission in 1994. The results of the simulations are compared with SEASAT and SIR-C/X-SAR imagery of ship wake patterns. It is shown that the dependence of the observed radar signatures of Kelvin arms on radar look direction is consistent with the proposed imaging theory and that the measured relative mean NRCS values induced by Kelvin arms can be fairly well reproduced by the proposed model. The simulations indicate that ship wake signatures should be more clearly visible on SEASAT L-band SAR than on ERS-1/-2 or RADARSAT C-band SAR images. The radar signatures of Kelvin arms are strongest at low wind speeds and are not very sensitive to wind direction.  相似文献   

3.
Abstract

A new era of remote sensing for coastal and oceanographic monitoring was born on 26 June 1978 with the launch of Seasat. Duck-X was a 2 month experiment conducted during August to October 1978 off the east coast of the U.S.A. for the validation of the Seasat synthetic aperture radar (SAR), During this field experiment, various oceanographic phenomena were monitored. Ground truth observations of these phenomena have been correlated with Seasat SAR imagery. The ground truth sensors included airborne photographic and radar imagery, meteorological satellite imagery, land based radars, and conventional wave gauges. This paper focuses on ocean surface waves, ocean currents and coastal inlet discharge

Specifically, the direction and length of the principal ocean wave trains are compared for the periods of Seasat overflight of the Duck-X area. During these overflights significant wave heights were 1.5 m and less and the maximum wave period was 14 s. The current correlations concentrate on the western boundary of the Gulf Stream and its associated eddy structure. Inlet outflow is shown for inlets on the east coast of the U.S.A

This ground truth study has indicated that the SAR imagery contains an unanticipated abundance of information on a variety of oceanographic and coastal phenomena.  相似文献   

4.
Scatterometer surface wind speed and direction observations in combination with radiometer wind speeds allow to generate surface wind analyses with high space and time resolutions over global as well as at regional scales. Regarding scatterometer sampling schemes and physics, the resulting surface wind analyses suffer from lack of accuracy in areas near coasts. The use of the synthetic aperture radar (SAR) onboard the Sentinel-1A satellite attempts to address the enhancement of surface wind analyses issues. In this study, SAR wind speeds and directions retrieved from backscatter coefficients acquired in interferometric wide (IW) swath mode are used. Their accuracy is determined through comprehensive comparisons with moored buoy wind measurements. SAR and buoy winds agree well at offshore and nearshore locations. The statistics characterizing the comparison of SAR and buoy wind speeds and directions are of the same order as those obtained from scatterometer (Advanced SCATterometer (ASCAT) and RapidScat) and buoy wind comparisons. The main discrepancy between SAR and buoy data are found for high wind speeds. SAR wind speeds exceeding 10 m s–1 tend to be underestimated. A similar conclusion is drawn from SAR and scatterometer wind speed comparisons. It is based on the underestimation of SAR backscatter coefficient (σ°) with respect to σ° estimated from scatterometer winds and the geophysical model function (GMF) named CMOD-IFR2 (Ifremer C band MODel). New SAR wind speeds are retrieved using CMOD-IFR2. The corrected SAR retrievals allow better determination of the spatial characteristics of surface wind speeds and of the related wind components in near-coast areas. They are used for enhancing the determination of the spatial structure function required for the estimation of wind fields gridded in space and time at the regional scale. The resulting wind fields are only determined from scatterometer wind observations in combination with radiometer retrievals. Their qualities are determined through comparisons with SAR wind speeds and directions, and through their application for determination of wind power off Brittany coasts.  相似文献   

5.
This paper discusses the occurrences of von Karman vortices in the lee of Cheju Island. In this area an isolated mountain often disturbs airflow over the Korea South Sea. The vortex street often occurs with a strong northerly flow with persistent patterns of cumulus and stratiform clouds. The generation of vortices and cloud streaks including enhancement of precipitation in the lee side observed with NOAA satellites are discussed in this paper. Observations of vortex and cloud formation by the Halla Mountain up to 800 km are recorded. Ten cases shown here are classic examples of lee disturbances in meso‐ and synoptic‐scales which are of interest for atmospheric dynamics research and tutorial aid.  相似文献   

6.
Images of rain events over the ocean acquired by a multi‐frequency/multi‐polarization Synthetic Aperture Radar (SAR) show different radar contrasts at different frequencies and polarizations. In order to better understand these effects, field and laboratory experiments were performed at different rain rates and wind speeds with scatterometers working at different radar frequencies, polarizations, and incidence angles. Our results show that the dominant scattering mechanism on a rain‐roughened water surface, observed at VV polarization, at all incidence angles is Bragg scattering from ring waves. At HH polarization the radar backscatter is caused by both ring waves and non‐propagating splash products, with the dominating effect depending on incidence angle. The reduction and enhancement of the surface roughness by ring waves and sub‐surface phenomena, respectively, result in a transition wavenumber between reduction of the radar backscattering and its enhancement of about 100 rad m?1. We assume that this transition wavenumber depends on the drop‐size distribution of the rain. Taking into consideration the different dependencies of the radar backscatter at different frequencies and polarizations on rain rate, we suggest a method to estimate rain rates by calculating the ratio of the radar cross‐sections at L band, VV polarization and at C band, HV polarization. Provided an availability of SAR data at the respective frequency–polarization combinations, this method allows for investigating the nature of small‐scale (convective) rain events over the ocean.  相似文献   

7.
Several algorithms have been proposed to retrieve near-surface wind fields from C-band synthetic aperture radar (SAR) images acquired over the ocean. They mainly differ in the way they retrieve the wind direction. Conventionally, the wind direction is taken from atmospheric models or is extracted from the linear features sometimes visible in SAR images. Recently, a new wind retrieval algorithm has been proposed, which also includes the Doppler shift induced by motions of the sea surface. In this article, we apply three wind retrieval algorithms, including the one using Doppler information, to three complex wind events encountered over the Black Sea and compare the SAR-derived wind fields with model wind fields calculated using the high-resolution weather research and forecasting (WRF) model. It is shown that the new algorithm is very efficient in resolving the 180° ambiguity in the wind direction, which is often a problem in the streak-based wind retrieval algorithms. However, the Doppler-based algorithm only yields good results for wind directions that have a significant component in the look direction of the SAR antenna. Furthermore, it is dependent on good separation of the contributions to the Doppler shift induced by surface currents and wind-related effects (wind drift and wind-sea components of the ocean wave spectrum). We conclude that an optimum wind retrieval algorithm should consist of a combination of the algorithms based on linear features and Doppler information.  相似文献   

8.
In this review, recent studies on the observations of typhoon eyes by images acquired by multiple sensors, including synthetic aperture radar (SAR), and infrared (IR) radiometer, are first summarized. Large horizontal distances between typhoon eyes on the ocean surface by SAR and those on the cloud top by IR sensors have been demonstrated; these have previously been ignored but should not be ignored in typhoon forecasts and numerical simulations. Then, based on nine published typhoon cases, the horizontal shifts and vertical tilt angles from the cloud-top typhoon eye locations by IR sensors on board the Feng-Yun 2 (FY-2) and Multi Functional Transport Satellite (MTSAT) to those at sea surface by SAR are further estimated. This shift difference between different sensors raises an issue on project distortion and navigation system errors for FY-2 and MTSAT satellites, which are of concern to both space agencies and data users. Finally, issues for current ongoing study and future research related to typhoon eyes are discussed, including rainband tracking between sensors for local wind speeds.  相似文献   

9.
The vertical and horizontal distributions of the cloud types across different seasons and over the contiguous USA and surrounding areas are studied. The study is performed by collecting two years (2007 and 2008) of data from the CloudSat 2B-CLDCLASS product that uses effective radar reflectivity factor Ze, the presence of precipitation and ancillary data such as surface topography and the model-predicted temperature profile to classify clouds into seven distinct types. Considerable seasonal variations of the horizontal distribution of the cloud-type fractions are observed in the study area among different seasons and for both daytime and night-time CloudSat observations. It was found that during spring and summer, deep convective (Dc) clouds are observed much more frequently during night‐time than during daytime over both the land and ocean. For the studied area and during daytime, low clouds were more frequent (up to ?50%) over the land and less frequent over the ocean compared with night-time observations. Analysis of the vertical distribution of cloud layers reveals that the fraction of cloudy scenes with two or more distinct cloud layers is the highest (up to 30%) over the northwest corner of the USA and the southwest corner of Canada and the nearby oceans. The southwest corner of the USA and the nearby east Pacific Ocean appeared to have the lowest fraction (<0.05%) of cloudy scenes with two or more distinct cloud layers. Over the land, approximately 18% of the total cloudy scenes are classified as two-layer clouds, whereas over the ocean, two-layer clouds are less frequent and range from 13% to 17% with a stronger seasonal dependency. Only about 2–3% of the total cloudy scenes are classified as multilayer clouds, with three or more distinct layers over both the land and ocean. The vertical distribution of cloud-top heights over both the land and ocean shows two distinct peaks. Over the land, the lower peak, at around 2 km, is almost independent of season, whereas the higher peak is seasonally dependent and varies between ?8 km (during winter) and ?11 km (during summer). Over the ocean, the lower peak is also observed near 2 km (or less), whereas the higher peak ranges approximately from 11 km (during winter) to 12 km (during summer).  相似文献   

10.
海洋锋区剪切流、辐聚流与海洋表面波之间存在相互作用,将这种相互作用处理为对海浪谱的扰动,来分析其对海浪谱密度和谱梯度的影响;根据两尺度模型,分析了海面小尺度波(厘米级)和大尺度波(海浪)与雷达后向散射系数的关系,从而说明了海洋锋的SAR成像原理。在海面SAR图像中,海洋锋的尺度比海浪的尺度大2~3个数量级,可以通过二维空间谱分析,将海浪信息的主要部分滤除,再利用数字图像处理技术提取海洋锋的特征信息,由此形成了一套用海面SAR图像提取海洋锋特征参数的信息处理方法。  相似文献   

11.

The results of two polarization airborne radar imagery tests of the ocean surface obtained during the JUSREX'92 experiment are presented. It is shown that the traditional composite surface model with small-scale 'Bragg' waves superposed over larger gravity waves can not explain either the contrasts of internal wave surface manifestations in conditions of a stable atmospheric boundary layer at low grazing angles (LGA), or the apparent difference between the images obtained at different polarizations in unstable atmospheric conditions. We attribute this discrepancy to the presence of mesoscale steep waves, which produce non-resonant scattering and make different relative contributions to the total cross sections for the two polarizations. The possibility of distinguishing between surface manifestations of atmospheric and oceanic origin is also discussed.  相似文献   

12.
ABSTRACT

Improving the understanding of cloud–radiation–monsoon interactions is difficult due to the limited knowledge regarding the impacts of vertical cloud radiative forcing on monsoon circulation. Here, we focus on the annual cycle of the vertical structure of cloud-induced radiative heating (CRH) to evaluate further their impacts on the East Asian monsoon circulation (100°–140° E, 20°–45° N) derived from satellite observations and reanalysis datasets. Entire troposphere and lower stratosphere are heated by vertical CRH, with the peak reaching 1 K day?1 at the mid-level troposphere (4–10 km) during summer. Although radiative warming occurs below 3 km from the prevailing stratocumulus, widespread weak radiative cooling (approximately ?0.2 K day?1) occurs at a wide vertical range above 3 km during winter. Consequently, the wind vector variations resulting from vertical CRH highly coincide with the monsoon circulation, leading to the increase in wind speeds by 1.8 and 0.5 m s?1 during summer and winter, respectively, while a weakly negative influence (about 0.3 m s?1) occurs at the low-level troposphere below 3 km during winter. Although high clouds, stratiform clouds, and stratocumulus dominate these wind vector variations, deep convective clouds generate the strongest updraft (up to 7 m s?1) amongst all cloud categories despite their low occurrence frequency. Results highlight the important enhancement of vertical CRH to East Asian monsoon circulation by perturbing the vertical structure of heating rate.  相似文献   

13.
SAR(Synthetic Aperture Radar,合成孔径雷达)作为一种现代高空间分辨率成像侧视雷达,对地球表面海洋所成的图像中蕴含了极为丰富的中尺度及亚中尺度海洋大气边界层的信息,因此对边界层气象学研究有着非常重要的意义。但是,使用SAR研究海气边界层这一涵盖微波遥感、气象学及海洋学等学科的科学前沿课题在国内却少有文献报道。在此背景下,首先介绍了SAR反演海洋大气边界层的研究概况,回顾了SAR反演海气边界层参数的原理和方法。然后以2002年5月7日当地时间10时53分ERS-2卫星获取的香港地区(22.097°N,E 114.300°E)SAR海洋图像为例,进行了反演风向风速的初步试验,最终获得了较高精度的风矢量。具体过程如下:先对SAR图像进行预处理,包括ADC(Analog Digital Converter,模数转换器)补偿、精确校准及斑点滤波等过程;然后利用经典的谱分析方法求得具有180°模糊度的风向,再用香港天文台气象浮标实测资料消除这一不确定性得到了真实的相对风向;紧接着利用CMOD4地球物理模式函数计算得到了海面上10 m高的风速。与气象浮标站所记录的平均风速和风向比较,两个20 km×20 km大小的试验区域求得的风向误差分别为23.71°和7.00°,平均风速误差分别为0.18 m/s和-0.12 m/s。结果表明,如果对SAR预先进行严格的预处理,结合经典的谱分析方法和CMOD4模型,即可获取高精度的风矢量。这一结果为今后海洋大气边界层的研究奠定了良好的基础。  相似文献   

14.
This paper outlines the development of a multi‐satellite precipitation estimation methodology that draws on techniques from machine learning and morphology to produce high‐resolution, short‐duration rainfall estimates in an automated fashion. First, cloud systems are identified from geostationary infrared imagery using morphology based watershed segmentation algorithm. Second, a novel pattern recognition technique, growing hierarchical self‐organizing map (GHSOM), is used to classify clouds into a number of clusters with hierarchical architecture. Finally, each cloud cluster is associated with co‐registered passive microwave rainfall observations through a cumulative histogram matching approach. The network was initially trained using remotely sensed geostationary infrared satellite imagery and hourly ground‐radar data in lieu of a dense constellation of polar‐orbiting spacecraft such as the proposed global precipitation measurement (GPM) mission. Ground‐radar and gauge rainfall measurements were used to evaluate this technique for both warm (June 2004) and cold seasons (December 2004–February 2005) at various temporal (daily and monthly) and spatial (0.04° and 0.25°) scales. Significant improvements of estimation accuracy are found classifying the clouds into hierarchical sub‐layers rather than a single layer. Furthermore, 2‐year (2003–2004) satellite rainfall estimates generated by the current algorithm were compared with gauge‐corrected Stage IV radar rainfall at various time scales over continental United States. This study demonstrates the usefulness of the watershed segmentation and the GHSOM in satellite‐based rainfall estimations.  相似文献   

15.
Under strong ocean surface wind conditions, the normalized radar cross section of synthetic aperture radar (SAR) is dampened at certain incident angles, compared with the signals under moderate winds. This causes a wind speed ambiguity problem in wind speed retrievals from SAR, because two solutions may exist for each backscattered signal. This study shows that the problem is ubiquitous in the images acquired by operational space‐borne SAR sensors. Moreover, the problem is more severe for the near range and range travelling winds. To remove this ambiguity, a method was developed based on characteristics of the hurricane wind structure. A SAR image of Hurricane Rita (2005) was analysed to demonstrate the wind speed ambiguity problem and the method to improve the wind speed retrievals. Our conclusions suggest that a speed ambiguity removal algorithm must be used for wind retrievals from SAR in intense storms and hurricanes.  相似文献   

16.
Alternative dark–bright patterns on two ENVISAT Advanced Synthetic Aperture Radar (ASAR) images of the east coast of the Korean Peninsula acquired on 18 and 19 May 2004 are interpreted as atmospheric gravity waves (AGWs) on the basis of simultaneous multi‐satellite observations and atmospheric gravity wave theory. The AGWs appeared in the form of a wave packet containing several waves located between 50 and 200 km offshore. The wavelengths were ranging from 13 to 20 km. The lengths of AGW crests were from 20 to 150 km. An NOAA‐17 pass was received about 30 min after the ASAR pass. Its channel 4 infrared (IR) image clearly shows wave‐like moisture patterns. However, the sea surface temperature (SST) image derived after applying nonlinear calibration and split‐window atmospheric correction shows no wave patterns. A daytime true‐colour MODIS image taken about 14 h later still shows a cloud band of same AGWs, indicating the lifespan of the standing AGWs can be over half a day. Although oceanic internal waves (IWs) may also cause similar wave patterns imaged by spaceborne SAR as they modulate the ocean surface roughness, we provide evidence to eliminate this possibility in this case. The characteristics of satellite observed AGWs are in good agreement with these simulated by a linear coastal AGW model.  相似文献   

17.
SCIAMACHY (SCanning Imaging Absorption spectroMeter for Atmospheric CartograpHY) is a passive remote sensing spectrometer observing backscattered radiation from the atmosphere and the Earth's surface, in the wavelength range between 240 and 2380 nm. The instrument is onboard ENVironmental SATellite (ENVISAT) which was launched on 1 March 2002. The Medium Resolution Imaging Spectrometer (MERIS) is also one of the 10 instruments onboard the ENVISAT satellite. MERIS is a 68.5° field-of-view nadir-pointing imaging spectrometer which measures the solar radiation reflected by the Earth in 15 spectral bands (visible and near-infrared). It obtains a global coverage of the Earth in three days. Its main objective is to measure sea colour and quantify ocean chlorophyll content and sediment, thus providing information on the ocean carbon cycle and thermal regime. It is also used to derive the cloud top height, aerosol and cloud optical thickness, and water vapour column. The ground spatial resolution of the instrument is 260 m × 290 m. This paper is aimed at determining the cloud fraction in SCIAMACHY pixels (typically, 30 km × 60 km ground scenes) using MERIS observations and number of thresholds for MERIS top-of-atmosphere reflectances and their ratios. Thresholds utilize the fact that clouds are bright white objects having similar reflectances in the blue and red. The MERIS cloud fraction has been derived for a number of SCIAMACHY states with area of 916 km × 400 km. The results are compared with correspondent cloud fractions obtained using SCIAMACHY polarization measurement devices (PMDs). Large differences are found between cloud fractions derived using SCIAMACHY and MERIS measurements. It is recommended to use highly spatially resolved MERIS observations instead of SCIAMACHY PMD measurements to retrieve cloud fractions in SCIAMACHY pixels. The improvements advised will enhance SCIAMACHY trace gas and cloud retrievals in the presence of broken cloud fields.  相似文献   

18.
In this article, the polarization ratio (PR) of TerraSAR-X (TS-X) vertical–vertical (VV) and horizontal–horizontal (HH) polarization data acquired over the ocean is investigated. Similar to the PR of C-band synthetic aperture radar (SAR), the PR of X-band SAR data also shows significant dependence on incidence angle. The normalized radar cross-section (NRCS) in VV polarization data is generally larger than that in HH polarization for incidence angles above 23°. Based on the analysis, two PR models proposed for C-band SAR were retuned using TS-X dual-polarization data. A new PR model, called X-PR hereafter, is proposed as well to convert the NRCS of TS-X in HH polarization to that in VV polarization. By using the developed geophysical model functions of XMOD1 and XMOD2 and the tuned PR models, the sea surface field is retrieved from the TS-X data in HH polarization. The comparisons with in situ buoy measurements show that the combination of XMOD2 and X-PR models yields a good retrieval with a root mean square error (RMSE) of 2.03 m s–1 and scatter index (SI) of 22.4%. A further comparison with a high-resolution analysis wind model in the North Sea is also presented, which shows better agreement with RMSE of 1.76 m s–1 and SI of 20.3%. We also find that the difference between the fitting of the X-PR model and the PR derived from TS-X dual-polarization data is close to a constant. By adding the constant to the X-PR model, the accuracy of HH polarization sea surface wind speed is further improved with the bias reduced by 0.3 m s–1. A case acquired at the offshore wind farm in the East China Sea further demonstrates that the improvement tends to be more effective for incidence angles above 40°.  相似文献   

19.
In this article, a method for the detection of wave field parameters from synthetic aperture radar (SAR) imagery in the fetch-limited Baltic Sea is presented. Over the Baltic Sea region, common southwest (SW) and west (W) winds induce steep waves with shorter wavelengths compared with ocean waves. Thus, with the use of previous SAR sensors (e.g. ENVISAT/ASAR), it was not possible to detect individual waves and retrieve image wave number spectra. Since the year 2007, when TerraSAR-X (TS-X) reached its orbit, high spatial resolution data is available for measuring the sea-state parameters: the individual waves up to 30 m wavelength and their refraction can be distinguished. The main objective of this work was to demonstrate the capability of detecting wave field parameter from (TS-X) imagery in the Baltic Sea. The wave field parameters obtained from the SAR imagery were compared with in situ measurements and the Simulating WAves Nearshore (SWAN) wave model. The comparison of SAR-based wave field information with buoy measurements showed high agreement in case of wave propagation direction (r = 0.95) and wavelength (r = 0.83). A significant correlation is also seen between SWAN- and SAR-derived wave propagation direction (r = 0.87) and wavelengths (r = 0.91). With the case studies, it is shown that SAR data enables one to detect land shadow effects and small-scale wave field variations in the coastal zone. It was shown that SAR data is also valuable for improving and interpreting the wave model results. In consequence of common slanting fetch cases over the Baltic Sea region, it was demonstrated that the peak wave directions differ from the mean wind directions up to 43°.  相似文献   

20.
合成孔径雷达海面风场反演研究进展   总被引:1,自引:0,他引:1       下载免费PDF全文
合成孔径雷达高分辨率、全天候的观测能力使其成为大面积获取高分辨海面风场的重要手段,这有助于我们理解各种海洋现象的物理过程,尤其是在海岸带,合成孔径雷达反演风场更具优势。以合成孔径雷达风向确定方式为主线,详细阐述了国内外学者在合成孔径雷达风场反演方面的研究进展,主要包括:利用影像线性纹理特征进行风场反演、基于外部初始风向的风场反演方法以及利用影像本身所包含的其它信息进行的风场反演工作,诸如距离向入射角差异、多极化后向散射等,并对合成孔径雷达风场反演的发展给出自己的观点。  相似文献   

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

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