CN109993691A - 基于分段式线性转换的非均匀图像转换方法及系统 - Google Patents
基于分段式线性转换的非均匀图像转换方法及系统 Download PDFInfo
- Publication number
- CN109993691A CN109993691A CN201910261626.0A CN201910261626A CN109993691A CN 109993691 A CN109993691 A CN 109993691A CN 201910261626 A CN201910261626 A CN 201910261626A CN 109993691 A CN109993691 A CN 109993691A
- Authority
- CN
- China
- Prior art keywords
- gray value
- uniform
- bitmap
- black region
- region
- 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
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 101
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000011426 transformation method Methods 0.000 claims abstract description 9
- 238000004590 computer program Methods 0.000 claims description 13
- CXENHBSYCFFKJS-OXYODPPFSA-N (Z,E)-alpha-farnesene Chemical compound CC(C)=CCC\C(C)=C\C\C=C(\C)C=C CXENHBSYCFFKJS-OXYODPPFSA-N 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000013507 mapping Methods 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 238000012545 processing Methods 0.000 description 6
- 230000009466 transformation Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 3
- 230000011218 segmentation Effects 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 241001270131 Agaricus moelleri Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008439 repair process Effects 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/04—Context-preserving transformations, e.g. by using an importance map
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本公开提供了一种基于分段式线性转换的非均匀图像转换方法及系统。其中,基于分段式线性转换的非均匀图像转换方法,包括:计算灰度值不均匀的非8位图的所有像素点的灰度值;按照像素点的灰度值大小,将灰度值不均匀的非8位图分成黑色区域和非黑色区域,同时将非黑色区域划分成若干区间区域;利用分段线性转换方法,计算黑色区域以及非黑色区域内各个区间区域转换为相应灰度值均匀的8位图的灰度转换曲线;利用各个线性分段的灰度转换曲线,得到灰度值均匀的8位图。本公开能够利用分段式线性转换,将灰度值不均匀的非8位图准确地转换为均匀的8位图像。
Description
技术领域
本公开属于图像处理领域,尤其涉及一种基于分段式线性转换的非均匀图像转换方法及系统。
背景技术
本部分的陈述仅仅是提供了与本公开相关的背景技术信息,不必然构成在先技术。
CT是用X射线束对人体某部一定厚度的层面进行扫描,扫描所得信息经计算而获得每个体素的X射线衰减系数或吸收系数,再排列成矩阵,即数字矩阵。经数字/模拟转换器把数字矩阵中的每个数字转为由黑到白不等灰度的小方块,即像素,并按矩阵排列,即构成CT图像。所以,CT图像是重建图像。每个体素的X射线吸收系数可以通过不同的数学方法算出。
目前,CT图像的存储方式为将其转换为特定位数的图像数据,比如:16位图像数据。以16位图像数据为例:按照CT图像成像规律,大部分CT图像都无法完美的分布到十六位数据(0~65535)上,根据现有的CT图像统计,大部分图像数字矩阵数据处在0~4000之间,而且大部分数据集中在0~2000上,在2000~4000上,像素点占很小的一部分。这样就带来极大的问题是:图像灰度不均匀。
发明人发现,在机器视觉领域,各种算法对8位图有更好的支持,对于不均匀的图像转换为8位数图像数据时,按照最值线性转换,对于低灰度值和高灰度值及其系数,转换后的图像仍然是不均匀的。例如:将CT图像的16位图像数据转换为8位图像数据,按照通用方法转换图像,CT图像灰度值大部分将分布到0~16之间,相对于0~255,这个数据是极其不均匀的。
发明内容
为了解决上述问题,本公开的第一个方面提供一种基于分段式线性转换的非均匀图像转换方法,其能够利用分段式线性转换,将灰度值不均匀的非8位图准确地转换为均匀的8位图像。
为了实现上述目的,本公开采用如下技术方案:
一种基于分段式线性转换的非均匀图像转换方法,包括:
计算灰度值不均匀的非8位图的所有像素点的灰度值;
按照像素点的灰度值大小,将灰度值不均匀的非8位图分成黑色区域和非黑色区域,同时将非黑色区域划分成若干区间区域;
利用分段线性转换方法,计算黑色区域以及非黑色区域内各个区间区域转换为相应灰度值均匀的8位图的灰度转换曲线;
利用各个线性分段的灰度转换曲线,得到灰度值均匀的8位图。
为了解决上述问题,本公开的第二个方面提供一种基于分段式线性转换的非均匀图像转换系统,其能够利用分段式线性转换,将灰度值不均匀的非8位图准确地转换为均匀的8位图像。
为了实现上述目的,本公开采用如下技术方案:
一种基于分段式线性转换的非均匀图像转换系统,包括:
灰度值计算模块,其用于计算灰度值不均匀的非8位图的所有像素点的灰度值;
图像划分模块,其用于按照像素点的灰度值大小,将灰度值不均匀的非8位图分成黑色区域和非黑色区域,同时将非黑色区域划分成若干区间区域;
灰度转换曲线计算模块,其用于利用分段线性转换方法,计算黑色区域以及非黑色区域内各个区间区域转换为相应灰度值均匀的8位图的灰度转换曲线;
图像转换模块,其用于利用各个线性分段的灰度转换曲线,得到灰度值均匀的8位图。
为了解决上述问题,本公开的第三个方面提供一种计算机可读存储介质,其能够利用分段式线性转换,将灰度值不均匀的非8位图准确地转换为均匀的8位图像。
为了实现上述目的,本公开采用如下技术方案:
一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现上述所述的基于分段式线性转换的非均匀图像转换方法中的步骤。
为了解决上述问题,本公开的第四个方面提供一种计算机设备,其能够利用分段式线性转换,将灰度值不均匀的非8位图准确地转换为均匀的8位图像。
为了实现上述目的,本公开采用如下技术方案:
一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现上述所述的基于分段式线性转换的非均匀图像转换方法中的步骤。
本公开的有益效果是:
本公开按照灰度值不均匀的非8位图的像素点的灰度值大小,将灰度值不均匀的非8位图分成黑色区域和非黑色区域,同时将非黑色区域划分成若干区间区域;利用分段线性转换方法,计算黑色区域以及非黑色区域内各个区间区域转换为相应灰度值均匀的8位图的灰度转换曲线;利用各个线性分段的灰度转换曲线,得到灰度值均匀的8位图,最终实现准确地将灰度值不均匀的非8位图转换为灰度均匀的8位图像,这样避免了图像转换过程中信息的畸变或丢失。
附图说明
构成本公开的一部分的说明书附图用来提供对本公开的进一步理解,本公开的示意性实施例及其说明用于解释本公开,并不构成对本公开的不当限定。
图1是本公开实施例提供的一种基于分段式线性转换的非均匀图像转换方法流程图。
图2是本公开实施例提供的灰度转换曲线。
图3是本公开实施例提供的一种基于分段式线性转换的非均匀图像转换系统结构示意图。
具体实施方式
下面结合附图与实施例对本公开作进一步说明。
应该指出,以下详细说明都是例示性的,旨在对本公开提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本公开所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本公开的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
本公开的灰度值不均匀的非8位图不限于十六位图,也适用于二十四位图,三十二位图和不均匀的八位图。
图1是本公开实施例提供的一种基于分段式线性转换的非均匀图像转换方法流程图。
图1是本公开实施例提供的一种基于分段式线性转换的非均匀图像转换方法流程图。
如图1所示,本实施例的一种基于分段式线性转换的非均匀图像转换方法,包括:
S101:计算灰度值不均匀的非8位图(如16位图)的所有像素点的灰度值。
具体地,设灰度值不均匀的非8位图(如16位图)的任意像素点的灰度值为x;
将各个像素点的灰度值存储至二维数组Xoriginal中;
其中,minoriginal≤x≤maxoriginal;maxoriginal为二维数组Xoriginal中的最大值,minoriginal为二维数组Xoriginal中的最小值;minoriginal=0,m为大于或等于3的正整数。
其中,对于灰度值不均匀的16位图,m=4。
灰度值不均匀的16位图以CT图像为例:
CT图像灰度值大部分将分布到0~16之间,相对于0~255,这个数据是极其不均匀的。
读取CT图像的DICOM格式文件,预处理数组。
其中,DICOM文件时包含图片矩阵、患者信息和序列信息,这里只需要读取图片矩阵,保存到二维数组Xoriginal。
S102:按照像素点的灰度值大小,将灰度值不均匀的非8位图(如16位图)分成黑色区域和非黑色区域,同时将非黑色区域划分成若干区间区域。
例如:CT图像中大部分是空气部分(黑色区域)。如果将低灰度纳入计算,会是整体图片偏暗,不能使图像均匀显示,因此要单独计算低灰度值像素点。
S103:利用分段线性转换方法,计算黑色区域以及非黑色区域内各个区间区域转换为相应灰度值均匀的8位图的灰度转换曲线。
具体地,利用分段线性转换方法,计算黑色区域以及非黑色区域内各个区间区域转换为相应灰度值均匀的8位图的灰度转换曲线的公式为:
其中,black2m为灰度值不均匀的非8位图的黑色区域灰度值的最大值,m为大于或等于3的正整数;black8为8位图的对应灰度值不均匀的非8位图的黑色区域的灰度值的最大值;step为步长,step为正数且大于非8位图的任意像素点的灰度值与255的比值;ki为非黑色区域的第i个区间区域的转换斜率,bi为非黑色区域的第i个区间区域的偏移量,i≥1。
对于公式f(x)=kix+bi来说:
统计灰度值处于(black2m+step×(i-1),black2m+step×i]范围内的像素点总数sumcur,总像素数为sum,则映射角度alph=sumcur/sum;
第i个区间的起点坐标是:
xi-1=black2m+step×(i-1),yi-1=f(black2m+step×(i-1));
第i个区间的终点坐标:xi=black2m+step×i,yi=yi-1+(255-black8)×alph;
ki和bi的计算公式分别为:bi=yi-kixi。
对于灰度值不均匀的16位图来说,定义16位图像黑色区域最大值black16,转换后的最大值black8,则黑色区域的转换曲线是y=black8/black16*x。
计算X16=[black16+1,max]之间的灰度转换曲线。
统计像素点灰度值在X16范围内的像素点总数sum,设当前映射后最大八位灰度值cur=black8,设迭代开始灰度值n=black16+1。以step为步长(step越小,图像转换越精细,但step必须大于max/255,max=maxoriginal-minoriginal),遍历X16,开始计算每个步长内的转换曲线:
1)统计灰度值处于[n,n+step]范围内的像素点总数sumcur;
2)计算映射角度alph=sumcur/sum,当前直线起点x0=n,y0=cur,当前直线终点x1=n+step,y1=cur+(255–black8)*alph。则转换斜率k1=(y1-y0)/(x1-x0),偏移量b1=y1-k1*x1,最终转换直线为:y1=k1*x1+b1;
3)迭代cur=y1,n=n+step。
通过此步骤计算,可以得到最终转换公式,如图2所示,
S104:利用各个线性分段的灰度转换曲线,得到灰度值均匀的8位图。
具体地,计算最终的8位灰度矩阵Y=f(X),得到灰度值均匀的8位图。
本实施例按照灰度值不均匀的非8位图的像素点的灰度值大小,将灰度值不均匀的非8位图分成黑色区域和非黑色区域,同时将非黑色区域划分成若干区间区域;利用分段线性转换方法,计算黑色区域以及非黑色区域内各个区间区域转换为相应灰度值均匀的8位图的灰度转换曲线;利用各个线性分段的灰度转换曲线,得到灰度值均匀的8位图,最终实现准确地将灰度值不均匀的非8位图转换为灰度均匀的8位图像,这样避免了图像转换过程中信息的畸变或丢失。
图3是本公开实施例提供的一种基于分段式线性转换的非均匀图像转换系统结构示意图。
如图3所示,本实施例的一种基于分段式线性转换的非均匀图像转换系统,包括:
(1)灰度值计算模块,其用于计算灰度值不均匀的非8位图的所有像素点的灰度值。
灰度值不均匀的16位图以CT图像为例:
CT图像灰度值大部分将分布到0~16之间,相对于0~255,这个数据是极其不均匀的。
读取CT图像的DICOM格式文件,预处理数组。
其中,DICOM文件时包含图片矩阵、患者信息和序列信息,这里只需要读取图片矩阵,保存到二维数组Xoriginal。
(2)图像划分模块,其用于按照像素点的灰度值大小,将灰度值不均匀的非8位图分成黑色区域和非黑色区域,同时将非黑色区域划分成若干区间区域。
例如:CT图像中大部分是空气部分(黑色区域)。如果将低灰度纳入计算,会是整体图片偏暗,不能使图像均匀显示,因此要单独计算低灰度值像素点。
(3)灰度转换曲线计算模块,其用于利用分段线性转换方法,计算黑色区域以及非黑色区域内各个区间区域转换为相应灰度值均匀的8位图的灰度转换曲线。
具体地,在所述灰度转换曲线计算模块中,计算黑色区域以及非黑色区域内各个区间区域转换为相应灰度值均匀的8位图的灰度转换曲线的公式为:
其中,black2m为灰度值不均匀的非8位图的黑色区域灰度值的最大值,m为大于或等于3的正整数;black8为8位图的对应灰度值不均匀的非8位图的黑色区域的灰度值的最大值;step为步长,step为正数且大于非8位图的任意像素点的灰度值与255的比值;ki为非黑色区域的第i个区间区域的转换斜率,bi为非黑色区域的第i个区间区域的偏移量,i≥1。
对于公式f(x)=kix+bi来说:
统计灰度值处于(black2m+step×(i-1),black2m+step×i]范围内的像素点总数sumcur,总像素数为sum,则映射角度alph=sumcur/sum;
第i个区间的起点坐标是:
xi-1=black2m+step×(i-1),yi-1=f(black2m+step×(i-1));
第i个区间的终点坐标:xi=black2m+step×i,yi=yi-1+(255-black8)×alph;
ki和bi的计算公式分别为:bi=yi-kixi。
对于灰度值不均匀的16位图来说,定义16位图像黑色区域最大值black16,转换后的最大值black8,则黑色区域的转换曲线是y=black8/black16*x。
计算X16=[black16+1,max]之间的灰度转换曲线。
统计像素点灰度值在X16范围内的像素点总数sum,设当前映射后最大八位灰度值cur=black8,设迭代开始灰度值n=black16+1。以step为步长(step越小,图像转换越精细,但step必须大于max/255,max=maxoriginal-minoriginal),遍历X16,开始计算每个步长内的转换曲线:
1)统计灰度值处于[n,n+step]范围内的像素点总数sumcur;
2)计算映射角度alph=sumcur/sum,当前直线起点x0=n,y0=cur,当前直线终点x1=n+step,y1=cur+(255–black8)*alph。则转换斜率k1=(y1-y0)/(x1-x0),偏移量b1=y1-k1*x1,最终转换直线为:y1=k1*x1+b1;
3)迭代cur=y1,n=n+step。
通过此步骤计算,可以得到最终转换公式,如图2所示,
(4)图像转换模块,其用于利用各个线性分段的灰度转换曲线,得到灰度值均匀的8位图。
具体地,计算最终的8位灰度矩阵Y=f(X),得到灰度值均匀的8位图。
本实施例按照灰度值不均匀的非8位图的像素点的灰度值大小,将灰度值不均匀的非8位图分成黑色区域和非黑色区域,同时将非黑色区域划分成若干区间区域;利用分段线性转换方法,计算黑色区域以及非黑色区域内各个区间区域转换为相应灰度值均匀的8位图的灰度转换曲线;利用各个线性分段的灰度转换曲线,得到灰度值均匀的8位图,最终实现准确地将灰度值不均匀的非8位图转换为灰度均匀的8位图像,这样避免了图像转换过程中信息的畸变或丢失。
在另一实施例中,该系统,还包括:
存储模块,其用于将各个像素点的灰度值存储至二维数组Xoriginal中;设灰度值不均匀的非8位图的任意像素点的灰度值为x;其中,minoriginal≤x≤maxoriginal;maxoriginal为二维数组Xoriginal中的最大值,minoriginal为二维数组Xoriginal中的最小值;minoriginal=0,m为大于或等于3的正整数。
在另一实施例中,还提供了一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现如图1所示的基于分段式线性转换的非均匀图像转换方法中的步骤。
在另一实施例中,还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如图1所示的基于分段式线性转换的非均匀图像转换方法中的步骤。
本领域内的技术人员应明白,本公开的实施例可提供为方法、系统、或计算机程序产品。因此,本公开可采用硬件实施例、软件实施例、或结合软件和硬件方面的实施例的形式。而且,本公开可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器和光学存储器等)上实施的计算机程序产品的形式。
本公开是参照根据本公开实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)或随机存储记忆体(RandomAccessMemory,RAM)等。
以上所述仅为本公开的优选实施例而已,并不用于限制本公开,对于本领域的技术人员来说,本公开可以有各种更改和变化。凡在本公开的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本公开的保护范围之内。
Claims (10)
1.一种基于分段式线性转换的非均匀图像转换方法,其特征在于,包括:
计算灰度值不均匀的非8位图的所有像素点的灰度值;
按照像素点的灰度值大小,将灰度值不均匀的非8位图分成黑色区域和非黑色区域,同时将非黑色区域划分成若干区间区域;
利用分段线性转换方法,计算黑色区域以及非黑色区域内各个区间区域转换为相应灰度值均匀的8位图的灰度转换曲线;
利用各个线性分段的灰度转换曲线,得到灰度值均匀的8位图。
2.如权利要求1所述的一种基于分段式线性转换的非均匀图像转换方法,其特征在于,该方法还包括:
设灰度值不均匀的非8位图的任意像素点的灰度值为x;
将各个像素点的灰度值存储至二维数组Xoriginal中;
其中,minoriginal≤x≤maxoriginal;maxoriginal为二维数组Xoriginal中的最大值,minoriginal为二维数组Xoriginal中的最小值;minoriginal=0,m为大于或等于3的正整数。
3.如权利要求1所述的一种基于分段式线性转换的非均匀图像转换方法,其特征在于,利用分段线性转换方法,计算黑色区域以及非黑色区域内各个区间区域转换为相应灰度值均匀的8位图的灰度转换曲线的公式为:
其中,black2m为灰度值不均匀的非8位图的黑色区域灰度值的最大值,m为大于或等于3的正整数;black8为8位图的对应灰度值不均匀的非8位图的黑色区域的灰度值的最大值;step为步长,step为正数且大于非8位图的任意像素点的灰度值与255的比值;ki为非黑色区域的第i个区间区域的转换斜率,bi为非黑色区域的第i个区间区域的偏移量,i≥1。
4.如权利要求3所述的一种基于分段式线性转换的非均匀图像转换方法,其特征在于,ki和bi的计算方法为:
统计灰度值处于(black2m+step×(i-1),black2m+step×i]范围内的像素点总数sumcur,总像素数为sum,则映射角度alph=sumcur/sum;
第i个区间的起点坐标是:
xi-1=black2m+step×(i-1),yi-1=f(black2m+step×(i-1));
第i个区间的终点坐标:xi=black2m+step×i,yi=yi-1+(255-black8)×alph;
ki和bi的计算公式分别为:bi=yi-kixi。
5.一种基于分段式线性转换的非均匀图像转换系统,其特征在于,包括:
灰度值计算模块,其用于计算灰度值不均匀的非8位图的所有像素点的灰度值;
图像划分模块,其用于按照像素点的灰度值大小,将灰度值不均匀的非8位图分成黑色区域和非黑色区域,同时将非黑色区域划分成若干区间区域;
灰度转换曲线计算模块,其用于利用分段线性转换方法,计算黑色区域以及非黑色区域内各个区间区域转换为相应灰度值均匀的8位图的灰度转换曲线;
图像转换模块,其用于利用各个线性分段的灰度转换曲线,得到灰度值均匀的8位图。
6.如权利要求5所述的一种基于分段式线性转换的非均匀图像转换系统,其特征在于,该系统,还包括:
存储模块,其用于将各个像素点的灰度值存储至二维数组Xoriginal中;设灰度值不均匀的非8位图的任意像素点的灰度值为x;其中,minoriginal≤x≤maxoriginal;maxoriginal为二维数组Xoriginal中的最大值,minoriginal为二维数组Xoriginal中的最小值;minoriginal=0,m为大于或等于3的正整数。
7.如权利要求5所述的一种基于分段式线性转换的非均匀图像转换系统,其特征在于,在所述灰度转换曲线计算模块中,计算黑色区域以及非黑色区域内各个区间区域转换为相应灰度值均匀的8位图的灰度转换曲线的公式为:
其中,black2m为灰度值不均匀的非8位图的黑色区域灰度值的最大值,m为大于或等于3的正整数;black8为8位图的对应灰度值不均匀的非8位图的黑色区域的灰度值的最大值;step为步长,step为正数且大于非8位图的任意像素点的灰度值与255的比值;ki为非黑色区域的第i个区间区域的转换斜率,bi为非黑色区域的第i个区间区域的偏移量,i≥1。
8.如权利要求6所述的一种基于分段式线性转换的非均匀图像转换系统,其特征在于,ki和bi的计算方法为:
统计灰度值处于(black2m+step×(i-1),black2m+step×i]范围内的像素点总数sumcur,总像素数为sum,则映射角度alph=sumcur/sum;
第i个区间的起点坐标是:
xi-1=black2m+step×(i-1),yi-1=f(black2m+step×(i-1));
第i个区间的终点坐标:xi=black2m+step×i,yi=yi-1+(255-black8)×alph;
ki和bi的计算公式分别为:bi=yi-kixi。
9.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现如权利要求1-4中任一项所述的基于分段式线性转换的非均匀图像转换方法中的步骤。
10.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现如权利要求1-4中任一项所述的基于分段式线性转换的非均匀图像转换方法中的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910261626.0A CN109993691B (zh) | 2019-04-02 | 2019-04-02 | 基于分段式线性转换的非均匀图像转换方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910261626.0A CN109993691B (zh) | 2019-04-02 | 2019-04-02 | 基于分段式线性转换的非均匀图像转换方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109993691A true CN109993691A (zh) | 2019-07-09 |
CN109993691B CN109993691B (zh) | 2022-11-11 |
Family
ID=67132193
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910261626.0A Expired - Fee Related CN109993691B (zh) | 2019-04-02 | 2019-04-02 | 基于分段式线性转换的非均匀图像转换方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109993691B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103530847A (zh) * | 2013-09-24 | 2014-01-22 | 电子科技大学 | 一种红外图像增强方法 |
CN105678750A (zh) * | 2015-12-31 | 2016-06-15 | 上海联影医疗科技有限公司 | 医学图像的灰度映射曲线生成方法及装置 |
WO2017096946A1 (zh) * | 2015-12-07 | 2017-06-15 | 乐视控股(北京)有限公司 | 图像高频信息的定位方法和装置 |
-
2019
- 2019-04-02 CN CN201910261626.0A patent/CN109993691B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103530847A (zh) * | 2013-09-24 | 2014-01-22 | 电子科技大学 | 一种红外图像增强方法 |
WO2017096946A1 (zh) * | 2015-12-07 | 2017-06-15 | 乐视控股(北京)有限公司 | 图像高频信息的定位方法和装置 |
CN105678750A (zh) * | 2015-12-31 | 2016-06-15 | 上海联影医疗科技有限公司 | 医学图像的灰度映射曲线生成方法及装置 |
Non-Patent Citations (2)
Title |
---|
史玉林等: "用于图像对比度增强的直方图锥形拉伸算法", 《计算机工程》 * |
江忠良: "一种非线性组合变换的图像增强技术", 《电子世界》 * |
Also Published As
Publication number | Publication date |
---|---|
CN109993691B (zh) | 2022-11-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110321920A (zh) | 图像分类方法、装置、计算机可读存储介质和计算机设备 | |
CN105354809B (zh) | 一种基于输出图像像素位置索引的预畸变方法和装置 | |
CN108460730B (zh) | 一种图像处理方法及其装置 | |
CN101779962B (zh) | 一种增强医学x射线影像显示效果的方法 | |
EP2169959A2 (en) | Lossless compression-encoding device | |
CN103914852B (zh) | 基于cuda的dicom医学影像动态非线性调窗方法 | |
CN104380340A (zh) | 基于邻域关系对像素进行评分和调整以揭示图像中的数据 | |
CN111754928B (zh) | 一种显示屏分层校正方法及装置 | |
CN107292819A (zh) | 一种基于边缘细节保护的红外图像超分辨率重建方法 | |
CN103281537B (zh) | 一种图像动态范围压缩方法及装置 | |
CN111489322A (zh) | 给静态图片加天空滤镜的方法及装置 | |
CN105338214B (zh) | 图像处理方法和图像处理装置 | |
CN109285527A (zh) | 显示装置及其图像处理方法 | |
CN108510562A (zh) | 基于图像分形纹理的数码迷彩图案生成方法 | |
CN108184122A (zh) | 基于多元线性回归的图像可逆信息隐藏方法和装置 | |
EP2169958A2 (en) | Lossless compression-encoding device and decoding device for image data | |
CN110047077B (zh) | 一种用于以太坊共识机制的图像处理方法 | |
CN106548223A (zh) | 图像条码编码方法、解码方法、编码装置及解码装置 | |
US20170195520A1 (en) | Image processing method, apparatus, and image forming device | |
CN113256657B (zh) | 一种高效医学图像分割方法及系统、终端、介质 | |
CN109993691A (zh) | 基于分段式线性转换的非均匀图像转换方法及系统 | |
US9159120B2 (en) | Image processing apparatus, image processing method, and storage medium | |
CN111862071A (zh) | 一种基于ct图像测量腰1椎体ct值的方法 | |
CN110458754B (zh) | 图像生成方法及终端设备 | |
CN112102216B (zh) | 自适应权重的全变差图像融合方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20221111 |