How to solve the problem?

1 次查看（过去 30 天）

Darsana P M 2018-4-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/397441-how-to-solve-the-problem

评论： Darsana P M 2018-4-27

I have a 198x198 matrix, an image. I need to get 3x3 matrices each from the above matrix. Total, 22 matrices, how to access each of them? can somebody help me with the code??

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Walter Roberson 2018-4-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/397441-how-to-solve-the-problem#answer_317252

在 MATLAB Online 中打开

blocks = mat2cell(YourMatrix, 3*ones(1,22), 3*ones(1,22));

then you would access (for example) blocks{7,4) to give you the 7th block of 3 down, 4'th block of 3 across.

17 个评论
显示 15更早的评论隐藏 15更早的评论

Darsana P M 2018-4-27

在 MATLAB Online 中打开

clc;
close all;
clear all;
k=0;
I = dicomread('C:\Users\Click Me\Desktop\NIELIT\Output images\dic5.dcm');
info = dicominfo('C:\Users\Click Me\Desktop\NIELIT\Output images\dic5.dcm');
figure(1);
image(I);
h=rgb2gray(I);
imhist(h);
N= 257;
l=zeros(200,400);
B_mat=h(1:200,:);
figure(2);
image(B_mat);
%Rot=[-15,-20,60,20;52,-39,0,25;36,48,32,13;12,43,23,76];
Rot=[-15,-20,60;52,-39,0;36,48,32]
o=[3 4 2 6];
q1=o;
q2=qInv(q1);
L0=B_mat(:,1:200);
R0=B_mat(:,201:400);
h=qGetR(R0);
[M,N]=size(L0);
b=3;
M=M-rem(M,b);
N=N-rem(N,b);
L0=L0;
%A=zeros(200,200);
lo=L0(1:200,1:200);
%s=size(lo);
% dim=198;
% [row column]=size(lo);
% count=0;
% Bigblock=zeros(198,198);
[r, c] = size(L0);
temp = L0;
temp( r+1:ceil(r/3)*3, c+1:ceil(c/3)*3 ) = 0;
[r, c] = size(temp);
blocks = mat2cell(temp, 3*ones(1,r/3), 3*ones(1,c/3));
%J=cell2mat(blocks);
clear temp;
this_block = blocks{7,4}; 
%extract 3x3 matrix
n1=(this_block).*Rot;
         L1=mod(n1,N);
        %B=cat(dim,Bigblock,block);
        R1=(L1 + this_block);
         m=mod(R1,N);
         %new_block=this_block;
%new_block=mat2cell(blocks);
%new_block = reshape(sort(this_block(:)), size(this_block));  %perform certain operations
new_blocks{7,4} = new_block;
new_matrix = cell2mat(new_blocks);
%new_blocks= cell2mat(blocks);
new_matrix = new_matrix(1:size(L0,1), 1:size(L0,2));  %trim out padding