comp.lang.idl-pvwave archive: archive » interpolate weirdness

Home » Public Forums » archive » interpolate weirdness

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: interpolate weirdness [message #87660 is a reply to message #87659]

Fri, 21 February 2014 04:33

Helder Marchetto
Messages: 520
Registered: November 2011

Senior Member

On Friday, February 21, 2014 1:00:36 PM UTC+1, CR wrote:
> Hi Folks,
>
>
>
> I know that there is no 'correct' way for interpolation, but maybe a 'best' IDL way. The following example easily demonstrates what I mean:
>
>
>
> IDL> d=dist(20) & print,interpolate(d,[.5],[.5],/grid)
>
> 0.853553
>
>
>
> IDL> print,interpolate(d,[.5],[.5],/grid,/cubic)
>
> 0.622236
>
>
>
> IDL> print,interpolate(d,[.400544],[.400544],/grid)
>
> 0.707107
>
>
>
> IDL> print,interpolate(d,[.581946],[.581946],/grid,/cubic)
>
> 0.707107
>
>
>
> IDL> print,'The result should be: ',sqrt(.5)
>
> The result should be: 0.707107
>
>
>
>
>
> It does not play a role which IDL version is used - for IDL 6.4 and for IDL 8.3 I got the same results. Is there any way to get 'closer' to a reasonable result?
>
>
>
> Thanks and Cheers
>
>
>
> Chris

I'm not sure what you're trying to do here, but interpolating near borders is kind of complicated because you can only have one side to rely on. However, if you're going to use bicubic interpolation, I would recommend using:
print,interpolate(d,[.5],[.5],/grid,cubic=-0.5)
then you get:
0.738791
Hope it helps,
h

Report message to a moderator

[Message index]

		interpolate weirdness By: rogass on Fri, 21 February 2014 04:00
		Re: interpolate weirdness By: Helder Marchetto on Fri, 21 February 2014 04:33
		Re: interpolate weirdness By: Helder Marchetto on Fri, 21 February 2014 05:59
		Re: interpolate weirdness By: Helder Marchetto on Fri, 21 February 2014 06:06
		Re: interpolate weirdness By: rogass on Fri, 28 February 2014 08:31

Previous Topic:	Behaviour of colorbar title in z-buffer?
Next Topic:	Orientation of Contour Labels

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sat Oct 11 13:23:48 PDT 2025

Total time taken to generate the page: 0.95943 seconds