clc;clear;close?all
im=imread（‘1.jpg‘）;
im=rgb2gray（im）;
figure
imshow（im）??????????%企鵝灰度圖像

[mn]=size（im）;?????%原圖的尺寸
blocksize=8;????????%分塊大小，此處塊為正方形，方便后面差值運算
%分塊處理
blocknum1?=floor（m/blocksize）;????%分塊數量（對應圖像長度，即m），取整后，原圖會剩下幾個像素的長或寬未做處理
blocknum2?=floor（n/blocksize）;????%分塊數量（對應圖像寬度，即n）
length?=blocknum1*blocksize;???%處理的長度像素數
height?=blocknum2*blocksize;???%處理的寬度像素數
A=zeros（blocknum1blocknum2）?;
ff?=?im（1:length1:height）;?????%生成一個用來原圖副本（去除多余像素的）
%開始分塊處理
for?k?=?1:blocknum2
????for?h?=?1:blocknum1
????
????????%生成模板
????????block?=?zeros（size（ff））;????%將模板初始化為0
????????lini?=?1?+?blocksize?*?（h?-?1）;
????????hini?=?1?+?blocksize?*?（k?-?1）;?????%分塊的第一個像素在原圖中的坐標
????????x?=?lini:（lini?+?blocksize?-?1）;
????????y?=?hini:（hini?+?blocksize?-?1）;????%生成分塊長、寬坐標序列
????????block（x?y）?=?1;?????????????????????%將模板上需要進行分塊的部分轉換成1，用來提取該分塊
????????
????????%提取分塊
????????ff?=?im2double（ff）;?????%原圖轉換