CN104224175B - Method of fusing two-dimensional magnetic resonance spectrum and three-dimensional magnetic resonance navigation image - Google Patents

Abstract

The invention belongs to the field of medical image postprocessing and relates to a method of fusing a two-dimensional magnetic resonance spectrum and a three-dimensional magnetic resonance navigation image. According to the method, three-dimensional magnetic resonance navigation information fused with a spectrum image is generated by means of spatial information and image registration by the aid of spatial information fields defined in header files of magnetic resonance original data and a non-rigid registration algorithm. The method has the advantages that the problem that the two-dimensional magnetic resonance spectrum image and the three-dimensional magnetic resonance navigation image cannot be fused is solved, metabolic information is increased for navigation, and the method is simple and convenient to operate, short in consumed time (20 minutes averagely), stable, reliable, free of extra introduction of hardware facilities and convenient for clinical practical application.

Description

A kind of method merging twodimensional magnetic resonance wave spectrum and three-dimensional magnetic resonance navigation image

Technical field

The invention belongs to medical image post processing field, it is related to a kind of twodimensional magnetic resonance wave spectrum and three-dimensional magnetic resonance navigation shadow Method as merging；The method is applied to all brain structures image and the post processing of metabolism image.

Background technology

The appearance of Use of Neuronavigation technology is the milestone in neurosurgery development history, and it makes neurosurgery from simple Rely on patient's subjective experience and enter into the epoch that objective image guides, greatly improve operation safety.Use of Neuronavigation is the most often adopted Guide image is nuclear magnetic resonance image, and three-dimensional nuclear magnetic resonance image can comprehensive display brain anatomy structure；But regrettably, Current three-dimensional magnetic resonance navigation image does not also include Magnetic Resonance Spectrum image.Described Magnetic Resonance Spectrum is for diagnostic imaging Mature technology, it is possible to provide the metabolic information of cerebral tissue or pathological changes, differentiating tumor property, grade malignancy, invading profit scope etc. side Face has substantial worth, if Magnetic Resonance Spectrum visual fusion can be entered navigation directions operation, will further improve operation safety.

The main cause that described twodimensional magnetic resonance wave spectrum image and three-dimensional magnetic resonance navigation image cannot merge is the two lattice Formula is different, and wave spectrum image is .rda form, and pixel format used by navigation sequence is DICOM format, there is no at present and merges the two Software.And in process of clinical application, spectral information can only be reduced in three-dimensional navigation image with subjective memory, method is coarse, Error is larger, and credibility is poor.

Therefore, how by Imaging processing method, twodimensional magnetic resonance wave spectrum accurately to be melted with three-dimensional magnetic resonance navigation image Close, be clinical at present demand urgently to be resolved hurrily.

Content of the invention

The purpose of the present invention is defect or the deficiency overcoming prior art, provides a kind of twodimensional magnetic resonance wave spectrum and Three-Dimensional Magnetic The method that resonance navigation image merges；The method adopts spatial information registration and non-rigid surface adjustment fusion evaluation, solves two Tie up the problem that Magnetic Resonance Spectrum image and three magnetic resonance navigation images can not merge, increase metabolic information, operation letter for navigation Just, time-consuming short, it is easy to Clinical practice.

The method that the twodimensional magnetic resonance wave spectrum of the present invention is merged with three-dimensional magnetic resonance navigation image is it is characterised in that adopt Spatial information field defined in the header file of magnetic resonance raw data and non-rigid surface adjustment algorithm, carry out spatial information and shadow As registration, generate the magnetic resonance three-dimensional navigation information merging wave spectrum image, concretely comprise the following steps：

（1）Obtain twodimensional magnetic resonance initial data, preliminary post processing carried out by the business software that magnetic resonance equipment carries, With Cho/NAA for value interested, generate twodimensional magnetic resonance wave spectrum image；

（2）Obtain three-dimensional magnetic resonance navigation DICOM data, with twodimensional magnetic resonance wave spectrum image as reference, with intracerebral lesion Centered on area, define region of interest；

（3）Above-mentioned data is imported work station, carries out wave spectrum image registering with the spatial information of navigation image；Two-dimensional magnetic is altogether Defined in the header file of vibration wave spectrum initial data, the field of spatial information is mainly:SliceThickness (thickness), PositionVector（Position vector）,RowVector（Row vector）,ColumnVector（Column vector）, PixelSpacingRow（Pixel between-line spacing）,PixelSpacingCol（Pixel column is spaced）；Three-dimensional magnetic resonance navigation image Defined in DICOM image header file, the field of spatial information is mainly:SliceThickness (thickness), ImagePostionPatient（Patient image position）,ImageOrientationPatient（Patient image positions）, PixelSpacing（Pixel separation）；Transformational relation between this two groups of data is（As shown in Figure 2）：

SliceThickness(DICOM)=SliceThickness(rda)

ImagePostionPatient(DICOM)=PositionVector(rda)

ImageOrientationPatient(DICOM)=[RowVector,ColumnVector](rda)

PixelSpacing(DICOM)=[PixelSpacingRow,PixelSpacingCol]；

（4）Further non-rigid method Image registration, specific algorithm is：

As with reference to image R (i, j, k) (i, j, k are body bundle index), its characteristic point is set to three-dimensional navigation imageAs floating image F (i, j, k) (i, j, k are body bundle index), its characteristic point is set to two-dimentional wave spectrum imageDeformation field U (x | D)；Due to point set T and S concordance it needs to be determined that, the energy function of therefore design is,

J (R, F, S, T | D)=IntensityMatching (R, F (U (P | D)))

+λ×PointMatching(T,U(S|D)

For measuring the similarity of above-mentioned function, using cross-correlation（Cross-Correlation）As module,

$\frac{{\mathrm{\Σ}}_{i,i,k}(R(i,j,k)-\stackrel{\‾}{R})(F\left(U\left(p_{\mathrm{ijk}}\right|D)\right)-\stackrel{\‾}{F})}{\sqrt{{\mathrm{\Σ}}_{i,j,k}{(R(i,j,k)-\stackrel{\‾}{R})}^2{\mathrm{\Σ}}_{i,j,k}{\left(F\left(U\left(p_{\mathrm{ijk}}\right|D)\right)\stackrel{\‾}{F}\right)}^2}}$ ；

（5）By setting alpha value, adjust the pixel weight of two images, if Alpha is 1, show navigation picture completely Picture；If 0.5, respectively show 50%；If 0, only display wave spectrum；Can arbitrarily adjust as needed；

（6）If stage construction two dimension wave spectrogram, repeat the above steps；

（7）The image that registration is merged ultimately generates DICOM file, imports navigator, is applied to clinic.

Compared with prior art, its advantage is the present invention：Solve twodimensional magnetic resonance wave spectrum image and three magnetic resonance The problem that navigation image can not merge, easy and simple to handle, take short（Average 20 minutes）, reliable and stable it is not necessary to additionally introduce hard Part facility, is easy to clinical practice.

Brief description

The present invention is further described with reference to the accompanying drawings and examples.

Fig. 1 is the overall flow chart of steps of the method for the invention.

Fig. 2 is twodimensional magnetic resonance wave spectrogram and three magnetic resonance navigation image data conversion relation schematic diagrams, wherein,

（1）Represent row pixel separation,（2）Represent row pixel separation.

Fig. 3 shows that twodimensional magnetic resonance wave spectrum selects region of interest and definition space information.

Twodimensional magnetic resonance wave spectrum interested position is merged after showing registration by Fig. 4 with three-dimensional magnetic resonance navigation data, and Realize three dimensional display, wherein, highlighted square is wave spectrum target spot interested.

Specific embodiment

Embodiment 1

The method that this twodimensional magnetic resonance wave spectrum is merged with three-dimensional magnetic resonance navigation image, concretely comprises the following steps：

（1）Obtain Magnetic Resonance Spectrum and magnetic resonance navigation picture initial data, preliminary treatment generates two dimension and 3-D view；

（2）Region of interest is set respectively on twodimensional magnetic resonance wave spectrum image and three-dimensional magnetic resonance navigation image；

（3）Carry out wave spectrum image registering with the spatial information of navigation image；

（4）Carry out non-rigid method and carry out the fusion registering with three-dimensional magnetic resonance navigation image of twodimensional magnetic resonance wave spectrum image, Generate the DICOM file merging, export to navigator and be used for clinic；

（5）Two kinds of image weights can be adjusted as needed at any time, reach optimum display effect.

Embodiment 2

The present invention adopts spatial information field and the calculation of non-rigid surface adjustment defined in the header file of magnetic resonance raw data Method, carries out spatial information and Image registration, generates the magnetic resonance three-dimensional navigation information merging wave spectrum image, concretely comprises the following steps：

（1）Obtain twodimensional magnetic resonance initial data, preliminary post processing carried out by the business software that magnetic resonance equipment carries, With Cho/NAA for value interested, generate twodimensional magnetic resonance wave spectrum image；

（2）Obtain three-dimensional magnetic resonance navigation DICOM data, with twodimensional magnetic resonance wave spectrum image as reference, with intracerebral lesion Centered on area, define region of interest；

（3）Above-mentioned data is imported work station, carries out wave spectrum image registering with the spatial information of navigation image；Two-dimensional magnetic is altogether Defined in the header file of vibration wave spectrum initial data, the field of spatial information is mainly:SliceThickness (thickness), PositionVector（Position vector）,RowVector（Row vector）,ColumnVector（Column vector）, PixelSpacingRow（Pixel between-line spacing）,PixelSpacingCol（Pixel column is spaced）；Three-dimensional magnetic resonance navigation image Defined in DICOM image header file, the field of spatial information is mainly:SliceThickness (thickness), ImagePostionPatient（Patient image position）,ImageOrientationPatient（Patient image positions）, PixelSpacing（Pixel separation）；Transformational relation between this two groups of data is（As shown in Figure 2）：

SliceThickness(DICOM)=SliceThickness(rda)

ImagePostionPatient(DICOM)=PositionVector(rda)

ImageOrientationPatient(DICOM)=[RowVector,ColumnVector](rda)

PixelSpacing(DICOM)=[PixelSpacingRow,PixelSpacingCol]；

（4）Further non-rigid method Image registration, specific algorithm is：

As with reference to image R (i, j, k) (i, j, k are body bundle index), its characteristic point is set to three-dimensional navigation imageAs floating image F (i, j, k) (i, j, k are body bundle index), its characteristic point is set to two-dimentional wave spectrum imageDeformation field U (x | D)；Due to point set T and S concordance it needs to be determined that, the energy function of therefore design is,

J (R, F, S, T | D)=IntensityMatching (R, F (U (P | D)))

+λ×PointMatching(T,U(S|D)

For measuring the similarity of above-mentioned function, using cross-correlation（Cross-Correlation）As module,

$\frac{{\mathrm{\Σ}}_{i,i,k}(R(i,j,k)-\stackrel{\‾}{R})(F\left(U\left(p_{\mathrm{ijk}}\right|D)\right)-\stackrel{\‾}{F})}{\sqrt{{\mathrm{\Σ}}_{i,j,k}{(R(i,j,k)-\stackrel{\‾}{R})}^2{\mathrm{\Σ}}_{i,j,k}{\left(F\left(U\left(p_{\mathrm{ijk}}\right|D)\right)\stackrel{\‾}{F}\right)}^2}}$ ；

（5）By setting alpha value, adjust the pixel weight of two images, if Alpha is 1, show navigation picture completely Picture；If 0.5, respectively show 50%；If 0, only display wave spectrum；Can arbitrarily adjust as needed；

（6）If stage construction two dimension wave spectrogram, repeat the above steps；

（7）The image that registration is merged ultimately generates DICOM file, imports navigator, application.

Claims (2)

1. a kind of fusion twodimensional magnetic resonance wave spectrum and the method for three-dimensional magnetic resonance navigation image are it is characterised in that adopt magnetic resonance Spatial information field defined in the header file of initial data and non-rigid surface adjustment algorithm, carry out spatial information and image are joined Standard, generates the magnetic resonance three-dimensional navigation information merging wave spectrum image, it includes step：

(1) obtain twodimensional magnetic resonance initial data, preliminary post processing is carried out by the business software of magnetic resonance equipment, with Cho/ NAA is value interested, generates twodimensional magnetic resonance wave spectrum image；

(2) obtain three-dimensional magnetic resonance navigation DICOM data, with twodimensional magnetic resonance wave spectrum image as reference, with intracerebral lesion area be Center, defines region of interest；

(3) the three-dimensional magnetic resonance navigation DICOM data of the two dimension generating resonance wave spectrum image and acquisition is imported work station, carry out Wave spectrum image is registering with the spatial information of navigation image；Space defined in the header file of described twodimensional magnetic resonance wave spectrum initial data The field of information is mainly:SliceThickness,PositionVector,RowVector,ColumnVector, PixelSpacingRow,PixelSpacingCol；Defined in the DICOM image header file of described three-dimensional magnetic resonance navigation image The field of spatial information is mainly：

SliceThickness,ImagePostionPatient,ImageOrientationPatient,PixelSpacing；

(4) non-rigid method Image registration further；

(5) pass through to set alpha value, adjust the pixel weight of wave spectrum image and navigation image, if Alpha is 1, show completely Navigation image；If 0.5, then wave spectrum image each display 50% with navigation image；If 0, only display wave spectrum image；As needed Arbitrarily adjust the size of alpha value；

(6) if stage construction two dimension wave spectrogram, repeat the above steps；

(7) image merging registration ultimately generates DICOM file, imports navigator.

2. the method as described in claim 1 is it is characterised in that in described step (3), described twodimensional magnetic resonance wave spectrum original number According to the transformational relation between the DICOM image head file data of three-dimensional magnetic resonance navigation image it is：

SliceThickness (DICOM)=SliceThickness (rda)

ImagePostionPatient (DICOM)=PositionVector (rda)

ImageOrientationPatient (DICOM)=[RowVector, ColumnVector] (rda)

PixelSpacing (DICOM)=[PixelSpacingRow, PixelSpacingCol].