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

Purpose

Patient-specific models of anatomical structures and pathologies generated from volumetric medical images play an increasingly central role in many aspects of patient care. A key task in generating these models is the segmentation of anatomical structures and pathologies of interest. Although numerous segmentation methods are available, they often produce erroneous delineations that require time-consuming modifications.

Methods

   We present a new geometry-based algorithm for the reliable detection and correction of segmentation errors in volumetric medical images. The method is applicable to anatomical structures consisting of a few 3D star-shaped components. First, it detects segmentation errors by casting rays from the initial segmentation interior to its outer surface. It then classifies the segmentation surface into correct and erroneous regions by minimizing an energy functional that incorporates first- and second-order properties of the rays lengths. Finally, it corrects the segmentation errors by computing new locations for the erroneous surface points by Laplace deformation so that the new surface has maximum smoothness with respect to the rays-length gradient magnitude.

Results

   Our evaluation on initial segmentations of 16 abdominal aortic aneurysm and 12 lung tumors in CT scans obtained by both adaptive region-growing and active contours level-set segmentation improved the volumetric overlap error by 66 and 70.5 % respectively, with respect to the ground-truth.

Conclusions

   The advantages of our method are that it is independent of the initial segmentation algorithm that covers a variety of anatomical structures and pathologies, that it does not require a shape prior, and that it requires minimal user interaction.
  相似文献   

2.

Purpose

To constrain the risk of severe toxicity in radiotherapy and radiosurgery, precise volume delineation of organs at risk is required. This task is still manually performed, which is time-consuming and prone to observer variability. To address these issues, and as alternative to atlas-based segmentation methods, machine learning techniques, such as support vector machines (SVM), have been recently presented to segment subcortical structures on magnetic resonance images (MRI).

Methods

SVM is proposed to segment the brainstem on MRI in multicenter brain cancer context. A dataset composed by 14 adult brain MRI scans is used to evaluate its performance. In addition to spatial and probabilistic information, five different image intensity values (IIVs) configurations are evaluated as features to train the SVM classifier. Segmentation accuracy is evaluated by computing the Dice similarity coefficient (DSC), absolute volumes difference (AVD) and percentage volume difference between automatic and manual contours.

Results

Mean DSC for all proposed IIVs configurations ranged from 0.89 to 0.90. Mean AVD values were below \(1.5\,\hbox {cm}^{3}\), where the value for best performing IIVs configuration was \(0.85\,\hbox {cm}^{3}\), representing an absolute mean difference of \(3.99\,\%\) with respect to the manual segmented volumes.

Conclusion

Results suggest consistent volume estimation and high spatial similarity with respect to expert delineations. The proposed approach outperformed presented methods to segment the brainstem, not only in volume similarity metrics, but also in segmentation time. Preliminary results showed that the approach might be promising for adoption in clinical use.
  相似文献   

3.

Purpose

In the current standard of care, real-time transrectal ultrasound (TRUS) is commonly used for prostate brachytherapy guidance. As TRUS provides limited soft tissue contrast, segmenting the prostate gland in TRUS images is often challenging and subject to inter-observer and intra-observer variability, especially at the base and apex where the gland boundary is hard to define. Magnetic resonance imaging (MRI) has higher soft tissue contrast allowing the prostate to be contoured easily. In this paper, we aim to show that prostate segmentation in TRUS images informed by MRI priors can improve on prostate segmentation that relies only on TRUS images.

Methods

First, we compare the TRUS-based prostate segmentation used in the treatment of 598 patients with a high-quality MRI prostate atlas and observe inconsistencies at the apex and base. Second, motivated by this finding, we propose an alternative TRUS segmentation technique that is fully automatic and uses MRI priors. The algorithm uses a convolutional neural network to segment the prostate in TRUS images at mid-gland, where the gland boundary can be clearly seen. It then reconstructs the gland boundary at the apex and base with the aid of a statistical shape model built from an MRI atlas of 78 patients.

Results

Compared to the clinical TRUS segmentation, our method achieves similar mid-gland segmentation results in the 598-patient database. For the seven patients who had both TRUS and MRI, our method achieved more accurate segmentation of the base and apex with the MRI segmentation used as ground truth.

Conclusion

Our results suggest that utilizing MRI priors in TRUS prostate segmentation could potentially improve the performance at base and apex.
  相似文献   

4.

Purpose

The aim of the study was to determine if increasing post-therapy calcification in peritoneal metastases in recurrent low-grade serous ovarian carcinomas indicated response to therapy.

Materials and methods

Retrospective analysis of patients with histologically confirmed, recurrent low-grade serous ovarian carcinoma who received treatment at our institution between 2000 and 2014 was performed. Only patients who had calcified tumor implants and showed either interval increase or decrease in tumor calcification following therapy were included in the study. Pre- and post-therapy CT scans of these patients were reviewed by 2 radiologists independently. Changes in the tumor calcification status and tumor deposits size were correlated with serum CA-125 levels. Fisher’s exact test was used to assess the association between peritoneal deposit and calcification status with serum CA-125 status.

Results

35 Patients were included in the study. Based on serial serum CA 125 levels, 22 patients (63%) had progressive disease, 12 (34%) had partial response and 1 (3%) had stable disease. Using RECIST 1.1, 16 had progressive disease, 3 had partial response and 16 had stable disease. In the patients with progressive disease, post-therapy tumor calcification increased in 77% and decreased in 23%. Fischer’s exact test showed that serum CA 125 change was significantly associated with change in size of peritoneal deposits and calcification change.

Conclusions

This preliminary study shows that post-therapy increase in peritoneal implant calcification in low-grade serous ovarian carcinomas is not an indicator of response to therapy.
  相似文献   

5.

Purpose

Simple renal cysts are a common benign finding in abdominal CT scans. However, since they may evolve in time, simple cysts need to be reported. With an ever-growing number of slices per CT scan, cysts are easily overlooked by the overloaded radiologist. In this paper, we address the detection of simple renal cysts as an incidental finding in a real clinical setting.

Methods

We propose a fully automatic framework for renal cyst detection, supported by a robust segmentation of the kidneys performed by a fully convolutional neural network. A combined 3D distance map of the kidneys and surrounding fluids provides initial candidates for cysts. Eventually, a second convolutional neural network classifies the candidates as cysts or non-cyst objects.

Results

Performance was evaluated on 52 abdominal CT scans selected at random in a real radiological workflow and containing over 70 cysts annotated by an experienced radiologist. Setting the minimal cyst diameter to 10 mm, the algorithm detected 59/70 cysts (true-positive rate = 84.3%) while producing an average of 1.6 false-positive per case.

Conclusions

The obtained results suggest the proposed framework is a promising approach for the automatic detection of renal cysts as incidental findings of abdominal CT scans.
  相似文献   

6.

Purpose

A new method for acetabular cartilage segmentation in both computed tomography (CT) arthrography and magnetic resonance imaging (MRI) datasets with leg tension is developed and tested.

Methods

The new segmentation method is based on the combination of shape and intensity information. Shape information is acquired according to the predictable nonlinear relationship between the U-shaped acetabulum region and acetabular cartilage. Intensity information is obtained from the acetabular cartilage region automatically to complete the segmentation procedures. This method is evaluated using 54 CT arthrography datasets with two different radiation doses and 20 MRI datasets. Additionally, the performance of this method in identifying acetabular cartilage is compared with four other acetabular cartilage segmentation methods.

Results

This method performed better than the comparison methods. Indeed, this method maintained good accuracy level for 74 datasets independent of the cartilage modality and with minimum user interaction in the bone segmentation procedures. In addition, this method was efficient in noisy conditions and in detection of the damaged cartilages with zero thickness, which confirmed its potential clinical usefulness.

Conclusions

Our new method proposes acetabular cartilage segmentation in three different datasets based on the combination of the shape and intensity information. This method executes well in situations where there are clear boundaries between the acetabular and femoral cartilages. However, the acetabular cartilage and pelvic bone information should be obtained from one dataset such as CT arthrography or MRI datasets with leg traction.
  相似文献   

7.

Purpose

Abdominal aortic aneurysm (AAA) is a localized, permanent and irreversible enlargement of the artery, with the formation of thrombus into the inner wall of the aneurysm. A precise patient-specific segmentation of the thrombus is useful for both the pre-operative planning to estimate the rupture risk, and for post-operative assessment to monitor the disease evolution. This paper presents a generic approach for 3D segmentation of thrombus from patients suffering from AAA using computed tomography angiography (CTA) scans.

Methods

A fast and versatile thrombus segmentation approach has been developed. It is composed of initial centerline detection and aorta lumen segmentation, an optimized pre-processing stage and the use of a 3D deformable model. The approach has been designed to be very generic and requires minimal user interaction. The proposed method was tested on different datasets with 145 patients overall, including pre- and post-operative CTAs, abdominal aorta and iliac artery sections, different calcification degrees, aneurysm sizes and contrast enhancement qualities.

Results

The thrombus segmentation approach showed very accurate results with respect to manual delineations for all datasets (\(\hbox {Dice} = 0.86 \pm 0.06, 0.81 \pm 0.06\) and \(0.87 \pm 0.03\) for abdominal aorta sections on pre-operative CTA, iliac artery sections on pre-operative CTAs and aorta sections on post-operative CTA, respectively). Experiments on the different patient and image conditions showed that the method was highly versatile, with no significant differences in term of precision. Comparison with the level-set algorithm also demonstrated the superiority of the 3D deformable model. Average processing time was \(8.2 \pm 3.5 \hbox { s}\).

Conclusion

We presented a near-automatic and generic thrombus segmentation algorithm applicable to a large variability of patient and imaging conditions. When integrated in an endovascular planning system, our segmentation algorithm shows its compatibility with clinical routine and could be used for pre-operative planning and post-operative assessment of endovascular procedures.
  相似文献   

8.

Purpose

The aim of this study is to introduce a fully automatic and reproducible short echo-time (STE) magnetic resonance imaging (MRI) segmentation approach for MR-based attenuation correction of positron emission tomography (PET) data in head region.

Procedures

Single STE-MR imaging was followed by generating attenuation correction maps (μ-maps) through exploiting an automated clustering-based level-set segmentation approach to classify head images into three regions of cortical bone, air, and soft tissue. Quantitative assessment was performed by comparing the STE-derived region classes with the corresponding regions extracted from X-ray computed tomography (CT) images.

Results

The proposed segmentation method returned accuracy and specificity values of over 90 % for cortical bone, air, and soft tissue regions. The MR- and CT-derived μ-maps were compared by quantitative histogram analysis.

Conclusions

The results suggest that the proposed automated segmentation approach can reliably discriminate bony structures from the proximal air and soft tissue in single STE-MR images, which is suitable for generating MR-based μ-maps for attenuation correction of PET data.
  相似文献   

9.

Purpose

Brain tumor segmentation is a required step before any radiation treatment or surgery. When performed manually, segmentation is time consuming and prone to human errors. Therefore, there have been significant efforts to automate the process. But, automatic tumor segmentation from MRI data is a particularly challenging task. Tumors have a large diversity in shape and appearance with intensities overlapping the normal brain tissues. In addition, an expanding tumor can also deflect and deform nearby tissue. In our work, we propose an automatic brain tumor segmentation method that addresses these last two difficult problems.

Methods

We use the available MRI modalities (T1, T1c, T2) and their texture characteristics to construct a multidimensional feature set. Then, we extract clusters which provide a compact representation of the essential information in these features. The main idea in this work is to incorporate these clustered features into the 3D variational segmentation framework. In contrast to previous variational approaches, we propose a segmentation method that evolves the contour in a supervised fashion. The segmentation boundary is driven by the learned region statistics in the cluster space. We incorporate prior knowledge about the normal brain tissue appearance during the estimation of these region statistics. In particular, we use a Dirichlet prior that discourages the clusters from the normal brain region to be in the tumor region. This leads to a better disambiguation of the tumor from brain tissue.

Results

We evaluated the performance of our automatic segmentation method on 15 real MRI scans of brain tumor patients, with tumors that are inhomogeneous in appearance, small in size and in proximity to the major structures in the brain. Validation with the expert segmentation labels yielded encouraging results: Jaccard (58%), Precision (81%), Recall (67%), Hausdorff distance (24 mm).

Conclusions

Using priors on the brain/tumor appearance, our proposed automatic 3D variational segmentation method was able to better disambiguate the tumor from the surrounding tissue.
  相似文献   

10.

Purpose

We present a platform, GRAphical Pipeline Environment (GRAPE), to facilitate the development of patient-adaptive magnetic resonance imaging (MRI) protocols.

Methods

GRAPE is an open-source project implemented in the Qt C++ framework to enable graphical creation, execution, and debugging of real-time image analysis algorithms integrated with the MRI scanner. The platform provides the tools and infrastructure to design new algorithms, and build and execute an array of image analysis routines, and provides a mechanism to include existing analysis libraries, all within a graphical environment. The application of GRAPE is demonstrated in multiple MRI applications, and the software is described in detail for both the user and the developer.

Results

GRAPE was successfully used to implement and execute three applications in MRI of the brain, performed on a 3.0-T MRI scanner: (i) a multi-parametric pipeline for segmenting the brain tissue and detecting lesions in multiple sclerosis (MS), (ii) patient-specific optimization of the 3D fluid-attenuated inversion recovery MRI scan parameters to enhance the contrast of brain lesions in MS, and (iii) an algebraic image method for combining two MR images for improved lesion contrast.

Conclusions

GRAPE allows graphical development and execution of image analysis algorithms for inline, real-time, and adaptive MRI applications.
  相似文献   

11.

Purpose

Patient-specific quantitative assessments of muscle mass and biomechanical musculoskeletal simulations require segmentation of the muscles from medical images. The objective of this work is to automate muscle segmentation from CT data of the hip and thigh.

Method

We propose a hierarchical multi-atlas method in which each hierarchy includes spatial normalization using simpler pre-segmented structures in order to reduce the inter-patient variability of more complex target structures.

Results

The proposed hierarchical method was evaluated with 19 muscles from 20 CT images of the hip and thigh using the manual segmentation by expert orthopedic surgeons as ground truth. The average symmetric surface distance was significantly reduced in the proposed method (1.53 mm) in comparison with the conventional method (2.65 mm).

Conclusion

We demonstrated that the proposed hierarchical multi-atlas method improved the accuracy of muscle segmentation from CT images, in which large inter-patient variability and insufficient contrast were involved.
  相似文献   

12.

Background

Resting state magnetic resonance imaging allows studying functionally interconnected brain networks. Here we were aimed to verify functional connectivity between brain networks at rest and its relationship with thalamic microstructure in migraine without aura (MO) patients between attacks.

Methods

Eighteen patients with untreated MO underwent 3 T MRI scans and were compared to a group of 19 healthy volunteers (HV). We used MRI to collect resting state data among two selected resting state networks, identified using group independent component (IC) analysis. Fractional anisotropy (FA) and mean diffusivity (MD) values of bilateral thalami were retrieved from a previous diffusion tensor imaging study on the same subjects and correlated with resting state ICs Z-scores.

Results

In comparison to HV, in MO we found significant reduced functional connectivity between the default mode network and the visuo-spatial system. Both HV and migraine patients selected ICs Z-scores correlated negatively with FA values of the thalamus bilaterally.

Conclusions

The present results are the first evidence supporting the hypothesis that an abnormal resting within networks connectivity associated with significant differences in baseline thalamic microstructure could contribute to interictal migraine pathophysiology.
  相似文献   

13.

Purpose

With the recent trend toward big data analysis, neuroimaging datasets have grown substantially in the past years. While larger datasets potentially offer important insights for medical research, one major bottleneck is the requirement for resources of medical experts needed to validate automatic processing results. To address this issue, the goal of this paper was to assess whether anonymous nonexperts from an online community can perform quality control of MR-based cortical surface delineations derived by an automatic algorithm.

Methods

So-called knowledge workers from an online crowdsourcing platform were asked to annotate errors in automatic cortical surface delineations on 100 central, coronal slices of MR images.

Results

On average, annotations for 100 images were obtained in less than an hour. When using expert annotations as reference, the crowd on average achieves a sensitivity of 82 % and a precision of 42 %. Merging multiple annotations per image significantly improves the sensitivity of the crowd (up to 95 %), but leads to a decrease in precision (as low as 22 %).

Conclusion

Our experiments show that the detection of errors in automatic cortical surface delineations generated by anonymous untrained workers is feasible. Future work will focus on increasing the sensitivity of our method further, such that the error detection tasks can be handled exclusively by the crowd and expert resources can be focused on error correction.
  相似文献   

14.

Purpose

Despite the success of total knee arthroplasty, there continues to be a significant proportion of patients who are dissatisfied. One explanation may be a shape mismatch between pre- and postoperative distal femurs. The purpose of this study was to investigate methods suitable for matching a statistical shape model (SSM) to intraoperatively acquired point cloud data from a surgical navigation system and to validate these against the preoperative magnetic resonance imaging (MRI) data from the same patients.

Methods

A total of 10 patients who underwent navigated total knee arthroplasty also had an MRI scan <2 months preoperatively. The standard surgical protocol was followed which included partial digitization of the distal femur. Two different methods were employed to fit the SSM to the digitized point cloud data, based on (1) iterative closest points and (2) Gaussian mixture models. The available MRI data were manually segmented and the reconstructed three-dimensional surfaces used as ground truth against which the SSM fit was compared.

Results

For both approaches, the difference between the SSM-generated femur and the surface generated from MRI segmentation averaged less than 1.7 mm, with maximum errors occurring in less clinically important areas.

Conclusion

The results demonstrated good correspondence with the distal femoral morphology even in cases of sparse datasets. Application of this technique will allow for measurement of mismatch between pre- and postoperative femurs retrospectively on any case done using the surgical navigation system and could be integrated into the surgical navigation unit to provide real-time feedback.
  相似文献   

15.

Purpose

We aimed to precisely estimate intra-tumoral heterogeneity using spatially regularized spectral clustering (SRSC) on multiparametric MRI data and compare the efficacy of SRSC with the previously reported segmentation techniques in MRI studies.

Procedures

Six NMRI nu/nu mice bearing subcutaneous human glioblastoma U87 MG tumors were scanned using a dedicated small animal 7T magnetic resonance imaging (MRI) scanner. The data consisted of T2 weighted images, apparent diffusion coefficient maps, and pre- and post-contrast T2 and T2* maps. Following each scan, the tumors were excised into 2–3-mm thin slices parallel to the axial field of view and processed for histological staining. The MRI data were segmented using SRSC, K-means, fuzzy C-means, and Gaussian mixture modeling to estimate the fractional population of necrotic, peri-necrotic, and viable regions and validated with the fractional population obtained from histology.

Results

While the aforementioned methods overestimated peri-necrotic and underestimated viable fractions, SRSC accurately predicted the fractional population of all three tumor tissue types and exhibited strong correlations (rnecrotic = 0.92, rperi-necrotic = 0.82 and rviable = 0.98) with the histology.

Conclusions

The precise identification of necrotic, peri-necrotic and viable areas using SRSC may greatly assist in cancer treatment planning and add a new dimension to MRI-guided tumor biopsy procedures.
  相似文献   

16.

Purpose

The aim of this study was to develop an efficient CT scan range estimation method that is based on the analysis of image data itself instead of metadata analysis. This makes it possible to quantitatively compare the scan range of two studies.

Methods

In our study, 3D stacks are first projected to 2D coronal images via a ray casting-like process. Trained 2D body part classifiers are then used to recognize different body parts in the projected image. The detected candidate regions go into a structure grouping process to eliminate false-positive detections. Finally, the scale and position of the patient relative to the projected figure are estimated based on the detected body parts via a structural voting. The start and end lines of the CT scan are projected to a standard human figure. The position readout is normalized so that the bottom of the feet represents 0.0, and the top of the head is 1.0.

Results

Classifiers for 18 body parts were trained using 184 CT scans. The final application was tested on 136 randomly selected heterogeneous CT scans. Ground truth was generated by asking two human observers to mark the start and end positions of each scan on the standard human figure. When compared with the human observers, the mean absolute error of the proposed method is 1.2 % (max: 3.5 %) and 1.6 % (max: 5.4 %) for the start and end positions, respectively.

Conclusion

We proposed a scan range estimation method using multiple body parts detection and relative structure position analysis. In our preliminary tests, the proposed method delivered promising results.
  相似文献   

17.

Purpose

Low-dose CT screening of the lungs is becoming a reality, triggering many more CT-guided lung biopsies. During these biopsies, the patient is submitted to repeated guiding scans with substantial cumulated radiation dose. Extension of the dose reduction to the biopsy procedure is therefore necessary. We propose an image denoising algorithm that specifically addresses the setup of CT-guided lung biopsies. It minimizes radiation exposure while keeping the image quality appropriate for navigation to the target lesion.

Methods

A database of high-SNR CT patches is used to filter noisy pixels in a non-local means framework, while explicitly enforcing local spatial consistency in order to preserve fine image details and structures. The patch database may be created from a multi-patient set of high-SNR lung scans. Alternatively, the first scan, acquired at high-SNR right before the needle insertion, can provide a convenient patient-specific patch database.

Results

The proposed algorithm is compared to state-of-the-art denoising algorithms for a dataset of 43 real CT-guided biopsy scans. Ultra-low-dose scans were simulated by synthetic noise addition to the sinogram, equivalent to a 96% reduction in radiation dose. The feature similarity score for the proposed algorithm outperformed the compared methods for all the scans in the dataset. The benefit of the patient-specific patch database over the multi-patient one is demonstrated in terms of recovered contrast for a tiny porcine lung nodule, following denoising with both approaches.

Conclusions

The proposed method provides a promising approach to the denoising of ultra-low-dose CT-guided biopsy images.
  相似文献   

18.

Purpose

A fully automatic multiatlas-based method for segmentation of the spine and pelvis in a torso CT volume is proposed. A novel landmark-guided diffeomorphic demons algorithm is used to register a given CT image to multiple atlas volumes. This algorithm can utilize both grayscale image information and given landmark coordinate information optimally.

Methods

The segmentation has four steps. Firstly, 170 bony landmarks are detected in the given volume. Using these landmark positions, an atlas selection procedure is performed to reduce the computational cost of the following registration. Then the chosen atlas volumes are registered to the given CT image. Finally, voxelwise label voting is performed to determine the final segmentation result.

Results

The proposed method was evaluated using 50 torso CT datasets as well as the public SpineWeb dataset. As a result, a mean distance error of \(0.59\pm 0.14\hbox { mm}\) and a mean Dice coefficient of \(0.90\pm 0.02\) were achieved for the whole spine and the pelvic bones, which are competitive with other state-of-the-art methods.

Conclusion

From the experimental results, the usefulness of the proposed segmentation method was validated.
  相似文献   

19.

Purpose

The goal of medical content-based image retrieval (M-CBIR) is to assist radiologists in the decision-making process by retrieving medical cases similar to a given image. One of the key interests of radiologists is lesions and their annotations, since the patient treatment depends on the lesion diagnosis. Therefore, a key feature of M-CBIR systems is the retrieval of scans with the most similar lesion annotations. To be of value, M-CBIR systems should be fully automatic to handle large case databases.

Methods

We present a fully automatic end-to-end method for the retrieval of CT scans with similar liver lesion annotations. The input is a database of abdominal CT scans labeled with liver lesions, a query CT scan, and optionally one radiologist-specified lesion annotation of interest. The output is an ordered list of the database CT scans with the most similar liver lesion annotations. The method starts by automatically segmenting the liver in the scan. It then extracts a histogram-based features vector from the segmented region, learns the features’ relative importance, and ranks the database scans according to the relative importance measure. The main advantages of our method are that it fully automates the end-to-end querying process, that it uses simple and efficient techniques that are scalable to large datasets, and that it produces quality retrieval results using an unannotated CT scan.

Results

Our experimental results on 9 CT queries on a dataset of 41 volumetric CT scans from the 2014 Image CLEF Liver Annotation Task yield an average retrieval accuracy (Normalized Discounted Cumulative Gain index) of 0.77 and 0.84 without/with annotation, respectively.

Conclusions

Fully automatic end-to-end retrieval of similar cases based on image information alone, rather that on disease diagnosis, may help radiologists to better diagnose liver lesions.
  相似文献   

20.

Background

Hip arthroscopy has been established over the past decade as a safe treatment of many hip diseases. Complications such as fractures of the femoral neck are very rare.

Case report

We report the case of a 54-year-old woman who developed a stress fracture of the femoral neck 5 weeks after arthroscopic femoral neck osteochondroplasty. To our knowledge this is the first reported case of a female patient.

Conclusions

Postoperative groin pain after load increase should the surgeon aware of a possible femoral neck stress fracture. For further diagnosis, physicians should not hesitate to ask for an MRI or a CT scan because conventional radiographs can be normal. Stable stress fractures can be treated conservatively with restricted weight bearing on 2 crutches for 6 weeks.
  相似文献   

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

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