| Interpolation: grid --> observation point [message #19103] |
Thu, 24 February 2000 00:00  |
Andy Loughe
Messages: 174
Registered: November 1995
|
Senior Member |
|
|
Hello,
I have inherited code that takes atmospheric model (gridded) data
and interpolates it to station observation points that report
hourly precipiation. The model points form a cartesian coordinate
system, while the station observations are spread randomly
within the boundaries of the full model domain.
Currently, a bi-linear approach is used to interpolate the
model data to each observation point. This method uses the
four closest model grid points which surround the observation
point.
I would like to use more than just the four model grid points
which surround the observation point, and weight the more distant
points appropriately.
Any suggestions on effective methods for doing this in IDL?
Thanks!
--------
Andrew F. Loughe email:loughe@fsl.noaa.gov phone:(303)497-6211
|
|
|
|
| Re: Interpolation: grid --> observation point [message #19174 is a reply to message #19103] |
Fri, 25 February 2000 00:00  |
wmc
Messages: 117
Registered: February 1995
|
Senior Member |
|
|
Andy Loughe <loughe@fsl.noaa.gov> wrote:
> I have inherited code that takes atmospheric model (gridded) data
> and interpolates it to station observation points that report
> hourly precipiation. The model points form a cartesian coordinate
> system, while the station observations are spread randomly
> within the boundaries of the full model domain.
> Currently, a bi-linear approach is used to interpolate the
> model data to each observation point. This method uses the
> four closest model grid points which surround the observation
> point.
> I would like to use more than just the four model grid points
> which surround the observation point, and weight the more distant
> points appropriately.
Is this a good idea? If there is any justification for
doing this, its that the model numerics means that more than just the
surrounding grid points influences a point. But in that case, your
interpolation should be based on the model numerics.
This might just be sensible for, say, surface pressure, but I'd have
thought that it would not be appropriate for ppn (by virtue of how its
generated).
In fact, one could probably make a reasonable case for saying that, for
convective ppn, even bilinear interpolation is inappropriate and you
you should just use the value from whichever grid square the station
happens to be in.
-W
--
William M Connolley | wmc@bas.ac.uk | http://www.nbs.ac.uk/public/icd/wmc/
Climate Modeller, British Antarctic Survey | Disclaimer: I speak for myself
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
