图像处理：马赛克与毛玻璃算法

马赛克：区域内随机取样覆盖整个区域
毛玻璃：区域内随机取样代替单个像素点

太久没碰c，忘得差不多了orz~

#define min(a,b) a<b?a:b
void MorA(int w, int h, int k1, int k2)
{
    //w:图片宽,h:图片高,k1：随机取样半径，k2像素覆盖半径
    for (int i = 0; i < w; i += k2)
    {
        for (int j = 0; j < h; j += k2)
        {
            int rx = i + rand() % k1;
            rx = min(rx, w - 1);
            int ry = j + rand() % k1;
            ry = min(ry, h - 1);

            int m = a[rx][ry];

            int lx = min(i + k2, w - 1);
            int ly = min(j + k2, h - 1);
            for (int x = i; x < lx; x++)
                for (int y = j; y < ly; y++)
                    b[x][y] = m;
        }
    }
}

为了不让指针溢出，写了四个三目运算符做安全检查，所以说溢出真是个很麻烦的东西。

在main函数内测试：

#define f(n,m) for(int i=0;i<n;i++)for(int j=0;j<m;j++)
#define pfln printf("\n")
#define pfA(n,m,a) for(int i=0;i<n;i++,pfln)for(int j=0;j<m;j++)printf("%d",a[i][j]);
int main()
{
    freopen("in.txt", "r", stdin);
    int n, m; scanf("%d%d", &n, &m);
    f(n, m)scanf("%d", &a[i][j]);
    pfA(n, m, a); pfln;
    MorA(n, m, 2, 4);
    pfA(n, m, b); pfln;
    MorA(n, m, 2, 1);
    pfA(n, m, b); pfln;
}

输出如下

可以发现雾化因为选择像素范围的原因会往左上角飘，还有待优化。
接下来开始处理真实图片
原图

k1=3,k2=3

k1=5,k2=5

k1=10,k2=10

k1=5,K2=1

k1=10,k2=1

后续更改：
关于雾化会往左上角偏移的优化：

#define max(a,b) a>b?a:b
#define min(a,b) a<b?a:b
void MorA(int w, int h, int k1, int k2)
{
    //w:图片宽,h:图片高,k1：随机取样半径，k2像素覆盖半径
    for (int i = 0; i < w; i += k2)
    {
        for (int j = 0; j < h; j += k2)
        {
            int rx = i - k1 / 2 + rand() % k1;
            rx = min(max(rx, 0), w - 1);
            int ry = j - k1 / 2 + rand() % k1;
            ry = min(max(ry, 0), h - 1);

            int m = a[rx][ry];

            int localx = min(i + k2, w - 1);
            int localy = min(j + k2, h - 1);
            for (int x = i; x < localx; x++)
                for (int y = j; y < localy; y++)
                    b[x][y] = m;
        }
    }
}

提升空间：
提高处理性能：
1.使用位运算比较坐标：

int a,b,p;
p=a>>31;
a=(a^p)-p;//取绝对值
p=-(a>b)?
a=(a&p)|(b&~p);//取大值 

2.图像长宽各扩充k1像素,一劳永逸取消越界判定
1适合小图像，2适合大图像

减少内存开销：
处理马赛克不需要b数组，减少内存开销
处理雾化时b数组为b[k1][w]，减少内存开销

提高处理效果：
使取样在以k1为半径的圆范围内
以后再实现吧。