GI 鬼成像

GI 鬼成像源代码，关联算法和一些数据读取存储等问题 的解决方案

%取一个缝
clear all;    %  清除变量
clc;    %  清除屏幕
close all;    %  关闭图片

Dir1 = ‘D:\backup\experdata\08122201-1\‘;    %  文件夹路径
Directory1 = ‘D:\backup\experdata\08122201-1\10000.jpeg‘;    %  第一幅图像的路径，以便测试图像大小
Dir2 = ‘D:\backup\experdata\08122201-2\‘;    %  文件夹路径
Directory2 = ‘D:\backup\experdata\08122201-2\10000.jpeg‘;    %  第一幅图像的路径，以便测试图像大小
Num =10000;    %  图片文件数

Image1 = imread（Directory1 ‘jpeg‘）;    %  读取第一幅图像
Image2 = imread（Directory2 ‘jpeg‘）;    %  读取第一幅图像
[Row Col] = size（Image1） %  得到图像的行列数
[Row Col] = size（Image2）
ImageSump = zeros（500500）;   %  生成一个零值矩阵，以便图像求和
Im = zeros（500500）;   %  生成一个零值矩阵，以便图像求和
Ima=zeros（450480）;
ImageSumpp = zeros（450480）;
Colg2 = zeros（500500）;    %  g2的存储矩阵
ColG2 = zeros（500500）;    %  G2的存储矩阵
I2 = zeros（1 1）;  
for I = 1 : Num    %  循环将NUM幅图像累加
    ImageI=9999+I;
    FileName =  [int2str（ImageI） ‘.jpeg‘];    %  形成文件名称，如Image1.jpg
    Image1 = imread（[Dir1 FileName] ‘jpeg‘）;    %  读取图像 
    Image2 = imread（[Dir2 FileName] ‘jpeg‘）;    %  读取图像 
    Im =  Image1（121:620321:820）;
    Ima = Image2（21:47081:560）;
    ImageSump = ImageSump + double（Im）; %  图像求和
    ImageSumpp = ImageSumpp + double（Ima）;
    SumI2 = sum（sum（double（Ima）‘））;
    I2 = I2 + SumI2;
    %ColG2 = ColG2 + double（Im） * SumI2 ;
end
ImageAvep = ImageSump / Num; %  求得平均图像
ImageAvepp = ImageSumpp / Num; %  求得平均图像

ImageNormalizep = ImageAvep / max（max（ImageAvep））; %  图像最大值
ImageNormalizepp = ImageAvepp / max（max（ImageAvepp））;
ImageGrayPicp = uint8（round（255 * ImageNormalizep））;%  归一化图像，并投影到0-255的灰度区间
ImageGrayPicpp = uint8（round（255 * ImageNormalizepp））; 
figure;
imshow（ImageGrayPicp）;    %  显示平均图像
figure;
imshow（ImageGrayPicpp）;    %  显示平均图像

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1926  2010-10-17 23:24  GI.m

----------- ---------  ---------- -----  ----

                 1926                    1