首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
We compare line-of-sight magnetograms from the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) and the Michelson Doppler Imager (MDI) onboard the Solar and Heliospheric Observatory (SOHO). The line-of-sight magnetic signal inferred from the calibrated MDI data is greater than that derived from the HMI data by a factor of 1.40. This factor varies somewhat with center-to-limb distance. An upper bound to the random noise for the 1′′ resolution HMI 720-second magnetograms is 6.3 Mx?cm?2, and 10.2 Mx?cm?2 for the 45-second magnetograms. Virtually no p-mode leakage is seen in the HMI magnetograms, but it is significant in the MDI magnetograms. 12-hour and 24-hour periodicities are detected in strong fields in the HMI magnetograms. The newly calibrated MDI full-disk magnetograms have been corrected for the zero-point offset and underestimation of the flux density. The noise is 26.4 Mx?cm?2 for the MDI one-minute full-disk magnetograms and 16.2 Mx?cm?2 for the five-minute full-disk magnetograms observed with four-arcsecond resolution. The variation of the noise over the Sun’s disk found in MDI magnetograms is likely due to the different optical distortions in the left- and right-circular analyzers, which allows the granulation and p-mode to leak in as noise. Saturation sometimes seen in sunspot umbrae in MDI magnetograms is caused by the low intensity and the limitation of the onboard computation. The noise in the HMI and MDI line-of-sight magnetic-field synoptic charts appears to be fairly uniform over the entire map. The noise is 2.3 Mx?cm?2 for HMI charts and 5.0 Mx?cm?2 for MDI charts. No evident periodicity is found in the HMI synoptic charts.  相似文献   

2.
Sixty days of Doppler images from the Solar and Heliospheric Observatory (SOHO) / Michelson Doppler Imager (MDI) investigation during the 1996 and 2008 solar minima have been analyzed to show that certain supergranule characteristics (size, size range, and horizontal velocity) exhibit fluctuations of three?to?five days. Cross-correlating parameters showed a good, positive correlation between supergranulation size and size range, and a moderate, negative correlation between size range and velocity. The size and velocity do exhibit a moderate, negative correlation, but with a small time lag (less than 12 hours). Supergranule sizes during five days of co-temporal data from MDI and the Solar Dynamics Observatory (SDO) / Helioseismic Magnetic Imager (HMI) exhibit similar fluctuations with a high level of correlation between them. This verifies the solar origin of the fluctuations, which cannot be caused by instrumental artifacts according to these observations. Similar fluctuations are also observed in data simulations that model the evolution of the MDI Doppler pattern over a 60-day period. Correlations between the supergranule size and size range time-series derived from the simulated data are similar to those seen in MDI data. A simple toy-model using cumulative, uncorrelated exponential growth and decay patterns at random emergence times produces a time-series similar to the data simulations. The qualitative similarities between the simulated and the observed time-series suggest that the fluctuations arise from stochastic processes occurring within the solar convection zone. This behavior, propagating to surface manifestations of supergranulation, may assist our understanding of magnetic-field-line advection, evolution, and interaction.  相似文献   

3.
We compare photospheric line-of-sight magnetograms from the Synoptic Optical Long-term Investigations of the Sun (SOLIS) Vector Spectro-Magnetograph (VSM) instrument with observations from the 150-foot Solar Tower at Mt. Wilson Observatory (MWO), the Helioseismic and Magnetic Imager (HMI) on the Solar Dynamics Observatory (SDO), and the Michelson Doppler Imager (MDI) on the Solar and Heliospheric Observatory (SOHO). We find very good agreement between VSM and the other data sources for both disk-averaged flux densities and pixel-by-pixel measurements. We show that the VSM mean flux density time series is of consistently high signal-to-noise ratio with no significant zero offsets. We discuss in detail some of the factors ?C spatial resolution, flux dependence, and position on the solar disk ?C affecting the determination of scaling between VSM and SOHO/MDI or SDO/HMI magnetograms. The VSM flux densities agree well with spatially smoothed data from MDI and HMI, although the scaling factors show a clear dependence on flux density. The factor to convert VSM to HMI increases with increasing flux density (from ??1 to ??1.5). The nonlinearity is smaller for the VSM vs. SOHO/MDI scaling factor (from ??1 to ??1.2).  相似文献   

4.
We present a new sigmoid catalog covering the duration of the Hinode mission and the Solar Dynamics Observatory (SDO) until the end of 2012. The catalog consists of 72 mostly long-lasting sigmoids. We collect and make available all X-ray and EUV data from Hinode, SDO, and the Solar TErrestrial RElations Observatory (STEREO), and we determine the sigmoid lifetimes, sizes, and aspect ratios. We also collect the line-of-sight magnetograms from the Helioseismic and Magnetic Imager (HMI) for SDO or the Michelson Doppler Imager (MDI) on the Solar and Heliospheric Observatory (SOHO) to measure flux versus time for the lifetime of each region. We determine that the development of a sigmoidal shape and eruptive activity is more strongly correlated with flux cancelation than with emergence. We find that the eruptive properties of the regions correlate well with the maximum flux, largest change, and net change in flux in the regions. These results have implications for constraining future flux-rope models of ARs and gaining insight into their evolutionary properties.  相似文献   

5.
We describe the imaging quality of the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) as measured during the ground calibration of the instrument. We describe the calibration techniques and report our results for the final configuration of HMI. We present the distortion, modulation transfer function, stray light, image shifts introduced by moving parts of the instrument, best focus, field curvature, and the relative alignment of the two cameras. We investigate the gain and linearity of the cameras, and present the measured flat field.  相似文献   

6.
L. Gy?ri 《Solar physics》2012,280(2):365-378
Sunspot and white light facular areas are important data for solar activity and are used, for example, in the study of the evolution of sunspots and their effect on solar irradiance. Solar Dynamic Observatory??s Helioseismic and Magnetic Imager (SDO/HMI) solar images have much higher resolution (??0.5????pixel?1) than Solar and Heliospheric Observatory??s Michelson Doppler Imager (SOHO/MDI) solar images (??2????pixel?1). This difference in image resolution has a significant impact on the sunspot and white light facular areas measured in the two datasets. We compare the area of sunspots and white light faculae derived from SDO/HMI and SOHO/MDI observations. This comparison helps the calibration of the SOHO sunspot and facular area to those in SDO observations. We also find a 0.22 degree difference between the North direction in SDO/HMI and SOHO/MDI images.  相似文献   

7.
The Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) provides a new tool for the systematic observation of white-light flares, including Doppler and magnetic information as well as continuum. In our initial analysis of the highly impulsive $\mathrm{\gamma}$ -ray flare SOL2010-06-12T00:57 (Martínez Oliveros et al., Solar Phys. 269, 269, 2011), we reported the signature of a strong blueshift in the two footpoint sources. Concerned that this might be an artifact due to aliasing peculiar to the HMI instrument, we undertook a comparative analysis of Global Oscillation Network Group (GONG++) observations of the same flare, using the PArametric Smearing Correction ALgorithm (PASCAL) algorithm to correct for artifacts caused by variations in atmospheric smearing. This analysis confirms the artifactual nature of the apparent blueshift in the HMI observations, finding weak redshifts at the footpoints instead. We describe the use of PASCAL with GONG++ observations as a complement to the SDO observations and discuss constraints imposed by the use of HMI far from its design conditions. With proper precautions, these data provide rich information on flares and transients.  相似文献   

8.
The Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) is designed to study oscillations and the magnetic field in the solar photosphere. It observes the full solar disk in the Fe?i absorption line at 6173 Å. We use the output of a high-resolution, 3D, time-dependent, radiation-hydrodynamic simulation based on the CO 5 BOLD code to calculate profiles F(??,x,y,t) for the Fe?i 6173 Å line. The emerging profiles F(??,x,y,t) are multiplied by a representative set of HMI filter-transmission profiles R i (??, 1??i??6) and filtergrams I i (x,y,t; 1??i??6) are constructed for six wavelengths. Doppler velocities V HMI(x,y,t) are determined from these filtergrams using a simplified version of the HMI pipeline. The Doppler velocities are correlated with the original velocities in the simulated atmosphere. The cross-correlation peaks near 100 km, suggesting that the HMI Doppler velocity signal is formed rather low in the solar atmosphere. The same analysis is performed for the SOHO/MDI Ni?i line at 6768 Å. The MDI Doppler signal is formed slightly higher at around 125 km. Taking into account the limited spatial resolution of the instruments, the apparent formation height of both the HMI and MDI Doppler signal increases by 40 to 50 km. We also study how uncertainties in the HMI filter-transmission profiles affect the calculated velocities.  相似文献   

9.
The Solar Dynamics Observatory (SDO) was launched on 11 February 2010 at 15:23 UT from Kennedy Space Center aboard an Atlas V 401 (AV-021) launch vehicle. A?series of apogee-motor firings lifted SDO from an initial geosynchronous transfer orbit into a circular geosynchronous orbit inclined by 28° about the longitude of the SDO-dedicated ground station in New Mexico. SDO began returning science data on 1 May 2010. SDO is the first space-weather mission in NASA’s Living With a Star (LWS) Program. SDO’s main goal is to understand, driving toward a predictive capability, those solar variations that influence life on Earth and humanity’s technological systems. The SDO science investigations will determine how the Sun’s magnetic field is generated and structured, how this stored magnetic energy is released into the heliosphere and geospace as the solar wind, energetic particles, and variations in the solar irradiance. Insights gained from SDO investigations will also lead to an increased understanding of the role that solar variability plays in changes in Earth’s atmospheric chemistry and climate. The SDO mission includes three scientific investigations (the Atmospheric Imaging Assembly (AIA), Extreme Ultraviolet Variability Experiment (EVE), and Helioseismic and Magnetic Imager (HMI)), a spacecraft bus, and a dedicated ground station to handle the telemetry. The Goddard Space Flight Center built and will operate the spacecraft during its planned five-year mission life; this includes: commanding the spacecraft, receiving the science data, and forwarding that data to the science teams. The science investigations teams at Stanford University, Lockheed Martin Solar Astrophysics Laboratory (LMSAL), and University of Colorado Laboratory for Atmospheric and Space Physics (LASP) will process, analyze, distribute, and archive the science data. We will describe the building of SDO and the science that it will provide to NASA.  相似文献   

10.
The Helioseismic and Magnetic Imager (HMI) instrument will produce Doppler-velocity and vector-magnetic-field maps of the solar surface, whose accuracy is dependent on a thorough knowledge of the transmission profiles of the components of the HMI optical-filter system. Here we present a series of wavelength-dependence calibration tests, performed on the instrument from 2005 onwards, to obtain these profiles. We obtained the transmittances as a function of wavelength for the tunable and non-tunable filter elements, as well as the variation of these transmittances with temperature and the angle of incidence of rays of light. We also established the presence of fringe patterns produced by interferences inside the blocking filter and the front window, as well as a change in transmitted intensity with the tuning position. This thorough characterization of the HMI-filter system confirmed the very high quality of the instrument, and showed that its properties are well within the required specifications to produce superior data with high spatial and temporal resolution.  相似文献   

11.
We report observations of a white-light solar flare (SOL2010-06-12T00:57, M2.0) observed by the Helioseismic Magnetic Imager (HMI) on the Solar Dynamics Observatory (SDO) and the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). The HMI data give us the first space-based high-resolution imaging spectroscopy of a white-light flare, including continuum, Doppler, and magnetic signatures for the photospheric Fe i line at 6173.34 Å and its neighboring continuum. In the impulsive phase of the flare, a bright white-light kernel appears in each of the two magnetic footpoints. When the flare occurred, the spectral coverage of the HMI filtergrams (six equidistant samples spanning ±172 mÅ around nominal line center) encompassed the line core and the blue continuum sufficiently far from the core to eliminate significant Doppler crosstalk in the latter, which is otherwise a possibility for the extreme conditions in a white-light flare. RHESSI obtained complete hard X-ray and γ-ray spectra (this was the first γ-ray flare of Cycle 24). The Fe i line appears to be shifted to the blue during the flare but does not go into emission; the contrast is nearly constant across the line profile. We did not detect a seismic wave from this event. The HMI data suggest stepwise changes of the line-of-sight magnetic field in the white-light footpoints.  相似文献   

12.
The Very Fast Inversion of the Stokes Vector (VFISV) is a Milne–Eddington spectral line inversion code used to determine the magnetic and thermodynamic parameters of the solar photosphere from observations of the Stokes vector in the 6173 Å Fe i line by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). We report on the modifications made to the original VFISV inversion code in order to optimize its operation within the HMI data pipeline and provide the smoothest solution in active regions. The changes either sped up the computation or reduced the frequency with which the algorithm failed to converge to a satisfactory solution. Additionally, coding bugs which were detected and fixed in the original VFISV release are reported here.  相似文献   

13.
In Fall 2008 NASA selected a large international consortium to produce a comprehensive automated feature-recognition system for the Solar Dynamics Observatory (SDO). The SDO data that we consider are all of the Atmospheric Imaging Assembly (AIA) images plus surface magnetic-field images from the Helioseismic and Magnetic Imager (HMI). We produce robust, very efficient, professionally coded software modules that can keep up with the SDO data stream and detect, trace, and analyze numerous phenomena, including flares, sigmoids, filaments, coronal dimmings, polarity inversion lines, sunspots, X-ray bright points, active regions, coronal holes, EIT waves, coronal mass ejections (CMEs), coronal oscillations, and jets. We also track the emergence and evolution of magnetic elements down to the smallest detectable features and will provide at least four full-disk, nonlinear, force-free magnetic field extrapolations per day. The detection of CMEs and filaments is accomplished with Solar and Heliospheric Observatory (SOHO)/Large Angle and Spectrometric Coronagraph (LASCO) and ground-based Hα data, respectively. A?completely new software element is a trainable feature-detection module based on a generalized image-classification algorithm. Such a trainable module can be used to find features that have not yet been discovered (as, for example, sigmoids were in the pre-Yohkoh era). Our codes will produce entries in the Heliophysics Events Knowledgebase (HEK) as well as produce complete catalogs for results that are too numerous for inclusion in the HEK, such as the X-ray bright-point metadata. This will permit users to locate data on individual events as well as carry out statistical studies on large numbers of events, using the interface provided by the Virtual Solar Observatory. The operations concept for our computer vision system is that the data will be analyzed in near real time as soon as they arrive at the SDO Joint Science Operations Center and have undergone basic processing. This will allow the system to produce timely space-weather alerts and to guide the selection and production of quicklook images and movies, in addition to its prime mission of enabling solar science. We briefly describe the complex and unique data-processing pipeline, consisting of the hardware and control software required to handle the SDO data stream and accommodate the computer-vision modules, which has been set up at the Lockheed-Martin Space Astrophysics Laboratory (LMSAL), with an identical copy at the Smithsonian Astrophysical Observatory (SAO).  相似文献   

14.
The Helioseismic and Magnetic Imager (HMI) began near-continuous full-disk solar measurements on 1 May 2010 from the Solar Dynamics Observatory (SDO). An automated processing pipeline keeps pace with observations to produce observable quantities, including the photospheric vector magnetic field, from sequences of filtergrams. The basic vector-field frame list cadence is 135 seconds, but to reduce noise the filtergrams are combined to derive data products every 720 seconds. The primary 720 s observables were released in mid-2010, including Stokes polarization parameters measured at six wavelengths, as well as intensity, Doppler velocity, and the line-of-sight magnetic field. More advanced products, including the full vector magnetic field, are now available. Automatically identified HMI Active Region Patches (HARPs) track the location and shape of magnetic regions throughout their lifetime. The vector field is computed using the Very Fast Inversion of the Stokes Vector (VFISV) code optimized for the HMI pipeline; the remaining 180° azimuth ambiguity is resolved with the Minimum Energy (ME0) code. The Milne–Eddington inversion is performed on all full-disk HMI observations. The disambiguation, until recently run only on HARP regions, is now implemented for the full disk. Vector and scalar quantities in the patches are used to derive active region indices potentially useful for forecasting; the data maps and indices are collected in the SHARP data series, hmi.sharp_720s. Definitive SHARP processing is completed only after the region rotates off the visible disk; quick-look products are produced in near real time. Patches are provided in both CCD and heliographic coordinates. HMI provides continuous coverage of the vector field, but has modest spatial, spectral, and temporal resolution. Coupled with limitations of the analysis and interpretation techniques, effects of the orbital velocity, and instrument performance, the resulting measurements have a certain dynamic range and sensitivity and are subject to systematic errors and uncertainties that are characterized in this report.  相似文献   

15.
In our previous articles (Chertok et al. in Solar Phys. 282, 175, 2013; Chertok et al. in Solar Phys. 290, 627, 2015), we presented a preliminary tool for the early diagnostics of the geoeffectiveness of solar eruptions based on the estimate of the total unsigned line-of-sight photospheric magnetic flux in accompanying extreme ultraviolet (EUV) arcades and dimmings. This tool was based on the analysis of eruptions observed during 1996?–?2005 with the Extreme-ultraviolet Imaging Telescope (EIT) and the Michelson Doppler Imager (MDI) onboard the Solar and Heliospheric Observatory (SOHO). Empirical relationships were obtained to estimate the probable importance of upcoming space weather disturbances caused by an eruption, which just occurred, without data on the associated coronal mass ejections. In particular, it was possible to estimate the intensity of a non-recurrent geomagnetic storm (GMS) and Forbush decrease (FD), as well as their onset and peak times. After 2010?–?2011, data on solar eruptions are obtained with the Atmospheric Imaging Assembly (AIA) and the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). We use relatively short intervals of overlapping EIT–AIA and MDI–HMI detailed observations, and additionally, a number of large eruptions over the next five years with the 12-hour cadence EIT images to adapt the SOHO diagnostic tool to SDO data. We show that the adopted brightness thresholds select practically the same areas of arcades and dimmings from the EIT 195 Å and AIA 193 Å image, with a cross-calibration factor of 3.6?–?5.8 (5.0?–?8.2) for the AIA exposure time of 2.0 s (2.9 s). We also find that for the same photospheric areas, the MDI line-of-sight magnetic flux systematically exceeds the HMI flux by a factor of 1.4. Based on these results, the empirical diagnostic relationships obtained from SOHO data are adjusted to SDO instruments. Examples of a post-diagnostics based on SDO data are presented. As before, the tool is applicable to non-recurrent GMSs and FDs caused by nearly central eruptions from active regions, provided that the southern component of the interplanetary magnetic field near the Earth is predominantly negative, which is not predicted by this tool.  相似文献   

16.
The Multi-Application Solar Telescope (MAST) is a 50 cm off-axis Gregorian telescope that has recently become operational at the Udaipur Solar Observatory (USO). An imaging spectropolarimeter is being developed as one of the back-end instruments of MAST to gain a better understanding of the evolution and dynamics of solar magnetic and velocity fields. This system consists of a narrow-band filter and a polarimeter. The polarimeter includes a linear polarizer and two sets of liquid crystal variable retarders (LCVRs). The instrument is intended for simultaneous observations in the spectral lines 6173 Å and 8542 Å, which are formed in the photosphere and chromosphere, respectively. In this article, we present results from the characterization of the LCVRs for the spectral lines of interest and the response matrix of the polarimeter. We also present preliminary observations of an active region obtained using the spectropolarimeter. For verification purposes, we compare the Stokes observations of the active region obtained from the Helioseismic Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) with that of MAST observations in the spectral line 6173 Å. We find good agreement between the two observations, considering the fact that MAST observations are limited by seeing.  相似文献   

17.
We test the reliability of helioseismic far-side active-region predictions, made using Dopplergrams from both the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) and the Global Oscillation Network Group (GONG), by comparison with far-side observation of solar activity from the Solar TErrestrial RElations Observatory (STEREO). Both GONG and HMI produce seismic Carrington maps that show strong magnetic-field regions, labeling predictions of far-side active regions that have a probability ≥?70 %. By visual comparison of these prediction maps with STEREO extreme ultraviolet (EUV) Carrington maps, we determine whether or not solar activity, as evidenced as brightness in EUV, is observed at the predicted locations. We analyzed nine months of data from 2011 and 2012. For both GONG and HMI, we find that for approximately 90 % of the active-region predictions, activity/brightness is observed in EUV at the predicted location. We also investigated the success of GONG and HMI at predicting large active regions before they appear at the east limb as viewed from Earth. Of the 27 identified large east-limb active regions in the nine months of data analyzed, GONG predicted 15 (55 %) at least once within the week prior to Earth-side appearance and HMI predicted 13 (48 %). Based on the STEREO far-side EUV observations, we suggest that 9 of the 27 active regions were probably too weak to be predicted while on the far side. Overall, we conclude that HMI and GONG have similar reliability using the current data-processing procedures.  相似文献   

18.
The Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO/HMI) provides continuous full-disk observations of solar oscillations. We develop a data-analysis pipeline based on the time–distance helioseismology method to measure acoustic travel times using HMI Doppler-shift observations, and infer solar interior properties by inverting these measurements. The pipeline is used for routine production of near-real-time full-disk maps of subsurface wave-speed perturbations and horizontal flow velocities for depths ranging from 0 to 20?Mm, every eight hours. In addition, Carrington synoptic maps for the subsurface properties are made from these full-disk maps. The pipeline can also be used for selected target areas and time periods. We explain details of the pipeline organization and procedures, including processing of the HMI Doppler observations, measurements of the travel times, inversions, and constructions of the full-disk and synoptic maps. Some initial results from the pipeline, including full-disk flow maps, sunspot subsurface flow fields, and the interior rotation and meridional flow speeds, are presented.  相似文献   

19.
Numerical MHD simulations of 3D reconnection events in the solar corona have improved enormously over the last few years, not only in resolution, but also in their complexity, enabling more and more realistic modeling. Various ways to obtain the initial magnetic field, different forms of solar atmospheric models as well as diverse driving speeds and patterns have been employed. This study considers differences between simulations with stratified and non-stratified solar atmospheres, addresses the influence of the driving speed on the plasma flow and energetics, and provides quantitative formulas for mapping electric fields and dissipation levels obtained in numerical simulations to the corresponding solar quantities. The simulations start out from a potential magnetic field containing a null-point, obtained from a Solar and Heliospheric Observatory (SOHO) Michelson Doppler Imager (MDI) magnetogram magnetogram extrapolation approximately 8?hours before a C-class flare was observed. The magnetic field is stressed with a boundary motion pattern similar to?–?although simpler than?–?horizontal motions observed by SOHO during the period preceding the flare. The general behavior is nearly independent of the driving speed, and is also very similar in stratified and non-stratified models, provided only that the boundary motions are slow enough. The boundary motions cause a build-up of current sheets, mainly in the fan-plane of the magnetic null-point, but do not result in a flare-like energy release. The additional free energy required for the flare could have been partly present in non-potential form at the initial state, with subsequent additions from magnetic flux emergence or from components of the boundary motion that were not represented by the idealized driving pattern.  相似文献   

20.
The Helioseismic and Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory produces line-of-sight (LOS) observables (Doppler velocity, magnetic-field strength, Fe i line width, line depth, and continuum intensity) as well as vector magnetic-field maps at the solar surface. The accuracy of LOS observables is dependent on the algorithm used to translate a sequence of HMI filtergrams into the corresponding observables. Using one hour of high-cadence imaging spectropolarimetric observations of a sunspot in the Fe i line at 6173 Å through the Interferometric Bidimensional Spectrometer installed at the Dunn Solar Telescope, and the Milne–Eddington inversion of the corresponding Stokes vectors, we test the accuracy of the observables algorithm currently implemented in the HMI data-analysis pipeline: the MDI-like algorithm. In an attempt to improve the accuracy of HMI observables, we also compare this algorithm to others that may be implemented in the future: a least-squares fit with a Gaussian profile, a least-squares fit with a Voigt profile, and the use of second Fourier coefficients in the MDI-like algorithm.  相似文献   

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

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