CN101322648A - NMR imaging equipment stability and method for measuring image-forming index - Google Patents
NMR imaging equipment stability and method for measuring image-forming index Download PDFInfo
- Publication number
- CN101322648A CN101322648A CNA2008100456644A CN200810045664A CN101322648A CN 101322648 A CN101322648 A CN 101322648A CN A2008100456644 A CNA2008100456644 A CN A2008100456644A CN 200810045664 A CN200810045664 A CN 200810045664A CN 101322648 A CN101322648 A CN 101322648A
- Authority
- CN
- China
- Prior art keywords
- image
- data
- magnetic resonance
- area
- nuclear magnetic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Abstract
The invention discloses a method for measuring the stability and imaging indexes of a nuclear magnetic resonance imaging device and is characterized in that: the images of a water phantom are obtained by adopting the nuclear magnetic resonance imaging device to scan the water phantom, the images are transmitted to a specific nuclear magnetic resonance imaging quality assurance and analysis work station by a DICOM protocol, and then the quality assurance and analysis work station is self started and enters a corresponding measurement process (referring to the figure 1 in the abstract of the process) according to DICOM control elements so as to process and calculate the images of the water phantom; the measuring content includes common nuclear magnetic resonance imaging quality assurance indexes, the shifting indexes which have important significance to Functional MRI, and other indexes; meanwhile, an automatic early warning mechanism is established for auto prompting when an abnormal condition occurs and gives the possible reason, and all measured results are sent into a database for long-term analysis of data and evaluation of running state of the device, so as to offer guidance to relevant technicians to master device state and carry out equipment testing.
Description
Technical field
The present invention relates to the measuring method of a kind of NMR (Nuclear Magnetic Resonance)-imaging stabilization of equipment performance and imaging index, belong to the application of nuclear magnetic resonance equipment.
Background technology
Since first nuclear magnetic resonance scanner in 1972 came out, (Magnetic ResonanceImage MRI) had obtained development at full speed as important clinical inspection means to NMR (Nuclear Magnetic Resonance)-imaging, is bringing into play important effect in clinical diagnosis work.Simultaneously, the imaging mode of nuclear magnetic resonance, NMR also all advances with day, and (functional Magnetic Resonance Image is fMRI) at some interior new imaging techniques functional mri to have occurred comprising.But as a kind of very complicated system, NMR (Nuclear Magnetic Resonance)-imaging is subjected to the influence of several factors and link, and there is the lot of unstable factor in equipment itself, is reflected as the image quality problem at last, causes mistaken diagnosis and wrong experimental result.Can the NMR (Nuclear Magnetic Resonance)-imaging equipment running status have determined the gained image accurately reflect patient's physiologic character, and the staff who therefore is engaged in NMR (Nuclear Magnetic Resonance)-imaging need understand the operation conditions of equipment usually.The MRI quality assurance then is to guarantee MRI clinical diagnosis quality and the prerequisite of carrying out smoothly based on the Neuroscience Research of MRI.MRI quality assurance (Quality Assurence, be called for short QA) be meant MRI equipment is carried out the problem of quality testing, performance evaluation discovering device existence in time and solves, make the performance of MRI equipment reach best, guarantee the correctness of diagnostic imaging and experiment.
Because the complexity of nuclear magnetic resonance, therefore, the quality control of magnetic resonance also relates to all many-sides.Examination after MRI equipment receives examination, Breakdown Maintenance, daily state detection and quality all relate to QA work when evaluating.Generally test being divided into three kinds uses: check and accept detection, Performance Detection, daily test.Different phase QA test job is different.Checking and accepting to detect with Performance Detection often organizes the professional and technical personnel to carry out in the specific time by the MRI equipment supplier.Institute's detection method is also provided by the equipment separate procedure.QA involved in the present invention partly is meant daily test, the detection that daily test is carried out before and after the start in every day by the applying unit of equipment, the project that detects is several crucial imaging parameters, daily test lays particular emphasis on index important and not consuming time, and be the trend analysis accumulation data, generally be arranged in and carry out before the formal inspection patient of start back morning, the water mould of simulated human tissue magnetics parameter is passed through in test, scanning back detection signal-to-noise ratio, bed thickness, spatial resolution and linearity.Resonant frequency and RF power also need record in addition; In logging test results, test parameter also needs record, has adopted special pulse train when particularly testing.
That abroad, be engaged in the magnetic resonance quality control is special magnetic resonance engineer.Because domestic conditions is limited, the responsibility of monitoring magnetic resonance operation conditions is served as by radiologic technologist usually.Show that according to statistics the operation conditions of the present magnetic resonance tool of China is not up to standard fully.Present existing problem is mainly as follows: (1) test job loaded down with trivial details: usually, be used for clinical NMR (Nuclear Magnetic Resonance)-imaging equipment before scanning patient every day, all to carry out the QA test of equipment, promptly the water mould is scanned with the MRI scanner unit, then scanogram is handled and measured according to correlation technique, obtain the QA index and analyze and record.This work is mainly finished by manual, and process is loaded down with trivial details; (2) test has very strong subjectivity, because manual ROI zone or the axis etc. selected lack unified standard, data are subjected to operator's influence; (3) record of parameter mainly is undertaken by craft with preserving, and is unfavorable for long preservation and to the statistical analysis of historical data.The state of equipment and its change are that the QA data statistics by the long period draws, and lack statistics and can't give full play to QA work importance, and this also is the reason that QA work is not fully carried out.Though have many relevant professional softwares to finish the QA process at present, these softwares are often provided by the MRI equipment supplier, do not possess general scan-data disposal ability and simple external data interface.Some specialized function is powerful on the other hand, but the implementation manual intervention is too much, complicated operation, even need carry out some Programming Design, be unfavorable for radiologic technologist work.This just presses for a cover at the magnetic resonance tool quality control, simple method.
Summary of the invention
The objective of the invention is to propose the measuring method of a kind of NMR (Nuclear Magnetic Resonance)-imaging stabilization of equipment performance and imaging index, be characterized in objective, automatically detect and the stability parameter of analytical equipment at the deficiencies in the prior art, avoid manual detection with subjectivity.
Purpose of the present invention is realized by following technical measures:
The measuring method of NMR (Nuclear Magnetic Resonance)-imaging stabilization of equipment performance and imaging index may further comprise the steps:
1 utilizes the DICOM agreement to finish (QA) process control of NMR (Nuclear Magnetic Resonance)-imaging quality assurance
(1) server end (SCP) keeps accepting state at any time in the mode of listening port, when receiving the association request of client (SCU is a MRI equipment work station) here, accept related request and judge whether COS is storage class, if then accept, otherwise refusal and continuation are monitored;
(2) receive message data,, judge whether to discharging related request to the message data that receives, if then release is related, and execution in step (3), otherwise execution in step (4);
(3) when judging message data when discharging related request, discharge related and judge whether view data (DQA) into the quality assurance, if execution in step 3 then continues to monitor, otherwise finish relatedly, continue to keep listening state;
(4) when the message data that receives is not the related request of release, execution in step 2, and repeating step (2);
2 extraction data also write the data base
(1) resolve the packet that receives, judge whether to accept data according to DICOM control label, selecting the control label here is patient name (PatientName), if patient name is " mriDQA " or " fmriDQA ", then accept data, and continue next step, otherwise finish;
(2) inspection is unified identifier (Study UID) and is judged whether to write the data base according to DICOM control label, if then continue next step, otherwise finishes;
(3) in the data base, set up new record, write date, scan protocols, mid frequency (AX) and amplification, and the data storage path is write the data base;
3 NMR (Nuclear Magnetic Resonance)-imaging quality assurance index analysis are measured
When server end receives the related request of release, discharge related also start-up routine and begin the detection and the calculating of NMR (Nuclear Magnetic Resonance)-imaging quality assurance (QA) index, the steps include:
(1) judges whether data according to label, if then calculate the QA index as follows into fMRI;
1. computed image signal to noise ratio (snr); With repeatedly being with threshold method to cut apart image segmentation is background area and prospect, and prospect is circular water mould image, circular water mould image with the center of circle be the center select 80% area of a circle zone for signaling zone area-of-interest (ROI); Circular water mould image is background area with the exterior domain, to the background area with repeatedly being with threshold method to cut apart once more, obtain pseudo-shadow zone of ghost and background noise district, choosing 2% image area in the left-right symmetric outside in background noise district is area-of-interest (ROI), extract signaling zone area-of-interest average gray value and noise range area-of-interest standard deviation, that calculates the two obtains signal noise ratio (snr) of image (SNR) than promptly;
2. the uniformity of measurement image; To cutting apart the signaling zone of back image, in selected area-of-interest (ROI) zone, obtain the signal value of bright spot of image and dim spot, by formula the uniformity of computed image;
3. measurement image torsion resistance; With repeatedly being with the threshold method split image, the prospect of image is asked its maximum gauge L with the searching algorithm of translation-angle
Max, minimum diameter L
MinAnd average diameter L
m, the torsion resistance D:D=(L of computed image
Max-L
Min)/L
m
4. the pseudo-shadow (GSR) of the ghost of image; With repeatedly being with the threshold method split image, image segmentation is once more carried out in the noise range of split image, prospect is the pseudo-shadow zone of ghost, background is the noise range, can try to achieve the pseudo-shadow of ghost of image:
GSR=(the pseudo-shadow meansigma methods of ghost-noise range meansigma methods)/signaling zone meansigma methods
5. the signal drift of sequence of computed images; A selection area to temporal sequence of images writes down its signal value one by one, represents its variation tendency with matched curve, and calculates drift value:
S
Max, S
MinMaximum, minima for selection area signal in the temporal sequence of images; Match value and actual value are subtracted each other, and its difference obtains its frequency spectrum with fast fourier transform, searches for maximum in frequency spectrum, and its range value and frequency are write the data base;
(2) if judge it is not functional mri (fMRI) data, then calculate MRI quality assurance (QA) index as follows:
1. computed image signal to noise ratio; In image sequence, select piece image, with repeatedly being with threshold method to cut apart image segmentation is background and prospect, in background and prospect, select the area-of-interest of noise range and signaling zone, extract signaling zone area-of-interest average gray value and noise range area-of-interest standard deviation, that calculates the two obtains signal noise ratio (snr) of image than promptly;
2. the uniformity of measurement image; To cutting apart the signaling zone of back image, in selected area-of-interest, obtain the signal value of bright spot of image and dim spot, according to the uniformity of formula computed image;
3. the linearity of measurement image; In the known water mould under the prerequisite of size,, calculate the linearity of image according to water mould length on respective frequencies and the phase-encoding direction in image;
(3) the QA indication information that calculates is write the data base;
(4) with result of calculation and normal data and historical average according to contrast, going beyond the scope then provides prompting.
Description of drawings:
Fig. 1 finishes QA process control flow figure for the DICOM agreement
Fig. 2 is QA data fetch and data base's the flow chart that writes
The flowchart that Fig. 3 measures for QA
The specific embodiment:
Below by embodiment the present invention is carried out concrete description; be necessary to be pointed out that at this present embodiment only is used for invention is further specified; but can not be interpreted as limiting the scope of the invention, the person skilled in the art in this field can do some nonessential improvement and adjustment according to the content of the invention described above.
Embodiment
The image sequence that client will scan the water mould utilizes the DICOM agreement to send to server end (server end is the QA data workstation, and client is the equipment work station), and server end determines whether to store and start the QA analysis process through judging.
Server end utilizes the DICOM agreement to finish the QA process control and realizes according to the following steps:
When (1) receiving the association request message of client, decode messages also judges whether to be storage class that if then make up response message A-AssociateAC PDU, client begins to transmit data; Otherwise refusal.
(2) server end reception message data, decoding and judging whether is to discharge association request message; And if judgment data discharges related beginning QA analysis process when being the DQA data.Then continue to connect data-message if not, execution in step 3 is until receiving the related request that discharges.(3) resolve the packet that receives, meet under the condition of acceptance prerequisite at DICOM control label PatientName and Study UID, in the data base, set up new record, write date, scan protocols, mid frequency and amplification AX, and the data storage path is write the data base.
When server end receives that client discharges related message, start the measurement that corresponding process begins the QA data according to the PatientName label.If the PatientName label value is " fmriDQA ", then be judged as the fMRI data, measure content and comprise: signal to noise ratio, uniformity degree of image, torsion resistance, the pseudo-shadow of ghost and signal drift.Introduce the measuring process of every index below:
(1) signal to noise ratio (snr); With repeatedly being with threshold method to cut apart the DICOM image, with image segmentation prospect and background, prospect is the part of middle section scanning water mould, and background is dry acoustical signal; Select the area-of-interest (ROI) of sense signal and dry sound, determine the central point of image according to the pixel of image, center with image is that circle is drawn in the center of circle, select suitable radius, make that the area of drawn circle is 70% of a foreground area, the zone that this circle is enclosed is the area-of-interest (ROI) of the signaling zone of selection.When choosing the area-of-interest of noise range,, should avoid pseudo-shadow zone owing to there is the pseudo-shadow of ghost.Image segmentation is once more carried out in the noise range, and the background after cutting apart is the noise range, and prospect is the pseudo-shadow zone of ghost, chooses noise range ROI and then can avoid pseudo-shadow zone in background.Extract signaling zone area-of-interest average gray value and noise range area-of-interest standard deviation, the ratio of the two is the signal to noise ratio of trying to achieve.
(2) torsion resistance; Torsion resistance has reflected the change of body form on image, for example may be ellipse with circular water mould scanning at image.Before the calculating of the torsion resistance of carrying out geometry deformation, cut apart image, threshold method repeatedly is with in the same employing of method.
What torsion resistance was paid close attention to is the prospect (water mould) of image segmentation, at first will find the geometric center of prospect.Here we are with the meansigma methods of the horizontal stroke of prospect, vertical coordinate horizontal stroke, the vertical coordinate as central point.Cross the axle of central point and the intersection point at prospect edge when asking angle to be zero, calculate and write down the distance of two intersection points.Translation-angle repeats the step, obtains the maximum gauge L through the axial length of central point
Max, minimum diameter L
MinAnd average diameter L
m, with formula D=(L
Max-L
Min)/L
mThe torsion resistance of computed image.
(3) uniformity degree of image (PIU); Uniformity degree of image refers to the uniformity coefficient of image intensity signal response, and it has described the reproduction of MRI system to same substance zone in the phantom.
The measurement image uniformity at first needs split image, and threshold method repeatedly is with in same here employing.In the water mould part of cutting apart the back image, be that circle is drawn in the center of circle with the centre coordinate, making round area is 70% of water mould image area, the ROI that this border circular areas is measured as the uniformity.Record the signal value I of bright spot of image and dim spot among the identical ROI respectively
2, I
1, by the uniformity of following formula computed image:
(4) the pseudo-shadow of ghost; The pseudo-shadow of bigger ghost inevitably can appear in the fMRI image sequence, and the signal that collects is impacted.The size of the pseudo-shadow of strict control ghost also is the importance of fMRI QA.In order to measure Ghost-to-signal ratio (GSR), the image of cutting apart in need measuring SNR carries out image segmentation once more, concrete grammar is: to the background of image segmentation gained image, adopt the iteration threshold method to cut apart once more, background after cutting apart is a noise region, and prospect is territory, the pseudo-shadow zone of ghost; Select suitable area-of-interest (ROI) in background and foreground area, as noise signal and the pseudo-shadow signal of ghost, the area-of-interest of water mould part is pressed the GSR of following formula computed image as signaling zone:
GSR＝(mean?of?ghost?ROI-mean?noise?ROI)/(mean?signal?ROI)
(5) signal drift: the fMRI sequence is the time dependent Dynamic Display of signal, in time series between image and image the stability of imaging seem very important.Therefore, its controlled factordisease also has signal drift except conventional factor.Signal drift is the important indicator of reflection stability.Because the unstability of MR imaging apparatus performance finally all can show as the unstability of gradation of image.Therefore the change procedure that detects same image-region signal is a kind of method of simple evaluation MR imaging apparatus stability.Obtain the central area of different time points moisture film image, calculate the drift of its image sequence:
S
Max, S
MinBe maximum, the minima on the signal time sequence of image sequence selection area; Come its variation tendency of Fitting Analysis with conic section.Match value and actual value are subtracted each other, and its difference obtains its frequency spectrum with fast fourier transform.In frequency spectrum, search for maximum, its range value and frequency are write the data base.
If if the PatientName label value is " mriDQA ", then be judged as common MRI data, then the QA indication information of Ce Lianging comprises: signal to noise ratio, the linearity, uniformity degree of image and signal drift.Wherein the metering system of uniformity degree of image and signal drift is identical with fMRI, introduces signal to noise ratio and the linearity measure method of common MRI here:
Signal to noise ratio: with repeatedly being with the threshold method split image, choose the signaling zone area-of-interest in prospect, method is identical with the fMRI choosing method; To background area, do not need to consider the pseudo-shadow of ghost this moment, the area-of-interest of noise signal is selected in zone on four angles: be radius, be that circle is drawn in the center of circle with the picture centre with half of the view picture figure image length of side, with drawn circle and intersection of diagonal is the center of circle, select suitable radius to determine area-of-interest, making its area is 2% of view picture image.
The linearity: the detection of the linearity needs that the known scaling point of a series of distances is arranged in the water mould, and the present invention utilizes the known dimensions of water mould to calibrate.The quality of image linearity is generally represented with the distortion percentage rate.The big more image linearity that distorts is poor more, and the more little image linearity that distorts is good more.The linearity computing formula of image is: L
F=(X-I)/I * 100%, L
P=(Y-I)/I * 100%; Wherein I is a physical length; L
FBe the linearity on the frequency coding direction, X is the length on the frequency coding direction; L
PThe linearity on the phase-encoding direction, Y is the length on phase-encoding direction.
After finishing above QA index calculating, result of calculation is write the data base preserve.Start relatively working of measurement result simultaneously:
(1) whether the pseudo-shadow of the image under consideration uniformity (PIU), signal drift (shift) and ghost (GSR) is at allowed band (PIU＞82%, shift＜1%, GSR＜3%);
(2) relatively, whether investigate set point (generally be set at difference and be no more than standard deviation) visual signal to noise ratio, torsion resistance and history average;
(3) frequency spectrum after signal drift is handled, whether the difference of investigating maximum and amplitude meansigma methods is 1.5 times of standard deviation scopes.
In the present embodiment, the common MRI of scanning water mould records signal to noise ratio 628 by above step, uniformity degree of image 91.70%, linearity L
F=0.87%, L
P=0.74%, signal drift 0.59%, all results are in allowed band.
In result of calculation was handled, any one index went beyond the scope, and then ejected the dialog box prompting, specification exception index item in the prompted dialog frame.
Claims (1)
1. the measuring method of NMR (Nuclear Magnetic Resonance)-imaging stabilization of equipment performance and imaging index is characterized in that this method may further comprise the steps:
(1) utilize the DICOM agreement to finish the process control of nuclear magnetic resonance, NMR quality assurance
1. server end keeps accepting state at any time in the mode of listening port, when receiving the association request of client, accepts related request and judges whether COS is storage class, if then accept, otherwise refuses and continue to monitor;
2. receive message data,, judge whether that if then discharge association, and execution in step 3. in order to discharge related request, otherwise execution in step 4. the message data that receives;
3. when judging message data when discharging related request, discharge related and according to the view data that judges whether to the quality assurance data, if execution in step (3) then and continues to monitor, otherwise finish relatedly, continue to keep listening state;
4. the message data when reception is not when discharging related request, execution in step (2), and repeating step is 2.;
(2) extraction data and write the data base
1. resolve the packet that receives, judge whether to accept data,, otherwise finish if then continue next step according to the DICOM label;
2. judge whether to write the data base according to the DICOM label,, otherwise finish if then continue next step;
3. in the data base, set up new record, write date, scan protocols, mid frequency and amplification, and the data storage path is write the data base;
(3) NMR (Nuclear Magnetic Resonance)-imaging quality assurance index analysis is measured
When server end receives the related request of release, discharge related also start-up routine and begin the detection and the calculating of nuclear magnetic resonance, NMR quality assurance index, the steps include:
1. judge whether to be functional MRI data according to label, if then calculate nuclear magnetic resonance, NMR quality assurance index as follows;
A computed image signal to noise ratio; With repeatedly being with threshold method to cut apart image segmentation is background area and prospect, selects the area-of-interest of signaling zone in prospect; The noise range area-of-interest with repeatedly being with threshold method to cut apart once more, is selected at background area in the noise range, extract the standard deviation of signaling zone area-of-interest average gray value and noise range area-of-interest, that calculates the two obtains signal noise ratio (snr) of image than promptly;
The uniformity of b measurement image; To cutting apart the signaling zone of back image, in selected area-of-interest, obtain the signal value of bright spot of image and dim spot, by formula the uniformity of computed image;
C measurement image torsion resistance; With repeatedly being with the threshold method split image, the prospect of image is asked its maximum gauge L with the varied angle searching algorithm
Max, minimum diameter L
MinAnd average diameter L
m, the torsion resistance D of computed image:
D＝(L
max-L
min)/L
m；
The pseudo-shadow of the ghost of d image; With repeatedly being with the threshold method split image, image segmentation is once more carried out in the noise range of split image, prospect is the pseudo-shadow zone of ghost, background is the noise range, can try to achieve the pseudo-shadow of ghost of image:
The pseudo-shadow of the ghost of image=(the pseudo-shadow meansigma methods of ghost-noise range meansigma methods)/signaling zone meansigma methods
The signal drift of e sequence of computed images; A selection area to temporal sequence of images writes down its signal value one by one, represents its variation tendency with matched curve, and calculates drift value shift:
S
Max, S
MinMaximum, minima for selection area signal in the temporal sequence of images; Match value and actual value are subtracted each other, and its difference obtains its frequency spectrum with fast fourier transform, searches for maximum in frequency spectrum, and its range value and frequency are write the data base;
2. if judge it is not functional MRI data, then calculate nuclear magnetic resonance, NMR quality assurance index according to the following steps:
A computed image signal to noise ratio; In image sequence, select piece image, with repeatedly being with threshold method to cut apart image segmentation is background and prospect, in background and prospect, select the area-of-interest of noise range and signaling zone, extract the standard deviation of signal area-of-interest average gray value and noise area-of-interest, that calculates the two obtains signal noise ratio (snr) of image than promptly;
The uniformity of b measurement image; To cutting apart the signaling zone of back image, in selected area-of-interest, obtain the signal value of bright spot of image and dim spot, according to the uniformity of formula computed image;
The linearity of c measurement image; In the known water mould under the prerequisite of size,, calculate the linearity of image according to water mould length on respective frequencies and the phase-encoding direction in image;
3. the nuclear magnetic resonance, NMR quality assurance indication information that calculates is write the data base;
4. with result of calculation and normal data and historical average according to contrast, going beyond the scope then provides prompting.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2008100456644A CN100571620C (en) | 2008-07-29 | 2008-07-29 | The measuring method of a kind of NMR (Nuclear Magnetic Resonance)-imaging stabilization of equipment performance and imaging index |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2008100456644A CN100571620C (en) | 2008-07-29 | 2008-07-29 | The measuring method of a kind of NMR (Nuclear Magnetic Resonance)-imaging stabilization of equipment performance and imaging index |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101322648A true CN101322648A (en) | 2008-12-17 |
| CN100571620C CN100571620C (en) | 2009-12-23 |
Family
ID=40186459
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB2008100456644A Expired - Fee Related CN100571620C (en) | 2008-07-29 | 2008-07-29 | The measuring method of a kind of NMR (Nuclear Magnetic Resonance)-imaging stabilization of equipment performance and imaging index |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN100571620C (en) |
Cited By (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103366082A (en) * | 2012-03-29 | 2013-10-23 | 通用电气公司 | System and method of managing technician review of medical test data |
| CN104161516A (en) * | 2014-01-09 | 2014-11-26 | 上海联影医疗科技有限公司 | Magnetic resonance imaging direction judging method and device |
| CN104700440A (en) * | 2015-03-17 | 2015-06-10 | 上海理工大学 | Magnetic resonance partial K-space image reconstruction method |
| CN105139410A (en) * | 2015-09-11 | 2015-12-09 | 浙江工商大学 | Brain tumor MRI image segmentation algorithm utilizing space section projection |
| CN105158818A (en) * | 2015-08-31 | 2015-12-16 | 合肥美亚光电技术股份有限公司 | Signal evaluation method of X-ray detection device |
| CN106575118A (en) * | 2014-08-22 | 2017-04-19 | 西门子股份公司 | Detection of consistent steady-state behavior in process plants |
| CN106657991A (en) * | 2017-01-25 | 2017-05-10 | 泰山医学院 | Magnetic resonance quality control method and system, and server |
| CN106859679A (en) * | 2017-01-25 | 2017-06-20 | 泰山医学院 | A kind of PET MR method for evaluating quality, server, terminal and system |
| CN108169081A (en) * | 2017-12-14 | 2018-06-15 | 四川大学华西医院 | The difference of blood cell analysis verifies model and its application process |
| CN109073723A (en) * | 2016-04-21 | 2018-12-21 | 皇家飞利浦有限公司 | Usage history database update MRI pulse train parameters |
| CN109801258A (en) * | 2018-12-18 | 2019-05-24 | 重庆市忠县人民医院 | A kind of medical imaging diagnosis quality control system and method |
| CN109830286A (en) * | 2019-02-13 | 2019-05-31 | 四川大学 | Brain function magnetic resonance code energy imaging method based on nonparametric statistics |
| CN109959886A (en) * | 2017-12-22 | 2019-07-02 | 西门子保健有限责任公司 | For determining the method, apparatus of the image quality information of MR imaging apparatus |
| CN110322408A (en) * | 2019-06-11 | 2019-10-11 | 浙江大学 | Multicenter magnetic resonance image automated quality control method based on cloud platform |
| CN110346397A (en) * | 2019-08-05 | 2019-10-18 | 固安县朝阳生物科技有限公司 | Nuclear magnetic resonance multimode detects die body |
| CN110850344A (en) * | 2018-08-21 | 2020-02-28 | 西门子医疗有限公司 | Method of operating an MRI apparatus |
| CN112690809A (en) * | 2020-02-04 | 2021-04-23 | 首都医科大学附属北京友谊医院 | Method, device, server and storage medium for determining equipment abnormality reason |
| CN113985334A (en) * | 2021-11-08 | 2022-01-28 | 电子科技大学 | Method for evaluating signal-to-noise ratio of magnetic resonance scanning image |
| CN117437321A (en) * | 2023-12-20 | 2024-01-23 | 清华大学 | Automatic evaluation method and device for arterial transit artifact through magnetic resonance spin labeling imaging |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6694047B1 (en) * | 1999-07-15 | 2004-02-17 | General Electric Company | Method and apparatus for automated image quality evaluation of X-ray systems using any of multiple phantoms |
| US6771067B2 (en) * | 2001-04-03 | 2004-08-03 | The United States Of America As Represented By The Department Of Health And Human Services | Ghost artifact cancellation using phased array processing |
| EP1728089A1 (en) * | 2004-03-12 | 2006-12-06 | Koninklijke Philips Electronics N.V. | Prescan for optimization of mri scan parameters |
| CN1799498B (en) * | 2004-12-31 | 2010-04-28 | 西门子（中国）有限公司 | Quick generalized partially parallel acquisition(GRAPA) image reconstruction algorithm for magnetic resonance imaging(MRI) |
| US7532942B2 (en) * | 2006-01-30 | 2009-05-12 | Bruce Reiner | Method and apparatus for generating a technologist quality assurance scorecard |
-
2008
- 2008-07-29 CN CNB2008100456644A patent/CN100571620C/en not_active Expired - Fee Related
Cited By (30)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103366082B (en) * | 2012-03-29 | 2018-01-26 | 通用电气公司 | The system and method that managerial and technical staff examines medical test data |
| CN103366082A (en) * | 2012-03-29 | 2013-10-23 | 通用电气公司 | System and method of managing technician review of medical test data |
| CN104161516A (en) * | 2014-01-09 | 2014-11-26 | 上海联影医疗科技有限公司 | Magnetic resonance imaging direction judging method and device |
| CN104161516B (en) * | 2014-01-09 | 2015-09-02 | 上海联影医疗科技有限公司 | Nuclear magnetic resonance orientation determination methods and device thereof |
| CN106575118A (en) * | 2014-08-22 | 2017-04-19 | 西门子股份公司 | Detection of consistent steady-state behavior in process plants |
| CN104700440A (en) * | 2015-03-17 | 2015-06-10 | 上海理工大学 | Magnetic resonance partial K-space image reconstruction method |
| CN104700440B (en) * | 2015-03-17 | 2017-10-13 | 上海理工大学 | Magnetic resonant part K spatial image reconstruction method |
| CN105158818A (en) * | 2015-08-31 | 2015-12-16 | 合肥美亚光电技术股份有限公司 | Signal evaluation method of X-ray detection device |
| CN105158818B (en) * | 2015-08-31 | 2017-10-10 | 合肥美亚光电技术股份有限公司 | A kind of method of X-ray detection equipment signal evaluation |
| CN105139410A (en) * | 2015-09-11 | 2015-12-09 | 浙江工商大学 | Brain tumor MRI image segmentation algorithm utilizing space section projection |
| CN105139410B (en) * | 2015-09-11 | 2017-11-10 | 浙江工商大学 | The brain tumor MRI image dividing method projected using aerial cross sectional |
| CN109073723A (en) * | 2016-04-21 | 2018-12-21 | 皇家飞利浦有限公司 | Usage history database update MRI pulse train parameters |
| WO2018137311A1 (en) * | 2017-01-25 | 2018-08-02 | 泰山医学院 | Magnetic resonance quality control method, server and system |
| CN106859679A (en) * | 2017-01-25 | 2017-06-20 | 泰山医学院 | A kind of PET MR method for evaluating quality, server, terminal and system |
| CN106657991A (en) * | 2017-01-25 | 2017-05-10 | 泰山医学院 | Magnetic resonance quality control method and system, and server |
| CN106859679B (en) * | 2017-01-25 | 2023-11-03 | 泰山医学院 | PET-MR quality evaluation method, server, terminal and system |
| CN108169081A (en) * | 2017-12-14 | 2018-06-15 | 四川大学华西医院 | The difference of blood cell analysis verifies model and its application process |
| CN109959886A (en) * | 2017-12-22 | 2019-07-02 | 西门子保健有限责任公司 | For determining the method, apparatus of the image quality information of MR imaging apparatus |
| CN109959886B (en) * | 2017-12-22 | 2021-09-10 | 西门子保健有限责任公司 | Method and device for determining imaging quality information of magnetic resonance imaging device |
| CN110850344A (en) * | 2018-08-21 | 2020-02-28 | 西门子医疗有限公司 | Method of operating an MRI apparatus |
| CN109801258A (en) * | 2018-12-18 | 2019-05-24 | 重庆市忠县人民医院 | A kind of medical imaging diagnosis quality control system and method |
| CN109830286A (en) * | 2019-02-13 | 2019-05-31 | 四川大学 | Brain function magnetic resonance code energy imaging method based on nonparametric statistics |
| CN109830286B (en) * | 2019-02-13 | 2022-09-30 | 四川大学 | Brain function magnetic resonance encoding energy imaging method based on nonparametric statistics |
| CN110322408A (en) * | 2019-06-11 | 2019-10-11 | 浙江大学 | Multicenter magnetic resonance image automated quality control method based on cloud platform |
| CN110346397A (en) * | 2019-08-05 | 2019-10-18 | 固安县朝阳生物科技有限公司 | Nuclear magnetic resonance multimode detects die body |
| CN112690809A (en) * | 2020-02-04 | 2021-04-23 | 首都医科大学附属北京友谊医院 | Method, device, server and storage medium for determining equipment abnormality reason |
| CN112690809B (en) * | 2020-02-04 | 2021-09-24 | 首都医科大学附属北京友谊医院 | Method, device, server and storage medium for determining equipment abnormality reason |
| CN113985334A (en) * | 2021-11-08 | 2022-01-28 | 电子科技大学 | Method for evaluating signal-to-noise ratio of magnetic resonance scanning image |
| CN117437321A (en) * | 2023-12-20 | 2024-01-23 | 清华大学 | Automatic evaluation method and device for arterial transit artifact through magnetic resonance spin labeling imaging |
| CN117437321B (en) * | 2023-12-20 | 2024-03-12 | 清华大学 | Automatic evaluation method and device for arterial transit artifact through magnetic resonance spin labeling imaging |
Also Published As
| Publication number | Publication date |
|---|---|
| CN100571620C (en) | 2009-12-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN100571620C (en) | The measuring method of a kind of NMR (Nuclear Magnetic Resonance)-imaging stabilization of equipment performance and imaging index | |
| US20200146648A1 (en) | Methods and apparatus for extended low contrast detectability for radiographic imaging systems | |
| CN101711670B (en) | Method for determining a kidney function parameter | |
| Mento et al. | On the influence of imaging parameters on lung ultrasound B-line artifacts, in vitro study | |
| Ruano et al. | Fetal lung volume estimated by 3‐dimensional ultrasonography and magnetic resonance imaging in cases with isolated congenital diaphragmatic hernia | |
| WO2018137311A1 (en) | Magnetic resonance quality control method, server and system | |
| WO2003063702A1 (en) | Magnetic resonance imaging method and apparatus | |
| EP3502727B1 (en) | Method for determining an imaging quality information for a magnetic resonance imaging apparatus from a signal-to-noise ratio at an anatomical landmark | |
| Alizadeh Dehnavi et al. | Assessment of the carotid artery by MRI at 3T: a study on reproducibility | |
| CN103857331A (en) | Method and system for quantifying hepatic fat in humans | |
| US10345411B2 (en) | Frequency monitoring of gradient pulses during magnetic resonance imaging | |
| US9636085B2 (en) | Method for determining, in real time, the probability that target biological tissue is opposite an ultrasonic transducer | |
| Wear et al. | US backscatter for liver fat quantification: an AIUM-RSNA QIBA pulse-echo quantitative ultrasound initiative | |
| Medvedofsky et al. | Multi-parametric quantification of tricuspid regurgitation using cardiovascular magnetic resonance: A comparison to echocardiography | |
| Di Marcello et al. | Assessment of mitral regurgitation severity by Doppler color flow mapping of the vena contracta in dogs | |
| CN110916663B (en) | Portable nuclear magnetic resonance organ elasticity noninvasive quantitative detection method | |
| CN106859679B (en) | PET-MR quality evaluation method, server, terminal and system | |
| US20230084562A1 (en) | Non-destructive inspection method and system based on artificial intelligence | |
| CN206517540U (en) | A kind of magnetic resonance quality control system | |
| Charan Lanjewar et al. | Planimetry of mitral valve stenosis in rheumatic heart disease by magnetic resonance imaging | |
| Fang et al. | Apparent diffusion coefficient reproducibility across 3 T scanners in a breast diffusion phantom | |
| Cheng et al. | Usefulness of contrast kinetics for predicting and monitoring tissue changes in muscle following thermal therapy in long survival studies | |
| Bonny et al. | Uncertainties of calculated Cramér‐Rao lower bounds: implications for quantitative MRS | |
| US7415371B2 (en) | Systems and methods for proactive detection of imaging chain problems during normal system operation | |
| CN113261940B (en) | Method and device for detecting magnetic resonance image quality |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20091223 Termination date: 20170729 |