why do I get a strange result with griddata cubic?

Question

Valeria Leto 2021-7-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/886424-why-do-i-get-a-strange-result-with-griddata-cubic

评论： Valeria Leto 2021-7-30

Hi! I used the griddata cubic interpolation for the values of an image. I get this image that has a much higher range, from -300 to 400 aproximately. The non processed image has a range from 90 to 130. Any ideas why this happens?Thanks!

3 个评论
显示 1更早的评论隐藏 1更早的评论

Valeria Leto 2021-7-27

在 MATLAB Online 中打开

u=[-20:0.05 :20];%5 cm

v=[0:-0.05:-20];

[Xq,Yq] = meshgrid(u,v);

Zq=zeros(length(v),length(u));

figure(30)

mesh(Xq,Yq,Zq)

view(2)

% Interpolazione con griddata: CUBIC
vq_4 = griddata(x,y,C,xq,yq,'cubic');
Unrecognized function or variable 'x'.
V_4=zeros(length(v),length(u));
for i=1:1:length(v)
V_4(i,:)=vq_4(1+(i-1)*801:1+(i-1)*801+800,1)';
end
%% PLOT
figure(34)
surf(Xq,Yq,Zq,V_4,'EdgeColor','none')
title('INTERPOLAZIONE CUBIC RISOLUZIONE 5 cm')
colormap gray
colorbar
axis equal
xlabel('est')
ylabel('nord')
% view(2)

Matt J 2021-7-27

在 MATLAB Online 中打开

The code you've posted does not run, I'm afraid:

Unrecognized function or variable 'x'.
Error in test (line 12)
vq_4 = griddata(x,y,C,xq,yq,'cubic');

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Walter Roberson 2021-7-27

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/886424-why-do-i-get-a-strange-result-with-griddata-cubic#answer_754904

在 MATLAB Online 中打开

This is expected behavior when there are sharp edges, especially with a surrounding flat area.

     B--
     |
   -A-

Especially if the region near B is dense, but the edge itself is sparse, then the polynomial generated near the edge needs to mathematically dip down below the base in order to gain the necessary steepness to match the flat control points to the left together with the high points to the right. Polynomials cannot just suddenly rise steeply without a leading edge, so the math needs to insert a downward edge to be able to rise from.

9 个评论
显示 7更早的评论隐藏 7更早的评论

Walter Roberson 2021-7-28

在 MATLAB Online 中打开

Your image is grayscale, so each location is a grayscale intensity.

Take a horizontal cross-section across your image. Moving left to right, you have dull gray, brighter color, narrow black, a bit brighter, then a bit duller.

Now imagine plotting that cross-section line with position on the x axes, and intensity on the y axes. Low, high, sudden very low, up to moderate again soon after, then medium.

         ----
 -------/    \--+ +-----
                | |
                 v
                 

So on the cross-section graph of intensity, there are edges.

And what you are finding is that near that sharp edge where it goes black, that you are getting interpolated values that are negative. This is not unexpected.

You are using cubic fit, so at each step there are four control values that a degree 3 polynomial is being fit to. Imagine you have some points

D

A B C

In order to fit a cubic to those four points, there has to be two reversal of direction, and the line has to exit on opposite sides of the baseline (if it starts from above it has to end below). In this case with D being above C, the most natural fit would be that the line goes up through C into D and exits above right, which would require that the line enters below left, going up through A, reversing direction between A and B, going down through B, reversing direction between B and C, and going up through C into D. Notice that this implied a minima between B and C -- a point below the baseline.

Now as you bring B and C closer together, the dip and rise has to get steeper and steeper. If your control points are close enough, you can end up with an arbitrarily low minima between B and C that is mathematically required to fit a cubic and yet is not present in any of the data.

/ \ /D

A B C

/ \ /

Likewise as you bring A and B together, you can end up with an implied peak that is higher than any of the datapoints you have.

Bruno Luong 2021-7-30

在 MATLAB Online 中打开

working with xs/ys instead of x/y

xs = sx*x
ys = sy*y

sx, sy are appropriate scaling factors that you have to figure it out.

Valeria Leto 2021-7-30

@Bruno Luong thank you

请先登录，再进行评论。

why do I get a strange result with griddata cubic?

3 个评论
显示 1更早的评论隐藏 1更早的评论

采纳的回答

9 个评论
显示 7更早的评论隐藏 7更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

Community Treasure Hunt

why do I get a strange result with griddata cubic?

3 个评论 显示 1更早的评论隐藏 1更早的评论

采纳的回答

9 个评论 显示 7更早的评论隐藏 7更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

Community Treasure Hunt

3 个评论
显示 1更早的评论隐藏 1更早的评论

9 个评论
显示 7更早的评论隐藏 7更早的评论