Float扩散算法数值异常暴涨问题咨询(网格智能体位置值扩散场景)
为啥你的扩散函数会让数值疯涨?
嘿,我一眼就瞅出你代码里的几个致命问题,直接导致数值越跑越高,来给你挨个捋清楚:
1. 权重算错了,总质量直接爆炸
你写的权重是(1 / blurRad * blurRad),这完全是个逻辑错误:
- 要是
blurRad是整数,1/blurRad会触发整数除法,比如blurRad=2时,1/2结果是0,乘完还是0——这时候等于没扩散,但你说数值暴涨,那大概率你用的是浮点数?那1/blurRad*blurRad结果就是1,相当于每个像素把自己的完整值加到周围所有格子里!总质量直接变成原来的(2*blurRad+1)^2倍,这不暴涨才怪! - 正确的权重应该是1除以扩散区域的总像素数,也就是
1.0 / ((2*blurRad + 1)*(2*blurRad + 1)),这样每次扩散后,总数值和之前保持一致,不会凭空变多。
2. 循环逻辑乱了,重复叠加还乱覆盖
看你内层循环里的这两行:
temp[i + k][j + l] += (1 / blurRad * blurRad) * grid[i][j]; temp[i][j] = (1 / blurRad * blurRad) * grid[i][j];
这俩操作全错了:
- 你每次处理
grid[i][j]的邻居时,给temp[i+k][j+l]加值,但转头又直接给temp[i][j]赋值,这会把之前其他像素扩散过来的值直接覆盖掉,逻辑彻底混乱。 - 更糟的是,每个像素会被多个来源重复加值:比如处理
grid[x][y]时会给temp[i][j]加一次,处理grid[i][j]的邻居时又会再加一次,加上权重是1(错误计算),总数值直接指数级增长。
3. 没初始化temp,旧值残留越积越多
你扩散前没把temp数组清零,这意味着上一帧的数值还留在里面,新的扩散值又叠加上去,数值自然越来越大。
给你改好的扩散函数
我把代码调整了,解决这些问题:
void diffuse() { // 先把temp数组清0,避免残留上一帧的旧值 // 要是二维数组用memset有问题,就用双重循环逐个赋值0 memset(temp, 0, sizeof(temp)); // 计算扩散区域的总像素数,得到归一化权重 int totalPixels = (2 * blurRad + 1) * (2 * blurRad + 1); float weight = 1.0f / totalPixels; for (int i = 0; i < width; i++) { for (int j = 0; j < height; j++) { float currentVal = grid[i][j]; // 遍历当前像素的所有扩散邻居 for (int k = -blurRad; k <= blurRad; k++) { int ni = i + k; if (ni < 0 || ni >= width) continue; // 跳过边界外的格子 for (int l = -blurRad; l <= blurRad; l++) { int nj = j + l; if (nj < 0 || nj >= height) continue; // 把当前像素的一部分值分配到对应的temp位置 temp[ni][nj] += currentVal * weight; } } } } // 最后把扩散后的temp值复制回grid,完成这一帧的扩散 memcpy(grid, temp, sizeof(grid)); // 同样,二维数组如果memcpy有问题,就用双重循环逐个复制 }
关于Float扩散算法的暴涨问题
如果换了浮点数版本还是疯涨,核心问题还是没做好归一化——也就是每次扩散时,所有分配出去的权重加起来得等于1,不然总质量就会越变越大。另外,边界处理如果截断了扩散范围,你可能需要用镜像或者边界重复的方式补全,避免总质量变化,但暴涨的话肯定是权重没算对或者重复叠加了。浮点数精度误差只会导致细微漂移,不会急剧暴涨。
内容的提问来源于stack exchange,提问作者Tristen Gordon
相关产品推荐
相关产品推荐

