Matlab绘图基础——图形绘制的插值 以及 图像大小的重采样

使用说明：图形绘制时的插值

interp1   %1-D data interpolation

interpft  %使用fft算法插值

    %将原数据x转换到频率域，再逆转换回来更密集的数据采样点

spline    %一维Cubic spline（三次样条） data interpolation

interp2   %2-D data interpolation

interp3   %3-D data interpolation (table lookup)

注：The interp2 command interpolates between data points. It finds values of a two-dimensional function f(x,y) underlying the data at intermediate points.

案例1：二维插值显示

[X,Y] =meshgrid(-3:.25:3);Z=peaks(X,Y);     %间隔为0.25

[XI,YI]=meshgrid(-3:.33:3,-3:0.05:3);          %XI间隔为0.33，YI间隔为0.05

ZI=interp2(X,Y,Z,XI,YI,'cubic');

                 %设定插值使用的算法。其他选项还有'nearest','linear,'spline'

mesh(X,Y,Z), hold on;axis([-3 3 -3 3 -5 20]);

mesh(XI,YI,ZI+15),hold off;   %绘制插值算得的图形

 

案例2：图像大小的重采样

M = imread('rice.png'); 

%% 使用插值的方法对图像进行重采样

% the original image

[p,q] = size(M);

figure;imshow(M);title('Original Image')

% original location

[X,Y] = meshgrid( (0:p-1)/(p-1), (0:q-1)/(q-1) );

 

% target image size

p1 = 100; q1 = 100;

% new sampling location

[XI,YI] = meshgrid( (0:p1-1)/(p1-1) , (0:q1-1)/(q1-1) );

% the new image 插值方法可以自己更改

M1 = interp2( X,Y,double(M), XI,YI ,'cubic');  % 注意这里的M原本是unit8类型的数据，这里必须转化为double类型的数据

figure;imshow(M1,[0,255]);title('Resized Image')  % 这里用imshow显示double类型的图像时要指定显示范围[0,255]

 

%% 简洁方法（IPT工具箱）

figure;imshow(M);title('Original Image')

% 重采样

p1 = 100; q1 = 100;

M1 = imresize(M,[p1,q1]);

figure;imshow(M1);title('Resized Image')

 

%案例3：已知点不是规则的的网格点，是无规则的TIN模型，如何插值获得规则网格点上数据

慢慢更新。。。