最快的3x3中值模糊
10.1国庆后,知名博主:laviewpbt http://www.cnblogs.com/Imageshop/
发起了一个优化3x3中值模糊的小活动。
俺也参加其中,今天博主laviewpbt 共享了一份不错的CLAHE代码。
free精神,真心为其点赞。
故俺也分享这份最快的3x3中值模糊的代码。
/// 编写者: laviewpbt, 编写时间: 2015.10.16, 联系QQ: 33184777 /// <summary>
/// 快速的实现3*3大小的中值模糊,边缘1像素未做处理。(2015.10.12日)
/// </summary>
/// <param name="Src">原始数据。</param>
/// <param name="Dest">目标数据。</param>
/// 百度: Fast median search_ an ANSI C implementation.html
/// 优化版本: 落羽の殇 联系QQ:200759103 static unsigned int* pixGreater = NULL;
static unsigned int* pixLess = NULL;
static unsigned int cmpTable[256 * 256 * 2] = { 0 };
class autoCmpTable
{
public:
autoCmpTable() {
unsigned int x, y, p;
unsigned int tableLength = 256 * 256;
pixGreater = cmpTable;
pixLess = cmpTable + tableLength;
for (x = 0; x < 256; x++) {
for (y = 0; y < 256; y++) {
p = x + y * 256;
if (x > y)
{
pixLess[p] = x;
pixGreater[p] = y;
}
else
{
pixGreater[p] = x;
pixLess[p] = y;
}
}
}
}
};
static autoCmpTable initCmpTable; IS_RET __stdcall Fastest33MedianBlur(TMatrix *Src, TMatrix *Dest)
{
if (Src == NULL || Dest == NULL) return IS_RET_ERR_NULLREFERENCE;
if (Src->Data == NULL || Dest->Data == NULL) return IS_RET_ERR_NULLREFERENCE;
if (Src->Width != Dest->Width || Src->Height != Dest->Height || Src->Channel != Dest->Channel || Src->Depth != Dest->Depth || Src->WidthStep != Dest->WidthStep) return IS_RET_ERR_PARAMISMATCH;
if (Src->Depth != IS_DEPTH_8U || Dest->Depth != IS_DEPTH_8U) return IS_RET_ERR_NOTSUPPORTED;
IS_RET Ret = IS_RET_OK; if (!pixGreater&&!pixLess)
{
return IS_RET_ERR_NOTSUPPORTED;
}
if (Src->Data == Dest->Data)
{
TMatrix *Clone = NULL;
Ret = IS_CloneMatrix(Src, &Clone);
if (Ret != IS_RET_OK) return Ret;
Ret = Fastest33MedianBlur(Clone, Dest);
IS_FreeMatrix(&Clone);
return Ret;
}
unsigned int X, Y, Width = Src->Width, Height = Src->Height;
unsigned char *LineP0, *LineP1, *LineP2, *LinePD;
unsigned int srcWidthStep = Src->WidthStep;
unsigned int dstWidthStep = Dest->WidthStep;
unsigned int srcChannel = Src->Channel;
unsigned int dstChannel = Dest->Channel;
unsigned char* SrcData = Src->Data;
unsigned char* DstData = Dest->Data; if (srcChannel == 1)
{
unsigned int Gray0, Gray1, Gray2, Gray3, Gray4, Gray5, Gray6, Gray7, Gray8, pos;
for (Y = 1; Y < Height - 1; Y++)
{
LineP0 = SrcData + (Y - 1) * srcWidthStep + 1;
LineP1 = LineP0 + srcWidthStep;
LineP2 = LineP1 + srcWidthStep;
LinePD = DstData + Y *dstWidthStep + 1;
for (X = 1; X < Width - 1; X++)
{
Gray0 = LineP0[X - 1];
Gray1 = LineP0[X];
Gray2 = LineP0[X + 1];
Gray3 = LineP1[X - 1];
Gray4 = LineP1[X];
Gray5 = LineP2[X + 1];
Gray6 = LineP2[X - 1];
Gray7 = LineP2[X];
Gray8 = LineP2[X + 1];
/* if (Gray1 > Gray2) Swap(Gray1, Gray2);*/
pos = Gray1 + Gray2 * 256;
Gray2 = pixLess[pos];
Gray1 = pixGreater[pos];
/* if (Gray4 > Gray5) Swap(Gray4, Gray5);*/
pos = Gray4 + Gray5 * 256;
Gray5 = pixLess[pos];
Gray4 = pixGreater[pos];
/* if (Gray7 > Gray8) Swap(Gray7, Gray8);*/
pos = Gray7 + Gray8 * 256;
Gray8 = pixLess[pos];
Gray7 = pixGreater[pos];
/* if (Gray0 > Gray1) Swap(Gray0, Gray1);*/
pos = Gray0 + Gray1 * 256;
Gray1 = pixLess[pos];
Gray0 = pixGreater[pos];
/* if (Gray3 > Gray4) Swap(Gray3, Gray4);*/
pos = Gray3 + Gray4 * 256;
Gray4 = pixLess[pos];
Gray3 = pixGreater[pos];
/* if (Gray6 > Gray7) Swap(Gray6, Gray7);*/
pos = Gray6 + Gray7 * 256;
Gray7 = pixLess[pos];
Gray6 = pixGreater[pos];
/* if (Gray1 > Gray2) Swap(Gray1, Gray2);*/
pos = Gray1 + Gray2 * 256;
Gray2 = pixLess[pos];
Gray1 = pixGreater[pos];
/* if (Gray4 > Gray5) Swap(Gray4, Gray5);*/
pos = Gray4 + Gray5 * 256;
Gray5 = pixLess[pos];
Gray4 = pixGreater[pos];
/* if (Gray7 > Gray8) Swap(Gray7, Gray8);*/
pos = Gray7 + Gray8 * 256;
Gray8 = pixLess[pos];
Gray7 = pixGreater[pos];
/* if (Gray0 > Gray3) Swap(Gray0, Gray3);*/
pos = Gray0 + Gray3 * 256;
Gray3 = pixLess[pos];
Gray0 = pixGreater[pos];
/* if (Gray5 > Gray8) Swap(Gray5, Gray8);*/
pos = Gray5 + Gray8 * 256;
Gray8 = pixLess[pos];
Gray5 = pixGreater[pos];
/* if (Gray4 > Gray7) Swap(Gray4, Gray7);*/
pos = Gray4 + Gray7 * 256;
Gray7 = pixLess[pos];
Gray4 = pixGreater[pos];
/* if (Gray3 > Gray6) Swap(Gray3, Gray6);*/
pos = Gray3 + Gray6 * 256;
Gray6 = pixLess[pos];
Gray3 = pixGreater[pos];
/* if (Gray1 > Gray4) Swap(Gray1, Gray4);*/
pos = Gray1 + Gray4 * 256;
Gray4 = pixLess[pos];
Gray1 = pixGreater[pos];
/* if (Gray2 > Gray5) Swap(Gray2, Gray5);*/
pos = Gray2 + Gray5 * 256;
Gray5 = pixLess[pos];
Gray2 = pixGreater[pos];
/* if (Gray4 > Gray7) Swap(Gray4, Gray7);*/
pos = Gray4 + Gray7 * 256;
Gray7 = pixLess[pos];
Gray4 = pixGreater[pos];
/* if (Gray4 > Gray2) Swap(Gray4, Gray2);*/
pos = Gray4 + Gray2 * 256;
Gray2 = pixLess[pos];
Gray4 = pixGreater[pos];
/* if (Gray6 > Gray4) Swap(Gray6, Gray4);*/
pos = Gray6 + Gray4 * 256;
Gray4 = pixLess[pos];
Gray6 = pixGreater[pos];
/* if (Gray4 > Gray2) Swap(Gray4, Gray2); */
pos = Gray4 + Gray2 * 256;
Gray2 = pixLess[pos];
Gray4 = pixGreater[pos]; LinePD[1] = Gray4;
LinePD++;
}
}
return Ret;
}
else
{
TMatrix *Blue = NULL, *Green = NULL, *Red = NULL, *Alpha = NULL; // 由于C变量如果不初始化,其值是随机值,可能会导致释放时的错误。
IS_RET Ret = SplitRGBA(Src, &Blue, &Green, &Red, &Alpha);
if (Ret != IS_RET_OK) goto Done24;
#pragma omp parallel num_threads(3)
{
#pragma omp sections
{
#pragma omp section
{
Ret = Fastest33MedianBlur(Blue, Blue);
}
#pragma omp section
{
Ret = Fastest33MedianBlur(Green, Green);
}
#pragma omp section
{
Ret = Fastest33MedianBlur(Red, Red);
} }
}
Ret = CombineRGBA(Dest, Blue, Green, Red, Alpha);
Done24:
IS_FreeMatrix(&Blue);
IS_FreeMatrix(&Green);
IS_FreeMatrix(&Red);
IS_FreeMatrix(&Alpha);
return Ret;
} return Ret;
}
关于交换法短值快速排序的参考资料见:https://github.com/afniHQ/AFNI/blob/ab8ca253b784ae71401927df88da2bc6e16d07c1/src/cs_qsort_small.h
最快的3x3中值模糊的更多相关文章
- 【算法随记三】小半径中值模糊的急速实现(16MB图7.5ms实现) + Photoshop中蒙尘和划痕算法解读。
在本人的博客里,分享了有关中值模糊的O(1)算法,详见:任意半径中值滤波(扩展至百分比滤波器)O(1)时间复杂度算法的原理.实现及效果 ,这里的算法的执行时间和参数是无关的.整体来说,虽然速度也很快, ...
- Emgu-WPF学习使用-中值模糊
原文:Emgu-WPF学习使用-中值模糊 实现效果: 实现途径: 前提:Image File-> System.Drawing.Bitmap->Image<Bgr, byte> ...
- OpenCV笔记(1)(图片读取与现实、色彩空间、基础运算、均值方差、逻辑运算、泛洪填充、均值中值及自定义平滑)
一.图片读取和显示 import cv2 as cv # 图片读取cv.imread(img_path) car_img = cv.imread("car1.png") # 图片显 ...
- opencv-11-中值滤波及自适应中值滤波
开始之前 在上一篇我们实现了读取噪声图像, 然后 进行三种形式的均值滤波得到结果, 由于我们自己写的均值滤波未作边缘处理, 所以效果有一定的下降, 但是总体来说, 我们得到的结果能够说明我们的算法执行 ...
- OpenCV计算机视觉学习(4)——图像平滑处理(均值滤波,高斯滤波,中值滤波,双边滤波)
如果需要处理的原图及代码,请移步小编的GitHub地址 传送门:请点击我 如果点击有误:https://github.com/LeBron-Jian/ComputerVisionPractice &q ...
- Atitit 图像处理 平滑 也称 模糊, 归一化块滤波、高斯滤波、中值滤波、双边滤波)
Atitit 图像处理 平滑 也称 模糊, 归一化块滤波.高斯滤波.中值滤波.双边滤波) 是一项简单且使用频率很高的图像处理方法 用途 去噪 去雾 各种线性滤波器对图像进行平滑处理,相关OpenC ...
- opencv3 图片模糊操作-均值滤波 高斯滤波 中值滤波 双边滤波
#include <iostream>#include <opencv2/opencv.hpp> using namespace std;using namespace cv; ...
- 学习 opencv---(8)非线性滤波:中值滤波,双边滤波
正如我们上一篇文章中讲到的,线性滤波可以实现很多种不同的图像变换.然而非线性滤波,如中值滤波器和双边滤波器,有时可以达到更好的实现效果. 邻域算子的其他一些例子还有对 二值图像进行操作的形态学算子,用 ...
- 基于FPGA的中值滤波算法实现
在这一篇开篇之前,我需要解决一个问题,上一篇我们实现了基于FPGA的均值滤波算法的实现,最后的显示效果图上发现有一些黑白色的斑点,我以为是椒盐噪声,然后在做基于FPGA的中值滤波算法的实验时,我发现黑 ...
随机推荐
- 以太坊开发DApp实战教程——用区块链、星际文件系统(IPFS)、Node.js和MongoDB来构建电商平台(一)
第一节 简介 欢迎和我们一起来用以太坊开发构建一个去中心化电商DApp!我们将用区块链.星际文件系统(IPFS).Node.js和MongoDB来构建电商平台类似淘宝的在线电商应用,卖家可以自由地出售 ...
- MyAdapter Andriod
private List<T> listdate;//定义数据对象 //为了获取item中的点击事件定义ViewHolderprivate static class ViewHolder ...
- 创建带缩进的XML
from xml.etree import ElementTree as ET from xml.dom import minidom root = ET.Element('}) son=ET.Sub ...
- iOS开发所有KeyboardType与图片对应展示
1.UIKeyboardTypeAlphabet 2.UIKeyboardTypeASCIICapable 3.UIKeyboardTypeDecimalPad 4.UIKeyboardTypeDe ...
- faster-rcnn 结构杂谈
faster-rcnn结构图: (只截取了最难理解的部分) 这个网络看似很复杂,但是理解了其中关键的层,就基本可以掌握这个结构了.要看源码!!要看源码!!要看源码 !!重要的事情说三遍. 关键的层: ...
- [知识梳理]课本3&9.1
函数:关键词:参数.返回值.函数返回类型.函数体. 函数按照返回类型,可以分为有参函数和无参函数. 函数根据是否有返回值,可以分为返回值函数和非返回值函数. 函数的定义:函数的定义可以放在任意 ...
- python入门(5)使用文件编辑器编写代码并保存执行
python入门(5)使用文件编辑器编写代码并保存执行 两款文本编辑器: 一个是Sublime Text,免费使用,但是不付费会弹出提示框: 一个是Notepad++,免费使用,有中文界面: 请注意, ...
- 关于阿里巴巴iconfont的使用方法
iconfont网址:http://www.iconfont.cn/ 说起iconfont,做前端开发的应该知道它的好处,图标库之丰富,只有你想不到的,没有你找不到的,而且轻量高清.用户在iconfo ...
- SQL Server数据库优化的10多种方法
巧妙优化sql server数据库的几种方法,在实际操作中导致查询速度慢的原因有很多,其中最为常见有以下的几种:没有索引或者没有用到索引(这是查询慢最常见的问题,是程序设计的缺陷). I/O吞吐量小, ...
- RxJava系列6(从微观角度解读RxJava源码)
RxJava系列1(简介) RxJava系列2(基本概念及使用介绍) RxJava系列3(转换操作符) RxJava系列4(过滤操作符) RxJava系列5(组合操作符) RxJava系列6(从微观角 ...