GDI+有Bitmap类。
应用例子:
有关GDI+对bmp的处理
数字图像处理算法实现
[ 作者:admin | 转贴自:本站原创 | 点击数:634 | 更新时间:2004-12-31 | 文章录入:admin ]
摘要: 关于空间域图像处理算法框架,直方图处理,空间域滤波器算法框架的编程心得,使用GDI+(C++)
一,图像文件的读取
初学数字图像处理时,图像文件的读取往往是一件麻烦的事情,我们要面对各种各样的图像文件格式,如果仅用C++的fstream库那就必须了解各种图像编码格式,这对于初学图像处理是不太现实的,需要一个能帮助轻松读取各类图像文件的库。在Win32平台上GDI+(C++)是不错的选择,不光使用上相对于Win32 GDI要容易得多,而且也容易移植到.Net平台上的GDI+。
Gdiplus::Bitmap类为我们提供了读取各类图像文件的接口,Bitmap::LockBits方法产生的BitmapData类也为我们提供了高速访问图像文件流的途径。这样我们就可以将精力集中于图像处理算法的实现,而不用关心各种图像编码。具体使用方式请参考MSDN中GDI+文档中关于Bitmap类和BitmapData类的说明。另外GDI+仅在Windows XP/2003上获得直接支持,对于Windows 2000必须安装相关DLL,或者安装有Office 2003,Visual Studio 2003 .Net等软件。
二,空间域图像处理算法框架
(1) 在空间域图像处理中,对于一个图像我们往往需要对其逐个像素的进行处理,对每个像素的处理使用相同的算法(或者是图像中的某个矩形部分)。即,对于图像f(x,y),其中0≤x≤M,0≤y≤N,图像为M*N大小,使用算法algo,则f(x,y) = algo(f(x,y))。事先实现一个算法框架,然后再以函数指针或函数对象(functor,即实现operator()的对象)传入算法,可以减轻编程的工作量。
如下代码便是一例:
#ifndef PROCESSALGO_H
#define PROCESSALGO_H
#include <windows.h>
#include <Gdiplus.h>
namespace nsimgtk
{
template <typename pixelType, Gdiplus::PixelFormat pixelFormat, class Processor>
bool ProcessPixelsOneByOne(Gdiplus::Bitmap* const p_bitmap, Processor processor, unsigned int x, unsigned int y,
unsigned int width, unsigned int height)
{
if (p_bitmap == NULL)
{
return false;
}
if ((width + x > p_bitmap->GetWidth()) || (height + y >p_bitmap->GetHeight()))
{
return false;
}
Gdiplus::BitmapData bitmapData;
Gdiplus::Rect rect(x, y, width,height);
if (p_bitmap->LockBits(&rect, Gdiplus::ImageLockModeWrite, pixelFormat, &bitmapData) != Gdiplus::Ok)
{
return false;
}
pixelType *pixels = (pixelType*)bitmapData.Scan0;
for (unsigned int row=0; row<height; ++row)
{
for (unsigned int col=0; col<width; ++col)
{
processor(&pixels[col+row*bitmapData.Stride/sizeof(pixelType)]);
}
}
if (p_bitmap->UnlockBits(&bitmapData) != Gdiplus::Ok)
{
return false;
}
return true;
}
}
#endif
ProcessPixelsOneByOne函数可以对图像中从(x,y)位置起始,width*height大小的区域进行处理。模板参数pixelType用于指定像素大小,例如在Win32平台上传入unsigned char即为8位,用于8阶灰度图。模板参数Processor为图像处理算法实现,可以定义类实现void operator(pixelType *)函数,或者传入同样接口的函数指针。
如下便是一些算法示例(说明见具体注释):
#ifndef SPATIALDOMAIN_H
#define SPATIALDOMAIN_H
#include <cmath>
#include <string>
namespace nsimgtk
{
// 8阶灰度图的灰度反转算法
class NegativeGray8
{
public:
void operator()(unsigned char *const p_value)
{
*p_value ^= 0xff;
}
};
// 8阶灰度图的Gamma校正算法
class GammaCorrectGray8
{
private:
unsigned char d_s[256];
public:
GammaCorrectGray8::GammaCorrectGray8(double c, double gamma);
void operator()(unsigned char*const p_value)
{
*p_value = d_s[*p_value];
}
};
// 8阶灰度图的饱和度拉伸算法
class ContrastStretchingGray8
{
private:
unsigned char d_s[256];
public:
ContrastStretchingGray8::ContrastStretchingGray8(double a1, double b1, unsigned int x1,
double a2, double b2, unsigned int x2, double a3, double b3);
void operator()(unsigned char*const p_value)
{
*p_value = d_s[*p_value];
}
};
// 8阶灰度图的位平面分割,构造函数指定位平面号
class BitPlaneSliceGray8
{
private:
unsigned char d_s[256];
public:
BitPlaneSliceGray8(unsigned char bitPlaneNum);
void operator()(unsigned char* const p_value)
{
*p_value = d_s[*p_value];
}
};
}
#endif
// 上述类中各构造函数的实现代码,应该分在另一个文件中,此处为说明方便,一并列出
#include "SpatialDomain/spatialDomain.h"
namespace nsimgtk
{
GammaCorrectGray8::GammaCorrectGray8(double c, double gamma)
{
double temp;
for (unsigned int i=0; i<256; ++i)
{
temp = ceil(c * 255.0 * pow(double(i)/255.0, gamma));
d_s[i] = unsigned char(temp);
}
}
ContrastStretchingGray8::ContrastStretchingGray8(double a1, double b1, unsigned int x1,
double a2, double b2, unsigned int x2, double a3, double b3)
{
if (x1 > 255 || x2 > 255 || x1 > x1)
{
for (unsigned int i=0; i<256; ++i)
d_s[i] = i;
}
else
{
double tmp;
for (unsigned int i=0; i<x1; ++i)
{
tmp = ceil(a1*double(i)+b1);
d_s[i] = (unsigned char)tmp;
}
for (unsigned int i=x1; i<x2; ++i)
{
tmp = ceil(a2*double(i)+b2);
d_s[i] = (unsigned char)tmp;
}
for (unsigned int i=x2; i<256; ++i)
{
tmp = ceil(a3*double(i)+b3);
d_s[i] = (unsigned char)tmp;
}
}
}
BitPlaneSliceGray8::BitPlaneSliceGray8(unsigned char bitPlaneNum)
{
unsigned char bitMaskArray[8] =
{
0x01, 0x02, 0x04, 0x08,
0x10, 0x20, 0x40, 0x80
};
for (unsigned int i=0; i<256; ++i)
{
unsigned char tmp = i;
tmp &= bitMaskArray[bitPlaneNum];
tmp = (tmp >> bitPlaneNum) * 255;
d_s[i] = tmp;
}
}
}
(2) 直方图在GDI+1.0中没有获得支持,我们必须自行实现。直方图相关的处理在数字图像处理中占有重要地位,可以通过它获取图像灰度级的统计信息,且可以通过直方图进行一些重要的图像增强技术,如直方图均衡化,直方图规定化,基本全局门限等。
下面是获取8阶图像直方图的算法实现:
namespace nsimgtk
{
bool GetHistogramNormalizeGray8(Gdiplus::Bitmap * const p_bitmap, float *histogramArray)
{
if (p_bitmap == NULL || histogramArray == NULL)
{
return false;
}
Gdiplus::BitmapData bitmapData;
Gdiplus::Rect rect(0, 0, p_bitmap->GetWidth(), p_bitmap->GetHeight());
if (p_bitmap->LockBits(&rect, Gdiplus::ImageLockModeRead, PixelFormat8bppIndexed, &bitmapData) != Gdiplus::Ok)
{
return false;
}
unsigned char *pixels = (unsigned char*)bitmapData.Scan0;
unsigned int histogram[256];
for (int i=0; i<256; ++i)
{
histogram[i] = 0;
}
for (unsigned int row=0; row<p_bitmap->GetWidth(); ++row)
{
for (unsigned int col=0; col<p_bitmap->GetHeight(); ++col)
{
++histogram[pixels[col+row*bitmapData.Stride]];
}
}
const unsigned int totalPixels = p_bitmap->GetWidth() * p_bitmap->GetHeight();
for (int i=0; i<256; ++i)
{
histogramArray[i] = float(histogram[i]) / float(totalPixels);
}
if (p_bitmap->UnlockBits(&bitmapData) != Gdiplus::Ok)
{
return false;
}
return true;
}
}
在获取直方图后(即上面算法的第二个参数),再将其作为参数传入下面的对象的构造函数,然后以该对象为仿函数传入ProcessPixelsOneByOne即可实现8阶图像直方图均衡化:
#ifndef SPATIALDOMAIN_H
#define SPATIALDOMAIN_H
#include <cmath>
#include <string>
namespace nsimgtk
{
// 8阶灰度图的直方图均衡化
class HistogramEqualizationGray8
{
private:
unsigned char d_s[256];
public:
HistogramEqualizationGray8(const float *const histogramArray);
void operator()(unsigned char *const p_value)
{
*p_value = d_s[*p_value];
}
};
}
#endif
//
#include "SpatialDomain/spatialDomain.h"
namespace nsimgtk
{
HistogramEqualizationGray8::HistogramEqualizationGray8(const float *const histogramArray)
{
if (histogramArray != NULL)
{
float sum = 0.0;
for (int i=0; i<256; ++i)
{
sum += histogramArray[i];
d_s[i] = unsigned char(sum * 255);
}
}
}
}
(3)空间域滤波器,滤波器是一个m*n大小的掩模,其中m,n均为大于1的奇数。滤波器逐像素地通过图像的全部或部分矩形区域,然后逐像素地对掩模覆盖下的像素使用滤波器算法获得响应,将响应赋值于当前像素即掩模中心像素,另外滤波器算法使用中将会涉及到图像边缘的问题,这可以对边缘部分掩模使用补零法补齐掩模下无像素值的区域,或者掩模的移动范围以不越出图像边缘的方式移动,当然这些处理方法都会给图像边缘部分带来不良效果,但是一般情况下,图像边缘部分往往不是我们关注的部分或者没有重要的信息。
下面的滤波器算法框架SpatialFilterAlgo即以补零法(zero-padding)实现:
#ifndef SPATIALFILTER_H
#define SPATIALFILTER_H
#include <vector>
#include <numeric>
#include <algorithm>
#include <gdiplus.h>
#include <fstream>
#include <cmath>
namespace nsimgtk
{
template <typename pixelType, Gdiplus::PixelFormat pixelFormat, class FilterMask>
bool SpatialFilterAlgo(Gdiplus::Bitmap* const p_bitmap, FilterMask filterMask, unsigned int x, unsigned int y,
unsigned int width, unsigned int height)
{
if (p_bitmap == NULL)
{
return false;
}
if ((width + x > p_bitmap->GetWidth()) || (height + y >p_bitmap->GetHeight()))
{
return false;
}
Gdiplus::BitmapData bitmapData;
Gdiplus::Rect rect(x, y, width,height);
if (p_bitmap->LockBits(&rect, Gdiplus::ImageLockModeWrite, pixelFormat, &bitmapData) != Gdiplus::Ok)
{
return false;
}
pixelType *pixels = (pixelType*)bitmapData.Scan0;
const unsigned int m = filterMask.d_m; // mask's width
const unsigned int n = filterMask.d_n; // mask's height
std::vector<pixelType> tmpImage((m-1+width)*(n-1+height)); // extend image to use zero-padding
// copy original bitmap to extended image with zero-padding method
for (unsigned int row=0; row<height; ++row)
{
for (unsigned int col=0; col<width; ++col)
{
tmpImage[(col+m/2)+(row+n/2)*(bitmapData.Stride/sizeof(pixelType)+m-1)] =
pixels[col+row*bitmapData.Stride/sizeof(pixelType)];
}
}
// process every pixel with filterMask
for (unsigned int row=0; row<height; ++row)
{
for (unsigned int col=0; col<width; ++col)
{
// fill the "m*n" mask with the current pixel's neighborhood
for (unsigned int i=0; i<n; ++i)
{
for (unsigned int j=0; j<m; ++j)
{
filterMask.d_mask[i*m+j] = tmpImage[(col+j)+(row+i)*(bitmapData.Stride/sizeof(pixelType)+m-1)];
}
}
// replace the current pixel with filter mask's response
pixels[col+row*bitmapData.Stride/sizeof(pixelType)] = filterMask.response();
}
}
if (p_bitmap->UnlockBits(&bitmapData) != Gdiplus::Ok)
{
return false;
}
return true;
}
}
#endif
其中模板参数FilterMask即为滤波掩模算法。通常的滤波算法有均值滤波器,可以模糊化图像,去除图形中的细节部分,使得我们可以关注图像中较为明显的部分,均值滤波器用于周期性噪声。中值滤波器用于图像中存在椒盐噪声也即脉冲噪声的情况下。另外有基于一阶微分的Sobel梯度算子和基于两阶微分的拉普拉斯算子,它们往往被用于边缘检测中。
下面是一些滤波器算法的具体实现,所以滤波器算法都应该实现pixelType response()函数以及有d_mask,d_m,d_n成员,这可以通过继承__filteMask类获得(不需要付出虚函数代价)。
#ifndef SPATIALFILTER_H
#define SPATIALFILTER_H
#include <vector>
#include <numeric>
#include <algorithm>
#include <gdiplus.h>
#include <fstream>
#include <cmath>
namespace nsimgtk
{
// 滤波器掩模的基类,提供掩模大小d_m,d_n,掩模覆盖下的m*n个像素值d_mask
// others filterMask should inherit it
template <typename pixelType>
struct __filterMask
{
const unsigned int d_m;
const unsigned int d_n;
// image's pixels under the m*n filter mask
std::vector<pixelType> d_mask;
// filter mask's width and heigh must be a odd, if not, it will plus one for the width or the height
__filterMask(unsigned int m, unsigned int n)
: d_m(m%2 ? m:m+1), d_n(n%2 ? n:n+1), d_mask(d_m*d_n)
{
}
};
// 掩模权值为全1的均值滤波器
template <typename pixelType>
class averagingFilterMaskSp
: public __filterMask<pixelType>
{
public:
averagingFilterMaskSp(unsigned int m, unsigned int n)
: __filterMask<pixelType>(m, n)
{ }
pixelType response()
{
return std::accumulate(d_mask.begin(), d_mask.end(), 0) / (d_m * d_n);
}
};
// 可自定义掩模权值的均值滤波器
template <typename pixelType>
class averagingFilterMask
: public __filterMask<pixelType>
{
private:
std::vector<pixelType> d_weight; // weights' vector(m*n)
int d_weight_sum; // all weights' sum
public:
averagingFilterMask(unsigned int m, unsigned int n, const std::vector<pixelType>& weightVec)
: __filterMask<pixelType>(m, n), d_weight(weightVec)
{
if (weightVec.size() != d_mask.size())
{
// if weight's size isn't equal to mask's size, it will change filter mask as a special filter mask
d_weight.resize(d_mask.size(), 1);
}
d_weight_sum = std::accumulate(d_weight.begin(), d_weight.end(), 0);
}
pixelType response()
{
return std::inner_product(d_mask.begin(), d_mask.end(), d_weight.begin(), 0) / d_weight_sum;
}
};
// 中值滤波器
template <typename pixelType>
class medianFilterMask
: public __filterMask<pixelType>
{
public:
medianFilterMask(unsigned int m, unsigned int n)
: __filterMask<pixelType>(m, n)
{ }
pixelType response()
{
std::sort(d_mask.begin(), d_mask.end());
return d_mask[d_mask.size()/2];
}
};
// 3*3拉普拉斯滤波器
// the mask is: [0 1 0 [0 -1 0
// 1 -5 1 or -1 5 -1
// 0 1 0] 0 -1 0]
// if pixel's brightness is less than min, set it to min
// if pixel's brightness is larger than max, set it to max
template <typename pixelType, pixelType min, pixelType max>
class laplacianFilter
: public __filterMask<pixelType>
{
public:
laplacianFilter()
: __filterMask<pixelType>(3, 3)
{ }
pixelType response()
{
int ret = (int)(5*(int)d_mask[4]) - ((int)d_mask[5]+d_mask[3]+d_mask[1]+d_mask[7]);
if (ret < min)
ret = min;
if (ret > max)
ret = max;
return ret;
}
};
// 3*3Sobel滤波器
// the mask is: [-1 -2 -1 [-1 0 1
// 0 0 0 and -2 0 2
// 1 2 1] -1 0 1]
// if pixel's brightness is larger than max, set it to max
template <typename pixelType, pixelType max>
class sobelFilter
: public __filterMask<pixelType>
{
public:
sobelFilter()
: __filterMask<pixelType>(3, 3)
{ }
pixelType response()
{
int ret = ::abs(d_mask[6]+2*d_mask[7]+d_mask[8]-d_mask[0]-2*d_mask[1]-d_mask[2])
+ ::abs(d_mask[2]+2*d_mask[5]+d_mask[8]-d_mask[0]-2*d_mask[3]-d_mask[6]);
if (ret > max)
ret = max;
return ret;
}
};
}
#endif
数字图像处理算法实现
------------编程心得(1)
2001414班 朱伟 20014123
摘要: 关于空间域图像处理算法框架,直方图处理,空间域滤波器算法框架的编程心得,使用GDI+(C++)
一,图像文件的读取
初学数字图像处理时,图像文件的读取往往是一件麻烦的事情,我们要面对各种各样的图像文件格式,如果仅用C++的fstream库那就必须了解各种图像编码格式,这对于初学图像处理是不太现实的,需要一个能帮助轻松读取各类图像文件的库。在Win32平台上GDI+(C++)是不错的选择,不光使用上相对于Win32 GDI要容易得多,而且也容易移植到.Net平台上的GDI+。
Gdiplus::Bitmap类为我们提供了读取各类图像文件的接口,Bitmap::LockBits方法产生的BitmapData类也为我们提供了高速访问图像文件流的途径。这样我们就可以将精力集中于图像处理算法的实现,而不用关心各种图像编码。具体使用方式请参考MSDN中GDI+文档中关于Bitmap类和BitmapData类的说明。另外GDI+仅在Windows XP/2003上获得直接支持,对于Windows 2000必须安装相关DLL,或者安装有Office 2003,Visual Studio 2003 .Net等软件。
二,空间域图像处理算法框架
(1) 在空间域图像处理中,对于一个图像我们往往需要对其逐个像素的进行处理,对每个像素的处理使用相同的算法(或者是图像中的某个矩形部分)。即,对于图像f(x,y),其中0≤x≤M,0≤y≤N,图像为M*N大小,使用算法algo,则f(x,y) = algo(f(x,y))。事先实现一个算法框架,然后再以函数指针或函数对象(functor,即实现operator()的对象)传入算法,可以减轻编程的工作量。
如下代码便是一例:
#ifndef PROCESSALGO_H
#define PROCESSALGO_H
#include <windows.h>
#include <Gdiplus.h>
namespace nsimgtk
{
template <typename pixelType, Gdiplus::PixelFormat pixelFormat, class Processor>
bool ProcessPixelsOneByOne(Gdiplus::Bitmap* const p_bitmap, Processor processor, unsigned int x, unsigned int y,
unsigned int width, unsigned int height)
{
if (p_bitmap == NULL)
{
return false;
}
if ((width + x > p_bitmap->GetWidth()) || (height + y >p_bitmap->GetHeight()))
{
return false;
}
Gdiplus::BitmapData bitmapData;
Gdiplus::Rect rect(x, y, width,height);
if (p_bitmap->LockBits(&rect, Gdiplus::ImageLockModeWrite, pixelFormat, &bitmapData) != Gdiplus::Ok)
{
return false;
}
pixelType *pixels = (pixelType*)bitmapData.Scan0;
for (unsigned int row=0; row<height; ++row)
{
for (unsigned int col=0; col<width; ++col)
{
processor(&pixels[col+row*bitmapData.Stride/sizeof(pixelType)]);
}
}
if (p_bitmap->UnlockBits(&bitmapData) != Gdiplus::Ok)
{
return false;
}
return true;
}
}
#endif
ProcessPixelsOneByOne函数可以对图像中从(x,y)位置起始,width*height大小的区域进行处理。模板参数pixelType用于指定像素大小,例如在Win32平台上传入unsigned char即为8位,用于8阶灰度图。模板参数Processor为图像处理算法实现,可以定义类实现void operator(pixelType *)函数,或者传入同样接口的函数指针。
如下便是一些算法示例(说明见具体注释):
#ifndef SPATIALDOMAIN_H
#define SPATIALDOMAIN_H
#include <cmath>
#include <string>
namespace nsimgtk
{
// 8阶灰度图的灰度反转算法
class NegativeGray8
{
public:
void operator()(unsigned char *const p_value)
{
*p_value ^= 0xff;
}
};
// 8阶灰度图的Gamma校正算法
class GammaCorrectGray8
{
private:
unsigned char d_s[256];
public:
GammaCorrectGray8::GammaCorrectGray8(double c, double gamma);
void operator()(unsigned char*const p_value)
{
*p_value = d_s[*p_value];
}
};
// 8阶灰度图的饱和度拉伸算法
class ContrastStretchingGray8
{
private:
unsigned char d_s[256];
public:
ContrastStretchingGray8::ContrastStretchingGray8(double a1, double b1, unsigned int x1,
double a2, double b2, unsigned int x2, double a3, double b3);
void operator()(unsigned char*const p_value)
{
*p_value = d_s[*p_value];
}
};
// 8阶灰度图的位平面分割,构造函数指定位平面号
class BitPlaneSliceGray8
{
private:
unsigned char d_s[256];
public:
BitPlaneSliceGray8(unsigned char bitPlaneNum);
void operator()(unsigned char* const p_value)
{
*p_value = d_s[*p_value];
}
};
}
#endif
// 上述类中各构造函数的实现代码,应该分在另一个文件中,此处为说明方便,一并列出
#include "SpatialDomain/spatialDomain.h"
namespace nsimgtk
{
GammaCorrectGray8::GammaCorrectGray8(double c, double gamma)
{
double temp;
for (unsigned int i=0; i<256; ++i)
{
temp = ceil(c * 255.0 * pow(double(i)/255.0, gamma));
d_s[i] = unsigned char(temp);
}
}
ContrastStretchingGray8::ContrastStretchingGray8(double a1, double b1, unsigned int x1,
double a2, double b2, unsigned int x2, double a3, double b3)
{
if (x1 > 255 || x2 > 255 || x1 > x1)
{
for (unsigned int i=0; i<256; ++i)
d_s[i] = i;
}
else
{
double tmp;
for (unsigned int i=0; i<x1; ++i)
{
tmp = ceil(a1*double(i)+b1);
d_s[i] = (unsigned char)tmp;
}
for (unsigned int i=x1; i<x2; ++i)
{
tmp = ceil(a2*double(i)+b2);
d_s[i] = (unsigned char)tmp;
}
for (unsigned int i=x2; i<256; ++i)
{
tmp = ceil(a3*double(i)+b3);
d_s[i] = (unsigned char)tmp;
}
}
}
BitPlaneSliceGray8::BitPlaneSliceGray8(unsigned char bitPlaneNum)
{
unsigned char bitMaskArray[8] =
{
0x01, 0x02, 0x04, 0x08,
0x10, 0x20, 0x40, 0x80
};
for (unsigned int i=0; i<256; ++i)
{
unsigned char tmp = i;
tmp &= bitMaskArray[bitPlaneNum];
tmp = (tmp >> bitPlaneNum) * 255;
d_s[i] = tmp;
}
}
}
(2) 直方图在GDI+1.0中没有获得支持,我们必须自行实现。直方图相关的处理在数字图像处理中占有重要地位,可以通过它获取图像灰度级的统计信息,且可以通过直方图进行一些重要的图像增强技术,如直方图均衡化,直方图规定化,基本全局门限等。
下面是获取8阶图像直方图的算法实现:
namespace nsimgtk
{
bool GetHistogramNormalizeGray8(Gdiplus::Bitmap * const p_bitmap, float *histogramArray)
{
if (p_bitmap == NULL || histogramArray == NULL)
{
return false;
}
Gdiplus::BitmapData bitmapData;
Gdiplus::Rect rect(0, 0, p_bitmap->GetWidth(), p_bitmap->GetHeight());
if (p_bitmap->LockBits(&rect, Gdiplus::ImageLockModeRead, PixelFormat8bppIndexed, &bitmapData) != Gdiplus::Ok)
{
return false;
}
unsigned char *pixels = (unsigned char*)bitmapData.Scan0;
unsigned int histogram[256];
for (int i=0; i<256; ++i)
{
histogram[i] = 0;
}
for (unsigned int row=0; row<p_bitmap->GetWidth(); ++row)
{
for (unsigned int col=0; col<p_bitmap->GetHeight(); ++col)
{
++histogram[pixels[col+row*bitmapData.Stride]];
}
}
const unsigned int totalPixels = p_bitmap->GetWidth() * p_bitmap->GetHeight();
for (int i=0; i<256; ++i)
{
histogramArray[i] = float(histogram[i]) / float(totalPixels);
}
if (p_bitmap->UnlockBits(&bitmapData) != Gdiplus::Ok)
{
return false;
}
return true;
}
}
在获取直方图后(即上面算法的第二个参数),再将其作为参数传入下面的对象的构造函数,然后以该对象为仿函数传入ProcessPixelsOneByOne即可实现8阶图像直方图均衡化:
#ifndef SPATIALDOMAIN_H
#define SPATIALDOMAIN_H
#include <cmath>
#include <string>
namespace nsimgtk
{
// 8阶灰度图的直方图均衡化
class HistogramEqualizationGray8
{
private:
unsigned char d_s[256];
public:
HistogramEqualizationGray8(const float *const histogramArray);
void operator()(unsigned char *const p_value)
{
*p_value = d_s[*p_value];
}
};
}
#endif
//
#include "SpatialDomain/spatialDomain.h"
namespace nsimgtk
{
HistogramEqualizationGray8::HistogramEqualizationGray8(const float *const histogramArray)
{
if (histogramArray != NULL)
{
float sum = 0.0;
for (int i=0; i<256; ++i)
{
sum += histogramArray[i];
d_s[i] = unsigned char(sum * 255);
}
}
}
}
(3)空间域滤波器,滤波器是一个m*n大小的掩模,其中m,n均为大于1的奇数。滤波器逐像素地通过图像的全部或部分矩形区域,然后逐像素地对掩模覆盖下的像素使用滤波器算法获得响应,将响应赋值于当前像素即掩模中心像素,另外滤波器算法使用中将会涉及到图像边缘的问题,这可以对边缘部分掩模使用补零法补齐掩模下无像素值的区域,或者掩模的移动范围以不越出图像边缘的方式移动,当然这些处理方法都会给图像边缘部分带来不良效果,但是一般情况下,图像边缘部分往往不是我们关注的部分或者没有重要的信息。
下面的滤波器算法框架SpatialFilterAlgo即以补零法(zero-padding)实现:
#ifndef SPATIALFILTER_H
#define SPATIALFILTER_H
#include <vector>
#include <numeric>
#include <algorithm>
#include <gdiplus.h>
#include <fstream>
#include <cmath>
namespace nsimgtk
{
template <typename pixelType, Gdiplus::PixelFormat pixelFormat, class FilterMask>
bool SpatialFilterAlgo(Gdiplus::Bitmap* const p_bitmap, FilterMask filterMask, unsigned int x, unsigned int y,
unsigned int width, unsigned int height)
{
if (p_bitmap == NULL)
{
return false;
}
if ((width + x > p_bitmap->GetWidth()) || (height + y >p_bitmap->GetHeight()))
{
return false;
}
Gdiplus::BitmapData bitmapData;
Gdiplus::Rect rect(x, y, width,height);
if (p_bitmap->LockBits(&rect, Gdiplus::ImageLockModeWrite, pixelFormat, &bitmapData) != Gdiplus::Ok)
{
return false;
}
pixelType *pixels = (pixelType*)bitmapData.Scan0;
const unsigned int m = filterMask.d_m; // mask's width
const unsigned int n = filterMask.d_n; // mask's height
std::vector<pixelType> tmpImage((m-1+width)*(n-1+height)); // extend image to use zero-padding
// copy original bitmap to extended image with zero-padding method
for (unsigned int row=0; row<height; ++row)
{
for (unsigned int col=0; col<width; ++col)
{
tmpImage[(col+m/2)+(row+n/2)*(bitmapData.Stride/sizeof(pixelType)+m-1)] =
pixels[col+row*bitmapData.Stride/sizeof(pixelType)];
}
}
// process every pixel with filterMask
for (unsigned int row=0; row<height; ++row)
{
for (unsigned int col=0; col<width; ++col)
{
// fill the "m*n" mask with the current pixel's neighborhood
for (unsigned int i=0; i<n; ++i)
{
for (unsigned int j=0; j<m; ++j)
{
filterMask.d_mask[i*m+j] = tmpImage[(col+j)+(row+i)*(bitmapData.Stride/sizeof(pixelType)+m-1)];
}
}
// replace the current pixel with filter mask's response
pixels[col+row*bitmapData.Stride/sizeof(pixelType)] = filterMask.response();
}
}
if (p_bitmap->UnlockBits(&bitmapData) != Gdiplus::Ok)
{
return false;
}
return true;
}
}
#endif
其中模板参数FilterMask即为滤波掩模算法。通常的滤波算法有均值滤波器,可以模糊化图像,去除图形中的细节部分,使得我们可以关注图像中较为明显的部分,均值滤波器用于周期性噪声。中值滤波器用于图像中存在椒盐噪声也即脉冲噪声的情况下。另外有基于一阶微分的Sobel梯度算子和基于两阶微分的拉普拉斯算子,它们往往被用于边缘检测中。
下面是一些滤波器算法的具体实现,所以滤波器算法都应该实现pixelType response()函数以及有d_mask,d_m,d_n成员,这可以通过继承__filteMask类获得(不需要付出虚函数代价)。
#ifndef SPATIALFILTER_H
#define SPATIALFILTER_H
#include <vector>
#include <numeric>
#include <algorithm>
#include <gdiplus.h>
#include <fstream>
#include <cmath>
namespace nsimgtk
{
// 滤波器掩模的基类,提供掩模大小d_m,d_n,掩模覆盖下的m*n个像素值d_mask
// others filterMask should inherit it
template <typename pixelType>
struct __filterMask
{
const unsigned int d_m;
const unsigned int d_n;
// image's pixels under the m*n filter mask
std::vector<pixelType> d_mask;
// filter mask's width and heigh must be a odd, if not, it will plus one for the width or the height
__filterMask(unsigned int m, unsigned int n)
: d_m(m%2 ? m:m+1), d_n(n%2 ? n:n+1), d_mask(d_m*d_n)
{
}
};
// 掩模权值为全1的均值滤波器
template <typename pixelType>
class averagingFilterMaskSp
: public __filterMask<pixelType>
{
public:
averagingFilterMaskSp(unsigned int m, unsigned int n)
: __filterMask<pixelType>(m, n)
{ }
pixelType response()
{
return std::accumulate(d_mask.begin(), d_mask.end(), 0) / (d_m * d_n);
}
};
// 可自定义掩模权值的均值滤波器
template <typename pixelType>
class averagingFilterMask
: public __filterMask<pixelType>
{
private:
std::vector<pixelType> d_weight; // weights' vector(m*n)
int d_weight_sum; // all weights' sum
public:
averagingFilterMask(unsigned int m, unsigned int n, const std::vector<pixelType>& weightVec)
: __filterMask<pixelType>(m, n), d_weight(weightVec)
{
if (weightVec.size() != d_mask.size())
{
// if weight's size isn't equal to mask's size, it will change filter mask as a special filter mask
d_weight.resize(d_mask.size(), 1);
}
d_weight_sum = std::accumulate(d_weight.begin(), d_weight.end(), 0);
}
pixelType response()
{
return std::inner_product(d_mask.begin(), d_mask.end(), d_weight.begin(), 0) / d_weight_sum;
}
};
// 中值滤波器
template <typename pixelType>
class medianFilterMask
: public __filterMask<pixelType>
{
public:
medianFilterMask(unsigned int m, unsigned int n)
: __filterMask<pixelType>(m, n)
{ }
pixelType response()
{
std::sort(d_mask.begin(), d_mask.end());
return d_mask[d_mask.size()/2];
}
};
// 3*3拉普拉斯滤波器
// the mask is: [0 1 0 [0 -1 0
// 1 -5 1 or -1 5 -1
// 0 1 0] 0 -1 0]
// if pixel's brightness is less than min, set it to min
// if pixel's brightness is larger than max, set it to max
template <typename pixelType, pixelType min, pixelType max>
class laplacianFilter
: public __filterMask<pixelType>
{
public:
laplacianFilter()
: __filterMask<pixelType>(3, 3)
{ }
pixelType response()
{
int ret = (int)(5*(int)d_mask[4]) - ((int)d_mask[5]+d_mask[3]+d_mask[1]+d_mask[7]);
if (ret < min)
ret = min;
if (ret > max)
ret = max;
return ret;
}
};
// 3*3Sobel滤波器
// the mask is: [-1 -2 -1 [-1 0 1
// 0 0 0 and -2 0 2
// 1 2 1] -1 0 1]
// if pixel's brightness is larger than max, set it to max
template <typename pixelType, pixelType max>
class sobelFilter
: public __filterMask<pixelType>
{
public:
sobelFilter()
: __filterMask<pixelType>(3, 3)
{ }
pixelType response()
{
int ret = ::abs(d_mask[6]+2*d_mask[7]+d_mask[8]-d_mask[0]-2*d_mask[1]-d_mask[2])
+ ::abs(d_mask[2]+2*d_mask[5]+d_mask[8]-d_mask[0]-2*d_mask[3]-d_mask[6]);
if (ret > max)
ret = max;
return ret;
}
};
}
#endif
转载于:https://www.cnblogs.com/songtzu/archive/2013/01/05/2845072.html
GDI+有Bitmap类。相关推荐
- c#图像处理入门(-bitmap类和图像像素值获取方法)
c#图像处理入门 -bitmap类和图像像素值获取方法 一.Bitmap类 Bitmap对象封装了GDI+中的一个位图,此位图由图形图像及其属性的像素数据组成.因此Bitmap是用于处理由像素数据定义 ...
- C#中的bitmap类和图像像素值获取方法
private void button1_Click_1(object sender, EventArgs e) { // 画直线 //Bitmap bitM ...
- 【转】c#数字图像处理(一)Bitmap类、 Bitmapdata类和 Graphics类
转自:https://www.cnblogs.com/dearzhoubi/p/8553763.html Bitmap类. Bitmapdata类和 Graphics类是C#图像处理中最重要的3个类, ...
- 在GDI+中用Mattix类对2D矢量图形进行平移、缩放操作
在GDI+中用Mattix类对2D矢量图形进行平移.缩放操作 1. GDI+中点坐标的格式及矩阵乘法的定义 GDI+中的的点按照1行3列的格式,即(x坐标,y坐标,1),其中1为哑元坐标.变换矩 ...
- C# 使用Bitmap类进行图片裁剪
各位朋友,大家好!前些天在博客园申请开通了自己的博客,主要是为了能和各位大侠和菜鸟分享软件开发这点乐趣,当然,也是自己经验积累和奋斗历程的简要版写照吧.好了,跑题了,下面就让我们一起去研究下C#中怎样 ...
- 用java写一个金字塔,一个Java写的用以构建影像金字塔的Bitmap类
一个Java写的用来构建影像金字塔的Bitmap类 一个Java写的用来构建金字塔影像的Bitmap类 cheungmine 2012 下面每个图像都是256x256像素.目的就是把这4幅影像合成一个 ...
- Bitmap类用法 详细说明
1. BitMap类 public void recycle()--回收位图占用的内存空间,把位图标记为Dead public final boolean isRecycled() --判断位图内 ...
- Bitmap类01_浅探
Bitmap类解析 Bitmap类是对图像进行处理的类,可以获取图像信息,进行图像颜色变换等操作 在安卓中,Bitmap指的是一张图片,可以是 .png 或者 .jpg 等其他常见的图像格式 参考AP ...
- GDI+——使用Graphics类绘制基本图形
目录 绘制基本图形 绘制直线 绘制矩形 绘制椭圆 绘制圆弧 绘制扇形 绘制多边形 绘制文本 综合案例 绘制波形曲线 使用双缓冲技术绘图 绘制基本图形 绘制直线 绘制直线用到了Graphics类的Dra ...
最新文章
- 2021年春季学期-信号与系统-第九次作业参考答案
- 单片机数码管从00到99C语言_MSP430单片机轻松入门与实践 — 畅学单片机
- oracle中if/else功能的实现的3种写法
- PHP魔法函数性能分析
- c语言vc数据类型长度,vc和gcc对C语言数据类型长度的定义
- 《ASP.NET Core 微服务实战》-- 读书笔记(第3章)
- C#学习之按钮点击事件
- STL容器之deque
- Vue商品添加到购物车
- 2019年最佳Python学习路线
- vue 文件转换二进制_vue项目将file转换成二进制流
- 【论文写作】招聘系统总体流程图如何画
- 【HDOJ6992】Lawn of the Dead(线段树×, 模拟大法好√)
- POSIX 与 CMSIS 标准
- 【Windows网络编程】完成端口IOCP原理及案例
- Flutter 与 Compose 应该怎么选择?它们冲突吗?
- java读取目录下所有txt_java读取同目录下的txt文件
- opencv人脸识别(一)调用笔记本摄像头
- React中setState() 函数的三种用法
- 【线性代数笔记】关于两个矩阵相乘等于零矩阵(AB=O)
热门文章
- python爬虫cookie池 与ip绑定_Python爬虫防封ip的一些技巧
- python在金融工程领域包括_金融工程专业对计算机的能力要求到底是什么?
- 计算机为什么启动二次才能打开,为什么电脑要2次重启才能启动显示器那
- python给出数据点进行插值_Python对数据进行插值和下采样的方法
- android多屏应用程序,微软也尝试“多屏协同” Windows系统可以运行安卓程序
- python函数调用的一般形式_Python的函数的定义与使用示例
- python实现弹幕_python实现b站直播自动发送弹幕功能
- MATLAB 长度和像素_MATLAB——单车道NaSch模型
- 砰的一声,实验室又炸鸡了
- 还没搞懂串口通信?一文带你读懂