首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
We present a novel method to generate parallax-free panoramic X-ray images during surgery by enabling the mobile C-arm to rotate around its X-ray source center, relative to the patient’s table. Rotating the mobile C-arm around its X-ray source center is impractical and sometimes impossible due to the mechanical design of mobile C-arm systems. In order to ensure that the C-arm motion is a relative pure rotation around its X-ray source center, we propose to move the table to compensate for the translational part of the motion based on C-arm pose estimation. For this we employ a visual marker pattern and a Camera Augmented Mobile C-arm system that is a standard mobile C-arm augmented by a video camera and mirror construction. We are able to produce a parallax-free panoramic X-ray image independent of the geometric configuration of imaged anatomical structures. Our method does not require a fronto-parallel setup or any overlap between the acquired X-ray images. This generated parallax-free panoramic X-ray image preserves the linear perspective projection property. It also presents a negligible difference (below 2 pixels) in the overlapping area between two consecutive individual X-ray images and has a high visual quality. This promises suitability for intra-operative clinical applications in orthopedic and trauma surgery. The experiments on phantoms and ex-vivo bone structure demonstrate both the functionality and accuracy of the method.  相似文献   

2.

Purpose

The camera-augmented mobile C-arm (CamC) augments any mobile C-arm by a video camera and mirror construction and provides a co-registration of X-ray with video images. The accurate overlay between these images is crucial to high-quality surgical outcomes. In this work, we propose a practical solution that improves the overlay accuracy for any C-arm orientation by: (i) improving the existing CamC calibration, (ii) removing distortion effects, and (iii) accounting for the mechanical sagging of the C-arm gantry due to gravity.

Methods

A planar phantom is constructed and placed at different distances to the image intensifier in order to obtain the optimal homography that co-registers X-ray and video with a minimum error. To alleviate distortion, both X-ray calibration based on equidistant grid model and Zhang’s camera calibration method are implemented for distortion correction. Lastly, the virtual detector plane (VDP) method is adapted and integrated to reduce errors due to the mechanical sagging of the C-arm gantry.

Results

The overlay errors are 0.38±0.06 mm when not correcting for distortion, 0.27±0.06 mm when applying Zhang’s camera calibration, and 0.27±0.05 mm when applying X-ray calibration. Lastly, when taking into account all angular and orbital rotations of the C-arm, as well as correcting for distortion, the overlay errors are 0.53±0.24 mm using VDP and 1.67±1.25 mm excluding VDP.

Conclusion

The augmented reality fluoroscope achieves an accurate video and X-ray overlay when applying the optimal homography calculated from distortion correction using X-ray calibration together with the VDP.  相似文献   

3.
This paper describes a method for tracking the camera motion of a flexible endoscope, in particular a bronchoscope, using epipolar geometry analysis and intensity-based image registration. The method proposed here does not use a positional sensor attached to the endoscope. Instead, it tracks camera motion using real endoscopic (RE) video images obtained at the time of the procedure and X-ray CT images acquired before the endoscopic examination. A virtual endoscope system (VES) is used for generating virtual endoscopic (VE) images. The basic idea of this tracking method is to find the viewpoint and view direction of the VES that maximizes a similarity measure between the VE and RE images. To assist the parameter search process, camera motion is also computed directly from epipolar geometry analysis of the RE video images. The complete method consists of two steps: (a) rough estimation using epipolar geometry analysis and (b) precise estimation using intensity-based image registration. In the rough registration process, the method computes camera motion from optical flow patterns between two consecutive RE video image frames using epipolar geometry analysis. In the image registration stage, we search for the VES viewing parameters that generate the VE image that is most similar to the current RE image. The correlation coefficient and the mean square intensity difference are used for measuring image similarity. The result obtained in the rough estimation process is used for restricting the parameter search area. We applied the method to bronchoscopic video image data from three patients who had chest CT images. The method successfully tracked camera motion for about 600 consecutive frames in the best case. Visual inspection suggests that the tracking is sufficiently accurate for clinical use. Tracking results obtained by performing the method without the epipolar geometry analysis step were substantially worse. Although the method required about 20 s to process one frame, the results demonstrate the potential of image-based tracking for use in an endoscope navigation system.  相似文献   

4.
背景:基于C型臂2D投影的3D模型重建是一种以XRII图像作为基础,经过校正后的运用一定数值函数进行3D模型重建的技术,可在手术过程中提供给手术者丰富的图像信息,方便手术的进行。目的:探讨基于C型臂2D投影的3D模型重建技术诸多方面的问题。方法:由第一作者检索1990/2010PubMed数据库、CNKI系列数据库及万方数据库有关图像引导手术技术、C型臂2D投影图像校正与重建、基于2D图像的3D模型重建以及医学图像配准等方面的文献。结果与结论:基于C型臂2D投影图像的3D模型重建是指以C型臂获取的2D投影图像为基础,实现骨骼3D模型的术中重建。3D重建模型不仅含有更为丰富的骨骼外部形状等解剖结构信息,而且还可包含骨密度及强度等骨骼内部多元有用信息。该技术可分为两条主线:有限角度锥形束X射线摄影合成方法;基于统计可变模型的非刚性配准方法。未来的研究可将该技术与手术导航相关技术进行结合从而建立手术导航系统。  相似文献   

5.

Purpose

Calibration and registration are the first steps for augmented reality and mixed reality applications. In the medical field, the calibration between an RGB-D camera and a C-arm fluoroscope is a new topic which introduces challenges.

Method

A convenient and efficient calibration phantom is designed by combining the traditional calibration object of X-ray images with a checkerboard plane. After the localization of the 2D marker points in the X-ray images and the corresponding 3D points from the RGB-D images, we calculate the projection matrix from the RGB-D sensor coordinates to the X-ray, instead of estimating the extrinsic and intrinsic parameters simultaneously.

Validation

In order to evaluate the effect of every step of our calibration process, we performed five experiments by combining different steps leading to the calibration. We also compared our calibration method to Tsai’s method to evaluate the advancement of our solution. At last, we simulated the process of estimating the rotation movement of the RGB-D camera using MATLAB and demonstrate that calculating the projection matrix can reduce the angle error of the rotation.

Results

A RMS reprojection error of 0.5 mm is achieved using our calibration method which is promising for surgical applications. Our calibration method is more accurate when compared to Tsai’s method. Lastly, the simulation result shows that using a projection matrix has a lower error than using intrinsic and extrinsic parameters in the rotation estimation.

Conclusions

We designed and evaluated a 3D/2D calibration method for the combination of a RGB-D camera and a C-arm fluoroscope.
  相似文献   

6.
C-arm fluoroscopy images are frequently used for qualitative assessment of prostate brachytherapy. Three-dimensional seed reconstruction from C-arm images is necessary for intraoperative dosimetry and quantitative assessment. Seed reconstruction requires accurately known C-arm poses. We propose to measure the C-arm rotation angles and computationally compensate for inevitable C-arm motion to compute the pose. We compensate the translational motions of a C-arm, such as oscillation, sagging and wheel motion using a three-level optimization algorithm and obviate the need for full pose tracking using external trackers or fiducials. We validated our approach on simulated and 100 clinical data sets from 10 patients and gained on average, a seed matching rate of 98.5%, projection error of 0.33 mm (STD = 0.21 mm) and computation time of 19.8 s per patient, which must be considered as clinically excellent results. We also show that without motion compensation the reconstruction is likely to fail.  相似文献   

7.
8.
Purpose  This paper is concerned with the reconstruction of vascular trees from few projections using discrete tomography. However, its computational cost is high and it lacks robustness when the data are inconsistent. We improve robustness by incorporating an intensity-based camera-correction method. The proposed approach is also capable of handling small motion artifacts by modeling them as repositionings of a virtual X-ray camera. We also present a parallel implementation which substantially reduces reconstruction time. Methods  We propose a data-driven reduction of positional inconsistencies by minimizing the reconstruction residual to increase the robustness. Inspired by motion compen-sation algorithms in SPECT imaging, we combine an intensity-based 2D/3D-registration method with itera-tive reconstruction methods. Our objective is the robust vascular-tree reconstruction from positionally inconsistent data. The speed of the reconstruction is substantially increased by a volume-splitting scheme that allows parallel processing. Results  Vascular trees in the liver can be accurately reconstructed from few positionally inconsistent projections using digitally reconstructed radiographs. We have tested the proposed method on synthetic projection data and on objects imaged with a new robotized C-arm. We measured a decrease in the average reconstruction residual of about 13% for real data compared to projection data without preprocessing. Over 4,600 reconstruction experiments were conducted to evaluate the speed-up obtained when employing the volume-splitting scheme. Reconstruction time decreased linearly with increased number of processor-cores, both for real and synthetic data. Conclusions  The proposed method reduces inconsistencies caused by positioning errors and small motion artifacts. No prior segmentation or detection of correspondences between projections is necessary, because all algorithms are intensity-based. As a result, the proposed method allows for robust, high-quality reconstructions, while reducing radiation dose substantially.  相似文献   

9.

Purpose

   Precise localization in bronchoscopy is challenging, particularly for peripheral lesions that cannot be reached by conventional bronchoscopes with a large working channel. Existing navigation methods are hampered by respiratory motion, e.g., in the lower lobes. We present an image-guided approach that considers respiratory motion and can localize instruments.

Methods

   We developed a rigid chest marker containing steel balls visible in X-ray images and a pattern for passive tracking with an optical camera system. An experimental setup to evaluate stereoscopic localization and to mimic chest motion was established in our interventional suite. The marker motion was recorded, and X-ray images were acquired from different angles using a standard C-arm. All coordinates were expressed with respect to the stationary tracking camera. The feasibility of motion-compensated stereoscopic localization was assessed.

Results

   The orientation of the C-arm could be established with a mean error of less than $1^{\circ }$ . Triangulation based on two different X-ray images from different angles resulted in a mean error of 1.8 ( $\pm $ 0.7) mm. A similar result was obtained when the marker was moved between X-ray acquisitions, and the mean error was 1.6 ( $\pm $ 1.4) mm. The latencies were approximately 80 and 380 ms for tracking camera and X-ray imaging, respectively. Stereoscopic localization of a moving target was feasible.

Conclusions

   The system presents a flexible alternative for precise stereoscopic localization of a bronchoscope or instruments using a standard C-arm. We demonstrated the ability to track multiple moving markers and to compensate for respiratory motion.  相似文献   

10.

Purpose

For augmented reality surgical navigation based on C-arm imaging, accuracy of the overlaid augmented reality onto the X-ray image is imperative. However, overlay displacement is generated when a conventional pinhole model describing a geometric relationship of a normal camera is adopted for C-arm calibration. Thus, a modified model for C-arm calibration is proposed to reduce this displacement, which is essential for accurate surgical navigation.

Method

Based on the analysis of displacement pattern generated for three-dimensional objects, we assumed that displacement originated by moving the X-ray source position according to the depth. In the proposed method, X-ray source movement was modeled as variable intrinsic parameters and represented in the pinhole model by replacing the point source with a planar source.

Results

The improvement which represents a reduced displacement was verified by comparing overlay accuracy for augmented reality surgical navigation between the conventional and proposed methods. The proposed method achieved more accurate overlay on the X-ray image in spatial position as well as depth of the object volume.

Conclusion

We validated that intrinsic parameters that describe the source position were dependent on depth for a three-dimensional object and showed that displacement can be reduced and become independent of depth by using the proposed planar source model.
  相似文献   

11.
A photoacoustic tomograph based on optical ultrasound detection is demonstrated, which is capable of high resolution real-time projection imaging and fast three-dimensional (3D) imaging. Snapshots of the pressure field outside the imaged object are taken at defined delay times after photoacoustic excitation by use of a charge coupled device (CCD) camera in combination with an optical phase contrast method. From the obtained wave patterns photoacoustic projection images are reconstructed using a back propagation Fourier domain reconstruction algorithm. Applying the inverse Radon transform to a set of projections recorded over a half rotation of the sample provides 3D photoacoustic tomography images in less than one minute with a resolution below 100 µm. The sensitivity of the device was experimentally determined to be 5.1 kPa over a projection length of 1 mm. In vivo images of the vasculature of a mouse demonstrate the potential of the developed method for biomedical applications.OCIS codes: (110.0110) Imaging systems, (170.5120) Photoacoustic imaging, (170.3880) Medical and biological imaging, (170.6960) Tomography  相似文献   

12.
Minimally invasive interventions are often performed under fluoroscopic guidance. Drawbacks of fluoroscopic guidance are the fact that the presented images are 2D projections and that both the patient and the clinician are exposed to radiation. Image-guided navigation using pre-interventionally acquired 3D MR or CT data is an alternative. However, this often requires invasive anatomical landmark-based, marker-based or surface-based image-to-patient registration. In this paper, a coupling between an image-guided navigation system and an intraoperative C-arm X-ray device with 3D imaging capabilities (3D rotational X-ray (3DRX) system) that enables direct navigation without invasive image-to-patient registration on 3DRX volumes, is described and evaluated. The coupling is established in a one-time preoperative calibration procedure. The individual steps in the registration procedure are explained and evaluated. The acquired navigation accuracy using this coupling is approximately one millimeter.  相似文献   

13.

Purpose

Conventional tracker configurations for surgical navigation carry a variety of limitations, including limited geometric accuracy, line-of-sight obstruction, and mismatch of the view angle with the surgeon??s-eye view. This paper presents the development and characterization of a novel tracker configuration (referred to as ??Tracker-on-C??) intended to address such limitations by incorporating the tracker directly on the gantry of a mobile C-arm for fluoroscopy and cone-beam CT (CBCT).

Methods

A video-based tracker (MicronTracker, Claron Technology Inc., Toronto, ON, Canada) was mounted on the gantry of a prototype mobile isocentric C-arm next to the flat-panel detector. To maintain registration within a dynamically moving reference frame (due to rotation of the C-arm), a reference marker consisting of 6 faces (referred to as a ??hex-face marker??) was developed to give visibility across the full range of C-arm rotation. Three primary functionalities were investigated: surgical tracking, generation of digitally reconstructed radiographs (DRRs) from the perspective of a tracked tool or the current C-arm angle, and augmentation of the tracker video scene with image, DRR, and planning data. Target registration error (TRE) was measured in comparison with the same tracker implemented in a conventional in-room configuration. Graphics processing unit (GPU)-accelerated DRRs were generated in real time as an assistant to C-arm positioning (i.e., positioning the C-arm such that target anatomy is in the field-of-view (FOV)), radiographic search (i.e., a virtual X-ray projection preview of target anatomy without X-ray exposure), and localization (i.e., visualizing the location of the surgical target or planning data). Video augmentation included superimposing tracker data, the X-ray FOV, DRRs, planning data, preoperative images, and/or intraoperative CBCT onto the video scene. Geometric accuracy was quantitatively evaluated in each case, and qualitative assessment of clinical feasibility was analyzed by an experienced and fellowship-trained orthopedic spine surgeon within a clinically realistic surgical setup of the Tracker-on-C.

Results

The Tracker-on-C configuration demonstrated improved TRE (0.87 ± 0.25)?mm in comparison with a conventional in-room tracker setup (1.92 ± 0.71)?mm (p Conclusions The proposed tracker configuration demonstrated sub-?mm TRE from the dynamic reference frame of a rotational C-arm through the use of the multi-face reference marker. Real-time DRRs and video augmentation from a natural perspective over the operating table assisted C-arm setup, simplified radiographic search and localization, and reduced fluoroscopy time. Incorporation of the proposed tracker configuration with C-arm CBCT guidance has the potential to simplify intraoperative registration, improve geometric accuracy, enhance visualization, and reduce radiation exposure.  相似文献   

14.
Three dimensional quantification of organ shape and structure during minimally invasive surgery (MIS) could enhance precision by allowing the registration of multi-modal or pre-operative image data (US/MRI/CT) with the live optical image. Structured illumination is one technique to obtain 3D information through the projection of a known pattern onto the tissue, although currently these systems tend to be used only for macroscopic imaging or open procedures rather than in endoscopy. To account for occlusions, where a projected feature may be hidden from view and/or confused with a neighboring point, a flexible multispectral structured illumination probe has been developed that labels each projected point with a specific wavelength using a supercontinuum laser. When imaged by a standard endoscope camera they can then be segmented using their RGB values, and their 3D coordinates calculated after camera calibration. The probe itself is sufficiently small (1.7 mm diameter) to allow it to be used in the biopsy channel of commonly used medical endoscopes. Surgical robots could therefore also employ this technology to solve navigation and visualization problems in MIS, and help to develop advanced surgical procedures such as natural orifice translumenal endoscopic surgery.  相似文献   

15.
Cardiologists use two-dimensional projection images in conventional X-ray coronary angiography for the assessment of three-dimensional structures. During minimally invasive interventions there is a need to clearly visualize and analyze contrast filled coronary arteries, surrounding tissue, and implanted devices. Three-dimensional reconstruction of these structures is challenging due to the cardiac and respiratory motion. In this paper we describe a method to automatically generate motion compensated reconstructions of various structures using rotational X-ray angiography.The method uses markers on a device or guide wire to identify and estimate the motion of an object or region of interest in order to register and motion compensate the projection images to generate a motion compensated reconstruction. The method is evaluated on 20 rotational acquisitions and the average marker couple detection rate is 84% for cardiac stents, 90% for closure devices and 20% for contrast filled coronaries. The projection images are motion compensated based on the semi-automatically detected markers and subsequently used for reconstruction. We conclude that it is feasible to reconstruct cardiac stents, closure devices, contrast filled coronaries, and calcified plaques using rotational X-ray angiography.  相似文献   

16.
This paper presents a new hybrid camera motion tracking method for bronchoscopic navigation combining SIFT, epipolar geometry analysis, Kalman filtering, and image registration. In a thorough evaluation, we compare it to state-of-the-art tracking methods. Our hybrid algorithm for predicting bronchoscope motion uses SIFT features and epipolar constraints to obtain an estimate for inter-frame pose displacements and Kalman filtering to find an estimate for the magnitude of the motion. We then execute bronchoscope tracking by performing image registration initialized by these estimates. This procedure registers the actual bronchoscopic video and the virtual camera images generated from 3D chest CT data taken prior to bronchoscopic examination for continuous bronchoscopic navigation. A comparative assessment of our new method and the state-of-the-art methods is performed on actual patient data and phantom data. Experimental results from both datasets demonstrate a significant performance boost of navigation using our new method. Our hybrid method is a promising means for bronchoscope tracking, and outperforms other methods based solely on Kalman filtering or image features and image registration.  相似文献   

17.
In many cases, radio-frequency catheter ablation of the pulmonary veins attached to the left atrium still involves fluoroscopic image guidance. Two-dimensional X-ray navigation may also take advantage of overlay images derived from static pre-operative 3D volumetric data to add anatomical details otherwise not visible under X-ray. Unfortunately, respiratory motion may impair the utility of static overlay images for catheter navigation. We developed a novel approach for image-based 3D motion estimation and compensation as a solution to this problem. It is based on 3D catheter tracking which, in turn, relies on 2D/3D registration. To this end, a bi-plane C-arm system is used to take X-ray images of a special circumferential mapping catheter from two directions. In the first step of the method, a 3D model of the device is reconstructed. Three-dimensional respiratory motion at the site of ablation is then estimated by tracking the reconstructed catheter model in 3D based on bi-plane fluoroscopy. Phantom data and clinical data were used to assess model-based catheter tracking. Our phantom experiments yielded an average 2D tracking error of 1.4 mm and an average 3D tracking error of 1.1 mm. Our evaluation of clinical data sets comprised 469 bi-plane fluoroscopy frames (938 monoplane fluoroscopy frames). We observed an average 2D tracking error of 1.0 ± 0.4 mm and an average 3D tracking error of 0.8 ± 0.5 mm. These results demonstrate that model-based motion-compensation based on 2D/3D registration is both feasible and accurate.  相似文献   

18.

Purpose

 X-ray fluoroscopy guidance is frequently used in medical interventions. Image-guided interventional procedures that employ localization for registration require accurate information about the C-arm’s rotation angle that provides the data externally in real time. Optical, electromagnetic, and image-based pose tracking systems have limited convenience and accuracy. An alternative method to recover C-arm orientation was developed using an accelerometer as tilt sensor.

Methods

    The fluoroscopic C-arm’s orientation was estimated using a tri-axial acceleration sensor mounted on the X-ray detector as a tilt sensor. When the C-arm is stationary, the measured acceleration direction corresponds to the gravitational force direction. The accelerometer was calibrated with respect to the C-arm’s rotation along its two axes, using a high-accuracy optical tracker as a reference. The scaling and offset error of the sensor was compensated using polynomial fitting. The system was evaluated on a GE OEC 9800 C-arm. Results obtained by accelerometer, built-in sensor, and image-based tracking were compared, using optical tracking as ground truth data.

Results

The accelerometer-based orientation measurement error for primary angle rotation was $-0.1\pm 0.0^{\circ }$ and for secondary angle rotation it was $0.1\pm 0.0^{\circ }$ . The built-in sensor orientation measurement error for primary angle rotation was $-0.1\pm 0.2^{\circ }$ , and for secondary angle rotation it was $0.1\pm 0.2^{\circ }$ . The image-based orientation measurement error for primary angle rotation was $-0.1\pm 1.3^{\circ }$ , and for secondary angle rotation it was $-1.3\pm 0.3^{\circ }$ .

Conclusion

The accelerometer provided better results than the built-in sensor and image-based tracking. The accelerometer sensor is small, inexpensive, covers the full rotation range of the C-arm, does not require line of sight, and can be easily installed to any mobile X-ray machine. Therefore, accelerometer tilt sensing is a very promising applicant for orientation angle tracking of C-arm fluoroscopes.  相似文献   

19.
背景:CT成像质量的优劣不仅取决于仪器的精密性和先进性,在很大程度上也取决于重建算法,由二维扇束扫描向三维锥束扫描是CT技术的发展方向,因此,寻找一种合适的锥束重建算法具有无法忽略的意义。目的:探讨基于C型臂超短扫描路径锥束投影的图像合成,为实现基于C型臂2D投影图像的3D模型重建提供算法支持。方法:由第一作者于2012年3至5月检索PubMed数据库、CNKI系列数据库及万方数据库1990年至2011年文献。检索词为"C型臂,超短扫描路径,FDK算法,有限角锥形束三维重建,超短扫描扇束重建算法",检索文章的语言种类为中文和英文。计算机初检得到58篇文献,其中19篇符合纳入标准被保留。结果与结论:基于C型臂2D投影图像的3D模型重建必须进行三维模型的重建,目前应用最为广泛的的三维图像重建方法仍然是FDK。但是FDK算法适用于全路径,对超短路径而言不能直接采用,而通过将二维扇束重建算法推广到三维空间中而获得的短扫描轨迹的FDK类型锥束重建算法可对采集到的锥束投影数据进行感兴趣区域重建。未来的研究可针对减少噪声等干扰数据对重建质量造成的影响进行探讨。  相似文献   

20.
OBJECTIVE: The objective of this study is to illustrate the low accuracy of two-dimensional (2-D) X-ray projection methods for the quantification of the three-dimensional (3-D) shoulder motions. BACKGROUND: The traditional method for the quantification of the gleno-humeral motion is by means of 2-D X-ray recording. The motion was characterized by the scapulo-humeral rhythm: the ratio of the nett humeral elevation over nett scapular rotation. The method was based on the quantification of the planar projection of the spatial positions of X-ray dense structures of the scapula. The deformations introduced by the central projection method, a feature of X-ray projection, cannot be compensated for by calibration: the position of the scapula with respect to the camera setting is unknown, and skeletal landmarks of the scapula cannot uniquely be identified. The transformation from 3-D orientations to 2-D angles will, therefore, be inaccurate. METHODS: A 2-D X-ray projection of the scapula during a typical arm abduction was simulated. The 3-D motion was obtained by means of palpation and subsequent digitization of skeletal landmarks of the scapula. The 3-D positions of the recorded landmarks were projected on a plane by a simulation based on the parameters of the X-ray equipment. The scapulo-humeral rhythm was calculated for the different scapular landmarks, and for the orientation of the subject with respect to the projection axis. The results were compared with previous published scapulo-humeral rhythms. RESULTS: The scapulo-humeral rhythm depends both on the choice of the skeletal landmarks, used to quantify the scapular rotations, and on the orientation of the subject in the X-ray setting. The full range of results obtained from earlier published experiments could be obtained from a simulation based on a single 3-D arm abduction. CONCLUSIONS: The 2-D scapulo-humeral rhythm, obtained from planar X-ray projection, is an inaccurate parameter to define the scapular motions. RELEVANCE: The 2-D scapulo-humeral rhythm is an insensitive parameter to identify clinical disorders in the gleno-humeral motions, 3-D motion recording should be applied. Only when stringent precautions are taken with respect to repeatability of positioning of the subjects, can the method be used to study intra-individual effects, e.g., the follow-up of patients during treatment.  相似文献   

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

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