Scatteredinterpolant (linear) from symmetric data does not produce symmetric isolines

I utilised the scatteredinterpolant function to generate an interpolation from which I derived isolines. My data are situated on a non-equispaced grid of points that are symmetric with respect to the midline. In fact, I only have data on one side, so I placed data on the other side by copying the values from the original data set, exploiting the symmetry.

When I generate iso-lines using the NATURAL method, the isolines have a symmetric trend. If I use the simplest method, i.e. LINEAR, the iso-lines are not symmetric, which does not make sense.

Do you know what may cause this issue and how to sort it out?

Answer 1

[Bruno Luong]

scatteredinterpolant all methods is based on triangulation of the data points. Even the data points are symmetric, the triangulation migh not.

[ES_Thorny]

Thanks for you answer. Is it possible to force a symmetric triangulation in some way?

[John D'Errico]

@ES_Thorny

No. You cannot enforce a symmetric triangulation in scatteredInterpolant. That would require you to be able to supply the triangulation itself to scatteredInterpolant, and this is not an option.

Could you create a triangulation that has the property of symmetry that you desire, and using other interpolation tools perform the interpolation? The answer is, well, yes, in theory. It would require some serious effort on your part, and it would force you to learn a lot about triangulation methods. Is it worth the effort? That is, just to make a subtly different picture that you prefer, when as you found, other methods were better able to handle your problem?

[Bruno Luong]

The easiest workaround is to symmetrize the interpolation result if you query points are symmetric

[ES_Thorny]

Yes, I thought about it, but how can it be done? I need to understand what is the output of scatteredinterpolant ... and how I can ducplicate it on the other side ...

[Walter Roberson]

Map the query location...

QX = X;
mask = QX > MidPoint;
QX(mask) = 2 * MidPoint - QX(mask);

[Bruno Luong]

@ES_Thorny post your code (that can run) if you want to get any specific help. The output is scattered is the z coordinates of the query points (xq,yq) assumed you want to interpolated 2D function. It's described in the doc page. If you have specific question just ask.