The Math Forum

Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Math Forum » Discussions » Software » comp.soft-sys.matlab

Topic: Bivariate Pchip
Replies: 9   Last Post: Jan 19, 2013 6:14 AM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
John D'Errico

Posts: 9,128
Registered: 12/7/04
Re: Bivariate Pchip
Posted: Jan 18, 2013 4:46 PM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

"Bruno Luong" <b.luong@fogale.findmycountry> wrote in message <kdccnj$kcc$>...
> "John D'Errico" <> wrote in message <kdbk43$afl$>...

> >
> > pchip in 2-d as a tensor product form has been shown
> > NOT to be adequate for the general desired behavior.
> > (It sometimes will produce an acceptable result, but in
> > general, it is not adequate.)
> >

> John, I must admit that I don't know what is the adequate definition of shape preserving interpolation in 2D.
> I just though at least the tensorial product would preserve at least the monotonic in any cut parallel to the two axis. Unless I'm mistaken, this also implies that the interpolation within a patch will necessary be bounded by the four corner data, thus never overshoot.
> May be these characteristics are not enough for OP, but it is better than nothing. No?
> Bruno

Yes, monotonicity is not a trivial thing to discuss
when you move to more than one independent

It has been many years since I looked at it, but I
recall the statement that a simple tensor product
version of pchip is not adequate here. And I don't
recall under which circumstances that interpolant

The case of a linear tensor product interpolant
(often known as a bilinear interplant, as used by
photoshop) is a good one to study, as that may
offer an idea. Consider the function f(x,y), where
it is defined at the 4 corners of the unit square,
and we will use bilinear interpolation over that

f((0,0) = 1
f(0,1) = 0
f(1,0) = 0
f(1,1) = 1

Within the unit square, the tensor product linear
interpolant reduces to

f(x,y) = 1 - y - x + 2*x*y

which is clearly not linear. If we hold either x or y
fixed, then of course it is again linear. But if we
interpolate along some other path, perhaps a
diagonal one from corner to corner, then the
interpolant will be a quadratic polynomial with
min or max at the center of the square. This is
the classic problem with a tensor product
interpolant when applied in three dimensions
for color science problems, and it is why that
method is avoided.

What I don't know (without some study) is if the
tensor product pchip can still have non-monotonic
behavior parallel to an axis. (Perhaps tomorrow
if I have some energy, as I think I know how to
show if this can happen.)


Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© The Math Forum at NCTM 1994-2017. All Rights Reserved.