《计算机断层图像的重建方法.pdf》由会员分享,可在线阅读,更多相关《计算机断层图像的重建方法.pdf(7页珍藏版)》请在专利查询网上搜索。
1、10申请公布号CN104207798A43申请公布日20141217CN104207798A21申请号201310365301X22申请日20130820201310211930720130531CNA61B6/03200601G06T11/0020060171申请人上海联影医疗科技有限公司地址201815上海市嘉定区嘉定工业区兴贤路1180号8幢72发明人张笛儿54发明名称计算机断层图像的重建方法57摘要一种计算机断层图像的重建方法,包括对待扫描对象进行CT扫描,并通过数据采集系统采集原始数据,且将其发送到数据测量系统;通过所述数据测量系统利用分段对数算法对所述原始数据处理,并获得返回值AD。
2、TC;将所述返回值ADTC通过射频端口发送至重建服务器进行重建图像。本发明提供的计算机断层图像的重建方法,利用分段对数算法对所述原始数据进行处理,保留了更多小数位数,从而可以获得能够产生高质量图像的高精确度的数据,解决了现有技术中由于受数据测量系统射频端口的限制,计算机断层扫描系统不能获取更高质量图像的问题。66本国优先权数据51INTCL权利要求书1页说明书4页附图1页19中华人民共和国国家知识产权局12发明专利申请权利要求书1页说明书4页附图1页10申请公布号CN104207798ACN104207798A1/1页21一种计算机断层图像的重建方法,其特征在于,包括如下步骤对待扫描对象进行C。
3、T扫描,并通过数据采集系统采集原始数据,且将其发送到数据测量系统;通过所述数据测量系统利用分段对数算法对所述原始数据处理,并获得返回值ADTC;将所述返回值ADTC通过射频端口发送至重建服务器进行图像重建。2如权利要求1所述的计算机断层图像的重建方法,其特征在于,所述原始数据包括计算机断层扫描设备中检测器上的暗电流IDTCDARK、所述检测器返回的积分强度值IDTC以及积分时间T。3如权利要求2所述的计算机断层图像的重建方法,其特征在于,所述返回值ADTC利用如下公式获得ADTCA0CLOG2IDTCIDTC_DARKCLOG2T其中CLOG2IDTCIDTCDARK以及CLOG2T为对数运算。
4、单元,C是位移运算系数,A0为平移量系数。4如权利要求3所述的计算机断层图像的重建方法,其特征在于,所述数据采集系统采集的原始数据位宽为24BIT,经过所述数据测量系统处理之后的数据位宽为16BIT。5如权利要求4所述的计算机断层图像的重建方法,其特征在于,所述位移运算系数C212。6如权利要求5所述的计算机断层图像的重建方法,其特征在于,所述对数运算单元CLOG2X采用如下分段算法进行运算7如权利要求6所述的计算机断层图像的重建方法,其特征在于,所述对数运算单元CLOG2X采用如下分段算法进行查表运算其中M为输出的整数位,U为查找对数表输入的索引值,为查找对数表输入的余数值。8如权利要求3所。
5、述的计算机断层图像的重建方法,其特征在于,所述返回值ADTC0,2161。9如权利要求8所述的计算机断层图像的重建方法,其特征在于,所述平移量系数A0利用所述公式A02161212LOG2TMAX获取,其中所述TMAX是单次采样中积分时间的最大值。权利要求书CN104207798A1/4页3计算机断层图像的重建方法技术领域0001本发明涉及计算机断层扫描领域,特别涉及一种计算机断层图像的重建方法。背景技术0002计算机断层扫描设备简称CT在进行临床扫描时,数据测量系统DMS通过射频端口传送CT扫描数据。为了节约带宽,一般CT射频端口的带宽只有16BIT,而典型的原始数据为24BIT,因此需要对。
6、原始数据进行压缩,通常采用对数压缩算法对数据进行处理。已有的对数压缩算法是对数据进行LOG2运算以2为底的对数运算以及位移运算处理。LOG2运算是一种有损压缩算法,会降低数据精度,影响图像质量。图像质量与图像损失的程度成负相关,图像损失的程度取决于保留的小数位位数,小数位位数越多,则数据越精确,对图像质量的影响越小。典型的24BIT的原始数据,经过LOG2压缩处理后转为16BIT,16BIT其中包括5BIT整数位、11BIT小数位,如果需要让压缩后的数据全部落入16BIT范围之内,则只能保留11位小数。发明内容0003基于此,有必要提供一种能够获得更高质量的计算机断层图像的重建方法。0004一。
7、种计算机断层图像的重建方法,包括如下步骤0005对待扫描对象进行CT扫描,并通过数据采集系统采集原始数据,且将其发送到数据测量系统;0006通过所述数据测量系统利用分段对数算法对所述原始数据处理,并获得返回值ADTC;0007将所述返回值ADTC通过射频端口发送至重建服务器进行图像重建。0008在一个实施例中,所述原始数据包括计算机断层扫描设备中的检测器上的暗电流IDTCDARK、所述检测器返回的积分强度值IDTC以及积分时间T。0009在一个实施例中,所述返回值ADTC利用如下公式获得0010ADTCA0CLOG2IDTCIDTC_DARKCLOG2T0011其中CLOG2IDTCIDTCD。
8、ARK以及CLOG2T为对数运算单元,C是位移运算系数,A0为平移量系数。0012在一个实施例中,所述数据采集系统采集的原始数据位宽为24BIT,经过所述数据测量系统处理之后的数据位宽为16BIT。0013在一个实施例中,所述位移运算系数C212。0014在一个实施例中,所述对数运算单元CLOG2X采用如下分段算法进行运算00150016在一个实施例中,所述对数运算单元CLOG2X采用如下分段算法进行查表运算说明书CN104207798A2/4页400170018其中M为输出的整数位,U为查找对数表输入的索引值,为查找对数表输入的余数值。0019在一个实施例中,所述返回值ADTC0,2161。。
9、0020在一个实施例中,所述平移量系数A0利用所述公式A02161212LOG2TMAX获取,其中所述TMAX是单次采样中积分时间的最大值。0021本发明提供的计算机断层图像的重建方法,利用分段对数算法对所述原始数据进行处理,保留了更多小数位数,从而可以获得能够产生高质量图像的高精确度数据,解决了现有技术中由于受数据测量系统射频端口限制,不能获取更高质量计算机断层图像的问题。附图说明0022图1是本发明中一种计算机断层图像的重建方法流程示意图;0023图2是采用分段对数运算算法与常用对数运算算法的值域范围的对比示意图。具体实施方式0024为使本发明的上述目的、特征和优点能够更为明显易懂,下面结。
10、合附图和实施例对本发明的具体实施方式做详细的说明。0025在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是本发明还可以采用其他不同于在此描述的其它方式来实施,因此本发明不受下面公开的具体实施例的限制。0026正如背景技术中所描述的,现有的对数压缩算法是对数据进行LOG2运算以及位移运算处理。LOG2运算是一种有损压缩算法,可降低数据精度,影响图像质量。图像质量与图像损失的程度成负相关,图像损失的程度取决于保留的小数位的位数,小数位位数越多,则数据越精确,对图像质量的影响越小。典型的原始数据为24BIT,经过LOG2压缩处理后转为16BIT,16BIT其中包括5BIT整数位、11BIT。
11、小数位,如果需要让压缩后的数据全部落入16BIT范围之内,则只能保留11位小数。为了得到更高质量的图像,有必要提供一种能够保留更多小数位数的图像重建方法。0027一种计算机断层图像的重建方法,包括如下步骤0028步骤S1对待扫描对象进行CT扫描,并通过数据采集系统采集原始数据,且将其发送至数据测量系统。0029数据采集系统通过多组数字积分电路,将每个探测器探测到的光电信号转换为24BIT位宽的逻辑信号,每个探测器对应一组数字积分电路。所述全部数字积分电路返回的逻辑信号通过特定的时序,汇总到数据测量系统的高速数据链路上。0030步骤S2通过所述数据测量系统利用分段对数算法对所述原始数据进行处理并。
12、获得返回值ADTC。0031具体地,所述逻辑信号进入所述数据测量系统的高速数据链路,在数据链路上,所说明书CN104207798A3/4页5述数据测量系统利用分段对数算法对所述数据进行同步处理并得到返回值。0032步骤S3将所述返回值ADTC通过射频端口发送至重建服务器进行图像重建。0033下面结合附图与具体实施方式对本发明的技术方案作进一步说明。0034如图1所示执行步骤S1对待扫描对象进行CT扫描,并通过数据采集系统采集原始数据,且将其发送到数据测量系统。0035所述数据采集系统返回的原始数据包括计算机断层扫描系统中的检测器上的暗电流IDTCDARK,所述检测器返回的积分强度值IDTC以及。
13、积分时间T。上述返回的原始数据的位宽为24BIT,射频端口的带宽无法满足传输24BIT位宽的数据,因此需要经过数据测量系统进行后处理。0036执行步骤S2通过所述数据测量系统利用分段对数算法对所述原始数据进行处理,并获得返回值ADTC。由于受射频端口带宽的影响,所述返回值ADTC只能是16BIT的数据,所述返回值ADTC利用如下公式获得0037ADTCA0CLOG2IDTCIDTC_DARKCLOG2T0038其中A0为平移量系数,要保证ADTC0,2161,A0可以通过如下公式获取,A02161212LOG2TMAX,其中TMAX是单次采样中积分时间的最大值。0039其中C是位移运算系数,当。
14、C2A时,即表示左移A位,保留A位小数精度,在本实施例中,优选地C212,保留12位小数精度。0040其中CLOG2IDTCIDTCDARK以及CLOG2T为对数运算单元,在本实施例中,所述分段对数运算单元CSEGLOG2X的算法如下00410042如上述算法,保留了12为的小数精度,即当输入数据的有效位数超过包含12BIT时,采用212LOG2X运算,并左移12位;当输入数据的有效位数不足12BIT时,直接加上11212后返回。0043反之数据所需要的最小带宽是由其最大值减去最小值决定的,当X1,224时,采用原有算法其最大值减去最小值为MAX212LOG2XMIN212LOG2X15216。
15、,超出了16BIT可表示的范围,而采用了本发明中技术方案中改进后的算法,MAX212LOG2XMIN11212X08125216,在16BIT可表示的范围之内。0044如图2所示,采用改进后的算法值域范围实线为0,6104,明显小于原有的算法值域范围虚线0,10104,从而达到更好的压缩效果。0045在另一个实施例中,所述分段对数运算单元CLOG2X采用如下公式所示的查表算法,所述公式如下00460047说明书CN104207798A4/4页60048表格10049如表格1所示,其中X是输入值;如表格1第1列所示,根据X的最高非零位判断X是否大于等于212;表格1第2行显示的是X小于212时的。
16、情形,第3到14行显示的是X大于等于212时的情形;如表格1第4列所示,M是输出的整数位。0050参数U为查找对数表输入的索引值,如表格1第3列所示,U1,XXXXXXXXXX,其中XXXXXXXXXX是X在最高位之后连续10个二进制位,每个X表示0或1;LOG2U210和LOG2U的值通过查表取得,该数值表包含210个值,对应每一个LOG21XXXXXXXXXX的值。0051参数为查找对数表输入的余数值,如表格1第3列所示0YYY,每个Y表示0或1。0052执行步骤S3将所述返回值ADTC通过射频端口发送至重建服务器进行重建图像。0053如上所述,本发明提供的计算机断层图像的重建方法,利用分段对数算法对所述数据进行处理,保留了更多小数位数,从而可以获得能够产生高质量图像的高精确度数据,解决了现有技术中由于受数据测量系统射频端口限制,计算机断层扫描系统不能获取更高质量图像的问题。0054本发明虽然已以较佳实施例公开如上,但其并不是用来限定本发明,任何本领域技术人员在不脱离本发明的精神和范围内,都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改,因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰,均属于本发明技术方案的保护范围。说明书CN104207798A1/1页7图1图2说明书附图CN104207798A。