首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 78 毫秒
1.
MRI-guided microwave ablation (MWA) is a minimally invasive treatment for localized cancer. MR thermometry has been shown to be able to provide vital information for monitoring the procedure in real-time. However, MRI during active MWA can suffer from image quality degradation due to intermittent electromagnetic interference (EMI). A novel approach to correct for EMI-contaminated images is presented here to improve MR thermometry during clinical hepatic MWA. The method was applied to MR-thermometry images acquired during four MR-guided hepatic MWA treatments using a commercially available MRI-configured microwave generator system. During the treatments MR thermometry data acquisition was synchronized to respiratory cycle to minimize the impact of motion. EMI was detected and corrected using uncontaminated k-space data from nearby frames in k-space. Substantially improved temperature and thermal damage maps have been obtained and the treatment zone can be better visualized. Our ex vivo tissue sample study shows the correction introduced minimal errors to the temperature maps and thermal damage maps.  相似文献   

2.
In some dynamic magnetic resonance imaging (MRI) applications, the sample is still, and only the signal intensity changes with time. For such cases, the keyhole imaging principle can be used. In standard keyhole imaging, a low-frequency image signal is acquired, using a limited number of phase-encoding steps, which correspond to the rectangular sampling region in the k-space center. However, such a region practically never coincides with the position of the k-space points, which carry the most relevant low-frequency image information. In this paper we propose an improved keyhole method, which allows dynamic acquisition of a low-frequency image signal from selected most relevant k-space points via fast imaging mechanisms. Dynamic data acquisition is executed in the presence of time-varying magnetic-field (MF) gradients after single sample excitation. Special care has been taken in the design of the gradient sequence to minimize gradient load. This improved keyhole imaging method has been considered theoretically and verified experimentally on a model system.  相似文献   

3.

Purpose

In cardiac MRI (cMRI) injection of contrast medium may be performed prior to the acquisition of cine steady-state free precession (SSFP) imaging to speed up the protocol and avoid delay before late Gadolinium enhancement (LGE) imaging. Aim of this study was to evaluate whether a condensed clinical protocol with contrast cine SSFP imaging is able to detect early microvascular obstruction (MO) and determine the infarct size compared to the findings of LGE inversion recovery sequences.

Materials and Methods

The study complies with the Declaration of Helsinki and was performed following approval by the ethic committee of the University of Erlangen-Nuremberg. Written informed consent was obtained from every patient. 68 consecutive patients (14 females/54 males) with acute ST-elevation myocardial infarction (STEMI) treated by percutaneous coronary revascularization were included in this study. CMRI was performed 6.6±2 days after symptom onset and MO and infarct size in early contrast SSFP cine imaging were compared to LGE imaging.

Results

MO was detected in 47/68 (69%) patients on cine SSFP and in 41/68 (60%) patients on LGE imaging. In 6 patients MO was found on cine SSFP imaging but was not detectable on LGE imaging. Infarct size on cine SSFP showed a strong agreement to LGE imaging (intraclass correlation coefficient [ICC] of 0.96 for enddiastolic, p<0.001 and 0.96 for endsystolic, p<0.001 respectively). Significant interobserver agreement was found measuring enddiastolic and endsystolic infarct size on cine SSFP imaging (p<0.01).

Conclusions

In patients after STEMI infarct size and presence of MO can be detected with contrast cine SSFP imaging. This could be an option in patients who are limited in their ability to comply with the demands of a cMRI protocol.  相似文献   

4.
The feasibility of using a commercial magnetic resonance imaging (MRI) scanner to do either imaging or hyperthermic treatment was demonstrated. Radiofrequency (RF) induced focal heating of phantoms and animal tissues was performed using a MRI scanner as the RF power source and a grounded interstitial probe as a device to produce hyperthermia via eddy current convergence. In the therapeutic mode, a pulse width of 900 microseconds and interval of 50 ms were used to give 2% duty cycle (closest simulation to continuous wave (CW) mode without bypassing imaging filters). Temperature in the vicinity of the grounded probe was measured with a field nonperturbing fluoroptic probe. Temperatures increased 4.5 degrees C in 5 minutes in a dielectrically uniform phantom, 3.1 degrees C in 6.7 minutes in rats' leg muscles, and 5.0 degrees C in 6.0 minutes in rats' peritoneum. The MRI of the phantom with the grounded probe and the fluoroptic probe was obtained using spin echo sequences. The potential advantage of this approach is visualization of deep-seated tumors and hyperthermic treatment with minimal modification of the MRI scanner.  相似文献   

5.
Sparse MRI has been introduced to reduce the acquisition time and raw data size by undersampling the k-space data. However, the image quality, particularly the contrast to noise ratio (CNR), decreases with the undersampling rate. In this work, we proposed an interpolated Compressed Sensing (iCS) method to further enhance the imaging speed or reduce data size without significant sacrifice of image quality and CNR for multi-slice two-dimensional sparse MR imaging in humans. This method utilizes the k-space data of the neighboring slice in the multi-slice acquisition. The missing k-space data of a highly undersampled slice are estimated by using the raw data of its neighboring slice multiplied by a weighting function generated from low resolution full k-space reference images. In-vivo MR imaging in human feet has been used to investigate the feasibility and the performance of the proposed iCS method. The results show that by using the proposed iCS reconstruction method, the average image error can be reduced and the average CNR can be improved, compared with the conventional sparse MRI reconstruction at the same undersampling rate.  相似文献   

6.
Cardiovascular transgenic mouse models with an early phenotype or even premature death require noninvasive imaging methods that allow for accurate visualization of cardiac morphology and function. Thus the purpose of our study was to assess the feasibility of magnetic resonance imaging (MRI) to characterize cardiac function and mass in newborn, juvenile, and adult mice. Forty-five C57bl/6 mice from seven age groups (3 days to 4 mo after birth) were studied by MRI under isoflurane anesthesia. Electrocardiogram-gated cine MRI was performed with an in-plane resolution of (78-117 microm)(2). Temporal resolution per cine frame was 8.6 ms. MRI revealed cardiac anatomy in mice from all age groups with high temporal and spatial resolution. There was close correlation between MRI- and autopsy-determined left ventricular (LV) mass (r = 0.95, SE of estimate = 9.5 mg). The increase of LV mass (range 9.6-101.3 mg), cardiac output (range 1.1-14.3 ml/min), and stroke volume (range 3. 2-40.2 microl) with age could be quantified by MRI measurements. Ejection fraction and cardiac index did not change with aging. However, LV mass index decreased with increasing age (P < 0.01). High-resolution MRI allows for accurate in vivo assessment of cardiac function in neonatal, juvenile, and adult mice. This method should be useful when applied in transgenic mouse models.  相似文献   

7.
The three-dimension (3D) magnetization-prepared rapid gradient-echo (MP-RAGE) sequence is one of the most popular sequences for structural brain imaging in clinical and research settings. The sequence captures high tissue contrast and provides high spatial resolution with whole brain coverage in a short scan time. In this paper, we first computed the optimal k-space sampling by optimizing the contrast of simulated images acquired with the MP-RAGE sequence at 3.0 Tesla using computer simulations. Because the software of our scanner has only limited settings for k-space sampling, we then determined the optimal k-space sampling for settings that can be realized on our scanner. Subsequently we optimized several major imaging parameters to maximize normal brain tissue contrasts under the optimal k-space sampling. The optimal parameters are flip angle of 12°, effective inversion time within 900 to 1100 ms, and delay time of 0 ms. In vivo experiments showed that the quality of images acquired with our optimal protocol was significantly higher than that of images obtained using recommended protocols in prior publications. The optimization of k-spacing sampling and imaging parameters significantly improved the quality and detection sensitivity of brain images acquired with MP-RAGE.  相似文献   

8.
The purpose of this project was to assess the feasibility of imaging the velopharynx of adult volunteers during repetitive speech, using gated magnetic resonance imaging (MRI). Although a number of investigators have used conventional MRI in the study of the human vocal tract, the mismatch between the lengthy time necessary to acquire sufficiently detailed images and the rapidity of movement of the vocal tract during speech has forced investigators to acquire images either while the subject is at rest or during sustained utterances. The technique used here acquired a portion of each image during repetitive utterances, building the full image over multiple utterance cycles. The velopharyngeal portal was imaged on a 1.5-Tesla GE Signa LX 8.2 platform with gated fast spoiled gradient echo protocol. An external 1-Hertz trigger was fed to the cardiac gate. Subjects synchronized utterance of consonant-vowel syllables to a flashing light synchronized with the external trigger. Each acquisition of 30 phases per second at a single-slice location took 22 to 29 seconds. Four consonant-vowel syllables (/pa/, /ma/, /sa/, and /ka/) were evaluated. Subjects vocalized throughout the acquisition, beginning 5 to 6 seconds beforehand to establish a regular rhythm. Imaging of the velopharyngeal portal was performed for sagittal, velopharyngeal axial (aligned perpendicular to the "knee" of the velum), axial, and coronal planes. Volumes were obtained by sequential acquisition of six to 10 slices (each with 30 phases) in the axial or sagittal planes during repetition of the /pa/ syllable. Spatiotemporal volumes of the single-slice data were sectioned to provide time-motion images (analogous to M-mode echocardiograms). Three-dimensional dynamic volume renderings of palate motion were displayed interactively (Vortex; CieMed, Singapore). A method suitable for the collection and visualization of four-dimensional information regarding monosyllabic speech using gated MRI was developed. These techniques were applied to a population of adult volunteer subjects with no history of speech problems and two patients with a history of cleft lip and palate. The techniques allowed good real-time visualization of velopharyngeal anatomy during its entire range of motion and was also able to image pathology-specific anatomic differences in the subjects with cleft lip and cleft palate. These methods may be applicable to a wide spectrum of problems in speech physiology research and for clinical decision-making regarding surgery for speech and outcomes analysis.  相似文献   

9.

Purpose

Respiratory motion correction remains a challenge in coronary magnetic resonance imaging (MRI) and current techniques, such as navigator gating, suffer from sub-optimal scan efficiency and ease-of-use. To overcome these limitations, an image-based self-navigation technique is proposed that uses “sub-images” and compressed sensing (CS) to obtain translational motion correction in 2D. The method was preliminarily implemented as a 2D technique and tested for feasibility for targeted coronary imaging.

Methods

During a 2D segmented radial k-space data acquisition, heavily undersampled sub-images were reconstructed from the readouts collected during each cardiac cycle. These sub-images may then be used for respiratory self-navigation. Alternatively, a CS reconstruction may be used to create these sub-images, so as to partially compensate for the heavy undersampling. Both approaches were quantitatively assessed using simulations and in vivo studies, and the resulting self-navigation strategies were then compared to conventional navigator gating.

Results

Sub-images reconstructed using CS showed a lower artifact level than sub-images reconstructed without CS. As a result, the final image quality was significantly better when using CS-assisted self-navigation as opposed to the non-CS approach. Moreover, while both self-navigation techniques led to a 69% scan time reduction (as compared to navigator gating), there was no significant difference in image quality between the CS-assisted self-navigation technique and conventional navigator gating, despite the significant decrease in scan time.

Conclusions

CS-assisted self-navigation using 2D translational motion correction demonstrated feasibility of producing coronary MRA data with image quality comparable to that obtained with conventional navigator gating, and does so without the use of additional acquisitions or motion modeling, while still allowing for 100% scan efficiency and an improved ease-of-use. In conclusion, compressed sensing may become a critical adjunct for 2D translational motion correction in free-breathing cardiac imaging with high spatial resolution. An expansion to modern 3D approaches is now warranted.  相似文献   

10.
Chan DD  Neu CP 《PloS one》2012,7(3):e33463
Characterization of spatiotemporal deformation dynamics and material properties requires non-destructive methods to visualize mechanics of materials and biological tissues. Displacement-encoded magnetic resonance imaging (MRI) has emerged as a noninvasive and non-destructive technique used to quantify deformation and strains. However, the techniques are not yet applicable to a broad range of materials and load-bearing tissues. In this paper, we visualize transient and internal material deformation through the novel synchrony of external mechanical loading with rapid displacement-encoded MRI. We achieved deformation measurements in silicone gel materials with a spatial resolution of 100 μm and a temporal resolution (of 2.25 ms), set by the repetition time (TR) of the rapid MRI acquisition. Displacement and strain precisions after smoothing were 11 μm and 0.1%, respectively, approaching cellular length scales. Short (1/2 TR) echo times enabled visualization of in situ deformation in a human tibiofemoral joint, inclusive of multiple variable T(2) biomaterials. Moreover, the MRI acquisitions achieved a fivefold improvement in imaging time over previous technology, setting the stage for mechanical imaging in vivo. Our results provide a general approach for noninvasive and non-destructive measurement, at high spatial and temporal resolution, of the dynamic mechanical response of a broad range of load-bearing materials and biological tissues.  相似文献   

11.
In this paper, we propose a method for non-invasively measuring three-dimensional in vivo kinematics of the ankle joint from a dynamic MRI acquisition of a single range-of-motion cycle. The proposed approach relies on an intensity-based registration method to estimate motion from multi-plane dynamic MRI data. Our approach recovers not only the movement of the skeleton, but also the possibly non-rigid temporal deformation of the joint. First, the rigid motion of each ankle bone is estimated. Second, a four-dimensional (3D+time) high-resolution dynamic MRI sequence is estimated through the use of the log-euclidean framework for the computation of temporal dense deformation fields. This approach has been then applied and evaluated on in vivo dynamic MRI data acquired for a pilot study on six healthy pediatric cohort in order to establish in vivo normative joint biomechanics. Results demonstrate the robustness of the proposed pipeline and very promising high resolution visualization of the ankle joint.  相似文献   

12.
A validation study and early results for non-invasive, in vivo measurement of coronary artery blood flow using phase contrast magnetic resonance imaging (PC-MRI) at 3.0T is presented. Accuracy of coronary artery blood flow measurements by phase contrast MRI is limited by heart and respiratory motion as well as the small size of the coronary arteries. In this study, a navigator echo gated, cine phase velocity mapping technique is described to obtain time-resolved velocity and flow waveforms of small diameter vessels at 3.0T. Phantom experiments using steady, laminar flow are presented to validate the technique and show flow rates measured by 3.0T phase contrast MRI to be accurate within 15% of true flow rates. Subsequently, in vivo scans on healthy volunteers yield velocity measurements for blood flow in the right, left anterior descending, and left circumflex arteries. Measurements of average, cross-sectional velocity were obtainable in 224/243 (92%) of the cardiac phases. Time-averaged, cross-sectional velocity of the blood flow was 6.8+/-4.3cm/s in the LAD, 8.0+/-3.8cm/s in the LCX, and 6.0+/-1.6cm/s in the RCA.  相似文献   

13.
To understand the structure-function relationship in the postinfarcted myocardium in rabbits, we induced cardiac ischemia by ligating the left circumflex coronary artery. Sham controls underwent thoracotomy only. At 7 and 30 d after ligation, cardiac MRI was conducted by using pulse-oxymetry-gated cine acquisition to provide complete phases of the heartbeat. The rabbits were anesthetized under 1.5% isoflurane ventilation, and ultrafast techniques made breath-hold 3D coverage in different cardiac axes feasible. Viability imaging was performed after intravenous injection of 0.15 mmol/kg gadolinium to assess the extent of infarction. Data (n ≥ 6) are presented as mean ± SEM and analyzed by ANOVA and ANCOVA. In postligation rabbits, end-systolic (mean ± SEM, 2.3 ± 0.3 mL) and end-diastolic (4.2 ± 0.4 mL) volumes were increased compared with preligation values (end-systolic, 1.1 ± 0.1 mL; end-diastolic, 2.98 ± 0.2 mL). Ejection fraction was influenced adversely by the presence of scar tissue at both 7 and 30 d after ligation and apparently nonlinear with the heart rate. Cardiac force was increased in the basal region in both end-systole and end-diastole in postligation hearts but progressively decreased toward the apex. Late gadolinium enhancement delineated 15.2 ± 5.8% myocardial infarction at 7 d after ligation and 14.5 ± 5.8% at 30 d, with limited wall motion and wall thinness. Compensatory wall thickening was present in the basal region when compared with that in preligation hearts. MRI offers detailed spatial resolution and tissue characterization after myocardial infarction.  相似文献   

14.
The objective of this study was to apply cine magnetic resonance imaging (MRI) using displacement encoding with stimulated echoes (DENSE) to measure the dynamic two-dimensional (2D) displacement and Lagrangian strain fields in the biceps brachii muscle. Six healthy volunteers underwent cine DENSE MRI during repeated elbow flexion against the load of gravity. Displacement encoded dynamic images of the upper arm were acquired with spatial and temporal resolutions of 1.9 x 1.9 mm(2) and 30 ms, respectively. Pixel-wise Lagrangian displacement and strain fields were calculated from the measured images. We extracted the first and second principal strains (E1 and E2) along the centerline and anterior regions of the muscle. E1 and E2 were relatively uniform along the anterior region. However, E1 and E2 were both non-uniform along the centerline region-normalized values for E1 and E2 varied over the ranges of 0.27-1.35, and 0.45-2.36, respectively. The directions of the first and second principal strains varied throughout the muscle and showed that the direction of principal shortening is not necessarily aligned with fascicle direction. This study demonstrates the utility of cine DENSE MRI for analyzing skeletal muscle mechanics and provides data describing the in vivo mechanics of muscle tissue to a level of detail that has not been previously possible.  相似文献   

15.
Blood oxygen level-dependent (BOLD) MRI is a non-invasive diagnostic method for assessing tissue oxygenation level, by changes in the transverse relaxation time T2*. 3D BOLD imaging of lung tumours is challenging, because respiratory motion can lead to significant image quality degradation. The purpose of this work was to explore the feasibility of a three dimensional (3D) Cartesian multi gradient echo (MGRE) sequence for T2* measurements of non-small cell lung tumours during free-breathing. A non-uniform quasi-random reordering of the pahse encoding lines that allocates more sampling points near the k-space origin resulting in efficient undersampling pattern for parallel imaging was combined with multi echo acquisition and self-gating. In a series of three patients 3D T2* maps of lung carcinomas were generated with isotropic spatial resolution and full tumour coverage at air inhalation and after hyperoxic gas challenge in arbitrary respiratory phases using the proposed self-gated MGRE acquisition. The changes in T2* on the inhalation of hyperoxic gas relative to air were quantified. Significant changes in T2* were observed following oxygen inhalation in the tumour (p < 0.02). Thus, the self-gated MGRE sequence can be used for assessment of BOLD signal with isotropic resolution and arbitrary respiratory phases in non-small cell lung cancer.  相似文献   

16.
PurposeOne of the challenges of cardiac MR imaging is the compensation of respiratory motion, which causes the heart and the surrounding tissues to move. Commonly-used methods to overcome this effect, breath-holding and MR navigation, present shortcomings in terms of available acquisition time or need to periodically interrupt the acquisition, respectively. In this work, an implementation of respiratory motion compensation that obtains information from abdominal ultrasound and continuously adapts the imaged slice position in real time is presented.MethodsA custom workflow was developed, comprising an MR-compatible ultrasound acquisition system, a feature-motion-tracking system with polynomial predictive capability, and a custom MR sequence that continuously adapts the position of the acquired slice according to the tracked position. The system was evaluated on a moving phantom by comparing sharpness and image blurring between static and moving conditions, and in vivo by tracking the motion of the blood vessels of the liver to estimate the cardiac motion. Cine images of the heart were acquired during free breathing.ResultsIn vitro, the predictive motion correction yielded significantly better results than non-predictive or non-corrected acquisitions (p ≪ 0.01). In vivo, the predictive correction resulted in an image quality very similar to the breath-hold acquisition, whereas the uncorrected images show noticeable blurring artifacts.ConclusionIn this work, the possibility of using ultrasound navigation with tracking for the real-time adaptation of MR imaging slices was demonstrated. The implemented technique enabled efficient imaging of the heart with resolutions that would not be feasible in a single breath-hold.  相似文献   

17.
This work’s aim was to minimize the acquisition time of a radial 3D ultra-short echo-time (UTE) sequence and to provide fully automated, gradient delay compensated, and therefore artifact free, reconstruction. The radial 3D UTE sequence (echo time 60 μs) was implemented as single echo acquisition with center-out readouts and improved time efficient spoiling on a clinical 3T scanner without hardware modifications. To assess the sequence parameter dependent gradient delays each acquisition contained a quick calibration scan and utilized the phase of the readouts to detect the actual k-space center. This calibration scan does not require any user interaction. To evaluate the robustness of this automatic delay estimation phantom experiments were performed and 19 in vivo imaging data of the head, tibial cortical bone, feet and lung were acquired from 6 volunteers. As clinical application of this fast 3D UTE acquisition single breath-hold lung imaging is demonstrated. The proposed sequence allowed very short repetition times (TR~1ms), thus reducing total acquisition time. The proposed, fully automated k-phase based gradient delay calibration resulted in accurate delay estimations (difference to manually determined optimal delay −0.13 ± 0.45 μs) and allowed unsupervised reconstruction of high quality images for both phantom and in vivo data. The employed fast spoiling scheme efficiently suppressed artifacts caused by incorrectly refocused echoes. The sequence proved to be quite insensitive to motion, flow and susceptibility artifacts and provides oversampling protection against aliasing foldovers in all directions. Due to the short TR, acquisition times are attractive for a wide range of clinical applications. For short T2* mapping this sequence provides free choice of the second TE, usually within less scan time as a comparable dual echo UTE sequence.  相似文献   

18.
Coronary bifurcations represent specific regions of the arterial tree that are susceptible to atherosclerotic lesions. While the effects of vessel compliance, curvature, pulsatile blood flow, and cardiac motion on coronary endothelial shear stress have been widely explored, the effects of myocardial contraction on arterial wall stress/strain (WS/S) and vessel stiffness distributions remain unclear. Local increase of vessel stiffness resulting from wall-strain stiffening phenomenon (a local process due to the nonlinear mechanical properties of the arterial wall) may be critical in the development of atherosclerotic lesions. Therefore, the aim of this study was to quantify WS/S and stiffness in coronary bifurcations and to investigate correlations with plaque sites. Anatomic coronary geometry and cardiac motion were generated based on both computed tomography and MRI examinations of eight patients with minimal coronary disease. Computational structural analyses using the finite element method were subsequently performed, and spatial luminal arterial wall stretch (LW(Stretch)) and stiffness (LW(Stiff)) distributions in the left main coronary bifurcations were calculated. Our results show that all plaque sites were concomitantly subject to high LW(Stretch) and high LW(Stiff), with mean amplitudes of 34.7 ± 1.6% and 442.4 ± 113.0 kPa, respectively. The mean LW(Stiff) amplitude was found slightly greater at the plaque sites on the left main coronary artery (mean value: 482.2 ± 88.1 kPa) compared with those computed on the left anterior descending and left circumflex coronary arteries (416.3 ± 61.5 and 428.7 ± 181.8 kPa, respectively). These findings suggest that local wall stiffness plays a role in the initiation of atherosclerotic lesions.  相似文献   

19.
《Médecine Nucléaire》2007,31(11):597-603
This article presents the current and future possibilities of the cardiac magnetic resonance imaging (MRI) in the ischemic heart disease. It is not an emergency technique, but is becoming an element of decision for an acute myocardial infarction. It makes it possible to visualize the extent of the myocardial lesions in post acute infarction, to appreciate functional recovery and to uncover a nonsymptomatic coronary stenosis or of atypical clinical presentation. The visualization of the coronary arteries in MRI is realizable, but remains experimentation field. The cardiac MRI place is to be defined by each medical team according to their diagram of ischemic heart disease follow-up.  相似文献   

20.
Imaging of cartilage has traditionally been achieved indirectly with conventional radiography. Loss of joint space width, or 'joint space narrowing', is considered a surrogate marker for cartilage thinning. However, radiography is severely limited by its inability to visualize cartilage, the difficulty of ascertaining the optimum and reproducible positioning of the joint in serial assessments, and the difficulty of grading joint space narrowing visually. With the availability of advanced magnetic resonance imaging (MRI) scanners, new pulse sequences, and imaging techniques, direct visualization of cartilage has become possible. MRI enables visualization not only of cartilage but also of other important features of osteoarthritis simultaneously. 'Pre-radiographic' cartilage changes depicted by MRI can be measured reliably by a semiquantitative or quantitative approach. MRI enables accurate measurement of longitudinal changes in quantitative cartilage morphology in knee osteoarthritis. Moreover, compositional MRI allows imaging of 'pre-morphologic' changes (that is, visualization of subtle intrasubstance matrix changes before any obvious morphologic alterations occur). Detection of joint space narrowing on radiography seems outdated now that it is possible to directly visualize morphologic and pre-morphologic changes of cartilage by using conventional as well as complex MRI techniques.  相似文献   

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

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