首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
PurposeTo compare, via Monte Carlo simulations, homogeneous and non-homogenous breast models adopted for mean glandular dose (MGD) estimates in mammography vs. patient specific digital breast phantoms.MethodsWe developed a GEANT4 Monte Carlo code simulating four homogenous cylindrical breast models featured as follows: (1) semi-cylindrical section enveloped in a 5-mm adipose layer; (2) semi-elliptical section with a 4-mm thick skin; (3) semi-cylindrical section with a 1.45-mm skin layer; (4) semi-cylindrical section in a 1.45-mm skin layer and 2-mm subcutaneous adipose layer. Twenty patient specific digital breast phantoms produced from a dedicated CT scanner were assumed as reference in the comparison. We simulated two spectra produced from two anode/filter combinations. An additional digital breast phantom was produced via BreastSimulator software.ResultsWith reference to the results for patient-specific breast phantoms and for W/Al spectra, models #1 and #3 showed higher MGD values by about 1% (ranges [–33%; +28%] and [−31%; +30%], respectively), while for model #4 it was 2% lower (range [−34%; +26%]) and for model #2 –11% (range [−39%; +14%]), on average. On the other hand, for W/Rh spectra, models #1 and #4 showed lower MGD values by 2% and 1%, while for model #2 and #3 it was 14% and 8% lower, respectively (ranges [−43%; +13%] and [−41%; +21%]). The simulation with the digital breast phantom produced with BreastSimulator showed a MGD overestimation of +33%.ConclusionsThe homogeneous breast models led to maximum MGD underestimation and overestimation of 43% and 28%, respectively, when compared to patient specific breast phantoms derived from clinical CT scans.  相似文献   

2.
PurposeTo provide mean glandular dose (MGD) estimates via Monte Carlo (MC) simulations as a function of the breast models and scan parameters in mammography, digital breast tomosynthesis (DBT) and dedicated breast CT (BCT).MethodsThe MC code was based on GEANT4 toolkit. The simulated compressed breast was either a cylinder with a semi-circular section or ad hoc shaped for oblique view (MLO). In DBT we studied the influence of breast models and exam parameters on the T-factors (i.e. the conversion factor for the calculation of the MGD in DBT from that for a 0-degree projection), and in BCT we investigated the influence on the MGD estimates of the ion chamber volume used for the air kerma measurements.ResultsIn mammography, a model representative of a breast undergoing an MLO view exam did not produce substantial differences (0.4%) in MGD estimates, when compared to a conventional cranio-caudal (CC) view breast model. The beam half value layer did not present a significant influence on T-factors in DBT (<0.8%), while the skin model presented significant influence on MGD estimates (up to 3.3% at 30 degrees scan angle), increasing for larger scan angles. We derived a correction factor for taking into account the different ion chamber volume used in MGD estimates in BCT.ConclusionsA series of MC code modules for MGD estimates in 2D and 3D breast imaging have been developed in order to take into account the most recent advances in breast models.  相似文献   

3.
Virtual clinical trials (VCT) are in-silico reproductions of medical examinations, which adopt digital models of patients and simulated devices. They are intended to produce clinically equivalent outcome data avoiding long execution times, ethical issues related to radiation induced risks and huge costs related to real clinical trials with a patient population. In this work, we present a platform for VCT in 2D and 3D X-ray breast imaging. The VCT platform uses Monte Carlo simulations based on the Geant4 toolkit and patient breast models derived from a cohort of high resolution dedicated breast CT (BCT) volume data sets. Projection images of the breast and three-dimensional glandular dose maps are generated for a given breast model, by simulating both 2D full-field digital mammography (DM) and 3D BCT examinations. Uncompressed voxelized breast models were derived from segmented patient images. Compressed versions of the digital breast phantoms for DM were generated using a previously published digital compression algorithm. The Monte Carlo simulation framework has the capability of generating and tracking ~105 photons/s using a server equipped with 16-cores and 3.0 GHz clock speed. The VCT platform will provide a framework for scanner design optimization, comparison between different scanner designs and between different modalities or protocols on computational breast models, without the need for scanning actual patients as in conventional clinical trials.  相似文献   

4.
PurposeThe purpose of this study was to develop and validate a Monte Carlo (MC) simulation tool for patient dose assessment for a 320 detector-row CT scanner, based on the recommendations of International Commission on Radiological Protection (ICRP). Additionally, the simulation was applied on four clinical acquisition protocols, with and without automatic tube current modulation (TCM).MethodsThe MC simulation was based on EGS4 code and was developed specifically for a 320 detector-row cone-beam CT scanner. The ICRP adult reference phantoms were used as patient models. Dose measurements were performed free-in-air and also in four CTDI phantoms: 150 mm and 350 mm long CT head and CT body phantoms. The MC program was validated by comparing simulations results with these actual measurements acquired under the same conditions. The measurements agreed with the simulations across all conditions within 5%. Patient dose assessment was performed for four clinical axial acquisitions using the ICRP adult reference phantoms, one of them using TCM.ResultsThe results were nearly always lower than those obtained from other dose calculator tools or published in other studies, which were obtained using mathematical phantoms in different CT systems. For the protocol with TCM organ doses were reduced by between 28 and 36%, compared to the results obtained using a fixed mA value.ConclusionsThe developed simulation program provides a useful tool for assessing doses in a 320 detector-row cone-beam CT scanner using ICRP adult reference computational phantoms and is ready to be applied to more complex protocols.  相似文献   

5.
PurposeEvaluation of the out-of-field dose is an important aspect in radiotherapy. Due to the fetus radiosensitivity, this evaluation becomes even more conclusive when the patient is pregnant. In this work, a linear accelerator Varian Clinac 2100c operating at 6 MV, a pregnant anthropomorphic phantom (Maria), and different shields added above the abdominal region of the phantom were used for the analysis based on MCNPX. Methods: The simulations were performed for the medial and lateral projections, using either an open field collimation (10×16 cm2) or a multileaf collimator. The added shields (M1 and M2) were designed based on models proposed by Stovall et al. [1], intending to reduce the deposited dose on the fetus and related structures. Results: The presence of the shields showed to be effective in reducing the doses on the fetus, amniotic sac, and placenta, for example. A reduction of about 43% was found in the dose on the fetus when M2 was added, using the open field collimation, in comparison with the situation with no shield, being the lateral projection the main responsible for the dose. The use of MLC significatively reduced the doses in different structures, including on the fetus and amniotic sac, for example, in comparison to the open field situation. A slight increment on the dose in organs such as the eyes, thyroid and brain was found in both collimation systems, due to the presence of the shields. The contribution of the leakage radiation from the tube head of the linear accelerator was found to be in the order of µGy, being reduced by the presence of the M2 shield. Conclusion: Using the shields showed to be an essential feature in order to reduce the dose not only on the fetus, but also in important structures responsible to its development.  相似文献   

6.
PurposeIn cerebral angiography, for diagnosis and interventional neuroradiology, cone-beam computed tomography (CBCT) scan is frequently performed for evaluating brain parenchyma, cerebral hemorrhage, and cerebral infarction. However, the patient’s eye lens is more frequently exposed to excessive doses in these scans than in the previous angiography and interventional neuroradiology (INR) procedures. Hence, radioprotection for the lenses is needed. This study selects the most suitable eye lens protection material for CBCT from among nine materials by evaluating the dose reduction rate and image quality.MethodsTo determine the dose reduction rate, the lens doses were measured using an anthropomorphic head phantom and a real-time dosimeter. For image quality assessment, the artifact index was calculated based on the pixel value and image noise within various regions of interest in a water phantom.ResultsThe protective materials exhibited dose reduction; however, streak artifacts were observed near the materials. The dose reduction rate and the degree of the artifact varied significantly depending on the protective material. The dose reduction rates were 14.6%, 14.2%, and 26.0% when bismuth shield: normal (bismuth shield in the shape of an eye mask), bismuth shield: separate (two separate bismuth shields), and lead goggles were used, respectively. The “separate” bismuth shield was found to be effective in dose reduction without lowering the image quality.ConclusionWe found that bismuth shields and lead goggles are suitable protective devices for the optimal reduction of lens doses.  相似文献   

7.
ObjectiveTo evaluate the effect of cone-beam computed tomography (CBCT) image acquisition protocols on image quality, lesion detection, delineation, and patient dose.Methods100-patients and a CTDI phantom combined with an electron density phantom were examined using four different CBCT-image acquisition protocols during image-guided transarterial chemoembolization (TACE). Protocol-1 (time: 6 s, tube rotation: 360°), protocol-2 (5 s, 300°), protocol-3 (4 s, 240°) and protocol-4 (3 s, 180°) were used. The protocols were first investigated using a phantom. The protocols that were found to be clinically appropriate in terms of image quality and radiation dose were then assessed on patients. A higher radiation dose and/or a poor image quality were inappropriate for the patient imaging. Patient dose (patient-entrance dose and dose-area product), image quality (Hounsfield Unit, noise, signal-to-noise ratio and contrast-to-noise ratio), and lesion delineation (tumor-liver contrast) were assessed and compared using appropriate statistical tests. Lesion detectability, sensitivity, and predictive values were estimated for CBCT-image data using pre-treatment patient magnetic resonance imaging.ResultsThe estimated patient dose showed no statistical significance (p > 0.05) between protocols-2 and -3; the assessed image quality between these protocols manifested insignificant difference (p > 0.05). Two other phantom protocols were not considered for patient imaging due to significantly higher dose (protocols-1) and poor image quality (protocol-4). Lesion delineation and detection were insignificant (p > 0.05) between protocols-2 and -3. Lesion sensitivities generated were 81–89% (protocol-2) and 81–85% (protocol-3) for different lesion types.ConclusionData acquisition using protocols-2 and -3 provided good image quality, lesion detection and delineation with acceptable patient dose during CBCT-imaging mainly due to similar frame numbers acquired.  相似文献   

8.
Rationale and objectivesDedicated breast CT and PET/CT scanners provide detailed 3D anatomical and functional imaging data sets and are currently being investigated for applications in breast cancer management such as diagnosis, monitoring response to therapy and radiation therapy planning. Our objective was to evaluate the performance of the diffeomorphic demons (DD) non-rigid image registration method to spatially align 3D serial (pre- and post-contrast) dedicated breast computed tomography (CT), and longitudinally-acquired dedicated 3D breast CT and positron emission tomography (PET)/CT images.MethodsThe algorithmic parameters of the DD method were optimized for the alignment of dedicated breast CT images using training data and fixed. The performance of the method for image alignment was quantitatively evaluated using three separate data sets; (1) serial breast CT pre- and post-contrast images of 20 women, (2) breast CT images of 20 women acquired before and after repositioning the subject on the scanner, and (3) dedicated breast PET/CT images of 7 women undergoing neo-adjuvant chemotherapy acquired pre-treatment and after 1 cycle of therapy.ResultsThe DD registration method outperformed no registration (p < 0.001) and conventional affine registration (p ≤ 0.002) for serial and longitudinal breast CT and PET/CT image alignment. In spite of the large size of the imaging data, the computational cost of the DD method was found to be reasonable (3–5 min).ConclusionsCo-registration of dedicated breast CT and PET/CT images can be performed rapidly and reliably using the DD method. This is the first study evaluating the DD registration method for the alignment of dedicated breast CT and PET/CT images.  相似文献   

9.
10.
The aim of this work was to research and evaluate the performance of three different digital breast tomosynthesis (DBT) systems in the clinical environment (Siemens Mammomat Inspiration, Hologic Selenia Dimensions, and Fujifilm Amulet Innovality). The characterization included the study of the detector, the automatic exposure control, and the resolution of DBT projections and reconstructed planes.The modulation transfer function (MTF) of the DBT projections was measured with a 1 mm thick steel edge, showing a strong anisotropy (30–40% lower MTF0.5 frequencies in the tube travel direction). The in-plane MTF0.5, measured with a 25 μm tungsten wire, ranges from 1.3 to 1.8 lp/mm in the tube-travel direction and between 2.4 and 3.7 lp/mm in the chest wall–nipple. In the latter direction, the MTF peak shift is more emphasized for large angular range systems (2.0 versus 1.0 lp/mm). In-depth resolution of the planes, via the full width at half maximum (FWHM) from the point spread function of a 25 μm tungsten wire, is not only influenced by angular range and yields 1.3–4.6 mm among systems. The artifact spread function from 1 mm diameter tungsten beads depends mainly on angular range, yielding two tendencies whether large (FWHM is 4.5 mm) or small (FWHM is 10 mm) angular range is used. DBT delivers per scan a mean glandular dose between 1.4 and 2.7 mGy for a 45 mm thick polymethyl methacrylate (PMMA) block.In conclusion, we have identified and analysed specific metrics that can be used for quality assurance of DBT systems.  相似文献   

11.
AimThe main purpose of the present study is assessment of skin dose in breast cancer radiotherapy.BackgroundAccurate assessment of skin dose in radiotherapy can provide useful information for clinical considerations.Materials and methodsA RANDO phantom was irradiated using a 6 MV Siemens Primus linac with medial and tangential radiotherapy fields for simulating breast cancer treatment. Dosimetry was also performed on various positions across the fields using an EBT3 radiochromic film. Similar conditions of measurement on the RANDO phantom including field size, irradiation angle, number of fields, etc. were subsequently simulated via the Monte Carlo N-Particle Transport code (MCNP). Ultimately, dose values for corresponding points from both methods were compared.ResultsConsidering dosimetry using radiochromic films on the RANDO phantom, there were points having underdose and overdose based on the prescribed dose and skin tolerance levels. In this respect, 81.25% and 18.75% of the points had underdose and overdose, respectively. In some cases, several differences were observed between the measurement and the MCNP simulation results associated with skin dose.ConclusionBased on the results of the points which had underdose, it was suggested that a bolus should be used for the given points. With regard to overdose points, it was advocated to consider skin tolerance dose in treatment planning. Differences between the measurement and the MCNP simulation results might be due to voxel size of tally cells in simulations, effect of beam’s angle of incidence, validation time of linac’s head, lack of electronic equilibrium in the build-up region, as well as MCNP tally type.  相似文献   

12.
PurposeThis study reports a sensitivity enhancement of gold-coated contact lens-type ocular in vivo dosimeters (CLODs) for low-dose measurements in computed tomography (CT).MethodsMonte Carlo (MC) simulations were conducted to evaluate the dose enhancement from the gold (Au) layers on the CLODs. The human eye and CLODs were modeled, and the X-ray tube voltages were defined as 80, 120, and 140 kVp. The thickness of the Au layer attached to a CLOD ranged from 100 nm to 10 μm. The thickness of the active layer ranged from 20 to 140 μm. The dose ratio between the active layer of the Au-coated CLOD and a CLOD without a layer, i.e., the dose enhancement factor (DEF), was calculated.ResultsThe DEFs of the first 20-μm thick active layer of the 5-μm thick Au-coated CLOD were 18.4, 19.7, 20.2 at 80, 120, and 140 kVp, respectively. The DEFs decreased as the thickness of the active layer increased. The DEFs of 100-nm to 5-μm thick Au layers increased from 1.7 to 5.4 for 120-kVp X-ray tube voltage when the thickness of the active layer was 140 μm.ConclusionsThe MC results presented a higher sensitivity of Au-coated CLODs (∼20-times higher than that of CLODs without a gold layer). Au-coated CLODs can be applied to an evaluation of very low doses (a few cGy) delivered to patients during CT imaging.  相似文献   

13.
Micro computed tomography (µCT) scanners are used to create high-resolution images and to quantify properties of the scanned objects. While modern µCT scanners benefit from the cone beam geometry, they are compromised by scatter radiation. This work aims to develop a Monte Carlo (MC) model of a µCT scanner in order to characterize the scatter radiation in the detector plane.The EGS++ framework with the MC code EGSnrc was used to simulate the particle transport through the main components of the XtremeCT (SCANCO Medical AG, Switzerland). The developed MC model was based on specific information of the manufacturer and was validated against measurements. The primary and the scatter radiation were analyzed and by implementing a dedicated tracing method, the scatter radiation was subdivided into different scatter components.The comparisons of measured and simulated transmission values for different absorber and filter combinations result in a mean difference of 0.2% ± 1.4%, with a maximal local difference of 3.4%. The reconstructed image of the phantom based on measurements agrees well with the image reconstructed using the MC model. The local contribution of scattered radiation is up to 10% of the total radiation in the detector plane and most of the scattered particles result from interactions in the scanned object. The MC simulations show that scatter radiation contains information about the structure of the object.In conclusion, a MC model for a µCT scanner was successfully validated and applied to analyze the characteristics of the scatter radiation for a µCT scanner.  相似文献   

14.
Obesity is associated with increased morbidity and mortality as well as reduced metrics in quality of life. Both environmental and genetic factors are associated with obesity, though the precise underlying mechanisms that contribute to the disease are currently being delineated. Several small animal models of obesity have been developed and are employed in a variety of studies. A critical component to these experiments involves the collection of regional and/or total animal fat content data under varied conditions. Traditional experimental methods available for measuring fat content in small animal models of obesity include invasive (e.g. ex vivo measurement of fat deposits) and non-invasive (e.g. Dual Energy X-ray Absorptiometry (DEXA), or Magnetic Resonance (MR)) protocols, each of which presents relative trade-offs. Current invasive methods for measuring fat content may provide details for organ and region specific fat distribution, but sacrificing the subjects will preclude longitudinal assessments. Conversely, current non-invasive strategies provide limited details for organ and region specific fat distribution, but enable valuable longitudinal assessment. With the advent of dedicated small animal X-ray computed tomography (CT) systems and customized analytical procedures, both organ and region specific analysis of fat distribution and longitudinal profiling may be possible. Recent reports have validated the use of CT for in vivo longitudinal imaging of adiposity in living mice. Here we provide a modified method that allows for fat/total volume measurement, analysis and visualization utilizing the Carestream Molecular Imaging Albira CT system in conjunction with PMOD and Volview software packages.  相似文献   

15.
16.
PurposeTo estimate organ dose and effective dose for patients for cardiac CT as applied in an international multicenter study (CORE320) with a 320-Detector row CT scanner using Monte Carlo (MC) simulations and voxelized phantoms. The effect of positioning of the arms, off-centering the patient and heart rate on patient dose was analyzed.MethodsA MC code was tailored to simulate the geometry and characteristics of the CT scanner. The phantoms representing the adult reference male and female were implemented according to ICRP 110. Effective dose and organ doses were obtained for CT acquisition protocols for calcium scoring, coronary angiography and myocardial perfusion.ResultsFor low heart rate, the normalized effective dose (E) for cardiac CT was higher for female (5.6 mSv/100 mAs) compared to male (2.2 mSv/100 mAs) due to the contribution of female breast tissue. Averaged E for female and male was 11.3 mSv for the comprehensive cardiac protocol consisting of calcium scoring (1.9 mSv); coronary angiography including rest cardiac perfusion (5.1 mSv) and stress cardiac perfusion (4.3 mSv). These values almost doubled at higher heart rates (20.1 mSv). Excluding the arms increased effective dose by 6–8%, centering the patient showed no significant effect. The k-factor (0.028 mSv/mGy.cm) derived from this study leads to effective doses up to 2–3 times higher than the values obtained using now outdated methodologies.ConclusionMC modeling of cardiac CT examinations on realistic voxelized phantoms allowed us to assess patient doses accurately and we derived k-factors that are well above those published previously.  相似文献   

17.
PurposeTo determine how the skin dose varies in patients receiving radiation treatment for breast cancer in the prone and supine positions.MethodsFifty patients were scanned in the prone and supine positions. A radiation treatment plan was created for the left breast using a 6-MV beam for a prescribed dose of 42.66 Gy in 16 fractions. The dose was calculated using 1- and 2.5-mm calculation grid sizes and the surface dose was compared in both techniques.ResultsThe median gantry angles relative to the skin surface at the central axis were 8 and 52 degrees for treatment in the prone and supine positions, respectively. The mean dose difference between the prone and supine techniques was statistically significant from 3- to 5-mm depth for both grid sizes. For the 1-mm calculation grid size, the doses at 3-, 4-, and 5-mm depths in the prone and supine techniques were 87.80% and 89.10% (P < 0.003), 91.92% and 94.50% (P < 0.00), and 95.30% and 98.20% (P < 0.00), respectively; for the 2.5-mm grid size, the respective doses were 87.10% and 88.59% (P < 0.00), 91.60% and 94.63% (P < 0.00), and 95.10% and 97.80% (P < 0.00), respectively.ConclusionsThis study demonstrates that the prone technique facilitates a relatively lower skin dose than the supine technique. This observation is probably due to the beam angle. The beam is more perpendicular to the skin surface in the prone technique, whereas it is more tangential in the supine technique, which may deliver a higher skin dose. Thus, the dose to the skin should be evaluated in the prone technique, and if desired, the skin dose could be carefully augmented via a bolus or beam spoiler.  相似文献   

18.
PurposeConverting the measurable quantities to patient organ doses in projection radiography is usually based on a standard-sized patient model and a specific radiation quality, which are likely to differ from the real situation. Large inaccuracies can therefore be obtained in organ doses, because organ doses are dependent on the exposure parameters, exposure geometry and patient anatomy. In this study, the effect of radiation quality and patient thickness on the organ dose conversion factors were determined.MethodsIn this study, the posterior–anterior projection radiograph of the thorax was selected in order to determine the effect of radiation quality (tube voltages of 70–130 kV and total filtrations of 3 mmAl to 4 mmAl + 0.2 mmCu) and patient thickness (anterior–posterior thicknesses of 19.4–30.8 cm) on the breast and lung dose conversion factors. For this purpose, Monte Carlo simulation programs ImpactMC and PCXMC were used with computed tomography examination data of adult male and female patients and mathematical hermaphrodite phantoms, respectively.ResultsCompared to the reference beam quality and patient thickness, the relative variation range in organ dose conversion factors was up to 74% for different radiation qualities and 122% for different patient thicknesses.ConclusionsConversion factors should only be used with comprehensive understanding of the exposure conditions, considering the exposure parameters, exposure geometry and patient anatomy they are valid for. This study demonstrates that patient thickness-specific and radiation quality-specific conversion factors are needed in projection radiography.  相似文献   

19.
Objective:  To analyse the correlation between cytomorphological criteria in smears with atypical glandular cells (AGC) or adenocarcinoma in situ (AIS) and human papillomavirus (HPV) reflex test results with different neoplastic histological diagnoses, particularly to distinguish between glandular and squamous neoplasia.
Methods:  A series of 155 women with glandular abnormalities in their conventional cervical smears was included: 106 with AGC, 35 with AGC associated with high-grade squamous intraepithelial lesion (HSIL) and 14 with AIS. Two reviewers evaluated 35 cytomorphological criteria and hybrid capture II (HCII) was performed in all cases. Colposcopy was carried out in all cases and biopsy in 126/155. For statistical purposes, predictive values and odds ratio (OR) were calculated, followed by chi-square automatic interaction detection.
Results:  Histology detected 56 cases of squamous and 17 of glandular intraepithelial or invasive neoplasia. Predictive values of the papillary groups and feathering criteria for glandular neoplasia were, respectively, 80.0% and 73.3%. Feathering was the criterion with the highest OR for distinguishing glandular from squamous neoplasia and also for distinguishing between glandular and non-neoplastic diagnosis. Rosettes and pseudostratified strips did not perform as well. Multivariant Classification and Regression Trees analysis identified feathering as the best criterion for distinguishing between glandular, squamous and non-neoplastic diagnoses regardless of HPV status.
Conclusions:  Feathering was the best criterion for predicting glandular neoplasia.  相似文献   

20.
ObjectiveThis paper aims to provide some practical recommendations to reduce eye lens dose for workers exposed to X-rays in interventional cardiology and radiology and also to propose an eye lens correction factor when lead glasses are used.MethodsMonte Carlo simulations are used to study the variation of eye lens exposure with operator position, height and body orientation with respect to the patient and the X-ray tube. The paper also looks into the efficiency of wraparound lead glasses using simulations. Computation results are compared with experimental measurements performed in Spanish hospitals using eye lens dosemeters as well as with data from available literature.ResultsSimulations showed that left eye exposure is generally higher than the right eye, when the operator stands on the right side of the patient. Operator height can induce a strong dose decrease by up to a factor of 2 for the left eye for 10-cm-taller operators. Body rotation of the operator away from the tube by 45°–60° reduces eye exposure by a factor of 2. The calculation-based correction factor of 0.3 for wraparound type lead glasses was found to agree reasonably well with experimental data.ConclusionsSimple precautions, such as the positioning of the image screen away from the X-ray source, lead to a significant reduction of the eye lens dose. Measurements and simulations performed in this work also show that a general eye lens correction factor of 0.5 can be used when lead glasses are worn regardless of operator position, height and body orientation.  相似文献   

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

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