Matlab计算热像图序列的特征粗糙度

特征粗糙度Ra是衡量粗糙度大小的量。公式如下。

其中，k为双对数坐标轴中的纵轴截距，D为分形维数。Ra将分形维数D和尺度系数k的作用同时进行了考虑，能够唯一表征热像温度场的粗糙度。当热像温度场温度分布越光滑，Ra越小，当温度场分布越粗糙，Ra越大。

1. 文件、数据介绍。

文件夹下有55个格式为mat的热像图数据。

每个mat文件包含一个200*300的juzhen变量。

2.运行程序，结果如下。

分形维数

特征粗糙度

3.部分源代码。

clc 
clear all

%%%%%%%%找到图片序列的最大值、最小值%%%%%%%%%%%%%%
%读取照片名
fid_file=fopen('E:\实用程序\热像\东北大学\k2-1-处理\k2-1-mat\k2-1-matname.txt');

%读新图像文件名
   
    %图像裁剪
    %picture=temp_picture;
    picture=temp_picture;
    %将差图像保存出来
    pic{k}= picture;
    
    %找每幅图像的最大值和最小值
%%%%%%%%%计算每张差图像的分形维数%%%%%%%%%%%%%%%%%%
for i=1:(k-1)
   photo=pic{i};
   %求图像大小
   [M N]=size(photo);
      %求最大灰度级G
   G=max_value-min_value;
   %确定不同子块长度L
   %拟合结果及相关系数
   %计算特征粗糙度
  %分形维数，相关系数，特征粗糙度
D=D'
R=R'
Ra=Ra'