首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 140 毫秒
1.
Ultrasound B-mode visualization of lesions produced in soft tissues using high-intensity focused ultrasound (HIFU) has been shown to be challenging when there is no cavitation activity and, therefore, no hyperechogenecity in the focal region. We investigated a method for the visualization and localization of HIFU-induced lesions after HIFU delivery was complete based on the change in backscattered radio-frequency (RF) signals. A HIFU transducer was used with focal dimension of 8 mm by 2 mm working at 5 MHz. HIFU was applied at different intensities to produce lesions in ex vivo chicken breast, with or without the generation of hyperecho in B-mode images. We compared lesion locations obtained from our RF-processing method, from measurement of physical lesions after exposure and from the B-mode images, if exposures had resulted in hyperecho. The results showed that the RF amplitude decreased as a function of time immediately after stopping the HIFU exposure. The lesions were clearly visualized in two-dimensional (2-D) images of the decay rate of RF amplitude, no matter with or without hyperecho. In experiments with hyperecho, when comparing to physical lesion locations, there was no statistically significant difference in the localization accuracy between the RF-based and the hyperecho-based method (p = 0.76). In cases without hyperecho, the distance between RF-based locations and measured lesion locations was 3.37 ± 1.59 mm (mean ± standard deviation). The axial and lateral difference were 2.00 ± 2.31 mm and 0.85 ± 2.15 mm, respectively, and no statistically significant difference was found between lesion coordinates (axial: p = 0.37 and lateral: p = 0.15). We demonstrated the feasibility of our proposed RF-based method for the localization of HIFU-induced lesions immediately after HIFU treatment. Using the decay rate in RF amplitude as the signature of lesion formation, our method can detect lesion locations even without the appearance of hyperecho. (E-mail: xlzheng@u.washington.edu)  相似文献   

2.
Adiposity accumulation in the liver is an early-stage indicator of non-alcoholic fatty liver disease. Analysis of ultrasound (US) backscatter echoes from liver parenchyma with deep learning (DL) may offer an affordable alternative for hepatic steatosis staging. The aim of this work was to compare DL classification scores for liver steatosis using different data representations constructed from raw US data. Steatosis in N = 31 patients with confirmed or suspected non-alcoholic fatty liver disease was stratified based on fat-fraction cutoff values using magnetic resonance imaging as a reference standard. US radiofrequency (RF) frames (raw data) and clinical B-mode images were acquired. Intermediate image formation stages were modeled from RF data. Power spectrum representations and phase representations were also calculated. Co-registered patches were used to independently train 1-, 2- and 3-D convolutional neural networks (CNNs), and classifications scores were compared with cross-validation. There were 67,800 patches available for 2-D/3-D classification and 1,830,600 patches for 1-D classification. The results were also compared with radiologist B-mode annotations and quantitative ultrasound (QUS) metrics. Patch classification scores (area under the receiver operating characteristic curve [AUROC]) revealed significant reductions along successive stages of the image formation process (p < 0.001). Patient AUROCs were 0.994 for RF data and 0.938 for clinical B-mode images. For all image formation stages, 2-D CNNs revealed higher patch and patient AUROCs than 1-D CNNs. CNNs trained with power spectrum representations converged faster than those trained with RF data. Phase information, which is usually discarded in the image formation process, provided a patient AUROC of 0.988. DL models trained with RF and power spectrum data (AUROC = 0.998) provided higher scores than conventional QUS metrics and multiparametric combinations thereof (AUROC = 0.986). Radiologist annotations indicated lower hepatic steatosis classification accuracies (Acc = 0.914) with respect to magnetic resonance imaging proton density fat fraction that DL models (Acc = 0.989). Access to raw ultrasound data combined with artificial intelligence techniques may offer superior opportunities for quantitative tissue diagnostics than conventional sonographic images.  相似文献   

3.
This paper simultaneously investigated the transient characteristics of integrated backscatter (IBS), attenuation coefficient and bubble activities as time traces before, during and after HIFU treatment, with different HIFU parameters (acoustic power and duty cycle) in both transparent tissue-mimicking phantoms and freshly excised bovine livers. These dynamic changes of acoustic parameters and bubble activities were correlated with the visualization of lesion development selected from photos, conventional B-mode ultrasound images and differential IBS images over the whole procedure of HIFU treatment. Two-dimensional radiofrequency (RF) data were acquired by a modified diagnostic ultrasound scanner to estimate the changes of mean IBS and attenuation coefficient averaged in the lesion region, and to construct the differential IBS images and B-mode ultrasound images simultaneously. Bubble activities over the whole procedure of HIFU treatment were investigated by the passive cavitation detection (PCD) method and the changes in subharmonic and broadband noise were correlated with the transient characteristics of IBS and attenuation coefficient. When HIFU was switched on, IBS and attenuation coefficient increased with the appearance of bubble clouds in the B-mode and differential IBS image. At the same time, the level of subharmonic and broadband noise rose abruptly. Then, there was an initial decrease in the attenuation coefficient, followed by an increase when at lower HIFU power. As the lesion appeared, IBS and attenuation coefficient both increased rapidly to a value twice that of normal. Then the changes in IBS and attenuation coefficient showed more complex patterns, but still showed a slower trend of increases with lesion development. Violent bubble activities were visible in the gel and were evident as strongly echogenic regions in the differential IBS images and B-mode images simultaneously. This was detected by a dramatic high level of subharmonic and broadband noise at the same time. These bubble activities caused fluctuations in IBS and attenuation coefficient during HIFU treatment. After HIFU, IBS and attenuation coefficient decreased gradually accompanied by the fadeout of bright hyperechoic spot in the B-mode and differential IBS image, but were still higher than normal when they were stable. The increases of IBS and attenuation coefficient were greater when using higher acoustic power or a higher duty cycle of the therapeutic emission. These experiments indicated that the bubble activities had the dominant effects on the transient characteristics of IBS and attenuation. This should be taken into consideration when using the dynamic acoustic-property changes for the potentially real-time monitoring imaging of HIFU treatment. (E-mail: mxwan@mail.xjtu.edu.cn)  相似文献   

4.
This paper presents an algorithm for segmenting left ventricular endocardial boundaries from RF ultrasound. Our method incorporates a computationally efficient linear predictor that exploits short-term spatio-temporal coherence in the RF data. Segmentation is achieved jointly using an independent identically distributed (i.i.d.) spatial model for RF intensity and a multiframe conditional model that relates neighboring frames in the image sequence. Segmentation using the RF data overcomes challenges due to image inhomogeneities often amplified in B-mode segmentation and provides geometric constraints for RF phase-based speckle tracking. The incorporation of multiple frames in the conditional model significantly increases the robustness and accuracy of the algorithm. Results are generated using between 2 and 5 frames of RF data for each segmentation and are validated by comparison with manual tracings and automated B-mode boundary detection using standard (Chan and Vese-based) level sets on echocardiographic images from 27 3D sequences acquired from six canine studies.  相似文献   

5.
In ultrasound-guided high-intensity focused ultrasound (HIFU) therapy, the changes observed on tissue are subtle during treatment; some ultrasound-guided HIFU protocols rely on the observation of significant brightness changes as the indicator of tissue lesions. The occurrence of a distinct hyperechogenic region (“bright-up”) around the focus is often associated with acoustic cavitation resulting in microbubble formation, but it may indicate different physical events such as larger bubbles from boiling (known to alter acoustic impedance) or sometimes lesion formation. A reliable method to distinguish and spatially localize these causes within the tissue would assist the control of HIFU delivery, which is the subject of this paper. Spectral analysis of the radio frequency (RF) signal underlying the B-mode image provides more information on the physical cause, but the usual techniques that are methods on the Fourier transform require a long series for good spectral resolution and so they give poor spatial resolution. This paper introduces an active spectral cavitation detection method to attain high spatial resolution (0.15 × 0.15 mm per pixel) through a parametric statistical method (ARMA modeling) used on finite-length data sets, which enables local changes to be identified more easily. This technique uses the characteristics of the signal itself to optimize the model parameters and structure. Its performance is assessed using synthesized cavitation RF data, and it is then demonstrated in ex vivo bovine liver during and after HIFU exposure. The results suggest that good spatial and spectral resolution can be obtained by the design of suitable algorithms. In ultrasound-guided HIFU, the technique provides a useful supplement to B-mode analysis, with no additional time penalty in data acquisition.  相似文献   

6.
High-intensity focused ultrasound (HIFU) is a promising method used to treat cardiac arrhythmias, as it can induce lesions at a distance throughout myocardium thickness. Numerical modeling is commonly used for ultrasound probe development and optimization of HIFU treatment strategies. This study was aimed at describing a numerical method to simulate HIFU thermal ablation in elastic and mobile heart models. The ultrasound pressure field is computed on a 3-D orthonormal grid using the Rayleigh integral method, and the attenuation is calculated step by step between cells. The temperature distribution is obtained by resolution of the bioheat transfer equation on a 3-D non-orthogonally structured curvilinear grid using the finite-volume method. The simulation method is applied on two regions of the heart (atrioventricular node and ventricular apex) to compare the thermal effects of HIFU ablation depending on deformation, motion type and amplitude. The atrioventricular node requires longer sonication than the ventricular apex to reach the same lesion volume. Motion considerably influences treatment duration, lesion shape and distribution in cardiac HIFU treatment. These results emphasize the importance of considering local motion and deformation in numerical studies to define efficient and accurate treatment strategies.  相似文献   

7.
A whole mount histology protocol for 3-D tissue reconstruction to compare the size and spatial location of tumors (and other components) identified in histology data with that from 3-D ultrasound (US) images is presented. Prostate specimens are imaged in 3-D using B-mode (US) and sonoelastography. The prostate surface is outlined in each B-mode image and a 3-D surface reconstruction is made. The specimen is then prepared for whole mount histology and the histology slides are digitally reconstructed to make a 3-D surface. These two surfaces are then aligned using a 3-D correlation algorithm, and the tumor boundary determined by the pathologist is compared with that using sonoelastography. 3-D images showing the overlapping histology and sonoelastography of prostate surface reconstructions for one prostate are presented to illustrate the technique; results for four prostates yielded an accuracy of 92% +/- 3%.  相似文献   

8.
Model-based elastography is an emerging technique with clinical applications in imaging vascular tissues, guiding minimally invasive therapies and diagnosing breast and prostate cancers. Its usage is limited because ultrasound can measure only the axial component of displacement with high precision. The goal of this study was to assess the effect of lateral sampling frequency, lateral beam-width and the number of active transmission elements on the quality of axial and lateral strain elastograms. Elastographic imaging was performed on gelatin-based phantoms with a modified commercial ultrasound scanner. Three groups of radio-frequency (RF) echo frames were reconstructed from fully synthetic aperture data. In the first group, all 128 transmission elements (corresponding to a lateral beamwidth of 0.22 mm at the center of the field of view) were used to reconstruct RF echo frames with A-line densities that varied from 6.4 lines/mm to 51.2 lines/mm. In the second group, the size of the aperture was varied to produce RF echo frames with lateral beamwidths ranging from 0.22 mm to 0.43 mm and a fixed A-line density of 25.6 lines/mm. In the third group, sparse arrays with varying number of active transmission elements (from 2 to 128) were used to reconstruct RF echo frames, whose A-line density and lateral beamwidth were fixed to 25.6 lines/mm and 0.22 mm, respectively. Applying a two-dimensional (2-D) displacement estimator to the pre- and post-deformed RF echo frames produced displacement elastograms. Axial and lateral strain elastograms were computed from displacement elastograms with a least squares strain estimator. The quality of axial and lateral strain elastograms improved with increasing applied strain and A-line density but decreased with increasing lateral beamwidth and deteriorated as the number of active transmission elements in the sparse arrays were reduced. This work demonstrated that the variance incurred when estimating the lateral component of displacement was reduced considerably when elastography was performed with a synthetic aperture ultrasound imaging system. Satisfactory axial and lateral strain elastograms were produced using a sparse array with as few as 16 active transmission elements.  相似文献   

9.
高强度聚焦超声原位灭活胸腹壁转移瘤的疗效观察   总被引:3,自引:0,他引:3  
目的 研究高强度聚焦超声 (high intensity focused ultrasound,HIFU)无侵袭性原位灭活胸腹壁转移瘤的疗效。方法  18例晚期恶性肿瘤胸腹壁转移患者接受 HIFU治疗 ,观察 HIFU治疗前和治疗后的临床症状和影像学 (B超、 MRI)变化情况。结果  HIFU治疗后 12例患者转移瘤疼痛不适症状完全消失 ,B超和 MRI显示 HIFU治疗转移瘤的靶区完全凝固性坏死 ,13例转移瘤完全坏死 ,5例大部分坏死 (范围 >80 % ) ,HIFU治疗的总有效率 10 0 %。结论  HIFU对胸腹壁转移瘤无侵袭性原位灭活的疗效是肯定和显著的  相似文献   

10.
Observing the morphology of human skin is important in the diagnosis of skin cancer and inflammation and in the assessment of skin aging. High-frequency ultrasound imaging provides high spatial resolution of the deep layers of the skin, which cannot be visualized by optical methods. The objectives of the present study were to develop a three-dimensional (3-D) ultrasound microscope and to observe the morphology of normal human skin in vivo. A concave polyvinylidene fluoride transducer with a central frequency of 120 MHz was excited using an electric pulse generated by semiconductor switching. The transducer was scanned two-dimensionally by using two linear motors on the region-of-interest and the ultrasonic reflection was digitized with 2-GHz sampling. Consecutive B-mode images perpendicular to the skin surface were reconstructed to generate multiplanar reconstructed images and 3-D volume-rendering images clearly showing microstructures such as sebaceous glands and hair follicles. The 3-D ultrasound microscope could be used to successfully image the morphology of human skin noninvasively and may provide important information on skin structure.  相似文献   

11.
We investigated the feasibility of two monitoring imaging methods to visualize and evaluate the high intensity focused ultrasound (HIFU) induced lesions in vitro during and after their formation, which were based on differential ultrasonic parameter estimation. Firstly, ultrasonic attenuation slope of tissue sample was estimated based on the spectral analysis of ultrasound RF backscattered signals. The differential attenuation slope maps were acquired, which were interpreted as the differences between the pretreatment image and those obtained in different stages during HIFU therapy. Secondly, ultrasonic integrated backscatter (IBS), defined as the frequency average of the backscatter transfer function over the useful bandwidth, was proposed quantitatively to evaluate the extent of lesions with the same RF signals as the first method. Differential IBS maps were also acquired to visualize temporal evolution of lesion formation. It was found in pig liver in vitro that more precise definition of the treated area was obtained from the differential IBS images than from differential attenuation slope images. Dramatic increase in both attenuation and IBS value was observed during the therapy, which may be related to dramatic enhancement of cavitation due to boiling and accompanying tissue damage. Two methods to obtain one differential image were compared and the cumulative differential image was found to be able to eliminate noises and artifacts to some extent, which was the cumulation of a series of differential images acquired from the differences between the temporally adjacent RF data frames. Moreover, we presented a bidirectional color code for identification of the artifacts due to tissue movements caused by HIFU radiation force. We conclude that cumulative differential IBS images have the potential to monitor the formation of HIFU-induced lesions.  相似文献   

12.
目的比较增强MRI与超声造影在高强度聚焦超声(HIFU)治疗子宫肌瘤疗效观察中的应用价值。方法 选取2019年2月~2020年2月我院收治102例行HIFU治疗的子宫肌瘤患者,患者治疗前后均进行超声造影和增强MRI检查,比较两种检测方法在疗效评价中的作用。结果两种检测方法治疗前后肌瘤均径、肌瘤体积比较差异无统计学意义(P>0.05),其中增强MRI可明显显示肌瘤的大小和边界,而超声造影对子宫肌瘤的边界图像显示相对不足,且增强MRI图像质量评分、消融率高于超声造影, 差异有统计学意义(P<0.05);增强MRI和超声造影对消融区、残余灶均能有效显示,两者评价HIFU治疗子宫肌瘤总有效率比 较差异无统计学意义(P>0.05)。结论增强MRI在子宫肌瘤定性诊断方面优于超声造影,而超声造影检查更具有可重复性、灵活性,两者均为HIFU治疗前后的必要检查。   相似文献   

13.
In this paper, a new positioning system is proposed for the 3-D ultrasound (US). This system combines the image registration technique and speckle decorrelation algorithm to accurately position sequential ultrasonic images without any additional positioning hardware. The speckle decorrelation algorithm estimates the relative distance of two neighboring frames and the image registration technique gets the range of the whole 3-D ultrasonic data set and makes slight modification on each frame's position. The image registration technique is based on the reference image, which is perpendicular to the 3-D ultrasonic data set. This reference image intersects each frame of the 3-D ultrasonic data set in a line. For each frame, the intersectional line is first found and then the location in the reference image can be used to estimate the position of this frame. This system uses the data set of consecutive 2-D freehand-scanned US B-mode images to construct the 3-D US volume data, and it can be integrated into the 3-D US volume rendering system.  相似文献   

14.
High-intensity focused ultrasound (HIFU) and conventional B-mode ultrasound (US) imaging were synchronized to develop a system for real-time visualization of HIFU treatment. The system was tested in vivo in pig liver. The HIFU application resulted in the appearance of a hyperechoic spot at the focus that faded gradually after cessation of HIFU exposure. The duration of HIFU exposure needed for a hyperechoic spot to appear, was inversely related to the HIFU intensity. The threshold intensity required to produce a hyperechoic spot in liver in < 1 s was 970 W/cm(2), in situ. At this HIFU dose, no immediate cellular damage was observed, providing a potential for pretreatment targeting. The real-time visualization method was used in hemostasis of actively bleeding internal pelvic vessels, allowing targeting and monitoring of successful treatment. Real-time US imaging may provide a useful tool for image-guided HIFU therapy.  相似文献   

15.
High-intensity focused ultrasound (HIFU) is a non-invasive tool that can be used for targeted thermal ablation treatments. Currently, HIFU is clinically approved for treatment of uterine fibroids, various cancers, and certain brain applications. However, for brain applications such as essential tremors, HIFU can only be used to treat limited areas confined to the center of the brain because of geometrical limitations (shape of the transducer and skull). A major obstacle to advancing this technology is the inability to treat non-central brain locations without causing damage to the skin and/or skull. Previous research has indicated that cavitation-induced bubbles or microbubble contrast agents can be used to enhance HIFU treatments by increasing ablation regions and shortening acoustic exposures at lower acoustic pressures. However, little research has been done to explore the interplay between microbubble concentration and pressure amplitude on HIFU treatments. We developed an in vitro experimental setup to study lesion formation at three different acoustic pressures and three microbubble concentrations. Real-time ultrasound imaging was integrated to monitor initial microbubble concentration and subsequent behavior during the HIFU treatments. Depending on the pressure used for the HIFU treatment, there was an optimal concentration of microbubbles that led to enhanced heating in the focal area. If the concentration of microbubbles was too high, the treatment was detrimentally affected because of non-linear attenuation by the pre-focal microbubbles. Additionally, the real-time ultrasound imaging provided a reliable method to monitor microbubble activity during the HIFU treatments, which is important for translation to in vivo HIFU applications with microbubbles.  相似文献   

16.
Assessing arterial-wall motion and deformations may reveal pathologic alterations in biomechanical properties of the parietal tissues and, thus, contribute to the detection of vascular disease onset. Ultrasound image sequences allow the observation of this motion and many methods have been developed to estimate temporal changes in artery diameter and wall thickness and to track 2-D displacements of selected points. Some methods enable the assessment of shearing or stretching within the wall, but none of them can estimate all these deformations simultaneously. The method herein proposed was devised to simultaneously estimate translation, compression, stretching and shearing of the arterial wall in ultrasound B-mode image sequences representing the carotid artery longitudinal section. Salient blob-like patterns, called key points, are automatically detected in each frame and matched between successive frames. A robust estimator based on an affine transformation model is then used to assess frame-to-frame motion explaining at best the key-point matches and to reject outliers. Realistic simulated image sequences were used to evaluate the accuracy and robustness of the method against ground truth. The method was also visually assessed on clinical image sequences, for which true deformations are unknown.  相似文献   

17.
Spectral parameter imaging in both the fundamental and harmonic of backscattered radio-frequency (RF) data were used for immediate visualization of high-intensity focused ultrasound (HIFU) lesion sites. A focused 5-MHz HIFU transducer with a coaxial 9-MHz focused single-element diagnostic transducer was used to create and scan lesions in chicken breast and freshly excised rabbit liver. B-mode images derived from the backscattered RF signal envelope were compared with midband fit (MBF) spectral parameter images in the fundamental (9-MHz) and harmonic (18-MHz) bands of the diagnostic probe. Images of HIFU-induced lesions derived from the MBF to the calibrated spectrum showed improved contrast (approximately 3 dB) of tumor margins versus surround compared with images produced from the conventional signal envelope. MBF parameter images produced from the harmonic band showed higher contrast in attenuated structures (core, shadow) compared with either the conventional envelope (3.3 dB core; 11.6 dB shadow) or MBF images of the fundamental band (4.4 dB core; 7.4 dB shadow). The gradient between the lesion and surround was 3.4 dB/mm, 6.9 dB/mm and 17.2 dB/mm for B-mode, MBF-fundamental mode and MBF-harmonic mode, respectively. Images of threshold and "popcorn" lesions produced in freshly excised rabbit liver were most easily visualized and boundaries best-defined using MBF-harmonic mode.  相似文献   

18.
One of the main problems encountered when using conventional B-mode ultrasound (US) for targeting and monitoring purposes during ablation therapies employing high-intensity focused US (HIFU) is the appearance of strong interference in the obtained diagnostic US images. In this study, instead of avoiding the interference noise, we demonstrate how we used it to locate the focus of the HIFU transducer in both in vitro tissue-mimicking phantoms and an ex vivo tissue block. We found that when the B-mode image plane coincided with the HIFU focal plane, the interference noise was maximally converged and enhanced compared with the off-focus situations. Stronger interference noise was recorded when the angle (alpha) between the US image plane and the HIFU axis was less than or equal to 90 degrees. By intentionally creating a target (group of bubbles) at the 3.5-MHz HIFU focus (7.1 mm in length and 0.7 mm in diameter), the position of the maximal noise convergence coincided well with the target. The difference between the predicted focus and the actual one (bubbles) on x and z axes (axes perpendicular to the HIFU central axis, Fig. 1) were both about 0.9 mm. For y axis (HIFU central axis), the precision was within 1.0 mm. For tissue block ablation, the interference noise concentrated at the position of maximal heating of the HIFU-induced lesions. The proposed method can also be used to predict the position of the HIFU focus by using a low intensity output scheme before permanent changes in the target tissue were made.  相似文献   

19.
The aim of this study was to establish a proof of concept for an ultrasonic image-guided high-intensity focused ultrasound (HIFU) breast treatment system. An ultrasonic computed tomography (UCT) scanner served as a platform for image acquisition and thermal mapping. With this system, images depicting the speed of sound (SOS) can be obtained. After HIFU activation, the resulting changes in SOS can be mapped. For the temperature range T ≤ 45°C, the temperature elevation is obtained directly from the change in SOS. For thermal ablation, monitoring is obtained by studying the temporal derivative of the SOS while continuous HIFU heating is applied. The method was implemented to in vitro tissue specimens and the ability for thermal monitoring was demonstrated. The results indicate that the suggested concept is indeed feasible.  相似文献   

20.
Patients with diabetic neuropathy often have neuropathic pain. The purpose of our work was to investigate the effects of high-intensity focused ultrasound (HIFU) on the conduction block of normal and neuropathic nerves for soothing pain. Adult male Sprague–Dawley rats were used, and diabetes was induced by streptozotocin injection. Diabetic neuropathy was evaluated with animal behavior tests. Sciatic nerves of both control and neuropathic rats were dissected from the starting point of the sciatic nerve to the point where the sural nerve ends near the ankle. The nerves were stored in Ringer's solution. The in vitro nerve was placed on a self-developed experimental platform for HIFU exposure. Stimulation and recording of the compound action potentials (CAPs) and sensory action potentials (SAPs) were performed. Control and neuropathic nerves exposed or not exposed to HIFU were submitted to histologic analysis. For the control and neuropathic nerves, suppression of CAPs and SAPs started 2 min post-HIFU treatment. Maximum suppression of SAPs was 34.4 ± 3.2% for the control rats and 11.6 ± 2.0% and 9.8 ± 3.0% for rats 4 wk post-injection and 8 wk post-injection, respectively. Time to full recovery was 25, 70 and 80 min, respectively. Histologic analysis revealed that the nerves in which CAPs and SAPs did not fully recover were damaged thermally or mechanically by HIFU. It is feasible to reversibly block nerves with appropriate HIFU treatment. Diabetic nerves were less suppressed by HIFU and were more vulnerable to permanent damage.  相似文献   

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

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