首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Quantification of 3-D intravascular flow is valuable for studying arterial wall diseases but currently there is a lack of effective clinical tools for this purpose. Divergence-free interpolation (DFI) using radial basis function (RBF) is an emerging approach for full-field flow reconstruction using experimental sparse flow field samples. Previous DFI reconstructs full-field flow from scattered 3-D velocity input obtained using phase-contrast magnetic resonance imaging with low temporal resolution. In this study, a new DFI algorithm is proposed to reconstruct full-field flow from scattered 2-D in-plane velocity vectors obtained using ultrafast contrast-enhanced ultrasound (>1000 fps) and particle imaging velocimetry. The full 3-D flow field is represented by a sum of weighted divergence-free RBFs in space. Because the acquired velocity vectors are only in 2-D and hence the problem is ill-conditioned, a regularized solution of the RBF weighting is achieved through singular value decomposition (SVD) and the L-curve method. The effectiveness of the algorithm is determined via numerical experiments for Poiseuille flow and helical flow with added noise, and it is found that an accuracy as high as 95.6% can be achieved for Poiseuille flow (with 5% input noise). Experimental feasibility is also determined by reconstructing full-field 3-D flow from experimental 2-D ultrasound image velocimetry measurements in a carotid bifurcation phantom. The method is typically faster for a range of problems compared with computational fluid dynamics, and has been found to be effective for the three flow cases.  相似文献   

2.
In this study, a technique for high-frame-rate ultrasound imaging velocimetry (UIV) is extended first to provide more robust quantitative flow velocity mapping using ensemble correlation of images without coherent compounding, and second to generate spatio-temporal wall shear stress (WSS) distribution. A simulation model, which couples the ultrasound simulator with analytical flow solution, was implemented to evaluate its accuracy. It is shown that the proposed approach can reduce errors in velocity estimation by up to 10-fold in comparison with the coherent correlation approach. Mean errors (ME) of 3.2% and 8.6% were estimated under a steady flow condition, while 3.0% and 10.6% were found under a pulsatile condition for the velocity and wall shear rate (WSR) measurement, respectively. Appropriate filter parameters were selected to constrain the velocity profiles before WSR estimations and the effects of incorrect wall tracking were quantified under a controlled environment. Although accurate wall tracking is found to be critical in WSR measurement (as a 200?µm deviation from the wall may yield up to a 60% error), this can be mitigated by HFR imaging (of up to 10?kHz) with contrast agents, which allow for improved differentiation of the wall-fluid boundaries. In vitro investigations on two carotid bifurcation phantoms, normal and diseased, were conducted, and their relative differences in terms of the flow patterns and WSR distribution were demonstrated. It is shown that high-frame-rate UIV technique can be a non-invasive tool to measure quantitatively the spatio-temporal velocity and WSS distribution.  相似文献   

3.
Risk stratification in coronary artery disease is an ongoing challenge for which few tools are available for quantifying physiology within coronary arteries. Recently, anatomy-driven computational fluid dynamic modeling has enabled the mapping of local flow dynamics in coronary stenoses, with derived parameters such as WSS exhibiting a strong capability for predicting adverse clinical events on a patient-specific basis. As cardiac catheterization is common in patients with coronary artery disease, minimally invasive technologies capable of identifying pathologic flow in situ in real time could have a significant impact on clinical decision- making. As a step toward in vivo quantification of slow flow near the arterial wall, proof-of-concept for 3-D intravascular imaging of blood flow dynamics is provided using a 118-element forward-viewing ring array transducer and a research ultrasound system. Blood flow velocity components are estimated in the direction of primary flow using an unfocused wave Doppler approach, and in the lateral and elevation directions, using a transverse oscillation approach. This intravascular 3-D vector velocity system is illustrated by acquiring real-time 3-D data sets in phantom experiments and in vivo in the femoral artery of a pig. The effect of the catheter on blood flow dynamics is also experimentally assessed in flow phantoms with both straight and stenotic vessels. Results indicate that 3-D flow dynamics can be measured using a small form factor device and that a hollow catheter design may provide minimal disturbance to flow measurements in a stenosis (peak velocity: 54.97 ± 2.13 cm/s without catheter vs. 51.37 ± 1.08 cm/s with hollow catheter, 6.5% error). In the future, such technologies could enable estimation of 3-D flow dynamics near the wall in patients already undergoing catheterization.  相似文献   

4.
This article reports the use of a two-dimensional (2-D) capacitive micro-machined ultrasound transducer (CMUT) to acquire radio-frequency (RF) echo data from relatively large volumes of a simple ultrasound phantom to compare three-dimensional (3-D) elasticity imaging methods. Typical 2-D motion tracking for elasticity image formation was compared with three different methods of 3-D motion tracking, with sum-squared difference (SSD) used as the similarity measure. Differences among the algorithms were the degree to which they tracked elevational motion: not at all (2-D search), planar search, combination of multiple planes and plane independent guided search. The cross-correlation between the predeformation and motion-compensated postdeformation RF echo fields was used to quantify motion tracking accuracy. The lesion contrast-to-noise ratio was used to quantify image quality. Tracking accuracy and strain image quality generally improved with increased tracking sophistication. When used as input for a 3-D modulus reconstruction, high quality 3-D displacement estimates yielded accurate and low noise modulus reconstruction. (E-mail: tgfisher@wisc.edu)  相似文献   

5.
Left ventricular (LV) blood flow is an inherently complex time-varying 3-D phenomenon, where 2-D quantification often ignores the effect of out-of-plane motion. In this study, we describe high frame rate 4-D echocardiographic particle image velocimetry (echo-PIV) using a prototype matrix transesophageal transducer and a dynamic LV phantom for testing the accuracy of echo-PIV in the presence of complex flow patterns. Optical time-resolved tomographic PIV (tomo-PIV) was used as a reference standard for comparison. Echo-PIV and tomo-PIV agreed on the general profile of the LV flow patterns, but echo-PIV smoothed out the smaller flow structures. Echo-PIV also underestimated the flow rates at greater imaging depths, where the PIV kernel size and transducer point spread function were large relative to the velocity gradients. We demonstrate that 4-D echo-PIV could be performed in just four heart cycles, which would require only a short breath-hold, providing promising results. However, methods for resolving high velocity gradients in regions of poor spatial resolution are required before clinical translation.  相似文献   

6.
Current Doppler echocardiography quantification of mitral regurgitation (MR) severity has shortcomings. Proximal isovelocity surface area (PISA)-based methods, for example, are unable to account for the fact that ultrasound Doppler can measure only one velocity component: toward or away from the transducer. In the present study, we used ultrasound-based computational fluid dynamics (Ub-CFD) to quantify mitral regurgitation and study its advantages and disadvantages compared with 2-D and 3-D PISA methods. For Ub-CFD, patient-specific mitral valve geometry and velocity data were obtained from clinical ultrasound followed by 3-D CFD simulations at an assumed flow rate. We then obtained the average ratio of the ultrasound Doppler velocities to CFD velocities in the flow convergence region, and scaled CFD flow rate with this ratio as the final measured flow rate. We evaluated Ub-CFD, 2-D PISA and 3-D PISA with an in vitro flow loop, which featured regurgitation flow through (i) a simplified flat plate with round orifice and (ii) a 3-D printed realistic mitral valve and regurgitation orifice. The Ub-CFD and 3-D PISA methods had higher precision than the 2-D PISA method. Ub-CFD had consistent accuracy under all conditions tested, whereas 2-D PISA had the lowest overall accuracy. In vitro investigations indicated that the accuracy of 2-D and 3-D PISA depended significantly on the choice of aliasing velocity. Evaluation of these techniques was also performed for two clinical cases, and the dependency of PISA on aliasing velocity was similarly observed. Ub-CFD was robustly accurate and precise and has promise for future translation to clinical practice.  相似文献   

7.
Blood flow velocity and wall shear stress (WSS) influence and are influenced by vascular disease. Their measurement is consequently useful in the laboratory and clinic. Contrast-enhanced ultrasound image velocimetry (UIV) can estimate them accurately but the need to inject contrast agents limits utility. Singular value decomposition and high-frame-rate imaging may render contrast agents dispensable. Here we determined whether contrast agent-free UIV can measure flow and WSS. In simulation, accurate measurements were achieved with a signal-to-noise ratio of 13.5 dB or higher. Signal intensity in the rabbit aorta increased monotonically with mechanical index; it was lowest during stagnant flow and uneven across the vessel. In vivo measurements with contrast-free and contrast-enhanced UIV differed by 4.4% and 1.9% for velocity magnitude and angle and by 9.47% for WSS. Bland–Altman analysis of waveforms revealed good agreement between contrast-free and contrast-enhanced UIV. In five rabbits, the root-mean-square errors were as low as 0.022 m/s (0.81%) and 0.11 Pa (1.7%). This study indicates that with an optimised protocol, UIV can assess flow and WSS without contrast agents. Unlike contrast-enhanced UIV, contrast-free UIV could be routinely employed.  相似文献   

8.
We developed a duplex ultrasound simulator for training and assessment of scanning skills. We used the simulator to test examiner performance in the measurement of flow velocities in dialysis access fistulas. Test cases were created from 3-D ultrasound scans of two dialysis access fistulas by reconstructing 3-D blood vessel models and simulating blood flow velocity fields within the lumens. The simulator displays a 2-D B-mode or color Doppler image corresponding to transducer position on a mannequin; a spectral waveform is generated according to Doppler sample volume location and system settings. Examiner performance was assessed by comparing the measured peak systolic velocity (PSV) with the true PSV provided by the computational flow model. The PSV measured by four expert examiners deviated from the true value by 7.8?±?6.1%. The results indicate the ability of the simulator to objectively assess an examiner's measurement accuracy in complex vascular targets.  相似文献   

9.
Ultrasound (US) examination of the common carotid artery was compared with a through-plane magnetic resonance imaging (MRI) sequence to validate a recently proposed technique for 3-D US vector flow imaging. Data from the first volunteer examined were used as the training set, before volume flow and peak velocities were calculated for the remaining eight volunteers. Peak systolic velocities (PSVs) and volume flow obtained with 3-D US were, on average, 34% higher and 24% lower than those obtained with MRI, respectively. A high correlation was observed for PSV (r = 0.79), whereas a lower correlation was observed for volume flow (r = 0.43). The overall standard deviations were ±5.7% and ±5.7% for volume flow and PSV with 3-D US, compared with ±2.7% and ±3.2% for MRI. Finally, the data were re-processed with a change in the parameter settings for the echo-canceling filter to investigate its influence on overall performance. PSV was less affected by the re-processing, whereas the difference in volume flow between 3-D vector flow imaging and MRI was reduced to ?9%, and with an improved overall standard deviation of ±4.7%. The results illustrate the feasibility of using 3-D US for precise and angle-independent volume flow and PSV estimation in vivo.  相似文献   

10.
A system for acquisition of 3-D arterial ultrasound geometries and integration with computational fluid dynamics (CFD) is described. The 3-D ultrasound is based on freehand B-mode imaging with positional information obtained using an optical tracking system. A processing chain was established, allowing acquisition of cardiac-gated 3-D data and segmentation of arterial geometries using a manual method and a semi-automated method, 3D meshing and CFD. The use of CFD allowed visualization of flow streamlines, 2-D velocity contours and 3-D wall shear stress. Three-dimensional positional accuracy was 0.17–1.8 mm, precision was 0.06–0.47 mm and volume accuracy was 4.4–15%. Patients with disease and volunteers were scanned, with data collection from one or more of the carotid bifurcation, femoral bifurcation and abdominal aorta. An initial comparison between a manual segmentation method and a semi-automated method suggested some advantages to the semi-automated method, including reduced operator time and the production of smooth surfaces suitable for CFD, but at the expense of over-smoothing in the diseased region. There were considerable difficulties with artefacts and poor image quality, resulting in 3-D geometry data that was unsuitable for CFD. These artefacts were exacerbated in disease, which may mean that future effort, in the integration of 3-D arterial geometry and CFD for clinical use, may best be served using alternative 3-D imaging modalities such as magnetic resonance imaging and computed tomography. (E-mail: P.Hoskins@ed.ac.uk)  相似文献   

11.
To develop an off-line system for three-dimensional (3-D) ultrasound (US) reconstruction of fetoplacental vasculature using colour segmentation and reconstruction software and to determine sources of error in fully freehand ultrasound image acquisition. US images were acquired freehand with the Acuson Sequoia (5C 2-MHz transducer) using power Doppler. After digital transfer to a personal computer, CQ Analysis software (Kinetic Imaging Ltd, Liverpool, UK) was used to segment the colour information from these images, and the resulting 8-bit grey-scale images were used for 3-D rendering using commercial software (VoxBlast, Vaytek Inc., Fairfield, IA, USA). 2-D scanning, software and freehand acquisition accuracy were assessed using a linear test rig and distance and volume phantoms (Dansk Phantom Service Ltd); 2-D scanning accuracy was within 1.3%, and software reconstruction accuracy within 1% for x and y planes and up to 3% for the z plane. Fully freehand acquisition was associated with a 12% to 18% mean percentage error in distance measurement in the plane of acquisition. Volumetric reconstruction inaccuracy was between 1.5% and 19.7% for precisely separated images and between 16.2% and 39.2% for fully freehand image acquisition. Rendered 3-D US vascular images clearly delineated vascular anatomy within the placenta and cord. Fully freehand 3-D US does have a role in off-line reconstruction of vascular anatomy, although variability in the z plane precludes its use for volumetric measurement. (E-mail: a.welsh@ic.ac.uk)  相似文献   

12.
Visualization and quantification of blood flow are considered important for early detection of atherosclerosis and patient-specific diagnosis and intervention. As conventional Doppler imaging is limited to 1-D velocity estimates, 2-D and 3-D techniques are being developed. We introduce an adaptive velocity compounding technique that estimates the 2-D velocity vector field using predominantly axial displacements estimated by speckle tracking from dual-angle plane wave acquisitions. Straight-vessel experiments with a 7.8-MHz linear array transducer connected to a Verasonics Vantage ultrasound system revealed that the technique performed with a maximum velocity magnitude bias and angle bias of –3.7% (2.8% standard deviation) and –0.16° (0.41° standard deviation), respectively. In vivo, complex flow patterns were visualized in two healthy and three diseased carotid arteries and quantified using a vector complexity measure that increased with increasing wall irregularity. This measure could potentially be a relevant clinical parameter which might aid in early detection of atherosclerosis.  相似文献   

13.
A vascular vector flow mapping (VFM) method visualizes 2-D cardiac flow dynamics by estimating the radial component of flow from the Doppler velocities and wall motion velocities using the mass conservation equation. Although VFM provides 2-D flow, the algorithm is applicable only to bounded regions. Here, a modified VFM algorithm, vascular VFM, is proposed so that the velocities are estimated regardless of the flow geometry. To validate the algorithm, a phantom mimicking a carotid artery was fabricated and VFM velocities were compared with optical particle image velocimetry (PIV) data acquired in the same imaged plane. The validation results indicate that given optimal beam angle condition, VFM velocitiy is fairly accurate, where the correlation coefficient R between VFM and PIV velocities is 0.95. The standard deviation of the total VFM error, normalized by the maximum velocity, ranged from 8.1% to 16.3%, whereas the standard deviation of the measured input errors ranged from 8.9% to 12.7% for color flow mapping and from 4.5% to 5.9% for subbeam calculation. These results indicate that vascular VFM is reliable as its accuracy is comparable to that of conventional Doppler-flow images.  相似文献   

14.
A new mathematical method for estimating velocity vectors from color Doppler datasets is proposed to image blood flow dynamics; this method has been called echodynamography or vector flow mapping (VFM). In this method, the concept of stream function is exploited to expand a 2-D distribution of radial velocities in polar coordinates, observed with color Doppler, to a 2-D distribution of velocity vectors. This study was designed to validate VFM using 3-D numerical flow models. Velocity fields were reconstructed from the virtual color Doppler datasets derived from the models. VFM captured the gross features of flow structures and produced comparable images of the distribution of vorticity, which correlated significantly with the original field (for velocity magnitudes, standard error of estimate = 0.003 to 0.007m/s; for vorticity, standard error of estimate = 0.35 to 2.01/s). VFM may be sensitive for depicting flow structures derived from color Doppler velocities with reasonable accuracy. (E-mail: t.uejima@nifty.ne.jp)  相似文献   

15.
We tested the feasibility of real-time, 3-D ultrasound (US) imaging in the brain. The 3-D scanner uses a matrix phased-array transducer of 512 transmit channels and 256 receive channels operating at 2.5 MHz with a 15-mm diameter footprint. The real-time system scans a 65 degrees pyramid, producing up to 30 volumetric scans per second, and features up to five image planes as well as 3-D rendering, 3-D pulsed-wave and color Doppler. In a human subject, the real-time 3-D scans produced simultaneous transcranial horizontal (axial), coronal and sagittal image planes and real-time volume-rendered images of the gross anatomy of the brain. In a transcranial sheep model, we obtained real-time 3-D color flow Doppler scans and perfusion images using bolus injection of contrast agents into the internal carotid artery.  相似文献   

16.
Wang J  Lu JY 《Ultrasonic imaging》2007,29(2):105-121
Based on the high frame rate (HFR) imaging theories, an extended HFR imaging method has been developed recently in our lab where multiple limited-diffraction array beams or steered plane waves are used in transmissions to reconstruct a high quality image of an equivalent dynamic focusing in both transmissions and receptions. The method has the potential to simplify imaging systems because the fast Fourier transform and square-wave aperture weightings can be used. The method is also flexible in using different numbers of transmissions for a continuous trade-off between image quality and frame rate. In this paper, we study the effects of phase aberration and noise on the extended HFR imaging method with in vitro experiments and compare the results with those obtained with a conventional delay-and-sum (D&S) method of a fixed-transmission focus and a dynamically-focused reception. In the experiments, an ATS539 tissue-mimicking phantom and an Acuson V2 phase array transducer (128 elements, 2.5 MHz, and 0.15-mm pitch) were used. The transducer was driven by a homemade general-purpose HFR imaging system that was capable of producing both the limited-diffraction array beams and steeredplane waves and echo data were acquired with the same system and then transferred to a personal computer via a universal serial bus (USB) 2.0 link for image reconstructions. The phase aberration was introduced by adding random phase shifts to both transmission and reception beams. The random noise was added to the received radiofrequency echo data. Results show that the phase aberration and noise degrade both the extended HFR and the conventional delay-and-sum (D&S) imaging method. However, images reconstructed with the extended HFR imaging method have an overall higher quality than those with the D&S method given the phase aberration and noise models studied.  相似文献   

17.
Three-dimensional (3-D) ultrasound imaging has improved evaluation of organ circulation and might contribute new information on maternal and fetal blood supply. Flow index (FI) of 3-D color images has been proposed as a measure of perfusion. The aim of this study was to evaluate whether the 3-D FI is a parameter of volume flow and flow velocity in a human vessel and in a flow phantom. A 1-cm-long strip of the uterine artery was recorded in 3-D power Doppler (3D-PD) mode in a cross-sectional study of 170 normal singleton pregnancies between 26 and 42 weeks' gestation. A fixed ultrasound system installation was used during the examination. The VOCAL software integrated in the ultrasound unit calculated vessel volume and FI. Reproducibility of the measurements was tested. The method was also tested on a commercially available flow phantom. Reproducibility measurements gave satisfactory results, both in terms of inter- and intraobserver variation. Unexpectedly, in normal pregnancy, the uterine artery FI decreased slightly with gestation. Uterine artery vessel volume increased, however, with gestational age. A poor correlation was found between the FI and both flow velocity and volume flow in the flow phantom. In conclusion, 3D-PD imaging can give impressive anatomical pictures of organ vascular tree. However, the new FI is poorly related to flow velocity or volume of flow.  相似文献   

18.
OBJECTIVES: To assess the ability of Doppler-gated 3-D fetal echocardiography to reconstruct and display specific cardiac structures routinely visualized during antenatal ultrasound in a population at low risk for cardiac anomalies. To determine whether any advantage is offered by 3-D sonographic cardiac examination over conventional sonographic fetal screening techniques. DESIGN: After routine two-dimensional sonographic examination, 3-D cardiac data were collected prospectively in 30 fetuses with gestational ages between 19 and 23 weeks from a low risk patient population. Basic echocardiographic key views were derived from 3-D data and selected for reconstruction and analysis. Four- and five-chamber views were rated and only those views judged to be well visualized were considered as positive results. RESULTS: The four- and five-chamber views were well visualized in all but one fetus using conventional 2-D imaging. Gated 3-D volume data sets enabled visualization of these structures in only 19 of 30 fetuses but provided additional structural depth and allowed a dynamic 3-D perspective of valvar morphology and ventricular wall motion. The right ventricular outflow tract was available from the 3-D volumes in 16 subjects. CONCLUSIONS: Considering the versatility of gated 3-D fetal cardiac imaging we believe that it may soon become an important component of fetal screening thus helping to retrieve standard cardiac cross sections when 2-D imaging is limited by lack of sonographer experience or sonographic windows. Diagnostically acceptable echocardiographic views were obtained more consistently with 2-D ultrasound than with 3-D volume data.  相似文献   

19.
Joint arthropathies often require continuous monitoring of the joint condition, typically performed using magnetic resonance (MR) or ultrasound (US) imaging. US imaging is often the preferred screening or diagnostic tool as it is fast and inexpensive. However, conventional 2-D US has limited capability to compare imaging results between examinations because of its operator dependence and challenges related to repeat imaging in the same location and orientation. Comparison between several imaging sessions is crucial to assess the interval progression of joint conditions. We propose a novel 3-D US scanner for ankle joint assessment that can partially overcome these issues by enabling 3-D imaging. Here, we (i) present the design of the 3-D US ankle scanner system, (ii) validate the geometric reconstruction accuracy of the system, (iii) provide preliminary images of healthy volunteer ankles and (iv) compare 3-D US imaging results with MR imaging. The 3-D ankle scanner consists of a tub filled with water, a linear US probe attached to the wall of the tub and a motorized unit that rotates the US probe 360° around the center of the tub. As the probe rotates, a 3-D US image is formed of the ankle of the patient positioned in the middle of the tub. US probe height, angle and distance from the tub center can be adjusted. The reconstruction accuracy of the system was validated in each of the coordinate directions at different probe angles using two test phantoms. A phantom consisting of numerous Ø200-µm nylon threads with known spacing and a metal rod with machined grooves was used for validation in the horizontal and vertical directions, respectively. The volumetric reconstruction accuracy validation was performed by imaging an agar phantom with two embedded spheres of known volumes and comparing the segmented sphere volume and surface area with the expected. Three-dimensional US and MR images of both ankles of five healthy volunteers were acquired. Distal tibia and proximal talus were segmented in both imaging modalities and the surfaces of these segmentations were compared using the 95% Hausdorff and mean surface distances. The observed mean linear measurement error in all the coordinate directions and over several probe angles was 2.98%. The mean measured volumetric measurement error was 3.45%. The volunteer study revealed useful features for joint assessment present in the 3-D ankle scanner images, such as joint spacing, distal tibia and proximal talus. The mean 95% Hausdorff and mean surface distances between segmentations in 3-D US and MR images were 5.68 ± 0.83 and 2.01 ± 0.30 mm, respectively. In this proof-of-concept study, the 3-D US ankle scanner enabled visualization of the ankle joint features that are useful for joint assessment.  相似文献   

20.
High-frequency 3-D color-flow imaging of the microcirculation   总被引:7,自引:0,他引:7  
High-frequency (> 20 MHz) ultrasound (US) flow imaging has the potential to be an important tool for assessing microvascular blood flow in superficial tissues noninvasively. In this paper, we describe the development and evaluation of a 3-D US flow imaging system capable of operating at center frequencies in the 20- to 50-MHz range. Flow images are made for tissue volumes of sizes up to 10 mm laterally and 5 mm in depth, permitting a range of scientific and clinical applications. To acquire data sets in a reasonable time, the 2-D sections were derived from data collected with a transducer that was scanning continuously in a direction perpendicular to the beam axis. Due to spectral broadening effects induced by scanning tissue, significant tradeoffs must be made between frame rate, lateral resolution and the minimum detectable blood velocity. 3-D flow images were reconstructed with flow data acquired from a series of adjacent planes. The system was evaluated at a center frequency of 50 MHz, using two PVDF transducers with lateral resolutions of 43 microm and 65 microm and axial resolutions of 66 microm to 72 microm, respectively. Velocity ranges were from below 1 mm/s to 25 mm/s. In vivo validation experiments using the mouse ear demonstrated the ability to follow branching patterns of closely spaced microvessels from 30 microm to 100 microm in diameter. Experiments conducted on mouse tumors successfully imaged microvessel morphology in the tumor microcirculation.  相似文献   

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

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