一、皮尔逊相关系数

在统计学中，皮尔逊相关系数( Pearson correlation coefficient），又称皮尔逊积矩相关系数（Pearson product-moment correlation coefficient，简称 PPMCC或PCCs），是用于度量两个变量X和Y之间的相关（线性相关），其值介于-1与1之间。

计算步骤参考下面文章的第七节相关性

机器学习笔记 - 数据和统计常用术语_bashendixie5的博客-CSDN博客从以下角度描述了数据及其统计方法。1、数据的类型定性与定量、离散或连续2、集中趋势的度量均值、中位数、众数、异常值、几何平均值、调和平均值、加权平均值。3、数据的采样随机采样、系统采样，分层采样，聚类样本等。4、数据分布的度量分位数、平均偏差、标准差和方差5、数据的比较单变量和双变量、离群值、相关性概念及其公式https://blog.csdn.net/bashendixie5/article/details/123613591

二、matlab计算相关性

1、示例1

按公式计算，A 和 B 高度相关。相关系数为1。

A= [1 4 7; 2 5 8; 3 6 9]
B = A*2%Find the average of the matrix A
meanA = mean2(A);%Find the average of the matrix B
meanB = mean2(B);%Subtract the average value from matrix A
Asub = A-meanA;%Subtract the average value from matrix B
Bsub = B-meanB;%Covariance of matrix A and matrix B
covAB = mean2(Asub.*Bsub);%Find the standard deviation of the matrix A
stdA = std(A(:),1);%Find the standard deviation of the matrix B
stdB = std(B(:),1);%Find the correlation Cofficient
Rho = covAB./(stdA*stdB)

2、示例2

使用corr2函数，C和D相关系数为 -1。

C= [1 4 7; 2 5 8; 3 6 9]
D = [9 6 3;8 5 2; 7 4 1];Rho = corr2(C,D)

三、模板匹配示例

1、使用normxcorr2函数

参考代码

%Read an Image A(Template)
A1 = imread('benten.jpg');%Read the Target Image
B1 = imread('watch.jpg');A = A1(:,:,1);
B = B1(:,:,1);normx_corrmap=normxcorr2(B(:,:,1),A(:,:,1));maxptx = max(normx_corrmap(:));
[x1,y1]=find(normx_corrmap==maxptx);
figure,
imagesc(A1(x1-size(B,1):x1,y1-size(B,2):y1,:));axis image

运行结果如下

2、使用corr2函数

参考代码

%Read an Image A(Template)A1 = imread('图片/benten.jpg');%Read the Target ImageB1 = imread('图片/watch.jpg');A = A1(:,:,1);B = B1(:,:,1);corr_map = zeros([size(A,1),size(A,2)]);for i = 1:size(A,1)-size(B,1)for j = 1:size(A,2)-size(B,2)%Construct the correlation mapcorr_map(i,j) = corr2(A(i:i+size(B,1)-1,j:j+size(B,2)-1),B);endendfigure,imagesc(corr_map);colorbar;%Find the maximum valuemaxpt = max(corr_map(:));[x,y]=find(corr_map==maxpt);%Display the image from the templatefigure,imagesc(B1);title('Target Image');colormap(gray);axis imagegrayA = rgb2gray(A1);Res   = A;Res(:,:,1)=grayA;Res(:,:,2)=grayA;Res(:,:,3)=grayA;Res(x:x+size(B,1)-1,y:y+size(B,2)-1,:)=A1(x:x+size(B,1)-1,y:y+size(B,2)-1,:);figure,imagesc(Res);

运行结果如下

3、频域模板匹配

基于傅里叶域中的归一化互相关。也称为相位相关。这里使用的两张图片是同一场景的不同快照。Image1.jpg用作模板图像，Image2.jpg的子图像用作目标图像。目标图像用零填充以匹配模板图像的大小。经过傅里叶变换后，模板信号与目标信号的共轭相乘并归一化。然后应用逆傅里叶，提取最大值对应的像素位置。

%Read two images of same scene
A = imread('Image1.jpg');
B = imread('Image2.jpg');figure,subplot(2,1,1);imagesc(A);title('Image 1');axis image
subplot(2,1,2);imagesc(B);title('Image 2');axis image

从图像矩阵B中裁剪一部分

B = imcrop(B,[58.5 49.5 226 102]);
figure,imagesc(B);title('sub Image - Image 2');axis image

应用傅里叶变换和逆变换

%Pad the image matrix B with zeros
B1 = zeros([size(A,1),size(A,2)]);
B1(1:size(B,1),1:size(B,2))=B(:,:,1);%Apply Fourier Transform
Signal1 = fftshift(fft2(A(:,:,1)));
Signal2 = fftshift(fft2(B1));%Mulitply Signal1 with the conjugate of Signal2
R = Signal1 .*conj(Signal2);%Normalize the result
Ph = R./abs(R);%Apply inverse fourier transform
IFT = ifft2(fftshift(Ph));figure,imagesc((abs((IFT))));colormap(gray);

查找最大值的像素位置，并从image1和image2中截取

%Find the maximum value
maxpt = max(real(IFT(:)));%Find the pixel position of the maximum value
[x,y]= find(real(IFT)==maxpt);figure,subplot(1,2,1);imagesc(A(x:x+size(B,1),y:y+size(B,2),:));axis image
subplot(1,2,2);imagesc(B);axis image

结果如下

MATLAB学习笔记皮尔逊相关系数和模板匹配相关推荐

OpenCV4学习笔记（19）——模板匹配
本次要整理记录的内容是:模板匹配. 模板匹配可以说是一种最简单的模式识别方法,它的实现主要是通过模板图像在被匹配图像中的平移,在被匹配图像中逐个区域寻找和模板图像相似的区域,如果存在某区域的相似度大于 ...
【opencv学习笔记】028之模板匹配——matchTemplate函数详解
目录一.前言二.模板匹配 1.模板匹配是个啥 2.常用匹配算法 1.平方差匹配-CV_TM_SQDIFF 2.标准平方差匹配-CV_TM_SQDIFF_NORMED 3.相关匹配-CV_TM_CC ...
opencv学习笔记十八：模板匹配(cv2.matchTemplate、cv2.minMaxLoc)
1.目标匹配函数:cv2.matchTemplate() res=cv2.matchTemplate(image, templ, method, result=None, mask=None) ima ...
数字图像处理MATLAB学习笔记（五）
数字图像处理MATLAB学习笔记(五) Color Image Processing 1 Color Image Representation in MATLAB 这里不多说了,彩色图片在计算机中以R ...
MATLAB学习笔记（二）
MATLAB学习笔记(二) 一.矩阵运算矩阵分析向量和矩阵的范数运算矩阵的秩矩阵的化零矩阵矩阵的化简rref()函数线性方程组超定线性方程组求解矩阵分解 1.对称正定矩阵的Choles ...
matlab 线性规划_从零开始的matlab学习笔记——（37）线性规划——后传
matlab应用--求极限,求导,求积分,解方程,概率统计,函数绘图,三维图像,拟合函数,动态图,傅里叶变换,随机数,优化问题....更多内容尽在个人专栏:matlab学习呼,终于将matlab的优 ...
MATLAB学习笔记（十八）
MATLAB学习笔记(十八) 一.Simulink仿真基础 1.1 Simulink的启动 1.2 系统仿真模型的创建 1.3 仿真参数的设置 1.4 总结二.子系统的创建与封装 2.1 子系统的创 ...
Matlab学习笔记5——fprintf
Matlab学习笔记5--fprintf 将数据写入文本文件语法 fprintf(fileID,formatSpec,A1,-,An) fprintf(formatSpec,A1,-,An) nby ...
Matlab学习笔记(3)—GUI程序设计与图像处理基本操作
Matlab学习笔记(3)-GUI程序设计与图像处理基本操作一.GUI程序设计如果想要开始使用Matlab进行GUI编程,首先需要进入Matlab的GUI界面对GUI有一个基本的认识. 1.gui ...

MATLAB学习笔记皮尔逊相关系数和模板匹配

一、皮尔逊相关系数

二、matlab计算相关性

1、示例1

2、示例2

三、模板匹配示例

1、使用normxcorr2函数

2、使用corr2函数

3、频域模板匹配

MATLAB学习笔记皮尔逊相关系数和模板匹配相关推荐

最新文章

热门文章

MATLAB学习笔记 皮尔逊相关系数和模板匹配

一、皮尔逊相关系数

二、matlab计算相关性

1、示例1

2、示例2

三、模板匹配示例

1、使用normxcorr2函数

2、使用corr2函数

3、频域模板匹配

MATLAB学习笔记 皮尔逊相关系数和模板匹配相关推荐

最新文章

热门文章

MATLAB学习笔记皮尔逊相关系数和模板匹配

MATLAB学习笔记皮尔逊相关系数和模板匹配相关推荐