JPH06274651A - Reconstituting method for density interpolation of three-dimensional image - Google Patents
Reconstituting method for density interpolation of three-dimensional imageInfo
- Publication number
- JPH06274651A JPH06274651A JP5064982A JP6498293A JPH06274651A JP H06274651 A JPH06274651 A JP H06274651A JP 5064982 A JP5064982 A JP 5064982A JP 6498293 A JP6498293 A JP 6498293A JP H06274651 A JPH06274651 A JP H06274651A
- Authority
- JP
- Japan
- Prior art keywords
- points
- density
- point
- slice images
- interpolation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Abstract
Description
【0001】[0001]
【産業上の利用分野】本発明は、互いに直交するx、
y、z座標系に置かれた3次元物体を、x、y、z方向
の各方向についてそれぞれのスライス間隔で撮影した画
像の並びである複数枚のスライス画像を受け、複数枚の
スライス画像から3次元物体を濃度値によって再構成す
る3次元画像の濃度補間による再構成方法に関し、特に
格子面画像を用いた3次元物体の濃度補間による再構成
方法に関する。BACKGROUND OF THE INVENTION The present invention relates to x, which are orthogonal to each other.
A three-dimensional object placed in the y, z coordinate system is captured at a slice interval in each of the x, y, and z directions, and a plurality of slice images, which are an array of images, are received. The present invention relates to a reconstruction method by density interpolation of a three-dimensional image for reconstructing a three-dimensional object with density values, and particularly to a reconstruction method by density interpolation of a three-dimensional object using a lattice plane image.
【0002】[0002]
【従来の技術】一般に、複数枚のスライス画像から元の
3次元物体の構造を再構成する際、スライス間隔が大き
い程、構造は粗くなる。そこで、滑らかな3次元構造を
形成するため、複数のスライス画像間において濃度補間
を行う。2. Description of the Related Art Generally, when a structure of an original three-dimensional object is reconstructed from a plurality of slice images, the larger the slice interval, the coarser the structure. Therefore, in order to form a smooth three-dimensional structure, density interpolation is performed between a plurality of slice images.
【0003】従来、複数枚のスライス画像から濃度値に
関して補間する方法は、大別して以下の2通り考えられ
た。1つ目は、図2に示すように単純に3方向の内スラ
イス間隔の一番小さい方向に着目し、その方向の2枚の
画像における濃度値を線形補間を行い、任意の注目点P
での濃度値を求める方法である。2つ目は、図3に示す
ように注目点Pを囲む格子面画像の8つの格子点に着目
し、8格子点での濃度値と、任意の注目点から各格子点
までの距離を考慮した補間を行い、点Pでの濃度値を求
める方法である。補間を周囲8個の格子点に制限するこ
とは、例えば、特開平2−58183号公報に従来技術
として開示されている。Conventionally, the following two methods are roughly considered as methods for interpolating density values from a plurality of slice images. First, as shown in FIG. 2, the focus is simply on the direction with the smallest inner slice interval in the three directions, and the density values in the two images in that direction are linearly interpolated to obtain an arbitrary target point P.
This is a method of obtaining the density value in. Secondly, as shown in FIG. 3, paying attention to eight grid points of the grid surface image surrounding the target point P, and taking into consideration the density value at the eight grid points and the distance from an arbitrary target point to each grid point. This is a method of performing the above interpolation to obtain the density value at the point P. Limiting the interpolation to eight surrounding grid points is disclosed as a conventional technique in, for example, Japanese Patent Application Laid-Open No. 2-58183.
【0004】[0004]
【発明が解決しようとする課題】この従来の3次元画像
の補間法の内、前者は1方向についての補間であるため
任意の注目点での濃度値の精度、更に再構成した3次元
物体の精度が低いという問題点があった。後者は3方向
を考慮しているため高精度ではあるが、格子点から注目
点までの距離が3次元的距離であるため、注目点を多く
取るほど距離計算時間が多くなり、更に3次元物体を再
構成するために多くの処理時間がかかるという問題があ
った。Among the conventional three-dimensional image interpolation methods, the former is interpolation in one direction, so that the accuracy of the density value at an arbitrary point of interest and the reconstructed three-dimensional object There was a problem that the accuracy was low. The latter is highly accurate because it considers three directions, but since the distance from the grid point to the point of interest is a three-dimensional distance, the more the number of points of interest, the longer the distance calculation time. There is a problem that it takes a lot of processing time to reconfigure.
【0005】本発明の課題は、高精度でかつ短縮された
処理時間で、3次元画像を濃度補間によって再構成でき
る、3次元画像の濃度補間による再構成方法を提供する
ことにある。An object of the present invention is to provide a reconstruction method by density interpolation of a three-dimensional image, which is capable of reconstructing a three-dimensional image by density interpolation with high accuracy and reduced processing time.
【0006】[0006]
【課題を解決するための手段】本発明による3次元画像
の濃度補間による再構成方法は、注目点の濃度値の補
間、更に3次元物体の再構成を行う際に、高精度かつ処
理時間低減のため、任意の注目点Pから格子面画像に引
いた、垂線が交わる6点に着目し、この6点での濃度値
と各点から注目点Pまでの距離(1次元的距離)を考慮
した補間法を用いることを特徴とする。A reconstruction method by density interpolation of a three-dimensional image according to the present invention is highly accurate and reduces the processing time when the density value of a point of interest is interpolated and the three-dimensional object is reconstructed. Therefore, pay attention to 6 points where perpendicular lines intersect from an arbitrary point of interest P drawn on the lattice plane image, and consider the density value at these 6 points and the distance (one-dimensional distance) from each point to the point of interest P. It is characterized by using the interpolation method described above.
【0007】即ち、本発明によれば、互いに直交する
x、y、z座標系に置かれた3次元物体を、x、y、z
方向の各方向についてそれぞれのスライス間隔で撮影し
た画像の並びである複数枚のスライス画像を受け、該複
数枚のスライス画像から前記3次元物体を濃度値によっ
て再構成する方法において、注目点Pを囲む格子面画像
を構成する、2枚のx方向のスライス画像、2枚のy方
向のスライス画像、及び2枚のz方向のスライス画像の
計6枚のスライス画像における濃度値より、注目点での
濃度値を補間して、3次元物体を再構成することを特徴
とする、3次元画像の濃度補間による再構成方法が得ら
れる。That is, according to the present invention, a three-dimensional object placed in an x, y, z coordinate system orthogonal to each other is converted into an x, y, z
In a method of receiving a plurality of slice images, which is an array of images taken at respective slice intervals in each direction, and reconstructing the three-dimensional object from the plurality of slice images by density values, From the density values of the two slice images in the x direction, the two slice images in the y direction, and the two slice images in the z direction, which constitute the surrounding lattice plane image, the density values in a total of six slice images There is obtained a reconstruction method by density interpolation of a three-dimensional image, which is characterized in that the three-dimensional object is reconstructed by interpolating the density values of.
【0008】更に本発明によれば、x方向に関して点P
から2枚のx方向のスライス画像に引いた垂線の交わる
点をx1 及びx2 、点x1 及びx2 における濃度値をf
x1及びfx2、点Pから点x1 及びx2 までの距離をdx1
及びdx2とし、y方向に関して点Pから2枚のy方向の
スライス画像に引いた垂線の交わる点をy1 及びy2、
点y1 及びy2 における濃度値をfy1及びfy2、点Pか
ら点y1 及びy2 までの距離をdy1及びdy2とし、z方
向に関して点Pから2枚のz方向のスライス画像に引い
た垂線の交わる点をz1 及びz2 、点z1 及びz2 にお
ける濃度値をfz1及びfz2、点Pから点z1 及びz2 ま
での距離をdz1及びdz2とすると、x方向、y方向、z
方向に関する点Pでの濃度値fpx,fpy,fpzを、各方
向の2枚の画像から fpi=(fi1×di2+fi2×di1)/Di Di =di1+di2 (ただし、i=x,y,zであり、Di は各方向の2枚
のスライス画像の間隔を表す。)に基づいて線形補間に
より求めることを特徴とする3次元画像の濃度補間によ
る再構成方法が得られる。Further in accordance with the present invention, the point P in the x direction is
X1 and x2 at the points where the perpendiculars drawn from the two slice images in the x direction intersect, and the density values at the points x1 and x2 are f
x1 and fx2, and the distance from point P to points x1 and x2 is dx1
And dx2, and the points where the perpendiculars drawn from the point P in the y direction to the two slice images in the y direction intersect are y1 and y2,
The density values at the points y1 and y2 are fy1 and fy2, the distances from the point P to the points y1 and y2 are dy1 and dy2, and the points where the perpendiculars drawn from the point P to the two z-direction slice images intersect with respect to the z direction are defined. If the density values at z1 and z2, the points z1 and z2 are fz1 and fz2, and the distances from the point P to the points z1 and z2 are dz1 and dz2, the x direction, the y direction, and the z direction.
The density values fpx, fpy, and fpz at the point P relating to the direction are calculated from two images in each direction by fpi = (fi1 × di2 + fi2 × di1) / Di Di = di1 + di2 (where i = x, y, z, Di represents a space between two slice images in each direction.), And a reconstruction method by density interpolation of a three-dimensional image is obtained which is obtained by linear interpolation.
【0009】又、本発明によれば、各方向の2枚のスラ
イス画像の間隔Di が大きくなるほど、各方向の濃度値
fpiが点Pへ影響を与える度合いが小さくなるのを、補
正するため、点Pでの濃度値fp を、 fp =(fpx/Dx +fpy/Dy +fpz/Dz )×{D
x Dy Dz /(Dx Dy +Dy Dz +Dz Dx )} に基づいて求めることを特徴とする請求項2に記載の3
次元画像の濃度補間による再構成方法が得られる。Further, according to the present invention, it is corrected that the degree to which the density value fpi in each direction influences the point P decreases as the distance Di between the two slice images in each direction increases. The density value fp at the point P is expressed by fp = (fpx / Dx + fpy / Dy + fpz / Dz) × {D
3. The method according to claim 2, wherein the value is determined based on x Dy Dz / (Dx Dy + Dy Dz + Dz Dx)}.
A reconstruction method by density interpolation of a three-dimensional image is obtained.
【0010】更に本発明によれば、3次元空間上の複数
の任意の点に注目して、各注目点Pでの濃度値fp を求
め、等濃度である点を連結することにより、滑らかな3
次元物体を再構成することを特徴とする3次元画像の濃
度補間による再構成方法が得られる。Further, according to the present invention, by paying attention to a plurality of arbitrary points in the three-dimensional space, the density value fp at each target point P is obtained, and points having the same density are connected to each other, thereby smoothing. Three
A reconstruction method by density interpolation of a three-dimensional image characterized by reconstructing a three-dimensional object is obtained.
【0011】[0011]
【実施例】以下に本発明について図面を参照して説明す
る。DESCRIPTION OF THE PREFERRED EMBODIMENTS The present invention will be described below with reference to the drawings.
【0012】図1は、本発明の一実施例による3次元画
像の濃度補間による再構成方法を説明するための図であ
る。図1は補間して求めるべき任意の注目点Pと格子面
画像との位置関係を示している。またx方向に関して点
Pから引いた垂線の交わる点がx1 ,x2 、その点にお
ける濃度値がfx1,fx2、更に点Pからx1 ,x2 まで
の距離がそれぞれdx1,dx2である。y,z方向も同様
である。FIG. 1 is a diagram for explaining a reconstruction method by density interpolation of a three-dimensional image according to an embodiment of the present invention. FIG. 1 shows the positional relationship between an arbitrary point of interest P to be obtained by interpolation and the lattice plane image. The points where the perpendiculars drawn from the point P intersect in the x direction are x1 and x2, the density values at those points are fx1 and fx2, and the distances from the point P to x1 and x2 are dx1 and dx2, respectively. The same applies to the y and z directions.
【0013】まずx方向、y方向、z方向に関する点P
での濃度値fpx,fpy,fpzを各方向の2枚の画像から
以下の式のように線形補間により求める。First, a point P concerning the x-direction, the y-direction and the z-direction.
The density values fpx, fpy, and fpz at are calculated from two images in each direction by linear interpolation as in the following equation.
【0014】fpi=(fi1×di2+fi2×di1)/Di Di =di1+di2 ただし、i=x,y,z この各方向の濃度値はスライス画像の間隔Di が大きく
なるほど、点Pへの影響の度合いは小さくなる。そこで
各方向の画像間隔を考慮すると、点Pでの濃度値fp は
次式で表される。Fpi = (fi1 × di2 + fi2 × di1) / Di Di = di1 + di2 where i = x, y, z The density value in each direction has a greater influence on the point P as the distance Di between slice images increases. Get smaller. Therefore, considering the image interval in each direction, the density value fp at the point P is expressed by the following equation.
【0015】 fp =(fpx/Dx +fpy/Dy +fpz/Dz )×{D
x Dy Dz /(Dx Dy +Dy Dz +Dz Dx )} 上式を適用して、3次元空間上の多数の任意の点につい
ての濃度値を求め、等濃度である点を連結することによ
り、滑らかな3次元物体が再構成出来る。Fp = (fpx / Dx + fpy / Dy + fpz / Dz) × {D
x Dy Dz / (Dx Dy + Dy Dz + Dz Dx)} By applying the above formula, the density values for a large number of arbitrary points in the three-dimensional space are obtained, and by connecting points of equal density, smooth points are connected. 3D objects can be reconstructed.
【0016】[0016]
【発明の効果】以上に説明したように本発明は、与えら
れたスライス画像の格子面画像の6点に着目したので、
高精度かつ少ない処理時間で元の3次元物体を再構成出
来るという効果を有する。As described above, the present invention focuses on the six points of the lattice plane image of the given slice image.
The effect is that the original three-dimensional object can be reconstructed with high accuracy and a short processing time.
【図1】本発明の一実施例による3次元画像の濃度補間
による再構成方法を説明するための図で、本発明に用い
る補間法における、3次元空間上の注目点Pと格子画面
像との位置関係を表す図である。FIG. 1 is a diagram for explaining a reconstruction method by density interpolation of a three-dimensional image according to an embodiment of the present invention, which shows an attention point P in a three-dimensional space and a grid screen image in the interpolation method used in the present invention. It is a figure showing the positional relationship of.
【図2】従来の1方向の補間法における、3次元空間上
の注目点Pとスライス画像との位置関係を表す図であ
る。FIG. 2 is a diagram showing a positional relationship between a target point P in a three-dimensional space and a slice image in a conventional one-direction interpolation method.
【図3】従来の8格子点による補間法における、3次元
空間上の注目点Pと8格子点との位置関係を表す図であ
る。FIG. 3 is a diagram showing a positional relationship between a target point P in a three-dimensional space and eight lattice points in a conventional interpolation method using eight lattice points.
P 注目点 Dx ,Dy ,Dz x,y,z各方向のスライス画像の
間隔P point of interest Dx, Dy, Dz x, y, z spacing of slice images in each direction
Claims (4)
れた3次元物体を、x、y、z方向の各方向についてそ
れぞれのスライス間隔で撮影した画像の並びである複数
枚のスライス画像を受け、該複数枚のスライス画像から
前記3次元物体を濃度値によって再構成する方法におい
て、注目点Pを囲む格子面画像を構成する、2枚のx方
向のスライス画像、2枚のy方向のスライス画像、及び
2枚のz方向のスライス画像の計6枚のスライス画像に
おける濃度値より、注目点での濃度値を補間して、3次
元物体を再構成することを特徴とする、3次元画像の濃
度補間による再構成方法。1. A plurality of slices, which are an array of images obtained by photographing a three-dimensional object placed in mutually orthogonal x, y, and z coordinate systems at respective slice intervals in each of the x, y, and z directions. In the method of receiving an image and reconstructing the three-dimensional object from the plurality of slice images by using density values, two slice images in the x direction and two y images forming a lattice plane image surrounding the target point P are formed. Characterized in that the three-dimensional object is reconstructed by interpolating the density value at the point of interest from the density values in a total of six slice images of the slice images in the direction and two slice images in the z direction. A reconstruction method of a three-dimensional image by density interpolation.
スライス画像に引いた垂線の交わる点をx1 及びx2 、
点x1 及びx2 における濃度値をfx1及びfx2、点Pか
ら点x1 及びx2 までの距離をdx1及びdx2とし、y方
向に関して点Pから2枚のy方向のスライス画像に引い
た垂線の交わる点をy1 及びy2 、点y1 及びy2 にお
ける濃度値をfy1及びfy2、点Pから点y1 及びy2 ま
での距離をdy1及びdy2とし、z方向に関して点Pから
2枚のz方向のスライス画像に引いた垂線の交わる点を
z1 及びz2 、点z1 及びz2 における濃度値をfz1及
びfz2、点Pから点z1 及びz2 までの距離をdz1及び
dz2とすると、x方向、y方向、z方向に関する点Pで
の濃度値fpx,fpy,fpzを、各方向の2枚の画像から fpi=(fi1×di2+fi2×di1)/Di Di =di1+di2 (ただし、i=x,y,zであり、Di は各方向の2枚
のスライス画像の間隔を表す。)に基づいて線形補間に
より求めることを特徴とする請求項1に記載の3次元画
像の濃度補間による再構成方法。2. The points where the perpendiculars drawn from the point P in the x direction to the two slice images in the x direction intersect are x1 and x2,
The density values at the points x1 and x2 are fx1 and fx2, the distances from the point P to the points x1 and x2 are dx1 and dx2, and the points where the perpendiculars drawn from the point P to the two y-direction slice images intersect in the y-direction Let y1 and y2 be the density values at points y1 and y2 be fy1 and fy2, and the distances from point P to points y1 and y2 be dy1 and dy2. Perpendicular lines drawn from the point P to two slice images in the z direction with respect to the z direction. Let z1 and z2 be the intersecting points, fz1 and fz2 be the density values at points z1 and z2, and dz1 and dz2 be the distances from point P to points z1 and z2. The density values fpx, fpy, and fpz are calculated from two images in each direction: fpi = (fi1 × di2 + fi2 × di1) / Di Di = di1 + di2 (where i = x, y, z, and Di is 2 in each direction). Based on the interval between slice images.) The method for reconstructing a three-dimensional image by density interpolation according to claim 1, wherein the reconstruction is performed by linear interpolation.
が大きくなるほど、各方向の濃度値fpiが点Pへ影響を
与える度合いが小さくなるのを、補正するため、点Pで
の濃度値fp を、 fp =(fpx/Dx +fpy/Dy +fpz/Dz )×{D
x Dy Dz /(Dx Dy +Dy Dz +Dz Dx )} に基づいて求めることを特徴とする請求項2に記載の3
次元画像の濃度補間による再構成方法。3. An interval Di between two slice images in each direction
In order to correct that the density value fpi in each direction has a smaller influence on the point P as becomes larger, the density value fp at the point P is calculated as fp = (fpx / Dx + fpy / Dy + fpz / Dz) × {D
3. The method according to claim 2, wherein the value is determined based on x Dy Dz / (Dx Dy + Dy Dz + Dz Dx)}.
Reconstruction method of three-dimensional image by density interpolation.
て、各注目点Pでの濃度値fp を求め、等濃度である点
を連結することにより、滑らかな3次元物体を再構成す
ることを特徴とする請求項3に記載の3次元画像の濃度
補間による再構成方法。4. A smooth three-dimensional object is reconstructed by paying attention to a plurality of arbitrary points in a three-dimensional space, obtaining a density value fp at each target point P, and connecting points having equal density. The reconstruction method by density interpolation of a three-dimensional image according to claim 3, wherein the reconstruction method is configured.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP5064982A JP2656707B2 (en) | 1993-03-24 | 1993-03-24 | Reconstruction method of 3D image by density interpolation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP5064982A JP2656707B2 (en) | 1993-03-24 | 1993-03-24 | Reconstruction method of 3D image by density interpolation |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH06274651A true JPH06274651A (en) | 1994-09-30 |
JP2656707B2 JP2656707B2 (en) | 1997-09-24 |
Family
ID=13273779
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP5064982A Expired - Lifetime JP2656707B2 (en) | 1993-03-24 | 1993-03-24 | Reconstruction method of 3D image by density interpolation |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2656707B2 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09293146A (en) * | 1996-04-26 | 1997-11-11 | Nec Corp | Visualization processing method |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0258183A (en) * | 1988-08-24 | 1990-02-27 | Gurafuika:Kk | Forming device for stereoscopic image |
-
1993
- 1993-03-24 JP JP5064982A patent/JP2656707B2/en not_active Expired - Lifetime
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0258183A (en) * | 1988-08-24 | 1990-02-27 | Gurafuika:Kk | Forming device for stereoscopic image |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09293146A (en) * | 1996-04-26 | 1997-11-11 | Nec Corp | Visualization processing method |
Also Published As
Publication number | Publication date |
---|---|
JP2656707B2 (en) | 1997-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10944960B2 (en) | Free-viewpoint video generating method and free-viewpoint video generating system | |
EP0216931B1 (en) | Image rotating system having an arbitrary angle | |
Kiryati et al. | Estimating shortest paths and minimal distances on digitized three-dimensional surfaces | |
US5309356A (en) | Three-dimensional reprojected image forming apparatus | |
US20020082800A1 (en) | Method of composing three-dimensional multi-viewpoints data | |
JPH0528243A (en) | Image-forming device | |
DE69721979T2 (en) | GRADIENT-BASED MOTION ESTIMATION | |
US20070206847A1 (en) | Correction of vibration-induced and random positioning errors in tomosynthesis | |
US5079697A (en) | Distortion reduction in projection imaging by manipulation of fourier transform of projection sample | |
JPH06274651A (en) | Reconstituting method for density interpolation of three-dimensional image | |
EP3573018A1 (en) | Image generation device, and image display control device | |
KR20090072030A (en) | An implicit geometric regularization of building polygon using lidar data | |
CN107798712B (en) | IDL polar region vector field visualization oriented algorithm | |
CN105913380B (en) | The method and apparatus of medical image conversion | |
JP2834560B2 (en) | Method for generating tomographic plane from three-dimensional image | |
JP2005312014A (en) | Resolution converting method | |
KR19980060719A (en) | Height map generation method and height map synthesis method from 2D video image | |
US6674922B1 (en) | Image processing method, image processing apparatus, and storage medium | |
JPH04106673A (en) | Moving image area dividing device | |
JP3783815B2 (en) | Image processing device | |
JPH05166706A (en) | Method for preparing data for charged particle beam lithography | |
JPS6188374A (en) | Picture data rotation processing device | |
JP2604711B2 (en) | Image conversion device | |
JP4156629B2 (en) | Three-dimensional backprojection method and X-ray CT apparatus | |
Zou et al. | Research on Surface Super Resolution for 3D optical scanning |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 19970506 |