تصحیح خطاIndex exceeds matrix dimensions

[sama.fa]

باسلام من میخوام تصویر به اندازه16*16بلوک بندی کنم ولی کدی که نوشتم خطامیده فقط به تصاویر خاص مثل تصویر لنا اجرا میشه برای تصویر دلخواه اجرا نمیشه کد و خطا را در اینجا قرار میدهم لطفا راهنمایی کنید ممنون میشم زودجواب بدهید.

کد:



clc
clear
close all
img = imread('1.png');

figure;
subplot(221)
imshow(img)
g = mean2(img);

k=16;
[h,w]=size(img);

n = 1;
for i=1:k:h
for j=1:k:w
imblock=img(i:i+k-1,j:j+k-1);
if nnz(imblock>g)>0
S(n) = sum(imblock(imblock>g));
else
S(n) = k * g;
end
n = n+1;
end
end
subplot(2,2,2)
imshow(reshape(S,h/16,w/16)',[],'InitialMagnification','fit')

خطا:

Index exceeds matrix dimensions

[sama.fa]

در ضمن من این کد برای یه پایگاه داده تصویر می خوام نه یه دونه تصویر.حالا چطور میتونم بلوک بندیش کنم؟؟؟؟

[rahnema1]

سلام
به این صورت تصویر به دست میاد
در نسخه های جدیدتر متلب به جای blkproc از تابع blockproc می شه استفاده کرد

clc
clear
close all
img = imread('1.png');

figure;
subplot(221)
imshow(img)
g = mean2(img);
img_gt_g = img > g;
k=16;
S = blkproc (img_gt_g,[k k], @(x) sum(x(:)));
S(S==0) = k *g;

subplot(2,2,2)
imshow(S,[],'InitialMagnification','fit')

[hhamidy]

سلام
به این صورت تصویر به دست میاد
در نسخه های جدیدتر متلب به جای blkproc از تابع blockproc می شه استفاده کرد

clc
clear
close all
img = imread('1.png');

figure;
subplot(221)
imshow(img)
g = mean2(img);
img_gt_g = img > g;
k=16;
S = blkproc (img_gt_g,[k k], @(x) sum(x(:)));
S(S==0) = k *g;

subplot(2,2,2)
imshow(S,[],'InitialMagnification','fit')


Subscripted assignment dimension mismatch.


Error in blkproc (line 92)
aa(border(1)+(1:ma),border(2)+(1:na)) = a;

[sama.fa]

باسلام ممنون از پاسختون ولی بازم ارور میده

Subscripted assignment dimension mismatch.


Error in blkproc (line 92)
aa(border(1)+(1:ma),border(2)+(1:na)) = a

وبرای بلوک های بعدی که من تو کد قبلی قراردادم بقیه کد قرارداده نشده

[rahnema1]

من که با چند تا عکس مختلف با سایزهای مختلف امتحان کردم ( البته در نرم افزار octave ) مشکلی وجود نداشت اگه بتونید عکس را بفرستید تا امتحان کنم
این برنامه برای تمام بلاکها انجام میده در واقع با حذف حلقه ها کار سریعتر شده

[sama.fa]

مثلا این تصویر ولی کلا رو چندتاتصاویرپایگاه داده میخوام کد اجرا شه


131775

[rahnema1]

در نرم افزار ocatve تست کردم بدون هیچ مشکلی اجرا شد و S هم به این شکل شد

6.0000 155.0000 246.0000 233.0000 72.0000 2003.4363
110.0000 256.0000 256.0000 256.0000 216.0000 2003.4363
128.0000 225.0000 247.0000 227.0000 162.0000 2003.4363
147.0000 140.0000 133.0000 112.0000 55.0000 2.0000
127.0000 191.0000 178.0000 175.0000 106.0000 3.0000
74.0000 139.0000 42.0000 78.0000 11.0000 2003.4363
50.0000 75.0000 92.0000 17.0000 35.0000 2003.4363

[sama.fa]

ممنون ولی من تو نرم افزار متلب میخوام که متاسفانه خطا میده.

ضمنا این کدی که نوشتم ببنین درسته برای بلوک بندی؟؟؟

img = imread('D:\image\roz.png',3);
img = imresize(img,[92,112]);
k=16;
L1 = size(img,1);
L2 = size(img,2);

L1 = k*fix(L1/k);
L2 = k*fix(L2/k);

L1_S = L1/k;
L2_S = L2/k;

img = imresize(img,[L1,L2]);
S = zeros(L1_S,L2_S);
figure(1);
subplot(121)
imshow(img)
g = mean2(img);
[h,w]=size(img);
n = 1;

i1 = 0;j1 = 0;


for i=1:L1_S
for j=1:L2_S
nn_i = (1:k)+(i-1)*k;
nn_j = (1:k)+(j-1)*k;
imblock=img(nn_i,nn_j);
if nnz(imblock>g)>0
S(i,j) = sum(imblock(imblock>g));
else
S(i,j) = k * g;
end
end
end
subplot(122)
imshow(S,[],'InitialMagnification','fit')

[rahnema1]

این کد blkproc در اکتاو هست شما این را در فایلی به نام بذارید و نام blkproc را در برنامه بالا به octblkproc تغییر بدید ببینیم اجرا می شه یا نه.

function B = octblkproc(A, varargin)
if(nargin<3)
error('blkproc: invalid number of parameters.');
end

%% check 'indexed' presence
indexed=false;
p=1;
if(ischar(varargin{1}) && strcmp(varargin{1}, 'indexed'))
indexed=true;
p = p + 1;
if(isa(A,'uint8') || isa(A,'uint16'))
padval=0;
else
padval=1;
end
else
padval=0;
end

%% check [m,n]
if(~isvector(varargin{p}))
error('blkproc: expected [m,n] but param is not a vector.');
end
if(length(varargin{p})~=2)
error('blkproc: expected [m,n] but param has wrong length.');
end
sblk=varargin{p}(:);
p= p + 1;

%% check [mborder,nborder]
if(nargin<p)
error('blkproc: required parameters haven't been supplied.');
end

if(isvector(varargin{p}) && isnumeric(varargin{p}))
if(length(varargin{p})~=2)
error('blkproc: expected [mborder,nborder] but param has wrong length.');
end
sborder=varargin{p}(:);
p = p + 1;
else
sborder=[0;0];
end

%% check fun
%% TODO: add proper checks for this one
if(nargin<p)
error('blkproc: required parameters haven't been supplied.');
end

fun=varargin{p};
if(~isa(fun,'function_handle') &&
~isa(fun,'inline function') &&
~ischar(fun))
error('blkproc: invalid fun parameter.');
end

%% remaining params are params to fun
%% extra params are p+1:nargin-1

%% First of all we calc needed padding which will be applied on bottom
%% and right borders
%% The '-' makes the function output needed elements to fill another
%% block directly
sp=mod(-size(A)',sblk);

if(any(sp))
A=padarray(A,sp,padval,'post');
end

%% we store A size without border padding to iterate later
soa=size(A);

%% If we have borders then we need more padding
if(any(sborder))
A=padarray(A,sborder,padval);
end

%% calculate end of block
eblk=sblk+sborder*2-1;

%% now we can process by blocks
%% we try to preserve fun return type by concatenating everything
for i=1:sblk(1):soa(1)
%% This assures r has the same class as returned by fun
r=feval(fun,A(i:i+eblk(1),1:1+eblk(2)),varargin{p+ 1:nargin-1});
for j=1+sblk(2):sblk(2):soa(2)
r=horzcat(r,feval(fun,A(i:i+eblk(1),j:j+eblk(2)),v arargin{p+1:nargin-1}));
end
if(i==1) %% this assures B has the same class as A
B=r;
else
B=vertcat(B,r);
end
end
end

[sama.fa]

من با نرم افزار ocatve کارر نکردم نمیدونم اصلا چیه؟میشه بیشتر راهنمایی کنید؟؟؟ممنون میشم

ضمنا من میخوام بلوک 16*16 هم پوشان باشه

[rahnema1]

نگفتم از ocatve استفاده کنید گفتم به جای blkproc از این تابعی که در بالا گذاشتم استفاده کنید که اسمش octblkproc هست . همپوشان هم میشه به راحتی انجام داد شما فعلا این تابع را اجرا کنید

[sama.fa]

این کد برای بلوک بندی 16*16هم پوشان صحیح نیست؟؟؟


img = imread('D:\image\roz.png',3);
img = imresize(img,[92,112]);
k=16;
L1 = size(img,1);
L2 = size(img,2);

L1 = k*fix(L1/k);
L2 = k*fix(L2/k);

L1_S = L1/k;
L2_S = L2/k;

img = imresize(img,[L1,L2]);
S = zeros(L1_S,L2_S);
figure(1);
subplot(121)
imshow(img)
g = mean2(img);
[h,w]=size(img);
n = 1;

i1 = 0;j1 = 0;


for i=1:L1_S
for j=1:L2_S
nn_i = (1:k)+(i-1)*k;
nn_j = (1:k)+(j-1)*k;
imblock=img(nn_i,nn_j);
if nnz(imblock>g)>0
S(i,j) = sum(imblock(imblock>g));
else
S(i,j) = k * g;
end
end
end
subplot(122)
imshow(S,[],'InitialMagnification','fit')

[sama.fa]

این کد blkproc در اکتاو هست شما این را در فایلی به نام بذارید و نام blkproc را در برنامه بالا به octblkproc تغییر بدید ببینیم اجرا می شه یا نه.

function B = octblkproc(A, varargin)
if(nargin<3)
error('blkproc: invalid number of parameters.');
end

%% check 'indexed' presence
indexed=false;
p=1;
if(ischar(varargin{1}) && strcmp(varargin{1}, 'indexed'))
indexed=true;
p = p + 1;
if(isa(A,'uint8') || isa(A,'uint16'))
padval=0;
else
padval=1;
end
else
padval=0;
end

%% check [m,n]
if(~isvector(varargin{p}))
error('blkproc: expected [m,n] but param is not a vector.');
end
if(length(varargin{p})~=2)
error('blkproc: expected [m,n] but param has wrong length.');
end
sblk=varargin{p}(:);
p= p + 1;

%% check [mborder,nborder]
if(nargin<p)
error('blkproc: required parameters haven't been supplied.');
end

if(isvector(varargin{p}) && isnumeric(varargin{p}))
if(length(varargin{p})~=2)
error('blkproc: expected [mborder,nborder] but param has wrong length.');
end
sborder=varargin{p}(:);
p = p + 1;
else
sborder=[0;0];
end

%% check fun
%% TODO: add proper checks for this one
if(nargin<p)
error('blkproc: required parameters haven't been supplied.');
end

fun=varargin{p};
if(~isa(fun,'function_handle') &&
~isa(fun,'inline function') &&
~ischar(fun))
error('blkproc: invalid fun parameter.');
end

%% remaining params are params to fun
%% extra params are p+1:nargin-1

%% First of all we calc needed padding which will be applied on bottom
%% and right borders
%% The '-' makes the function output needed elements to fill another
%% block directly
sp=mod(-size(A)',sblk);

if(any(sp))
A=padarray(A,sp,padval,'post');
end

%% we store A size without border padding to iterate later
soa=size(A);

%% If we have borders then we need more padding
if(any(sborder))
A=padarray(A,sborder,padval);
end

%% calculate end of block
eblk=sblk+sborder*2-1;

%% now we can process by blocks
%% we try to preserve fun return type by concatenating everything
for i=1:sblk(1):soa(1)
%% This assures r has the same class as returned by fun
r=feval(fun,A(i:i+eblk(1),1:1+eblk(2)),varargin{p+ 1:nargin-1});
for j=1+sblk(2):sblk(2):soa(2)
r=horzcat(r,feval(fun,A(i:i+eblk(1),j:j+eblk(2)),v arargin{p+1:nargin-1}));
end
if(i==1) %% this assures B has the same class as A
B=r;
else
B=vertcat(B,r);
end
end
end


این تابع در متلب خطا میده

[rahnema1]

این دستور ها را بزنید چه نتیجه ای چاپ می شه؟
img = imread('1.png');
disp(size(img))
disp(class(img))

[sama.fa]

این دستور ها را بزنید چه نتیجه ای چاپ می شه؟
img = imread('1.png');
disp(size(img))
disp(class(img))

من این دستور اجرا کردم اینم نتیجه بعدش چیکار کنم؟؟؟

img = imread('D:\image\roz.png');
disp(size(img))
disp(class(img))
113 91 3

uint8

[rahnema1]

عکس 1 را گفتم همون مرد کچل عینکی

[sama.fa]

اینم مرد کچل خب بعدش چیکار کنم:


img = imread('D:\image\1.png');

disp(size(img))
disp(class(img))
112 92

uint8
>>

[rahnema1]

شما همون برنامه را فکر کنم با roz.png اجرا کرده بودید که خطای Subscripted assignment dimension mismatch. میداد. لطفا با عکس همون مرد کچل ااجراکنید فکر نکنم خطا بده حالا اگه خطا داد احتمالا تصویر rgb بوده که باید به تصویر تک رنگ خاکستری تبدیل بشه
حالا این را امتحان کنید دیگه فکر نکنم خطا بده

clc
clear
close all
img = imread('1.png');
if size(img , 3) == 3
img = rgb2gray(img);
end
img = double(img);
figure;
subplot(221)
imshow(img)
g = mean2(img);
img_gt_g = img > g;
k=16;
S = blkproc (img_gt_g,[k k], @(x) sum(x(:)));
S(S==0) = k *g;

subplot(2,2,2)
imshow(S,[],'InitialMagnification','fit')

[sama.fa]

ممنون از راهنمایی .فقط یه سوال منظور این چند سطر چیه؟؟؟؟

if size(img , 3) == 3
img = rgb2gray(img);
end
img = double(img);

واینکه گفته::

1.میانگین تصویرG
2.برای هر بلوک B از اندازه K*K پیکسل های بزرگتر از عدد میانگین بدست آمده که اول میانگین را حساب بعد جمع میکند و درداخل S قرار می دهد
3.اگر پیکسلی وجود نداشت S=K*G که در آن K ثابت Kدر اینجا بعنوان اندازه بلوک در نطر گرفته میشود
4.S را بعنوان یه ویژگی توی بردار ستونی ذخیره می کنیم
5.مرحله 2تا4را برای بلوک های دیگر هم تکرار میشود


کلیه مراحل که من تو کد قبلی نوشته بودم تواین کد همه این مراحل ذکر نشده ؟؟؟!!

[rahnema1]

تابع blkproc جوری تنظیم شده که فقط با تصاویر خاکستری کار می کنه یعنی تصاویری که دو بعدی هستند
اما تصاویر رنگی در واقع سه بعدی هستند یعنی برای هر رنگ ما تصویر جداگانه داریم
شما سایز را برای هم تصویر خانم و هم مرد کچل زدید اطلاعاتی که چاپ شد برای عکس زن این بود:
113 91 3
یعنی یک عکس 113 در 91 که برای هر رنگ 3 تصویر داریم
اما عکس مرد کچل این جور بود:
112 92
یعنی شامل فقط یک تصویر بود ( یعنی همون چیزی که blkproc انتظار اون را داشت)
حالا با تابع rgb2gray می تونیم تصویر خانم را مثل مرد کچل کنیم که ترکیب سه تایی را به یک عکس تبدیل کنه که متناسب با ورودی blkproc بشه
اون double که زدم شاید اینجا لازم نباشه. ولی در کل بد نیست وقتی می خواهیم روی تصویر عملیات ریاضی انجام بدیم این کار صورت بگیره تا محاسبات دقیق تر بشه
تمام این مراحل که ذکر کردید در این برنامه پیاده کردم. حلقه ها را برداشتم به جای اون blkproc گذاشتم

[sama.fa]

ممنون و سپاس از شما

فقط یه سوال مرحله 4و5 یعنی بردار ستون وتکرار بلوک های دیگر منظوراینه؟؟؟ ((:)sum(x

[rahnema1]

S در واقع در یک ماتریس ذخیره شده نه در یک بردار
بنابراین در پایان همه کدها می تونیم بذاریم که شکل بردار ستونی به خودش بگیره
S_bordar= S(:);

[sama.fa]

در این کد پس مرحله 5 یعنی بلوک های دیگر هم ذکر شده؟؟؟؟؟؟

[rahnema1]

بله محاسبه برای تمام بلوکها صورت گرفته اون هم از طریق blkproc

[sama.fa]

ممنون در مورد فیلتر بالاگذر گوسی سوالمو اینجا مطرح میکنم برای اینکه از تصویر بردارویژگی فرکانس بالا استخراج کنیم من از فیلتر بالاگذر گوسی بردار استخراج کردم اینم کدش :

میشه بگین آیا کد وطرز استفاده ای که کردم درسته یانه؟؟؟؟؟؟؟؟؟

clc;
clear;

img_orginal=im2double(imread('D:\image\roz.png'));
a=img_orginal;
a = rgb2gray(a);

figure,
subplot(221);
imshow(a),title('main image','fontsize',15);

H =fspecial('gaussian' ,[5 5],2);
hpg= imfilter(a,H);

subplot(222);
imshow(hpg),title('matlab gausian filter','fontsize',15);

[m,n]=size(a);
high_filter=zeros(m,n);
f_transform=fft2(a);
f_shift=fftshift(f_transform);
p=m/2;
q=n/2;
d0=70;
for i=1:m
for j=1:n
distance=sqrt((i-p)^2+(j-q)^2);
high_filter(i,j)=exp(-(distance)^2/(2*(d0^2)));
end
end
temp=high_filter.*f_shift;
x=real(ifft2(double(temp)));

subplot(223);
imshow(x,[]),title('my high pass gaussian filter','fontsize',15);

A=conv2(H,x,'same');
[U R]=eig(A)
B=eig(A)

[rahnema1]

این کد را از جایی گرفتید؟ اونجا چی گفته بود؟

[sama.fa]

کد ازیه جایی گرفتم ولی سه سطر اخرشو خودم اضافه کردم برای استخراج بردار ویژگی ولی نمیدونم درست یانه ولی یه کد هم خودم نوشتم که خطا میده ولی اونم دقیق نمیدونم برای استخراج بردار ویژگی هست یانه میشه شما راهنمایی کنید

f=im2double(imread('D:\image\roz.png'));
f = rgb2gray(f);
f=imresize(f,[92 112]);
PQ=paddedsize(size(f));
D0=0.5*PQ(1);
PQ(1)=PQ(2);
H=hpfilter('gaussian',PQ(1),PQ(2),D0);
g=dftfilt(f,H);
figure,imshow(g,[])

[rahnema1]

حالا من یه بررسی می کنم

[sama.fa]

ممنون ازشما

[rahnema1]

چه کسی گفته لازمه از فیلتر بالاگذر گوسی استفاده بشه