首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
This work aimed at developing a rapid and clinically applicable method for the assessment of left atrial size and function using magnetic resonance imaging (MRI). We studied 17 healthy subjects and 26 cardiac patients. Left atrial cine MRI with 50 ms phases was made in 6–12 contiguous long-axis sections encompassing the entire atrial cavity. A volume-time curve was reconstructed to measure the minimum and maximum volumes as well as the fractional volume change, reservoir function, ejection fraction, and mean filling and emptying rates of the left atrium. The image section with the largest left atrial area was then selected and a comparable area-time curve was reconstructed. The atrial phasic areas and functional indices were determined analogously to the volume-based assessment. The contours of atrial area-time and volume-time curves agreed closely in individual subjects. All area-based left atrial measurements distinguished cardiac patients as a group from healthy persons. The combined specificity of the area-based analyses was 92% and the sensitivity, 65%, in identifying abnormal results in individual patients. The accuracy of the area-based data was best for the atrial minimum size, fractional change, reservoir function, and mean filling rate. The estimated time savings with the simplified method were 5 to 6 h per patient. Left atrial size and function can be studied by reconstructing a phasic atrial area-time curve with cine MRI. Atrial enlargement and abnormalities of filling and reservoir function can be reliably identified, but if data on conduit or stroke function are crucial the three-dimensional MRI technique is still recommended.  相似文献   

2.

Purpose

The objective of this study is to determine regional left ventricle (LV) function and temporal heterogeneity of LV wall contraction by analyzing regional time–volume curve (TVC) after Fourier fitting and to assess altered systolic and diastolic functions and temporal indices of myocardial contraction in infarcted segments in comparison with noninfarcted myocardium in patients with myocardial infarction (MI).

Methods

Steady-state cine magnetic resonance (MR) and late gadolinium-enhanced (LGE) MR images were acquired using a 1.5-T MR system in 60 patients with MI. Regional LV function was determined by analyzing regional TVC in 16 segments. The fitted regional TVC was generated by Fourier curve fitting with five harmonics. Regional LV ejection fraction (EF), peak ejection rate (PER), peak filling rate (PFR), time to end-systole and time to peak filling (TPF) were determined from TVC and the first derivative curve.

Results

On LGE MR imaging (MRI), MI was observed in 307 of 960 segments (32.0%). Regional EF and PER averaged in LGE segments were 49.3±14.5% and 2.83±0.65 end-diastolic volume (EDV)/s, significantly lower than those in normal segments (66.7±11.9% and 3.63±0.60 EDV/s, P<.001 and P<.01, respectively). In addition, regional PFR, an index of diastolic function, was significantly reduced in LGE segments (1.94±0.54 vs. 2.86±0.68 EDV/s, P<.01). Time to end-systole and TPF were significantly greater in LGE segments (380.2±57.6 and 169.3±45.4 ms) than in normal segments (300.9±55.1 and 132.3±43.0 ms, P<.01 and P<.01, respectively).

Conclusions

Analysis of regional TVC on cine MRI after Fourier fitting allows quantitative assessment of regional systolic and diastolic LV functions and temporal heterogeneity of LV wall contraction in patients with MI.  相似文献   

3.
BackgroundMRI assessment of aortic pulse wave velocity (PWV) helps predict the risk of vascular events, but the recommended phase contrast sampling rate is faster than what is utilized in most clinical sequences. There are many existing MRI databases obtained for assessment of cardiac output using lower temporal frequency sampling where information might be obtained about aortic stiffness (PWV). In this work, we sought to evaluate whether the Group Delay (GD) method can generate a reproducible measure of stiffness and describe expected age-related stiffening of the aortic arch using lower sampling rates in standard clinical sequences.MethodsPhase contrast (PC) MRI was obtained on the ascending and descending aortic arch in a heterogeneous adult cohort (n = 23; 9 women) spanning over a wide range of ages (ages 24–89, mean 49.4 ± 18.4). Data was collected with standard cardiac MRI protocols for cardiac output evaluation (repetition time = 7.8 ms, views-per-segment = 4, encoding velocity = 200 cm/s). Pulse wave transit times (TT) were computed using the GD method, two other validated automated approaches (cross correlation TT Algorithm by Gaddum and Segment by Medviso), and the manual tangent method. Pressure waveforms from tonometry and flow waveforms from PC MRI were used to assess wave reflections.ResultsGroup Delay and TT-Algorithm showed significant and high retest reproducibility (r = 0.86 for both) as well as high PWV correlation with age (r = 0.93, P-value < 0.00005 and r = 0.96, P-value < 0.00005 respectively) and with each other (r = 0.94, P-value < 0.00001, RMSE = 0.94 m/s). Arbitrary altering of the image acquisition trigger in the GD method introduced error of 10%–13%, but the TT-algorithm error range was 11%–25%.ConclusionGroup Delay enables reproducible assessment of transit time to derive PWV from low temporal resolution clinical cardiac MRI sequences that can also identify age-related stiffening.  相似文献   

4.

Purpose

The objective of this study was to compare multiple methods for estimation of PWV from 4D flow MRI velocity data and to investigate if 4D flow MRI-based PWV estimation with piecewise linear regression modeling of travel-distance vs. travel time is sufficient to discern age-related regional differences in PWV.

Methods

4D flow MRI velocity data were acquired in 8 young and 8 older (age: 23 ± 2 vs. 58 ± 2 years old) normal volunteers. Travel-time and travel-distance were measured throughout the aorta and piecewise linear regression was used to measure global PWV in the descending aorta and regional PWV in three equally sized segments between the top of the aortic arch and the renal arteries. Six different methods for extracting travel-time were compared.

Results

Methods for estimation of travel-time that use information about the whole flow waveform systematically overestimate PWV when compared to methods restricted to the upslope-portion of the waveforms (p < 0.05). In terms of regional PWV, a significant interaction was found between age and location (p < 0.05). The age-related differences in regional PWV were greater in the proximal compared to distal descending aorta.

Conclusion

Care must be taken as different classes of methods for the estimation of travel-time produce different results. 4D flow MRI-based PWV estimation with piecewise linear regression modeling of travel-distance vs. travel time can discern age-related differences in regional PWV well in line with previously reported data.  相似文献   

5.
Obstructive sleep apnea (OSA) is a sleep-related breathing disorder that can cause left ventricular (LV) dysfunction. In patients with OSA, the LV dysfunction is usually evaluated by echocardiography. The purpose of this study was to evaluate whether the use of breathhold cine MRI for the study of LV dysfunction would be feasible and well tolerated by patients with OSA. Six volunteers and five patients underwent a breathhold cine MRI study of the LV using a 1.5 Tesla MR imager. Cine MRI was performed using a breathhold k-space segmented TurboFLASH technique during end-expiration. Systolic thickening of the LV septal wall was 49% +/- 16% in normals vs. 25% +/- 10.5% in patients (p < 0.05). Systolic thickening of the LV free wall was 42% +/- 12% in normals vs. 22% +/- 9% in patients (p < 0.05). There was a significant difference in end-diastolic wall thickness between the two groups. All patients tolerated the procedure well. The total duration of each study was relatively short (less than 11 min). Breathhold MRI techniques can be used to study LV dysfunction in patients with respiratory disability such as OSA.  相似文献   

6.
The loss of pulmonary artery (PA) compliance has significant pathophysiological effect on the right ventricle. Noninvasive and reliable assessment of PA wall stiffness would be an essential determiner of right heart load and a clinically useful factor to assess cardiovascular risk. Two MRI techniques have been proposed for assessing PA stiffness by measuring pulse wave velocity (PWV): transit time (TT) and flow area (QA). However, no data are available that compares the two techniques and evaluates their performance, especially over a wide range of PWV values or at 3.0-T, which is the purpose of the present study. Thirty-three patients with different heart conditions were imaged using optimized high-temporal resolution and high-spatial resolution velocity-encoding MRI sequences. Statistical analysis was conducted to study intermethod, interobserver and intraobserver variabilities. The PWV measurements using TT and QA techniques showed good agreement (P>0.1). The Bland-Altman analysis showed negligible differences between the two methods (mean±S.D.=0.11±0.35 m/s, correlation coefficient r=0.94). The repeated measurements showed low interobserver and intraobserver variabilities, although the S.D. of the differences was larger in the QA technique. The mean±S.D. of the TT/QA measurement differences were −0.05±0.2/0.0±0.36 m/s and 0.02±0.26/0.02±0.39 m/s for the interobserver and intraobserver differences, respectively. In conclusion, each technique has its own advantages and disadvantages. The two techniques result in similar measurements, although the QA method is more subjective due to its dependency on operator intervention.  相似文献   

7.
The aim of this study was to evaluate the capability of using micro-magnetic resonance imaging (MRI) to visualize and characterize atherosclerotic plaques of mouse models. Twenty five apolipoprotein E-knockout mice were fed atherogenic diet, which enabled creation of aortic atherosclerotic plaques. Aortic plaques were examined in vivo by 4.7 T MRI and then characterized ex vivo by 11.7 T three-dimensional MRI. MR images were correlated with subsequent histological confirmation. In vivo 4.7-T MRI demonstrated unevenly thickened aortic walls due to formation of atherosclerotic plaques. Ex vivo 11.7-T MRI enabled not only to acquire full volume-rendered images of the entire vessels but also to characterize plaque components (such as lipid cores and fibrous caps) at any level and any projection, which were confirmed by histological correlation. Micro-MRI provides an excellent imaging tool for basic science to investigate atherosclerosis in small animal models, which may become a supplement to histopathology of atherosclerotic cardiovascular disease.  相似文献   

8.
Magnetic resonance imaging (MRI) is particularly well adapted to the evaluation of aortic distensibility. The calculation of this parameter, based on the change in vessel cross-sectional area per unit change in blood pressure, requires precise delineation of the aortic wall on a series of cine-MR images. Firstly, the study consisted in validating a new automatic method to assess aortic elasticity. Secondly, aortic distensibility was studied for the ascending and descending thoracic aortas in 26 healthy subjects. Two homogeneous groups were available to evaluate the influence of sex and age (with an age limit value of 35 years). The automatic postprocessing method proved to be robust and reliable enough to automatically determine aortic distensibility, even on artefacted images. In the 26 healthy volunteers, a marked decrease in distensibility appears with age, although this decrease is only significant for the ascending aorta (8.97±2.69 10−3 mmHg−1 vs. 5.97±2.02 10−3 mmHg−1). Women have a higher aortic distensibility than men but only significantly at the level of the descending aorta (7.20±1.61 10−3 mmHg−1 vs. 5.05±2.40 10−3 mmHg−1). Through our automatic contouring method, the aortic distensibility from routine cine-MRI has been studied on a healthy subject population providing reference values of aortic stiffness. The aortic distensibility calculation shows that age and sex are causes of aortic stiffness variations in healthy subjects.  相似文献   

9.

Purpose

To assess the feasibility of measuring pulmonary artery (PA) pulse wave velocity (PWV) in children breathing ambient air and 12% oxygen.

Methods

Velocity-encoded phase-contrast MR images of the PA were acquired in 15 children, aged 9–12 years, without evidence of cardiac or pulmonary diseases. PWV was derived as the ratio of flow to area changes during early systole. Each child was scanned twice, in air and after at least 20 minutes into inspiratory hypoxic challenge. Intra-observer and inter-observer variability and repeatability were also compared.

Results

PA PWV, which was successfully measured in all subjects, increased from 1.31 ± 0.32 m/s in air to 1.61 ± 0.58 m/s under hypoxic challenge (p = 0.03). Intra- and inter-observer coefficients of variations were 9.0% and 15.6% respectively. Good correlation within and between observers of r = 0.92 and r = 0.72 respectively was noted for PA PWV measurements. Mean (95% limit of agreement) intra- and inter-observer agreement on Bland–Altman analysis were − 0.02 m/s (− 0.41–0.38 m/s) and -0.28 m/s (− 1.06–0.49 m/s).

Conclusion

PA PWV measurement in children using velocity-encoded MRI is feasible, reproducible and sufficiently sensitive to detect differences in PA compliance between normoxia and hypoxia. This technique can be used to detect early changes of PA compliance and monitor PAH in children.  相似文献   

10.
The purpose of this study was to investigate the feasibility of a noninvasive approach that combines magnetic resonance imaging (MRI) oximetry and flow measurement to obtain the oxygen consumption in the myocardium and in the whole body. Thirteen healthy male volunteers [mean (+/-S.D.) age: 35+/-7 years] underwent this MR study, which included myocardial oxygen consumption (MVO(2)) measurements in 11 subjects and whole-body oxygen consumption (VO(2)) measurements in 8 subjects. In six subjects, both measurements were obtained. Five subjects had repeated MRI measurements of global MVO(2) in order to verify the reproducibility of this approach. The protocol included in vitro blood sample T(2)-%O(2) calibration, coronary sinus (CS) and main pulmonary artery (MPA) T(2) and phase contrast flow measurement and left ventricular (LV) mass calculation. Based on Fick's law, a global measurement of LV MVO(2) and whole-body VO(2) using MRI was feasible. The MVO(2) values were 11+/-3 ml/min per 100 g LV mass. For repeated measurements, differences in MVO(2) of 1 ml/min per 100 g LV mass appear detectable. The whole-body VO(2) values were 3.8+/-0.8 ml/min/kg body weight. MRI techniques that combine CS and MPA T(2), flow and LV mass measurements to quantify MVO(2) and whole-body VO(2) noninvasively in healthy subjects appear feasible, based on their correspondence to previously published work.  相似文献   

11.
Zhang X  Greenleaf JF 《Ultrasonics》2006,44(Z1):e165-e168
Increased stiffness of the arteries has recently gained acceptance as a potential risk for cardiovascular and many other diseases. Pulse wave velocity (PWV) is widely used for estimating the stiffness of an artery. However, PWV is an indicator of average artery stiffness between two measuring points. In addition to measured PWV, the diameter and thickness are needed to calculate the elastic modulus of the artery. We present a new method to generate a torsion wave in the artery wall and measure its propagation speed in the circumferential direction of the artery. The elastic modulus of the artery can be calculated from the torsion wave velocity without the knowledge of the thickness of the artery, which is difficult to measure with accuracy.  相似文献   

12.
Arterial wall stiffness can be associated with various diseases. Arteriosclerosis involves the buildup of plaques within artery walls that stiffen the arteries. The stiffness of an artery can be assessed by measurement of the pulse wave velocity (PWV). Usually, PWV is estimated using the foot-to-foot method. However, the foot of the pressure wave is not very clear due to reflected waves. Also, the blood pressure wave generated by the heart is normally a low frequency wave, hence the time resolution is low. PWV is an average indicator of artery stiffness between the two measuring positions, and therefore cannot easily identify local stiffness. In this paper a force on the arterial wall is generated noninvasively by the radiation force of ultrasound. Techniques for measuring the propagating wave due to this localized force are studied. The excitation force can be either a very short pulse or a modulated cw signal of a few hundred hertz. The temporal resolution of this method, which is in the range of microseconds, is much higher than the conventional pressure PWV method, and therefore allows the wave velocity to be measured accurately over short distances of a few millimeters.  相似文献   

13.
Left ventricle (LV) segmentation in cardiac MRI is an essential procedure for quantitative diagnosis of various cardiovascular diseases. In this paper, we present a novel fully automatic left ventricle segmentation approach based on convolutional neural networks. The proposed network fully takes advantages of the hierarchical architecture and integrate the multi-scale feature together for segmenting the myocardial region of LV. Moreover, we put forward a dynamic pixel-wise weighting strategy, which can dynamically adjust the weight of each pixel according to the segmentation accuracy of upper layer and force the pixel classifier to take more attention on the misclassified ones. By this way, the LV segmentation performance of our method can be improved a lot especially for the apical and basal slices in cine MR images. The experiments on the CAP database demonstrate that our method achieves a substantial improvement compared with other well-know deep learning methods. Beside these, we discussed two major limitations in convolutional neural networks-based semantic segmentation methods for LV segmentation.  相似文献   

14.

Introduction

Cardiac magnetic resonance (CMR) is a unique method to determine regional and local aortic stiffness parameters. Although various methods have been validated, there are no data in patients after acute ST-segment elevation myocardial infarction (STEMI). In the present study we assessed the feasibility of different CMR derived measures of aortic stiffness in patients after first acute STEMI for the first time.

Methods

CMR derived aortic pulse wave velocity (PWV) determined by the regional transit-time (PWVTT) and local flow-area (PWVQA) method as well as local distensibility coefficients (DCs) was analyzed in 22 healthy young volunteers and 28 patients with recent acute STEMI.

Results

PWVTT and DC of the ascending aorta differed significantly between healthy subjects and STEMI patients (all p < 0.001). PWVQA at thoracic levels of aorta was not different between groups (p > 0.520) and did not correlate with age (p > 0.149) and PWVTT (p > 0.310). Intra- and interobserver variability was high for PWVTT (r = 0.970, p < 0.001 and r = 0.920, p < 0.001), acceptable for DC (all r > 0. 809, p < 0.001 and all r > 0.510, p < 0.001) but low for thoracic PWVQA (all r < 0.330 and all r < 0.372).

Conclusion

PWVTT and local DC are robust methods for the assessment of aortic stiffness in patients after acute STEMI.  相似文献   

15.
This study describes early diastolic inflow dynamics based on three-directional magnetic resonance velocity data and investigates age-dependent changes in early diastolic inflow characteristics. We examined 26 young healthy volunteers age 25 (3) years (mean, SD), and 23 healthy older volunteers age 63 (8) years. Three-directional magnetic resonance velocity mapping was performed in a long axis plane through the heart. Transverse velocity profiles were read in five different positions in the early diastolic inflow stream of the left ventricle. The size and timing of the maximum velocities at each level were recorded and the repeatability of the method was tested. Compared with the younger group, the older group was characterized by: 1) lower maximum velocity in all positions, 2) increased deceleration of blood downstream from the mitral leaflet tips, and 3) delayed velocity propagation. The described method was repeatable and enabled detection of the age-dependent differences between groups of normal subjects. In conclusion, the early diastolic inflow pattern changes with age, probably reflecting changes in the diastolic function of the myocardium.  相似文献   

16.
This study was done to improve the ability of magnetic resonance (MR) imaging to provide clear cross-sectional images of the ascending and descending aorta in diastole. The study was motivated by interest in measuring the regional compliance of the ascending aorta, which requires determination of the change in cross sectional area of the vessel between systole and diastole. In diastolic images, residual signal from slow flowing blood and flow artifact consistently obscured the inner boundary of the aortic wall and precluded tracing and measurement of the cross sectional area. We concluded that cross sectional area measurement of the ascending aorta was impossible on our system using standard spin echo sequences. To improve wall delineation in diastolic images, SAT pulses were optimized with respect to pulse timing, slice thickness, and gap. Optimized SAT pulses greatly improved the delineation of the vessel wall by removing unwanted signal from flowing spins. Measurement precision was vastly improved by running two scans with and without flow compensation, and correlating visually and numerically the area measurements from each. We established that each image should be measured by two independent observers and traced three times by each. Using these procedures, diastolic cross-sectional areas of the mid-ascending aorta could be measured with a precision of 2.5%, and the change of cross-sectional area between systole and diastole could be measured with a precision of 10.8%. These measurements were precise enough to detect CAD patients with low aortic compliance from the age-matched controls previously reported in one study. The test based on cross sectional area measurement, with a false positive detection rate of 5%, had a false negative rate of 58%. Compliance measurements by MR at 1.5 T could become clinically useful if normal and abnormal populations are sufficiently separated.  相似文献   

17.
To determine whether differences exist between neurofibromatosis type 1 (NF1) patients with or without focal lesions and healthy normal volunteers in the metabolite ratios of normal appearing white matter, 27 patients with NF1 (with parenchymal lesion, MR positive, n: 17; without parenchymal lesions, MR negative, n: 10) and 20 healthy volunteers underwent MRI and short TE (31 ms) proton MR spectroscopy (MRS). In 17 patients with parenchymal lesions, 61 focal lesions were detected by MRI. MRS was performed from normal appearing frontal and posterior parietal white matter (FWM and PWM) in NF1 and from control groups. NAA/Cr, Cho/Cr and MI/Cr ratios were calculated. Significant increase in Cho/Cr and MI/Cr ratios were found in FWM and PWM in MR negative and positive groups when compared to control group. NAA/Cr ratio in MR positive group was significantly decreased in FWM compared to control group. There were no significant differences between FWM and PWM in all metabolite ratios of MR negative group. MI/Cr ratio in MR positive group was significantly elevated in PWM compared to FWM. Metabolite changes detected by MRS could indicate demyelination and gliosis in normal appearing white matter in all NF1 patients, and additionally neuroaxonal damage in the FWM of NF1 patients with focal lesions. For that reason, in the clinical evaluation and follow-up of these patients MRS features of normal appearing white matter should be considered in addition to focal lesions.  相似文献   

18.
OBJECTIVE: Monitoring of articular cartilage repair after matrix-associated autologous chondrocyte implantation (MACI) by a new grading and point-scale system based on noninvasive cartilage-specific magnetic resonance imaging (MRI) protocol. PATIENTS AND METHODS: In 20 patients, postoperative high-resolution MRI follow-up examinations at 4, 12, 24 and 52 weeks after matrix-based ACI for cartilage repair were initiated. The repair tissue was described with separate variables: degree of defect repair in width and length, surface, structure and signal intensity of the repair tissue, and status of the subchondral lamina and bone. For these variables, a grading system with point-scale evaluation was applied, and the mean average values were calculated for every follow-up MR exam of each patient. RESULTS: In 10 patients, an incomplete filling of the defect improved to complete filling (6 patients) or less incomplete (4 patients) filling of the defect. Three cases of implant hypertrophy returned to normal within 1 year. Complete filling of the defect by repair tissue was found in 2 patients from the beginning. Integration was complete in 10 cases. Improvement of incomplete to complete integration was found in 3 patients. The signal intensity of the implant developed to native cartilage signal in 13 patients. The mean average values increased from the 4th to the 52 nd week in 17 of 20 patients and decreased in 3 of 20 patients. CONCLUSION: High-resolution MRI provides a noninvasive tool for monitoring the development of cartilage repair tissue in MACI over time and helps to differentiate abnormal repair tissue from a normal maturation process.  相似文献   

19.

Purpose

The purpose of this study was to investigate the combined effect of hypertension and type 2 diabetes mellitus (DM2) on aortic stiffness and endothelial dysfunction by using an integrated MRI approach.

Materials and Methods

A total of 31 non-hypertensive DM2 patients and 31 hypertensive DM2 patients underwent 3.0-T MRI. Aortic distensibility (AD), pulse wave velocity (PWV) and brachial artery flow-mediated dilation (FMD) were assessed. Student's t-test, Mann–Whitney U test, chi-squared test, Pearson correlation analysis, and univariable and multiple linear regression analyses were used for statistical analyses.

Results

The hypertensive patients showed lower AD at multiple levels (ascending aorta [AA]: 2.07 ± 0.98 × 10− 3 mm Hg− 1 vs. 3.21 ± 1.70 × 10− 3 mm Hg− 1, p < 0.01; proximal thoracic descending aorta [PDA]: 2.58 ± 0.72 × 10− 3 mm Hg− 1 vs. 3.58 ± 1.47 × 10− 3 mm Hg− 1, p < 0.01; distal descending aorta [DDA]: 3.11 ± 1.84 × 10− 3 mm Hg− 1 vs. 4.27 ± 1.75 × 10− 3 mm Hg− 1, p < 0.01); faster PWV (7.46 ± 2.28 m/s vs. 5.82 ± 1.12 m/s, p < 0.05) and lower FMD (12.67% ± 6.49% vs. 20.66% ± 9.7%; p < 0.01). Systolic blood pressure was an independent predictor of PWV, AA-AD, DDA-AD and FMD. FMD was statistically significantly associated with PWV (r = − 0.37, p < 0.01) and AD (p < 0.01).

Conclusions

Hypertension has a contributive effect on aortic stiffness and endothelial dysfunction in DM2 patients.  相似文献   

20.
左心室心肌最为发达,心肌收缩产生的高压将动脉血泵入全身,集中体现了心脏的泵血能力.定量分析左心室收缩运动是诊断心血管疾病(如心肌梗死)的重要途径.本文采用描述左心室心肌材质的生物力学模型重建左心室位移场.该力学模型作为插值项,与心脏电影磁共振图像的观测位移场共同纳入贝叶斯估计框架,并采用有限元法求解位移场方程.实验比较了左心室射血无力组(46例)与正常组(55例)的左心室功能参数,发现两组在径向和圆周方向的位移、速度、应变和应变率都具有非常显著的差异(p < 0.001),这证明本文方法能够有效区别左心室运动正常与否.实验结果还与CVI软件测量的左心室功能参数具有较高的相关性,说明本文方法有望辅助心血管疾病的临床诊断.  相似文献   

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

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