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 
Switch to threaded view of this topic Create a new topic Submit Reply
Re: Vertical profiles [message #1819] Mon, 07 March 1994 23:53
roe is currently offline  roe
Messages: 4
Registered: March 1994
Junior Member
Hello!

I'm analysing a couple of SAR images from the ERS-1 satellite and am wondering if anybody has IDL routines for speckle filtering.

Thanks a lot.

---
| HansPeter Roesli Internet: roe@otl.sma.ch |
| Swiss Meteorological Institute Telex: 846 007 |
| Osservatorio Ticinese Fax: +4193 326 310 |
| CH-6605 Locarno-Monti, Switzerland Phone: +4193 326 319 |
Re: Vertical profiles [message #1822 is a reply to message #1819] Mon, 07 March 1994 20:35 Go 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
Re: Vertical profiles [message #1823 is a reply to message #1822] Mon, 07 March 1994 00:22 Go to previous message
knipp is currently offline  knipp
Messages: 68
Registered: January 1993
Member
In article 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 :-
>
> ____ ____ ___ y (long)
> /| /| /| /| /| | /
> / | / | / | / | / | | /
> / |/ |/ |/ |/ | | /
> / /___/ /___/ | | --------------> x(lat)
> | | | | | | | |
> | | | | | | | |
> | | | | | |__| |
> | | | | | / |
> | /| | /| | / |
> | / | | / | | / z(depth)
> |/ |___|/ |___|/
>
...
>
> thanks
> Tricia Matthews
> University of Reading
> England
> sgsmatpa@uk.ac.reading


Hi

I guess there's a routine existing in PvWave called PLOT3D
(..../lib/user/plot3d.pro).

It may be the solution to your problem.


Karl


____________________________________________________________ __________________
__ ____ __
/ // _ \ / / Karlheinz Knipp phone: +49 511 - 762 4922
/ // /_/ // / University of Hannover fax: +49 511 - 762 2483
/ // ____// / Institute for Photogrammetry
/ // / / / Nienburger Str.1
/_//_/ /_/ FRG 30167 Hannover 1 email: knipp@ipi.uni-hannover.de
  Switch to threaded view of this topic Create a new topic Submit Reply
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: Wed Oct 08 19:51:32 PDT 2025

Total time taken to generate the page: 0.00679 seconds