CN102512140B - Method for locating optical projection tomography imaging rotation center - Google Patents

Method for locating optical projection tomography imaging rotation center Download PDF

Info

Publication number
CN102512140B
CN102512140B CN 201110386219 CN201110386219A CN102512140B CN 102512140 B CN102512140 B CN 102512140B CN 201110386219 CN201110386219 CN 201110386219 CN 201110386219 A CN201110386219 A CN 201110386219A CN 102512140 B CN102512140 B CN 102512140B
Authority
CN
China
Prior art keywords
center
rotation
sinogram
cloth
width
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.)
Active
Application number
CN 201110386219
Other languages
Chinese (zh)
Other versions
CN102512140A (en
Inventor
田捷
董迪
秦承虎
杨鑫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Automation of Chinese Academy of Science
Original Assignee
Institute of Automation of Chinese Academy of Science
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Automation of Chinese Academy of Science filed Critical Institute of Automation of Chinese Academy of Science
Priority to CN 201110386219 priority Critical patent/CN102512140B/en
Publication of CN102512140A publication Critical patent/CN102512140A/en
Application granted granted Critical
Publication of CN102512140B publication Critical patent/CN102512140B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

The invention discloses a locating optical projection coherence tomography rotation center method, which introduces a rotation center computing method based on mass center track in the field of X-ray computed tomography imaging into the field of optical projection coherence tomography imaging and simultaneously is combined with a rotation center fine location method of rebuilding image mean-square deviation. The locating optical projection coherence tomography rotation center method can achieve automatic, quick and accurate rotation center location according to different experimental conditions.

Description

A kind of method of locating optical projection fault imaging center of rotation
Technical field
The present invention relates to optical projection fault imaging (Optical Projection Tomography is called for short OPT) technology, specially refer to a kind of method of locating optical projection fault imaging center of rotation.
Background technology
The optical projection tomography technology is to utilize light characteristics along straightline propagation in the small size organism, and the visible emitting line penetrates sample, gathers the sample projection view of a plurality of angles then with camera, carries out three-dimensional imaging.Specifically, when carrying out the optical projection fault imaging, need carry out multi-angle scanning to sample, generally adopt automatically controlled turntable that sample is carried out the step-by-step movement rotation, whenever rotate to an angle acquisition one or more projected image, when sample rotates a circle, stop scanning.The data that the optical projection computed tomography (SPECT) system finally collects are two-dimensional projection image that light passes sample under a series of different angles, if a certain provisional capital of all projected images is extracted, be superposed to piece image by row successively according to scanning sequency, just can obtain a similar sinusoidal sinogram, each width of cloth sinogram correspondence level of sample rebuild the cross section, all sinograms are the just corresponding three-dimension disclocation rebuilding body of sample, the process from sinogram to sample three-dimension disclocation structure is called optical projection fault imaging three-dimensional reconstruction.
The optical projection tomography technology can be realized the structure and the molecular specificity functional imaging of 1-10mm yardstick biological sample, have resolution height, structure-function integration, radiationless, plurality of advantages such as cost is low, it can carry out the qualitative and quantitative study of cellular level to living body biological in small scale, realize organism in real time, noinvasive, dynamically, imaging in vivo.But the optical projection tomography technology is needing the position of pre-aligned turntable rotating shaft on detector before the three-dimensional imaging, promptly high-quality three-dimensional imaging can be carried out in the center of rotation position.
The center of rotation computational methods of main flow have four kinds at present: first method is that the sample experiment is before by the artificial turntable of adjusting, the turntable rotary middle spindle is alignd with the detector centrage as far as possible, its shortcoming is to need artificial the participation, and accuracy is difficult to guarantee that time is long; Second method is to scan known imitative body in advance before the sample experiment, utilizes the property calculation center of rotation of imitative body, the scanning of the sample that experimentizes again, and its shortcoming is a complex operation, time is long; The third method is the highlighted curve calculation center of rotation that utilizes on the sinogram that sample scanning obtains, and its shortcoming is that these highlighted curves are not all to exist on each sample, and the method does not have versatility; The 4th kind of method is the sinogram that scanning samples is obtained, an artificial given initial center of rotation, repeatedly rebuild in its vicinity then, hand picking or procedure identification go out to rebuild effect position preferably, its shortcoming is to need artificial given initial value, operation inconvenience, if initial value is far away excessively apart from actual value simultaneously, maximum error can appear in the method.In sum, existing center of rotation computational methods need artificial the participation, complicated operation, and precision and speed all need be improved.
Summary of the invention
(1) technical problem that will solve
For solving above-mentioned one or more problems, the invention provides a kind of localization method of optical projection fault imaging center of rotation, to realize the automatic location of center of rotation, improve localized speed and precision.
(2) technical scheme
A kind of localization method of optical projection fault imaging center of rotation is provided according to an aspect of the present invention.This localization method comprises: each width of cloth in several selected sinograms of optical projection fault imaging, utilize the centroid trajectory method to obtain the center of rotation of this width of cloth sinogram; The center of rotation of several selected sinograms is averaged, obtain the coarse positioning value of optical projection fault imaging center of rotation; Center of rotation with the coarse positioning value is an initial value, center of rotation fine search interval is set, in this fine search interval each width of cloth in the selected sinogram is repeatedly rebuild, in a series of reconstructed images, center of rotation is decided again as this width of cloth sinogram in the center of rotation position of choosing reconstructed image mean square deviation maximum, and the center of rotation of deciding again of all selected sinograms is averaged the meticulous position rotating of the optical projection fault imaging center that obtains.
(3) beneficial effect
The localization method of optical projection fault imaging center of rotation of the present invention has fully ensured localized rapidity, automatic type, accuracy, specifically:
(1) on sinogram is selected, fully take into account the influence of noise of optical projection fault imaging scan-data, only select to carry out on the strong sinogram of sample signal the center of rotation search, can guarantee center of rotation search institute based on data have than high s/n ratio, improved search precision of the present invention from the data selection;
(2) in the coarse positioning process, directly utilize centroid trajectory to carry out automatic center of rotation and estimate that the barycenter method of coarse positioning process only needs data for projection, does not need human intervention, has fully guaranteed automaticity of the present invention;
(3) in thin position fixing process, result based on coarse positioning searches for once more, and employing parallel processing technique, meticulously search near being positioned at the coarse positioning value, reduce mistake of the present invention and located probability, ensured final positioning accuracy, also reduced simultaneously the hunting zone, the integrating parallel technical guarantee rapidity of the present invention.
Description of drawings
Fig. 1 is high sample signal intensity sinogram screening process figure in the localization method of embodiment of the invention optical projection fault imaging center of rotation;
Fig. 2 is based on the center of rotation coarse positioning flow process of centroid trajectory in the localization method of embodiment of the invention optical projection fault imaging center of rotation;
Fig. 3 is based on the meticulous positioning flow of the center of rotation of reconstructed image mean square deviation in the localization method of embodiment of the invention optical projection fault imaging center of rotation;
Fig. 4 is at the center of rotation computational process sketch map of fruit bat pupa experiment in the localization method of embodiment of the invention optical projection fault imaging center of rotation; Wherein, Fig. 4 a has shown center of rotation and the average coarse positioning center of rotation that each selected sinogram calculates with the centroid trajectory method; Fig. 4 b has shown with the coarse positioning center of rotation and has rebuild the cross section that has obtained in the three-dimensional reconstruction body; Fig. 4 c has shown near the coarse positioning center with what the reconstructed image average variance method calculated and has decided center of rotation and meticulous position rotating center again; Fig. 4 d has shown and utilizes meticulous position rotating center to rebuild the three-dimensional reconstruction body section that obtains;
Fig. 5 is at the reconstructed results of fruit bat pupa experiment in the localization method of embodiment of the invention optical projection fault imaging center of rotation.
The specific embodiment
For making the purpose, technical solutions and advantages of the present invention clearer, below in conjunction with specific embodiment, and with reference to accompanying drawing, the present invention is described in more detail.Though this paper can provide the demonstration of the parameter that comprises particular value, should be appreciated that parameter need not definitely to equal corresponding value, but can in acceptable error margin or design constraint, be similar to described value.
The present invention is by introducing optical projection fault imaging field with the center of rotation computational methods based on centroid trajectory in the X ray computer fault imaging field, combination is based on the meticulous localization method of the center of rotation of reconstructed image mean square deviation simultaneously, design a kind of localization method of optical projection fault imaging center of rotation, can be at different experiment situations, realization is automatically, center of rotation is located fast, accurately.
Implementation of the present invention is divided into three key steps: high sample signal intensity sinogram screening, based on the center of rotation coarse positioning of centroid trajectory, based on the meticulous location of the center of rotation of reconstructed image mean square deviation, wherein coarse positioning and thin location all are to carry out on the sinogram that the first step is selected, to improve positioning accuracy.Utilize the experiment of fruit bat pupa to describe step of the present invention below, adopt the live body pupa to carry out the optical projection fault imaging in the experiment, the projection view pixel count of imaging system collection is 500*500, each pixel size is 16 microns, the sample angle number of evenly gathering that rotates a circle is 330, with the experiment of live body pupa is example, and detailed step of the present invention is as follows:
Step S1: for improving the center of rotation positioning accuracy, the present invention screens the sinogram that the optical projection fault imaging obtains in advance, extracts the sinogram that noise is little, signal intensity is good and carries out subsequent calculations.
Whether have higher signal to noise ratio in order to weigh sinogram, the present invention utilizes sinogram sample number of pixels to weigh its signal to noise ratio relatively, and the sinogram that the sample number of pixels is many more, its signal intensity are good more, and noise is low more, and it is good more to be rotated the center search effect.If but calculate the sample number of pixels of all sinograms, calculation of complex, the present invention handles by simplifying, only add up two row sample pixel averages of mutual upright projection angle in each sinogram, the sample signal intensity of relative quantification sinogram, under the situation of less error, significantly reduced the processing time.Certainly, all sinograms are carried out subsequent calculations also be fine, the software and hardware resources that only expends is bigger.
Fig. 1 is high sample signal intensity sinogram screening process figure in the localization method of embodiment of the invention optical projection fault imaging center of rotation.Below in conjunction with Fig. 1, screening process is elaborated.
Step S1-1: shown in Fig. 1 (a row), the present invention at first filters out the perpendicular projection view of two width of cloth projection ray directions, the present invention selects 0 ° and 90 ° of projection views, the colleague mutually of these two projection views can be mapped to two perpendicular row of projection angle in the corresponding sinogram, just can assess the sample signal intensity of all sinograms with these two projections; Those skilled in the art should understand, select 0 ° and 90 ° of schemes that projection views are optimum herein, in fact, also can select three width of cloth or several the projection sample signal intensity of coming estimation of sinusoidal figure more, and projection angle also can be selected arbitrarily, certainly, preferably select representative several projections.
Step S1-2: shown in Fig. 1 (b row), carry out threshold value and cut apart selecting two projection views, distinguish sample and background pixel, and then assess the sample signal intensity of each sinogram, in order to simplify calculating, it is the pixel grey scale average of corresponding projection view that segmentation threshold is set, and can distinguish sample and background so simple and effectively;
Step S1-3: the number of each row sample pixel of statistics on two width of cloth projections, it represents the relative signal intensity of this row, with the sample number of pixels is transverse axis, capable number is the longitudinal axis, can obtain each the row sample signal intensity curve of projection view shown in Fig. 1 (c row), the height of curve is represented the power of this row sample signal;
Step S1-4: the intensity curves of two width of cloth projections is averaged, obtains the average signal strength curve shown in Fig. 1 (d row), the statement by before as can be known, this curve correspondence the sample signal intensity of each sinogram;
Step S1-5: shown in Fig. 1 (d row), filter out the sinogram of signal intensity pro-25%, gray area among the figure is the sinogram that filters out, its correspondence is the 216th to the 340th sinogram in the experiment of fruit bat pupa, then, these sinograms are carried out follow-up center of rotation calculate, shown in Fig. 1 (e row).
Step S2: at each selected sinogram of step S1, the calculating barycenter is projected in a series of positions on the detector, and calculate its average, promptly obtain the center of rotation of sinogram, with the center of rotation that calculates of useful centroid trajectory average and be coarse positioning center of rotation of the present invention.Fig. 2 is based on the center of rotation coarse positioning flow process of centroid trajectory in the localization method of embodiment of the invention optical projection fault imaging center of rotation.As shown in Figure 2, at each sinogram, calculate its coarse positioning center of rotation position according to following formula
Figure GDA0000130608070000051
s c ‾ = 1 M Σ i = 0 M - 1 s m ( γ i )
Wherein M is the projection view number of optical projection fault imaging, i.e. the projected angle number of degrees; γ iBe the projection angle of optical projection fault imaging, its computing formula is
Figure GDA0000130608070000053
I=0 ..., (M-1); s mi) be that sinogram is at γ iThe angle barycenter is projected in the position on the detector, and its computing formula is as follows:
s m ( γ i ) = Σ j = 1 N s j g ( s j , γ i ) Σ j = 1 N g ( s j , γ i )
Wherein, s jBe the center of each row of detector, N is the columns of detector, g (s j, γ i) be that sinogram is at γ iAngle s jThe pixel value at row place, this step utilize the centroid method method of loci to calculate the center of rotation position of each sinogram
Figure GDA0000130608070000061
With the sinogram center of rotation that calculates of useful centroid trajectory method average, obtain the coarse positioning center of rotation of present embodiment.
Step S3: the center of rotation with step S2 coarse positioning is an initial value, center of rotation fine search interval is set, in this interval, selected sinogram is repeatedly rebuild, in a series of reconstructed images, choose the center of rotation position of reconstructed image mean square deviation maximum, what be this width of cloth sinogram decides center of rotation again, and the center of rotation of deciding again of all sinograms is averaged the meticulous position rotating of the present invention center that is.
Fig. 3 is based on the meticulous positioning flow of the center of rotation of reconstructed image mean square deviation in the localization method of embodiment of the invention optical projection fault imaging center of rotation.The first half of Fig. 3 is based on the different pairing reconstructed images in alternative center of rotation position; The latter half of Fig. 3 is the mean square deviation at the alternative center of rotation of difference position reconstructed image.From the position, the reconstruction figure of the first half, there is corresponding relation in the value of pointing to mean square deviation with its below dotted line.
As shown in Figure 3, at the center of rotation position of coarse positioning, it is the center of rotation region of search that each 20 pixel coverage is set about it, and it is 1/8 pixel value that the search stepping is set, and filters out 320 alternative center of rotation positions.At each sinogram, utilize 320 alternative center of rotation respectively it to be rebuild, obtain 320 reconstructed images, calculate the mean square deviation F of each reconstructed image Variance, its computing formula is:
F variance = Σ y = 0 N - 1 Σ x = 0 N - 1 ( f ( x , y ) - f ‾ ) 2 N 2
Wherein f (x, y) be reconstructed image coordinate (x, the pixel value of y) locating, N are the line number and the columns of reconstructed image,
Figure GDA0000130608070000063
Be the pixel average of reconstructed image, its computing formula is
Figure GDA0000130608070000064
Calculate the center of rotation position of each sinogram reconstructed image mean square deviation maximum by above-mentioned formula, what be this width of cloth sinogram decides center of rotation again, as seen is better than the reconstructed results of step S2 coarse positioning center of rotation based on the reconstructed results of deciding center of rotation again as the latter half of Fig. 3.At last, the center of rotation position that relocates of all sinograms is averaged, obtain meticulous position rotating of the present invention center.
At center of rotation computational methods of the present invention, Fig. 4 is at the center of rotation computational process sketch map of fruit bat pupa experiment in the localization method of embodiment of the invention optical projection fault imaging center of rotation.Wherein, Fig. 4 a has shown center of rotation and the average coarse positioning center of rotation that each selected sinogram calculates with the centroid trajectory method; Fig. 4 b has shown with the coarse positioning center of rotation and has rebuild the cross section that has obtained in the three-dimensional reconstruction body that only there is faint pseudo-shadow in visual picture edge, illustrates that the coarse positioning center of rotation is near true center of rotation position; Fig. 4 c has shown near the coarse positioning center with what the reconstructed image average variance method calculated and has decided center of rotation and meticulous position rotating center again, as seen meticulous position rotating center and coarse positioning center of rotation only differ from two pixels, have proved localized effectiveness of two steps from the side; Fig. 4 d has shown and utilizes meticulous position rotating center to rebuild the three-dimensional reconstruction body section that obtains that it is compared Fig. 4 b definition and is greatly improved.
Fig. 5 is at the reconstructed results of fruit bat pupa experiment in the localization method of embodiment of the invention optical projection fault imaging center of rotation.Wherein, shown partly (a)-(d) and utilized the present invention some results in the experiment of fruit bat pupa that wherein (a) shown the projection view of a width of cloth pupa; (b) for utilizing the three-dimensional reconstruction body section after center of rotation computational methods of the present invention are proofreaied and correct; (c) carry out visual result for the three-dimensional reconstruction body; (d) be the structure collection of illustrative plates of fruit bat pupa, from Fig. 5 as seen, calibrate three-dimensional reconstruction result behind the center of rotation by the present invention, the structure collection of illustrative plates of its internal structure and pupa fits like a glove, from experimentally having verified effectiveness of the present invention.
In sum, the present invention proposes a kind of localization method of optical projection fault imaging center of rotation from coarse to fine, fully ensured localized rapidity, automatic type, accuracy, specifically:
(1) on sinogram is selected, fully take into account the influence of noise of optical projection fault imaging scan-data, only select to carry out on the strong sinogram of sample signal the center of rotation search, can guarantee center of rotation search institute based on data have than high s/n ratio, improved search precision of the present invention from the data selection;
(2) in the coarse positioning process, directly utilize centroid trajectory to carry out automatic center of rotation and estimate that the barycenter method of coarse positioning process only needs data for projection, does not need human intervention, has fully guaranteed automaticity of the present invention;
(3) in thin position fixing process, result based on coarse positioning searches for once more, and employing parallel processing technique, meticulously search near being positioned at the coarse positioning value, reduce mistake of the present invention and located probability, ensured final positioning accuracy, also reduced simultaneously the hunting zone, the integrating parallel technical guarantee rapidity of the present invention.
Above-described specific embodiment; purpose of the present invention, technical scheme and beneficial effect are further described; institute is understood that; the above only is specific embodiments of the invention; be not limited to the present invention; within the spirit and principles in the present invention all, any modification of being made, be equal to replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (8)

1. the localization method of an optical projection fault imaging center of rotation is characterized in that, comprising:
From whole sinograms of optical projection fault imaging, choose the strong sinogram of signal intensity as selected sinogram, this step comprises:
Extract the signal intensity of the whole sinograms of described optical projection fault imaging, comprising:
From described optical projection fault imaging, filter out at least two width of cloth projection views;
Described at least two width of cloth projection views are carried out threshold value cut apart, remove the pixel of background area;
The sample number of pixels of every row in described at least two width of cloth projections after the calculating removal background area pixels is with the signal intensity of described sample number of pixels as this row;
The signal intensity of corresponding row in described at least two width of cloth projections is average, with described meansigma methods as signal intensity to the sinogram that should go;
By signal intensity to weak order described whole sinograms being sorted by force;
The sinogram that ordering is forward is as described selected sinogram;
Each width of cloth in several selected sinograms of optical projection fault imaging utilizes the centroid trajectory method to obtain the center of rotation of this width of cloth sinogram; The center of rotation of described several selected sinograms is averaged, obtain the coarse positioning value of optical projection fault imaging center of rotation; Wherein, described is initial value with the coarse positioning value, and the step that center of rotation fine search interval is set comprises:
With described coarse positioning value is mid point, is provided with that each N pixel coverage is center of rotation fine search interval about it;
It is the M pixel value that the search stepping is set, and filters out [2N/M] individual fine search point, and described N is greater than M, and [] is for rounding symbol;
Coarse positioning value with described center of rotation is an initial value, center of rotation fine search interval is set, in this fine search interval each width of cloth in the selected sinogram is repeatedly rebuild, in a series of reconstructed images, center of rotation is decided again as this width of cloth sinogram in the center of rotation position of choosing reconstructed image mean square deviation maximum, and the center of rotation of deciding again of all selected sinograms is averaged the meticulous position rotating of the optical projection fault imaging center that obtains.
2. localization method according to claim 1 is characterized in that, in the described step of utilizing the center of rotation that the centroid trajectory method obtains this width of cloth sinogram, and the center of rotation position of calculating sinogram according to following formula
Figure FDA00003157075900011
s c ‾ = 1 M Σ i = 0 M - 1 s m ( γ i ) ,
Wherein M is the projection view number of optical projection fault imaging, i.e. the projected angle number of degrees; γ iBe the projection angle of optical projection fault imaging, its computing formula is
Figure FDA00003157075900022
I=0 ...., (M-1); s mi) be that sinogram is at γ iThe angle barycenter is projected in the position on the detector.
3. localization method according to claim 2 is characterized in that, the center of rotation position of described calculating sinogram
Figure FDA00003157075900023
Step in, calculate described s according to following formula mi) value:
s m ( γ i ) = Σ j = 1 N s j g ( s j , γ i ) Σ j = 1 N g ( s j , γ i )
Wherein, s jBe the center of each row of detector, N is the columns of detector, g (s j, r i) be that sinogram is at γ iAngle s jThe pixel value at row place.
4. localization method according to claim 1 is characterized in that, in the described step of deciding center of rotation again of center of rotation position as this width of cloth sinogram of choosing reconstructed image mean square deviation maximum, calculates the mean square deviation of reconstructed image according to following formula:
F variance = Σ y = 0 N - 1 Σ x = 0 N - 1 ( f ( x , y ) - f ‾ ) 2 N 2
Wherein f (x, y) be reconstructed image coordinate (x, the pixel value of y) locating, N are the line number and the columns of reconstructed image,
Figure FDA00003157075900026
Be the pixel average of reconstructed image, its computing formula is
5. localization method according to claim 4, it is characterized in that, describedly selected sinogram is repeatedly rebuild in this fine search interval, in a series of reconstructed images, choose in the step of deciding center of rotation again of center of rotation position as this width of cloth sinogram of reconstructed image mean square deviation maximum, adopt parallel mode to carry out sinogram and rebuild, and obtain the mean square deviation of each reconstructed image.
6. localization method according to claim 1 is characterized in that, described N=20, described M=1/8.
7. localization method according to claim 1 is characterized in that, described at least two width of cloth projection views are two width of cloth projection views, and projection ray's direction of described two width of cloth projection views is vertical mutually.
8. localization method according to claim 1 is characterized in that, described at least two width of cloth projection views is carried out in the step that threshold value cuts apart, and described threshold value is the pixel grey scale average of projection view.
CN 201110386219 2011-11-29 2011-11-29 Method for locating optical projection tomography imaging rotation center Active CN102512140B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110386219 CN102512140B (en) 2011-11-29 2011-11-29 Method for locating optical projection tomography imaging rotation center

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110386219 CN102512140B (en) 2011-11-29 2011-11-29 Method for locating optical projection tomography imaging rotation center

Publications (2)

Publication Number Publication Date
CN102512140A CN102512140A (en) 2012-06-27
CN102512140B true CN102512140B (en) 2013-07-24

Family

ID=46283441

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110386219 Active CN102512140B (en) 2011-11-29 2011-11-29 Method for locating optical projection tomography imaging rotation center

Country Status (1)

Country Link
CN (1) CN102512140B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103033809B (en) * 2012-09-25 2014-10-15 中山大学 Method of implementing mobile separated management of compression radio frequency tomography
CN102973246B (en) * 2012-12-03 2014-12-03 中国科学院自动化研究所 Optical projection tomography method capable of automatically processing depth of field
TWI481392B (en) * 2013-02-21 2015-04-21 Univ Nat Yang Ming A dental calculus imaging methods
CN104713864B (en) * 2015-03-20 2017-06-30 中国科学院自动化研究所 The geometric correction method of the optical projection computed tomography (SPECT) system based on imitative body
CN111256663B (en) * 2018-12-03 2021-10-19 北京世纪朝阳科技发展有限公司 Centering calibration method and device
CN110068572B (en) * 2019-04-02 2021-06-29 佛山科学技术学院 Method and device for calibrating rotation center based on optical projection chromatography sinogram

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001314397A (en) * 2000-05-10 2001-11-13 Ge Medical Systems Global Technology Co Llc Method for detecting attached attitude of phantom and its x-ray ct device
CN101584587A (en) * 2008-05-21 2009-11-25 中国科学院过程工程研究所 Automatic calibration method for CT projection center
CN101825433A (en) * 2010-03-31 2010-09-08 北京航空航天大学 Measuring method of offset of rotating center of rotating table of fan beam 2D-CT scanning system
US7912174B2 (en) * 2005-10-13 2011-03-22 Agency For Science, Technology And Research Computed tomography system and method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001314397A (en) * 2000-05-10 2001-11-13 Ge Medical Systems Global Technology Co Llc Method for detecting attached attitude of phantom and its x-ray ct device
US7912174B2 (en) * 2005-10-13 2011-03-22 Agency For Science, Technology And Research Computed tomography system and method
CN101584587A (en) * 2008-05-21 2009-11-25 中国科学院过程工程研究所 Automatic calibration method for CT projection center
CN101825433A (en) * 2010-03-31 2010-09-08 北京航空航天大学 Measuring method of offset of rotating center of rotating table of fan beam 2D-CT scanning system

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Micro-CT系统中对投影图像旋转中心的校正;张蔚等;《医疗卫生装备》;20090331;第30卷(第3期);第4-6页 *
X射线CT系统投影旋转中心的测量;李保磊等;《光学精密工程》;20110531;第19卷(第5期);第967-971页 *
张蔚等.Micro-CT系统中对投影图像旋转中心的校正.《医疗卫生装备》.2009,第30卷(第3期),第4-6页.
李保磊等.X射线CT系统投影旋转中心的测量.《光学精密工程》.2011,第19卷(第5期),第967-971页.

Also Published As

Publication number Publication date
CN102512140A (en) 2012-06-27

Similar Documents

Publication Publication Date Title
CN102512140B (en) Method for locating optical projection tomography imaging rotation center
CN102753097B (en) Method for x-ray marker localization in 3d space in the presence of motion
CN102754128B (en) Landmark identification in radioscopic image and process
JP5891560B2 (en) Identification-only optronic system and method for forming three-dimensional images
US10217293B2 (en) Depth camera-based human-body model acquisition method and network virtual fitting system
CN101133431B (en) Method for registering biomedical images with reduced imaging artifacts caused by object movement
CN103900494B (en) For the homologous points fast matching method of binocular vision 3 D measurement
RU2523929C1 (en) System and method for automated planning of views in 3d brain images
CN104504710B (en) Moore stripe recognition method and device for X-ray grating phase-contrast imaging
CN103308452B (en) Optical projection tomography image capturing method based on depth-of-field fusion
CN109754445A (en) Method and apparatus for generating X-ray tomographic image data
CN109643458A (en) The system and method that sinogram for automation is completed, combines and completed by combination
CN105931202A (en) Geometric correction die body correction method and system
CN104567758B (en) Stereo imaging system and its method
CN101454801A (en) Local motion compensation based on list mode data
CN105046743A (en) Super-high-resolution three dimensional reconstruction method based on global variation technology
CN103279929A (en) CT (Computerized Tomography) image metal track prediction and artifact reduction method based on integral cosine
EP2401719B1 (en) Methods for segmenting images and detecting specific structures
CN102499701B (en) Geometrical calibrating method for X-ray and fluorescent double-mode living body imaging system
CN108186037A (en) A kind of determining method and apparatus of area-of-interest
CN104268893B (en) Method for segmenting and denoising lung parenchyma through lateral scanning and four-corner rotary scanning
CN104156977B (en) Point target movement velocity detection method based on multiple linear moveout scanning, extending and sampling
CN102488528B (en) Correcting method for geometric parameters of tomography
CN106056646B (en) CT image rebuilding methods and system
CN105869130B (en) The method of inspection and device of CT system geometric correction effect

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