comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Re: Interpolation on a sphere
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Switch to threaded view of this topic Create a new topic Submit Reply
Re: Interpolation on a sphere [message #2925] Sun, 09 October 1994 05:53
cooper is currently offline  cooper
Messages: 1
Registered: October 1994
Junior Member
> In article <373pgv$gv@danberg.llnl.gov>, dan@danberg.llnl.gov (Dan Bergmann) writes:
> |>
> |> Using [simple interpolation], a data value at longitude=0 latitude=89
> |> will have very little influence on the interpolated value at
> |> longitude=180 latitude=89,
> |> but in reality, these points are right next to each other.


How about copying data near longitudes 0 and 360 (or 180 E and
180 W, of course) to "ghost" zones _above_ 360 and _below_ 0?
That is, map longitude = 0-2 (or whatever) onto longitude 360-362,
and 358-360 onto (-2)-0. Then interpolate as usual; then strip
off the ghost zones before using the result.

It isn't a pretty solution, but I thought it might work...

Glenn Cooper
Re: Interpolation on a sphere [message #2929 is a reply to message #2925] Fri, 07 October 1994 16:54 Go to previous message
todd is currently offline  todd
Messages: 25
Registered: February 1994
Junior Member
In article <373pgv$gv@danberg.llnl.gov>, dan@danberg.llnl.gov (Dan Bergmann) writes:
|>
|> In article <372kp6$54d@news.mic.ucla.edu>, todd@artemis.ess.ucla.edu (Todd Ratcliff) writes:
|> |> Dan,
|> |>
|> |> If you just want to see your data wrapped onto a sphere you can
|> |> simply interpolate your data into a 180X360 array using whatever
|> |> interpolation that you normally use, generate a
|>
|> Using this method, a data value at longitude=0 latitude=89 will have
|> very little influence on the interpolated value at longitude=180 latitude=89,
|> but in reality, these points are right next to each other. I'm not so interested
|> in viewing my data as I am in doing a good interpolation on the surface of
|> a sphere.


Ah, good point! In that case if you find a good spherical interpolator, please
let me know.


Todd
Re: Interpolation on a sphere [message #2932 is a reply to message #2929] Fri, 07 October 1994 08:29 Go to previous message
dan is currently offline  dan
Messages: 27
Registered: March 1993
Junior Member
In article <372kp6$54d@news.mic.ucla.edu>, todd@artemis.ess.ucla.edu (Todd Ratcliff) writes:
|> Dan,
|>
|> If you just want to see your data wrapped onto a sphere you can
|> simply interpolate your data into a 180X360 array using whatever
|> interpolation that you normally use, generate a

Using this method, a data value at longitude=0 latitude=89 will have
very little influence on the interpolated value at longitude=180 latitude=89,
but in reality, these points are right next to each other. I'm not so interested
in viewing my data as I am in doing a good interpolation on the surface of
a sphere.

|> spherical polygon (PV~Wave has Poly_Sphere and I've seen an IDL
|> routine called Mesh_Obj by Dan Carr which can do the same thing)
|> with 180 polygons around equator and 360 polygons around the
|> meridian then use polyshade with the Shades=Bytscl(your_array)
|> keyword.
|>
|> You'll have your data on sphere.
|>
|>
|> Todd
--
************************************************************ ***
** Dan Bergmann dbergmann@llnl.gov **
** Global Climate Research fax (510) 422-5844 **
** Lawrence Livermore National Lab human (510) 423-6765 **
Re: Interpolation on a sphere [message #2941 is a reply to message #2932] Thu, 06 October 1994 22:02 Go to previous message
todd is currently offline  todd
Messages: 25
Registered: February 1994
Junior Member
Dan,

If you just want to see your data wrapped onto a sphere you can
simply interpolate your data into a 180X360 array using whatever
interpolation that you normally use, generate a
spherical polygon (PV~Wave has Poly_Sphere and I've seen an IDL
routine called Mesh_Obj by Dan Carr which can do the same thing)
with 180 polygons around equator and 360 polygons around the
meridian then use polyshade with the Shades=Bytscl(your_array)
keyword.

You'll have your data on sphere.


Todd
  Switch to threaded view of this topic Create a new topic Submit Reply
Previous Topic: CALL_EXTERNAL-->Fortran in Solaris
Next Topic: The Game of LIFE

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Wed Oct 08 19:13:41 PDT 2025

Total time taken to generate the page: 0.00534 seconds