(2)计算atmospheric light

之前做了一个铺垫,拿到了一张图片的rgb中最小的值,接下来还会涉及到几行代码,其中有些语法不太懂,但是发现去百度搜索的一些博客说的不太对,还是谨慎一点,不找野路子,再看看官方文档吧。

首先是:

 darkc = min(im_min,[],3);   %得到暗通道——im_min[][]r,g,b空间中的的最小值, darkc是一个270*480*1的矩阵,代表每个点的3维的最小值

得到rgb矩阵的最小值——三维缩成了一维。还是举个例子吧——所有的人都不喜欢看文字,看例子,一眼明了:

% 创建三维数组:
d1 = [1,2,3;4,5,6;7,8,9;];
d2 = [4,5,6;7,8,9;1,2,3;];
d3 = [7,8,9;1,2,3;4,5,6;];
im_min(:,:,1) = d1;
im_min(:,:,2) = d2;
im_min(:,:,3) = d3;
im_min
im_min = min(im_min,[],3)

运行结果如下:



接下来,第二行代码

[~, idx] = sort(darkc(:),1,'descend');

首先,darkc(:):   Use A(:), the column representation of

C++实现大气校正代码可以使用OpenCV库来进行图像处理和计算。下面是一个简单的示例代码,用于实现大气校正: ```cpp #include <iostream> #include <opencv2/opencv.hpp> using namespace std; using namespace cv; Mat atmosphericCorrection(Mat inputImage, float A, float t0, float p) { Mat outputImage = inputImage.clone(); // 计算大气光照 Scalar atmosphericLight(A, A, A); // 对每个像素进行校正 for (int i = 0; i < inputImage.rows; i++) { for (int j = 0; j < inputImage.cols; j++) { Vec3f pixel = inputImage.at<Vec3f>(i, j); // 计算透射率 float t = 1 - p * min(pixel[0] / atmosphericLight[0], min(pixel[1] / atmosphericLight[1], pixel[2] / atmosphericLight[2])); // 校正像素值 for (int k = 0; k < 3; k++) { pixel[k] = (pixel[k] - atmosphericLight[k]) / max(t, t0) + atmosphericLight[k]; } outputImage.at<Vec3f>(i, j) = pixel; } } return outputImage; } int main() { // 读取输入图像 Mat inputImage = imread("input.jpg", IMREAD_COLOR); // 设置参数 float A = 220; // 大气光照强度 float t0 = 0.1; // 最小透射率 float p = 0.1; // 大气光照比例 // 进行大气校正 Mat outputImage = atmosphericCorrection(inputImage, A, t0, p); // 显示结果 imshow("Input Image", inputImage); imshow("Output Image", outputImage); waitKey(0); return 0; } ``` 上述代码中,`atmosphericCorrection`函数实现了大气校正的算法。它接受输入图像、大气光照强度A、最小透射率t0和大气光照比例p作为参数,并返回校正后的图像。在`main`函数中,我们读取输入图像,设置参数,然后调用`atmosphericCorrection`函数进行大气校正,并显示结果图像。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值