CN112308806B - 图像处理方法、装置、电子设备以及可读存储介质 - Google Patents
图像处理方法、装置、电子设备以及可读存储介质 Download PDFInfo
- Publication number
- CN112308806B CN112308806B CN202011256082.8A CN202011256082A CN112308806B CN 112308806 B CN112308806 B CN 112308806B CN 202011256082 A CN202011256082 A CN 202011256082A CN 112308806 B CN112308806 B CN 112308806B
- Authority
- CN
- China
- Prior art keywords
- pixel point
- filter window
- image processing
- image
- processed
- 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
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 25
- 238000001914 filtration Methods 0.000 claims abstract description 44
- 238000012545 processing Methods 0.000 claims abstract description 23
- 238000010276 construction Methods 0.000 claims description 4
- 230000006870 function Effects 0.000 description 11
- 238000000034 method Methods 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 7
- 238000004891 communication Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 230000002146 bilateral effect Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 238000010295 mobile communication Methods 0.000 description 4
- 238000005192 partition Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000005236 sound signal Effects 0.000 description 2
- 230000001413 cellular effect Effects 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010897 surface acoustic wave method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明提供一种图像处理方法、装置、电子设备以及存储介质。本发明通过以待处理图像的每个待处理像素点为中心像素点来构建滤波窗口，并计算滤波窗口内每个像素点与中心像素点的角度值，从而获得滤波窗口的方差，进而获取滤波窗口内每个像素点的权重系数，最终根据权重系数加权滤波窗口内的所有像素点，避免进行图像增强处理时造成梯度翻转等伪影的问题，实现滤波后的图像在降噪的同时保持良好的细节，增强图像显示质量。
Description
技术领域
本发明涉及图像处理技术领域，尤指涉及一种图像处理方法、装置、电子设备以及存储介质。
背景技术
目前，在降噪(Noise Reduction)、图像风格转换(Photographic StyleTransfer)、色调映射(Tone Mapping)、卡通演示(Cartoon Rendition)以及图像增强(Image Enhancement)领域广泛应用双边滤波器(Bilateral Filter)，虽然双边滤波器是一种采用非线性滤波方法的保边滤波器，但是它在应用于图像增强领域进行图像处理时会造成梯度翻转(Gradient Reversal)等伪影(Artifacts)。
发明内容
本发明提供一种图像处理方法、装置、电子设备以及存储介质，用以解决在图像增强领域进行图像处理时会造成梯度翻转等伪影的问题。
根据本发明的第一方面，本发明提供一种图像处理方法，所述方法包括：获取待处理图像，针对所述待处理图像中的每个待处理像素点，执行以下操作，以得到滤波处理后的输出图像；以所述待处理像素点作为中心像素点构建滤波窗口，其中，所述中心像素点为I(y,x)，所述滤波窗口的大小为M×N，M，N均为大于1的奇数；计算所述滤波窗口的每个像素点与所述中心像素点的角度值，所述滤波窗口的像素点为I(i,j)，其中i大于等于1且小于等于M，j大于等于1且小于等于N；根据所述角度值计算所述滤波窗口内的方差；根据所述方差计算所述滤波窗口内每个像素点的权重系数；以及根据所述权重系数加权所述滤波窗口的所有像素点，以获取加权结果。
在一些实施例中，所述角度值的计算公式为：
其中，theta(i,j)为所述角度值，sim(i,j)为所述滤波窗口的像素点与所述中心像素点的像素相似度，dis(i,j)为所述滤波窗口的像素点与所述中心像素点的像素距离，a为第一调整参数。
在一些实施例中，所述像素相似度的计算公式为：
sim(i,j)＝I(i,j)-I(y,x)。
在一些实施例中，所述像素距离的计算公式为：
在一些实施例中，所述方差的计算公式为：
其中，sigma为所述方差，b为预设目标均值，c为第二调整参数。
在一些实施例中，所述权重系数的计算公式为：
其中，weight(i,j)为所述权重系数。
在一些实施例中，所述加权结果的计算公式为：
其中，out(i,j)为所述加权结果。
根据本发明的第二方面，本发明提供一种图像处理装置，所述图像处理装置包括：获取模块，用于获取待处理图像，所述待处理图像中包括多个待处理像素点；构建模块，用于以所述待处理像素点作为中心像素点构建滤波窗口，其中，所述中心像素点为I(y,x)，所述滤波窗口的大小为M×N，M，N均为大于1的奇数；角度值计算模块，用于计算所述滤波窗口的每个像素点与所述中心像素点的角度值，所述滤波窗口的像素点为I(i,j)，其中i大于等于1且小于等于M，j大于等于1且小于等于N；方差计算模块，用于根据所述角度值计算所述滤波窗口内的方差；权重系数计算模块，用于根据所述方差计算所述滤波窗口内每个像素点的权重系数；以及加权模块，用于根据所述权重系数加权所述滤波窗口的所有像素点，以获取加权结果。
根据本发明的第三方面，本发明提供一种电子设备，包括处理器，存储器及存储在所述存储器上并可在所述处理器上运行的程序或指令，所述程序或指令被所述处理器执行时实现如上述的图像处理方法的步骤。
根据本发明的第四方面，本发明提供一种可读存储介质，所述可读存储介质上存储程序或指令，所述程序或指令被处理器执行时实现如上述的图像处理方法。
与现有技术相比，本发明的有益效果：通过以待处理图像的每个待处理像素点为中心像素点来构建滤波窗口，并计算滤波窗口内每个像素点与中心像素点的角度值，从而获得滤波窗口内的方差，进而获取滤波窗口内每个像素点的权重系数，最终根据权重系数加权滤波窗口内的所有像素点，避免进行图像增强处理时造成梯度翻转等伪影的问题，实现滤波后的图像在降噪的同时保持良好的细节，增强图像显示质量。
附图说明
图1为本发明实施例提供的一种图像处理方法的步骤流程示意图。
图2为本发明实施例提供的一种图像增强的框架结构示意图。
图3为本发明实施例提供的一种图像处理装置的结构示意图。
图4为本发明实施例提供的电子设备的结构示意图。
图5为本发明实施例提供的电子设备的具体结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，本发明实施例提供一种图像处理方法，应用于电子设备，该方法包括步骤S11至步骤S16。
步骤S11，获取待处理图像，所述待处理图像中包括多个待处理像素点，执行以下操作，以得到滤波处理后的输出图像。
在本发明实施例中，可以将原图片根据大小进行拆分处理，也可以直接将一张图片作为待处理图像进行处理，本发明实施例在此不做具体限定。
步骤S12，以所述待处理像素点作为中心像素点构建滤波窗口，其中，所述中心像素点为I(y,x)，所述滤波窗口的大小为M×N，M，N均为大于1的奇数。
在本发明实施例中，依次选择一个像素点作为中心像素点，以该中心像素作为滤波窗口的中心，向该中心像素的左上、正上、右上，正左、正右、左下、正下、右下八个方向拓展，其中，一对相反方向扩展的像素个数相同，选取该八个方向的像素组成M×N维矩阵作为滤波窗口。例如M＝5，N＝5时，即中心像素为I(3,3)。
步骤S13，计算所述滤波窗口的每个像素点与所述中心像素点的角度值，所述滤波窗口的像素点为I(i,j)，其中i大于等于1且小于等于M，j大于等于1且小于等于N。
在本发明实施例中，角度值的计算公式为：
其中，theta(i,j)为角度值，sim(i,j)为滤波窗口的像素点与中心像素点的像素相似度，dis(i,j)为滤波窗口的像素点与中心像素点的像素距离，a为第一调整参数，a用于对的结果进行调整，a越大，theta(i,j)的绝对值越大，滤波程度越大，其中a为大于0的实数。
像素相似度的计算公式为：
sim(i,j)＝I(i,j)-I(y,x)。
像素距离的计算公式为：
步骤S14,根据所述角度值计算所述滤波窗口内的方差。
在本发明实施例中，方差的计算公式为：
其中，sigma为所述方差，b为预设目标均值，b可以根据滤波窗口内的角度值加权平均得到，b的取值范围为c为第二调整参数，c用于对sigma的结果进行调整，c越大，sigma越大，滤波程度越大，其中c为大于0的实数。
步骤S15，根据所述方差计算所述滤波窗口内每个像素点的权重系数。
在本发明实施例中，权重系数的计算公式为：
其中，weight(i,j)为所述权重系数。
步骤S16，根据所述权重系数加权所述滤波窗口的所有像素点，以获取加权结果。
在本发明实施例中，加权结果的计算公式为：
其中，out(i,j)为所述加权结果。
通过对每个像素进行上述图像处理，以得到滤波后的图像。
如图2所示，本发明实施例提供一种图像增强的框架结构示意图。输入图像后，首先通过一个低通滤波器，得到低通分量和高频成分(High Frequency)；然后，对高频成分进行增强(Enhance)，如对高频成分乘增益值(Gain)以进行放大，最后将增强后的高频成分加回至输入图像的低通分量，从而得到输出的增强图像。
如图3所示，本发明实施例提供一种图像处理装置，应用于电子设备，该装置包括获取模块301、构建模块302、角度值计算模块303、方差计算模块304、权重系数计算模块305以及加权模块306。
获取模块301用于获取待处理图像，所述待处理图像中包括多个待处理像素点。
在本发明实施例中，可以将原图片根据大小进行拆分处理，也可以直接将一张图片作为待处理图像进行处理，本发明实施例在此不做具体限定。
构建模块302用于以所述待处理像素点作为中心像素点构建滤波窗口，其中，所述中心像素点为I(y,x)，所述滤波窗口的大小为M×N，M，N均为大于1的奇数。
在本发明实施例中，依次选择一个像素点作为中心像素点，以该中心像素作为滤波窗口的中心，向该中心像素的左上、正上、右上，正左、正右、左下、正下、右下八个方向拓展，其中，一对相反方向扩展的像素个数相同，选取该八个方向的像素组成M×N维矩阵作为滤波窗口。例如M＝5，N＝5时，即中心像素为I(3,3)。
角度值计算模块303用于计算所述滤波窗口的每个像素点与所述中心像素点的角度值，所述滤波窗口的像素点为I(i,j)，其中i大于等于1且小于等于M，j大于等于1且小于等于N。
在本发明实施例中，角度值的计算公式为：
其中，theta(i,j)为角度值，sim(i,j)为滤波窗口的像素点与中心像素点的像素相似度，dis(i,j)为滤波窗口的像素点与中心像素点的像素距离，a为第一调整参数，a用于对的结果进行调整，a越大，theta(i,j)的绝对值越大，滤波程度越大，其中a为大于0的实数。
像素相似度的计算公式为：
sim(i,j)＝I(i,j)-I(y,x)。
像素距离的计算公式为：
方差计算模块304用于根据所述角度值计算所述滤波窗口内的方差。
在本发明实施例中，方差的计算公式为：
其中，sigma为所述方差，b为预设目标均值，b可以根据滤波窗口内的角度值加权平均得到，b的取值范围为c为第二调整参数，c用于对sigma的结果进行调整，c越大，sigma越大，滤波程度越大，其中c为大于0的实数。
权重系数计算模块305用于根据所述方差计算所述滤波窗口内每个像素点的权重系数。
在本发明实施例中，权重系数的计算公式为：
其中，weight(i,j)为所述权重系数。
加权模块306用于根据所述权重系数加权所述滤波窗口的所有像素点，以获取加权结果。
在本发明实施例中，加权结果的计算公式为：
其中，out(i,j)为所述加权结果。
通过以上模块对待处理图像处理，以得到滤波后的图像。
参阅图4，本发明实施例还提供一种电子设备400，该电子设备400可以是手机、平板以及电脑等设备。如图4所示，电子设备400包括处理器401、存储器402。其中，处理器401与存储器402电性连接。
处理器401是电子设备400的控制中心，利用各种接口和线路连接整个电子设备的各个部分，通过运行或加载存储在存储器402内的应用程序，以及调用存储在存储器402内的数据，执行电子设备的各种功能和处理数据，从而对电子设备进行整体监控。
在本实施例中，该电子设备400设有多个存储分区，该多个存储分区包括系统分区和目标分区，电子设备400中的处理器401会按照如下的步骤，将一个或一个以上的应用程序的进程对应的指令加载到存储器402中，并由处理器401来运行存储在存储器402中的应用程序，从而实现各种功能：
获取待处理图像，针对所述待处理图像中的每个待处理像素点，执行以下操作，以得到滤波处理后的输出图像；
以所述待处理像素点作为中心像素点构建滤波窗口，其中，所述中心像素点为I(y,x)，所述滤波窗口的大小为M×N，M，N均为大于1的奇数；
计算所述滤波窗口的每个像素点与所述中心像素点的角度值，所述滤波窗口的像素点为I(i,j)，其中i大于等于1且小于等于M，j大于等于1且小于等于N；
根据所述角度值计算所述滤波窗口内的方差；
根据所述方差计算所述滤波窗口内每个像素点的权重系数；以及
根据所述权重系数加权所述滤波窗口的所有像素点，以获取加权结果。
参阅图5，图5示出了本发明实施例提供的电子设备500的具体结构框图，该电子设备500可以用于实施上述实施例中提供的图像处理方法。该电子设备500可以为手机或平板。电子设备500还包括以下部件。
RF电路510用于接收以及发送电磁波，实现电磁波与电信号的相互转换，从而与通讯网络或者其他设备进行通讯。RF电路510可包括各种现有的用于执行这些功能的电路元件，例如，天线、射频收发器、数字信号处理器、加密/解密芯片、用户身份模块(SIM)卡、存储器等等。RF电路510可与各种网络如互联网、企业内部网、无线网络进行通讯或者通过无线网络与其他设备进行通讯。上述的无线网络可包括蜂窝式电话网、无线局域网或者城域网。上述的无线网络可以使用各种通信标准、协议及技术，包括但并不限于全球移动通信系统(Global System for Mobile Communication,GSM)、增强型移动通信技术(Enhanced DataGSM Environment,EDGE)，宽带码分多址技术(Wideband Code Division MultipleAccess,WCDMA)，码分多址技术(Code Division Access,CDMA)、时分多址技术(TimeDivision Multiple Access,TDMA)，无线保真技术(Wireless Fidelity，Wi-Fi)(如美国电气和电子工程师协会标准IEEE 802.11a，IEEE 802.11b,IEEE802.11g和/或IEEE802.11n)、网络电话(Voice over Internet Protocol,VoIP)、全球微波互联接入(Worldwide Interoperability for Microwave Access，Wi-Max)、其他用于邮件、即时通讯及短消息的协议，以及任何其他合适的通讯协议，甚至可包括那些当前仍未被开发出来的协议。
存储器520可用于存储软件程序以及模块，如上述实施例中图像处理方法对应的程序指令/模块，处理器580通过运行存储在存储器520内的软件程序以及模块，从而执行各种功能应用以及数据处理，即实现图像处理方法的功能。存储器520可包括高速随机存储器，还可包括非易失性存储器，如一个或者多个磁性存储装置、闪存、或者其他非易失性固态存储器。在一些实例中，存储器520可进一步包括相对于处理器580远程设置的存储器，这些远程存储器可以通过网络连接至电子设备500。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
输入单元530可用于接收输入的数字或字符信息，以及产生与用户设置以及功能控制有关的键盘、鼠标、操作杆、光学或者轨迹球信号输入。具体地，输入单元530可包括触敏表面531以及其他输入设备532。触敏表面531，也称为触摸显示屏或者触控板，可收集用户在其上或附近的触摸操作(比如用户使用手指、触笔等任何适合的物体或附件在触敏表面531上或在触敏表面531附近的操作)，并根据预先设定的程式驱动相应的连接装置。可选的，触敏表面531可包括触摸检测装置和触摸控制器两个部分。其中，触摸检测装置检测用户的触摸方位，并检测触摸操作带来的信号，将信号传送给触摸控制器；触摸控制器从触摸检测装置上接收触摸信息，并将它转换成触点坐标，再送给处理器580，并能接收处理器580发来的命令并加以执行。此外，可以采用电阻式、电容式、红外线以及表面声波等多种类型实现触敏表面531。除了触敏表面531，输入单元530还可以包括其他输入设备532。具体地，其他输入设备532可以包括但不限于物理键盘、功能键(比如音量控制按键、开关按键等)、轨迹球、鼠标、操作杆等中的一种或多种。
显示单元540可用于显示由用户输入的信息或提供给用户的信息以及电子设备500的各种图形用户接口，这些图形用户接口可以由图形、文本、图标、视频和其任意组合来构成。显示单元540可包括显示面板541，可选的，可以采用LCD(Liquid Crystal Display，液晶显示器)、OLED(Organic Light-Emitting Diode，有机发光二极管)等形式来配置显示面板541。进一步的，触敏表面531可覆盖显示面板541，当触敏表面531检测到在其上或附近的触摸操作后，传送给处理器580以确定触摸事件的类型，随后处理器580根据触摸事件的类型在显示面板541上提供相应的视觉输出。虽然在图5中，触敏表面531与显示面板541是作为两个独立的部件来实现输入和输出功能，但是在某些实施例中，可以将触敏表面531与显示面板541集成而实现输入和输出功能。
电子设备500还可包括至少一种传感器550，比如光传感器、运动传感器以及其他传感器。具体地，光传感器可包括环境光传感器及接近传感器，其中，环境光传感器可根据环境光线的明暗来调节显示面板541的亮度，接近传感器可在电子设备500移动到耳边时，关闭显示面板541和/或背光。作为运动传感器的一种，重力加速度传感器可检测各个方向上(一般为三轴)加速度的大小，静止时可检测出重力的大小及方向，可用于识别手机姿态的应用(比如横竖屏切换、相关游戏、磁力计姿态校准)、振动识别相关功能(比如计步器、敲击)等；至于电子设备500还可配置的陀螺仪、气压计、湿度计、温度计、红外线传感器等其他传感器，在此不再赘述。
音频电路560、扬声器561，传声器562可提供用户与电子设备500之间的音频接口。音频电路560可将接收到的音频数据转换后的电信号，传输到扬声器561，由扬声器561转换为声音信号输出；另一方面，传声器562将收集的声音信号转换为电信号，由音频电路560接收后转换为音频数据，再将音频数据输出处理器580处理后，经RF电路510以发送给比如另一终端，或者将音频数据输出至存储器520以便进一步处理。音频电路560还可能包括耳塞插孔，以提供外设耳机与电子设备500的通信。
电子设备500通过传输模块570(例如Wi-Fi模块)可以帮助用户收发电子邮件、浏览网页和访问流式媒体等，它为用户提供了无线的宽带互联网访问。虽然图5示出了传输模块570，但是可以理解的是，其并不属于电子设备500的必须构成，完全可以根据需要在不改变发明的本质的范围内而省略。
处理器580是电子设备500的控制中心，利用各种接口和线路连接整个手机的各个部分，通过运行或执行存储在存储器520内的软件程序和/或模块，以及调用存储在存储器520内的数据，执行电子设备500的各种功能和处理数据，从而对手机进行整体监控。可选的，处理器580可包括一个或多个处理核心；在一些实施例中，处理器580可集成应用处理器和调制解调处理器，其中，应用处理器主要处理操作系统、用户界面和应用程序等，调制解调处理器主要处理无线通信。可以理解的是，上述调制解调处理器也可以不集成到处理器580中。
电子设备500还包括给各个部件供电的电源590(比如电池)，在一些实施例中，电源可以通过电源管理系统与处理器580逻辑相连，从而通过电源管理系统实现管理充电、放电、以及功耗管理等功能。电源590还可以包括一个或一个以上的直流或交流电源、再充电系统、电源故障检测电路、电源转换器或者逆变器、电源状态指示器等任意组件。
尽管未示出，电子设备500还可以包括摄像头(如前置摄像头、后置摄像头)、蓝牙模块等，在此不再赘述。具体在本实施例中，电子设备的显示单元是触摸屏显示器，电子设备还包括有存储器，以及一个或者一个以上的程序，其中一个或者一个以上程序存储于存储器中，且经配置以由一个或者一个以上处理器执行一个或者一个以上程序包含用于进行以下操作的指令：
获取待处理图像，针对所述待处理图像中的每个待处理像素点，执行以下操作，以得到滤波处理后的输出图像；
以所述待处理像素点作为中心像素点构建滤波窗口，其中，所述中心像素点为I(y,x)，所述滤波窗口的大小为M×N，M，N均为大于1的奇数；
计算所述滤波窗口的每个像素点与所述中心像素点的角度值，所述滤波窗口的像素点为I(i,j)，其中i大于等于1且小于等于M，j大于等于1且小于等于N；
根据所述角度值计算所述滤波窗口内的方差；
根据所述方差计算所述滤波窗口内每个像素点的权重系数；以及
根据所述权重系数加权所述滤波窗口的所有像素点，以获取加权结果。
具体实施时，以上各个模块可以作为独立的实体来实现，也可以进行任意组合，作为同一或若干个实体来实现，以上各个模块的具体实施可参见前面的方法实施例，在此不再赘述。
本领域普通技术人员可以理解，上述实施例的各种方法中的全部或部分步骤可以通过指令来完成，或通过指令控制相关的硬件来完成，该指令可以存储于一计算机可读的可读存储介质中，并由处理器进行加载和执行。为此，本发明实施例提供一种可读存储介质，其中存储有多条指令，该指令能够被处理器进行加载，以执行本发明实施例所提供的任一种图像处理方法中的步骤。
其中，该可读存储介质可以包括：只读存储器(ROM，Read Only Memory)、随机存取记忆体(RAM，Random Access Memory)、磁盘或光盘等。
由于该可读存储介质中所存储的指令，可以执行本发明实施例所提供的任一种图像处理方法中的步骤，因此，可以实现本发明实施例所提供的任一种图像处理方法所能实现的有益效果，详见前面的实施例，在此不再赘述。以上各个操作的具体实施可参见前面的实施例，在此不再赘述。
本发明的有益效果在于：通过以待处理图像的每个待处理像素点为中心像素点来构建滤波窗口，并计算滤波窗口内每个像素点与中心像素点的角度值，从而获得滤波窗口的方差，进而获取滤波窗口内每个像素点的权重系数，最终根据权重系数加权滤波窗口内的所有像素点，避免进行图像增强处理时造成梯度翻转等伪影的问题，实现滤波后的图像在降噪的同时保持良好的细节，增强图像显示质量。
在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述的部分，可以参见其他实施例的相关描述。
以上对本发明实施例所提供的一种图像处理方法、系统、可读存储介质及电子设备进行了详细介绍，本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的技术方案及其核心思想；本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例的技术方案的范围。
Claims (9)
1.一种图像处理方法，其特征在于，包括：
获取待处理图像，针对所述待处理图像中的每个待处理像素点，执行以下操作，以得到滤波处理后的输出图像：
以所述待处理像素点作为中心像素点构建滤波窗口，其中，所述中心像素点为I(y,x)，所述滤波窗口的大小为M×N，M，N均为大于1的奇数；
计算所述滤波窗口的每个像素点与所述中心像素点的角度值，所述滤波窗口的像素点为I(i,j)，其中i大于等于1且小于等于M，j大于等于1且小于等于N；
根据所述角度值计算所述滤波窗口内的方差；
根据所述方差计算所述滤波窗口内每个像素点的权重系数；以及
根据所述权重系数加权所述滤波窗口的所有像素点，以获取加权结果；
所述角度值的计算公式为：
其中，theta(i,j)为所述角度值，sim(i,j)为所述滤波窗口的像素点与所述中心像素点的像素相似度，dis(i,j)为所述滤波窗口的像素点与所述中心像素点的像素距离，a为第一调整参数。
2.如权利要求1所述的图像处理方法，其特征在于，所述像素相似度的计算公式为：sim(i,j)＝I(i,j)-I(y,x)。
3.如权利要求1所述的图像处理方法，其特征在于，所述像素距离的计算公式为：
4.如权利要求1所述的图像处理方法，其特征在于，所述方差的计算公式为：
其中，sigma为所述方差，theta(i,j)为角度值，b为预设目标均值，c为第二调整参数。
5.如权利要求1所述的图像处理方法，其特征在于，所述权重系数的计算公式为：
其中，weight(i,j)为所述权重系数，theta(i,j)为角度值，sigma为所述方差，b为预设目标均值。
6.如权利要求1所述的图像处理方法，其特征在于，所述加权结果的计算公式为：
其中，out(i,j)为所述加权结果，weight(i,j)为权重系数。
7.一种图像处理装置，其特征在于，包括：
获取模块，用于获取待处理图像，所述待处理图像中包括多个待处理像素点；
构建模块，用于以所述待处理像素点作为中心像素点构建滤波窗口，其中，所述中心像素点为I(y,x)，所述滤波窗口的大小为M×N，M，N均为大于1的奇数；
角度值计算模块，用于计算所述滤波窗口的每个像素点与所述中心像素点的角度值，所述滤波窗口的像素点为I(i,j)，其中i大于等于1且小于等于M，j大于等于1且小于等于N；
方差计算模块，用于根据所述角度值计算所述滤波窗口内的方差；
权重系数计算模块，用于根据所述方差计算所述滤波窗口内每个像素点的权重系数；以及
加权模块，用于根据所述权重系数加权所述滤波窗口的所有像素点，以获取加权结果；
所述角度值的计算公式为：
其中，theta(i,j)为所述角度值，sim(i,j)为所述滤波窗口的像素点与所述中心像素点的像素相似度，dis(i,j)为所述滤波窗口的像素点与所述中心像素点的像素距离，a为第一调整参数。
8.一种电子设备，其特征在于，包括处理器，存储器及存储在所述存储器上并可在所述处理器上运行的程序或指令，所述程序或指令被所述处理器执行时实现如权利要求1至6任一项所述的图像处理方法的步骤。
9.一种可读存储介质，其特征在于，所述可读存储介质上存储程序或指令，所述程序或指令被处理器执行时实现如权利要求1至6任一项所述的图像处理方法。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011256082.8A CN112308806B (zh) | 2020-11-11 | 2020-11-11 | 图像处理方法、装置、电子设备以及可读存储介质 |
| PCT/CN2021/094568 WO2022100048A1 (zh) | 2020-11-11 | 2021-05-19 | 图像处理方法、装置、电子设备以及可读存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011256082.8A CN112308806B (zh) | 2020-11-11 | 2020-11-11 | 图像处理方法、装置、电子设备以及可读存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112308806A CN112308806A (zh) | 2021-02-02 |
| CN112308806B true CN112308806B (zh) | 2024-06-18 |
Family
ID=74325846
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011256082.8A Active CN112308806B (zh) | 2020-11-11 | 2020-11-11 | 图像处理方法、装置、电子设备以及可读存储介质 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN112308806B (zh) |
| WO (1) | WO2022100048A1 (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112308806B (zh) * | 2020-11-11 | 2024-06-18 | 海宁奕斯伟集成电路设计有限公司 | 图像处理方法、装置、电子设备以及可读存储介质 |
| CN113034396B (zh) * | 2021-03-26 | 2023-05-12 | 龙芯中科(合肥)技术有限公司 | 图像滤波方法、装置、电子设备及存储介质 |
| CN114363477A (zh) * | 2021-12-30 | 2022-04-15 | 上海网达软件股份有限公司 | 基于滑窗权重回归的视频自适应锐化的方法及系统 |
| CN114972354B (zh) * | 2022-08-02 | 2022-10-28 | 济宁金筑新型建材科技有限公司 | 基于图像处理的蒸压加气混凝土砌块生产控制方法及系统 |
| CN116091405B (zh) * | 2022-11-30 | 2024-01-26 | 武汉极动智能科技有限公司 | 图像处理方法及装置、计算机设备、存储介质 |
| CN115840828B (zh) * | 2023-02-13 | 2023-06-06 | 湖北芯擎科技有限公司 | 一种图像比对显示方法、装置、设备及介质 |
| CN116071220B (zh) * | 2023-03-06 | 2023-06-20 | 浙江华感科技有限公司 | 一种图像窗口数据处理方法、装置、设备及介质 |
| CN116453029B (zh) * | 2023-06-16 | 2023-08-29 | 济南东庆软件技术有限公司 | 基于图像数据的楼宇火灾环境检测方法 |
| CN116503404B (zh) * | 2023-06-27 | 2023-09-01 | 梁山县创新工艺品股份有限公司 | 塑料玩具质量检测方法、装置、电子设备及存储介质 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104978715A (zh) * | 2015-05-11 | 2015-10-14 | 中国科学院光电技术研究所 | 一种基于滤波窗口及参数自适应的非局部均值图像去噪方法 |
Family Cites Families (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101877124A (zh) * | 2009-04-28 | 2010-11-03 | 北京捷科惠康科技有限公司 | 一种医学图像的滤波方法及系统 |
| JP5289288B2 (ja) * | 2009-11-27 | 2013-09-11 | 三菱電機株式会社 | ノイズ除去装置及びノイズ除去方法 |
| WO2013038333A2 (en) * | 2011-09-18 | 2013-03-21 | Forus Health Pvt. Ltd. | Method and system for enhancing image quality |
| CN102663708B (zh) * | 2012-04-27 | 2015-06-10 | 飞依诺科技(苏州)有限公司 | 基于方向加权中值滤波的超声图像处理方法 |
| CN103778613B (zh) * | 2014-02-21 | 2016-06-08 | 武汉大学 | 一种窗口自适应的极化sar影像滤波方法 |
| CN104331863B (zh) * | 2014-10-15 | 2018-01-09 | 中国航空工业集团公司洛阳电光设备研究所 | 一种图像滤波去噪方法 |
| US9792671B2 (en) * | 2015-12-22 | 2017-10-17 | Intel Corporation | Code filters for coded light depth acquisition in depth images |
| CN105894464B (zh) * | 2016-03-28 | 2018-08-31 | 福州瑞芯微电子股份有限公司 | 一种中值滤波图像处理方法和装置 |
| CN106373095B (zh) * | 2016-08-29 | 2017-11-17 | 广东欧珀移动通信有限公司 | 一种图像处理方法及终端 |
| CN106874912A (zh) * | 2016-12-20 | 2017-06-20 | 银江股份有限公司 | 一种基于改进lbp算子的图像目标检测方法 |
| CN108230255A (zh) * | 2017-09-19 | 2018-06-29 | 北京市商汤科技开发有限公司 | 用于实现图像增强的方法、装置和电子设备 |
| CN109087257B (zh) * | 2018-07-25 | 2020-08-07 | 武汉科技大学 | 一种基于参数估计框架的空域增量图像滤波方法 |
| CN109389560B (zh) * | 2018-09-27 | 2022-07-01 | 深圳开阳电子股份有限公司 | 一种自适应加权滤波图像降噪方法、装置及图像处理设备 |
| CN109741287B (zh) * | 2018-12-27 | 2021-01-01 | 湖南国科微电子股份有限公司 | 图像导向滤波方法及装置 |
| CN110211082B (zh) * | 2019-05-31 | 2021-09-21 | 浙江大华技术股份有限公司 | 一种图像融合方法、装置、电子设备及存储介质 |
| CN110619647B (zh) * | 2019-09-16 | 2022-12-30 | 中山大学 | 基于边缘点频域空域特征结合图像模糊区域定位方法 |
| CN112308806B (zh) * | 2020-11-11 | 2024-06-18 | 海宁奕斯伟集成电路设计有限公司 | 图像处理方法、装置、电子设备以及可读存储介质 |
-
2020
- 2020-11-11 CN CN202011256082.8A patent/CN112308806B/zh active Active
-
2021
- 2021-05-19 WO PCT/CN2021/094568 patent/WO2022100048A1/zh active Application Filing
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104978715A (zh) * | 2015-05-11 | 2015-10-14 | 中国科学院光电技术研究所 | 一种基于滤波窗口及参数自适应的非局部均值图像去噪方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2022100048A1 (zh) | 2022-05-19 |
| CN112308806A (zh) | 2021-02-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112308806B (zh) | 图像处理方法、装置、电子设备以及可读存储介质 | |
| US11363196B2 (en) | Image selection method and related product | |
| CN112532958B (zh) | 图像处理方法、装置、电子设备以及可读存储介质 | |
| CN107749046B (zh) | 一种图像处理方法及移动终端 | |
| CN112703530B (zh) | 图像处理方法、装置、计算机设备和存储介质 | |
| US20220319378A1 (en) | Screen color temperature control method, apparatus, storage medium, and mobile terminal | |
| CN111182236A (zh) | 一种图像合成方法、装置、存储介质及终端设备 | |
| JP6862564B2 (ja) | 画像合成のための方法、装置および不揮発性コンピュータ可読媒体 | |
| CN112330564B (zh) | 图像处理方法、装置、电子设备以及可读存储介质 | |
| CN108259808B (zh) | 一种视频帧压缩方法和移动终端 | |
| CN111556248B (zh) | 拍摄方法、装置、存储介质及移动终端 | |
| CN110855896A (zh) | 拍摄方法及装置、存储介质、终端设备 | |
| CN117116203A (zh) | 屏幕亮度调整方法、装置、存储介质及电子设备 | |
| CN107742279B (zh) | 一种图像处理方法、装置及存储介质 | |
| CN110706158A (zh) | 图像处理方法、图像处理装置及终端设备 | |
| CN116092434B (zh) | 调光方法、装置、电子设备及计算机可读存储介质 | |
| CN114125151B (zh) | 图像处理方法、移动终端及存储介质 | |
| CN112468725B (zh) | 照片拍摄方法、装置、存储介质及移动终端 | |
| CN108280816B (zh) | 一种高斯滤波方法及移动终端 | |
| CN113763230B (zh) | 图像风格迁移模型训练方法、风格迁移方法及装置 | |
| CN110990606B (zh) | 图片存储方法、装置、存储介质及电子设备 | |
| CN111343335B (zh) | 图像显示处理方法、系统、存储介质及移动终端 | |
| CN117392035A (zh) | 图像处理方法、装置、电子设备及计算机可读存储介质 | |
| CN109165173B (zh) | 一种握手处理方法及移动终端 | |
| CN107682631B (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 | ||
| CB02 | Change of applicant information |
Address after: Room 263, block B, science and technology innovation center, 128 Shuanglian Road, Haining Economic Development Zone, Haining City, Jiaxing City, Zhejiang Province, 314400 Applicant after: Haining yisiwei IC Design Co.,Ltd. Applicant after: Beijing ESWIN Computing Technology Co.,Ltd. Address before: Room 263, block B, science and technology innovation center, 128 Shuanglian Road, Haining Economic Development Zone, Haining City, Jiaxing City, Zhejiang Province, 314400 Applicant before: Haining yisiwei IC Design Co.,Ltd. Applicant before: Beijing yisiwei Computing Technology Co.,Ltd. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant |