CN101627918A - Ct影像压缩方法和装置 - Google Patents
Ct影像压缩方法和装置 Download PDFInfo
- Publication number
- CN101627918A CN101627918A CN200810137772A CN200810137772A CN101627918A CN 101627918 A CN101627918 A CN 101627918A CN 200810137772 A CN200810137772 A CN 200810137772A CN 200810137772 A CN200810137772 A CN 200810137772A CN 101627918 A CN101627918 A CN 101627918A
- Authority
- CN
- China
- Prior art keywords
- image
- compression
- classification model
- data
- value
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 66
- 238000007906 compression Methods 0.000 claims abstract description 144
- 230000006835 compression Effects 0.000 claims abstract description 144
- 230000000007 visual effect Effects 0.000 claims abstract description 53
- 238000013145 classification model Methods 0.000 claims description 96
- 230000006837 decompression Effects 0.000 claims description 20
- 238000012217 deletion Methods 0.000 claims description 10
- 230000037430 deletion Effects 0.000 claims description 10
- 238000013480 data collection Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 6
- 229940079593 drug Drugs 0.000 claims description 5
- 239000003814 drug Substances 0.000 claims description 5
- 238000011084 recovery Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000013500 data storage Methods 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 3
- 238000002591 computed tomography Methods 0.000 description 108
- 238000003759 clinical diagnosis Methods 0.000 description 9
- 210000001519 tissue Anatomy 0.000 description 9
- 241001269238 Data Species 0.000 description 8
- 210000004369 blood Anatomy 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 210000003710 cerebral cortex Anatomy 0.000 description 2
- 238000005056 compaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 210000004185 liver Anatomy 0.000 description 2
- 210000003205 muscle Anatomy 0.000 description 2
- 210000000496 pancreas Anatomy 0.000 description 2
- 210000002381 plasma Anatomy 0.000 description 2
- 230000002685 pulmonary effect Effects 0.000 description 2
- 210000000952 spleen Anatomy 0.000 description 2
- 210000001685 thyroid gland Anatomy 0.000 description 2
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 1
- CPGKMLVTFNUAHL-UHFFFAOYSA-N [Ca].[Ca] Chemical compound [Ca].[Ca] CPGKMLVTFNUAHL-UHFFFAOYSA-N 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 210000001185 bone marrow Anatomy 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 229910052791 calcium Inorganic materials 0.000 description 1
- 239000011575 calcium Substances 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 230000002490 cerebral effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M7/00—Conversion of a code where information is represented by a given sequence or number of digits to a code where the same, similar or subset of information is represented by a different sequence or number of digits
- H03M7/30—Compression; Expansion; Suppression of unnecessary data, e.g. redundancy reduction
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
- H04N19/60—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using transform coding
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Compression Of Band Width Or Redundancy In Fax (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
本发明公开了CT影像压缩方法和装置。本发明提供的对CT重建后的图像进行压缩的方法,包括如下的步骤:获得要压缩的重建后的图像;所述的重建后的图像包括视野外区域和视野内区域;根据预设的CT值分类模板,判断图像中的视野内区域各个像素的CT值属于哪类CT值分类模板;根据CT值分类模板和各分类模板中所确定的压缩方法对落于其中的像素数据进行相应的压缩。本发明通过利用CT值的分类的方式,对于重要的图像数据进行无损压缩,对那些对临床诊断不是很敏感的图像数据进行不同程度的有损压缩,不仅有效地保证了压缩后的数据的有效性,还大大地提高的压缩比例,增强了CT系统的影像存储能力。
Description
技术领域
本发明涉及影像压缩技术,特别是关于计算机断层摄影(CT)的影像压缩。
背景技术
在CT系统中,影像储存容量是一个非常重要的系统性能参数。为了提高影像储存容量,现有的方法是通过增加硬盘,但这同时也会增加成本。于是出现了通过压缩医学图像文件格式(Digital Imaging andCommunications in Medicine,简称DICOM)图像的方法来提高CT系统的影像储存容量。典型的影像压缩方法有两种:一种是无损压缩,另外一种是有损压缩。无损压缩包括Huffman(霍夫曼)编码、algorithm(算术)编码等压缩方法。这种压缩率比较低。有损压缩包括DCT(离散余弦)编码、Predictive(预测)编码、vector quantization(量化)等等压缩方法。有损压缩可以得到较高的压缩比率。现有技术如西门子公司申请的CN200610089179所揭示的压缩方法,其对图像的编码需要用到相邻图像的数据,如果相邻图像相关性不大时,编码效果就会较差。
发明内容
本发明的目的在于提出一种CT影像的压缩方法及装置,用于压缩存储该CT图像所需的信息量,从而提高CT系统的影像储存容量。
本发明提出一种对CT重建后的图像进行压缩的方法,包括如下的步骤:
步骤10)获得要压缩的重建后的图像;所述的重建后的图像包括视野外区域和视野内区域;
步骤11)根据预设的CT值分类模板,判断图像中的视野内区域各个像素的CT值属于哪类CT值分类模板;
步骤12)根据CT值分类模板和各分类模板中所确定的压缩方法对落于其中的像素数据进行相应的压缩。
其中,还包括如下的步骤:除去视野外区域数据,并建立相应的记录信息。
所述的重建后的图像为医学图像文件格式(Digital Imaging andCommunications in Medicine,简称DICOM)图像,其为512X512的矩阵。
所述的预设的CT值分类模板包括A[-1024~-111],B[-110~145],C[146~657],D[658~3071]四类;对CT值分类模板B[-110~145]内的图像数据采用无损压缩,对CT值分类模板A[-1024~-111],C[146~657],D[658~3071]的图像数据采用有损压缩。
所述的A、B、C、D四类的CT值分类模板中的所对应的压缩方式是按照如下的公式进行: 其中,y为压缩结果,x为各像素的CT值,intercept为每一分类模板中的下限值;slope为每一分类模板中的压缩比率,为常数,用户可以根据自己的要求定义压缩比率;其中,对应于CT值分类模板B[-110~145]的压缩比率slope为1;对应于CT值分类模板A[-1024~-111],C[146~657],D[658~3071]压缩比率slope分别为4,2,9。
A、B、C、D四类CT值分类模板分别采用二进制表示,A用00,B用01,C用10,D用11表示。
视野外区域像素的CT值为-3024。
本发明提出一种对上述压缩后的CT重建图像进行解压缩的方法,包括如下的步骤:
21)获得要解压的影像数据;
22)根据每个像素数据所在的CT值分类模板,应用与之相对应的解压算法对压缩的影像数据进行解压缩;
23)恢复视野外区域的图像数据;
24)将视野外区域数据加入到解压缩后的影像中,得到完整的图像。
其中,应用如下的公式对压缩的影像数据进行解压缩x=y·slope+intercept,其中x表示解压缩后的影像数据,y为被压缩的影像数据,压缩比率slope和intercept由各自所处的CT值分类模板决定。
本发明还提供一种CT装置,包括:扫描台,其对受检者进行X射线扫描;数据收集单元,其对扫描台输出的扫描数据进行收集、模数转换;图像重建单元,其根据数据收集单元发送来的扫描数据进行图像重建,并将所重构建得到的图像数据存储到存储单元和/或通过中央控制单元在用户操作界面上显示;供用户操作的用户操作界面;用于存储数据和信息的存储单元;和中央控制单元,其连接至所述图像重建单元的输出端,并控制扫描台、数据收集单元和用户操作界面、存储单元;所述的存储单元中存储有CT值分类模板及应用于CT值分类模板的压缩算法和解压缩算法;该CT装置还包括图像压缩/解压计算单元,其在中央控制单元的控制下调用存储在存储单元的CT值分类模板及相应的压缩算法和解压缩算法,对重建后的图像进行压缩或对压缩后的图像进行解压缩。
所述的图像压缩/解压计算单元包括压缩计算单元和解压计算单元;中央控制单元可以自动的控制图像压缩/解压计算单元中的压缩计算单元对图像重建单元重建后的图像进行压缩,也可以通过用户操作用户操作界面向中央控制单元发出相应的压缩指令。
所述的所述重建单元重建后的图像为医学图像文件格式(DigitalImaging and Communications in Medicine,简称DICOM)图像,其为512×512的矩阵,包括视野内区域和视野外区域。
所述的CT值分类模板包括A[-1024~-111],B[-110~145],C[146~657],D[658~3071]四类;对CT值分类模板B[-110~145]内的图像数据采用无损压缩,对CT值分类模板A[-1024~-111],C[146~657],D[658~3071]的图像数据采用有损压缩。
A、B、C、D四类的CT值分类模板中的所对应的压缩方式是按照如下的公式进行: 其中,y为压缩结果,x为各像素的CT值,intercept为每一分类模板中的下限值;slope为每一分类模板中的压缩比率,为常数,用户可以根据自己的要求定义压缩比率;其中,对应于CT值分类模板B[-110~145]的压缩比率slope为1;与之相对应的对压缩的影像数据进行解压缩的方法为x=y·slope+intercept,其中x表示解压缩后的影像数据,y为被压缩的影像数据,压缩比率slope和intercept由各自所处的CT值分类模板决定;y·slope表示y与slope的乘积。
对应于CT值分类模板A[-1024~-111],C[146~657],D[658~3071]的压缩比率slope分别为4,2,9。
A、B、C、D四类CT值分类模板分别采用二进制表示,A用00、B用01、C用10、D用11表示。
当压缩计算单元收到中央控制单元发送来的对重建图像进行压缩的指令时,从存储单元中调用要压缩的重建图像数据,并将图像中视野外区域数据删除,判断视野内的图像数据的各个像素CT值属于那个CT值分类模板,并根据各像素所落入CT值分类模板的压缩算法对像素数据进行压缩;所述的用户操作界面上设置有“删除空气数据”的操作键。
当解压计算单元收到中央控制单元发送来的对压缩后的重建图像进行解压缩的指令时,从存储单元中调用要解压缩的图像数据,根据各像素所落入CT值分类模板的解压缩算法对压缩的像素数据进行解压缩;并恢复视野外区域的图像数据;将视野外区域数据加入到解压缩后的影像中,得到完整的图像。
当压缩计算单元删除视野外区域数据时,其同时在存储单元中建立、存储一条对应的记录信息;压缩计算单元压缩后的图像数据或解压计算单元解压后的图像数据被存储到存储单元中。
本发明通过利用CT值分类的方式,用户可以根据自己需求,对于重要的图像数据进行无损压缩,对那些对临床诊断不是很敏感的图像数据进行不同程度的有损压缩,确保重要组织信息的数据是没有丢失的,对于那些丢失部分数据的包括空气、骨骼,其对诊断是不会产生多大的影响;本发明通过按照CT值分类模板的方式对不同的CT值像素数据进行的不同的压缩,不仅有效的保证了压缩后的数据的有效性,还大大的提高了压缩比例,增强了CT系统影像存储能力,提高了影像传送的能力。同时,本发明的压缩方法可以和其他方法联合使用,达到更高的压缩比,兼容性好。
附图说明
图1为本发明提出的一种压缩图像数据的流程图;
图2为CT装置重建后得到的一种DICOM影像的示意图;
图3a为人类组织CT值的范围的示意图;
图3b为表示图2中DICOM影像的各个像素所在的CT值分类模板的一实施例;
图4a为对C分类模板中的像素数据进行压缩的示意图;
图4b为对A,B,D分类模板中的像素数据分别进行压缩的示意图;
图5为对利用图1中的压缩方法压缩后的图像数据进行解压缩的方法流程图;
图6为通过用户操作后的DICOM图像示意图;
图7为具有压缩CT重建图像功能的CT装置的功能模块示意图。
具体实施方式
下面,参照附图来详细说明本发明的实施形态。本发明不限于实施形态。
本发明的CT影像压缩方法是基于影像中物体组织的CT值如人体的器官、组织的CT值对DICOM图像进行有损或无损压缩,从而不仅保证重要图像数据不会丢失,还较大的提高了CT系统的影像储存容量。
请参照图1所示,为本发明CT影像压缩方法的流程图,结合人体为例说明本发明的压缩方法:
步骤10)获得要压缩的DICOM图像数据。
请同时参阅图2,为一个CT DICOM图像;DICOM图像数据通常包括两部分的内容:头信息和影像数据。头信息记录了病人的信息、扫描协议等。一般的CT影像数据是512×512像素,每一像素通常采用2个字节(16位)的方式表示。图2所示的是一个典型的512×512 CTDICOM图像,它的尺寸是527672字节,包括3384个字节头信息和524288个字节影像数据。本发明中的压缩仅对影像数据进行压缩。
步骤11)除去视野外区域数据。
从图2可以看出,DICOM图像数据是512×512的矩阵,其包括两部分:视野外部分和视野内部分。视野内的数据才是真正的影像数据,其为影像矩阵的内切圆。视野外的数据对医生的临床诊断是没有用的。通常,这些视野外的数据被设置特殊的CT值,例如-3024,如图2中的黑色部分,以使从图像中可以明显的区别出视野内的图像。既然视野外部分不是真正的影像数据,可以先删除这部分数据以减少DICOM图像的大小。这种删除视野外部分数据可以将影像数据压缩掉21.5%即(1-π/4)。同时对删除的信息建立一个标记,如在头信息中包含了“CT值为-3024的视野外区域数据”这一信息。
步骤12)根据预设的CT值分类模板,判断视野内图像各个像素的CT值属于哪类CT值分类模板。
对于一般的CT图像,CT值的范围是从-1024到3071(总共4096),每一像素的CT值是以2字节16位的方式存储。典型的人类组织的CT值请参阅表1和图3a所示,图3a为人类组织CT值的范围。
表1
组织 | CT值范围 |
血浆(Blood plasma) | 3~14 |
血液(Blood) | 13~32 |
钙(Calcium) | 60~300 |
甲状腺(Thyroid gland) | 50~80 |
肝脏(Liver) | 45~75 |
脾(Spleen) | 35~55 |
肌肉(Muscle) | 35~50 |
胰腺(Pancreas) | 25~55 |
大脑皮层(Cerebral cortex) | 32~40 |
大脑骨髓(Cerebral medulla) | 28~32 |
脂肪(Fat) | -50~-100 |
骨骼(Bone) | 400~1000 |
空气(Air) | -1010~-990 |
为了压缩影像数据,可以将CT值范围分成4类:A[-1024~-111],B[-110~145],C[146~657],D[658~3071]。从上面的表1和CT值划分中,可以看出,最重要的CT值范围是B[-110~145]类模板,它包含了人体中大多数的组织,对于这部分,即使是一点小小的误差,都有可能影响临床诊断的结果,所以本发明中,对此部分进行无损压缩。对于C[146~657]类模板,它包含了钙和骨骼,1个CT值的误差是可以接受的。对于A[-1024~-111]类模板,它主要是空气,所以2个CT值的误差是可以接受的。对于D[658~3071]类模板,它主要包括骨骼和金属,4个CT值的误差是可以接受的。这种分法,可以保证B类模板是无损压缩,而对于C类模板的压缩,最大误差只有1个CT值,对于A和D类模板,最大误差分别是2个CT值和4个CT值。
本发明中可以通过如下的方法标注每一像素数据是属于哪个分类模板中,由于此处总共有4类,所以采用2进制可以区别表述他们,如A用00表示,B用01表示,C用10表示,D用11表示,如图3b所示,落在A范围内的所有像素数据用00表示,落在B范围内的所有像素数据用01,落在C范围内的所有像素数据用10表示,落在D范围内的所有像素数据用11表示。
步骤13)根据分类模板和各分类模板中所确定的压缩方法对落于其中的图像数据进行相应的压缩。
本实施例中应用如下的公式(1)压缩图像中各像素数据, 其中,y为压缩结果,x为各像素的CT值,intercept表示截距,即为每一分类模板中的下限值;slope表示斜率,即为每一分类模板中的压缩比率,为常数,用户可以根据自己的要求定义压缩比率。
每一分类模板所对应的intercept和slope压缩比率是不同的,intercept由该分类模板的下限值决定,如B分类模板intercept为-110,C分类模板intercept为146;slope压缩比率可以由用户根据所要求的压缩比例自己设定,但slope的应为大于或等于1的整数。
下面通过对C分类部分进行压缩来举例说明。如图4a所示,这是一个线性压缩,压缩结果满足如上的公式(1)。在此,像素的数据范围从原来的[146~675]被压缩为[0~255],原始的CT值x的范围是[146~675],Intercept是146,压缩比率是2,压缩误差是slope/2,即1个CT值。压缩结果y的范围是[0~255],正好满足了从2进制16位长度压缩到8位的长度。从而建立了应用于C分类模板的公式(1)中的intercept和slope的参数,分别为146和2。
请同时参阅图4b,通过相同的压缩方法,对其他分类部分进行压缩,并获得应用于A、B、D分类模板的公式(1)中的intercept和slope的参数:对于分类模板A[-1024~-111],intercept是-1024,slope是4;对于分类模板B[-110~145],intercept是-110,slope是1;对于分类模板D[658~3071],intercept是658,slope是9。
步骤14)得到压缩后的图像数据。
为了还原被压缩前的图像,需要对压缩的影像数据进行解码。请参阅图5,为本发明提供的对应于压缩图像进行解压的流程图,包括如下的步骤:
21)获得要解压的影像数据;
22)根据每个像素数据所在的CT值分类模板,应用与之相对应的解压算法对压缩的影像数据进行解压缩。
本实施例中对应的应用如下的公式(2)对压缩的影像数据进行解压缩
x=y·slope+intercept(2),其中x表示解压缩后的影像数据,y为通过上面的压缩方法获得的被压缩的影像数据,压缩比率slope和intercept与压缩时的值一样,由各自所处的CT值分类模板决定;y·slope表示y与slope的乘积。
步骤23)根据头信息中的记载,恢复视野外区域数据。
步骤24)将视野外区域数据加入到解压缩后的影像中,得到完整的图像。
本实施例中的对CT的DICOM图像进行压缩的方法,这种方法的压缩比例是 接近于2∶1,压缩比例大,能大幅的提高CT系统影像储存容量,其中,是通过删除视野外区域数据获得的,是通过将影像数据从16位压缩为10位(8位的影像数据和2位的分类模板数据)。通过使用CT值分类模板,可以确保重要的组织信息如本发明中所列举的B分类模板中的部分不会受到损害、丢失。对于丢失部分如空气或骨骼,其对临床诊断是没有多大影响的,另外,丢失的程度也是可以接受的。本发明的压缩方法并不影响其他压缩的方法的使用,对该CT的DICOM图像进行压缩,其他的有损、无损的压缩方法还可以使用,以提高影像的压缩比率。
在本发明的压缩方法中,客户还可以根据自己的要求,如对损害或病变的区域,不想对其进行压缩,其可以通过选择该区域将其保护起来,如图6所示,不进行本发明的压缩。
本发明的压缩方法中,基于CT值分类模板,许多现有的压缩技术可以应用,如对于B分类部分,由于该部分的数据很重要,对该部分的影像数据采用无损编码的方式压缩(如Huffman编码),而对于另外的A、C、D分类部分由于其数据对临床诊断不是很敏感,所以可以采用有损编码的方式压缩(如DCT编码)以提高压缩比率。当然,为了能够获得有效的解压图像数据,解压的解码算法应该与压缩的编码算法相对应。
通常,在CT影像中包括很多的空气数据,除了在肺部影像中有用外,基本上在其他的影像中对于临床诊断是没有用的,所以为了提高压缩比率,也可以将此部分的影像数据象视野外区域数据那样去掉,当然,应该有信息记录被删除的空气的信息,如在头信息中包含有“CT值为-1000的空气数据”,以便恢复图像数据。
应用本发明的压缩方法实现对CT DICOM图像压缩的CT装置,如图7所示,包括:扫描台30,其对受检者进行包括平扫、轴扫、螺旋扫或某一位置的X射线扫描;数据收集单元32,其对扫描台30输出的扫描数据包括由平扫、轴扫或螺旋扫或某一位置的位置扫描所获得的扫描数据进行收集、模数转换;图像重建单元33,其根据数据收集单元发送来的扫描数据包括进行图像重建,并将所重构建得到的图像数据存储到存储单元34并通过中央控制单元31在用户操作界面36上显示;供用户操作的用户操作界面36;和中央控制单元31,其连接至所述图像重建单元33的输出端,并控制扫描台30、数据收集单元32和用户操作界面36、存储单元34。该CT装置还包括图像压缩/解压计算单元35,其在中央控制单元31的控制下对由图像重建单元33重建后的DICOM图像进行压缩或对压缩后的DICOM图像进行解压缩。
所述的图像压缩/解压计算单元35包括压缩计算单元351和解压计算单元352。中央控制单元31可以自动的控制图像压缩/解压计算单元35中的压缩计算单元351对图像重建单元35重建后的DICOM图像进行压缩,也可以通过用户操作用户操作界面36向中央控制单元31发出相应的压缩指令。请再参阅图2所示,图像压缩/解压计算单元35在获得中央控制单元31发出的压缩指令后,从存储单元34中获取要压缩的DICOM图像数据,并对DICOM图像数据进行如上所述的压缩步骤:
首先,根据DICOM图像数据的特定CT值如视野外区域数据的CT值为-3024和/或空气的CT值为-1000,找出视野外区域图像数据和/或空气图像数据;压缩计算单元根据判断单元所找出的要删除的、在临床诊断中不需要的图像数据如视野外区域数据和/或空气数据进行删除,并在存储单元34中建立、存储相应的记录信息;
其次,根据CT值分类模板和DICOM图像数据中视野内图像的各像素CT值,判断各像素CT值属于哪个CT值分类模版,然后根据各分类模板,运用与之相对应的压缩算法对落入其中的像素数据进行压缩,如应用上述的公式(1)对各像素数据进行压缩;
最后得到压缩后的图像数据,并存储于存储单元34中。
其中,所述的CT值分类模板可如上述的压缩方法中的步骤12中所述,分成A、B、C、D四类,也可以根据客户自己的要求,建立符合自身要求的CT值分类模板,并对每一分类模板设置如上述压缩方法中的公式(1)的压缩编码方法,也可以是其他的如上所述的无损或有损压缩编码方法。
为了显示准确完整的影像,需要对压缩的图像数据进行解压缩。中央控制单元31控制图像压缩/解压计算单元35对压缩图像进行解压缩。图像压缩/解压计算单元35中的解压计算单元352在中央控制单元31的控制下,获取存储单元34中需要解压缩的图像数据,根据该图像数据所属的各CT值分类模板,应用与压缩编码所对应的解压方法对各个像素数据进行解压解码,如当压缩编码采用的是如公式(1)所述的压缩编码方法,则解压缩时,采用与之对应的公式(2)的解压方法,恢复图像数据;同时,该解压计算单元352恢复被删除的图像数据如视野外区域数据和空气数据,其通过将存储单元34中保存的视野外区域数据添到恢复图像视野外区域。
所述的CT值分类模板和每一分类模板的压缩编码算法及与之对应的解压缩算法预存于存储单元34中。
请再参阅图6,用户可以通过用户操作界面36对图像进行操作,对不需要进行压缩、重要的部位进行选中保护起来,则中央控制单元31控制图像压缩/解压计算单元35对该保护起来的图像数据不进行压缩。
为了提高压缩比率,用户可以根据自己的需要,对那些对临床诊断没有用的图像数据,可以象视野外区域数据那样去掉,如在CT影像中包括很多的空气数据,除了在肺部影像中有用外,基本上在其他的影像中对于临床诊断是没有用的,所以,该部分的影像数据可以删掉。CT装置可以预设的各物质、组织的CT值来做出判断,如空气的CT值为-1000,如当用户选中用户操作界面36上的“删除空气数据”的操作键时,则图像压缩/解压计算单元35中的压缩计算单元对DICOM影像中的空气数据进行删除,同时在存储单元34建立一条相应的记录信息。在解压恢复图像数据时根据相应的记录信息,对空气数据进行恢复。
Claims (20)
1、一种对CT重建后的图像进行压缩的方法,包括如下的步骤:
步骤10)获得要压缩的重建后的图像;所述的重建后的图像包括视野外区域和视野内区域;
步骤11)根据预设的CT值分类模板,判断图像中的视野内区域各个像素的CT值属于哪类CT值分类模板;
步骤12)根据CT值分类模板和各分类模板中所确定的压缩方法对落于其中的像素数据进行相应的压缩。
2、如权利要求1所述的对CT重建后的图像进行压缩的方法,其特征在于:还包括如下的步骤:除去视野外区域数据,并建立相应的记录信息。
3、如权利要求2所述的对CT重建后的图像进行压缩的方法,其特征在于:所述的重建后的图像为医学图像文件格式(Digital Imagingand Communications in Medicine,简称DICOM)图像,其为512×512的矩阵。
4、如权利要求3所述的对CT重建后的图像进行压缩的方法,其特征在于:所述的预设的CT值分类模板包括A[-1024~-111],B[-110~145],C[146~657],D[658~3071]四类;对CT值分类模板B[-110~145]内的图像数据采用无损压缩,对CT值分类模板A[-1024~-111],C[146~657],D[658~3071]的图像数据采用有损压缩。
5、如权利要求4所述的对CT重建后的图像进行压缩的方法,其特征在于:A、B、C、D四类的CT值分类模板中的所对应的压缩方式是按照如下的公式进行: 其中,y为压缩结果,x为各像素的CT值,intercept为每一分类模板中的下限值;slope为每一分类模板中的压缩比率,为常数,用户可以根据自己的要求定义压缩比率;其中,对应于CT值分类模板B[-110~145]的压缩比率slope为1。
6、如权利要求5所述的对CT重建后的图像进行压缩的方法,其特征在于:对应于CT值分类模板A[-1024~-111],C[146~657],D[658~3071]压缩比率slope分别为4,2,9。
7、如权利要求4-6中任一项所述的对CT重建后的图像进行压缩的方法,其特征在于:A、B、C、D四类CT值分类模板分别采用二进制表示,A用00,B用01,C用10,D用11表示。
8、如权利要求7所述的对CT重建后的图像进行压缩的方法,其特征在于:视野外区域像素的CT值为-3024。
9、对权利要求8压缩后的CT重建图像进行解压缩的方法,包括如下的步骤:
21)获得要解压的影像数据;
22)根据每个像素数据所在的CT值分类模板,应用与之相对应的解压算法对压缩的影像数据进行解压缩;
23)恢复视野外区域的图像数据;
24)将视野外区域数据加入到解压缩后的影像中,得到完整的图像。
10、如权利要求9所述的对压缩后的CT重建图像进行解压缩的方法,其特征在于应用如下的公式对压缩的影像数据进行解压缩x=y·slope+intercept,其中x表示解压缩后的影像数据,y为被压缩的影像数据,压缩比率slope和intercept由各自所处的CT值分类模板决定。
11、一种CT装置,包括:
扫描台,其对受检者进行X射线扫描;
数据收集单元,其对扫描台输出的扫描数据进行收集、模数转换;
图像重建单元,其根据数据收集单元发送来的扫描数据进行图像重建,并将所重构建得到的图像数据存储到存储单元和/或通过中央控制单元在用户操作界面上显示;
供用户操作的用户操作界面;
用于存储数据和信息的存储单元;
和中央控制单元,其连接至所述图像重建单元的输出端,并控制扫描台、数据收集单元和用户操作界面、存储单元;
其特征在于,所述的存储单元中存储有CT值分类模板及应用于CT值分类模板的压缩算法和解压缩算法;该CT装置还包括图像压缩/解压计算单元,其在中央控制单元的控制下调用存储在存储单元的CT值分类模板及相应的压缩算法和解压缩算法,对重建后的图像进行压缩或对压缩后的图像进行解压缩。
12、如权利要求11所述的CT装置,其特征在于所述的图像压缩/解压计算单元包括压缩计算单元和解压计算单元;中央控制单元可以自动的控制图像压缩/解压计算单元中的压缩计算单元对图像重建单元重建后的图像进行压缩,也可以通过用户操作用户操作界面向中央控制单元发出相应的压缩指令。
13、如权利要求12所述的CT装置,其特征在于所述的所述重建单元重建后的图像为医学图像文件格式(Digital Imaging andCommunications in Medicine,简称DICOM)图像,其为512×512的矩阵,包括视野内区域和视野外区域。
14、如权利要求13所述的CT装置,其特征在于:所述的CT值分类模板包括A[-1024~-111],B[-110~145],C[146~657],D[658~3071]四类;对CT值分类模板B[-110~145]内的图像数据采用无损压缩,对CT值分类模板A[-1024~-111],C[146~657],D[658~3071]的图像数据采用有损压缩。
15、如权利要求14所述的CT装置,其特征在于:A、B、C、D四类的CT值分类模板中的所对应的压缩方式是按照如下的公式进行: 其中,y为压缩结果,x为各像素的CT值,intercept为每一分类模板中的下限值;slope为每一分类模板中的压缩比率,为常数,用户可以根据自己的要求定义压缩比率;其中,对应于CT值分类模板B[-110~145]的压缩比率slope为1;与之相对应的对压缩的影像数据进行解压缩的方法为x=y·slope+intercept,其中x表示解压缩后的影像数据,y为被压缩的影像数据,压缩比率slope和intercept由各自所处的CT值分类模板决定。
16、如权利要求15所述的CT装置,其特征在于:对应于CT值分类模板A[-1024~-111],C[146~657],D[658~3071]的压缩比率slope分别为4,2,9。
17、如权利要求14-16中任一项所述的CT装置,其特征在于:A、B、C、D四类CT值分类模板分别采用二进制表示,A用00、B用01、C用10、D用11表示。
18、如权利要求13所述的CT装置,其特征在于:当压缩计算单元收到中央控制单元发送来的对重建图像进行压缩的指令时,从存储单元中调用要压缩的重建图像数据,并将图像中视野外区域数据删除,判断视野内的图像数据的各个像素CT值属于那个CT值分类模板,并根据各像素所落入CT值分类模板的压缩算法对像素数据进行压缩;所述的用户操作界面上设置有“删除空气数据”的操作键。
19、如权利要求18所述的CT装置,其特征在于:当解压计算单元收到中央控制单元发送来的对压缩后的重建图像进行解压缩的指令时,从存储单元中调用要解压缩的图像数据,根据各像素所落入CT值分类模板的解压缩算法对压缩的像素数据进行解压缩;并恢复视野外区域的图像数据;将视野外区域数据加入到解压缩后的影像中,得到完整的图像。
20、如权利要求19所述的CT装置,其特征在于:当压缩计算单元删除视野外区域数据时,其同时在存储单元中建立、存储一条对应的记录信息;压缩计算单元压缩后的图像数据或解压计算单元解压后的图像数据被存储到存储单元中。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200810137772A CN101627918A (zh) | 2008-07-18 | 2008-07-18 | Ct影像压缩方法和装置 |
JP2009160066A JP2010022821A (ja) | 2008-07-18 | 2009-07-06 | Ct画像の圧縮方法及びその装置 |
US12/502,867 US20100014627A1 (en) | 2008-07-18 | 2009-07-14 | Method and apparatus for ct image compression |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200810137772A CN101627918A (zh) | 2008-07-18 | 2008-07-18 | Ct影像压缩方法和装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101627918A true CN101627918A (zh) | 2010-01-20 |
Family
ID=41530285
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200810137772A Pending CN101627918A (zh) | 2008-07-18 | 2008-07-18 | Ct影像压缩方法和装置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20100014627A1 (zh) |
JP (1) | JP2010022821A (zh) |
CN (1) | CN101627918A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102655592A (zh) * | 2011-03-04 | 2012-09-05 | 索尼公司 | 信息处理装置、信息处理方法和程序 |
CN103999087A (zh) * | 2011-12-15 | 2014-08-20 | 皇家飞利浦有限公司 | 针对接收方优化的医学成像重建 |
CN108572183A (zh) * | 2017-03-08 | 2018-09-25 | 清华大学 | 检查设备和分割车辆图像的方法 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2663074B1 (en) * | 2012-05-11 | 2016-04-13 | XRing Technologies GmbH | Method and system for lossy compression and decompression of computed tomography data |
EP3404919A1 (en) * | 2017-05-17 | 2018-11-21 | Koninklijke Philips N.V. | Vector-valued diagnostic image encoding |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1512853A (zh) * | 2001-06-01 | 2004-07-14 | ��ʽ��������ҽҩ | X射线ct装置、图像处理装置以及x射线ct装置的图像处理方法 |
CN1934587A (zh) * | 2004-03-22 | 2007-03-21 | 康泰克斯特视觉公司 | 增强计算机化x射线断层摄影图像的方法、计算机程序产品以及设备 |
US20070065032A1 (en) * | 2005-09-22 | 2007-03-22 | Hernandez Albert A | Method and apparatus for boundary-based image compression |
Family Cites Families (42)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6275773A (ja) * | 1985-09-30 | 1987-04-07 | Toshiba Corp | 画像圧縮装置 |
US4788706A (en) * | 1987-12-17 | 1988-11-29 | General Electric Company | Method of measurement of x-ray energy |
JPH01270175A (ja) * | 1988-04-21 | 1989-10-27 | Toshiba Corp | 画像処理装置 |
EP0381785B1 (de) * | 1989-02-06 | 1993-11-18 | Siemens Aktiengesellschaft | Rotierende Datenübertragungsvorrichtung |
US5173787A (en) * | 1990-03-07 | 1992-12-22 | Fuji Photo Film Co., Ltd. | Method for storing image signals |
US5703965A (en) * | 1992-06-05 | 1997-12-30 | The Regents Of The University Of California | Image compression/decompression based on mathematical transform, reduction/expansion, and image sharpening |
US5594767A (en) * | 1995-11-02 | 1997-01-14 | General Electric Company | Methods and apparatus for enhancing image sharpness |
US5825909A (en) * | 1996-02-29 | 1998-10-20 | Eastman Kodak Company | Automated method and system for image segmentation in digital radiographic images |
IL119283A0 (en) * | 1996-09-19 | 1996-12-05 | Elscint Ltd | Adaptive filtering |
JPH1127680A (ja) * | 1997-07-01 | 1999-01-29 | Access:Kk | データ処理システム |
JP4035216B2 (ja) * | 1997-12-25 | 2008-01-16 | キヤノン株式会社 | 画像処理方法及び画像処理装置 |
US6035012A (en) * | 1998-05-14 | 2000-03-07 | Gen Electric | Artifact correction for highly attenuating objects |
US6198850B1 (en) * | 1998-06-12 | 2001-03-06 | Xerox Corporation | System and method for segmentation dependent lossy and lossless compression for higher quality |
US6301324B1 (en) * | 1999-03-31 | 2001-10-09 | General Electric Company | RF slipring receiver for a computerized tomography system |
US6327327B1 (en) * | 1999-09-27 | 2001-12-04 | Picker International, Inc. | Multi-channel segmented slip ring |
JP2001144966A (ja) * | 1999-11-11 | 2001-05-25 | Canon Inc | 画像処理装置及びその方法 |
US6795583B1 (en) * | 1999-11-24 | 2004-09-21 | General Electric Company | Image data compression employing embedded compression and decompression codes |
US7421136B2 (en) * | 1999-11-24 | 2008-09-02 | Ge Medical Systems Information Technologies Inc. | Image tessellation for region-specific coefficient access |
JP2002123816A (ja) * | 2000-10-18 | 2002-04-26 | Fuji Photo Film Co Ltd | 断層画像処理方法および装置並びに記録媒体 |
US6754376B1 (en) * | 2000-11-22 | 2004-06-22 | General Electric Company | Method for automatic segmentation of medical images |
US6501819B2 (en) * | 2000-12-18 | 2002-12-31 | Ge Medical Systems Global Technology Company, Llc | Medical diagnostic method and apparatus to control dual energy exposure techniques based on image information |
AU2002243602A1 (en) * | 2001-01-19 | 2002-07-30 | U.S. Army Medical Research And Materiel Command | A method and apparatus for generating two-dimensional images of cervical tissue from three-dimensional hyperspectral cubes |
US6480560B2 (en) * | 2001-03-16 | 2002-11-12 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for motion gating using CT projection data |
JP3766608B2 (ja) * | 2001-05-02 | 2006-04-12 | テラリコン・インコーポレイテッド | ネットワーク環境における三次元画像表示装置 |
IL152520A0 (en) * | 2001-12-19 | 2003-12-23 | Philips Medical Systems Techno | Wireless data transmission in ct-scanners |
JP2003290216A (ja) * | 2002-03-27 | 2003-10-14 | Ge Medical Systems Global Technology Co Llc | X線ct装置、システムおよび投影データ圧縮復元方法 |
US20040022447A1 (en) * | 2002-07-31 | 2004-02-05 | General Electric Company | Method and system for image compression and decompression using span of interest of an imaging sequence |
US6999549B2 (en) * | 2002-11-27 | 2006-02-14 | Ge Medical Systems Global Technology, Llc | Method and apparatus for quantifying tissue fat content |
US6944333B2 (en) * | 2003-04-30 | 2005-09-13 | Ventana Medical Systems, Inc. | Color image compression via spectral decorrelation and elimination of spatial redundancy |
DE10322138B3 (de) * | 2003-05-16 | 2004-09-30 | Siemens Ag | Datenübertragungssystem mit integriertem Fehlerüberwachungs- und Diagnosewerkzeug für die Datenübertragung in einem Computertomographen sowie zugehöriges Verfahren |
JP2005185788A (ja) * | 2003-12-26 | 2005-07-14 | Shimadzu Corp | コンピュータ画像診断装置 |
US7684629B2 (en) * | 2004-01-26 | 2010-03-23 | Fujifilm Corporation | Data compression apparatus, and data compression program storage medium |
US7664184B2 (en) * | 2004-07-21 | 2010-02-16 | Amimon Ltd. | Interpolation image compression |
JP4533043B2 (ja) * | 2004-08-25 | 2010-08-25 | キヤノン株式会社 | 画像符号化装置及び方法、並びに、コンピュータプログラム及びコンピュータ可読記憶媒体 |
US7826670B2 (en) * | 2005-06-15 | 2010-11-02 | Fujifilm Corporation | Data compression apparatus and data compression program storage medium |
JP2007097977A (ja) * | 2005-10-07 | 2007-04-19 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
DE602007000299D1 (de) * | 2006-01-18 | 2009-01-15 | Toshiba Kk | Röntgen-CT-Gerät und Verknüpfungssystem zur Kommunikation von medizinischen Daten |
JP4697967B2 (ja) * | 2006-03-02 | 2011-06-08 | キヤノン株式会社 | 画像符号化装置及びその制御方法 |
JP2008079760A (ja) * | 2006-09-27 | 2008-04-10 | Fujifilm Corp | 画像圧縮処理方法及び画像圧縮処理装置並びに医用ネットワークシステム |
JP2008148970A (ja) * | 2006-12-19 | 2008-07-03 | Ge Medical Systems Global Technology Co Llc | X線ctシステム |
US8121417B2 (en) * | 2007-03-19 | 2012-02-21 | General Electric Company | Processing of content-based compressed images |
DE102007018324B3 (de) * | 2007-04-18 | 2008-06-05 | Siemens Ag | Bilddatenerfassungssystem eines Röntgen-, CT- oder MRT-Geräts mit integriertem Datenkomprimierungsmodul zur Datenreduktion akquirierter Bilddaten |
-
2008
- 2008-07-18 CN CN200810137772A patent/CN101627918A/zh active Pending
-
2009
- 2009-07-06 JP JP2009160066A patent/JP2010022821A/ja active Pending
- 2009-07-14 US US12/502,867 patent/US20100014627A1/en not_active Abandoned
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1512853A (zh) * | 2001-06-01 | 2004-07-14 | ��ʽ��������ҽҩ | X射线ct装置、图像处理装置以及x射线ct装置的图像处理方法 |
CN1934587A (zh) * | 2004-03-22 | 2007-03-21 | 康泰克斯特视觉公司 | 增强计算机化x射线断层摄影图像的方法、计算机程序产品以及设备 |
US20070065032A1 (en) * | 2005-09-22 | 2007-03-22 | Hernandez Albert A | Method and apparatus for boundary-based image compression |
Non-Patent Citations (1)
Title |
---|
蒋宏晖 等: "DICOM格式CT图像的数据库存储技术", 《计算机应用》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102655592A (zh) * | 2011-03-04 | 2012-09-05 | 索尼公司 | 信息处理装置、信息处理方法和程序 |
CN103999087A (zh) * | 2011-12-15 | 2014-08-20 | 皇家飞利浦有限公司 | 针对接收方优化的医学成像重建 |
CN103999087B (zh) * | 2011-12-15 | 2020-06-30 | 皇家飞利浦有限公司 | 针对接收方优化的医学成像重建 |
CN108572183A (zh) * | 2017-03-08 | 2018-09-25 | 清华大学 | 检查设备和分割车辆图像的方法 |
CN108572183B (zh) * | 2017-03-08 | 2021-11-30 | 清华大学 | 检查设备和分割车辆图像的方法 |
Also Published As
Publication number | Publication date |
---|---|
JP2010022821A (ja) | 2010-02-04 |
US20100014627A1 (en) | 2010-01-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wong et al. | Radiologic image compression-a review | |
Kivijärvi et al. | A comparison of lossless compression methods for medical images | |
JP5394598B2 (ja) | 医用画像圧縮装置、医用画像圧縮方法、及び予測知識データベース作成装置 | |
Karadimitriou et al. | The centroid method for compressing sets of similar images | |
CN101627918A (zh) | Ct影像压缩方法和装置 | |
Karadimitriou | Set redundancy, the enhanced compression model, and methods for compressing sets of similar images | |
Gowda et al. | A novel method of data compression using ROI for biomedical 2D images | |
Karadimitriou et al. | Min-max compression methods for medical image databases | |
US6832001B1 (en) | Image processing apparatus and method | |
Thanki et al. | Hybrid and advanced compression techniques for medical images | |
Kaur et al. | Image compression using HAAR wavelet transform and discrete cosine transform | |
JP2010022821A6 (ja) | Ct画像の圧縮方法及びその装置 | |
Špelič et al. | Lossless compression of threshold-segmented medical images | |
Ho et al. | Dynamic image data compression in spatial and temporal domains: Theory and algorithm | |
Kil et al. | Lossless medical image compression using redundancy analysis | |
Aldemir et al. | Chain code strategy for lossless storage and transfer of segmented binary medical data | |
CN113538253B (zh) | 一种医疗影像的处理方法、介质及ct装置 | |
Komma et al. | Lossless Volume Data Compression Schemes. | |
Menegaz | Trends in medical image compression | |
JP2001218062A (ja) | 画像処理装置及び方法と記憶媒体 | |
Naaz et al. | Implementation of hybrid algorithm for image compression and decompression | |
Kumar et al. | Lossless compression of digital mammography using fixed block segmentation and pixel grouping | |
Akhtar et al. | Significance of ROI coding using MAXSHIFT scaling applied on MRI images in teleradiology-telemedicine | |
Ali et al. | Significance of region of interest applied on MRI and CT images in teleradiology-telemedicine | |
Nasifoglu et al. | Multi-regional Adaptive Image Compression (AIC) for hip fractures in pelvis radiography |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20100120 |