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

特征粗糙度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');

%读新图像文件名
   
    %图像裁剪
   ******
    
    %找每幅图像的最大值和最小值
%%%%%%%%%计算每张差图像的分形维数%%%%%%%%%%%%%%%%%%
******
D=D'
R=R'
Ra=Ra'