Re: Vertical profiles [message #1819] |
Mon, 07 March 1994 23:53 |
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  |
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  |
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
|
|
|