基于MATALB编程的异物检测（代码完整，数据齐全）

function main() clc close all [filename,pathname]=uigetfile(... {'*.jpg;*.png;*.bmp;*.tif;*.jpeg',... 'Image Files(*.jpg,*.png,*.bmp,*.tif,*.jpeg)'; '*.jpg','JPEG(*.jpg;*.jpeg)';... '*.png','PNG(*.png)';... '*.bmp','BMP(*.bmp)';... '*.tif','TIF(*.tif)';... '*.*','All Files(*.*)'},... 'Pick a file'); img=imread([pathname,filename]); %m:行 n：列 [m,n,~]=size(img); %得到三个矩阵B，C,D B=img(:,:,1); C=img(:,:,2); D=img(:,:,3); % 为了确保输入的数值a能同时被m和n整除，需先求取m和n的公约数 % findCommomDivisior : 求公约数 % commonDivisor ：m和n的公约数 disp('公约数如下，请输入其中一个公约数完成后续处理') commonDivisor = findCommonDivisor( m, n ) a=input('输入的数值为：'); % 判断输入的数值a是否是公约数中的一个 % 若不是，提示输入错误；若是，则进行后续处理 tmp=any( a == commonDivisor(:) ); if (~tmp) disp('输入了错误的数值'); else % constructMat: 构造所需要的E，F，G矩阵 [b, c, E, F, G] = constructMat(m, n, a, B, C, D); % E % F % G % 保存E，F，G矩阵为.mat格式，利于后续分析(相对路径) save ('.\保存的数据\E,F,G矩阵\myDateE.mat','E') save ('.\保存的数据\E,F,G矩阵\myDateF.mat','F') save ('.\保存的数据\E,F,G矩阵\myDateG.mat','G') t=input('辨识参数(该参数为下限值0-100之间)：'); % analyse:分析函数 [indexE , indexF , indexG] = analyse(t,b,c,E,F,G); disp('E矩阵中异常元素行（row）列（col）值'); indexE disp('F矩阵中异常元素行（row）列（col）值'); indexF disp('G矩阵中异常元素行（row）列（col）值'); indexG % 保存E,F,G矩阵中异常元素的行列值(相对路径) save ('.\保存的数据\E,F,G矩阵中异常元素的行列值\myIndexE.mat','indexE') save ('.\保存的数据\E,F,G矩阵中异常元素的行列值\myIndexF.mat','indexF') save ('.\保存的数据\E,F,G矩阵中异常元素的行列值\myIndexG.mat','indexG') end