i-frame

[zeinab 12]

سلام.من دارم روی پروژه i-frame کار میکنم که درواقع هم روی خود کد و هم بهینه سازی کد مشکل دارم.اگر کسی بتونه کمکم کنه واقعا ممنون میشم.
این i-frame شامل مراحل زیر است:
1.خواندن یک تصویر رنگی و تبدیل به grayscale
2.شکستن تصویر grayscale به بلوکهای 8*8
3.تبدیل dct روی بلوکها
4.کوانتایز کردن ضرایب dct طبق جول کوانتایزسیون استاندارد با ضریب q=1
5.کوانتایز معکوس
6.تبدیل dct معکوس
7.تشکیل تصویر grayscale
8.نمایش تصویر grayscale اولیه به دست آمده
9.محاسبه MSE بین تصویر grayscale به دست آمده و اولیه

code:
close all
clear all
I=imread ('cameraman.tif');
J=rgb2gray(I);
I2=blkproc(J,[8 8],dct2);
Q=1;
F1=Q*round(I2,Q);
F1(1,1)=F(1,1);
G=idct2(F1);
figure;
subplot(2,2,1),imshow(J);
subplot(2,2,2),imshow(G,[0,255]);
subplot(2,2,3),imhist(J);
subplot(2,2,4),imhist(G);

[rahnema1]

close all

clear all

img = imread ('cameraman.tif');

grayImg =double(rgb2gray(img));

DC = blkproc(grayImg, [8 8], @dct2);

zarib = 1;

Quantization_Matrix=...

	[ 16  11  10  16  24  40  51  61; 

	  12  12  14  19  26  58  60  55; 

	  14  13  16  24  40  57  69  56; 

	  14  17  22  29  51  87  80  62; 

	  18  22  37  56  68 109 103  77; 

	  24  35  55  64  81 104 113  92; 

	  49  64  78  87 103 121 120 101; 

	  72  92  95  98 112 100 103  99 ];

if zarib < 50

	Quantization_Matrix = Quantization_Matrix * 50 / zarib ;

else 

	Quantization_Matrix = Quantization_Matrix * (100 - zarib) / 50;

end

Quantization_Matrix = double(uint8((Quantization_Matrix - .2) / 254.8 * 254 + 1));

quantized_DC = round(blkproc(DC, [8 8], @(x, y) x ./ y , Quantization_Matrix));

reconstructed_DC = blkproc(quantized_DC, [8 8], @(x, y) x .* y , Quantization_Matrix);

compressed_image = round(blkproc(reconstructed_DC, [8 8], @idct2));

mse =sum(( grayImg(:) - compressed_image(:)) .^ 2) / numel(grayImg);

printf("MSE = %f\n", mse)

figure;

subplot(2,2,1),imshow(grayImg,[0,255]);

subplot(2,2,2),imshow(compressed_image,[0,255]);

subplot(2,2,3),imhist(grayImg);

subplot(2,2,4),imhist(compressed_image);

[zeinab 12]

سلام ببخشید اما میشه راجع به کدتون توضیح بدین چون من تازه شروع کردم خیلی درک بالایی از کدهاتون ندارم.ممنون

[zeinab 12]

سلام میشه قسمت های زیر رو توضیح بدید؟
if zarib < 50
Quantization_Matrix = Quantization_Matrix * 50 / zarib ;
else
Quantization_Matrix = Quantization_Matrix * (100 - zarib) / 50;
end
Quantization_Matrix = double(uint8((Quantization_Matrix - .2) / 254.8 * 254 + 1));
quantized_DC = round(blkproc(DC, [8 8], @(x, y) x ./ y , Quantization_Matrix));
reconstructed_DC = blkproc(quantized_DC, [8 8], @(x, y) x .* y , Quantization_Matrix);
compressed_image = round(blkproc(reconstructed_DC, [8 8], @idct2));
mse =sum(( grayImg(:) - compressed_image(:)) .^ 2) / numel(grayImg);

[rahnema1]

این pdf را بخونید داخلش توضیحات خوبی ارائه داده
http://www.lokminglui.com/dct.pdf
اون عبارت = quantized_DC یعنی DC / Quantization_Matrix
و عبارت = reconstructed_DC هم یعنی quantized_DC * Quantization_Matrix البته به صورت خانه های 8*8
اون = compressed_image هم که تبدیل معکوس dct هست
mse که MSE را به دست میاره!

[zeinab 12]

ممنون.یک سوال دیگه هم دارم چرا وقتی تصویر رنگی رو به خاکستری تبدیل کردید با عنوان داده ای double در grayImg ذخیره کردید؟

[rahnema1]

سلام
بهتره برای انجلم محاسبات ابتدا تصویر را از نوع unit8 به double تبدیل کنیم چون محدوده اعدادی که در uint8 جا می گیره بین 0 تا 255 هست که ما در محاسباتمون به عدد اعشار نیاز داریم مثل محاسبه mse یا فکر کنم یه جای دیگه

[zeinab 12]

سلام.اگر من بخوام برای ضرایب کوانتایز متفاوت mse و تصویر فشرده شده ام رو ببینم یک راهش اینه که برای تک تک ضرایب کوانتایز این کار رو انجام بدم که بهینه نیست و یکیشم اینه که یک functionتعریف کنم.
من میدونم چطور function را با یک ورودی و یک خروجی تنظیم کنم اما چطور یک ورودی (ضریب کوانتایز) بدم و دو خروجی (mseو تصویر کوانتایزشده)بخوام.
فرم تعریف functionرو ممنون میشم برام بگید.

[rahnema1]

سلام، به این صورت می تونید بنویسید البته می تونید از printf به بعد را از داخل تابع حذف کنید

function [mse compressed_image] = iframe(img, zarib)

	grayImg =double(rgb2gray(img));

	DC = blkproc(grayImg, [8 8], @dct2);

	Quantization_Matrix=...

		[ 16  11  10  16  24  40  51  61; 

		  12  12  14  19  26  58  60  55; 

		  14  13  16  24  40  57  69  56; 

		  14  17  22  29  51  87  80  62; 

		  18  22  37  56  68 109 103  77; 

		  24  35  55  64  81 104 113  92; 

		  49  64  78  87 103 121 120 101; 

		  72  92  95  98 112 100 103  99 ];

	if zarib < 50

		Quantization_Matrix = Quantization_Matrix * 50 / zarib ;

	else 

		Quantization_Matrix = Quantization_Matrix * (100 - zarib) / 50;

	end

	Quantization_Matrix = double(uint8((Quantization_Matrix - .2) / 254.8 * 254 + 1));

	quantized_DC = round(blkproc(DC, [8 8], @(x, y) x ./ y , Quantization_Matrix));

	reconstructed_DC = blkproc(quantized_DC, [8 8], @(x, y) x .* y , Quantization_Matrix);

	compressed_image = round(blkproc(reconstructed_DC, [8 8], @idct2));

	mse =sum(( grayImg(:) - compressed_image(:)) .^ 2) / numel(grayImg);

	printf("MSE = %f\n", mse)

	figure;

	subplot(2,2,1),imshow(grayImg,[0,255]);

	subplot(2,2,2),imshow(compressed_image,[0,255]);

	subplot(2,2,3),imhist(grayImg);

	subplot(2,2,4),imhist(compressed_image);

end

این هم نحوه استفاده

img = imread ('cameraman.tif');

[mse compressed_image] = iframe(img, 1);

[zeinab 12]

یعنی اون mse ای که در کنار compressed_image آوردین هم mse رو میده و هم تصویر compressed_imageرو؟!!!
نباید برای برگرداندن این دومقدار دستور imshow رو بیاریم؟

[rahnema1]

هر کدم از mse و یا compressed_image را جداگانه می تونیم استفاده کنیم
مثلا mse را چاپ کنیم یا یا compressed_image را imshow کنیم

[zeinab 12]

این اروری که در زمان اجرای برنامه متلب بهم میده چیه؟

??? Error using ==> rgb2gray>parse_inputs at 82MAP must be a m x 3 array.





Error in ==> rgb2gray at 35

X = parse_inputs(varargin{:});

[rahnema1]

rgb2gray را بردارید . احتمالا وقتی تصویر سیاه و سفید باشه ایراد می گیره