ISP图像信号处理——平场校正介绍以及C++实现

news2024/11/26 7:25:38

参考文章1:http://t.csdn.cn/h8TBy

参考文章2:http://t.csdn.cn/6nmsT

参考网址3:opencv平场定标 - CSDN文库

平场校正一般先用FPN(Fixed Pattern Noise)固定图像噪声校正,即暗场校正;再用PRNU(Photo Response Non Uniformity)图像非均匀性响应校正,即明场校正;进行图像或者相机的校正。

平场校正就是以整帧图像的均值或者中值为目标图像,根据每个像素的特性不同,以响应增益gain为系数K和偏置offset为系数B,对每个像素进行校正,从而使得整幅图像看起来很均匀平滑,当然,为了简化计算,节省资源,也可以使用多个相邻的像素使用同一组校准参数。

一.按照算法的理论步骤

校正时,可以采用多段校正法将响应曲线分段校正,而更多的是简单暴力的两点校正法,也就是看成线性响应。首先相机对暗场进行一次曝光,得到每个像元的偏移(Offset);接下来对均匀光照条件下的灰度均匀物体进行一次成像,得到均匀场图像,最好能够使图像中所有的点都接近最大的灰度值;最后用均匀光场图像减去暗场图像,用相对标定的方法对图像增益(Gain)进行校正。

在暗场校正中,可以得到均值Vavgb,像素值Vinb;

在明场校正中,可以得到均值Vavgr,像素值Vinr;

可以得到增益响应系数K=(Vavgr-Vavgb)/(Vinr-Vinb);

偏置B=Vavgb-Vinb*K;

将n个像素的每一个像素点的K[n]和B[n]写入一个ram表中,当读取一帧图像时,对每一个像素进行校正Vout=Vin*K+B。必须首先做暗场再做明场,计算PRNU因子时要用到FPN因子的值。FPN校正时,像素值必须在1DN与127DN(可理解为灰度值)之间;PRNU校正时,像素值必须在128DN与254DN之间。

以上是针对相机或者图像偏数学算法原理的介绍。

上面这种方法感觉步骤不够具体,个人能力对于很多参数计算及过程难以理解,无法实现完整代码的撰写,有能力者希望可以指点一二。

二.利用平场定标技术校正(一组图像)

下面是我在另外一个博主这里看到的另一种方法平场定标,跟平场校正的功能比较像,这是针对一组拍摄好的图像做处理,我觉得很有启发:

平场定标(Flat-field calibration))是指通过对图像进行校正,消除由于光照不均匀或传感器噪声等因素引起的图像亮度和颜色的不均匀性。在OpenCV中,可以使用平场定标技术来校正图像。具体步骤如下:

1.收集平场图像:在相机的工作环境中,拍摄一系列完全均匀且没有目标物体的图像,这些图像被称为平场图像。

⒉.计算平场图像的平均值:将所有平场图像叠加在一起,并计算每个像素位置的平均值。这将得到一个平均平场图像。(三中用到)

3.计算平场图像的标准差:将所有平场图像与平均平场图像进行逐像素的差异计算,并计算每个像素位置的标准差。这将得到一个标准差平场图像。

4.校正图像:对于要校正的图像,将其与平均平场图像相减,并除以标准差平场图像。这将消除图像中的亮度和颜色不均匀性。

这种方法写了代码,但是最后得到的显示结果存疑,这里就不贴出来了。

三.常见的平场校正流程(两个输入图,一个输出图)

这个方法步骤是通过AI问答得到的,也具备一定的可借鉴性。

  1. 读取平场图像和待校正图像。(个人认为,平场图像可以用上面二中的平均平场图像。)
//个人根据现有资料理解所写,不确保正确性!
int main()
{
	//读入图片路径下的所有图像文件
	string imgdirpath="D:\\flatimages\\";
	string imgpath = imgdirpath + "*.bmp";  // test_images
	vector<String> filesVec;
	//遍历文件夹,两个参数分别表示文件的路径和保存图像具体路径的Vector容器。
	glob(imgpath, filesVec);
	if (filesVec.size() == 0)
	{
		cout << "File is empty!" << endl;
		return -1;
	}
	
	//初始化平场图像,利用一组亮度均匀符合要求的图像的平均图像做平场图像
	Mat flatsum = Mat::zeros(height, width, CV_8UC3);
	int num = 0;
	for (int i = 0; i < filesVec.size(); i++)
	{
		//叠加所有平场图像,并求平均值图像
		add(flatsum, ReadRaw10(filesVec[i].c_str()), flatsum);
		num++;
	}
	Mat flatmean3 = flatsum / num;//得到平均平场图像
    Mat inimg=imread("D:\\imgtest\\01.bmp");//读取待处理的输入图像
    Mat img_flat = flatfield2(inimg, flatmean3);//平场校正,
	imshow("img_flat", img_flat);
	waitKey(0);
  1. 如果没有一张符合要求的平场图像,需要利用一组图像计算平均图像用作平场图像。
  2. 对平场图像和待校正图像进行预处理,如图像增强、降噪等。
  3. 计算平场图像和待校正图像的平均值。
  4. 计算非均匀性模式:
  5. 将平场图像减去其平均值,得到平场图像的差异图像。
  6. 将待校正图像减去其平均值,得到待校正图像的差异图像。
  7. 将平场图像的差异图像除以待校正图像的差异图像,得到非均匀性模式。
  8. 校正待校正图像:
  9. 将待校正图像减去非均匀性模式,得到校正后的图像。
  10. 可选步骤:根据需要,进行其他的图像处理,如对比度调整、增强等。
  11. 输出校正后的图像。

以下是本人根据理解写的,代码下方的图像监视(VS插件,自行安装)得到的放大像素值截图(不是一一对应):

//(3)平场校正2。两个输入,一个输出,需要提前知道平场图像,用平均值
Mat flatfield2(Mat inImg,Mat flatImg)
{
	Mat grayflatImg;
	Mat grayinImg;
	cvtColor(flatImg, grayflatImg, COLOR_BGR2GRAY);
	cvtColor(inImg, grayinImg, COLOR_BGR2GRAY);
	//计算非均匀性模式
	Mat nonUniformityPattern = (grayflatImg - mean(grayflatImg))/ (grayinImg-mean(grayinImg));
	//待校正图像减去非均匀性模式,得到校正后的图像
	Mat flatcorrected1 = grayinImg - nonUniformityPattern;
	//转回彩色图像
	Mat flatcorrected3 = Mat::zeros(height, width, CV_8UC3);
	//cvtColor(flatcorrected1, flatcorrected3, COLOR_GRAY2BGR);
	cvtColor(flatcorrected1, flatcorrected3, COLOR_BayerBG2RGB);
	return flatcorrected3;
}

四.较为简单的步骤(一个输入图像,一个输出图像)

还有另一种方法,使用了OpenCV库来实现平场校正:

首先,将彩色图像转换为灰度图像,然后计算图像的平均值来获取亮度分布图。

接下来,将原始图像除以亮度分布图,得到校正后的图像。

最后,将校正后的图像转换回彩色图像并显示。

这种方法也写了代码,但是最后得到的显示结果也存疑,这里就不贴出来了。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1048285.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

c++颜色空间转换

c颜色空间转换 cvtColor函数标识符颜色空间转换的类型 cvtColor函数使用不同的标识符来表示颜色空间转换的类型。这些标识符通常位于OpenCV的cv命名空间中&#xff0c;并以cv::COLOR_作为前缀。以下是一些常用的cvtColor函数标识符&#xff1a; BGR到灰度&#xff1a; cv::CO…

JEP290深入学习

关键类的学习 ObjectInputFilter serialFilter字段的类型就是ObjectInputFilter&#xff0c;这个接口是一个函数接口&#xff08;可以对方法进行赋值&#xff09;&#xff0c;其中定义的抽象方法是&#xff1a; Status checkInput(FilterInfo filterInfo); 接受的参数类型F…

JavaSE 栈和队列

目录 1 栈 Stack1.1 概念1.2 实现1.3 Stack常见方法的使用1.4 前、中、后缀表达式 2 队列 Queue2.1 概念2.2 实现2.3 循环队列2.4 双端队列 3 Java中的栈和队列4 练习题4.1 括号匹配问题4.2 用队列实现栈4.3 用栈实现队列4.4 实现一个最小栈4.5 设计循环队列 1 栈 Stack 1.1 概…

LIMS实验室信息管理系统源码 基于计算机的数据处理技术、数据存储技术、网络传输技术、自动化仪器分析技术于一体

LIMS 是一个集现代化管理思想与基于计算机的数据处理技术、数据存储技术、网络传输技术、自动化仪器分析技术于一体&#xff0c;以实验室业务和管理工作为核心&#xff0c;遵循实验室管理国际规范&#xff0c;实现对实验室全方位管理的信息管理系统。 LIMS将样品管理、数据管理…

前端利用策略模式优化if else

策略模式 策略模式定义了一系列算法,并将每种算法封装起来,使它们可以相互替换,且算法的变化不会影响使用算法的客户。 策略模式优化if else 策略模式可以用来代替多重if else语句。 Vue中使用策略模式 我们以一个表单验证的例子来看看如何使用策略模式优化if else。 methods:…

一文让你搞懂async/await存在时的执行顺序

最近在写项目&#xff0c;发现代码中总出现这样的内容&#xff0c;一直百思不得其解&#xff0c;后来同事说这段代码的意义是把后边的任务变成宏任务执行&#xff0c;这我更迷惑了&#xff0c;我发现自己的Promise基础还是很差&#xff0c;因此在摸索了几天后&#xff0c;感觉好…

SAP ABAP内部单位外部单位转换问题

之前没有留意&#xff0c;用户用了那么久也没有发现&#xff0c;今天用户找过来说他在前台维护的是PC&#xff0c;但是报表里面显示是ST。 看了一下原来这里的“内部度量单位”就是内部的、与语言无关的计量单位格式&#xff08;ST&#xff09;。该内部格式在输出到屏幕之前&a…

完整指南:如何使用 Node.js 复制文件

文件拷贝指的是将一个文件的数据复制到另一个文件中&#xff0c;使目标文件与源文件内容一致。Node.js 提供了文件系统模块 fs&#xff0c;通过该模块可以访问文件系统&#xff0c;实现文件操作&#xff0c;包括拷贝文件。 Node.js 中文件拷贝方法 在 Node.js 中&#xff0c;有…

竞赛选题 大数据商城人流数据分析与可视化 - python 大数据分析

0 前言 &#x1f525; 优质竞赛项目系列&#xff0c;今天要分享的是 &#x1f6a9; 基于大数据的基站数据分析与可视化 该项目较为新颖&#xff0c;适合作为竞赛课题方向&#xff0c;学长非常推荐&#xff01; &#x1f947;学长这里给一个题目综合评分(每项满分5分) 难度…

Postgresql源码(114)视图权限授予逻辑

0 速查 被授权的对象在系统表中记录授权信息&#xff0c;例如pg_namespace中的nspacl列&#xff1a; {mingjieUC/mingjie,UC/mingjie,pusr1UC/mingjie}pusr1UC/mingjie的含义&#xff1a; mingjie是赋予者pusr1是被赋予者UC是权限&#xff0c;表示USAGE和CREATE 1 视图权限…

@SpringBootApplication注解的理解——如何排除自动装配 分布式情况下如何自动加载 nacos是怎么被发现的

前言 spring作为主流的 Java Web 开发的开源框架&#xff0c;是Java 世界最为成功的框架&#xff0c;持续不断深入认识spring框架是Java程序员不变的追求。 本篇博客介绍SpringBootApplicant注解的自动加载相关内容 其他相关的Spring博客文章列表如下&#xff1a; Spring基…

ubuntu20.04 jammy 安装ros2

ubunut22.04 jammy&#xff08;5.15&#xff09; ros2版本: humble 安装参考&#xff1a; Ubuntu (Debian packages) — ROS 2 Documentation: Humble documentationl 按照官方给的操作指南进行操作即可&#xff0c;到安装软件包的时候&#xff0c;若只为开发&#xff0…

C++中实现雪花算法来在秒级以及毫秒及时间内生成唯一id

1、雪花算法原理 雪花算法&#xff08;Snowflake Algorithm&#xff09;是一种用于生成唯一ID的算法&#xff0c;通常用于分布式系统中&#xff0c;以确保生成的ID在整个分布式系统中具有唯一性。它的名称来源于雪花的形状&#xff0c;因为生成的ID通常是64位的整数&#xff0…

【小白专属02】SpringBoot集成MybatisPlus

目录 前言 一、添加依赖项 二、配置数据库连接信息 三、数据库创建表 四、创建实体表 五、编写Mapper接口 六、功能测试 前言 上节回顾 上一节我们对SpringBoot框架进行了搭建&#xff0c;其实按照流程走问题不大。最主要的部分我们把下载SpringBoot demo的网址换成了…

编程每日一练(多语言实现)基础篇:百元买百鸡

文章目录 一、实例描述二、技术要点三、代码实现3.1 C 语言实现3.2 Python 语言实现3.3 Java 语言实现3.4 JavaScript 语言实现 一、实例描述 中国古代数学家张丘建在他的 《算经》 中提出了一个著名的 "百钱买百鸡问题"&#xff0c; 鸡翁一&#xff0c;值钱五&…

重新定义旅行,AnyGo for Mac让你畅享虚拟GPS位置的无限可能!

旅行是一种令人兴奋和激动人心的体验&#xff0c;而现在有了AnyGo for Mac这款强大的虚拟GPS位置工具&#xff0c;你将能够重新定义旅行&#xff0c;尽情探索全新的地理环境和体验&#xff01; AnyGo for Mac为你带来了全新的旅行方式。无需离开家门&#xff0c;你就可以通过模…

基于php的物流信息公共平台设计与实现

项目描述 临近学期结束&#xff0c;还是毕业设计&#xff0c;你还在做java程序网络编程&#xff0c;期末作业&#xff0c;老师的作业要求觉得大了吗?不知道毕业设计该怎么办?网页功能的数量是否太多?没有合适的类型或系统?等等。今天给大家介绍一篇基于php的物流信息公共平…

面对海量数据挑战,企业怎样进行数据处理?

海量数据给企业带来机遇的同时&#xff0c;也给企业带来了一定的挑战。怎样高效处理海量数据也是成为了企业面临的重要问题。 一&#xff1a;数据采集 数据采集是数据处理的第一步。企业需要进行数据采集&#xff0c;但是要保证数据的准确性以及安全性。同时需要注意&#xf…

PgSQL-向量数据库插件-lantern

PgSQL-向量数据库插件-lantern 即pgvector、pg_embedding 后又一向量数据库扩展Lantern问世了。当然也为向量列提供了hnsw索引以加速ORDER BY... LIMIT查询。Lantern使用usearch实现hnsw。 使用方法 保留了标准PgSQL接口&#xff0c;兼容其生态工具。首先需要安装该插件&#x…

服务器基础知识:raid卡有什么优势?

当今直播电商、短视频呈现热火朝天的趋势&#xff0c;是近两年的商业热门议题&#xff0c;以淘宝天猫、拼多多、抖音、小红书、视频号为代表的平台数据量呈现爆炸性的增长。据IDC预测&#xff0c;2025年全球数据总量将增至175ZB&#xff08;1ZB1万亿GB&#xff09;&#xff0c;…