CN105046728A - Visualized organ based virtual tomographic image reconstruction method - Google Patents

Abstract

The present invention relates to the field of medical tomographic image processing and aims to provide a visualized organ based virtual tomographic image reconstruction method. The visualized organ based virtual tomographic image reconstruction method comprises: adopting an inertial coordinate system of a visualized organ as a quantitative positioning reference; choosing a new coordinate system according to needs, and reconstructing the visualized organ in the new coordinate system; then sectioning the visualized organ along a coordinate axis of the new coordinate system to obtain a tomographic image so as to obtain a tomographic image with an arbitrary angle. The method has the beneficial effects that: with respect to the visualized organ, virtual tomography performed along any axis relative to an inertial reference system can be realized so as to expand application of the visualized tomographic image to medical education; and with respect to a living organ, unlimited virtual scanning performed along any axis relative to the inertial reference system can be realized by scanning once only, thereby improving the reliability of a tomographic image diagnosis result.

Description

A kind of virtual fault image method for reconstructing based on visual organ

Technical field

The present invention relates to the process field of medical science fault image, particularly a kind of virtual fault image method for reconstructing based on visual organ.

Background technology

National Library of Medicine comes into effect " Visual Human Body (VisibleHuman, HV) plan " in 1989.Within 1994, complete the first male sex's visible human in the world, various countries follow up one after another afterwards.Advocate as Zhong Shi town academician and opened up Chinese Digital Human and digital medical research, completing the structure of the first Chinese Digital Human in 2003.VH has caused the revolution of a medical education in the world.But there is a great defect at the fault image cut based on organ: the cross section perpendicular or parallel with cutting axis can only be observed, because organ can only be cut once.

Fault image field is so important, to such an extent as to has only successively just been awarded six Nobel Prizes in nuclear magnetic resonance field.In clinical medicine, fault image is still the gold criterion of diagnosis various diseases.In preclinical medicine, organ cross-section analysis is very general, and fault image, to description organ morphosis and the change of prediction organ, has the effect do not replaced.But the same with visible human, scanning device also can only observe the organ fault image perpendicular or parallel with scan axis.

Although the functions such as (Reposition), again section (Reslice) reorientated by the medical science fault image post-processed softwares such as Mimics, AnalyzeDirect, but these functions are not the fault image reconstructions along appointment axle of relatively unique, constant inertial coordinates system, because we cannot allow experimenter accomplish, twice detection is all in duplicate attitude.Different scanning attitudes causes scan axis to be different.Like this, just as not having " first meridian ", " longitude " is just nonsensical the same, not fixing, unique inertial coordinates system organ fault image does not just have quantitative description, because different coordinate systems, coordinate axis is just different, and square with the axis or parallel fault image is also just different.

Based on the fault image that cutting or the scanning of organ are formed, storehouse forms visual organ.The relatively unique inertial coordinates system of the information realization of visual organ how is utilized to be a gordian technique urgently broken through along the virtual cross sectional reconstruction of arbitrary coordinate axle.Setting up a kind of relative inertness coordinate system and carry out the method that coordinate transform realizes carrying out along any direction virtual cross sectional reconstruction, is all have important theory and practice meaning to fields such as medical education, clinical diagnose or basic medical researches.

Summary of the invention

The object of the invention is to overcome above-mentioned defect, provide a kind of inertia coordinate system unique based on visual organ visual organ to be carried out to the virtual fault image method for reconstructing based on visual organ of Virtual cropping and scanning along arbitrary axis.

In order to solve the problems of the technologies described above, the technical solution used in the present invention is:

Based on a virtual fault image method for reconstructing for visual organ, comprise the following steps:

Step 1, with the principal axis of inertia of visual organ for coordinate axis builds inertial coordinates system;

Step 2, build new coordinate system, visual organ is got a some i, and acquisition point i is at the position coordinates of inertial coordinates system obtain the centroid position coordinate (x of visual organ relative inertness coordinate system _c, y _c, z _c), when acquisition is rebuild, visual organ is around x, y, the angular displacement alpha that z-axis is rotated, and beta, gamma, substitutes into formula by above-mentioned parameter:

$\left(\begin{array}{c}{x}_i^{new}+x_c\\ y_i^{new}+y_c\\ z_i^{new}+z_c\end{array}\right)=\left(\begin{array}{ccc}1& 0& 0\\ 0& \cos \left(\mathrm{\α}\right)& -\sin \left(\mathrm{\α}\right)\\ 0& \sin \left(\mathrm{\α}\right)& \cos \left(\mathrm{\α}\right)\end{array}\right)\left(\begin{array}{ccc}\cos \left(\mathrm{\β}\right)& 0& \sin \left(\mathrm{\β}\right)\\ 0& 1& 0\\ -\sin \left(\mathrm{\β}\right)& 0& \cos \left(\mathrm{\β}\right)\end{array}\right)\left(\begin{array}{ccc}\cos \left(\mathrm{\γ}\right)& -\sin \left(\mathrm{\γ}\right)& 0\\ \sin \left(\mathrm{\γ}\right)& \cos \left(\mathrm{\γ}\right)& 0\\ 0& 0& 1\end{array}\right)\left(\begin{array}{c}{x}_i^o-x_c\\ y_i^o-y_c\\ z_i^o-z_c\end{array}\right)$

Obtain visual organ being put the position coordinates of i in new coordinate system each point on visual organ is calculated described visual organ is rebuild in new coordinate system;

Step 3, by what obtain in step 2 substitute into formula:

$\left\{\begin{array}{c}{x}_i^{\′}=Round\left(\left(\frac{\left(x_i^{new}-x_c\right)}{R},0\right)R\right)\\ y_i^{\′}=Round\left(\left(\frac{\left(y_i^{new}-y_c\right)}{R},0\right)R\right)\\ z_i^{\′}=Round\left(\left(\frac{\left(z_i^{new}-z_c\right)}{R},0\right)R\right)\end{array}\right.$

Image procossing is carried out to the visual organ in new coordinate system, obtains as the position coordinates (x' of the some i on visual organ in new coordinate system after isotropy reconstruction _i, y' _i, z' _i), wherein R is isotropy parameter;

Step 4, along the coordinate axis of new coordinate system, virtual cross sectional reconstruction is carried out to visual organ, comprising:

By (the x' obtained in step 3 _i, y' _i, z' _i) substitute into formula:

Obtain the image of visual organ virtual tomography along the x-axis direction, wherein, p is the some points in coordinate axis, $p\&Element;(0,\&PlusMinus;R,\&PlusMinus;2R......),0<\mathrm{\&Delta;}<\frac{1}{2}R;$

By (the x' obtained in step 3 _i, y' _i, z' _i) substitute into formula:

Obtain the image of visual organ virtual tomography along the y-axis direction, wherein, p is the some points in coordinate axis, $p\&Element;(0,\&PlusMinus;R,\&PlusMinus;2R......),0<\mathrm{\&Delta;}<\frac{1}{2}R;$

By (the x' obtained in step 3 _i, y' _i, z' _i) substitute into formula:

Obtain the image of visual organ virtual tomography along the z-axis direction, wherein, p is the some points in coordinate axis, $p\&Element;(0,\&PlusMinus;R,\&PlusMinus;2R......),0<\mathrm{\&Delta;}<\frac{1}{2}R.$

Beneficial effect of the present invention is: the present invention creatively proposes to carry out virtual fault image method for reconstructing based on unique principal axis of inertia of visual organ, compared to prior art, the present invention can realize Relative inertia reference system for visual organ and carry out virtual cross sectional reconstruction along arbitrary axis, thus expands the application of visual fault image in medical education; As long as run-down can be realized for organs of living beings can carry out unlimited virtual scan reconstruction along arbitrary axis by relative inertness coordinate system, based on uniqueness and the unchangeability of the inertial reference system of principal axis of inertia structure, realize the quantification of visual organ Comparison of Morphological, basic medical research, clinical diagnose are all had great significance.Unique, constant visual organ principal axis of inertia is once by the inertial coordinates system as visual organ, and morphological analysis is just become quantitatively from sxemiquantitative.Set up the prerequisite that coordinate system is fault image quantitative test, there is no coordinate system, change just nonsensical with fault image to description organ morphosis and prediction organ.And visual organ unique, constant inertial coordinate is that the analysis of organ fault image provides quantitative location basis of reference.

Accompanying drawing explanation

Fig. 1 is the process flow diagram of the embodiment of the present invention based on the virtual fault image method for reconstructing of visual organ;

Fig. 2 is the stereographic map of the embodiment of the present invention one Cerebrum;

Fig. 3 is the front view of the embodiment of the present invention one Cerebrum;

Fig. 4 is the side view of the embodiment of the present invention one Cerebrum;

Fig. 5 is the vertical view of the embodiment of the present invention one Cerebrum;

Fig. 6 is the stereographic map that the embodiment of the present invention one is positioned at inertial coordinates system deutocerebrum organ;

Fig. 7 is the front view that the embodiment of the present invention one is positioned at inertial coordinates system deutocerebrum organ;

Fig. 8 is the side view that the embodiment of the present invention one is positioned at inertial coordinates system deutocerebrum organ;

Fig. 9 is the vertical view that the embodiment of the present invention one is positioned at inertial coordinates system deutocerebrum organ;

Figure 10 is the stereographic map of visual 1st metatarsal of the embodiment of the present invention 2 first;

Figure 11 is the positive apparent direction sectional view of visual 1st metatarsal of the embodiment of the present invention 2 first;

Figure 12 is the side-looking direction sectional view of visual 1st metatarsal of the embodiment of the present invention 2 first;

Figure 13 be visual 1st metatarsal of the embodiment of the present invention 2 first overlook direction sectional view;

Figure 14 is the stereographic map of visual 1st metatarsal of the embodiment of the present invention 2 second;

Figure 15 is the positive apparent direction sectional view of visual 1st metatarsal of the embodiment of the present invention 2 second;

Figure 16 is the side-looking direction sectional view of visual 1st metatarsal of the embodiment of the present invention 2 second;

Figure 17 be visual 1st metatarsal of the embodiment of the present invention 2 second overlook direction sectional view;

Figure 18 is the stereographic map of visual 1st metatarsal of the embodiment of the present invention 2 the 3rd;

Figure 19 is the positive apparent direction sectional view of visual 1st metatarsal of the embodiment of the present invention 2 the 3rd;

Figure 20 is the side-looking direction sectional view of visual 1st metatarsal of the embodiment of the present invention 2 the 3rd;

Figure 21 be visual 1st metatarsal of the embodiment of the present invention 2 the 3rd overlook direction sectional view;

Figure 22 is the stereographic map of visual 1st metatarsal of the embodiment of the present invention 2 the 4th;

Figure 23 is the positive apparent direction sectional view of visual 1st metatarsal of the embodiment of the present invention 2 the 4th;

Figure 24 is the side-looking direction sectional view of visual 1st metatarsal of the embodiment of the present invention 2 the 4th;

Figure 25 be visual 1st metatarsal of the embodiment of the present invention 2 the 4th overlook direction sectional view.

Embodiment

By describing technology contents of the present invention, structural attitude in detail, realized object and effect, accompanying drawing is coordinated to be explained in detail below in conjunction with embodiment.

The design of most critical of the present invention is: adopt the inertial coordinates system of visual organ as location basis of reference, new coordinate system can be selected according to demand, and visual organ is rebuild in new coordinate system, visual organ section is obtained fault image by the coordinate axis again along new coordinate system, thus fault image at any angle can be obtained, as long as run-down can be realized for organs of living beings can carry out unlimited virtual scan along arbitrary axis by relative inertness coordinate system.

Refer to shown in Fig. 1, the virtual fault image method for reconstructing based on visual organ of the present embodiment, comprises the following steps:

Step 1, with the principal axis of inertia of visual organ for coordinate axis builds inertial coordinates system;

Step 2, build new coordinate system, visual organ is got a some i, and acquisition point i is at the position coordinates of inertial coordinates system obtain the centroid position coordinate (x of visual organ relative inertness coordinate system _c, y _c, z _c), when acquisition is rebuild, visual organ is around x, y, the angular displacement alpha that z-axis is rotated, and beta, gamma, substitutes into formula by above-mentioned parameter:

$\left(\begin{array}{c}{x}_i^{new}+x_c\\ y_i^{new}+y_c\\ z_i^{new}+z_c\end{array}\right)=\left(\begin{array}{ccc}1& 0& 0\\ 0& \cos \left(\mathrm{\&alpha;}\right)& -\sin \left(\mathrm{\&alpha;}\right)\\ 0& \sin \left(\mathrm{\&alpha;}\right)& \cos \left(\mathrm{\&alpha;}\right)\end{array}\right)\left(\begin{array}{ccc}\cos \left(\mathrm{\&beta;}\right)& 0& \sin \left(\mathrm{\&beta;}\right)\\ 0& 1& 0\\ -\sin \left(\mathrm{\&beta;}\right)& 0& \cos \left(\mathrm{\&beta;}\right)\end{array}\right)\left(\begin{array}{ccc}\cos \left(\mathrm{\&gamma;}\right)& -\sin \left(\mathrm{\&gamma;}\right)& 0\\ \sin \left(\mathrm{\&gamma;}\right)& \cos \left(\mathrm{\&gamma;}\right)& 0\\ 0& 0& 1\end{array}\right)\left(\begin{array}{c}{x}_i^o-x_c\\ y_i^o-y_c\\ z_i^o-z_c\end{array}\right)$

Obtain visual organ being put the position coordinates of i in new coordinate system each point on visual organ is calculated described visual organ is rebuild in new coordinate system;

First can fasten at inertial coordinate and build arbitrary new coordinate system, fasten at inertial coordinate and choose arbitrfary point i, adopt above-mentioned computing method just can obtain the coordinate of i in new coordinate system, by each point that computer calculate inertial coordinate is fastened, thus visual organ is rebuild in new coordinate system.

Step 3, by what obtain in step 2 substitute into formula:

$\left\{\begin{array}{c}{x}_i^{\&prime;}=Round\left(\left(\frac{\left(x_i^{new}-x_c\right)}{R},0\right)R\right)\\ y_i^{\&prime;}=Round\left(\left(\frac{\left(y_i^{new}-y_c\right)}{R},0\right)R\right)\\ z_i^{\&prime;}=Round\left(\left(\frac{\left(z_i^{new}-z_c\right)}{R},0\right)R\right)\end{array}\right.$

Image procossing is carried out to the visual organ in new coordinate system, obtains as the position coordinates (x' of the some i on visual organ in new coordinate system after isotropy reconstruction _i, y' _i, z' _i), wherein R is isotropy parameter;

Round is bracket function, each point in new coordinate system can be arranged, will put corresponding arrangement in corresponding layer, and improve image effect after above-mentioned calculating.

Step 4, along the coordinate axis of new coordinate system, virtual tomography is carried out to visual organ, comprising:

By (the x' obtained in step 3 _i, y' _i, z' _i) substitute into formula:

Obtain the image of visual organ virtual tomography along the x-axis direction, wherein, p is the some points in coordinate axis, $p\&Element;(0,\&PlusMinus;R,\&PlusMinus;2R......),0<\mathrm{\&Delta;}<\frac{1}{2}R;$

By (the x' obtained in step 3 _i, y' _i, z' _i) substitute into formula:

Obtain the image of visual organ virtual tomography along the y-axis direction, wherein, p is the some points in coordinate axis, $p\&Element;(0,\&PlusMinus;R,\&PlusMinus;2R......),0<\mathrm{\&Delta;}<\frac{1}{2}R;$

By (the x' obtained in step 3 _i, y' _i, z' _i) substitute into formula:

Obtain the image of visual organ virtual tomography along the z-axis direction, wherein, p is the some points in coordinate axis, $p\&Element;(0,\&PlusMinus;R,\&PlusMinus;2R......),0<\mathrm{\&Delta;}<\frac{1}{2}R.$

After above-mentioned calculating, just can draw final fault image, make the visual organ obtained after a section can carry out arbitrarily time virtual tomography, make visual organ can efficiency utilization, expand the range of application of visual fault image in medical education.

From foregoing description, beneficial effect of the present invention is: the present invention creatively proposes to carry out virtual fault image method for reconstructing based on unique principal axis of inertia of visual organ, compared to prior art, the present invention can realize Relative inertia reference system for visual organ and carry out virtual cross sectional reconstruction along arbitrary axis, thus expands the application of visual fault image in medical education; As long as run-down can be realized for organs of living beings can carry out unlimited virtual scan along arbitrary axis by relative inertness coordinate system, based on uniqueness and the unchangeability of the inertial reference system of principal axis of inertia structure, realize the quantification of visual organ Comparison of Morphological, basic medical research, clinical diagnose are all had great significance.The principal axis of inertia of unique, constant visual organ is once by the inertial coordinates system as visual organ, and morphological analysis is just developed into quantitatively by sxemiquantitative.Set up the prerequisite that coordinate system is fault image quantitative test, there is no coordinate system, change just nonsensical with fault image to description organ morphosis and prediction organ.And visual organ unique, constant inertial coordinate is that the analysis of organ fault image provides quantitative location basis of reference.

Further, before described " obtaining visual organ ", comprise further:

Visual organ is formed by the section broken layer image storehouse of organ.

Seen from the above description, organs of living beings is successively scanned, after arrangement storehouse, just form visual organ.

Please refer to shown in Fig. 2 to Fig. 9, embodiments of the invention one are:

(1) visualizing cortico organ:

Based on the fault image that Cerebrum scanning is formed, storehouse forms visualizing cortico model;

(2) inertial coordinates system: form a coordinate system by the principal axis of inertia of Cerebrum, using the inertial coordinates system of this coordinate system as Cerebrum.

Please refer to shown in Figure 10 to Figure 25, embodiments of the invention two are:

(1) visual 1st metatarsal (BigBrain):

Carry out twice sweep to the 1st metatarsal (BigBrain) respectively with two different angles, the fault image of formation respectively storehouse forms the first visual 1st metatarsal (BigBrain) and the second visual 1st metatarsal (BigBrain);

(2) inertial coordinates system: form a coordinate system by the principal axis of inertia of the 1st metatarsal (BigBrain), using the inertial coordinates system of this coordinate system as the 1st metatarsal (BigBrain), the 1st metatarsal (BigBrain) in inertial coordinates system is defined as the 3rd visual 1st metatarsal (BigBrain);

(3) Virtual cropping: according to the needs such as research, diagnosis, analysis of the 1st metatarsal (BigBrain) morphosis, the inertial coordinates system of relative 1st metatarsal (BigBrain) carries out coordinate conversion, forms new coordinate system.In new coordinate system, rebuild along coordinate axis, form the virtual scan fault image of the 1st metatarsal (BigBrain).

(4) coordinate transform carried out the 1st metatarsal (BigBrain) is as follows:

$\left(\begin{array}{c}{x}_i^{new}+x_c\\ y_i^{new}+y_c\\ z_i^{new}+z_c\end{array}\right)=\left(\begin{array}{ccc}1& 0& 0\\ 0& \cos \left(\mathrm{\&alpha;}\right)& -\sin \left(\mathrm{\&alpha;}\right)\\ 0& \sin \left(\mathrm{\&alpha;}\right)& \cos \left(\mathrm{\&alpha;}\right)\end{array}\right)\left(\begin{array}{ccc}\cos \left(\mathrm{\&beta;}\right)& 0& \sin \left(\mathrm{\&beta;}\right)\\ 0& 1& 0\\ -\sin \left(\mathrm{\&beta;}\right)& 0& \cos \left(\mathrm{\&beta;}\right)\end{array}\right)\left(\begin{array}{ccc}\cos \left(\mathrm{\&gamma;}\right)& -\sin \left(\mathrm{\&gamma;}\right)& 0\\ \sin \left(\mathrm{\&gamma;}\right)& \cos \left(\mathrm{\&gamma;}\right)& 0\\ 0& 0& 1\end{array}\right)\left(\begin{array}{c}{x}_i^o-x_c\\ y_i^o-y_c\\ z_i^o-z_c\end{array}\right)$

In formula, be the position coordinates of the 1st metatarsal (BigBrain) relative inertness coordinate system, represent the 1st metatarsal (BigBrain) position coordinates in new coordinate system, (x _c, y _c, z _c) representing the centroid position coordinate of the 1st metatarsal (BigBrain) relative inertness coordinate system, α, beta, gamma represents respectively around x, y, the angular displacement that z-axis is rotated, and forms the 4th visual 1st metatarsal (BigBrain) after coordinate transform.

(5) in new coordinate system, along carrying out the process of fault image isotropy to the 4th visual 1st metatarsal (BigBrain):

$\left\{\begin{array}{c}{x}_i^{\&prime;}=Round\left(\left(\frac{\left(x_i^{new}-x_c\right)}{R},0\right)R\right)\\ y_i^{\&prime;}=Round\left(\left(\frac{\left(y_i^{new}-y_c\right)}{R},0\right)R\right)\\ z_i^{\&prime;}=Round\left(\left(\frac{\left(z_i^{new}-z_c\right)}{R},0\right)R\right)\end{array}\right.$

(x' in formula _i, y' _i, z' _i) be some i on the 1st metatarsal (BigBrain) in a coordinate system isotropy rebuild after position coordinates, R is isotropy parameter.

(6) along the axis of new coordinate system, virtual tomography is carried out to the 4th visual 1st metatarsal (BigBrain):

Along transverse axis:

Along Z-axis:

Along scan axis:

In formula, p is the some points in coordinate axis,

By the present embodiment, be no matter the cutting to organ or scanning, after they are become fault image, through the configuration of surface of all reducible organ of storehouse and the visual organ of inner structure.Once through overscanning, square with the axis or parallel cross section can only be observed cutting on scan axis.But the present invention is based on visual organ unique, constant inertial coordinates system can observe the arbitrary section of Relative inertia reference system along section axle quantitatively.

As Socrates said " people can not step into same river for twice ", we do not accomplish to make experimenter be scanned twice with identical attitude yet, organ can not be cut twice in tomography mode.This method is using unique for visual organ, the changeless principal axis of inertia inertial coordinates system as visual organ, according to the unchangeability of visual organ to coordinate transform, Virtual Reconstruction fault image, can carry out not by the virtual tomography of number of times restriction, because the virtual fault image of organ and the section attitude of organ have nothing to do, realize the quantitative test of visual organ fault image thus, this is innovation of the present invention just.Just as Ke Kendelu (Kendrew), Lu Ge (Klug) scan protein from different perspectives by x-ray, obtaining the three-dimensional structure of protein by diffraction analysis, the same (the former 1962 obtains the Nobel Prizes, the latter's nineteen eighty-two obtains the Nobel Prize), the present invention is based on that visual organ is unique, fixing inertial coordinates system, achieve by coordinate transform the quantification that fault image describes organ morphosis, the change of prediction organ.

In sum, unique principal axis of inertia based on visual organ provided by the invention carries out virtual fault image method for reconstructing, compared to prior art, the present invention can realize Relative inertia reference system for visual organ and carry out virtual cross sectional reconstruction along arbitrary axis, thus expands the application of visual fault image in medical education; As long as run-down can be realized for organs of living beings can carry out unlimited virtual scan along arbitrary axis by relative inertness coordinate system, based on uniqueness and the unchangeability of the inertial reference system of principal axis of inertia structure, realize the quantification of visual organ Comparison of Morphological, basic medical research, clinical diagnose are all had great significance.The principal axis of inertia of unique, constant visual organ is once by the inertial coordinates system as visual organ, and morphological analysis is just developed into quantitatively by sxemiquantitative.Set up the prerequisite that coordinate system is fault image quantitative test, there is no coordinate system, change just nonsensical with fault image to description organ morphosis and prediction organ.And visual organ unique, constant inertial coordinate is that the analysis of organ fault image provides quantitative location basis of reference.

The foregoing is only embodiments of the invention; not thereby the scope of the claims of the present invention is limited; every utilize instructions of the present invention and accompanying drawing content to do equivalent structure or equivalent flow process conversion; or be directly or indirectly used in other relevant technical fields, be all in like manner included in scope of patent protection of the present invention.

Claims (2)

1., based on a virtual fault image method for reconstructing for visual organ, it is characterized in that, comprising:

Step 1, with the principal axis of inertia of visual organ for coordinate axis builds inertial coordinates system;

Step 2, build new coordinate system, visual organ is got a some i, and acquisition point i is at the position coordinates of inertial coordinates system obtain the centroid position coordinate (x of visual organ relative inertness coordinate system _c, y _c, z _c), when acquisition is rebuild, visual organ is around x, y, the angular displacement alpha that z-axis is rotated, and beta, gamma, substitutes into formula by above-mentioned parameter:

$\left(\begin{array}{c}{x}_i^{new}+x_c\\ y_i^{new}+y_c\\ z_i^{new}+z_c\end{array}\right)=\left(\begin{array}{ccc}1& 0& 0\\ 0& \cos \left(\mathrm{\&alpha;}\right)& -\sin \left(\mathrm{\&alpha;}\right)\\ 0& \sin \left(\mathrm{\&alpha;}\right)& \cos \left(\mathrm{\&alpha;}\right)\end{array}\right)\left(\begin{array}{ccc}\cos \left(\mathrm{\&beta;}\right)& 0& \sin \left(\mathrm{\&beta;}\right)\\ 0& 1& 0\\ -\sin \left(\mathrm{\&beta;}\right)& 0& \cos \left(\mathrm{\&beta;}\right)\end{array}\right)\left(\begin{array}{ccc}\cos \left(\mathrm{\&gamma;}\right)& -\sin \left(\mathrm{\&gamma;}\right)& 1\\ \sin \left(\mathrm{\&gamma;}\right)& \cos \left(\mathrm{\&gamma;}\right)& 0\\ 0& 0& 1\end{array}\right)\left(\begin{array}{c}{x}_i^o-x_c\\ y_i^o-y_c\\ z_i^o-z_c\end{array}\right)$

Obtain visual organ being put the position coordinates of i in new coordinate system after each point on organ visual in inertial coordinates system is calculated, described visual organ can be rebuild in new coordinate system;

Step 3, by what obtain in step 2 substitute into formula:

$\left\{\begin{array}{c}{x}_i^{\&prime;}=Round\left(\left(\frac{\left(x_i^{new}-x_c\right)}{R},0\right)R\right)\\ y_i^{\&prime;}=Round\left(\left(\frac{\left(y_i^{new}-y_c\right)}{R},0\right)R\right)\\ z_i^{\&prime;}=Round\left(\left(\frac{\left(z_i^{new}-z_c\right)}{R},0\right)R\right)\end{array}\right.$

Image procossing is carried out to the visual organ in new coordinate system, obtain into the some i on visual organ in new coordinate system isotropy rebuild after position coordinates (x ' _i, y ' _i, z ' _i), wherein R is isotropy parameter;

Step 4, along the coordinate axis of new coordinate system, virtual tomography is carried out to visual organ, comprising:

By obtain in step 3 (x ' _i, y ' _i, z ' _i) substitute into formula:

Obtain the image of visual organ virtual tomography along the x-axis direction, wherein, p is the some points in coordinate axis, p ∈ (0, ± R, ± 2R ...),

By obtain in step 3 (x ' _i, y ' _i, z ' _i) substitute into formula:

Obtain the image of visual organ virtual tomography along the y-axis direction, wherein, p is the some points in coordinate axis, p ∈ (0, ± R, ± 2R ...),

By obtain in step 3 (x ' _i, y ' _i, z ' _i) substitute into formula:

Obtain the image of visual organ virtual tomography along the z-axis direction, wherein, p is the some points in coordinate axis, p ∈ (0, ± R, ± 2R ...),

2. the virtual fault image method for reconstructing based on visual organ according to claim 1, is characterized in that, before described " obtaining visual organ ", comprises further:

Cut or scan the section broken layer image storehouse formed by organ and form visual organ.