首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 156 毫秒
1.
A new algorithm for automatic detection of prominences on the solar limb in 304 Å EUV images is presented, and results of its application to SOHO/EIT data discussed. The detection is based on the method of moments combined with a classifier analysis aimed at discriminating between limb prominences, active regions, and the quiet corona. This classifier analysis is based on a Support Vector Machine (SVM). Using a set of 12 moments of the radial intensity profiles, the algorithm performs well in discriminating between the above three categories of limb structures, with a misclassification rate of 7%. Pixels detected as belonging to a prominence are then used as the starting point to reconstruct the whole prominence by morphological image-processing techniques. It is planned that a catalogue of limb prominences identified in SOHO and STEREO data using this method will be made publicly available to the scientific community.  相似文献   

2.
The SPIRIT complex onboard the CORONAS-F satellite has routinely imaged the Sun in the 171, 175, 195, 284, and 304 Å spectral bands since August 2001. The complex incorporates two telescopes. The Ritchey-Chretien telescope operates in the 171, 195, 284, and 304 Å bands and has an objective similar to that of the SOHO/EIT instrument. The Herschel telescope obtains solar images synchronously in the 175 and 304 Å bands with two multilayer-coated parabolic mirrors. The SPIRIT program includes synoptic observations, studies of the dynamics of various structures on the solar disk and in the corona up to 5 solar radii, and coordinated observations with other spaceborne and ground-based telescopes. In particular, in the period 2002–2003, synoptic observations with the SPIRIT Ritchey-Chretien telescope were coordinated with regular 6-hour SOHO/EIT observations. Since June 2003, when EIT data were temporarily absent (SOHO keyholes), the SPIRIT telescope has performed synoptic observations at a wavelength of 175 A. These data were used by the Solar Influence Data Analysis Center (SIDC) at the Royal Observatory of Belgium for an early space weather forecast. We analyze the photometric and spectral parameters of the SPIRIT and EIT instruments and compare the integrated (over the solar disk) EUV fluxes using solar images obtained with these instruments during the CORONAS-F flight from August 2001 through December 2003.  相似文献   

3.
The Solar Dynamics Observatory provides multiwavelength imagery from extreme ultraviolet (EUV) to visible light as well as magnetic-field measurements. These data enable us to study the nature of solar activity in different regions of the Sun, from the interior to the corona. For solar-cycle studies, synoptic maps provide a useful way to represent global activity and evolution by extracting a central meridian band from sequences of full-disk images over a full solar Carrington rotation (≈?27.3 days). We present the global evolution during Solar Cycle 24 from 20 May 2010 to 31 August 2013 (CR?2097?–?CR?2140), using synoptic maps constructed from full-disk, line-of-sight magnetic-field imagery and EUV imagery (171 Å, 193 Å, 211 Å, 304 Å, and 335 Å). The synoptic maps have a resolution of 0.1 degree in longitude and steps of 0.001 in sine of latitude. We studied the axisymmetric and non-axisymmetric structures of solar activity using these synoptic maps. To visualize the axisymmetric development of Cycle 24, we generated time–latitude (also called butterfly) images of the solar cycle in all of the wavelengths, by averaging each synoptic map over all longitudes, thus compressing it to a single vertical strip, and then assembling these strips in time order. From these time–latitude images we observe that during the ascending phase of Cycle 24 there is a very good relationship between the integrated magnetic flux and the EUV intensity inside the zone of sunspot activities. We observe a North–South asymmetry of the EUV intensity in high-latitudes. The North–South asymmetry of the emerging magnetic flux developed and resulted in a consequential asymmetry in the timing of the polar magnetic-field reversals.  相似文献   

4.
Coronal holes (CHs) are regions of open magnetic field lines in the solar corona and the source of the fast solar wind. Understanding the evolution of coronal holes is critical for solar magnetism as well as for accurate space weather forecasts. We study the extreme ultraviolet (EUV) synoptic maps at three wavelengths (195 Å/193 Å, 171 Å and 304 Å) measured by the Solar and Heliospheric Observatory/Extreme Ultraviolet Imaging Telescope (SOHO/EIT) and the Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA) instruments. The two datasets are first homogenized by scaling the SDO/AIA data to the SOHO/EIT level by means of histogram equalization. We then develop a novel automated method to identify CHs from these homogenized maps by determining the intensity threshold of CH regions separately for each synoptic map. This is done by identifying the best location and size of an image segment, which optimally contains portions of coronal holes and the surrounding quiet Sun allowing us to detect the momentary intensity threshold. Our method is thus able to adjust itself to the changing scale size of coronal holes and to temporally varying intensities. To make full use of the information in the three wavelengths we construct a composite CH distribution, which is more robust than distributions based on one wavelength. Using the composite CH dataset we discuss the temporal evolution of CHs during the Solar Cycles 23 and 24.  相似文献   

5.
Polar crown prominences, that partially circle the Sun’s poles between 60° and 70° latitude, are made of chromospheric plasma. We aim to diagnose the 3D dynamics of a polar crown prominence using high-cadence EUV images from the Solar Dynamics Observatory (SDO)/AIA at 304, 171, and 193 Å and the Ahead spacecraft of the Solar Terrestrial Relations Observatory (STEREO-A)/EUVI at 195 Å. Using time series across specific structures, we compare flows across the disk in 195 Å with the prominence dynamics seen on the limb. The densest prominence material forms vertical columns that are separated by many tens of Mm and connected by dynamic bridges of plasma that are clearly visible in 304/171 Å two-colour images. We also observe intermittent but repetitious flows with velocity 15 km?s?1 in the prominence that appear to be associated with EUV bright points on the solar disk. The boundary between the prominence and the overlying cavity appears as a sharp edge. We discuss the structure of the coronal cavity seen both above and around the prominence. SDO/HMI and GONG magnetograms are used to infer the underlying magnetic topology. The evolution and structure of the prominence with respect to the magnetic field seems to agree with the filament-linkage model.  相似文献   

6.
A new method for the automated detection of coronal holes and filaments on the solar disk is presented. The starting point is coronal images taken by the Extreme Ultraviolet Telescope on the Solar and Heliospheric Observatory (SOHO/EIT) in the Fe ix/x 171 Å, Fe xii 195 Å, and He ii 304 Å extreme ultraviolet (EUV) lines and the corresponding full-disk magnetograms from the Michelson Doppler Imager (SOHO/MDI) from different phases of the solar cycle. The images are processed to enhance their contrast and to enable the automatic detection of the two candidate features, which are visually indistinguishable in these images. Comparisons are made with existing databases, such as the He i 10830 Å NSO/Kitt Peak coronal-hole maps and the Solar Feature Catalog (SFC) from the European Grid of Solar Observations (EGSO), to discriminate between the two features. By mapping the features onto the corresponding magnetograms, distinct magnetic signatures are then derived. Coronal holes are found to have a skewed distribution of magnetic-field intensities, with values often reaching 100?–?200 gauss, and a relative magnetic-flux imbalance. Filaments, in contrast, have a symmetric distribution of field intensity values around zero, have smaller magnetic-field intensity than coronal holes, and lie along a magnetic-field reversal line. The identification of candidate features from the processed images and the determination of their distinct magnetic signatures are then combined to achieve the automated detection of coronal holes and filaments from EUV images of the solar disk. Application of this technique to all three wavelengths does not yield identical results. Furthermore, the best agreement among all three wavelengths and NSO/Kitt Peak coronal-hole maps occurs during the declining phase of solar activity. The He ii data mostly fail to yield the location of filaments at solar minimum and provide only a subset at the declining phase or peak of the solar cycle. However, the Fe ix/x 171 Å and Fe xii 195 Å data yield a larger number of filaments than the Hα data of the SFC.  相似文献   

7.
We briefly overview results of our study of the large-scale solar activity associated with coronal mass ejections (CMEs). The observational material is constituted with data of the SPIRIT telescope aboard the CORONAS-F satellite in the three EUV channels 175, 284, and 304 Å. In particular, we consider a powerful geoeffective event of November 4, 2003, which was not observed by the SOHO/EIT telescope, a series of extremely powerful events of October 2003, and an event of November 18, 2003, with filament eruption. The efficiency of combined analysis of the SPIRIT and EIT data is demonstrated. The analysis confirms the coincidence of many dimmings in different spectral channels, including coronal lines with different excitation temperatures and the transition-region line, as well as the global character and homology of dimmings in recurrent events. The higher cadence SPIRIT observations at 304 Å reveal a slowly propagating large-scale darkening probably caused by absorption of emission in the dense, cold plasma of an eruptive filament.  相似文献   

8.
We analyze multiwavelength observations of an M2.9/1N flare that occurred in AR NOAA 11112 on 16 October 2010. AIA 211 Å EUV images reveal the presence of a faster coronal wave (decelerating from ≈?1390 to ≈?830 km?s?1) propagating ahead of a slower wave (decelerating from ≈?416 to ≈?166 km?s?1) towards the western limb. The dynamic radio spectrum from Sagamore Hill radio telescope shows the presence of a metric type II radio burst, which reveals the presence of a coronal shock wave (speed ≈?800 km?s?1). The speed of the faster coronal wave, derived from AIA 211 Å images, is found to be comparable to the coronal shock speed. AIA 171 Å high-cadence observations showed that a coronal loop, which was located at a distance of ≈?0.32R ⊙ to the west of the flaring region, started to oscillate by the end of the impulsive phase of the flare. The results indicate that the faster coronal wave may be the first driver of the transversal oscillations of coronal loop. As the slower wave passed through the coronal loop, the oscillations became even stronger. There was a plasmoid eruption observed in EUV and a white-light CME was recorded, having velocity of ≈?340?–?350 km?s?1. STEREO 195 Å images show an EIT wave, propagating in the same direction as the lower-speed coronal wave observed in AIA, but decelerating from ≈?320 to ≈?254 km?s?1. These observations reveal the co-existence of both waves (i.e. coronal Moreton and EIT waves), and the type II radio burst seems to be associated with the coronal Moreton wave.  相似文献   

9.
We describe the automated extraction of active regions (ARs) or plages from the European Grid of Solar Observations (EGSO) Solar Feature Catalogue using a region-growing technique. In this work, Hα and Ca ii K3 solar images from the Meudon Observatory and EUV solar images from the SOHO/EIT instrument were used. For better detection accuracy, the statistical properties of each quarter of a full disk solar image are used to define local intensity thresholds for an initial segmentation that helps to define AR seeds. Median filtering and morphological operations are applied to the resulting binary image in order to remove noise and to merge broken regions. The centroids of each labelled region are used as seeds, from which a region-growing procedure starts. Statistics-based local thresholding is also applied to compute upper- and lower- threshold intensity values defining the spatial extents of the regions. The detection results obtained with the resulting automated thresholding and region-growing (ATRG) procedure are compared day-by-day with the synoptic maps manually generated by the Meudon Observatory and NOAA for 2 months in 2002 and more coarsely over a 5-year period. The moderate correlation found between our detection results and those produced manually on the other data sets reveals a need for a unified active region definition. As an application of the SFC for ARs we present the tracking of the active region AR NOAA 10484 during its appearance on the solar disk from 19–26 October 2003 and compare its intensity variations for Hα and Fe xii 195 Å wavelengths.  相似文献   

10.
We have analyzed dimmings, i.e., regions of temporarily reduced brightness, and manifestations of a coronal wave in the famous event of 14 July 2000 using images produced with the EUV telescope SOHO/EIT. Our analysis was inspired by a paper by Andrews (2001, Solar Phys. 204, 181 (Paper I)), in which this event was studied using running-difference EIT images at 195 Å formed by subtraction of a previous image from each current one. Such images emphasize changes of the brightness, location, and configuration of observed structures occurring during the 12-min interval between two subsequent heliograms. However, they distort the picture of large-scale disturbances caused by a CME, particularly, dimmings. A real picture of dimmings can be obtained from fixed-base difference ‘de-rotated’ images. The latter are formed in two stages: first, the solar rotation is compensated using three-dimensional rotation of all images (‘de-rotation’) to the time of a pre-event heliogram, here 10:00 UT, and then the base heliogram is subtracted from all others. We show real dimmings to be essentially different from those described by Andrews (Paper I). The restructuring of large-scale magnetic fields in the corona in connection with the CME was accompanied by the appearance and growth of two large dimmings. One of them was located along the central meridian, southward of the eruption center, at the place of the pre-eruption arcade. Another dimming occupied the space between the flare region and a remote western active region. Several smaller dimmings were observed virtually over the whole solar disk, especially, within the northwest quadrant. We have also revealed a propagating disturbance with properties of a coronal wave in the northern polar sector, where no dimmings were observed. This fact is discussed in the context of probable association between dimmings and coronal waves. Having suppressed the ‘snowstorm’ produced in the EIT images by energetic particles, we have considered dimming manifestations in all four EIT pass bands of 171, 195, 284, and 304 Å as well as the light curves of the main dimmings including several later images at 195 Å. Our analysis shows that the major cause of the dimmings was density depletion that reached up to 30% in this event. The picture of dimmings implies that the CME in the Bastille Day event was an octopus-like bundle of some magnetic ropes, with the ‘arms’ being connected to several active regions disposed over almost the whole visible solar surface.  相似文献   

11.
We present identifications of coronal holes (CHs) from observations in the He?i 10?830 Å line made at Kitt Peak Observatory (from 1975 to 2003) and in the EUV 195 Å wavelength with SOHO/EIT (from 1996 to 2012). To determine whether a feature is a CH we have developed semi-automatic techniques for delineating CH borders on synoptic charts and for subsequent mapping of these borders on magnetic-field charts. Using these techniques, we superimposed CH borders on magnetic-field charts over the time interval from 1975 to 2012. A major contribution to the total area was made by high-latitude CHs, but in the declining phase of solar cycle 23, the contribution from low-latitude CHs increased substantially. Variations in the flux of Galactic cosmic rays and those in the inclination angle of the heliospheric current sheet followed the cyclic variations of CH areas. High-latitude CHs affect the properties of the solar wind in the ecliptic plane.  相似文献   

12.
In this study, we present the three-dimensional (3D) configuration of a filament observed by STEREO and the Global High Resolution H-alpha Network (GHN) in EUV 304 Å and Hα line center, respectively. This was the largest filament located close to the active region NOAA 10956 that produced a small B9.6 flare and two Coronal Mass Ejections (CMEs) on 19 May 2007. The 3D coordinates of multiple points traced along this filament were reconstructed by triangulation from two different aspect angles. The two STEREO (A and B) spacecraft had a separation angle α of 8.6 degree on 19 May 2007. The “true” heights of the filament were estimated using STEREO images in EUV 304 and Hα images, respectively. Our results show that EUV emission of the filament originates from higher locations than the Hα emission. We also compare the measured reconstructed heights of the filaments in EUV with those reported in previous studies.  相似文献   

13.
We discuss early results derived from an algorithm that automates the detection, cataloging, and analysis of extreme-ultraviolet (EUV) “bright points” (BP) from 9 years of data acquired by the Extreme-ultraviolet Imaging Telescope (EIT) on the Solar and Heliospheric Observatory (SOHO). The algorithm relies upon the computation of a map of “intensity significance”; this then contains the location of the EUV BPs. By mapping the location of BPs in each image and linking them through long sequences of EIT images we can describe the temporal and spatial variations of the 1.3× 108 EUV BPs observed by SOHO to date. We suggest that there is a considerable amount of physical information about the solar coronal plasma that can be readily gleamed from the BP detection database. In this paper we discuss only a small portion of the possible correlations, but we point to the possibility of BP lifetime distributions that are well described by modified power-laws; the components of which vary with the different temperature filters and with time over the present solar cycle. Dedicated to the memory of John (Ian) Hamilton (1938–2004).  相似文献   

14.
One of the major discoveries of the Extreme ultraviolet Imaging Telescope (EIT) on SOHO was the intensity enhancements propagating over a large fraction of the solar surface. The physical origin(s) of the so-called EIT waves is still strongly debated with either wave (primarily fast-mode MHD waves) or nonwave (pseudo-wave) interpretations. The difficulty in understanding the nature of EUV waves lies in the limitations of the EIT observations that have been used almost exclusively for their study. They suffer from low cadence and single temperature and viewpoint coverage. These limitations are largely overcome by the SECCHI/EUVI observations onboard the STEREO mission. The EUVI telescopes provide high-cadence, simultaneous multitemperature coverage and two well-separated viewpoints. We present here the first detailed analysis of an EUV wave observed by the EUVI disk imagers on 7 December 2007 when the STEREO spacecraft separation was ≈?45°. Both a small flare and a coronal mass ejection (CME) were associated with the wave. We also offer the first comprehensive comparison of the various wave interpretations against the observations. Our major findings are as follows: (1) High-cadence (2.5-minute) 171 Å? images showed a strong association between expanding loops and the wave onset and significant differences in the wave appearance between the two STEREO viewpoints during its early stages; these differences largely disappeared later; (2) the wave appears at the active region periphery when an abrupt disappearance of the expanding loops occurs within an interval of 2.5 minutes; (3) almost simultaneous images at different temperatures showed that the wave was most visible in the 1?–?2 MK range and almost invisible in chromospheric/transition region temperatures; (4) triangulations of the wave indicate it was rather low lying (≈?90 Mm above the surface); (5) forward-fitting of the corresponding CME as seen by the COR1 coronagraphs showed that the projection of the best-fit model on the solar surface was inconsistent with the location and size of the co-temporal EUV wave; and (6) simulations of a fast-mode wave were found in good agreement with the overall shape and location of the observed wave. Our findings give significant support for a fast-mode interpretation of EUV waves and indicate that they are probably triggered by the rapid expansion of the loops associated with the CME.  相似文献   

15.
The increase in the amount of solar data provided by new satellites makes it necessary to develop methods to automate the detection of solar features. Here we present a method for automatically detecting active regions in solar extreme ultraviolet (EUV) images using a series of steps. Initially, the bright regions in the image are segmented using seeded region growing. In a second phase these bright regions are clustered into active regions. Partition-based clustering (both hard and fuzzy) and hierarchical clustering are compared in this work. The aim of the clustering phase is to associate a group to each segmented region in order to reduce the total number of active regions. This facilitates the documentation or subsequent monitoring of these regions. We use two indicators to validate the partitioning: i) the number of detected clusters approximates the number of active regions reported by the National Oceanic and Atmospheric Administration (NOAA) and ii) the area that defines each cluster overlaps with the area of an active region of NOAA. Experiments have been performed on over 6000 images from SOHO/EIT (195 Å). The best results were obtained using hierarchical clustering. The method detects a set of active regions in an image of the solar corona that successfully matches the number of NOAA regions. We will use these regions to perform real-time monitoring and flare detection.  相似文献   

16.
The STEREO mission has been providing a stereoscopic view of filament eruptions in the EUV. The clearest view during a filament eruption is seen in He ii 304 Å observations. One of the main problems in visualizing filament dynamics in He ii 304 Å is the strong background contrast due to surface features. We present a technique that removes background features and leaves behind only the filamentary structure, as seen by STEREO-A and -B. The technique uses a pair of STEREO He ii 304 Å images observed simultaneously. The STEREO-B image is geometrically transformed to a STEREO-A view so that the background images appear similar. Filaments, being elevated structures, i.e., not lying on the same spherical surface as background features, do not appear similar in the transformed view. Thus, subtracting the two images cancels the background but leaves behind the filament structure. We apply this technique to study the dynamics of the filament-eruption event of 22 May 2008, which was observed by STEREO and followed by several ground-based observatories participating in the Joint Observing Programme (JOP 178).  相似文献   

17.
We generated an event catalog with an automated detection algorithm based on the entire EUVI image database observed with the two Solar Terrestrial Relations Observatory (STEREO)-A and -B spacecraft over the first six years of the mission (2006?–?2012). The event catalog includes the heliographic positions of some 20?000 EUV events, transformed from spacecraft coordinates to Earth-based coordinates, and information on associated GOES flare events (down to the level of GOES A5-class flares). The 304 Å wavelength turns out to be the most efficient channel for flare detection (79?% of all EUVI event detections), while the 171 Å (4?%), 195 Å (10?%), and the 284 Å channel (7?%) retrieve substantially fewer flare events, partially due to the suppressing effect of EUV dimming, and partially due to the lower cadence in the later years of the mission. Due to the Sun-circling orbits of STEREO-A and -B, a large number of flares have been detected on the farside of the Sun, invisible from Earth, or seen as partially occulted events. The statistical size distributions of EUV peak fluxes (with a power-law slope of α P =2.5±0.2) and event durations (with a power-law slope of α T =2.4±0.3) are found to be consistent with the fractal-diffusive self-organized criticality model. The EUVI event catalog is available on-line at secchi.lmsal.com/EUVI/euvi_autodetection/euvi_events.txt and may serve as a comprehensive tool to identify stereoscopically observed flare events for 3D reconstruction and to study occulted flare events.  相似文献   

18.
Explosive events appear as broad non-Gaussian wings in the line profiles of small transition-region phenomena. Images from the Solar Dynamics Observatory (SDO) give a first view of the plasma dynamics at the sites of explosive events seen simultaneously in O?vi spectra of a region of quiet Sun, taken with the ultraviolet spectrometer Solar Ultraviolet Measurements of Emitted Radiation (SUMER) onboard the Solar and Heliospheric Observatory (SOHO). Distinct event bursts were seen either at the junction of supergranular network cells or near emerging flux. Three are described in the context of their surrounding transition region (304 Å) and coronal (171 Å) activity. One showed plasma ejections from an isolated pair of sites, with a time lag of 50 seconds between events. At the site where the later explosive event was seen, the extreme ultraviolet (EUV) images show a hot core surrounded by a small, expanding ring of chromospheric emission, which we interpret as a “splash.” The second explosive-event burst was related to flux cancellation, inferred from Helioseismic and Magnetic Imager (HMI) magnetograms, and a coronal dimming surrounded by a ring of bright EUV emission with explosive events at positions where the spectrometer slit crossed the bright ring. The third series of events occurred at the base of a slow, small coronal mass ejection (mini-CME). All events studied here imply jet-like flows probably triggered by magnetic reconnection at supergranular junctions. Events come from sites close to the footpoints of jets seen in Atmospheric Imaging Assembly (AIA) images, and possibly from the landing site of high-velocity flows. They are not caused by rapid rotation in spicules.  相似文献   

19.
The Sun Watcher using Active Pixel system detector and Image Processing (SWAP) onboard the PRoject for OnBoard Autonomy-2 (PROBA2) spacecraft provides images of the solar corona in EUV channel centered at 174 Å. These data, together with the Atmospheric Imaging Assembly (AIA) and the Helioseismic and Magnetic Imager (HMI) onboard Solar Dynamics Observatory (SDO), are used to study the dynamics of coronal bright points. The evolution of the magnetic polarities and associated changes in morphology are studied using magnetograms and multi-wavelength imaging. The morphology of the bright points seen in low-resolution SWAP images and high-resolution AIA images show different structures, whereas the intensity variations with time show similar trends in both SWAP 174 Å and AIA 171 Å channels. We observe that bright points are seen in EUV channels corresponding to a magnetic flux of the order of 1018 Mx. We find that there exists a good correlation between total emission from the bright point in several UV–EUV channels and total unsigned photospheric magnetic flux above certain thresholds. The bright points also show periodic brightenings, and we have attempted to find the oscillation periods in bright points and their connection to magnetic-flux changes. The observed periods are generally long (10?–?25 minutes) and there is an indication that the intensity oscillations may be generated by repeated magnetic reconnection.  相似文献   

20.
A filament disappearance event was observed on 22 May 2008 during our recent campaign JOP 178. The filament, situated in the Southern Hemisphere, showed sinistral chirality consistent with the hemispheric rule. The event was well observed by several observatories, in particular by THEMIS. One day, before the disappearance, Hα observations showed up- and down-flows in adjacent locations along the filament, which suggest plasma motions along twisted flux rope. THEMIS and GONG observations show shearing photospheric motions leading to magnetic flux canceling around barbs. STEREO A, B spacecraft with separation angle 52.4°, showed quite different views of this untwisting flux rope in He ii 304 Å images. Here, we reconstruct the three-dimensional geometry of the filament during its eruption phase using STEREO EUV He ii 304 Å images and find that the filament was highly inclined to the solar normal. The He ii 304 Å movies show individual threads, which oscillate and rise to an altitude of about 120 Mm with apparent velocities of about 100 km?s?1 during the rapid evolution phase. Finally, as the flux rope expands into the corona, the filament disappears by becoming optically thin to undetectable levels. No CME was detected by STEREO, only a faint CME was recorded by LASCO at the beginning of the disappearance phase at 02:00 UT, which could be due to partial filament eruption. Further, STEREO Fe xii 195 Å images showed bright loops beneath the filament prior to the disappearance phase, suggesting magnetic reconnection below the flux rope.  相似文献   

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

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