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

Home » Public Forums » archive » Reconstruct surface from gradient field?
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Reconstruct surface from gradient field? [message #87532] Thu, 13 February 2014 14:24 Go to previous message
dg86 is currently offline  dg86
Messages: 118
Registered: September 2012
Senior Member
Dear Folks,

I have measurements of the gradient field (tangent vectors) associated with
a two-dimensional surface, and would like to compute the height of the surface
itself. I understand that there are standard algorithms to do this, and wonder if
any of them have been implemented in IDL. If so, I'd be grateful for a pointer.

More specifically, I have values grad_x(x,y) and grad_y(x,y) that represent the
gradient of an unknown function z(x,y) measured at irregularly gridded locations
x and y in the plane. I can use griddata() to interpolate the gradient
functions smoothly across the plane. I want to calculate z(x,y) from grad_x and grad_y.
Moreover, I want to do this as accurately as possible.

Simple line integration is OK, but yields errors in the estimate for z(x,y) that vary
widely across the field of view. I've also taken a first shot at using FFTs to solve
Poisson's equation for z(x,y). Tests with simulated data were not encouraging.
These disappointing first efforts make me very hopeful indeed that somebody else
has a good solution already coded up, and is willing to share!

All the best,

David
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Escape "/" for NG LEGEND
Next Topic: MPFIT2DFUN- use for a time series of images?

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

Current Time: Wed Oct 08 15:06:34 PDT 2025

Total time taken to generate the page: 0.00387 seconds