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

Home » Public Forums » archive » Re: Vertical profiles
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: Vertical profiles [message #1822 is a reply to message #1819] Mon, 07 March 1994 20:35 Go to previous messageGo to previous message
dogru is currently offline  dogru
Messages: 1
Registered: March 1994
Junior Member
In <19574.9403041804@sgsscsc2> sgsmatpa@reading.ac.uk (sgsmatpa) writes:


> Hello!
> Another plea for help here.

> I have some oceangraphic data that I would like to display
> using Pvwave but I'm not having much success. The data basically
> consists of a series of almost parallel, vertical slices
> through the ocean :-

> ... [picture deleted.]

> Data is in the form of x, y, z, temperature. I want to be able to display
> the 3D data and roatate it around the z axis to look from different angles
> and to contour or shade the surface with the temperature values. Later
> on I want to interpolate the data to give a full 3D volume but thats a
> different problem!

> SHADE_SURF doesnt seem to be able to cope with this kind of 'vertical'
> surface. I can cheat and map longditude to z and depth to y but this
> means I can't rotate the data around the depth axis.

> Any suggestions as to alternative ways of displaying it? Can I convert it
> to a volume and make most of the volume transparent except where I have
> data points? And if I do that can I still display axis etc?

I faced this same problem a short while ago. The trick lies in
using the IDL routines TRIANGULATE and TRIGRID. The first routine
converts the given points to a set of triangles in the space.
The second routine then takes this set and maps them to a volume.
If, for example, x, y, and z are each vectors of size N, you can
obtain the (equivalent) 3xN set of triangles by

TRIANGULATE, x, y, triangles

The variable 'triangles' will hold the set of triangles, which you
can use in TRIGRID to obtain a volume (or, a 2D array accepted by
SURFACE and SHADE_SURF) of 51x51, which is the default:

vol=TRIGRID(x, y, z, triangles)
SHADE_SURF, vol

> Also is there a 'null' value for data points (e.g. -9999.000) ?

Yes, one of the above routines accepts a keyword that lets you
change the default for the missing values. Usually, zero works
much better, however, if you later wish to plot the graph.

A couple minor annoyances: Even replacing the null values with
zeros resulted in graphs that the users did not want to see.
And second, after the above transformations, you lose the actual
axis ranges of your data. Mail me or post if you want more info.


Sait Dogru
dogru@cs.umn.edu


> Help! All advice very welcome.

> thanks
> Tricia Matthews
> University of Reading
> England
> sgsmatpa@uk.ac.reading
[Message index]
 
Read Message
Read Message
Read Message
Previous Topic: Empirical Orthogonal Function Analysis?
Next Topic: binary data files with ascii header

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

Current Time: Fri Oct 10 19:34:06 PDT 2025

Total time taken to generate the page: 0.23971 seconds