Matlab图像处理与分析之压力分布图像获得

（若有侵权，请立刻联系笔者删除）
各位兄弟姐妹们好久不见，由于老板一直在催促着画图所以一直没有更新相关学习内容，此次带来的偏向于对课程实验结果的一次图像处理，测量的是压力分布，谈到压力测量可能我们第一想到的是水银计，高中还是初中学过水银的高度反应气压大小的实验，当然此次实验用的是压敏漆，简称PSP。此次会大体介绍一下PSP并重点介绍一下结果图像在MATLAB中的大致处理过程。
1.什么是PSP？
PSP是压敏漆英文单词的缩写，我们知道飞机飞行时需要得到表面所受力的分布，根据所得结果可以进一步改进飞机结构提高安全性和稳定性，但在高速飞行中或者说实验室气流特别猛烈地情况下，直接接触式测量就不太能用了，且测量范围受限，因此需要一种非接触式测量，由此，PSP就出现了。
将粘合剂涂抹到表面后，会产生一定的光物理现象，总的而言就是气压增加时，发光强度会降低，具体的光学原理可以探索一下Stern-Volmer原理，这个原理笔者也没有太多研究，此次实验更偏向于对装置的运用和结果的分析，大体运用的公式便是：

这个公式中间项没有考虑，左侧Iref是指未通风的结果，I是指通风后的结果，左边就是用通风后的结果除以未通风的结果，右边的AB都是和不同的压敏漆相关的系数，右边的pref则是测量时的标准气压，而P就是我们要求的值。获得左边的两幅图像不是难事，做好加工调整好光源就获取图像就好。
2.图像压力测量
本次的实验条件是环境20度，气压10万Pa左右，所获的参数在程序中显示，先看一下获得的原图：(下面的变量名有些尽管名称相同，但不一定会有相同的作用，希望各位主要理解函数的功能和一般图像处理方向)

opentif=imread('open.tif');%读取图像

（1）图像预处理：
最左面的是没有风和光源加入的情况下拍摄的图片，为了获得相机本身因素所造成的噪声，因此用后两张图，减掉这张背景图即可去除背景噪声。所用的函数为imsubtract。

Iref=imsubtract(off,background);%两者图像相减

得到的结果为：

其实对于上述的图片，我们可能肉眼看不出变化，但是有一个工具可以获得，就是直方图，直接上代码：

[n1,histon]=imhist(opentif);%opentif是定义在程序中的图像名称，n1和histon分别是图像的像素值和分布大小

最左边的直方图是背景噪音图的直方图，中间的图片是未去噪的图，右边的图是去噪后的直方图，可以明显看到相减后的效果是很好的。下面进行阈值分割，根据直方图选取阈值可以轻松的分离三角区和背景图像。结果不再展示因为背景分割后都归0也都是黑色的。
（2）图像配准

[optimizer,metric]=imregconfig('multimodal');                            %图像对齐的信息和处理方式
optimizer.InitialRadius = optimizer.InitialRadius/2;                     %改变步长
optimizer.MaximumIterations=300;                                         %改变迭代次数
movingRegisteredDefault=imregister(I,Iref,'affine',optimizer,metric);    %对准函数，采用的仿射变换模式
figure; imshowpair(movingRegisteredDefault, Iref);                       %展示两幅图像不同

代码给出后其实整个的过程也就很理解了，首先就是第一行代码对应的函数选取一种两幅图像的比较方式和对齐方式，然后通过改变步长和迭代次数后，通过第四行的代码进行配准，图像配准，其实也就是让两幅图像在大体轮廓位置上更加相似。最后一行代码是比较两幅图像不同，可以修改的参数都是单引号括起来的内容以及文字信息带改变二字 的内容。读者可以通过help查看不同参数的意思来尽情修改尝试。

不知道各位有没有看出来有没有什么变化，反正将第一幅图照着第二幅图配准后的结果肉眼我是看不太出来的。
（3）压力分布测试
这个过程跟PSP原理有关，将所提及的公式写成代码即可。然后得到的结果可能还需要进行中值滤波：

opentif=medfilt2(opentif,[5,5]);%opentif代表的是计算后的压力分布图像，后面的是滤波模板的大小

对于展示的结果可以通过颜色可视化：

figure;imagesc(opentif);
caxis([1.005e+05,1.03e+05]); %加入颜色条并限制范围
colorbar;

所得到的结果：

还是比较具有观赏性的吧，好了，今天就到这啦，有关图像处理的内容我可能了解的不是特别多，也就是传统的算法不断寻找然后集合起来，有需要的可以一起交流。最后还是那句话：如果作者有任何侵权的现象或者需要删除的内容请立刻联系，谢谢！

参考文献
1.MATLAB图像配准相关知识