红外非均匀矫正matlab实现,红外图像非均匀矫正——两点矫正
红外芯片由于工艺问题存在严重的分均匀性,所以非均匀矫正一直是影响红外图像质量的第一因素。分均匀矫正的算法也是红外图像处理研究的重点区域,建立了一些矫正的方式方法。其中最常用最简单的就应该算是两点温度定标算法。
应用两点法校正有两个前提条件,第一,探测器的响应在所关注的温度范围内是线性变化的,第二,探测器的响应具有时间的稳定性,并且其受随机噪声的影响较小,则非均匀性引入固定模式的乘性和加性噪声。
具体过程基本分为三步:
一、两温度下采集均匀原始数据
二、根据原始数据计算校正参数
三、将矫正参数带入实时原始图进行校正
matlab实现的的.m文件内容如下:
clear;
clc;
%%%%%%%%%%%%%%低温下对9帧图像取平均值%%%%%%%%%%%%%
save1=zeros(3,2,9,'uint16');
for i = 1:9; %连续读M帧图像,存入三维数组
fid = textread(strcat('300',num2str(i),'b.txt'),'%s');
fid1 = hex2dec(fid);
save1(:,:,i) = reshape(fid1,3,2);
end
%求M帧的平均值,消除偶然误差
sample1 = zeros(3,2,'double');
sample1 = double(save1(:,:,1)+save1(:,:,2)+save1(:,:,3)+save1(:,:,4)+save1(:,:,5)+save1(:,:,6)+save1(:,:,7)+save1(:,:,8)+save1(:,:,9))/9;
subplot(2,2,1), imshow(uint8(sample1));
title('T1温度原始图');
%%%%%%%%%%%%%低温下对9帧图像取平均值%%%%%%%%%%%%%%%
save2=zeros(3,2,9,'uint16');
for i = 1:9;
fid = textread(strcat('700',num2str(i),'b.txt'),'%s');
fid1 = hex2dec(fid);
save2(:,:,i) = reshape(fid1,3,2);
end
sample2 = zeros(3,2,'double');
sample2= double(save2(:,:,1)+save2(:,:,2)+save2(:,:,3)+save2(:,:,4)+save2(:,:,5)+save2(:,:,6)+save2(:,:,7)+save2(:,:,8)+save2(:,:,9))/9;
subplot(2,2,2),imshow(uint8(sample2));
title('T2温度原始图');
%%%%%%%%%%%%%%%计算校正因子%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
avr1=sum(sample1(:))/(3*2);
avr2=sum(sample2(:))/(3*2);
G = zeros(3,2,'double');
O = zeros(3,2,'double');
M = zeros(3,2,'double');
for i=1:3;
for j=1:2;
G(i,j) = (avr2-avr1)/(sample2(i,j)-sample1(i,j)); %计算增益
O(i,j) = avr1-G(i,j)*sample1(i,j); %计算偏移量
%O(i,j) =(avr2*double(sample1(i,j))-avr1*double(sample2(i,j)))/double(sample1(i,j)-sample2(i,j));
end
end
file = fopen('2dian4p000s17.txt','wt');
count = fwrite(file,G,'uint16');
%%%%%%%%%%%%%%校正场景图片%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fid = textread('4p000s.txt','%s');
fid1 = hex2dec(fid);
test2=reshape(fid1,3,2);
%%%%%%%%%%%%%%%%%%%%%%%%%校正得到全部结果%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:3;
for j=1:2;
M(i,j) = G(i,j)*double(test2(i,j))+O(i,j);
end
end
for i=1:3;
for j=1:2;
final(i,j) = uint16(M(i,j));
end
end
subplot(2,2,3),imshow(uint8(final));title('校正后 ');
subplot(2,2,4),imshow(uint8(test2));title('非均匀原始图');
file = fopen('2dian4p000s16.txt','wt');
count = fwrite(file,final,'double');
运行结果如下所示:
:
可以自己造图才跑,数据文件上传太繁琐了就没上传,如果有兴趣可以找我来要。
QQ:356636122
红外非均匀矫正matlab实现,红外图像非均匀矫正——两点矫正相关推荐
- matlab 刻度非均匀控制,MATLAB 出一张好看的图
1.坐标轴的视点(viewpoint):从哪个方向看整个坐标系统,这决定了坐标轴的方向和位置,通过view函数实现视点的设置: view([z y ]):(将坐标系统想象为一座房子,而自己是个会飞的天 ...
- matlab中predictor怎么填,在MATLAB中求解非線性有限元
我嘗試在MATLAB中求解帶有節點熱源的四面體有限元的熱擴散問題,這個節點取決於解矢量.非線性方程系統如下:在MATLAB中求解非線性有限元 乙U」 + A U = Q(T) 與B是熱capactiy ...
- 2ASK非相干解调matlab
2ASK非相干解调matlab 2ASK信号产生 振幅键控是利用载波的幅度变化来传递数字信息,而其频率和初始相位保持不变.在2ASK中,载波的幅度只有两种变化状态,分别对应二进制信息"0&q ...
- MATLAB程序采用非支配排序遗传算法(NSGA2)求解分布式电源选址定容问题
MATLAB程序采用非支配排序遗传算法(NSGA2)求解分布式电源选址定容问题,可作为一个有用的参考,程序注释明确,算法原理可以自己搜. 现有:6920651507678049浪迹天涯
- MATLAB程序采用非支配排序遗传算法(NSGA2)求解分布式电源选址定容问题,可作为一个有用的参考,程序注释明确,算法原理可以自己搜。
MATLAB程序采用非支配排序遗传算法(NSGA2)求解分布式电源选址定容问题,可作为一个有用的参考,程序注释明确,算法原理可以自己搜. :8620651507678049浪迹天涯
- 【MATLAB】直方图均衡化Matlab实现(非histeq内置函数)
直方图均衡化Matlab实现(非histeq内置函数) 直方图均衡化常用于图像增强,可以提高图像对比度,扩展图像动态范围. 算法步骤: 统计直方图 累计直方图 映射 close all clear a ...
- 【组合数学】非降路径问题 ( 非降路径问题概要说明 | 非降路径问题基本模型 | 非降路径问题拓展模型 1 非原点起点 | 非降路径问题拓展模型 2 有途经点 )
文章目录 一.非降路径问题 概要说明 二.非降路径问题 基本模型 二.非降路径问题 拓展模型 1 三.非降路径问题 拓展模型 2 组合恒等式参考博客 : [组合数学]二项式定理与组合恒等式 ( 二项式 ...
- 精算模型10 非参数模型0 精算数据、非参数建模的思路概述
精算模型10 非参数模型0 精算数据.非参数建模的目标概述 数据类型 完整数据 非完整数据 数据分析目标 这是精算模型的第十部分,精算的非参数统计模型,从这一部分开始我们讨论精算建模的统计方法.在讨论 ...
- 同步阻塞的BIO、同步非阻塞的NIO、异步非阻塞的AIO
IO的方式通常分为几种,同步阻塞的BIO.同步非阻塞的NIO.异步非阻塞的AIO. 一.BIO 在JDK1.4出来之前,我们建立网络连接的时候采用BIO模式,需要先在服务端启动一个ServerSock ...
最新文章
- 工作随笔——Intellij_idea-14官方快捷键中文版
- QT实现统计历年的太阳黑子数
- aws eks_在生产中配置和使用AWS EKS
- ACM竞赛、数论内容常用的定理(求解(a/b)%c,乘法逆元,费马小定理)
- .net string format
- conceptd什么时候上市_阳山水蜜桃多少钱一斤?什么时候成熟上市?
- 5G手机是否值得买?北京5G基站将超1.4万个
- yum yum doesn‘t match version of Python 终极解决方案
- wireshark网络分析笔记
- DoTween动画系统
- 大型电子病历系统【整套源代码分享】
- 女孩,请把第一胎留给丈夫!!!
- error: Apostrophe not preceded by \ (in XXX.)
- 用python,重温小时候猜数字大小游戏
- 面试官:大量请求 Redis 不存在的数据,从而影响数据库,该如何解决?
- mysql varchar字段自动忽略尾部空格
- 频数直方图的步骤_《频数直方图》PPT课件
- 羽毛球比赛中羽毛球擦网和卡网(挂网)怎么判罚?
- 悦读 | 公布你的原则,读瑞.达利欧的《原则》
- 视频聊天软件开发技术