专利名称：调频调幅混合网点层次连续调控制方法
技术领域：
本发明属于图像硬拷贝复制领域的半色调网点生成技术领域，涉及一种调频调幅混合网点层次连续调控制方法，具体涉及一种在多位成像深度设备上的调频调幅混合网点层次连续调控制方法。
背景技术：
图象的硬拷贝复制主要涉及打印机及高档印刷制版设备的挂网制版技术。用于图像硬拷贝复制的挂网技术又称为数字图象半色调技术。数字图象半色调技术可以分为两类，分别是调幅挂网和调频挂网。调幅挂网又称为聚集点有序抖动技术，其特征是所生成的半色调图象的染色点在几何位置上是两两相邻地聚合在一起的，从而形成一簇一簇的染色区域，这些染色区域又被称为网点。由于聚集点有序抖动技术采用控制网点面积的方法再现原稿图象灰度，因而被称为调幅网点。
基于调频网与调幅网多位混合挂网的技术，现今已有技术不多，典型的可以参见已有的专利《一种基于双反馈的误差扩散调频挂网方法》(申请号200510068127.8)，《在多位成像深度设备上的调频调幅混合网点网型控制方法》(申请号200510116636.3)，《多位成像深度设备上的调频调幅混合半色调图像处理方法》(申请号200510116635.9)。已有技术中混合网点发展的轨迹主要遵循着误差扩散基本理论主线展开，首先出现了基于双反馈的误差扩散调频挂网方法，该方法在挂网时，通过采用基于误差扩散的基本算法基础上的双反馈技术，实现了将传统的调频网与调幅网进行了有机的结合，达到了混合加网的效果。采用该方法，可以在输出设备分辨率较低的情况下，有效地解决了网点层次容易丢失的问题，并能同时降低因随机分布特性所带来的纹理，因而可以输出高质量层次饱满的半色调图像。紧接着该技术的出现，又出现了多位成像深度设备上的调频调幅混合半色调图像处理方法，该技术采用了多位混合网点的动态控制算法，可以在多位成像深度的输出设备上输出高质量层次饱满的多位半色调图像，并有效地解决原有的一位深度设备下输出的混合网点边沿锯齿现象，保证在低分辨率下输出高分辨率效果的层次阶调连续的调频调幅混合网点。为完善上述的多位网点，又出现了在多位成像深度设备上的调频调幅混合网点网型控制方法，引入相邻方向输出灰度动态统计算法，保证了网点形状的可控性，解决了因误差扩散的随机性而带来的网型变化的随意性，可以在原有技术基础上，根据设备成像对网点的要求，充分发挥多位成像设备的特性，能够在低分辨率下输出网点形状变化可控的调频调幅半色调混合挂网效果，解决了混合半色调网点在实际输出时的颗粒感问题，保证了层次阶调的平滑效果。
已有的上述技术主要是基于一位或多位成像深度设备基础上的网点生成通用算法，在实际设备输出网点时，由于已有技术中的多位网点在不同层次阶调上的网型再现完全基于误差随机扩散理论，采用了不同层次阶调的网型动态控制算法，使得输出的多位网点形状和特性具备了多位网点所应具有的网点层次递增均匀性等特点，在实际设备输出时主要表现为网点大小一致，并使得输出图像的质量有所提高，大大减少了混合网点固有的颗粒感问题，输出图像平滑细腻。
但是，多位混合网点网型动态控制算法中由于对网点输出层次进行了强制性的改变，使得网型在一定范围内得以控制，但这种处理或多或少破坏了原有误差扩散算法，加之该控制算法实施的有效区域主要在低密度和中间密度区域内分别进行操作，相互之间没有任何算法上的联系，同时在不同密度区域划分上，分界点采用固定取值的机制，最终导致在低密度区与中间密度区以及中间密度区与高密度区分界点上产生由于不同的处理手段导致的网点灰度层次跳变，影响了半色调图像的输出质量，使得输出图像在层次阶调上存在跳变反差，在多位成像设备输出时尤为明显。
发明内容
针对现有技术中存在的缺陷，本发明的目的在于提供一种在多位成像深度设备上的调频调幅混合网点层次连续调控制方法。该方法在多位成像深度设备上的调频调幅混合网点层次连续调控制技术及原有的多位成像深度的混合网点生成技术及网型控制等多项技术基础上，能够解决现有技术中因不同层次阶调所采用各异的动态网型控制手段而造成的网点灰度层次跳变，从而有效地发挥了多位成像深度设备的层次阶调再现效果，提高了半色调混合网点层次连续调输出的质量。
为达到以上目的，本发明的技术方案如下一种在多位成像深度设备上的调频调幅混合网点层次连续调控制方法，包括以下步骤步骤一首先根据设备的图像位深度n，将0到255的层次阶调范围均分为2n-1个不同区域 ，(R1，R2]，…，(Ri-1，Ri]，…(R2n-2，255]与其对应的位输出点阵范围如下(0，Out1)，(Out1，Out2)，…，(Outi-1，Outi)，…(Out2n-2，11…1)其中Ri为动态区域分界点，其大小依据随后的抖动算法而动态取值。
Outi为n位深度二进制表示；在每个区域范围内取中间点的阈值Mi作为该区域的阈值比较参数；更进一步，Ri的取值方法采用了层次阶调可调抖动算法，方法如下参见图例1，给定某一分界点的抖动基准值Ri，建立该分界点抖动查找表TDi，其中各表项值的获得可通过在Ri值附近随机取值得到，随机方式可以有所不同，只要保证各表项值在给定Ri值附近随机均匀的分布即可。
更进一步，本专利采用下面的计算公式得到抖动查找表各表项值TDi
＝Ri-1；TDi[1]＝Ri-2TDi[2]＝Ri+1；TDi[3]＝Ri-3TDi[4]＝Ri+4；TDi[5]＝Ri-3TDi[6]＝Ri-4；TDi[7]＝Ri+2 公式1TDi[8]＝Ri-2；TDi[9]＝Ri+4TDi[10]＝Ri+2；TDi[11]＝Ri+1TDi[12]＝Ri-1；TDi[13]＝Ri+3TDi[14]＝Ri+3；TDi[15]＝Ri-4步骤二在
的层次阶调范围内，设定n位成像深度输出概率阈值；更进一步，在
的层次阶调范围内，设定2n-2个整数作为n位成像深度输出概率阈值L1，L2，…，Li，…，L2n-2，该概率阈值的大小可调，可根据实际设备输出要求作相应改变，以达到最佳输出质量；步骤三在抖动表TDi-1，TDi中，针对每一个处理点分别随机取出两个抖动表项赋值给Ri-1，Ri作为处理层次区域分界点，在2n-1个不同区域(Ri-1，Ri]中，基于误差扩散双反馈的调频调幅混合网挂网方法，分别作相应运算处理，其具体处理方法包括以下步骤(1)首先扫描并输入原稿图像，对原稿图像的最终象素值g”(m，t)进行域值比较运算T，运算的结果被转换成为半色调图像的输出象素值b(m，t)；
更进一步，扫描并输入原稿图像时采用双向扫描处理的方法，即针对连续扫描的原稿每一行，扫描首先是从左到右扫描一行，紧接着下一行的扫描采用由右到左进行，依次类推交错进行，直到所有行扫描完毕。
(2)将步骤(1)中得到的输出象素值b(m，t)和被求域值的原稿图像的中间象素值g’(m，t))进行比较，并计算两者的差值即误差值e(m，t)；(3)将误差值e(m，t)通过一个扩散滤波器e与一定的权重分配值作乘积运算后扩散到当前处理象素周围未处理的象素上，被扩散位置的原稿象素值g(m，t)将和被扩散到此象素上的误差值加全求和得到新的原稿图像的中间象素值g’(m，t)；更进一步，这里的误差扩散滤波器e，采用下表的扩散原理和权重分配系数** d5d3d2d4d5d4d2d1d2d3d2d1其中**代表了当前点像素位置，其它位置的算术比值代表了与当前像素相对位置上的扩散权重值，取值范围在
，且满足2×d1+4×d2+2×d3+2×d4+2×d5∈
。
(4)与上面步骤(2)、(3)操作并行，将当前处理后的输出象素值b(m，t)与第二个扩散滤波器w进行乘积运算并进行抖动处理后扩散到周围相应的未处理象素上并与上面两步中误差扩散得到的的原稿图像的中间象素值g’(m，t)作最终的加全求和作为原稿图像的最终象素值g”(m，t)；更进一步，这里扩散滤波器w的扩散方向设置如下** w0w3w2w1其中扫描方向由左到右，**代表了当前点像素位置，其它位置的参数代表了与当前像素相对位置上的扩散权重值，取值范围在
，且满足w0+w1+w2+w3∈
。
更进一步，对扩散滤波器w的抖动算法采用如下的方法进行处理fRand＝(R(m，t)/R_MAX-0.5)×cDitherdw0＝w0-fRanddw2＝w2+fRanddw1＝w1+fRand
dw3＝w3-fRand上述公式中fRand为抖动微调参数；R(m，t)为扫描当前点的随机取值参数；R_MAX为随机参数R(i)的最大值；cDither为抖动幅度调整参数，该参数决定了调幅特性质量优劣；dw0～dw3为抖动后扩散滤波器w不同方向上的权重扩散值。
(5)循环上面步骤1)到步骤4)，直到将所有输入像素g(m，t)处理完成。
步骤四在抖动表TDi-1，TDi中，针对每一个处理点分别随机取出两个抖动表项赋值给Ri-1，Ri作为处理层次区域分界点，在区域(Ri-1，Ri)内实现混合挂网的同时，根据概率阈值Fi及当前点的网型控制累积值S，结合原有的动态可变层次输出机制采用相邻方向输出灰度动态统计算方法其动态可变层次输出方法包括如下计算公式根据当前点的网型控制累积值S产生当前点的伪随机值Fi＝random(S)， 公式2其中伪随机函数random可在编译环境下自动产生，Fi∈
；采用如下的公式动态计算输出点阵数据OUT=I1=0--------Fi∈
Out1------Fi∈(L1,L2]···Outi-1-----Fi∈(Li-2,Lo-1]I2=Outi------Fi∈(Li-1,Li]Outi+1-----Fi∈(L1,Li+1]···11...1-----Fi∈(L2n-2,255)]]>公式3相邻方向输出灰度动态计算输出点阵数据控制网型具体步骤如下(1)根据多位成像深度设备稳定性及输出线性化特征，确定网型控制调整密度区域；(2)根据多位成像深度设备脉冲宽度调节技术参数要求，结合不同密度区已有随机网型及设备输出稳定性对不同密度区域网点形状的要求，确定需要调整的网点形状；
(3)根据步骤(2)确定的需要调整的网点形状，结合输出设备的成像特点，确定目标网点形状及应该避免出现的网点形状；(4)最终依据网点密度区及已有的输出灰度层次，采用相邻方向输出灰度动态统计算法，确定多位混合半色调二进制网点输出范围参数，算法中的控制点设定方法包括以下步骤1)设定当前点前一个点的输出灰度值a，a前一个点的输出灰度值为b；2)设定当前点前一行中的当前点位置的点的灰度值为c，c前一个点灰度值为d，c后一个点灰度值为e；3)设定输出灰度值为c的前一行相同位置的点的灰度值为f；4)设定当前点周围输出灰度统计参数Sum1＝a+c+dSum2＝a+bSum3＝c+f其中，Sum1表示当前点周围2×2矩形区域灰度加全和，Sum2表示当前点水平方向上前两个输出灰度加全和，Sum3表示当前点垂直方向上前两个输出灰度加全和。
本发明的效果在于使用本发明提供的方法，可以在原有误差扩散双反馈的多位调频调幅混合网挂网方法的基础上，根据设备成像对网点的要求，充分发挥了多位成像设备的特性，在保证不同层次阶调网型的同时，消除了各阶调之间的层次跳变。采用本发明所述的方法可以有效地消除由于网型控制带来的层次阶调跳变问题，从而保证半色调图像输出的层次平滑度。
图1是多位输出设备层次阶调可调抖动算法原理图；图2是2位输出设备层次阶调可调抖动算法原理图；图3是调整前低密度区的网点形状示意图；图4是调整前中密度区的网点形状示意图；图5是调整后的低密度区网点形状示意图；图6是调整后的中间密度区网点形状示意图；图7是相邻输出灰度值动态统计图；图8是层次抖动前分界点连续调网点输出效果图；图9是层次抖动后分界点连续调网点输出效果图。
具体实施例方式
下面结合附图和实施例对本发明作进一步的描述
本实施例基于已有的误差扩散双反馈的多位成像深度的调频调幅混合网挂网方法及网型控制等多项技术基础上，引入动态层次阶调划分的抖动控制算法，实现了在多位成像深度设备上的调频调幅混合网点层次连续调控制技术，其具体实施方式
如下如图2所示，本实施例中，设定输出设备的成像深度为两位，即n＝2，层次分界点R1和R2的抖动基准值分别设定为84和171，抖动表的生成将根据两个抖动基准值计算，设备已有的脉冲宽度调节技术对网点的要求是(设定纯白为255，纯黑为0)低密度区(171～255)中由于网点没有搭接，相互独立，要求单象素的网点输出值为3(二进制11)，两象素网点输出组合为3与1的组合，三个象素或更多的网点组合要求不能存在3与3的水平或纵向排列。
中密度区(84～171)中由于独立的网点于与搭接后的网点同时存在，因此除了低密度区的要求外，要遵循网点大小同层次多数一致性原则和网点象素排列层次的递增性原则，从而使得网点的输出更接近于2位成像深度网点的特征。
高密度区(0～84)中由于本实施例中输出设备为低分辩激光打印机，因此，其打印线性特征对该密度区的要求不设限制，该密度区仍按原有多位网点的输出原则输出。
根据上面输出设备的要求，其具体实际解决方案如下步骤一首先根据设备的图像位深度n＝2，将0到255的层次阶调范围均分为2n-1＝3个不同区域
，(R1，R2]，(R2，255]，其中，R1＝84，R2＝171与其对应的位输出点阵范围如下(11，10)，(10，01)，(01，00)R1与R2取值的获得，通过两个抖动查找表的得到，本专利中针对R1与R2的两个抖动查找表TD1与TD2的设置通过公式1计算如下TD1＝{83，82，85，81，88，81，80，86，82，88，86，85，83，87，87，80}TD2＝{170，169，172，168，175，168，167，173，169，175，173，172，170，174，174，167}步骤二在
的层次阶调范围内，设定2n-2＝2个整数作为n位成像深度输出概率阈值F1＝8，F2＝24，通过该阈值概率控制了11，10，01，00的输出概率，从而达到2位混合网点的输出效果；
步骤三在上面3个不同区域中，根据原稿图像输入象素值确定操作区域范围，动态随机从抖动表TD1与TD2取出相应的R1和R2，基于已有的误差扩散双反馈的调频调幅混合网挂网方法，分别作相应运算处理，其具体处理方法如下(1)首先扫描并输入原稿图像，对原稿图像的最终象素值g”(m，t)进行域值比较运算T，运算的结果被转换成为半色调图像的输出象素值b(m，t)；这里扫描原稿象素数据时，为避免由于扫描方向与随机分布的网点频率产生干涉现象，俗称“撞网”，本实施例中采取双向扫描的方法。即针对连续扫描的原稿每一行，扫描首先是从左到右扫描一行，紧接着下一行的扫描采用由右到左进行，依次类推交错进行，直到所有行扫描完毕。
(2)将此步骤(1)中得到的输出象素值b(m，t)和被求域值的原稿图像的中间象素值g’(m，t))进行比较，并计算两者的差值即误差值e(m，t)；(3)将误差值e(m，t)通过一个扩散滤波器e与一定的权重分配值作乘积运算后扩散到当前处理象素周围未处理的象素上，被扩散位置的原稿象素值g(m，t)将和被扩散到此象素上的误差值加全求和得到新的原稿图像的中间象素值g’(m，t)。所述的误差扩散滤波器e，采用下表的扩散原理和权重分配系数** d5d3d2d4d5d4d2d1d2d3d2d1其中**代表了当前点象素位置，其它位置的算术比值代表了与当前象素相对位置上的扩散权重值，在本实施例中，采用如下参数取值d1＝1/44 d2＝2/44 d3＝5/44 d4＝4/44 d5＝8/44上面步骤2，3完成了第一次扩散反馈操作，实现了误差扩散的基本原理。
(4)与上面步骤(2)、(3)操作并行，将当前处理后的输出象素值b(m，t)与第二个扩散滤波器w进行乘积运算并进行抖动处理后扩散到周围相应的未处理象素上并与上面两步中误差扩散得到的原稿图像的中间象素值g’(m，t)作最终的加全求和作为最终的输入象素值g”(m，t)。这里第二个滤波器w的扩散方向设置如下** w0w3w2w1
其中扫描方向由左到右，**代表了当前点象素位置，其它位置的参数代表了与当前象素相对位置上的扩散权重值，取值范围在
，且满足wsum＝(w0+w1+w2+w3)∈
。
本实施例中的抖动算法采用如下的算法进行处理fRand＝(R(m，t)/R_MAX-0.5)×cDitherdw0＝w0-fRanddw2＝w2+fRanddw1＝w1+fRanddw3＝w3-fRand上述公式中fRand为抖动微调参数；R(m，t)为扫描当前点的随机取值参数；R_MAX为随机参数R(i)的最大值；cDither为抖动幅度调整参数，该参数决定了调幅特性质量优劣；dw0～dw3为抖动后滤波器w不同方向上的权重扩散值。
上述步骤(4)完成第二次扩散反馈操作，实现了调频网的调幅特性。
在本实施例中，设定上面的参数如下w0＝w2＝0.175，w1＝w3＝0.025，则wsum＝0.4cDither＝0.2挂网过程中，通过调整wsum，改变调频网点大小，通过w0～w3取值分配，来调整调频网点的形状。
(5)循环上面步骤(1)到步骤(4)，直到将所有输入象素g(m，t)处理完成。
步骤四在3个均等区域内实现混合挂网的同时，根据概率阈值L1＝8，L2＝24及当前点的网型控制累积值S，动态可变层次输出机制，算法如下公式根据当前点的网型控制累积值S，结合公式1产生当前点的伪随机值Fi＝random(S)公式2其中伪随机函数random可在编译环境下自动产生。F∈
根据公式2，动态计算输出点阵数据。
相邻方向输出灰度动态计算输出点阵数据控制网型实现如下1)根据多位成像深度设备稳定性及输出线性化特征，确定网型控制调整密度区域。
本实施例中确定的调整密度区域是低密度区(R2～255)和中密度区(R1～R2)。
2)根据多位成像深度设备的脉冲宽度调节技术参数要求，结合不同密度区已有随机网型及设备输出稳定性对不同密度区域网点形状的要求，确定需要调整的网点形状；本实施例中，低密度区(R2～255)需要调整的网点形状可参考图3，中密度区(R1～R2)需要调整的网点形状参考图4。
3)根据步骤2)确定的需要调整的网点形状，结合输出设备的成像特点，确定目标网点形状及应该避免出现的网点形状。
本实施例中，低密度区目标网点形状参见图5，图中符号*代表需要调整的输出灰度，其值为1或0；中密度区目标网点形状参见图6。图5与图6中，理想状态下，2位网点的中心灰度值为3。
4)最终依据网点密度区及已有的输出灰度层次，采用相邻方向输出灰度动态统计算法，确定多位混合半色调二进制网点输出范围参数。根据步骤二，步骤三所确定的目标，本实施例采用如下的伪代码实现如上步骤根据相邻方向输出灰度动态统计算法中，控制点设定的方法如下1)设定当前点前一个点的输出灰度值a，a前一个点的输出灰度值为b；2)设定当前点前一行中的当前点位置的点的灰度值为c，c前一个点灰度值为d，c后一个点灰度值为e；3)设定输出灰度值为c的前一行相同位置的点的灰度值为f。
4)设定当前点周围输出灰度统计参数Sum1＝a+c+dSum2＝a+bSum3＝c+f上面a，b，c，d，e，f∈
(实例中，二位网点具有4个输出层次)，其具体分布图参见附图7(相邻输出灰度动态统计原理图)。
二位网点统计伪码实现逻辑如下R1＝TD1[random()％15]R2＝TD2[random()％15](其中伪随机函数random()可在编译环境下自动产生)IF(输入象素值∈[R1～R2])中密度区{IF(Sum1＝＝9)输出灰度值二进制范围[10，01]
ELSE IF(a！＝0&&c！＝0&&d！＝0&&输入象素值∈[123，R2])输出灰度值二进制范围
ELSE IF(Sum3＞＝5)&&输入象素值∈[123，R2])输出灰度值二进制范围
ELSE IF(Sum1＜3&&输入象素值∈[139，R2])输出灰度值二进制范围[11，00]ELSE IF(Sum2＝＝6&&输入象素值∈[123，R2])输出灰度值二进制范围
ELSE IF(Sum1＞＝7&&输入象素值∈[139，R2])输出灰度值二进制范围
ELSE依照“动态可变层次机制”输出}IF(输入象素值∈(R2，255])低密度区{IF(a＝＝3||c＝＝3){IF((a＝＝3&&c＝＝3)||(a＝＝3&&(c！＝3&&c！＝0))||(a！＝3&&a！＝0)&&c＝＝3))输出灰度值二进制范围
ELSE输出灰度值二进制范围
}ELSE IF(((a！＝3&&a！＝0)&&(c！＝3&&c！＝0))||(c！＝3&&c！＝0&&d＝＝0&&e＝＝0)||a！＝3&&a！＝0&&b＝＝3)输出灰度值二进制范围
ELSE输出灰度值二进制范围[11，00]}基于2位深度成像深度设备，结合上面的具体实施例的各个步骤，可以实现2位调频调幅混合挂网，其挂网效果图与原有的2位混合网点对比可参见图8，图9，其中图中左图为低密度区与中密度区层次分界点连续调网点样例，右图为中密度区与高密度区层次分界点连续调网点样例。通过比较可以明显看出，采用本发明所述的方法在多位成像深度设备上输出的调频调幅混合网点的大小、形状基本趋于统一，而且网点的输出效果更加平整、干净，从而有效地消除现有技术中由于网型控制带来的层次阶调跳变问题。
本发明所述的方法并不限于具体实施方式
中所述的实施例，本领域技术人员根据本发明的技术方案得出其他的实施方式，同样属于本发明的技术创新范围。
权利要求
1.一种调频调幅混合网点层次连续调控制方法，包括以下步骤步骤一首先根据设备的图像位深度n，将0到255的层次阶调范围均分为2n-1个不同区域
，(R1，R2]，...，(Ri-1，Ri]，...(R2n-2，255]与其对应的位输出点阵范围如下(0，Out1)，(Out1，Out2)，...，(Outi-1，Outi)，...(Out2n-2，11...1)其中Ri为动态区域分界点，Ri的取值方法采用了层次阶调可调抖动算法，方法如下给定某一分界点的抖动基准值Ri，建立该分界点抖动查找表TDi，其中各表项值的获得可通过在Ri值附近随机取值得到；Outi为n位深度二进制表示；在每个区域范围内取中间点的阈值Mi作为该区域的阈值比较参数；步骤二在
的层次阶调范围内，设定n位成像深度输出概率阈值；步骤三在抖动表TDi-1，TDi中，针对每一个处理点分别随机取出两个抖动表项赋值给Ri-1，Ri作为处理层次区域分界点，在2n-1个不同区域(Ri-1，Ri]中，基于误差扩散双反馈的调频调幅混合网挂网方法，分别作相应运算处理；步骤四在抖动表TDi-1，TDi中，针对每一个处理点分别随机取出两个抖动表项赋值给Ri-1，Ri作为处理层次区域分界点，在区域(Ri-1，Ri)内实现混合挂网的同时，根据概率阈值Fi及当前点的网型控制累积值S，结合原有的动态可变层次输出机制采用相邻方向输出灰度动态统计计算方法其动态可变层次输出方法包括如下计算公式根据当前点的网型控制累积值S产生当前点的伪随机值Fi＝random(S)， 公式2其中伪随机函数random可在编译环境下自动产生，Fi∈
；采用如下公式动态计算输出点阵数据OUT=l1=0--------Fi∈
Out1------F1∈(L1,L2]···Outi-1-----Fi∈(Li-2,Li-1]l2=Outi------Fi∈(Li-1,Li]Outi-----Fi∈(Li,Li+1]···11...1-----Fi∈(L2n-2,255]]]>公式3。
2.如权利要求1所述的一种调频调幅混合网点层次连续调控制方法，其特征是步骤一中，采用下面的计算公式得到抖动查找表TDi各表项值TDi
＝Ri-1；TDi[1]＝Ri-2TDi[2]＝Ri+1；TDi[3]＝Ri-3TDi[4]＝Ri+4；TDi[5]＝Ri-3TDi[6]＝Ri-4；TDi[7]＝Ri+2TDi[8]＝Ri-2；TDi[9]＝Ri+4TDi[10]＝Ri+2；TDi[11]＝Ri+1TDi[12]＝Ri-1；TDi[13]＝Ri+3TDi[14]＝Ri+3；TDi[15]＝Ri-4公式1。
3.如权利要求1所述的一种调频调幅混合网点层次连续调控制方法，其特征是步骤二中，在
的层次阶调范围内，设定2n-2个整数作为n位成像深度输出概率阈值L1，L2，...，Li，...，L2n-2，该概率阈值的大小可调，可根据实际设备输出要求作相应改变，以达到最佳输出质量。
4.如权利要求1所述的一种调频调幅混合网点层次连续调控制方法，其特征是步骤三中，具体处理方法包括以下步骤(1)首先扫描并输入原稿图像，对原稿图像的最终象素值g”(m，t)进行域值比较运算T，运算的结果被转换成为半色调图像的输出象素值b(m，t)；(2)将步骤(1)中得到的输出象素值b(m，t)和被求域值的原稿图像的中间象素值g’(m，t))进行比较，并计算两者的差值即误差值e(m，t)；(3)将误差值e(m，t)通过一个扩散滤波器e与一定的权重分配值作乘积运算后扩散到当前处理象素周围未处理的象素上，被扩散位置的原稿象素值g(m，t)将和被扩散到此象素上的误差值加全求和得到新的原稿图像的中间象素值g’(m，t)；(4)与上面步骤(2)、(3)操作并行，将当前处理后的输出象素值b(m，t)与第二个扩散滤波器w进行乘积运算并进行抖动处理后扩散到周围相应的未处理象素上并与上面两步中误差扩散得到的的原稿图像的中间象素值g’(m，t)作最终的加全求和作为原稿图像的最终象素值g”(m，t)；(5)循环上面步骤(1)到步骤(4)，直到将存储单元M中的所有输入像素g(m，t)处理完成。
5.如权利要求4所述的一种调频调幅混合网点层次连续调控制方法，其特征是步骤(1)中，扫描并输入原稿图像时采用双向扫描的方法，即针对连续扫描的原稿每一行，扫描首先是从左到右扫描一行，紧接着下一行的扫描采用由右到左进行，依次类推交错进行，直到所有行扫描完毕。
6.如权利要求4或5所述的一种调频调幅混合网点层次连续调控制方法，其特征是步骤(3)中，所述的误差扩散滤波器e，采用下表的扩散原理和权重分配系数** d5d3d2d4d5d4d2d1d2d3d2d1其中**代表了当前点像素位置，其它位置的算术比值代表了与当前像素相对位置上的扩散权重值，取值范围在
，且满足2×d1+4×d2+2×d3+2×d4+2×d5∈
。
7.如权利要求4或5所述的一种调频调幅混合网点层次连续调控制方法，其特征是步骤(4)中，所述的扩散滤波器w的扩散方向设置如下** w0w3w2w1其中扫描方向由左到右，**代表了当前点像素位置，其它位置的参数代表了与当前像素相对位置上的扩散权重值，取值范围在
，且满足w0+w1+w2+w3∈
。
8.如权利要求7所述的一种调频调幅混合网点层次连续调控制方法，其特征是对所述的扩散滤波器w的抖动算法采用如下的方法进行处理fRand＝(R(m，t)/R_MAX-0.5)×cDitherdw0＝w0-fRanddw2＝w2+fRanddw1＝w1+fRanddw3＝w3-fRand上述公式中fRand为抖动微调参数；R(m，t)为扫描当前点的随机取值参数；R_MAX为随机参数R(i)的最大值；cDither为抖动幅度调整参数，该参数决定了调幅特性质量优劣；dw0～dw3为抖动后扩散滤波器w不同方向上的权重扩散值。
9.如权利要求1、2、3、4或5所述的一种调频调幅混合网点层次连续调控制方法，其特征是步骤四中，相邻方向输出灰度动态计算输出点阵数据控制网型具体步骤如下(1)根据多位成像深度设备稳定性及输出线性化特征，确定网型控制调整密度区域；(2)根据多位成像深度设备脉冲宽度调节技术参数要求，结合不同密度区已有随机网型及设备输出稳定性对不同密度区域网点形状的要求，确定需要调整的网点形状；(3)根据步骤(2)确定的需要调整的网点形状，结合输出设备的成像特点，确定目标网点形状及应该避免出现的网点形状；(4)最终依据网点密度区及已有的输出灰度层次，采用相邻方向输出灰度动态统计算法，确定多位混合半色调二进制网点输出范围参数。
10.如权利要求9所述的一种调频调幅混合网点层次连续调控制方法，其特征是步骤(4)中所述相邻方向输出灰度动态统计算法中的控制点设定方法包括以下步骤1)设定当前点前一个点的输出灰度值a，a前一个点的输出灰度值为b；2)设定当前点前一行中的当前点位置的点的灰度值为c，c前一个点灰度值为d，c后一个点灰度值为e；3)设定输出灰度值为c的前一行相同位置的点的灰度值为f；4)设定当前点周围输出灰度统计参数Sum1＝a+c+dSum2＝a+bSum3＝c+f其中，Sum1表示当前点周围2×2矩形区域灰度加全和，Sum2表示当前点水平方向上前两个输出灰度加全和，Sum3表示当前点垂直方向上前两个输出灰度加全和。
全文摘要
本发明涉及一种在多位成像深度设备上的调频调幅混合网点层次连续调控制方法，属于图像硬拷贝复制领域的半色调网点生成技术领域。该技术在原有的多位成像深度的混合网点生成技术及网型控制等多项技术基础上，引入动态层次阶调划分的抖动控制算法，保证了原有的调频调幅多位混合网点技术具有了层次阶调划分动态调整机制，解决了因不同层次阶调所采用各异的动态网型控制手段而造成的网点灰度层次跳变，从而有效地发挥了多位成像深度设备的层次阶调再现效果，提高了半色调混合网点层次连续调输出的质量和多位深度设备输出图像的质量。
文档编号H04N1/405GK1913574SQ200610112740
公开日2007年2月14日 申请日期2006年8月31日 优先权日2006年8月31日
发明者李海峰, 杨斌 申请人:北京北大方正电子有限公司, 北京大学