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 #87863 is a reply to message #87662]

Fri, 28 February 2014 08:31

rogass
Messages: 200
Registered: April 2008

Senior Member

Dear Helder,
thank you for clarification. The situation has changed. Desired values are now far away from boundaries and I can use something like this:

function cr_interpol,dat,mat_old,mat_new,cubic=cubic

ndim = size(dat,/ndimensions)
dims = size(dat,/dimensions)
case (1b) of
ndim eq 1 : begin
newx = (interpol(ulindgen(dims[0]),(mat_old[0,*])[*],(mat_new[0,*]) [*]))[*]
result = interpolate(dat[*],newx,/grid,cubic=cubic)
end
ndim eq 2 : begin
newx = (interpol(ulindgen(dims[0]),(mat_old[0,*])[*],(mat_new[0,*]) [*]))[*]
newy = (interpol(ulindgen(dims[1]),(mat_old[1,*])[*],(mat_new[1,*]) [*]))[*]
result = interpolate(dat,newx,newy,/grid,cubic=cubic)
end
ndim eq 3 : begin
newx = (interpol(ulindgen(dims[0]),(mat_old[0,*])[*],(mat_new[0,*]) [*]))[*]
newy = (interpol(ulindgen(dims[1]),(mat_old[1,*])[*],(mat_new[1,*]) [*]))[*]
newz = (interpol(ulindgen(dims[2]),(mat_old[2,*])[*],(mat_new[2,*]) [*]))[*]
result = interpolate(dat,newx,newy,newz,/grid,cubic=cubic)
end
endcase
return,result
end

Cheers

Chris

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: Fri Oct 17 21:30:53 PDT 2025

Total time taken to generate the page: 1.35705 seconds