CN110415175B - 一种快速去除平坦区域编码马赛克的方法 - Google Patents

一种快速去除平坦区域编码马赛克的方法 Download PDF

Info

Publication number
CN110415175B
CN110415175B CN201910716238.7A CN201910716238A CN110415175B CN 110415175 B CN110415175 B CN 110415175B CN 201910716238 A CN201910716238 A CN 201910716238A CN 110415175 B CN110415175 B CN 110415175B
Authority
CN
China
Prior art keywords
pixel
image
mosaic
input image
mask
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.)
Active
Application number
CN201910716238.7A
Other languages
English (en)
Other versions
CN110415175A (zh
Inventor
陈宇
许祝登
刘宇新
朱政
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Microframe Information Technology Co ltd
Original Assignee
Hangzhou Microframe Information Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hangzhou Microframe Information Technology Co ltd filed Critical Hangzhou Microframe Information Technology Co ltd
Priority to CN201910716238.7A priority Critical patent/CN110415175B/zh
Publication of CN110415175A publication Critical patent/CN110415175A/zh
Application granted granted Critical
Publication of CN110415175B publication Critical patent/CN110415175B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4015Image demosaicing, e.g. colour filter arrays [CFA] or Bayer patterns

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)
  • Compression Or Coding Systems Of Tv Signals (AREA)

Abstract

本发明提供一种快速去除平坦区域编码马赛克的方法,通过利用图像的局部变化特性、形态学和快速导向图滤波来实现去除平坦区域的编码马赛克,基于图像的局部特征检测出马赛克的编码块,通过形态学操作去掉孤立的马赛克,并使用快速导向图滤波算法对马赛克区域进行滤波,对滤波后的图像和原图进行图像融合获得去马赛克的图像。本发明提供的所述方法计算复杂度低,简单易实现,可以非常快速有效的去除编码马赛克。

Description

一种快速去除平坦区域编码马赛克的方法
技术领域
本发明涉及图像处理技术领域,尤其涉及去除平坦区域编码马赛克。
背景技术
马赛克主要是指画面中出现的多处类似小方块的图像,导致画面的一些区域的马赛克效应。视频图像中的马赛克严重影响视频的质量。在视频编码的过程中,基于块的编码和量化的失真会导致某些区域出现马赛克,尤其在压缩率高的情况下。现有的去除编码马赛克方法的计算复杂度高,如何快速有效地去除编码马赛克是一个很有挑战的课题
发明内容
本发明提供的快速去除平坦区域编码马赛克的方法具体包括如下步骤:
步骤1,检测输入图像I中的马赛克像素块;
步骤1.1,把所述输入图像I划分为8x8的像素块;
步骤1.2,对所述输入图像I进行8倍的下采样得到缩略图I',即计算所述输入图像的I每个8x8像素块的平均值作为所述缩略图I'的像素值;
步骤1.3,计算所述输入图像I的平方图,所述平方图为所述输入图像的每个像素的像素值进行平方运算得到的图像,对所述平方图进行比例为8倍的下采样得到所述平方图的缩略图I'2,用于快速导向图滤波中的系数计算;
步骤1.4,计算所述输入图像I的每个8x8像素块的方差σ2
步骤1.5,对于所述缩略图I'上的每个像素,如果对应的所述输入图像的8x8像素块的方差σ2小于预设阈值t1,并且所述缩略图I'上的各个像素与周围4邻域的至少一个的差值小于阈值t2,那么就将所述差值小于所述阈值t2的所述像素标记为马赛克,并将对应的蒙板图mask'上相同位置的像素记为255,否则记0;
步骤2,去除孤立的马赛克:通过形态学去除孤立的马赛克,来提高马赛克标记的精确度;
步骤2.1,对所述蒙板图mask'上的每个像素8邻域的像素平均值
Figure BDA0002154967760000021
如果所述蒙版图mask'的像素点的平均值/>
Figure BDA0002154967760000022
等于0或者255,则把所述蒙版图mask'的像素点的像素值修改为/>
Figure BDA0002154967760000023
否则保持不变;
步骤2.2,将所述蒙版图mask'通过线性插值获得所述输入图像I尺寸大小的马赛克蒙板图;
步骤3,对马赛克区域进行快速导向图滤波,通过使用快速导向图滤波方法,对所述输入图像I进行滤波处理得到滤波后的图像I1
步骤4,通过马赛克蒙板图mask,将所述输入图像I和所述滤波后的图像I1进行融合,记α=mask/255,融合图像O的计算方式如下:
O=α*I1+(1-α)*I。
步骤5,输出去马赛克后的所述融合图像O。
本发明提出的去除平坦区域编码马赛克的快速方法,通过利用图像的局部变化特性、形态学和快速导向图滤波来实现去除平坦区域的编码马赛克。所述方法计算复杂度低,简单易实现,可以非常快速有效的去除编码马赛克。
附图说明
为了更清楚地说明本说明书实施例或现有技术中的技术方案,下面将对实施例或现有技术中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本说明书中记载的一些实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本说明书实施例提供的一种去除平坦区域编码马赛克的方法的流程图;
图2是本说明书实施例提供的3x3像素块;
图3是存在编码马赛克的输入图像;
图4是滤波后的图像。
具体实施方式
为了使本技术领域的人员更好地理解本说明书中的技术方案,下面将结合本说明书一个或多个实施例中的附图,对本说明书实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本说明书一部分实施例,而不是全部的实施例。基于本说明书实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本说明书保护的范围。
以下结合附图,详细说明本说明书实施例提供的技术方案。
如图1所示,本发明提供了一种去除平坦区域编码马赛克的快速方法,包括以下步骤:
步骤1,检测输入图像I中的马赛克像素块。所述步骤1的步骤具体包括如下子步骤:
步骤1.1,把所述输入图像I划分为8x8的像素块;
步骤1.2,对所述输入图像I进行8倍的下采样得到缩略图I',即计算所述输入图像的I每个8x8像素块的平均值作为所述缩略图I'的像素值;
步骤1.3,计算所述输入图像I的平方图,所述平方图为所述输入图像的每个像素的像素值进行平方运算得到的图像,对所述平方图进行比例为8倍的下采样得到所述平方图的缩略图I'2,用于快速导向图滤波中的系数计算;
步骤1.4,计算所述输入图像I的每个8x8像素块的方差;
步骤1.5,对于所述缩略图I'上的每个像素,如果对应的所述输入图像的8x8像素块的方差σ2小于某个预设的阈值t1,并且所述缩略图I'上的各个像素与周围4邻域的至少一个的差值小于阈值t2,那么就将所述差值小于阈值t2的所述像素标记为马赛克,并将对应的蒙板图mask'上相同位置的像素记为255,否则记0,其中,如图2所示,以像素x4为例,像素x4的所述周围4邻域为x1、x3、x5和x7。所述阈值参数t1优选为3.2,所述阈值参数t2优选为2。
步骤2,去除孤立的马赛克:通过形态学去除孤立的马赛克,来提高马赛克标记的精确度。
步骤2.1,对所述蒙板图mask'上的每个像素8邻域的像素平均值。如图2所示,像素点x4的8邻域的像素平均为
Figure BDA0002154967760000041
的取值为(x0+x1+x2+x3+x5+x6+x7+x8)/8,如果蒙版图mask'的像素点(比如蒙版图mask'的边缘像素点)的8邻域像素中有不可用的像素,则计算所述8邻域像素中所有可用的像素的像素值的均值,如图2所示中x0表示一帧图像的第一个像素点p(0,0)时,所述x0的/>
Figure BDA0002154967760000042
的取值为(x1+x3+x4)/3。如果蒙版图mask'的像素点的平均值/>
Figure BDA0002154967760000043
等于0或者255,则把所述蒙版图mask'的像素点的像素值修改为/>
Figure BDA0002154967760000044
否则保持不变。即如果一个8x8的像素块被标记为马赛克,而它周围所有的8x8像素块都被标记为非马赛克,则把这个像素块标记为非马赛克;如果一个8x8的像素块被标记为了非马赛克块,而它周围的所有8x8像素块都被标记为马赛克,则把这个像素块标记为马赛克。
步骤2.2,将所述蒙版图mask'通过线性插值获得所述输入图像I尺寸大小的马赛克蒙板图。
步骤3,对马赛克区域进行快速导向图滤波。通过使用快速导向图滤波方法,对所述输入图像I进行滤波处理得到滤波后的图像I1。其中所述快速导向图滤波的半径和正则化参数的设置优选为半径r=16,正则化参数ε=52
步骤4,通过马赛克蒙板图mask,将所述输入图像I和所述滤波后的图像I1进行融合,记α=mask/255,最后融合的图像O的计算方式如下:
O=α*I1+(1-α)*I。
步骤5,输出去马赛克后的图像O。
本发明提出的去除平坦区域编码马赛克的快速方法,可以非常快速的去除马赛克,实测在i7-8700k上单线程处理1080p灰度图可以达到400fps。图3所示是存在编码马赛克的输入图像。图4所示是经过本发明提出的去除平坦区域编码马赛克的快速方法处理后的滤波后的图像。图3的右上方有一个明显的马赛克区域,经过去编码马赛克处理后,从图4上可以看出,该区域的马赛克已经去除。
本申请可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
上述具体实施方式,并不构成对本发明保护范围的限制。本领域技术人员应该明白的是,取决于设计要求和其他因素,可以发生各种各样的修改、组合、子组合和替代。任何在本发明的精神和原则之内所作的修改、等同替换和改进等,均应包含在本发明保护范围之内。

Claims (3)

1.一种快速去除平坦区域编码马赛克的方法,其特征在于,包括如下步骤:
步骤1,检测输入图像I中的马赛克像素块;
步骤1.1,把所述输入图像I划分为8x8的像素块;
步骤1.2,对所述输入图像I进行8倍的下采样得到缩略图I',即计算所述输入图像的I每个8x8像素块的平均值作为所述缩略图I'的像素值;
步骤1.3,计算所述输入图像I的平方图,所述平方图为所述输入图像的每个像素的像素值进行平方运算得到的图像,对所述平方图进行比例为8倍的下采样得到所述平方图的
缩略图I′2,用于快速导向图滤波中的系数计算;
步骤1.4,计算所述输入图像I的每个8x8像素块的方差σ2
步骤1.5,对于所述缩略图I'上的每个像素,如果对应的所述输入图像的8x8像素块的方差σ2小于预设阈值t1,并且所述缩略图I'上的各个像素与周围4邻域的至少一个的差值小
于阈值t2,则将所述差值小于所述阈值t2的所述像素标记为马赛克,并将对应的蒙版图mask'上相同位置的像素记为255,否则记0;
步骤2,去除孤立的马赛克:通过形态学去除孤立的马赛克,来提高马赛克标记的精确度;
步骤2.1,对所述蒙版图mask'上的每个像素8邻域的像素平均值
Figure QLYQS_1
,如果所述蒙版图mask'的像素点的平均值/>
Figure QLYQS_2
等于0或者255,则把所述蒙版图mask'的像素点的像素值修改为/>
Figure QLYQS_3
,否则保持不变;
步骤2.2,将所述蒙版图mask'通过线性插值获得所述输入图像I尺寸大小的马赛克蒙版图mask;
步骤3,对马赛克区域进行快速导向图滤波,通过使用快速导向图滤波方法,对所述输入图像I进行滤波处理得到滤波后的图像I1
步骤4,通过马赛克蒙版图mask,将所述输入图像I和所述滤波后的图像I1进行融合,记α=mask/255,融合图像O的计算方式如下:
Figure QLYQS_4
步骤5,输出去马赛克后的所述融合图像O。
2.根据权利要求1所述的快速去除平坦区域编码马赛克的方法,其特征在于,所述阈值t1为3 .2,所述阈值t2为2。
3.根据权利要求1所述的快速去除平坦区域编码马赛克的方法,其特征在于,所述快速导向图滤波的半径和正则化参数的设置为半径r=16,正则化参数ε=52
CN201910716238.7A 2019-08-07 2019-08-07 一种快速去除平坦区域编码马赛克的方法 Active CN110415175B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910716238.7A CN110415175B (zh) 2019-08-07 2019-08-07 一种快速去除平坦区域编码马赛克的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910716238.7A CN110415175B (zh) 2019-08-07 2019-08-07 一种快速去除平坦区域编码马赛克的方法

Publications (2)

Publication Number Publication Date
CN110415175A CN110415175A (zh) 2019-11-05
CN110415175B true CN110415175B (zh) 2023-06-02

Family

ID=68365680

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910716238.7A Active CN110415175B (zh) 2019-08-07 2019-08-07 一种快速去除平坦区域编码马赛克的方法

Country Status (1)

Country Link
CN (1) CN110415175B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112070671A (zh) * 2020-09-04 2020-12-11 平安科技(深圳)有限公司 一种基于频谱分析的马赛克去除方法、系统、终端及存储介质
CN114202486B (zh) * 2022-02-16 2022-05-20 安翰科技(武汉)股份有限公司 胶囊内窥镜图像的马赛克去除方法和系统

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07282260A (ja) * 1994-04-13 1995-10-27 Daewoo Electron Co Ltd 時系列顔画像処理による個人の認識方法
US8131110B2 (en) * 2008-07-03 2012-03-06 Seiko Epson Corporation Reducing signal overshoots and undershoots in demosaicking
CN108805884A (zh) * 2018-06-13 2018-11-13 北京搜狐新媒体信息技术有限公司 一种马赛克区域检测方法、装置及设备
CN109410123B (zh) * 2018-10-15 2023-08-18 深圳市能信安科技股份有限公司 基于深度学习的去除马赛克的方法、装置及电子设备

Also Published As

Publication number Publication date
CN110415175A (zh) 2019-11-05

Similar Documents

Publication Publication Date Title
CN109194923B (zh) 基于局部非均匀分辨率的视频图像处理设备、系统及方法
CN108694705B (zh) 一种多帧图像配准与融合去噪的方法
WO2022257759A1 (zh) 一种图像带状伪影去除方法、装置、设备和介质
WO2018082185A1 (zh) 图像处理方法和装置
CN111275626A (zh) 一种基于模糊度的视频去模糊方法、装置及设备
CN110428382B (zh) 一种用于移动终端的高效视频增强方法、装置和存储介质
EP2670142A1 (en) Image processing method, image processing device, video encoding/decoding method, video encoding/decoding device, and programs therefor
CN110415175B (zh) 一种快速去除平坦区域编码马赛克的方法
CN109688413B (zh) 以支持辅助帧的视频编码格式编码视频流的方法和编码器
CN108234826B (zh) 图像处理方法及装置
US8295647B2 (en) Compressibility-aware media retargeting with structure preserving
CN105631854A (zh) 一种基于fpga平台的自适应图像清晰度评价算法
KR102327758B1 (ko) 다중 도메인 딥러닝을 이용한 영상 내 모아레 무늬 제거 장치 및 방법
Yae et al. Inverted residual fourier transformation for lightweight single image deblurring
CN113506305A (zh) 三维点云数据的图像增强方法、语义分割方法及装置
JP6403401B2 (ja) 画像処理装置、画像処理方法、及び、プログラム
JP5042251B2 (ja) 画像処理装置および画像処理方法
CN112929562B (zh) 视频抖动的处理方法、装置、设备以及存储介质
KR101667306B1 (ko) Topological Derivative 기반 저 복잡도 영상 세그먼테이션 방법 및 시스템
US9336573B2 (en) Self-adaptive image edge correction device and method thereof
US8369637B2 (en) Image processing apparatus, image processing method, and program
CN114596210A (zh) 噪声估计方法、装置、终端设备及计算机可读存储介质
KR101650897B1 (ko) 저해상도 컨텐츠를 위한 윈도우 사이즈 주밍 방법 및 그 장치
CN103414843A (zh) 基于监控应用的视频缩放算法及实现方案
KR102337745B1 (ko) 이중 도메인 뉴럴 네트워크 기반 영상 내 모아레 무늬 제거 장치 및 방법

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