非常抱歉，我的回答中存在错误。
这里是一个更正的代码：
- % 读取待校准的图像
- img = imread('test.jpg');
- % 转换为灰度图像
- if size(img, 3) == 3
- img = rgb2gray(img);
- end
- % 读取灰度参考目标
- ref = imread('reference.jpg');
- % 转换为灰度图像
- if size(ref, 3) == 3
- ref = rgb2gray(ref);
- end
- % 最大值归一化
- img = double(img) / double(max(img(:)));
- ref = double(ref) / double(max(ref(:)));
- % 计算平均灰度值
- mean_ref = mean(ref(:));
- mean_img = mean(img(:));
- % 校准图像
- max_val = double(intmax(class(img)));
- img_calibrated = img .* (max_val./double(img));
- % 显示结果
- figure;
- subplot(1, 2, 1); imshow(img); title('Original');
- subplot(1, 2, 2); imshow(img_calibrated); title('Calibrated');
请注意，我在此处添加了 `max_val` 变量，并且使用了点除号来执行归一化。此外，由于 `intmax` 函数返回的值的类型与 `img` 的类型相同，因此可以使用 `class(img)` 来动态计算最大值。
--- 光学专家Gpt |
