%function?aftans?=?image_rotate（imagepath?theta）
%?兩維圖像旋轉(zhuǎn)以及雙線性灰度差值算法的實現(xiàn)
%?Ref.?章毓晉.?圖像工程（上冊）——圖像處理.?清華大學(xué)出版社
%?Author:?lskyp????Date:?2009.08.12
im_init?=?imread（‘desert.jpg‘）;
%im_init?=?imread（imagepath）;
im_init?=?double（im_init）;?%對于灰度圖像，uint8表示范圍[0，255]，double型表示范圍[0，1]
im_height?=?size（im_init1）;
im_width?=?size（im_init2）;
%?分別處理灰度圖像和RGB圖像
if?ndims（im_init）?==?3
????im_final?=?zeros（im_heightim_width3）;
????R?=?im_init（::1）;
????G?=?im_init（::2）;
????B?=?im_init（::3）;
????R_final?=?im_final（::1）;
????G_final?=?im_final（::2）;
????B_final?=?im_final（::3）;
else
????im_final?=?zeros（im_heightim_width）;
end
?theta?=?pi/5;
rot_matrix?=?[cos（theta）?-sin（theta）;sin（theta）?cos（theta）];
tranh?=?（im_height+1）/2;
tranw?=?（im_width+1）/2;
for?h?=?1:im_height
????for?w?=?1:im_width
????????%?平移至原點，旋轉(zhuǎn)，然后再平移回去
????????new_position?=?rot_matrix*[h?-?tranh;w?-?tranw]+?[tranh;?tranw];
????????%?這里是認為圖像中心處為原點，建立局部坐標(biāo)系，方向還是向下向右
????????%?旋轉(zhuǎn)后，再平移回去
????????new_position（1）?=?mod（new_position（1）im_height）;
????????new_position（2）?=?mod（new_position（2）im_width）;
????????%?由于旋轉(zhuǎn)有可能會跑出圖像邊界，進行了模操作
????????%?這是控制操作范圍的好方法
????????if?new_position（1）?==?0
????????????new_position（1）?=?1;
????????end
????????if?new_position（2）?==?0
????????????new_position（2）?=?1;
????????end
????????%?當(dāng)然，圖像陣默認是從1開始計量，所以不能出現(xiàn)0
????????%?如果新位置為整數(shù)，那么直接賦予灰度值或者RGB值
????????%?否則，按照雙線性插值計算，使用后向映射
????????if?new_position?==?round（new_position）
????????????if?ndims（im_init）?==?3
????????????????R_final（hw）?=?R（new_position（1）new_position（2））;
????????????????G_final（hw）?=?G（new_position（1）new_position（2））;
????????????????B_final（hw）?=?B（new_position（1）new_position（2））;
????????????else
????????????????im_final（hw）?=?im_init（new_position（1）new_position（2））;
????????????end
????????????%?這里可以出對應(yīng)的final位置的像素是旋轉(zhuǎn)后所在位置對應(yīng)的原圖像像素值
????????else
????????????h_new?=?floor（new_position（1））;
????????????w_new?=?floor（new_position（2））;
????????????%?先取整，向0方向取整
????????????if?h_new?==?0
????????????????h_new?=?1;
????????????end
????????????if?w_new?==?0
????????????????w_new?=?1;
????????????end
????????????%?消除0位置
????????????if?ndims（im_init）?==?3
????????????????%?雙線性插值的實現(xiàn)過程?Ref.?章毓晉.?圖像工程（上冊）
????????????????%?——圖像處理.?清華大學(xué)出版社
????????????????R_temp1?=?R（h_new?+?1w_new）*（new_position（1）?-?h_new）+?R（h_neww_new）*（h_new?+?1?-?new_position（1））;
????????????????R_temp2?=?R（h_new?+?1w_new?+?1）*（new_position（1）?-?h_new）...
????????????????????????+?R（h_neww_new?+?1）*（h_new?+?1?-?new_position（1））;
????????????????%?先進行h方向的線性插值
????????????????R_final（hw）?=?R_temp1*（w_new?+?1?-?new_position（2））...
?????????????????????????????+?R_temp2*（new_position（2）?-?w_new）;
????????????????%?再進行w方向的線性插值
????????????????G_temp1?=?G（h_new?+?1w_new）*（new_position（1）?-?h_new）...
????????????????????????+?G（h_neww_new）*（h_new?+?1?-?new_position（1））;
????????????????G_temp2?=?G（h_new?+?1w_new?+?1）*（new_position（1）?-?h_new）...
???

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????目錄??????????0??2010-03-23?23:43??圖像縮放

?????文件???????5195??2009-11-28?12:38??圖像縮放\11111.asv

?????文件???????5198??2009-11-28?13:15??圖像縮放\11111.m

?????文件??????15026??2009-11-28?09:30??圖像縮放\cup.jpg

?????文件?????845941??2009-04-22?16:55??圖像縮放\Desert.jpg

?????文件???????2278??2009-11-28?09:29??圖像縮放\images.jpg

?????文件???????1616??2009-11-28?19:22??圖像縮放\im_f.asv

?????文件???????1613??2009-11-28?19:27??圖像縮放\im_f.m

?????文件????????739??2009-11-28?18:01??圖像縮放\im_K.asv

?????文件????????739??2009-11-28?18:01??圖像縮放\im_K.m

?????文件???????1074??2009-11-28?12:46??圖像縮放\im_R.m

?????文件???????1392??2009-11-28?20:30??圖像縮放\im_rotate.asv

?????文件???????1428??2009-11-28?20:44??圖像縮放\im_rotate.m

?????文件???????3224??2009-11-28?19:26??圖像縮放\man.jpg

-----------?---------??----------?-----??----

???????????????885463????????????????????14