趁热打铁！
学习图像剪切。


matlab源码
function [imout] = ImageCrop( imin,x,y,width,height )

%UNTITLED13 此处显示有关此函数的摘要

%   此处显示详细说明



[m,n,g] = size(imin);





imin = im2double(imin);

imout = zeros(width,height,g);





for i= 1:width

    for j = 1:height

        if (i+x <= m && j+y <=n)

        imout(i,j,:) = imin(i+x,j+y,:);

        end

    end

end



end




图像剪切函数测试代码


clc;

clear;



i = imread('g:\cool\data\1.bmp');

figure('name','图像剪切测试');

subplot(121);imshow(i);title('原始图片');



out = ImageCrop(i,100,100,200,250);



subplot(122);imshow(out);title('剪切后图片');

Matlab图像剪切，超出图像区域指定颜色填充
  【尊重原创，转载请注明出处】http://blog.csdn.net/guyuealian/article/details/77507969
   Matlab自带图像剪切（裁剪）函数，即imcrop函数，但该函数对于超出图像边界区域是停止裁剪的，即对超出图像边界的部分不作裁剪处理。正常情况下，imcrop函数足够使用了。但有时，我们希望对于超出图像边界的区域，可以指定颜色填充，以便我们可以观察裁剪情况。
   若想在OpenCV中实现图像剪切，对超出图像区域指定颜色填充，可参考本人的博客：http://blog.csdn.net/guyuealian/article/details/77486096
   本人实现方法如下：
function destImg=ImageCropPadding(gray_img,rect,paddingValue)
% ImageCropPadding函数功能说明：
% 该函数可实现根据rect(x,y,w,h)剪切指定的灰度图像，若剪切超出图像边界，则用paddingValue的值填充
% 输入参数：gray_img：必须输入灰度图像，暂时不支持RGB图像；
%              rect：指定剪切区域，类似于Matlab的imcrop函数用法
%      paddingValue：填充值（1-255）
%% For Example:
%  SrcImage=imread('images/liuyifei.jpg');
%  grabImage=rgb2gray(SrcImage);
%  figure,imshow(grabImage),title('grabImage');
%  rect=[-50,-50,500,500];
%  paddingValue=1;
%  destImg=ImageCropPadding(SrcImage,rect,paddingValue);
%  figure,imshow(destImg),title('destImg');
%%
x=rect(1);
y=rect(2);
w=rect(3);
h=rect(4);
destImg = im2uint8(zeros(h, w)+paddingValue/255); 
crop_x1 = max(1, x);
crop_y1 = max(1, y);
[rows,cols]=size(gray_img);
crop_x2 = min(cols, x + w); 
crop_y2 = min(rows, y + h);
roi_img = gray_img(crop_y1:crop_y2,crop_x1:crop_x2);
x1 = crop_x1 - x+1;
y1 = crop_y1 - y+1;
x2 = crop_x2 - x+1;
y2 = crop_y2 - y+1;
% I(y1:y2, x1:x2 ) = I(y1:y2 , x1:x2) + roi_img;
destImg(y1:y2, x1:x2 ) =  roi_img;
end
   调用方法：
SrcImage=imread('images/liuyifei.jpg');
grabImage=rgb2gray(SrcImage);
figure,imshow(grabImage),title('grabImage');
rect=[-50,-50,500,500];
paddingValue=1;%黑色填充
destImg=ImageCropPadding(SrcImage,rect,paddingValue);
figure,imshow(destImg),title('destImg');
   运行效果：

matlab按步长剪切图片：
clear;close all;
%% settings
folder = 'E:\train\';
savepath = 'F:\matlab\';
size_input = 60;
stride = 40;
%% initialization
data = zeros(size_input, size_input, 1, 1);
count = 0;
%% generate data
filepaths = dir(fullfile(folder,'*.bmp'));
for i = 1 : length(filepaths)
image = imread(fullfile(folder,filepaths(i).name));
image = rgb2ycbcr(image);
image = im2double(image(:, :, 1));
[hei,wid] = size(image);
for x = 1 : stride : hei-size_input+1
for y = 1 :stride : wid-size_input+1
count=count+1;
subim_input = image(x : x+size_input-1, y : y+size_input-1);
imwrite(subim_input,strcat(savepath,num2str(count),'.bmp'),'bmp');
end
end
end
重命名：
clear;close all; %% settings folder = 'D:\SYM\Caffe\caffe-master\mypro\GAN\data\train\'; savepath = 'F:\matlab\'; %% initialization count = 2386; %% generate data filepaths = dir(fullfile(folder,'*.jpg'));    for i = 1 : length(filepaths)      image = imread(fullfile(folder,filepaths(i).name));     count=count+1;             imwrite(image,strcat(savepath,num2str(count),'.jpg'),'jpg'); end

读入实验图像EXP1A.tif和EXP1B.bmp，分别剪切出上述图像右上角四分之一的内容，并要求在同一个图形界面上显示。

要求：

在同一个图形界面上，于第1行第1列位置上显示图像EXP1A.tif的原始内容，于第1行第2列位置上显示EXP1A.tif右上角四分之一的内容，于第2行第1列位置上显示图像EXP1B.bmp的原始内容，于第2行第2列位置上显示EXP1B.bmp右上角四分之一的内容。

实现思路：

1.使用imread函数依次读取图像'EXP1A.tif'和图像'EXP1B.bmp'

2.使用size函数读取图像的行数r和列数c，将两幅图像的行数和列数分别除以2：r/2、c/2，因为运算结果可能不是整数，所以需要使用ceil函数来对结果进行取整得到二分之一行数half_x和列数half_y

3.分别取两幅图像的1至half_x行和half_y至最后一行

4.按题目要求显示图像

实现效果：



实现源代码：

clc;
%读取图像
img_a = imread('EXP1A.tif');
img_b = imread('EXP1B.bmp');


[ar,ac,ak] = size(img_a);%取图像的行数r和列数c
%取图像右上角四分之一
half_ax = ceil(ar/2);%ceil是对r/2的结果取整(向正无穷方向)
half_ay = ceil(ac/2);
img_a1 = img_a(1:half_ax,half_ay:end);

[br,bc,bk] = size(img_b);%取图像的行数r和列数c
%取图像右上角四分之一
half_bx = ceil(br/2);%ceil是对r/2的结果取整(向正无穷方向)
half_by = ceil(bc/2);
img_b1 = img_b(1:half_bx,half_by:end,:);

%显示图像
figure('NumberTitle', 'off', 'Name', '实验一第二题'); 
subplot(2,2,1);
imshow(img_a);
title('EXP1A.tif 原图');

subplot(2,2,2);
imshow(img_a1);
title('EXP1A.tif 四分之一右上角图');

subplot(2,2,3);
imshow(img_b);
title('EXP1B.bmp 原图');

subplot(2,2,4);
imshow(img_b1);
title('EXP1B.bmp 四分之一右上角图');