CN108629735A - 检测图像插值及估计ct图像分辨率的方法和装置 - Google Patents
检测图像插值及估计ct图像分辨率的方法和装置 Download PDFInfo
- Publication number
- CN108629735A CN108629735A CN201710183393.8A CN201710183393A CN108629735A CN 108629735 A CN108629735 A CN 108629735A CN 201710183393 A CN201710183393 A CN 201710183393A CN 108629735 A CN108629735 A CN 108629735A
- Authority
- CN
- China
- Prior art keywords
- image
- images
- high pass
- band images
- specific direction
- 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
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000001514 detection method Methods 0.000 title claims abstract description 36
- 238000013139 quantization Methods 0.000 claims abstract description 50
- 230000009466 transformation Effects 0.000 claims abstract description 32
- 230000008859 change Effects 0.000 claims abstract description 11
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000002592 echocardiography Methods 0.000 claims 4
- 238000010586 diagram Methods 0.000 description 12
- 238000003325 tomography Methods 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000686 essence Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4007—Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种检测图像线性插值的方法和装置以及估计CT机的图像分辨率的方法和装置。根据本发明的第一方面，提供一种检测图像线性插值的方法，包括：S102，加载复数个图像；S108，在一特定方向上以一量化长度n和一预设滤波器库对各图像作小波变换，获得一相应的高通子带图像，n为一正整数；S110，检测各高通子带图像的连续性，或者显示各高通子带图像并接收对各高通子带图像的连续性的判断，若各高通子带图像存在周期性的不连续，则判断所述复数个图像在所述特定方向上有线性插值，否则改变所述量化长度n并转至步骤S108。
Description
技术领域
本发明涉及计算机断层扫描，特别是图像分辨率。
背景技术
空间分辨率(spatial resolution)又称高对比度分辨率(high contrastresolution)，它是衡量CT图像质量的一个重要参数，是测试一幅图像的量化指标，是指在高对比度(密度分辨率大10％)的情况下鉴别细微的能力，即显示最小体积病灶或结构的能力。它的定义是在两种物质CT值相差100HU以上时，能分辨最小的圆形孔径或是黑白相间(密度差相同)的线对数，单位是mm或lp/cm。其换算关系为：5÷lp/cm＝可分辨的最小物体的直径(mm)。
许多因素影响空间分辨率。探测器孔径越窄，孔径转移函数越宽，空间分辨率就越高。因焦点小的X线管产生窄的X射线，可获得较高的空间分辨率。探测器之间的距离它决定了采样间隔，间隔越小空间分辨率越高。在图像重建中选用的卷积滤波器的形式不同，空间分辨率也不同。射线剂量、矩阵、层厚、像素大小、扫描装置噪声等对空间分辨率均有影响。层厚越薄，空间分辨率越高；但层厚越薄，噪声就越大，低对比分辨率就会降低。
目前常用的检测CT空间分辨率的方法有调制传递函数(MTF)的截止频率法、分辨成排圆孔大小法和分辨线对数法。
发明内容
有鉴于此，本发明提出了一种检测图像线性插值的方法和装置以及估计CT机的图像分辨率的方法和装置。
根据本发明的第一方面，提供一种检测图像线性插值的方法，包括：S102，加载复数个图像；S108，在一特定方向上以一量化长度n和一预设滤波器库对各图像作小波变换，获得一相应的高通子带图像，n为一正整数；S110，检测各高通子带图像的连续性，或者显示各高通子带图像并接收对各高通子带图像的连续性的判断，若各高通子带图像存在周期性的不连续，则判断所述复数个图像在所述特定方向上有线性插值，否则改变所述量化长度n并转至步骤S108。
在一实施例中，所述预设滤波器库是Daubechies3小波滤波器库。
根据本发明的第二方面，提供一种估计CT机的图像分辨率的方法，包括：S102，加载复数个图像，所述图像是所述CT机的断层图像，其像素在一特定方向上占据的名义长度是Z；S108，在所述特定方向上以一量化长度n和一预设滤波器库对各图像作小波变换，获得一相应的高通子带图像，n为一正整数；S110，检测各高通子带图像的连续性，或者显示各高通子带图像并接收对各高通子带图像的连续性的判断，若各高通子带图像存在周期性的不连续，则判断所述复数个图像在所述特定方向上有线性插值，否则改变所述量化长度n并转至步骤S108；S112，估计所述CT机的图像分辨率在Z*(n-1)与Z*n之间。
在一实施例中，所述预设滤波器库是Daubechies3小波滤波器库。
根据本发明的第三方面，提供一种估计CT机的图像分辨率的方法，包括：S202，加载复数个图像，所述图像是所述CT机的断层图像，其像素在一特定方向上占据的名义长度是Z；S204，在一第一方向上以一第一量化长度m和一第一预设滤波器库对各图像作小波变换，获得一相应的第一低通子带图像，m为一正整数；S206，在一第二方向上以一第二量化长度k和一第二预设滤波器库对各第一低通子带图像作小波变换，获得一相应的第二低通子带图像，k为一正整数；S208，在一特定方向上以一量化长度n和一预设滤波器库对各第二低通子带图像作小波变换，获得一相应的高通子带图像，n为一正整数；S210，检测各高通子带图像的连续性，或者显示各高通子带图像并接收对各高通子带图像的连续性的判断，若各高通子带图像存在周期性的不连续，则判断所述复数个图像在所述特定方向上有线性插值，否则改变所述量化长度n并转至步骤S208；S212，估计所述CT机的图像分辨率在Z*(n-1)与Z*n之间。
在一实施例中，所述预设滤波器库、第一预设滤波器库和第二预设滤波器库中的至少一个是Daubechies3小波滤波器库。
在一实施例中，所述第一方向以及第二方向与所述特定方向正交。
根据本发明的第四方面，提供一种检测图像线性插值的装置，包括：一图像加载单元，其加载复数个图像；一特定方向小波变换单元，其在一特定方向上以一量化长度n和一预设滤波器库对各图像作小波变换，获得一相应的高通子带图像，n为一正整数；一检测单元，其检测各高通子带图像的连续性，或者显示各高通子带图像并接收对各高通子带图像的连续性的判断，若各高通子带图像存在周期性的不连续，则判断所述复数个图像在所述特定方向上有线性插值，否则改变所述量化长度n并转至所述特定方向小波变换单元。
在一实施例中，所述预设滤波器库是Daubechies3小波滤波器库。
根据本发明的第五方面，提供一种估计CT机的图像分辨率的装置，包括：一图像加载单元，加载复数个图像，所述图像是所述CT机的断层图像，其像素在一特定方向上占据的名义长度是Z；一特定方向小波变换单元，其在所述特定方向上以一量化长度n和一预设滤波器库对各图像作小波变换，获得一相应的高通子带图像，n为一正整数；一检测单元，其检测各高通子带图像的连续性，或者显示各高通子带图像并接收对各高通子带图像的连续性的判断，若各高通子带图像存在周期性的不连续，则判断所述复数个图像在所述特定方向上有线性插值，否则改变所述量化长度n并转至所述特定方向小波变换单元；一估计单元，其估计所述CT机的图像分辨率在Z*(n-1)与Z*n之间。
在一实施例中，所述预设滤波器库是Daubechies3小波滤波器库。
根据本发明的第六方面，提供一种估计CT机的图像分辨率的装置，包括：一图像加载单元，其加载复数个图像，所述图像是所述CT机的断层图像，其像素在一特定方向上占据的名义长度是Z；一第一方向小波变换单元，其在一第一方向上以一第一量化长度m和一第一预设滤波器库对各图像作小波变换，获得一相应的第一低通子带图像，m为一正整数；一第二方向小波变换单元，其在一第二方向上以一第二量化长度k和一第二预设滤波器库对各第一低通子带图像作小波变换，获得一相应的第二低通子带图像，k为一正整数；一特定方向小波变换单元，其在一特定方向上以一量化长度n和一预设滤波器库对各第二低通子带图像作小波变换，获得一相应的高通子带图像，n为一正整数；一检测单元，其检测各高通子带图像的连续性，或者显示各高通子带图像并接收对各高通子带图像的连续性的判断，若各高通子带图像存在周期性的不连续，则判断所述复数个图像在所述特定方向上有线性插值，否则改变所述量化长度n并转至特定方向小波变换单元；一估计单元，其估计所述CT机的图像分辨率在Z*(n-1)与Z*n之间。
在一实施例中，所述预设滤波器库、第一预设滤波器库和第二预设滤波器库中的至少一个是Daubechies3小波滤波器库。
在一实施例中，所述第一方向以及第二方向与所述特定方向正交。
本发明的检测图像线性插值的方法和装置可直接通过图像本身判断其是否存在线性插值，本发明的估计CT机的图像分辨率的方法和装置可直接通过CT断层图像本身估计CT机的图像分辨率，而无须以模体测试。
附图说明
下面将通过参照附图详细描述本发明的优选实施例，使本领域的普通技术人员更清楚本发明的上述及其它特征和优点，附图中：
图1为根据本发明的第一实施例的检测图像线性插值的方法和估计CT机的图像分辨率的方法的流程示意图。
图2为根据本发明的第二实施例的估计CT机的图像分辨率的方法的流程示意图。
图3为图2所示的方法获得的同一断层的图像的示意图。
图4A至图4F为图2所示的方法获得的几个连续断层的图像的示意图。
图5为根据本发明的第三实施例的检测图像线性插值的装置和估计CT机的图像分辨率的装置的结构示意图。
图6为根据本发明的第四实施例的估计CT机的图像分辨率的装置的结构示意图。
在上述附图中，所采用的附图标记如下：
100、200 估计CT机的图像分辨率的方法 404 第一方向小波变换单元
101 检测图像线性插值的方法 406 第二方向小波变换单元
300、400 估计CT机的图像分辨率的装置 308、408 特定方向小波变换单元
301 检测图像线性插值的装置 310、410 检测单元
302、402 图像加载单元 312、412 估计单元
S102、S108、
S110、S112、
S202、S204、 步骤
S206、S208、
S210、S212
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，以下举实施例对本发明进一步详细说明。
图1为根据本发明的第一实施例的检测图像线性插值的方法101和估计CT机的图像分辨率的方法100的流程示意图。估计CT机的图像分辨率的方法100包括步骤S102、步骤S108、步骤S110和步骤S112。
本实验采用原始图像名义长度为5mm(z方向即CT机的检查床进出机架孔的方向，也即断层方向)，0.759765625mm(x方向)，0.759765625mm(y方向)，x方向为水平方向，y方向为垂直方向，它们与z方向正交。为验证方法，采用线性插值将图像插值为各向同性，即三个方向的图像名义长度为1.5mm。实验采用提出的方法对插值后的图像检测。
在步骤S102中，加载复数个图像，这些图像是CT机的断层图像，其像素在一特定方向上占据的名义长度是Z。在本实施例中插值后图像为各向同性Z＝1.5mm。
在步骤S108中，在该特定方向上以一量化长度n和一预设滤波器库对各图像作小波变换，获得一相应的高通子带图像和一相应的低通子带图像，n为一正整数。在本实施例中，预设滤波器库是Daubechies3小波滤波器库，量化长度n＝4，因此相应的低通滤波器和高通滤波器分别为：
h[4]＝{0.4829629131445341,0.8365163037378077,0.2241438680420134,-0.1294095225512603}
g[4]＝{0.1294095225512603,0.2241438680420134,-0.8365163037378077,0.4829629131445341}。
在步骤S110中，检测各高通子带图像的连续性，或者显示各高通子带图像并接收对各高通子带图像的连续性的判断，若各高通子带图像存在周期性的不连续，则判断该复数个图像在该特定方向上有线性插值，否则改变量化长度n(例如变为n+2或n-2)并转至步骤S108。
在步骤S112中，估计CT机的图像分辨率在Z*(n-1)与Z*n之间。
若步骤S102中的图像是CT机断层图像意外的其他图像，上述部分步骤也可用于检测一般图像的线性插值。根据本发明的第一实施例，检测图像线性插值的方法101包括步骤S102、步骤S108和步骤S110。
在步骤S102中，加载复数个图像。
在步骤S108中，在一特定方向上以一量化长度n和一预设滤波器库对各图像作小波变换，获得一相应的高通子带图像和一相应的低通子带图像，n为一正整数。在本实施例中，预设滤波器库是Daubechies3小波滤波器库，量化长度n＝4，因此相应的低通滤波器和高通滤波器分别为：
h[4]＝{0.4829629131445341,0.8365163037378077,0.2241438680420134,-0.1294095225512603}
g[4]＝{0.1294095225512603,0.2241438680420134,-0.8365163037378077,0.4829629131445341}。
在步骤S110中，检测各高通子带图像的连续性，或者显示各高通子带图像并接收对各高通子带图像的连续性的判断，若各高通子带图像存在周期性的不连续，则判断该复数个图像在该特定方向上有线性插值，否则改变量化长度n(例如变为n+2或n-2)并转至步骤S108。
图2为根据本发明的第二实施例的估计CT机的图像分辨率的方法200的流程示意图，图3为图2所示的方法200获得的同一断层的图像的示意图。估计CT机的图像分辨率的方法200包括步骤S202、步骤S204、步骤S206、步骤S208、步骤S210和步骤S212。
在步骤S202中，加载复数个图像，这些图像是CT机的断层图像，其像素在一特定方向上占据的名义长度是Z。在本实施例中，该特定方向为z方向，即CT机的检查床进出机架孔的方向，也即断层方向。在本实施例中，加载153个插值后的断层图像，并且上述断层图像的像素在x和y方向上占据的名义长度分别是1.5mm和1.5mm。x方向为水平方向，y方向为垂直方向，它们与z方向正交，Z＝1.5mm。图3中的第一列示出了153个断层图像中的第38个。
在步骤S204中，在一第一方向上以一第一量化长度m和一第一预设滤波器库对各图像作小波变换，获得一相应的第一低通子带图像(L图像)和一相应的第一高通子带图像(H图像)，m为一正整数。在本实施例中，第一方向即x方向，第一量化长度m＝4。图3中的第二列示出了153个断层图像中的第38个经小波变换后获得的第一低通子带图像(L图像)和第一高通子带图像(H图像)。
在步骤S206中，在一第二方向上以一第二量化长度k和一第二预设滤波器库对各第一低通子带图像(L图像)作小波变换，获得一相应的第二低通子带图像(LL图像)和一相应的第二高通子带图像(LH图像)，k为一正整数。在本实施例中，第二方向即y方向，第二量化长度k＝4。图3中的第三列中的上两行示出了153个断层图像中的第38个的第一低通子带图像(L图像)经小波变换后获得的第二低通子带图像(LL图像)和一相应的第二高通子带图像(LH图像)。
在步骤S208中，在一特定方向上以一量化长度n和一预设滤波器库对各第二低通子带图像(LL图像)作小波变换，获得一相应的高通子带图像(LLH图像)和一相应的低通子带图像(LLL图像)，n为一正整数。在本实施例中，该特定方向为z方向，量化长度n＝4。图3中的第四、五列中的第一行示出了153个断层图像中的第38个的第二低通子带图像(LL图像)经小波变换后获得的低通子带图像(LLL图像)和高通子带图像(LLH图像)。153个断层图像可获得77个低通子带图像(LLL图像)和77个高通子带图像(LLH图像)，第38个断层图像对应的低通子带图像(LLL图像)和高通子带图像(LLH图像)排列第19。
在步骤S210中，检测各高通子带图像(LLH图像)的连续性，或者显示各高通子带图像(LLH图像)并接收对各高通子带图像(LLH图像)的连续性的判断，若各高通子带图像存在周期性的不连续，则判断复数个图像在该特定方向上有线性插值，否则改变量化长度n(例如变为n+2或n-2)并转至步骤S208。图4A至图4F为图2所示的方法100获得的几个连续断层的图像的示意图。图4A至图4F分别是77个高通子带图像(LLH图像)中的第23至28个。若观察全部77个高通子带图像，可发现每5个高通子带图像出现一个低响应的切片，即各高通子带图像存在周期性的不连续，因此可判断这些图像在该特定方向(z方向)上有线性插值。
在步骤S212中，估计CT机的图像分辨率在Z*(n-1)与Z*n之间。在本实施例中，估计z方向的图像分辨率在1.5*(4-1)mm与1.5*4mm之间，即在4.5mm至6mm之间。而实际上插值前的z方向图像分辨率为5mm。
在本实施例中，预设滤波器库、第一预设滤波器库和第二预设滤波器库中是Daubechies3小波滤波器库，因此相应的低通滤波器和高通滤波器分别为：
h[4]＝{0.4829629131445341,0.8365163037378077,0.2241438680420134,-0.1294095225512603}
g[4]＝{0.1294095225512603,0.2241438680420134,-0.8365163037378077,0.4829629131445341}。
图5为根据本发明的第三实施例的检测图像线性插值的装置301和估计CT机的图像分辨率的装置300的结构示意图。如图5所示，估计CT机的图像分辨率的装置300包括一图像加载单元302、一特定方向小波变换单元308、一检测单元310和一估计单元312。
本实验采用原始图像名义长度为5mm(z方向即CT机的检查床进出机架孔的方向，也即断层方向)，0.759765625mm(x方向)，0.759765625mm(y方向)，x方向为水平方向，y方向为垂直方向，它们与z方向正交。为验证估计CT机的图像分辨率的装置300，采用线性插值将图像插值为各向同性，即三个方向的图像名义长度为1.5mm。实验装置300对插值后的图像检测。
图像加载单元302加载复数个图像，这些图像是CT机的断层图像，其像素在一特定方向上占据的名义长度是Z。在本实施例中插值后图像为各向同性Z＝1.5mm。
特定方向小波变换单元308在该特定方向上以一量化长度n和一预设滤波器库对各图像作小波变换，获得一相应的高通子带图像和一相应的低通子带图像，n为一正整数。在本实施例中，预设滤波器库是Daubechies3小波滤波器库，量化长度n＝4，因此相应的低通滤波器和高通滤波器分别为：
h[4]＝{0.4829629131445341,0.8365163037378077,0.2241438680420134,-0.1294095225512603}
g[4]＝{0.1294095225512603,0.2241438680420134,-0.8365163037378077,0.4829629131445341}。
检测单元310检测各高通子带图像的连续性，或者显示各高通子带图像并接收对各高通子带图像的连续性的判断，若各高通子带图像存在周期性的不连续，则判断该复数个图像在该特定方向上有线性插值，否则改变量化长度n(例如变为n+2或n-2)并转至特定方向小波变换单元308。
估计单元312估计CT机的图像分辨率在Z*(n-1)与Z*n之间。
若图像是CT机断层图像意外的其他图像，上述部分模块也可用于检测一般图像的线性插值。根据本发明的第三实施例，检测图像线性插值的装置301包括一图像加载单元302、一特定方向小波变换单元308和一检测单元310。：
图像加载单元302加载复数个图像。
特定方向小波变换单元308在一特定方向上以一量化长度n和一预设滤波器库对各图像作小波变换，获得一相应的高通子带图像和一相应的低通子带图像，n为一正整数。在本实施例中，预设滤波器库是Daubechies3小波滤波器库，量化长度n＝4，因此相应的低通滤波器和高通滤波器分别为：
h[4]＝{0.4829629131445341,0.8365163037378077,0.2241438680420134,-0.1294095225512603}
g[4]＝{0.1294095225512603,0.2241438680420134,-0.8365163037378077,0.4829629131445341}。
检测单元310检测各高通子带图像的连续性，或者显示各高通子带图像并接收对各高通子带图像的连续性的判断，若各高通子带图像存在周期性的不连续，则判断该复数个图像在该特定方向上有线性插值，否则改变量化长度n(例如变为n+2或n-2)并转至特定方向小波变换单元308。
图6为根据本发明的第四实施例的估计CT机的图像分辨率的装置400的结构示意图，图3为图6所示的装置400获得的同一断层的图像的示意图。估计CT机的图像分辨率的装置400包括一图像加载单元402、一第一方向小波变换单元404、一第二方向小波变换单元406、一特定方向小波变换单元408、一检测单元410和一估计单元412。
图像加载单元402加载复数个图像，这些图像是CT机的断层图像，其像素在一特定方向上占据的名义长度是Z。在本实施例中，该特定方向为z方向，即CT机的检查床进出机架孔的方向，也即断层方向。在本实施例中，加载153个插值后的断层图像，并且上述断层图像的像素在x和y方向上占据的名义长度分别是1.5mm和1.5mm。x方向为水平方向，y方向为垂直方向，它们与z方向正交，Z＝1.5mm。图3中的第一列示出了153个断层图像中的第38个。
第一方向小波变换单元404在一第一方向上以一第一量化长度m和一第一预设滤波器库对各图像作小波变换，获得一相应的第一低通子带图像(L图像)和一相应的第一高通子带图像(H图像)，m为一正整数。在本实施例中，第一方向即x方向，第一量化长度m＝4。图3中的第二列示出了153个断层图像中的第38个经小波变换后获得的第一低通子带图像(L图像)和第一高通子带图像(H图像)。
第二方向小波变换单元406在一第二方向上以一第二量化长度k和一第二预设滤波器库对各第一低通子带图像(L图像)作小波变换，获得一相应的第二低通子带图像(LL图像)和一相应的第二高通子带图像(LH图像)，k为一正整数。在本实施例中，第二方向即y方向，第二量化长度k＝4。图3中的第三列中的上两行示出了153个断层图像中的第38个的第一低通子带图像(L图像)经小波变换后获得的第二低通子带图像(LL图像)和一相应的第二高通子带图像(LH图像)。
特定方向小波变换单元408在一特定方向上以一量化长度n和一预设滤波器库对各第二低通子带图像(LL图像)作小波变换，获得一相应的高通子带图像(LLH图像)和一相应的低通子带图像(LLL图像)，n为一正整数。在本实施例中，该特定方向为z方向，量化长度n＝4。图3中的第四、五列中的第一行示出了153个断层图像中的第38个的第二低通子带图像(LL图像)经小波变换后获得的低通子带图像(LLL图像)和高通子带图像(LLH图像)。153个断层图像可获得77个低通子带图像(LLL图像)和77个高通子带图像(LLH图像)，第38个断层图像对应的低通子带图像(LLL图像)和高通子带图像(LLH图像)排列第19。
检测单元410检测各高通子带图像(LLH图像)的连续性，或者显示各高通子带图像(LLH图像)并接收对各高通子带图像(LLH图像)的连续性的判断，若各高通子带图像存在周期性的不连续，则判断复数个图像在该特定方向上有线性插值，否则改变量化长度n(例如变为n+2或n-2)并转至步骤S208。图4A至图4F为图2所示的方法100获得的几个连续断层的图像的示意图。图4A至图4F分别是77个高通子带图像(LLH图像)中的第23至28个。若观察全部77个高通子带图像，可发现每5个高通子带图像出现一个低响应的切片，即各高通子带图像存在周期性的不连续，因此可判断这些图像在该特定方向(z方向)上有线性插值。
估计单元412估计CT机的图像分辨率在Z*(n-1)与Z*n之间。在本实施例中，估计Z方向的图像分辨率在1.5*(4-1)mm与1.5*4mm之间，即在4.5mm至6mm之间。而实际上插值前的z方向图像分辨率为5mm。
在本实施例中，预设滤波器库、第一预设滤波器库和第二预设滤波器库中是Daubechies3小波滤波器库，因此相应的低通滤波器和高通滤波器分别为：
h[4]＝{0.4829629131445341,0.8365163037378077,0.2241438680420134,-0.1294095225512603}
g[4]＝{0.1294095225512603,0.2241438680420134,-0.8365163037378077,0.4829629131445341}。
本发明的检测图像线性插值的方法和装置可直接通过图像本身判断其是否存在线性插值，本发明的估计CT机的图像分辨率的方法和装置可直接通过CT断层图像本身估计CT机的图像分辨率，而无须以模体测试。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (14)
1.一种检测图像线性插值的方法(101)，包括：
S102，加载复数个图像；
S108，在一特定方向上以一量化长度n和一预设滤波器库对各图像作小波变换，获得一相应的高通子带图像，n为一正整数；
S110，检测各高通子带图像的连续性，或者显示各高通子带图像并接收对各高通子带图像的连续性的判断，若各高通子带图像存在周期性的不连续，则判断所述复数个图像在所述特定方向上有线性插值，否则改变所述量化长度n并转至步骤S108。
2.如权利要求1所述的方法(101)，其特征是，所述预设滤波器库是Daubechies3小波滤波器库。
3.一种估计CT机的图像分辨率的方法(100)，包括：
S102，加载复数个图像，所述图像是所述CT机的断层图像，其像素在一特定方向上占据的名义长度是Z；
S108，在所述特定方向上以一量化长度n和一预设滤波器库对各图像作小波变换，获得一相应的高通子带图像，n为一正整数；
S110，检测各高通子带图像的连续性，或者显示各高通子带图像并接收对各高通子带图像的连续性的判断，若各高通子带图像存在周期性的不连续，则判断所述复数个图像在所述特定方向上有线性插值，否则改变所述量化长度n并转至步骤S108；
S112，估计所述CT机的图像分辨率在Z*(n-1)与Z*n之间。
4.如权利要求3所述的方法(100)，其特征是，所述预设滤波器库是Daubechies3小波滤波器库。
5.一种估计CT机的图像分辨率的方法(200)，包括：
S202，加载复数个图像，所述图像是所述CT机的断层图像，其像素在一特定方向上占据的名义长度是Z；
S204，在一第一方向上以一第一量化长度m和一第一预设滤波器库对各图像作小波变换，获得一相应的第一低通子带图像，m为一正整数；
S206，在一第二方向上以一第二量化长度k和一第二预设滤波器库对各第一低通子带图像作小波变换，获得一相应的第二低通子带图像，k为一正整数；
S208，在一特定方向上以一量化长度n和一预设滤波器库对各第二低通子带图像作小波变换，获得一相应的高通子带图像，n为一正整数；
S210，检测各高通子带图像的连续性，或者显示各高通子带图像并接收对各高通子带图像的连续性的判断，若各高通子带图像存在周期性的不连续，则判断所述复数个图像在所述特定方向上有线性插值，否则改变所述量化长度n并转至步骤S208；
S212，估计所述CT机的图像分辨率在Z*(n-1)与Z*n之间。
6.如权利要求5所述的方法(200)，其特征是，所述预设滤波器库、第一预设滤波器库和第二预设滤波器库中的至少一个是Daubechies3小波滤波器库。
7.如权利要求5所述的方法(200)，其特征是，所述第一方向以及第二方向与所述特定方向正交。
8.一种检测图像线性插值的装置(301)，包括：
一图像加载单元(302)，其加载复数个图像；
一特定方向小波变换单元(308)，其在一特定方向上以一量化长度n和一预设滤波器库对各图像作小波变换，获得一相应的高通子带图像，n为一正整数；
一检测单元(310)，其检测各高通子带图像的连续性，或者显示各高通子带图像并接收对各高通子带图像的连续性的判断，若各高通子带图像存在周期性的不连续，则判断所述复数个图像在所述特定方向上有线性插值，否则改变所述量化长度n并转至所述特定方向小波变换单元(308)。
9.如权利要求8所述的装置(301)，其特征是，所述预设滤波器库是Daubechies3小波滤波器库。
10.一种估计CT机的图像分辨率的装置(300)，包括：
一图像加载单元(302)，加载复数个图像，所述图像是所述CT机的断层图像，其像素在一特定方向上占据的名义长度是Z；
一特定方向小波变换单元(308)，其在所述特定方向上以一量化长度n和一预设滤波器库对各图像作小波变换，获得一相应的高通子带图像，n为一正整数；
一检测单元(310)，其检测各高通子带图像的连续性，或者显示各高通子带图像并接收对各高通子带图像的连续性的判断，若各高通子带图像存在周期性的不连续，则判断所述复数个图像在所述特定方向上有线性插值，否则改变所述量化长度n并转至所述特定方向小波变换单元(308)；
一估计单元(312)，其估计所述CT机的图像分辨率在Z*(n-1)与Z*n之间。
11.如权利要求10所述的装置(300)，其特征是，所述预设滤波器库是Daubechies3小波滤波器库。
12.一种估计CT机的图像分辨率的装置(400)，包括：
一图像加载单元(402)，其加载复数个图像，所述图像是所述CT机的断层图像，其像素在一特定方向上占据的名义长度是Z；
一第一方向小波变换单元(404)，其在一第一方向上以一第一量化长度m和一第一预设滤波器库对各图像作小波变换，获得一相应的第一低通子带图像，m为一正整数；
一第二方向小波变换单元(406)，其在一第二方向上以一第二量化长度k和一第二预设滤波器库对各第一低通子带图像作小波变换，获得一相应的第二低通子带图像，k为一正整数；
一特定方向小波变换单元(408)，其在一特定方向上以一量化长度n和一预设滤波器库对各第二低通子带图像作小波变换，获得一相应的高通子带图像，n为一正整数；
一检测单元(410)，其检测各高通子带图像的连续性，或者显示各高通子带图像并接收对各高通子带图像的连续性的判断，若各高通子带图像存在周期性的不连续，则判断所述复数个图像在所述特定方向上有线性插值，否则改变所述量化长度n并转至特定方向小波变换单元(408)；
一估计单元(412)，其估计所述CT机的图像分辨率在Z*(n-1)与Z*n之间。
13.如权利要求12所述的装置(400)，其特征是，所述预设滤波器库、第一预设滤波器库和第二预设滤波器库中的至少一个是Daubechies3小波滤波器库。
14.如权利要求12所述的装置(400)，其特征是，所述第一方向以及第二方向与所述特定方向正交。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710183393.8A CN108629735B (zh) | 2017-03-24 | 2017-03-24 | 检测图像插值及估计ct图像分辨率的方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710183393.8A CN108629735B (zh) | 2017-03-24 | 2017-03-24 | 检测图像插值及估计ct图像分辨率的方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108629735A true CN108629735A (zh) | 2018-10-09 |
| CN108629735B CN108629735B (zh) | 2021-10-22 |
Family
ID=63707740
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710183393.8A Active CN108629735B (zh) | 2017-03-24 | 2017-03-24 | 检测图像插值及估计ct图像分辨率的方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108629735B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000216981A (ja) * | 1999-01-25 | 2000-08-04 | Sony Corp | 電子透かし埋め込み方法および電子透かし埋め込み装置 |
| CN101082983A (zh) * | 2007-06-22 | 2007-12-05 | 崔志明 | 基于感兴趣区域的自适应医学序列图像插值方法 |
| US8744205B2 (en) * | 2011-06-16 | 2014-06-03 | Fraunhofer Gesellschaft zur Förderung der angewandten Forschung e.V. | Apparatus for decomposing images and recomposing adjusted images |
| CN103841421A (zh) * | 2014-02-26 | 2014-06-04 | 吉林省大为科技有限公司 | 一种usdc数据处理方法 |
| US20140233820A1 (en) * | 2012-11-01 | 2014-08-21 | Virginia Commonweath University | Segmentation and Fracture Detection in CT Images |
| US20150348288A1 (en) * | 2014-05-30 | 2015-12-03 | Ge Medical Systems Global Technology Company, Llc | Image producing apparatus, radiation tomography apparatus, and program |
-
2017
- 2017-03-24 CN CN201710183393.8A patent/CN108629735B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000216981A (ja) * | 1999-01-25 | 2000-08-04 | Sony Corp | 電子透かし埋め込み方法および電子透かし埋め込み装置 |
| CN101082983A (zh) * | 2007-06-22 | 2007-12-05 | 崔志明 | 基于感兴趣区域的自适应医学序列图像插值方法 |
| US8744205B2 (en) * | 2011-06-16 | 2014-06-03 | Fraunhofer Gesellschaft zur Förderung der angewandten Forschung e.V. | Apparatus for decomposing images and recomposing adjusted images |
| US20140233820A1 (en) * | 2012-11-01 | 2014-08-21 | Virginia Commonweath University | Segmentation and Fracture Detection in CT Images |
| CN103841421A (zh) * | 2014-02-26 | 2014-06-04 | 吉林省大为科技有限公司 | 一种usdc数据处理方法 |
| US20150348288A1 (en) * | 2014-05-30 | 2015-12-03 | Ge Medical Systems Global Technology Company, Llc | Image producing apparatus, radiation tomography apparatus, and program |
Non-Patent Citations (1)
| Title |
|---|
| SUPRIYA S. JARANDE ETC.: "Comparative Analysis of Image Enhancement", 《INTERNATIONAL CONFERENCE ON ELECTRICAL, ELECTRONICS, AND OPTIMIZATION TECHNIQUES》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108629735B (zh) | 2021-10-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Sharma et al. | Clinical image quality assessment of accelerated magnetic resonance neuroimaging using compressed sensing | |
| US8208709B2 (en) | System and method for improved real-time cine imaging | |
| EP1388817A1 (en) | Three-dimensional spatial filtering apparatus and method | |
| Li et al. | An effective noise reduction method for multi‐energy CT images that exploit spatio‐spectral features | |
| Pineda et al. | Beyond noise power in 3D computed tomography: the local NPS and off‐diagonal elements of the Fourier domain covariance matrix | |
| CN107656224B (zh) | 一种磁共振成像方法、装置及系统 | |
| Goerner et al. | A comparison of five standard methods for evaluating image intensity uniformity in partially parallel imaging MRI | |
| CN102525496A (zh) | 确定管电压值的方法、x射线装置、计算程序和数据载体 | |
| US5722410A (en) | Apparatus and method for automated determination of surface coil coordinates to correct signal inhomogeneities in magnetic resonance signals | |
| CN108629735A (zh) | 检测图像插值及估计ct图像分辨率的方法和装置 | |
| JP5819675B2 (ja) | 磁気共鳴イメージング装置、およびプログラム | |
| US9355437B2 (en) | Image position compensation method and apparatus | |
| KR101431646B1 (ko) | 데이터처리장치, 데이터처리방법 및 컴퓨터 판독가능 기록매체 | |
| Ji et al. | Compressed sensing parallel magnetic resonance imaging | |
| JP2003159246A (ja) | 医用診断撮像システムの像表示の質検査方法及び装置 | |
| Villanueva et al. | Medical diagnosis improvement through image quality enhancement based on super-resolution | |
| CN111145285B (zh) | Mri动态增强图像重建方法、装置、计算机设备和存储介质 | |
| CN113255515A (zh) | 发动机振动形态可视化方法、存储介质和电子设备 | |
| Pal et al. | Analysis of noise power spectrum for linear and non-linear reconstruction algorithms for CT | |
| Körner et al. | Sinogram interpolation to decrease acquisition time in X-ray computed tomography measurement of surface topography | |
| CN112712572A (zh) | Ct扫描设备的低信号噪声的抑制方法、系统及计算机可读存储介质 | |
| Miéville et al. | 3D noise power spectrum applied on clinical MDCT scanners: Effects of reconstruction algorithms and reconstruction filters | |
| Paskaš et al. | Edge examination using Hölder exponent and image statistics | |
| Kacperski | Attenuation correction in SPECT without attenuation map | |
| CN114879107B (zh) | 一种fMRI中射频接收线圈本征时域稳定性参数的测量方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |