202104-2邻域均值---前缀和

试题链接：首页 - 计算机软件能力认证考试系统

【试题背景】

顿顿在学习了数字图像处理后，想要对手上的一副灰度图像进行降噪处理。不过该图像仅在较暗区域有很多噪点，如果贸然对全图进行降噪，会在抹去噪点的同时也模糊了原有图像。因此顿顿打算先使用邻域均值来判断一个像素是否处于较暗区域，然后仅对处于较暗区域的像素进行降噪处理。

【问题描述】

待处理的灰度图像长宽皆为 n 个像素，可以表示为一个 n×n 大小的矩阵 A，其中每个元素是一个 [0,L) 范围内的整数，表示对应位置像素的灰度值。
对于矩阵中任意一个元素 Aij（0≤i,j<n），其邻域定义为附近若干元素的集和：

Neighbor(i,j,r)={Axy|0≤x,y<n and |x−i|≤r and |y−j|≤r}

这里使用了一个额外的参数 r 来指明 Aij 附近元素的具体范围。根据定义，易知 Neighbor(i,j,r) 最多有 (2r+1)2 个元素。

如果元素 Aij 邻域中所有元素的平均值小于或等于一个给定的阈值 t，我们就认为该元素对应位置的像素处于较暗区域。
下图给出了两个例子，左侧图像的较暗区域在右侧图像中展示为黑色，其余区域展示为白色。

现给定邻域参数 r 和阈值 t，试统计输入灰度图像中有多少像素处于较暗区域。

【输入形式】

输入共 n+1 行。

输入的第一行包含四个用空格分隔的正整数 n、L、r 和 t，含义如前文所述。

第二到第 n+1 行输入矩阵 A。
第 i+2（0≤i<n）行包含用空格分隔的 n 个整数，依次为 Ai0,Ai1,⋯,Ai(n−1)。

【输出形式】

输出一个整数，表示输入灰度图像中处于较暗区域的像素总数。

【样例1输入】

4 16 1 6

0 1 2 3

4 5 6 7

8 9 10 11

12 13 14 15

【样例1输出】

7

【样例2输入】

11 8 2 2

0 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 0 0 0 0

0 7 0 0 0 7 0 0 7 7 0

7 0 7 0 7 0 7 0 7 0 7

7 0 0 0 7 0 0 0 7 0 7

7 0 0 0 0 7 0 0 7 7 0

7 0 0 0 0 0 7 0 7 0 0

7 0 7 0 7 0 7 0 7 0 0

0 7 0 0 0 7 0 0 7 0 0

0 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 0 0 0 0

【样例2输出】

83

这道题就是用二维前缀和

        首先初始化二维前缀和数组

int b[n+1][n+1];
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            int cur;
            cin>>cur;
            b[i][j]=b[i-1][j]+b[i][j-1]-b[i-1][j-1]+cur;
        }
    }

之后就是遍历每一个像素点，计算该像素点所在邻域内，所有点的像素之和，并计算平均值；

然后将平均值（aver）与t进行比对，如果aver小于等于t，则处于暗区域的像素点个数（num）++；

for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            //求邻域左右边界值
            int lefti=max(1,i-r),righti=min(n,i+r),leftj=max(1,j-r),rightj=min(n,j+r);
            //求邻域内值
            int sum=b[righti][rightj]-b[righti][leftj-1]-b[lefti-1][rightj]+b[lefti-1][leftj-1];
            //求平均值
            double aver=sum/((righti-lefti+1)*(rightj-leftj+1)*1.0);
            if(aver<=t) num++;
        }
    }

最后输出num即可。

总的代码

#include <iostream>

using namespace std;

int main()
{
    int n,L,r,t,num=0;
    cin>>n>>L>>r>>t;
    int b[n+1][n+1];
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            int cur;
            cin>>cur;
            b[i][j]=b[i-1][j]+b[i][j-1]-b[i-1][j-1]+cur;
        }
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            //求邻域左右边界值
            int lefti=max(1,i-r),righti=min(n,i+r),leftj=max(1,j-r),rightj=min(n,j+r);
            //求邻域内值
            int sum=b[righti][rightj]-b[righti][leftj-1]-b[lefti-1][rightj]+b[lefti-1][leftj-1];
            //求平均值
            double aver=sum/((righti-lefti+1)*(rightj-leftj+1)*1.0);
            if(aver<=t) num++;
        }
    }
    cout<<num;
    return 0;
}

参考大佬文章：

CSP202104-2邻域均值（二维前缀和）_疼疼疼疼疼的博客-CSDN博客

前缀和应用总结——前缀和及矩阵前缀和_CAccept的博客-CSDN博客_矩阵前缀和