ISP图像处理中大数乘法溢出处理(一)

[二、方案选型（按 ISP 嵌入式工程推荐优先级排序）](#二、方案选型（按 ISP 嵌入式工程推荐优先级排序）)

[方案 1：64 位中间变量中转](#方案 1：64 位中间变量中转)

[方案 2：无 64‑bit 硬件（极老旧 MCU，无uint64_t/long long）------ 带余数补偿拆分算法](#方案 2：无 64‑bit 硬件（极老旧 MCU，无uint64_t/long long）—— 带余数补偿拆分算法)

[方案 3：溢出预判 + 分支分流（不推荐主路径，适合异常保护）](#方案 3：溢出预判 + 分支分流（不推荐主路径，适合异常保护）)

[方案 4：定点 Q 格式缩放（适合无除法硬件、硬件加速 ISP 算子）](#方案 4：定点 Q 格式缩放（适合无除法硬件、硬件加速 ISP 算子）)

各方案对比表

[ISP 直方图场景额外工程建议](#ISP 直方图场景额外工程建议)

一、问题

isp 图像处理中，需要统计0-255所有bin d的直方图信息之和。一般ISP会设置不同区域的权重再统计每一个像素进行构建直方图。如果分辨率比较大的时候对其做加权处理后直方图所有bin之和会是一个比较大的数值all_total，现在我想对取所有像素的百分比的时候，比如 all_total*ratio/100, 设定其其像素进行统计百分比阈值，当大于某一阈值进行一些逻辑操作。在32bit系统中，当 all_total比较大的时候all_total*ratio对 unsigned long 范围可能溢出，需要做溢出处理。

复制代码

// 危险代码，极易溢出
unsigned long total;
uint16_t ratio;
unsigned long thresh = (total * ratio) / 100U;

total是直方图像素总和（32‑bit unsigned long，最大值 0xFFFFFFFF = 4294967295）；先执行乘法 total * ratio，乘法结果超出 32 位范围就发生无符号溢出回绕，数值错乱，之后再除以 100 得到完全错误阈值。

溢出临界点：

例：ratio=100，安全上限 = 42,949,672；只要直方图总像素＞该值，乘法直接溢出。

二、方案选型（按 ISP 嵌入式工程推荐优先级排序）

方案 1：64 位中间变量中转【首选！如果系统为64bit优选】

核心技巧：先把其中一个操作数强制提升到 64bit 再相乘，乘法全程在 64 位空间运算，永远不会溢出 32bit 像素总和场景。

⚠️致命坑：不能 (uint64_t)(total*ratio)，此时乘法已经 32 位溢出后再强转，为时已晚！

复制代码

#include <stdint.h>

uint32_t CalcPercentThresh(uint32_t total, uint16_t ratio)
{
    /* (uint64_t)total 先提升位宽；ratio自动提升为64bit相乘 */
    uint64_t temp = (uint64_t)total * ratio;
    temp = temp / 100U;
    /* 结果一定 <= total，安全转回32bit */
    return (uint32_t)temp;
}

适用平台：所有带 C99 编译器嵌入式64bit Linux；
优点：零精度损失，逻辑最简单，无误差，无分支判断；
上限能力：total ≤ 4G、ratio ≤ 100，乘积最大 400G，远小于 64 位上限 (1.8e19)，绝对不会溢出。

ISP 直方图百分比阈值判断示例（累计直方图，大于阈值触发逻辑）

复制代码

uint32_t total_sum = ...; //直方图所有bin总和
uint32_t percent_thresh;
percent_thresh = (uint32_t)((uint64_t)total_sum * 95U / 100U); 

//遍历累计直方图，找到第一个累计>percent_thresh的灰度位置做阈值

方案 2：无 64‑bit 硬件（极老旧 MCU，无`uint64_t/long long`）------ 带余数补偿拆分算法

直接 (total / 100)*ratio 会产生截断精度丢失；利用恒等式做精确拆分：

A=total，B=ratio，C=100。

复制代码

uint32_t MulDiv32_NoOverflow(uint32_t A, uint32_t B, uint32_t C)
{
    uint32_t div  = A / C;
    uint32_t rem  = A % C;
    uint32_t part1= div * B;
    uint32_t part2= (rem * B)/ C;
    return part1 + part2;
}
//调用：thresh = MulDiv32_NoOverflow(total, ratio, 100U);

约束：rem < C，本例 C=100，rem*ratio最大 99*100=9900，永远不会溢出 32 位，完美规避大数乘法溢出。

优点：纯 32 位运算，无 64 位依赖，数学结果完全等价 (A*B)/C，精度无损；
缺点：多一次除法、取余运算，少量耗时；适合无 64‑bit 的极简内核。

方案 3：溢出预判 + 分支分流（不推荐主路径，适合异常保护）

先预判乘法会不会溢出，溢出时走备选计算路径。溢出判定条件：

复制代码

/* 判断 total * ratio 是否溢出32bit无符号 */
if (total > UINT32_MAX / ratio)
{
    //溢出分支：使用拆分除法方案
}
else
{
    //安全分支：直接运算
    thresh = (total*ratio)/100U;
}

缺点：代码双分支，逻辑复杂，ISP 流水线不喜欢条件分支，容易引入跳转延迟。

方案 4：定点 Q 格式缩放（适合无除法硬件、硬件加速 ISP 算子）

如果你不想运行时除法，可提前转为定点运算，例如 Q8（放大 256 倍）：

同样乘法必须 64 位中转防止溢出：

复制代码

uint32_t q_factor = (ratio * 256U)/100U;
uint32_t thresh = (uint32_t)(((uint64_t)total * q_factor) >> 8U);

适合硬件 ISP 内核，移位代替除法，适合做 LCE/AE 实时直方图统计。

各方案对比表

方案	是否需要 64 位	精度损失	ISP 工程推荐度
64bit 中转 (uint64_t)	是	无损失	⭐⭐⭐⭐⭐首选
带余数拆分 MulDiv32	否	无损失	⭐⭐⭐⭐ 无 64bit 时使用
先除后乘无余数补偿	否	有截断误差	❌禁止使用
溢出预判分支	是 / 否	无损失	⭐⭐ 仅做防护
Q 定点缩放	是	微小舍入误差	⭐⭐⭐适合硬件算子

ISP 直方图场景额外工程建议

永远优先用 (uint64_t) 强转方案，如果图像信号处理器 ARM 内核都支持 64 位，优先使用；
不要为了省一点点周期牺牲溢出安全性，直方图总像素经常可达 200 万、800 万、1200 万像素，极易触发 32‑bit 乘法溢出；
百分比阈值逻辑（90%、95% 累计直方图截断，DRC / 阴影抑制）的工业标准代码范式：

//累计直方图百分比阈值
uint32_t hist_sum = GetHistTotalSum();
uint32_t cumu_th = (uint32_t)((uint64_t)hist_sum * 95U / 100U);

uint32_t cum_sum = 0;
uint32_t gray_th = 255;
for(int bin=0;bin<=255;bin++)
{
cum_sum += hist[bin];
if(cum_sum > cumu_th)
{
gray_th = bin;
break;
}
}
//使用gray_th做后续图像亮度分割逻辑