Re: Interpolation: grid --> observation point

["J.D. Smith" <jdsmith(at)astro.cornell.edu>]

Andy Loughe wrote:
> 
> 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

See the CUBIC keyword to either poly_2d or interpolate.	 Full n-point order
polynomial interpolation is not available, but might not be what you want
either. You might also look at krig2d.

JD

-- 
 J.D. Smith                             |*|      WORK: (607) 255-5842    
 Cornell University Dept. of Astronomy  |*|            (607) 255-6263
 304 Space Sciences Bldg.               |*|       FAX: (607) 255-5875 
 Ithaca, NY 14853                       |*|