3 篇文章 0 订阅

# 3.插值处理

     插值的任务是由已知的观测点为物理量，建立一个简单、连续的解析模型，以便能根据该模型推测物理量在非观测点处的特性。插值包括多项式插值、艾尔米特插值、分段插值和样条插值等等，其中最常见的有双线性插值和双三次插值。这里选用常见的双线性插值以及双三次插值。


## 3.1 双线性插值

	clear;
clc;

a = [27.50, 28.25, 29.75, 30.50, 28.50, 28.00, 28.25, 28.00
27.75, 27.25, 27.50, 29.00, 31.00, 29.25, 30.00, 29.25
26.50, 27.50, 26.75, 27.00, 30.75, 34.00, 46.50, 33.25
26.50, 26.75, 27.00, 27.00, 31.00, 39.00, 59.00, 40.75
26.00, 26.25, 25.75, 27.25, 29.50, 40.75, 59.50, 40.50
26.25, 25.50, 21.75, 23.75, 28.00, 41.00, 59.50, 37.50
26.50, 25.75, 17.50, 20.75, 28.00, 41.25, 53.25, 34.50
27.25, 25.75, 24.75, 25.00, 28.00, 33.75, 36.25, 33.00               %%存数据组1到a矩阵中
]
b = [28.00, 31.25, 33.00, 29.75, 31.75, 29.25, 27.25, 27.25
29.75, 34.25, 38.25, 35.25, 36.75, 31.75, 27.50, 26.50
30.00, 37.00, 39.25, 38.25, 37.75, 33.00, 28.25, 26.75
29.50, 34.50, 39.25, 38.75, 37.00, 31.75, 27.25, 26.75
28.00, 31.00, 33.75, 34.25, 32.25, 28.75, 26.25, 26.25
26.75, 27.25, 28.25, 28.50, 27.25, 26.00, 26.75, 26.25
26.25, 26.50, 26.25, 26.00, 26.25, 26.50, 26.00, 26.00
27.50, 26.00, 26.50, 25.75, 26.00, 26.50, 26.25, 27.25            %%存数据组2到b矩阵中
]

c = [];
d = [];                                  %构建插值后的矩阵c、d
i=0;
for hang=1:36                            %双线性插值
if hang==1
i=i+1;
end
if hang>2

if rem(hang,5)==2
i=i+1;

end
end

if rem(hang,5)==1
ku=0;
end

if rem(hang,5)==2
ku=0.2;
end

if rem(hang,5)==3
ku=0.4;
end

if rem(hang,5)==4
ku=0.8;
end

if rem(hang,5)==0
ku=1;
end
if hang==6 || hang==11 || ...
hang==16 || hang==21 || ...
hang==26 || hang==31 || ...
hang==36
ku=1;
end

j=0;
for lie=1:36
if lie==1
j=j+1;
end
if lie>2

if rem(lie,5)==2
j=j+1;
end
end
if rem(lie,5)==1
kv=0;
end
if rem(lie,5)==2
kv=0.2;
end
if rem(lie,5)==3
kv=0.4;
end
if rem(lie,5)==4
kv=0.6;
end
if rem(lie,5)==0
kv=0.8;
end
if lie==6 || lie==11 ||...
lie==16 || lie==21 ||...
lie==26 || lie==31 || ...
lie==36
kv=1;
end
u=ku;
v=kv;

c(hang,lie)=(1-u)*(1-v)*b(i,j)...       %%利用权重公式进行赋值
+(1-u)*v*b(i,j+1)...
+u*(1-v)*b(i+1,j)...
+u*v*b(i+1,j+1);

d(hang,lie)=(1-u)*(1-v)*a(i,j)...
+(1-u)*v*a(i,j+1)...
+u*(1-v)*a(i+1,j)...
+u*v*a(i+1,j+1);

end

end

subplot(2,2,1)                          %显示原有像素图以及插值后的像素图
imagesc(b);
subplot(2,2,2)
imagesc(c);
subplot(2,2,3)
imagesc(a);
subplot(2,2,4)
imagesc(d);


## 3.2 双三次插值

	clc;
clear;

b = [28.00, 31.25, 33.00, 29.75, 31.75, 29.25, 27.25, 27.25
29.75, 34.25, 38.25, 35.25, 36.75, 31.75, 27.50, 26.50
30.00, 37.00, 39.25, 38.25, 37.75, 33.00, 28.25, 26.75
29.50, 34.50, 39.25, 38.75, 37.00, 31.75, 27.25, 26.75
28.00, 31.00, 33.75, 34.25, 32.25, 28.75, 26.25, 26.25
26.75, 27.25, 28.25, 28.50, 27.25, 26.00, 26.75, 26.25
26.25, 26.50, 26.25, 26.00, 26.25, 26.50, 26.00, 26.00
27.50, 26.00, 26.50, 25.75, 26.00, 26.50, 26.25, 27.25
]

t=imresize(b,45,'bicubic');   %45倍插值

imagesc(t);


# 4.期望

• 6
点赞
• 31
收藏
• 打赏
• 3
评论
02-05
05-11 3232
10-15 2543
06-22
03-31 4700
01-12 3360
04-02
08-12 3940

### “相关推荐”对你有帮助么？

• 非常没帮助
• 没帮助
• 一般
• 有帮助
• 非常有帮助

©️2022 CSDN 皮肤主题：游动-白 设计师：我叫白小胖

¥2 ¥4 ¥6 ¥10 ¥20

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、C币套餐、付费专栏及课程。