CN102074026A - Method for displaying vector graph on medicinal image - Google Patents

Abstract

The invention provides a method for displaying a vector graph on a medicinal image. According to the relative change of boundary rectangles of the medicinal image and the vector graph, the method enables the vector graph to change together with the medicinal image. After the method is used, a medicinal digital imaging and communications in medicine (DICOM) image and an artificial joint industrial template file of a graphics interchange format (DXF) in accordance with industrial standards can be displayed on the same coordinate system through a doctor, and accurate matching of type and size is performed by the operations such as movement and rotation of an artificial joint graph and the like. Therefore, the method is convenient and quick, is greatly convenient for the doctor to select the type of the artificial joint, and plays a role in assisting and planning the operation of the doctor.

Description

A kind of method that on medical image, shows vector graphics

Technical field

The invention provides a kind of method that on medical image, shows vector graphics, relate in particular to and a kind ofly on medical image, show vector graphics and can move method with convergent-divergent together.

Background technology

The definition of vector graphics: vector image, be also referred to as OO image or drawing image, on mathematics, be defined as a series of points that connect by line.Graphic element in the vector file is called object.Each object all is an entity that has a style of one's own, and it has attributes such as color, shape, profile, size and screen position.Since each object all is an entity that has a style of one's own, just can when keeping its original sharpness and flexibility, repeatedly moves and change its attribute, and can not influence other object in the legend.These features make the program based on vector be specially adapted to legend and three-dimensional modeling, because they require to create and to operate single object usually.Drawing based on vector is irrelevant with resolution.This means that they can be shown on the output device by highest resolution.

The file layout of polar plot has SWF, SVG, WMF, EMF, EPS, DXF totally six kinds of file layouts, this paper example uses the DXF file layout, the DXF file layout is the vector file form among the AutoCAD, it is the actual industrial standard of polar plot form now, it is with ASCII character mode storage file (format standard of redaction is supported simultaneously with the binary mode storage file), and is very accurate aspect performance figure big or small.The DXF file can be by many software transfers or output, and the expansion of its file is called " .dxf ".

The problem that on medical image, shows vector graphics.Because the medical science DICOM image of standard is a vector graphics, and the joint prosthesis of orthopaedics also is the vector graphics of the industrial DXF form of standard.If two kinds of images are simultaneously displayed on a coordinate system the inside, just can help the doctor to determine that before operation operation needs the size and the model of the joint prosthesis of use, accelerate the speed of operation, shorten the time of operation, reduce patient simultaneously because the sum of errors risk that the displacement joint prosthesis brings.

When carrying out joint prosthesis implantation or replacement operation, the doctor often will push bloody musculature aside, and dependence experience is then judged the size and the model of required joint prosthesis.The problem that such operation brings is exactly that the doctor has extended the time of performing the operation in the selection of constantly carrying out joint prosthesis model and size, has also increased the risk of operation simultaneously.In order to reduce this time of by virtue of experience changing joint prosthesis, the orthopaedics DICOM image of taking before patient can being performed the operation is put together with the joint prosthesis vector graphics that meets the DXF form of industrial standard and is shown, helps the doctor to determine the size and the model of joint prosthesis accurately before operation.

Existing and the most akin implementation of the present invention: the doctor becomes film with patient's orthopaedics DICOM image print, is placed on the special-purpose lamp box, uses size and model that survey instruments such as ruler and set square determine to need to replace joint prosthesis simultaneously.

The shortcoming of prior art is conspicuous, doctor's trouble that operates not only, but also need to calculate a lot of parameters.

Summary of the invention

The invention provides a kind of method that on medical image, shows vector graphics, after using method of the present invention, the doctor can be presented under the same coordinate system by the joint prosthesis industry template file with medical science DICOM image and the DXF form that meets industrial standard, and by the operations such as mobile and rotation of joint prosthesis figure being carried out the coupling of accurate model and size.So not only convenient but also quick, greatly facilitate the type selecting of doctor to joint prosthesis, doctor's operation has been played the effect of auxiliary and plan.

The present invention solves the problems of the technologies described above the technical scheme that is adopted to be:

A kind of method that shows vector graphics on medical image is characterized in that may further comprise the steps:

A. set up unified coordinate system;

B. in coordinate system, load medical image;

C. load vector graphics in coordinate system, vector graphics is presented on the medical image;

D. medical image pixel f (x _i, y _i) (0＜i≤N) coordinate in coordinate system is (x _i, y _i), wherein N is the number of pixel; Vector graphics is at the image pixel f of coordinate system ₁(x1 _i, y1 _i) (coordinate of 0＜i≤N1) is (x1 _i, y1 _i), wherein N1 is the number of pixel; According to the coordinate of medical image pixel, obtain the size of the border rectangle in medical image zone, it is Far Left left, rightmost right, top top and the bottommost bottom of the image-region of reference that the border rectangle has with coordinate; Four attributes of this border rectangle are obtained by following formula respectively:

left＝min{x ₀，x ₁，...，x _N-1}；

right＝max{x ₀，x ₁，...，x _N-1}；

top＝max{y ₀，y ₁，...，y _N-1}；

bottom＝min{y ₀，y ₁，...，y _N-1}；

Wherein N is the number of the image pixel in the rectangle of border, the border rectangle of medical image be (left, right, top, bottom), the border rectangle of vector graphics is (left ₁, right ₁, top ₁, bottom ₁);

E. the border rectangle of medical image obtains the center of this rectangle and is:

$(\frac{(\mathrm{left}+\mathrm{right})}{2},\frac{\mathrm{top}+\mathrm{bottom}}{2})$

The center that the border rectangle of vector graphics obtains this rectangle is:

$(\frac{({\mathrm{left}}_1+{\mathrm{right}}_1)}{2},\frac{{\mathrm{<figure-callout\; id="1"\; label="top"\; filenames="HSA00000298909500012.png"\; state="\{\{state\}\}">top</figure-callout>}}_1+{\mathrm{<figure-callout\; id="1"\; label="bottom"\; filenames="HSA00000298909500012.png"\; state="\{\{state\}\}">bottom</figure-callout>}}_1}{2})$

After medical image was moved, the moving method step of vector graphics was as follows:

E1. the distance that move in coordinate system relative to former rectangular centre at the center of obtaining the new border rectangle of medical image for (x, y),

E2. when the distance that moves is not 0, the border rectangle mobile phase of vector graphics with distance,

E3. in new border rectangle, show new vector graphics;

F. during the medical image convergent-divergent, the center of medical image border rectangle does not change, and medical image convergent-divergent rate is α, and the distance of the central point of the border rectangle of vector graphics and the border rectangle of medical image is:

$s={({(\frac{({\mathrm{left}}_1+{\mathrm{right}}_1)}{2}-\frac{(\mathrm{left}+\mathrm{right})}{2})}^2+{(\frac{\mathrm{top}+\mathrm{bottom}}{2}-\frac{{\mathrm{<figure-callout\; id="1"\; label="top"\; filenames="HSA00000298909500012.png"\; state="\{\{state\}\}">top</figure-callout>}}_1+{\mathrm{<figure-callout\; id="1"\; label="bottom"\; filenames="HSA00000298909500012.png"\; state="\{\{state\}\}">bottom</figure-callout>}}_1}{2})}^2)}^{1/2}$

New distance is after the convergent-divergent:

s ₁＝α*s

The angle β at the center of former vector graphics and the center of medical image is:

$\mathrm{\&beta;}=\mathrm{arctg}\left(\right|\frac{\frac{\mathrm{top}+\mathrm{bottom}}{2}-\frac{{\mathrm{<figure-callout\; id="1"\; label="top"\; filenames="HSA00000298909500012.png"\; state="\{\{state\}\}">top</figure-callout>}}_1+{\mathrm{<figure-callout\; id="1"\; label="bottom"\; filenames="HSA00000298909500012.png"\; state="\{\{state\}\}">bottom</figure-callout>}}_1}{2}}{\frac{{\mathrm{left}}_1+{\mathrm{right}}_1}{2}-\frac{\mathrm{left}+\mathrm{right}}{2}}\left|\right),$

The center of the border rectangle of then new vector graphics is:

$(\frac{(\mathrm{left}+\mathrm{right})}{2}-{\mathrm{<figure-callout\; id="1"\; label="s"\; filenames="HSA00000298909500012.png"\; state="\{\{state\}\}">s</figure-callout>}}_1*\cos \mathrm{\&beta;},\frac{\mathrm{top}+\mathrm{bottom}}{2}+{\mathrm{<figure-callout\; id="1"\; label="s"\; filenames="HSA00000298909500012.png"\; state="\{\{state\}\}">s</figure-callout>}}_1*\sin \mathrm{\&beta;})$

Can obtain new border rectangle after former vector graphics border rectangle calculates in the same way, draw new vector graphics in new border rectangle class.

Description of drawings

Fig. 1 is the border rectangle synoptic diagram of embodiment of the invention medical image in coordinate system;

Fig. 2 is the angle synoptic diagram at the center of the center of embodiment of the invention medical image and vector graphics;

Fig. 3 for the embodiment of the invention load vector graphics after main window show synoptic diagram;

Fig. 4 is an embodiment of the invention operating vector graphics class structural drawing;

Fig. 5 presses left mouse button for the embodiment of the invention and enters mobile status demonstration synoptic diagram in vector graphics;

Fig. 6 enters rotation status when pressing right mouse button for the embodiment of the invention in vector graphics;

Fig. 7 is an embodiment of the invention mouse roller amplification vector figure synoptic diagram;

Fig. 8 dwindles the polar plot synoptic diagram for embodiment of the invention mouse roller;

Fig. 9 is an embodiment of the invention operating process synoptic diagram.

Embodiment

With embodiment the present invention is described in further detail with reference to the accompanying drawings below:

The polar plot of DICOM image and DXF form be based upon with the current window be under the reference frame two not under the jurisdiction of each other with two width of cloth images that intersect.The DXF image is presented on the DICOM image.

In order to measure the size of bone, we are by loading polar plot, and the polar plot that loads moved and rotate mate with the bone that will undergo surgery.For the result who makes coupling is convenient to observe with more accurate, we have increased the same moved further and the synchronous zoom function of DICOM image and polar plot.

We suppose DICOM image pixel f (x _i, y _i) (0＜i≤N) coordinate in the client area is (x _i, y _i), wherein N is the number of pixel.Suppose that the DXF polar plot that loads is presented at the image pixel f of client area ₁(x1 _i, y1 _i) (coordinate of 0＜i≤N1) is (x1 _i, y1 _i), wherein N1 is the number of pixel.

According to the coordinate of image pixel, we can obtain the size of the border rectangle of image-region.The hypothetical boundary rectangle has left, right, top and four attributes of bottom, and representing respectively with the client area coordinate is Far Left, rightmost, top and the bottommost of the image-region of reference.As shown in Figure 1.

Four attributes of this border rectangle are obtained by following formula respectively:

left＝min{x ₀，x ₁，...，x _N-1}；

right＝max{x ₀，x ₁，...，x _N-1}；

top＝max{y ₀，y ₁，...，y _N-1}；

bottom＝min{y ₀，y ₁，...，y _N-1}；

Wherein N is the number of the image pixel in the rectangle of border.The border rectangle of DICOM image be (left, right, top, bottom), the border rectangle of DXF vector image is (left ₁, right ₁, top ₁, bottom ₁).

We by the center that the border rectangle of original DICOM image obtains this rectangle are:

$(\frac{(\mathrm{left}+\mathrm{right})}{2},\frac{\mathrm{top}+\mathrm{bottom}}{2})$

We by the center that the border rectangle of original DFX image obtains this rectangle are:

$(\frac{({\mathrm{left}}_1+{\mathrm{right}}_1)}{2},\frac{{\mathrm{top}}_1+{\mathrm{<figure-callout\; id="1"\; label="bottom"\; filenames="HSA00000298909500012.png"\; state="\{\{state\}\}">bottom</figure-callout>}}_1}{2})$

After mouse drag DICOM image is moved, the distance that move relative to former rectangular centre at the center of new border rectangle for (x, y).When we judge mobile distance is not after 0, the border rectangle mobile phase of DXF image distance (direction is also identical, and we distinguish direction with sign) together.By calculating, we show new DXF figure in new border rectangle.

In the time of image zoom, the center of image boundary rectangle does not change.Suppose that DICOM image zoom rate is α, the distance of the central point of the border rectangle of DFX image and the border rectangle of DICOM image is:

$s={({(\frac{({\mathrm{left}}_1+{\mathrm{right}}_1)}{2}-\frac{(\mathrm{left}+\mathrm{right})}{2})}^2+{(\frac{\mathrm{top}+\mathrm{bottom}}{2}-\frac{{\mathrm{<figure-callout\; id="1"\; label="top"\; filenames="HSA00000298909500012.png"\; state="\{\{state\}\}">top</figure-callout>}}_1+{\mathrm{<figure-callout\; id="1"\; label="bottom"\; filenames="HSA00000298909500012.png"\; state="\{\{state\}\}">bottom</figure-callout>}}_1}{2})}^2)}^{1/2}$

New distance is after the convergent-divergent: s ₁=α * s

As shown below, the angle β at the center of the center of former DXF image and DICOM image is:

$\mathrm{\&beta;}=\mathrm{arctg}\left(\right|\frac{\frac{\mathrm{top}+\mathrm{bottom}}{2}-\frac{{\mathrm{<figure-callout\; id="1"\; label="top"\; filenames="HSA00000298909500012.png"\; state="\{\{state\}\}">top</figure-callout>}}_1+{\mathrm{<figure-callout\; id="1"\; label="bottom"\; filenames="HSA00000298909500012.png"\; state="\{\{state\}\}">bottom</figure-callout>}}_1}{2}}{\frac{{\mathrm{left}}_1+{\mathrm{right}}_1}{2}-\frac{\mathrm{left}+\mathrm{right}}{2}}\left|\right)$

The center of the border rectangle of then new DXF image is:

$(\frac{(\mathrm{left}+\mathrm{right})}{2}-{\mathrm{<figure-callout\; id="1"\; label="s"\; filenames="HSA00000298909500012.png"\; state="\{\{state\}\}">s</figure-callout>}}_1*\cos \mathrm{\&beta;},\frac{\mathrm{top}+\mathrm{bottom}}{2}+{\mathrm{<figure-callout\; id="1"\; label="s"\; filenames="HSA00000298909500012.png"\; state="\{\{state\}\}">s</figure-callout>}}_1*\sin \mathrm{\&beta;})$

Can obtain new border rectangle after former DXF border rectangle calculates in the same way, draw new DXF image in new border rectangle class.Then we can realize the synchronous convergent-divergent of DXF image and DICOM image.

The orthopaedics image file is read in the internal memory from disk, in the main display window of screen, show the image that loads.In this process, calling the LoadFromFile function is loaded into the orthopaedics image the internal memory from disk, and the image that loading is come in is presented at the screen main window, when main window shows, employed function is the OnPaint among the VC, promptly in this function, the dc (contextual devices) that calls main window comes display image.After being loaded into internal memory from disk image, keep the pointer pImg of load image, use to treat the back.

The operating process of instantiation:

The image of needs is loaded into, be presented in the main window.

After showing the orthopaedics image, select corresponding orthopaedics template according to the different parts of health, this example uses the DXF file layout.

Selected with LoadFromFile vector file to be loaded into the internal memory from disk after the template, then the view data in the internal memory is presented on the basic display unit, in procedure for displaying, employed function is the OnPaint among the VC, promptly in this function, the dc (contextual devices) that calls the thumbnail window shows the template of loading.Load the back and preserve the address pDxfImg of view data in internal memory, use for later step.

Basic display unit shows as shown in Figure 3 after having loaded polar plot.

After having loaded polar plot, the user just can begin to adjust the position of polar plot in the orthopaedics image, mates by operations such as moving, rotate, amplification is dwindled.If the user thinks that the polar plot that loads can't well mate with the orthopaedics image all the time, can earlier former polar plot be deleted from display, reselect other polar plot then, till obtaining satisfied result.The employed function of delete procedure is the OnPaint among the VC, in this function, calls the dc (contextual devices) of thumbnail window, wipes former polar plot.In this example, we have defined class LAutoCadDXF, this class definition some interfaces, as SetPosition, SetRatio, finish operation, as shown in Figure 4 to vector graphics.

The pDxfImg vector image data pointer of preserving is previously passed to the LAutoCadDXF class, can operate the polar plot object of loading.

For can be on the orthopaedics image operating vector figure easily, also to respond window message, just can make polar plot along with corresponding adjustment is made in user's operation, the concrete operations flow process is as follows:

Response left mouse button is pressed message OnLButtonDown in polar plot, when left mouse button is pressed, this moment, polar plot entered mobile status (as shown in Figure 5), and (ShowMousePosX ShowMousePosY) uses for later step with polar plot viewing area oldRect to write down mouse coordinates.

Response mouse mobile messaging OnMouseMove in polar plot, the mouse coordinates parameter x, the y that pass over by OnMouseMove, with the mouse coordinates (ShowMousePosX that preserves in the step 1, ShowMousePosY) calculate displacement offsetX, the offsetY of mouse in directions X and Y direction, calculate the new viewing area newRect of polar plot by the oldRect in the step a) and the offsetX, the offsetY that calculate, the interface SetPosition that calls in the class LAutoCadDXF class is provided with the new display position of polar plot.

The OnRButtonDown message that response right mouse button is pressed in polar plot, this moment, polar plot entered rotation status (as shown in Figure 6), calculated angle of inclination npreAngle by mouse coordinates and polar plot centre coordinate, kept this value and used for later step.The concrete method of calculating angle is:

Mouse coordinates is oldPtX when supposing to press right mouse button, oldPtY, this coordinate is the coordinate in the polar plot coordinate system, be converted into screen coordinate so also will call ClientToScreen, obtain polar plot zone pRect by GetWindowRect after the coordinate conversion, obtain the coordinate centerX of polar plot center in screen thus, centerY, centerX=in the Windows system (pRect.right+pRect.left)/2, centerY=(pRect.top+pRect.bottom)/2 can obtain npreAngle=arctan ((oldPtY-centerY)/oldPtX-centerX)) by trigonometric function then

Response OnMouseMove message in the mouse moving process, calculate angle of inclination ncurAngle by mouse coordinates and polar plot centre coordinate, calculate the same step a) of method of angle, ncurAngle and the npreAngle that preserves previously subtracted each other obtain anglec of rotation angle, call the SetAngle clock-face diagram then, polar plot can move and rotates along with mouse, decontrols right mouse button, and polar plot stops the rotation.

Response mouse roller message OnMouseWheel in polar plot, rolling distance by mouse roller calculates amplification (dwindling) coefficient nRatio, concrete computing method genju passes to orthopaedics image and vector image simultaneously with coefficient, call the image after explicit function shows amplification (dwindling), as for polar plot, calling interface SetRatio gets final product, and this moment, polar plot dwindled (shown in Fig. 7,8) along with the orthopaedics image amplifies together.

Those skilled in the art do not break away from essence of the present invention and spirit, can there be the various deformation scheme to realize the present invention, the above only is the preferable feasible embodiment of the present invention, be not so limit to interest field of the present invention, the equivalent structure that all utilizations instructions of the present invention and accompanying drawing content are done changes, and all is contained within the interest field of the present invention.

Claims (1)

1. method that shows vector graphics on medical image is characterized in that may further comprise the steps:

A. set up unified coordinate system;

B. in coordinate system, load medical image;

C. load vector graphics in coordinate system, vector graphics is presented on the medical image;

D. medical image pixel f (x _i, y _i) (0＜i≤N) coordinate in coordinate system is (x _i, y _i), wherein N is the number of pixel; Vector graphics is at the image pixel f of coordinate system ₁(x1 _i, y1 _i) (coordinate of 0＜i≤N1) is (x1 _i, y1 _i), wherein N1 is the number of pixel; According to the coordinate of medical image pixel, obtain the size of the border rectangle in medical image zone, it is Far Left left, rightmost right, top top and the bottommost bottom of the image-region of reference that the border rectangle has with coordinate; Four attributes of this border rectangle are obtained by following formula respectively:

left＝min{x ₀，x ₁，...，x _N-1}；

right＝max{x ₀，x ₁，...，x _N-1}；

top＝max{y ₀，y ₁，...，y _N-1}；

bottom＝min{y ₀，y ₁，...，y _N-1}；

Wherein N is the number of the image pixel in the rectangle of border, the border rectangle of medical image be (left, right, top, bottom), the border rectangle of vector graphics is (left ₁, right ₁, top ₁, bottom ₁);

E. the border rectangle of medical image obtains the center of this rectangle and is:

$(\frac{(\mathrm{left}+\mathrm{right})}{2},\frac{\mathrm{top}+\mathrm{bottom}}{2})$

The center that the border rectangle of vector graphics obtains this rectangle is:

$(\frac{({\mathrm{left}}_1+{\mathrm{right}}_1)}{2},\frac{{\mathrm{top}}_1+{\mathrm{bottom}}_1}{2})$

After medical image was moved, the moving method step of vector graphics was as follows:

E1. the distance that move in coordinate system relative to former rectangular centre at the center of obtaining the new border rectangle of medical image for (x, y),

E2. when the distance that moves is not 0, the border rectangle mobile phase of vector graphics with distance,

E3. in new border rectangle, show new vector graphics;

F. during the medical image convergent-divergent, the center of medical image border rectangle does not change, and medical image convergent-divergent rate is α, and the distance of the central point of the border rectangle of vector graphics and the border rectangle of medical image is:

$s={({(\frac{({\mathrm{left}}_1+{\mathrm{right}}_1)}{2}-\frac{(\mathrm{left}+\mathrm{right})}{2})}^2+{(\frac{\mathrm{top}+\mathrm{bottom}}{2}-\frac{{\mathrm{top}}_1+{\mathrm{bottom}}_1}{2})}^2)}^{1/2}$

New distance is after the convergent-divergent:

s ₁＝α*s

The angle β at the center of former vector graphics and the center of medical image is:

$\mathrm{\&beta;}=\mathrm{arctg}\left(\right|\frac{\frac{\mathrm{top}+\mathrm{bottom}}{2}-\frac{{\mathrm{top}}_1+{\mathrm{bottom}}_1}{2}}{\frac{{\mathrm{left}}_1+{\mathrm{right}}_1}{2}-\frac{\mathrm{left}+\mathrm{right}}{2}}\left|\right),$

The center of the border rectangle of then new vector graphics is:

$(\frac{(\mathrm{left}+\mathrm{right})}{2}-s_1*\cos \mathrm{\&beta;},\frac{\mathrm{top}+\mathrm{bottom}}{2}+s_1*\sin \mathrm{\&beta;})$

Can obtain new border rectangle after former vector graphics border rectangle calculates in the same way, draw new vector graphics in new border rectangle class.

Images