首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.

Background

In recent years, new microscopic imaging techniques have evolved to allow us to visualize several different proteins (or other biomolecules) in a visual field. Analysis of protein co-localization becomes viable because molecules can interact only when they are located close to each other. We present a novel approach to align images in a multi-tag fluorescence image stack. The proposed approach is applicable to multi-tag bioimaging systems which (a) acquire fluorescence images by sequential staining and (b) simultaneously capture a phase contrast image corresponding to each of the fluorescence images. To the best of our knowledge, there is no existing method in the literature, which addresses simultaneous registration of multi-tag bioimages and selection of the reference image in order to maximize the overall overlap between the images.

Methodology/Principal Findings

We employ a block-based method for registration, which yields a confidence measure to indicate the accuracy of our registration results. We derive a shift metric in order to select the Reference Image with Maximal Overlap (RIMO), in turn minimizing the total amount of non-overlapping signal for a given number of tags. Experimental results show that the Robust Alignment of Multi-Tag Bioimages (RAMTaB) framework is robust to variations in contrast and illumination, yields sub-pixel accuracy, and successfully selects the reference image resulting in maximum overlap. The registration results are also shown to significantly improve any follow-up protein co-localization studies.

Conclusions

For the discovery of protein complexes and of functional protein networks within a cell, alignment of the tag images in a multi-tag fluorescence image stack is a key pre-processing step. The proposed framework is shown to produce accurate alignment results on both real and synthetic data. Our future work will use the aligned multi-channel fluorescence image data for normal and diseased tissue specimens to analyze molecular co-expression patterns and functional protein networks.  相似文献   

2.

Objectives

To evaluate the accuracy of advanced non-linear registration of serial lung Computed Tomography (CT) images using Large Deformation Diffeomorphic Metric Mapping (LDDMM).

Methods

Fifteen cases of lung cancer with serial lung CT images (interval: 62.2±26.9 days) were used. After affine transformation, three dimensional, non-linear volume registration was conducted using LDDMM with or without cascading elasticity control. Registration accuracy was evaluated by measuring the displacement of landmarks placed on vessel bifurcations for each lung segment. Subtraction images and Jacobian color maps, calculated from the transformation matrix derived from image warping, were generated, which were used to evaluate time-course changes of the tumors.

Results

The average displacement of landmarks was 0.02±0.16 mm and 0.12±0.60 mm for proximal and distal landmarks after LDDMM transformation with cascading elasticity control, which was significantly smaller than 3.11±2.47 mm and 3.99±3.05 mm, respectively, after affine transformation. Emerged or vanished nodules were visualized on subtraction images, and enlarging or shrinking nodules were displayed on Jacobian maps enabled by highly accurate registration of the nodules using LDDMM. However, some residual misalignments were observed, even with non-linear transformation when substantial changes existed between the image pairs.

Conclusions

LDDMM provides accurate registration of serial lung CT images, and temporal subtraction images with Jacobian maps help radiologists to find changes in pulmonary nodules.  相似文献   

3.

Background

Image registration is to produce an entire scene by aligning all the acquired image sequences. A registration algorithm is necessary to tolerance as much as possible for intensity and geometric variation among images. However, captured image views of real scene usually produce unexpected distortions. They are generally derived from the optic characteristics of image sensors or caused by the specific scenes and objects.

Methods and Findings

An analytic registration algorithm considering the deformation is proposed for scenic image applications in this study. After extracting important features by the wavelet-based edge correlation method, an analytic registration approach is then proposed to achieve deformable and accurate matching of point sets. Finally, the registration accuracy is further refined to obtain subpixel precision by a feature-based Levenberg-Marquardt (FLM) method. It converges evidently faster than most other methods because of its feature-based characteristic.

Conclusions

We validate the performance of proposed method by testing with synthetic and real image sequences acquired by a hand-held digital still camera (DSC) and in comparison with an optical flow-based motion technique in terms of the squared sum of intensity differences (SSD) and correlation coefficient (CC). The results indicate that the proposed method is satisfactory in the registration accuracy and quality of DSC images.  相似文献   

4.
WY Hsu 《PloS one》2012,7(7):e40558

Background

A common registration problem for the application of consumer device is to align all the acquired image sequences into a complete scene. Image alignment requires a registration algorithm that will compensate as much as possible for geometric variability among images. However, images captured views from a real scene usually produce different distortions. Some are derived from the optic characteristics of image sensors, and others are caused by the specific scenes and objects.

Methodology/Principal Findings

An image registration algorithm considering the perspective projection is proposed for the application of consumer devices in this study. It exploits a multiresolution wavelet-based method to extract significant features. An analytic differential approach is then proposed to achieve fast convergence of point matching. Finally, the registration accuracy is further refined to obtain subpixel precision by a feature-based modified Levenberg-Marquardt method. Due to its feature-based and nonlinear characteristic, it converges considerably faster than most other methods. In addition, vignette compensation and color difference adjustment are also performed to further improve the quality of registration results.

Conclusions/Significance

The performance of the proposed method is evaluated by testing the synthetic and real images acquired by a hand-held digital still camera and in comparison with two registration techniques in terms of the squared sum of intensity differences (SSD) and correlation coefficient (CC). The results indicate that the proposed method is promising in registration accuracy and quality, which are statistically significantly better than other two approaches.  相似文献   

5.

Purpose

To investigate the effect of B-spline-based elastic image registration on adaptive optics scanning laser ophthalmoscopy (AO-SLO)-assisted capillary visualization.

Methods

AO-SLO videos were acquired from parafoveal areas in the eyes of healthy subjects and patients with various diseases. After nonlinear image registration, the image quality of capillary images constructed from AO-SLO videos using motion contrast enhancement was compared before and after B-spline-based elastic (nonlinear) image registration performed using ImageJ. For objective comparison of image quality, contrast-to-noise ratios (CNRS) for vessel images were calculated. For subjective comparison, experienced ophthalmologists ranked images on a 5-point scale.

Results

All AO-SLO videos were successfully stabilized by elastic image registration. CNR was significantly higher in capillary images stabilized by elastic image registration than in those stabilized without registration. The average ratio of CNR in images with elastic image registration to CNR in images without elastic image registration was 2.10 ± 1.73, with no significant difference in the ratio between patients and healthy subjects. Improvement of image quality was also supported by expert comparison.

Conclusions

Use of B-spline-based elastic image registration in AO-SLO-assisted capillary visualization was effective for enhancing image quality both objectively and subjectively.  相似文献   

6.

Purpose

Superimposition of two dimensional preoperative and postoperative facial images, including radiographs and photographs, are used to evaluate the surgical changes after orthognathic surgery. Recently, three dimensional (3D) imaging has been introduced allowing more accurate analysis of surgical changes. Surface based registration and voxel based registration are commonly used methods for 3D superimposition. The aim of this study was to evaluate and compare the accuracy of the two methods.

Materials and methods

Pre-operative and 6 months post-operative cone beam CT scan (CBCT) images of 31 patients were randomly selected from the orthognathic patient database at the Dental Hospital and School, University of Glasgow, UK. Voxel based registration was performed on the DICOM images (Digital Imaging Communication in Medicine) using Maxilim software (Medicim-Medical Image Computing, Belgium). Surface based registration was performed on the soft and hard tissue 3D models using VRMesh (VirtualGrid, Bellevue City, WA). The accuracy of the superimposition was evaluated by measuring the mean value of the absolute distance between the two 3D image surfaces. The results were statistically analysed using a paired Student t-test, ANOVA with post-hoc Duncan test, a one sample t-test and Pearson correlation coefficient test.

Results

The results showed no significant statistical difference between the two superimposition methods (p<0.05). However surface based registration showed a high variability in the mean distances between the corresponding surfaces compared to voxel based registration, especially for soft tissue. Within each method there was a significant difference between superimposition of the soft and hard tissue models.

Conclusions

There were no significant statistical differences between the two registration methods and it was unlikely to have any clinical significance. Voxel based registration was associated with less variability. Registering on the soft tissue in isolation from the hard tissue may not be a true reflection of the surgical change.  相似文献   

7.

Aim

To investigate the clinical application of a technique for patient set-up verification in breast cancer radiotherapy based on a 3D surface image registration system.

Background

Accurate and reproducible patient set-up is a prerequisite to correctly deliver fractionated radiotherapy. Various approaches are available to verify and correct patient setup for 3D image acquisition in a radiation treatment room.

Materials and methods

The study analyzed the setup reproducibility of 15 patients affected by breast cancer and candidates for conformal radiotherapy by using the AlignRT system (VisionRT, London, UK). At the initial setup, electronic portal imaging device (EPID) images were compared with Digitally Reconstructed Radiographs (DRRs) and a reference three-dimensional (3D) surface image was obtained by AlignRT. Surface images were acquired prior to every subsequent setup procedure. The systematic and random errors along longitudinal and vertical directions were measured and compared for the two systems.

Results

The procedure for surface registration, image acquisition and comparison with the reference image took less than 1 min on average. The T test for systematic error showed no significant difference between the 2 verification systems along the longitudinal (p = 0.69) and vertical (p = 0.67) axes. The T-test for random error showed a significant difference between the 2 systems along the vertical axis (p = 0.05).

Conclusion

AlignRT is fast, simple, non-invasive and seems to be reliable in detecting patient setup errors. Our results suggest that it could be used to assess the setup reproducibility for breast cancer patients.  相似文献   

8.

Background and Aims

Although many methods exist for quantifying the number of pollen grains in a sample, there are few standard methods that are user-friendly, inexpensive and reliable. The present contribution describes a new method of counting pollen using readily available, free image processing and analysis software.

Methods

Pollen was collected from anthers of two species, Carduus acanthoides and C. nutans (Asteraceae), then illuminated on slides and digitally photographed through a stereomicroscope. Using ImageJ (NIH), these digital images were processed to remove noise and sharpen individual pollen grains, then analysed to obtain a reliable total count of the number of grains present in the image. A macro was developed to analyse multiple images together. To assess the accuracy and consistency of pollen counting by ImageJ analysis, counts were compared with those made by the human eye.

Key Results and Conclusions

Image analysis produced pollen counts in 60 s or less per image, considerably faster than counting with the human eye (5–68 min). In addition, counts produced with the ImageJ procedure were similar to those obtained by eye. Because count parameters are adjustable, this image analysis protocol may be used for many other plant species. Thus, the method provides a quick, inexpensive and reliable solution to counting pollen from digital images, not only reducing the chance of error but also substantially lowering labour requirements.  相似文献   

9.

Background

Cone-beam Computed Tomography (CBCT) and stereophotography are two of the latest imaging modalities available for three-dimensional (3-D) visualization of craniofacial structures. However, CBCT provides only limited information on surface texture. This can be overcome by combining the bone images derived from CBCT with 3-D photographs. The objectives of this study were 1) to evaluate the feasibility of integrating 3-D Photos and CBCT images 2) to assess degree of error that may occur during the above processes and 3) to identify facial regions that would be most appropriate for 3-D image registration.

Methodology

CBCT scans and stereophotographic images from 29 patients were used for this study. Two 3-D images corresponding to the skin and bone were extracted from the CBCT data. The 3-D photo was superimposed on the CBCT skin image using relatively immobile areas of the face as a reference. 3-D colour maps were used to assess the accuracy of superimposition were distance differences between the CBCT and 3-D photo were recorded as the signed average and the Root Mean Square (RMS) error.

Principal Findings:

The signed average and RMS of the distance differences between the registered surfaces were −0.018 (±0.129) mm and 0.739 (±0.239) mm respectively. The most errors were found in areas surrounding the lips and the eyes, while minimal errors were noted in the forehead, root of the nose and zygoma.

Conclusions

CBCT and 3-D photographic data can be successfully fused with minimal errors. When compared to RMS, the signed average was found to under-represent the registration error. The virtual 3-D composite craniofacial models permit concurrent assessment of bone and soft tissues during diagnosis and treatment planning.  相似文献   

10.
11.
Huang CB  Zhou J  Zhou Y  Lu ZL 《PloS one》2010,5(12):e15075

Background

How the visual system combines information from the two eyes to form a unitary binocular representation of the external world is a fundamental question in vision science that has been the focus of many psychophysical and physiological investigations. Ding & Sperling (2006) measured perceived phase of the cyclopean image, and developed a binocular combination model in which each eye exerts gain control on the other eye''s signal and over the other eye''s gain control. Critically, the relative phase of the monocular sine-waves plays a central role.

Methodology/Principal Findings

We used the Ding-Sperling paradigm but measured both the perceived contrast and phase of cyclopean images in three hundred and eighty combinations of base contrast, interocular contrast ratio, eye origin of the probe, and interocular phase difference. We found that the perceived contrast of the cyclopean image was independent of the relative phase of the two monocular gratings, although the perceived phase depended on the relative phase and contrast ratio of the monocular images. We developed a new multi-pathway contrast-gain control model (MCM) that elaborates the Ding-Sperling binocular combination model in two ways: (1) phase and contrast of the cyclopean images are computed in separate pathways, although with shared cross-eye contrast-gain control; and (2) phase-independent local energy from the two monocular images are used in binocular contrast combination. With three free parameters, the model yielded an excellent account of data from all the experimental conditions.

Conclusions/Significance

Binocular phase combination depends on the relative phase and contrast ratio of the monocular images but binocular contrast combination is phase-invariant. Our findings suggest the involvement of at least two separate pathways in binocular combination.  相似文献   

12.
Single Echo MRI     

Purpose

Previous nonlinear gradient research has focused on trajectories that reconstruct images with a minimum number of echoes. Here we describe sequences where the nonlinear gradients vary in time to acquire the image in a single readout. The readout is designed to be very smooth so that it can be compressed to minimal time without violating peripheral nerve stimulation limits, yielding an image from a single 4 ms echo.

Theory and Methods

This sequence was inspired by considering the code of each voxel, i.e. the phase accumulation that a voxel follows through the readout, an approach connected to traditional encoding theory. We present simulations for the initial sequence, a low slew rate analog, and higher resolution reconstructions.

Results

Extremely fast acquisitions are achievable, though as one would expect, SNR is reduced relative to the slower Cartesian sampling schemes because of the high gradient strengths.

Conclusions

The prospect that nonlinear gradients can acquire images in a single <10 ms echo makes this a novel and interesting approach to image encoding.  相似文献   

13.

Aim

To analyze intrafraction movement in patients undergoing frameless robotic radiosurgery and evaluate the influence of image acquisition frequency on global accuracy.

Background

Stereotactic radiosurgery requires high spatial accuracy in dose delivery. In conventional radiosurgery, a rigid frame is used to guarantee a correct target alignment and no subsequent movement. Frameless radiosurgery with thermoplastic mask for immobilization cannot completely eliminate intrafraction patient movement. In such cases, it is necessary to evaluate its influence on global treatment accuracy.

Materials and methods

We analyzed the intrafraction motion of the first 15 patients undergoing intracranial radiosurgery (39 fractions) with the CyberKnife VSI system at our institution. Patient position was measured at a 15–90-s interval and was used to estimate intrafraction patient movement.

Results

With our acquisition image protocol and immobilization device, the 99% displacement error was lower than 0.85 mm. The systematic movement components were lower than 0.05 mm and the random component was lower than 0.3 mm in the 3 translational axes. Clear linear time dependence was found in the random component.

Conclusions

Selection of the X-ray image acquisition time is necessary to meet the accuracy required for radiosurgery procedures with the CyberKnife VSI system. We verified that our image acquisition protocol met the 1-mm criterion.  相似文献   

14.

Purpose

The aim of this study was to develop a robust post-processing workflow for motion-corrupted datasets in diffusion kurtosis imaging (DKI).

Materials and methods

The proposed workflow consisted of brain extraction, rigid registration, distortion correction, artifacts rejection, spatial smoothing and tensor estimation. Rigid registration was utilized to correct misalignments. Motion artifacts were rejected by using local Pearson correlation coefficient (LPCC). The performance of LPCC in characterizing relative differences between artifacts and artifact-free images was compared with that of the conventional correlation coefficient in 10 randomly selected DKI datasets. The influence of rejected artifacts with information of gradient directions and b values for the parameter estimation was investigated by using mean square error (MSE). The variance of noise was used as the criterion for MSEs. The clinical practicality of the proposed workflow was evaluated by the image quality and measurements in regions of interest on 36 DKI datasets, including 18 artifact-free (18 pediatric subjects) and 18 motion-corrupted datasets (15 pediatric subjects and 3 essential tremor patients).

Results

The relative difference between artifacts and artifact-free images calculated by LPCC was larger than that of the conventional correlation coefficient (p<0.05). It indicated that LPCC was more sensitive in detecting motion artifacts. MSEs of all derived parameters from the reserved data after the artifacts rejection were smaller than the variance of the noise. It suggested that influence of rejected artifacts was less than influence of noise on the precision of derived parameters. The proposed workflow improved the image quality and reduced the measurement biases significantly on motion-corrupted datasets (p<0.05).

Conclusion

The proposed post-processing workflow was reliable to improve the image quality and the measurement precision of the derived parameters on motion-corrupted DKI datasets. The workflow provided an effective post-processing method for clinical applications of DKI in subjects with involuntary movements.  相似文献   

15.

Background

The eye gaze of other individuals conveys important social information and can trigger multiple psychological activities; some of which, such as emotional reactions and attention orienting, occur very rapidly. Although some neuroscientific evidence has suggested that the amygdala may be involved in such rapid gaze processing, no evidence has been reported concerning the speed at which the amygdala responds to eye gaze.

Methodology/Principal Findings

To investigate this issue, we recorded electrical activity within the amygdala of six subjects using intracranial electrodes. Subjects observed images of eyes and mosaics pointing in averted and straight directions. The amygdala showed higher gamma-band oscillations for eye gaze than for mosaics, which peaked at 200 ms regardless of the direction of the gaze.

Conclusion

These results indicate that the human amygdala rapidly processes eye gaze.  相似文献   

16.

Purpose

To investigate the feasibility of an intravascular imaging antenna to image abdominal aorta atherosclerotic plaque in swine using 3.0T magnetic resonance imaging (MRI).

Methods

Atherosclerotic model was established in 6 swine. After 8 months, swine underwent an MR examination, which was performed using an intravascular imaging guide-wire, and images of the common iliac artery and the abdominal aorta were acquired. Intravascular ultrasound (IVUS) was performed in the right femoral artery; images at the same position as for the MR examination were obtained. The luminal border and external elastic membrane of the targeted arteries were individually drawn in the MR and IVUS images. After co-registering these images, the vessel, lumen, and vessel wall areas and the plaque burden in the same lesions imaged using different modalities were calculated and compared. The diagnostic accuracy of intravascular MR examination in delineating the vessel wall and detecting plaques were analyzed and compared using IVUS.

Results

Compared with IVUS, good agreement was found between MRI and IVUS for delineating vessel, lumen, and vessel wall areas and plaque burden (r value: 0.98, 0.95, 0.96 and 0.91, respectively; P<0.001).

Conclusion

Compared with IVUS, using an intravascular imaging guide-wire to image deep seated arteries allowed determination of the vessel, lumen and vessel wall areas and plaque size and burden. This may provide an alternative method for detecting atherosclerotic plaques in the future.  相似文献   

17.

Objectives

Magnifying narrow-band imaging (M-NBI) is more accurate than white-light imaging for diagnosing small gastric cancers. However, it is uncertain whether moving M-NBI images have additional effects in the diagnosis of gastric cancers compared with still images.

Design

A prospective multicenter cohort study.

Methods

To identify the additional benefits of moving M-NBI images by comparing the diagnostic accuracy of still images only with that of both still and moving images. Still and moving M-NBI images of 40 gastric lesions were obtained by an expert endoscopist prior to this prospective multicenter cohort study. Thirty-four endoscopists from ten different Japanese institutions participated in the prospective multicenter cohort study. Each study participant was first tested using only still M-NBI images (still image test), then tested 1 month later using both still and moving M-NBI images (moving image test). The main outcome was a difference in the diagnostic accuracy of cancerous versus noncancerous lesions between the still image test and the moving image test.

Results

Thirty-four endoscopists were analysed. There were no significant difference of cancerous versus noncancerous lesions between still and moving image tests in the diagnostic accuracy (59.9% versus 61.5%), sensitivity (53.4% versus 55.9%), and specificity (67.0% versus 67.6%). And there were no significant difference in the diagnostic accuracy between still and moving image tests of demarcation line (65.4% versus 65.5%), microvascular pattern (56.7% versus 56.9%), and microsurface pattern (48.1% versus 50.9%). Diagnostic accuracy showed no significant difference between the still and moving image tests in the subgroups of endoscopic findings of the lesions.

Conclusions

The addition of moving M-NBI images to still M-NBI images does not improve the diagnostic accuracy for gastric lesions. It is reasonable to concentrate on taking sharp still M-NBI images during endoscopic observation and use them for diagnosis.

Trial registration

Umin.ac.jp UMIN-CTR000008048  相似文献   

18.

Introduction

In recent years, there has been an increased demand for computer-aided diagnosis (CAD) tools to support clinicians in the field of indirect immunofluorescence. To this aim, academic and industrial research is focusing on detecting antinuclear, anti-neutrophil, and anti-double-stranded (anti-dsDNA) antibodies. Within this framework, we present a CAD system for automatic analysis of dsDNA antibody images using a multi-step classification approach. The final classification of a well is based on the classification of all its images, and each image is classified on the basis of the labeling of its cells.

Methods

We populated a database of 342 images—74 positive (21.6%) and 268 negative (78.4%)— belonging to 63 consecutive sera: 15 positive (23.8%) and 48 negative (76.2%). We assessed system performance by using k-fold cross-validation. Furthermore, we successfully validated the recognition system on 83 consecutive sera, collected by using different equipment in a referral center, counting 279 images: 92 positive (33.0%) and 187 negative (67.0%).

Results

With respect to well classification, the system correctly classified 98.4% of wells (62 out of 63). Integrating information from multiple images of the same wells recovers the possible misclassifications that occurred at the previous steps (cell and image classification). This system, validated in a clinical routine fashion, provides recognition accuracy equal to 100%.

Conclusion

The data obtained show that automation is a viable alternative for Crithidia luciliae immunofluorescence test analysis.  相似文献   

19.

Purpose

To propose a new objective scatter index (OSI) based in the analysis of double-pass images of a point source to rank and classify cataract patients. This classification scheme is compared with a current subjective system.

Methods

We selected a population including a group of normal young eyes as control and patients diagnosed with cataract (grades NO2, NO3 and NO4) according to the Lens Opacities Classification System (LOCS III). For each eye, we recorded double-pass retinal images of a point source. In each patient, we determined an objective scatter index (OSI) as the ratio of the intensity at an eccentric location in the image and the central part. This index provides information on the relevant forward scatter affecting vision. Since the double-pass retinal images are affected by both ocular aberrations and intraocular scattering, an analysis was performed to show the ranges of contributions of aberrations to the OSI.

Results

We used the OSI values to classify each eye according to the degree of scatter. The young normal eyes of the control group had OSI values below 1, while the OSI for subjects in LOCS grade II were around 1 to 2. The use of the objective index showed some of the weakness of subjective classification schemes. In particular, several subjects initially classified independently as grade NO2 or NO3 had similar OSI values, and in some cases even higher than subjects classified as grade NO4. A new classification scheme based in OSI is proposed.

Conclusions

We introduced an objective index based in the analysis of double-pass retinal images to classify cataract patients. The method is robust and fully based in objective measurements; i.e., not depending on subjective decisions. This procedure could be used in combination with standard current methods to improve cataract patient surgery scheduling.  相似文献   

20.

Background

Circulating endothelial cells (CECs) are markers of vascular damage that have clinical relevance in many diseases, including acute myocardial infarction (AMI), and may be predictors of treatment responses. Herein, we investigated the diagnostic and prognostic value of CEC monitoring in AMI patients and a murine model.

Methodology/Principal Findings

CECs were defined as Hoechst 33342+/CD45−/CD31+/CD146+/CD133 in human blood samples and Hoechst 33342+/CD45−/CD31+/KDR+/CD117 in murine samples. To evaluate the validity and variability of our CEC detection system, peripheral blood samples of vascular endothelial growth factor-treated athymic nude mice and AMI patients were collected and subjected to intra-assay analysis. CEC detection by flow cytometry and real-time PCR were compared. Blood samples were obtained from 61 AMI patients, 45 healthy volunteers and 19 samples of the original AMI patients accepted one month treatment, via flow cytometry and expressed as a percentage of peripheral blood mononuclear cells.

Results

Our CEC detection method was validated and had limited variability. CEC concentrations were higher in AMI patients compared to healthy controls. One month post-treatment, CECs levels decreased significantly.

Conclusions/Significance

CEC levels may be useful as a diagnostic and prognostic biomarker in AMI patients.  相似文献   

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

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