实验一
一、插值和采样
1 （a）读入图像head.jpg并显示。

>> A=imread('C:\Documents and Settings\Administrator\桌面\head.jpg'); imshow(A)
（b）计算图像维度。

>> size(A)
ans =
256 256
（c）此图像大小为40cm*40cm，计算图像的采样距离。

x=40cm/256=0.15625cm=1.5625mm
同理，y=1.5625mm
（d）逻辑坐标（图像坐标）为（22, 54）、（126, 241）的点，其空间坐标是多少？
逻辑坐标（图像坐标）为（22, 54）的点，空间坐标为
(22*1.5625mm,54*1.5625mm)=(34.375mm,84.375mm)；
同理，逻辑坐标（图像坐标）为（126, 241）的点，空间坐标为(196.875mm,376.05625mm)。

（e）求空间坐标为（14.2188, 5.3125）、（21.4063,34.5313）处的像素值。

>>b=40/256;
b =
0.1563
>> x3=(14.2188/b)-1
x3 =
90.0003
>> y3=(5.3125/b)-1
y3 =
33
>> x4=floor(x3)
x4 =
90
>> x5=ceil(x3)
x5 =
91
>> C1=A([x4],[y3])
C1 =
115
>> C2=A([x5],[y3])
C2 =
108
>> C3=0.9997*C1+0.0003*C2
C3 =
115
>> x6=(21.4063/b)-1
x6 =
136.0003
>> y6=(34.5313/b)-1 y6 =
220.0003
>> x7=floor(x6)
x7 =
136
>> x8=ceil(x6)
x8 =
137
>> y7=floor(y6)
y7 =
220
>> y8=ceil(y6)
y8 =
221
>> C4=A([x7],[y7]) C4 =
128
>> C4=A([x7],[y8]) C4 =
98
>> C4=A([x7],[y7]) C4 =
128
>> C5=A([x7],[y8]) C5 =
98
>> C6=A([x8],[y7])
C6 =
121
>> C7=A([x8],[y8])
C7 =
127
>> C8=0.9997*C4+0.0003*C6
C8 =
128
>> C9=0.9997*C5+0.0003*C7
C9 =
98
>> C10=0.9997*C8+0.0003*C9
C10 =
128
综上所述，空间坐标为（14.2188, 5.3125）处的像素值为115，（21.4063,34.5313）处的像素值为128。

（f）使用最邻近插值法，求在逻辑坐标（65.2, 193.7）和空间坐标（22.58, 7.24）处的像素值。

>> C11=A([65],[194])
C11 =
91
>> x9=(22.58/b)-1
x9 =
143.5120
>> y9=(7.24/b)-1
y9 =
45.3360
>> x10=ceil(x9)
x10 =
144
>> y10=floor(y9)
y10 =
45
>> C12=A([x10],[y10])
C12 =
94
综上所述，逻辑坐标为（65.2, 193.7）的点像素值为91，空间坐标（22.58, 7.24）的点像素值为94。

（g）使用双线性插值法，求逻辑坐标为（1.5,1.5）处的像素值。

>> C13=A([1],[1])
C13 =
4
>> C14=A([1],[2])
C14 =
3
>> C15=A([2],[1])
C15 =
3
>> C16=A([2],[2])
C16 =
2
>> C17=0.5*C13+0.5*C14
C17 =
4
>> C18=0.5*C15+0.5*C16
C18 =
3
>> C19=0.5*C17+0.5*C18
C19 =
4
综上所述，逻辑坐标为（1.5,1.5）的点像素值为4。