首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
This paper describes the development and application of a new rapid, full k-space acquisition imaging pulse sequence based on the rapid acquisition with relaxation enhancement (RARE) methodology. We have termed this pulse sequence single excitation multiple image RARE (SEMI-RARE). We demonstrate the application of SEMI-RARE to the visualisation of a static liquid phantom and it is shown that up to 120 images can be acquired from a single excitation. By exploiting the inherent relaxation and diffusion contrast within the series of images, the SEMI-RARE provides an ultra-fast method for characterising the spatial distribution of chemical species and phases within complex systems. The pulse sequence is then applied to the study of single- and two-phase flow in a single narrow tube of inner diameter 2.9mm. In particular, it is shown that 8 two-dimensional slice-selective images for two-phase bubble-train flow in a single tube can be acquired from a single excitation at effective echo-times of 37, 109, 181, 253, 325, 397, 469, and 541ms. The visualisation enables the determination of gas/liquid bubble sizes and velocities during two-phase flow. We also report the first direct evidence, obtained from magnetic resonance measurements, of liquid re-circulation zones associated with bubble-train flow. The robustness of the SEMI-RARE technique makes it an attractive fast imaging technique for the study of multi-phase flow phenomena, which are often characterised by large variations in magnetic susceptibility, and are of widespread interest in chemical engineering.  相似文献   

2.
A treatment planning system based on magnetic resonance (MR) angiographic imaging data for the radiosurgery of inoperable cerebral arteriovenous malformations is reported. MR angiography was performed using a three-dimensional (3D) velocity-compensated fast imaging with steady-state precession (FISP) sequence. Depending on the individual MR system, inhomogeneities and nonlinearities induced by eddy currents during the pulse sequence can distort the images and produce spurious displacements of the stereotactic coordinates in both the x-y plane and the z axis. If necessary, these errors in position can be assessed by means of two phantoms placed within the stereotactic guidance system--a "2D-phantom" displaying "pincushion" distortion in the image, and a "3D-phantom" displaying displacement, warp, and tilt of the image plane itself. The pincushion distortion can be "corrected" (reducing displacements from 2-3 mm to 1 mm) by calculations based on modeling the distortion as a fourth order 2D polynomial. Displacement, warp, and tilt of the image plane may be corrected by adjustment of the gradient shimming currents. After correction, the accuracy of the geometric information is limited only by the pixel resolution of the image (= 1 mm). Precise definition of the target volume could be performed by the therapist either directly in the MR images or in calculated projection MR angiograms obtained by a maximum intensity projection algorithm. MR angiography provides a sensitive, noninvasive 3D method for defining target volume and critical structures, and for calculating precise dose distributions for radiosurgery of cerebral arteriovenous malformations.  相似文献   

3.
In this paper, we present the correction of the geometric distortion measured in the clinical magnetic resonance imaging (MRI) systems reported in the preceding paper (Part I) using a 3D method based on the phantom-mapped geometric distortion data. This method allows the correction to be made on phantom images acquired without or with the vendor correction applied. With the vendor's 2D correction applied, the method corrects for both the "residual" geometric distortion still present in the plane in which the correction method was applied (the axial plane) and the uncorrected geometric distortion along the axis normal to the plane. The evaluation of the effectiveness of the correction using this new method was carried out through analyzing the residual geometric distortion in the corrected phantom images. The results show that the new method can restore the distorted images in 3D nearly to perfection. For all the MRI systems investigated, the mean absolute deviations in the positions of the control points (along x-, y- and z-axes) measured on the corrected phantom images were all less than 0.2 mm. The maximum absolute deviations were all below approximately 0.8 mm. As expected, the correction of the phantom images acquired with the vendor's correction applied in the axial plane performed equally well. Both the geometric distortion still present in the axial plane after applying the vendor's correction and the uncorrected distortion along the z-axis have all been "restored."  相似文献   

4.
A phantom that can be used for mapping geometric distortion in magnetic resonance imaging (MRI) is described. This phantom provides an array of densely distributed control points in three-dimensional (3D) space. These points form the basis of a comprehensive measurement method to correct for geometric distortion in MR images arising principally from gradient field non-linearity and magnet field inhomogeneity. The phantom was designed based on the concept that a point in space can be defined using three orthogonal planes. This novel design approach allows for as many control points as desired. Employing this novel design, a highly accurate method has been developed that enables the positions of the control points to be measured to sub-voxel accuracy. The phantom described in this paper was constructed to fit into a body coil of a MRI scanner, (external dimensions of the phantom were: 310 mm x 310 mm x 310 mm), and it contained 10,830 control points. With this phantom, the mean errors in the measured coordinates of the control points were on the order of 0.1 mm or less, which were less than one tenth of the voxel's dimensions of the phantom image. The calculated three-dimensional distortion map, i.e., the differences between the image positions and true positions of the control points, can then be used to compensate for geometric distortion for a full image restoration. It is anticipated that this novel method will have an impact on the applicability of MRI in both clinical and research settings, especially in areas where geometric accuracy is highly required, such as in MR neuro-imaging.  相似文献   

5.
Magnetic resonance imaging sequences utilizing limited flip angles and gradient echoes yield rapid (less than 2 min) dynamic images of the cardiovascular system. These images contain both accurate anatomical and functional information. Using a gradient refocused acquisition in the steady state (GRASS) in the CINE mode, we studied the relationship between gradient echo signal intensity and velocity of steady and pulsatile flow in a phantom simulating medium to large vessels. Images were acquired on a 1.5 Tesla system (repetition TIME = 21 ms, echo TIME = 12 ms, flip ANGLE = 30 degrees). Data from each pulse interval were sorted in 16 images. Signal intensities from flow tube lumina and surrounding stationary water jacket were used to calculate contrast ratios which were compared to velocity measurements made with electromagnetic (EM) flow probes outside the magnet room. During steady flow, signal intensity contrast ratios increased with increasing flow and in a 10 mm thick slice, reached a peak at 48 cm/s, and declined for velocities up to 90 cm/s. Changes in instantaneous velocity during pulsatile flow correlated well (r > .88) with signal intensity changes up to a maximum mean velocity of 17 cm/s. Total signal intensity from the lumen for an “R to R” interval correlated extremely well (r > .97) with mean pulsatile flow velocities up to 30 cm/s. The excellent correlation between gradient echo signal intensity and actual flow velocities suggests that this imaging sequence might be useful for evaluating normal and pathologic flow phenomena.  相似文献   

6.
Although it is thought that perfluoro-2,4-dimethyl-3-isopropyl-3-pentyl (PFR-2) is a candidate for electron paramagnetic resonance (EPR) imaging agents because of its high stability, no study has been made yet on the EPR imaging of PFR-2. In this study, EPR imaging of a phantom including PFR-2 and mice that had received PFR-2 was performed by an in vivo EPR imaging system operating at an EPR frequency of 700 MHz equipped with a bridged loop-gap resonator (inner diameter, 41 mm; axial length, 10 mm). Because PFR-2 is insoluble in water, it was dissolved in perfluorocarbon. The PFR-2 solution was put in cylindrical sample tubes with various inner diameters, and these sample tubes were placed together in a larger cylindrical sample tube filled with a physiological saline solution, which was used as a phantom. The spatial resolution was estimated to be about 3 mm on the basis of EPR imaging of the phantom. EPR images of mice that had received a PFR-2 injection via the intraperitoneal route indicated that PFR-2 remained in the peritoneal cavity even 2 days after the injection. This finding suggests that it is possible to perform EPR imaging of experimental animals using PFR-2 as an imaging agent which persists in a biological system. Authors' address: Hidekatsu Yokoyama, National Institute of Advanced Industrial Science and Technology, 2266-98 Anagahora, Shimoshidami, Moriyama-ku, Nagoya, Aichi 463-8560, Japan  相似文献   

7.
A low-cost 35 mm PIV stereoscopic system for liquid flows is presented which has an imaging component cost under US$9000. The system uses an angular configuration, rotating mirror image shifting and in-situ calibration techniques. Image processing algorithms based on cross correlation and bicubic interpolation are also used to calculate the 3D data from the PIV images. Results from an error analysis have shown the system to have in plane errors ranging from 4.15 to 5.95% and out of plane errors of 7.01% providing an f-number of f2 is fixed for all imaging. Subsequent application of the system to a flow field generated by a free falling sphere in wheat syrup have produced results which when compared to previous flow visualisation give good qualitative agreement. Suggested improvements to the PIV system costing US$1300 would allow operation at f-numbers down to f by modifying the cameras for the Scheimpflug condition and using a corrective liquid prism.  相似文献   

8.
Geometric distortion in MR imaging predominantly arises from the inhomogeneity of the static field and the nonlinearity of the gradients. It is the purpose of this paper to analyse the object and machine related contributions to geometric distortion in order to determine which corrections are necessary for attaining a specified precision. System related imperfections were measured by systematic variation of the strength, direction, and polarity of the read-out gradient in imaging experiments on a grid of cylindrical sample tubes. For the 1.5-T system used in this study, static field related errors up to 7 mm and gradient related errors up to 4 mm were observed (midcoronal plane, FOV 400-mm, G-read between 0.5 and 3.0 mT/m). Field related errors were shown to be inversely proportional to gradient strength, whereas gradient related errors turned out to be virtually independent of gradient strength. It therefore seems recommendable to always apply the strongest available selection and read-out gradients when geometric fidelity is given preference to signal-to-noise considerations. Correction of system related geometric distortions in MR images can readily be performed by table lookup. Object-induced distortions of the gradient fields were studied by experiments on a grid of sample tubes immersed into a cylindrical water bath of variable saline concentration. These experiments revealed a negligible influence of the object on the gradient error distribution, and lead to the conclusion that correction for the nonlinearity of the gradients only requires the application of system dependent correction factors. Object-related distortions of B0 were studied by conventional SE and fat-suppressed IR experiments on phantoms and human subjects. In these experiments the polarity of the read-out gradient was reversed. Subtraction images showed significant object-induced inhomogeneities of the static field at tissue-air interfaces and in the immediate vicinity of the object being imaged. A first attempt to correct for object related B0 inhomogeneities was made by contour analysis of the source images. At present this correction still has to be done manually.  相似文献   

9.
The purpose of this investigation was to analyze the normal anatomy and leiomyomas of the uterus with an ultra low field (0.02 T) magnetic resonance imaging (MRI) device. MR imaging was performed on 18 uteri, 11 of which were imaged both preoperatively (in vivo) and as an operative specimen (in vitro), 6 only as an operative specimen, and 1 only preoperatively. All uteri were examined histologically after imaging. The junctional zone was much better delineated in vivo than in vitro, indicating that its appearance on MR is partly due to blood flow. No structures contributing to its visibility in vitro could be demonstrated histologically. Twenty leiomyomas (size range 7–79 mm) in 12 uteri were found with MRI. They were slightly better discerned in vivo than in vitro. The leiomyomas, having no degenerative changes, had a signal intensity which was the same or lower than that of the myometrium. On images obtained in vitro the signal intensity of these leiomyomas relative to that of myometrium correlated directly with their muscular content (R = 0.74, p = .002). The authors conclude that the junctional zone is a sum of physiological and structural factors, the latter being responsible for its in vitro delineation. MR imaging of the uterus in vitro did not give more information than MR imaging in vivo. All leiomyomas larger than 10 mm could be detected, indicating that MR imaging at 0.02 T is an accurate method for the imaging of the uterine leiomyomas.  相似文献   

10.
PurposeObjects falling outside of the true elliptical field-of-view (FOV) in Propeller imaging show unique aliasing artifacts. This study proposes a de-aliasing approach to restore the signal intensities in Propeller images without extra data acquisition.Materials and methodsComputer simulation was performed on the Shepp-Logan head phantom deliberately placed obliquely to examine the signal aliasing. In addition, phantom and human imaging experiments were performed using Propeller imaging with various readouts on a 3.0 Tesla MR scanner. De-aliasing using the proposed method was then performed, with the first low-resolution single-blade image used to find out the aliasing patterns in all the single-blade images, followed by standard Propeller reconstruction. The Propeller images without and with de-aliasing were compared.ResultsComputer simulations showed signal loss at the image corners along with aliasing artifacts distributed along directions corresponding to the rotational blades, consistent with clinical observations. The proposed de-aliasing operation successfully restored the correct images in both phantom and human experiments.ConclusionThe de-aliasing operation is an effective adjunct to Propeller MR image reconstruction for retrospective restoration of aliased signals.  相似文献   

11.
田丽萍  李立立  温文龙  王兴  陈萍  卢裕  王俊锋  赵卫  田进寿 《物理学报》2018,67(18):188501-188501
针对无人机载及星载激光成像雷达系统对条纹管的小型化、高空间分辨率与大探测面积的应用需求,研制了一台具有高边缘空间分辨能力、高亮度增益的小型条纹相机.采用球面光电阴极、球面荧光屏技术提高了条纹相机的边缘空间分辨率和探测面积,有利于增大激光成像雷达的探测视场;采用狭缝型加速电极代替传统栅网电极,有利于提高条纹相机的电耐性和可靠性;设计了加载高达-15 kV工作电压的像缩小型条纹管,增大了条纹管的亮度增益,有助于增大激光雷达系统的探测距离.测试结果显示:在有效工作面积16 mm×2 mm内,条纹管静态空间分辨率高于29.3 lp/mm@MTF=5%(MTF表示调制传递函数),亮度增益高达39.4.条纹相机光电阴极处静态空间分辨率高于15 lp/mm@CTF=11.64%(CTF表示对比度传递函数);边缘动态空间分辨率高于9.8 lp/mm@CTF=5.51%;时间分辨率优于54.6 ps@Tscreen=4.3 ns(Tscreen为全屏时间)且在整个工作面积内具有较高的一致性;动态范围为345:1@54.6 ps.同时,为满足不同的景深及探测精度需求,相机设置六个扫描档位,可以实现不同扫速下的超快速目标诊断.该条纹相机在无人机载及星载激光成像雷达探测中具有潜在的实用价值.  相似文献   

12.
A major determinant of the success of surgical vascular modifications, such as the total cavopulmonary connection (TCPC), is the energetic efficiency that is assessed by calculating the mechanical energy loss of blood flow through the new connection. Currently, however, to determine the energy loss, invasive pressure measurements are necessary. Therefore, this study evaluated the feasibility of the viscous dissipation (VD) method, which has the potential to provide the energy loss without the need for invasive pressure measurements. Two experimental phantoms, a U-shaped tube and a glass TCPC, were scanned in a magnetic resonance (MR) imaging scanner and the images were used to construct computational models of both geometries. MR phase velocity mapping (PVM) acquisitions of all three spatial components of the fluid velocity were made in both phantoms and the VD was calculated. VD results from MR PVM experiments were compared with VD results from computational fluid dynamics (CFD) simulations on the image-based computational models. The results showed an overall agreement between MR PVM and CFD. There was a similar ascending tendency in the VD values as the image spatial resolution increased. The most accurate computations of the energy loss were achieved for a CFD grid density that was too high for MR to achieve under current MR system capabilities (in-plane pixel size of less than 0.4 mm). Nevertheless, the agreement between the MR PVM and the CFD VD results under the same resolution settings suggests that the VD method implemented with a clinical imaging modality such as MR has good potential to quantify the energy loss in vascular geometries such as the TCPC.  相似文献   

13.
The amplitudes of gradient-echoes produced using static field gradients are sensitive to diffusion of tissue water during the echo evolution time. Gradient-echoes have been used to produce MR images in which image intensity is proportional to the self-diffusion coefficient of water. However, such measurements are subject to error due to the presence of background magnetic field gradients caused by variations in local magnetic susceptibility. These local gradients add to the applied gradients. The use of radiofrequency (RF) gradients to produce gradient-echoes may avoid this problem. The RF magnetic field is orthogonal to the offset field produced by local magnetic susceptibility gradients. Thus, the effect of the local gradients on RF gradient-echo amplitude is small if the RF field is strong enough to minimize resonance offset effects. The effects of susceptibility gradients can be further reduced by storing magnetization longitudinally during the echo evolution period. A water phantom was used to evaluate the effects of background gradients on the amplitudes of RF gradient-echoes. A surface coil was used to produce an RF gradient of between 1.3 and 1.6 gauss/cm. Gradient-echoes were detected with and without a 0.16 gauss/cm static magnetic field gradient applied along the same direction as the RF gradient. The background static field gradient had no significant effect on the decay of RF gradient-echo amplitude as a function of echo evolution time. In contrast, the effect of the background gradient on echoes produced using a 1.6 gauss/cm static field gradient is calculated to be significant. This analysis suggests that RF gradient-echoes can produce MR images in which signal intensity is a function of the self-diffusion coefficient of water, but is not significantly affected by background gradients.  相似文献   

14.
In this work, computer modeling based on a finite element method is used to simulate the T2* relaxation of hyperpolarized noble gases (HNG) in the lungs. A physical model of lung airways consisting of a phantom constructed from micro-capillary fibers of diameters similar to the size of lung airways with semi-permeable walls is also presented. The fibers are surrounded by a liquid medium (water) of magnetic susceptibility similar to lung tissue. Theoretical predictions of the field strength dependence of T2* for 129Xe in the phantom and in vivo rat lung are presented. These predictions are in good agreement with experimental T2* values obtained from the phantoms and in vivo rat lungs (160, 19 and 8 ms) at three different field strengths (0.074, 1.89 and 3T, respectively) using hyperpolarized 129Xe. The strong dependence of T2* on field strength is consistent with the theoretical prediction that low fields may be optimal for HNG MR imaging of the lungs as the decreased T2* at high fields necessitates an increase in bandwidth for conventional MR imaging.  相似文献   

15.
Diffusion-weighted echo planar imaging (DW-EPI) suffers from geometric distortion due to low phase-encoding bandwidth. Read-out segmented echo planar imaging (RS-EPI) reduces distortion but residual distortion remains in extreme cases. Additional corrections need to be applied, especially for radiotherapy applications where a high degree of accuracy is needed. In this study the use of magnetic field map corrections are assessed in DW-EPI and RS-EPI, to reduce geometric uncertainty for MRI-guided radiotherapy applications. Magnetic field maps were calculated from gradient echo images and distortion corrections were applied to RS-EPI images. Distortions were assessed in a prostate phantom by comparing to the known geometry, and in vivo using a modified Hausdorff distance metric using a T2-weighted spin echo as ground truth. Across 10 patients, field map-corrected RS-EPI reduced maximum distortion by 5 mm on average compared to DW-EPI (σ = 1.9 mm). Geometric distortions were also reduced significantly using field mapping with RS-EPI, compared to RS-EPI alone (p ≤ 0.05). The increased geometric accuracy of these techniques can potentially allow diffusion-weighted images to be fused with other MR or CT images for radiotherapy treatment purposes.  相似文献   

16.
The purpose was to analyse magnetic susceptibility effects on accuracy of point-wise measurements of signal profiles in the assessment of MRS volume selection performance. An existing phantom design consisting of a sphere with a movable signal source was used for the investigation. The influence from the phantom on magnetic field homogeneity was measured with phase sensitive 1H imaging and 31P spectroscopy on a 1.5 T whole body MR system. The susceptibility effects for such a phantom design can be separated in 1/ A variation in the background magnetic field, which is caused by the stationary structures and has a significant influence on spatial accuracy. 2/ A magnetic field distortion, which is caused by the movable signal source and has very little influence on accuracy. The spatial inaccuracy due to susceptibility effects in this phantom, was 0.03 mm for positions of the signal source covering a 40-mm VOI. Susceptibility effects from the movable signal source were substantial but had very little influence on spatial accuracy. Still, improvements of this phantom design are possible. Point-wise measurements using a phantom with a movable signal source is inherently insensitive to susceptibility effects from the signal source and permits accurate signal profile measurements of high spatial (sub-mm) resolution.  相似文献   

17.
Full field laser Doppler imaging (LDI) and single exposure laser speckle contrast imaging (LSCI) are directly compared using a novel instrument which can concurrently image blood flow using both LDI and LSCI signal processing. Incorporating a commercial CMOS camera chip and a field programmable gate array (FPGA) the flow images of LDI and the contrast maps of LSCI are simultaneously processed by utilizing the same detected optical signals. The comparison was carried out by imaging a rotating diffuser. LDI has a linear response to the velocity. In contrast, LSCI is exposure time dependent and does not provide a linear response in the presence of static speckle. It is also demonstrated that the relationship between LDI and LSCI can be related through a power law which depends on the exposure time of LSCI.  相似文献   

18.
《Ultrasonics》2013,53(1):36-44
Vibro-acoustography (VA) is a medical imaging method based on the difference-frequency generation produced by the mixture of two focused ultrasound beams. VA has been applied to different problems in medical imaging such as imaging bones, microcalcifications in the breast, mass lesions, and calcified arteries. The obtained images may have a resolution of 0.7–0.8 mm. Current VA systems based on confocal or linear array transducers generate C-scan images at the beam focal plane. Images on the axial plane are also possible, however the system resolution along depth worsens when compared to the lateral one. Typical axial resolution is about 1.0 cm. Furthermore, the elevation resolution of linear array systems is larger than that in lateral direction. This asymmetry degrades C-scan images obtained using linear arrays. The purpose of this article is to study VA image restoration based on a 3D point spread function (PSF) using classical deconvolution algorithms: Wiener, constrained least-squares (CLSs), and geometric mean filters. To assess the filters’ performance on the restored images, we use an image quality index that accounts for correlation loss, luminance and contrast distortion. Results for simulated VA images show that the quality index achieved with the Wiener filter is 0.9 (when the index is 1.0 this indicates perfect restoration). This filter yielded the best result in comparison with the other ones. Moreover, the deconvolution algorithms were applied to an experimental VA image of a phantom composed of three stretched 0.5 mm wires. Experiments were performed using transducer driven at two frequencies, 3075 kHz and 3125 kHz, which resulted in the difference-frequency of 50 kHz. Restorations with the theoretical line spread function (LSF) did not recover sufficient information to identify the wires in the images. However, using an estimated LSF the obtained results displayed enough information to spot the wires in the images. It is demonstrated that the phase of the theoretical and the experimental PSFs are dissimilar. This fact prevents VA image restoration with the current theoretical PSF. This study is a preliminary step towards understanding the restoration of VA images through the application of deconvolution filters.  相似文献   

19.
吴海兵  陶声祥  张良  张静 《应用光学》2016,37(5):673-679
鉴于低照度条件下彩色成像都是采用微光、红外、紫外等波段图像融合成假彩色,提出一种利用滤光片过滤和像增强器增强的三基色获取及真彩色融合方法。采用F P滤光片设计出了透射中心在三基色光中心波长的三基色滤光片,对其光谱透过性进行分析;通过旋转三基色滤色轮将目标反射光过滤,使用像增强器对过滤后反射光增强,借助FPGA编写的控制程序实现滤色轮转速和CCD摄像机成像时间的精确同步,利用CCD摄像机获取经增强后的三基色图像,实现目标三基色图像的动态采集。研制的样机在微光实验室进行了三基色采集及融合实验,实验时光照度分别1×10-1 lx (等效于月光下)和1×10-3 lx(等效于星光下),采集速度设置为60 f/s,对获取的图像质量进行了评价。结果表明:在照度为1×10-1 lx时,融合后的真彩色图像在灰度均值、灰度方差和信息熵3项指标方面,比3幅单色图像的平均值分别提高了5.06%、5.97%和1.08%;在照度为1×10-3 lx时,融合后的真彩色图像与3幅单色图像的平均值分别提高了13.18%、-9.86和8.65%。  相似文献   

20.
The principal objective of this study was to verify the utility of the reconstruction imaging technique in the brain perfusion computed tomography (PCT) scan by assessing reductions in the radiation dose and analyzing the generated images. The setting used for image acquisition had a detector coverage of 40 mm, a helical thickness of 0.625 mm, a helical shuttle mode scan type and a rotation time of 0.5 s as the image parameters used for the brain PCT scan. Additionally, a phantom experiment and an animal experiment were carried out. In the phantom and animal experiments, noise was measured in the scanning with the tube voltage fixed at 80 kVp (kilovolt peak) and the level of the adaptive statistical iterative reconstruction (ASIR) was changed from 0% to 100% at 10% intervals. The standard deviation of the CT coefficient was measured three times to calculate the mean value. In the phantom and animal experiments, the absorbed dose was measured 10 times under the same conditions as the ones for noise measurement before the mean value was calculated. In the animal experiment, pencil-type and CT-dedicated ionization chambers were inserted into the central portion of pig heads for measurement. In the phantom study, as the level of the ASIR changed from 0% to 100% under identical scanning conditions, the noise value and dose were proportionally reduced. In our animal experiment, the noise value was lowest when the ASIR level was 50%, unlike in the phantom study. The dose was reduced as in the phantom study.  相似文献   

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

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