Line Integral Convolution

Hi,

First of all, matplotlib is a very nice piece of software - I use it
all the time for interactive data analysis, for posters, and for
published papers.

For visualizing low-resolution vector fields, quiver() is quite nice.
But for high-resolution vector fields there is a technique called
"line integral convolution" that can be quite effective. Matplotlib
doesn't seem to implement it, but I have written a simple line
integral convolution operator I'd be happy to contribute. It's written
in cython, and is quite fast and fairly general. There is no
high-level interface like quiver() yet, in part because of the many
ways you might want to use the operator. In fact I don't know for sure
that it belongs here rather than (say) scipy.

Are you interested in including this in matplotlib?

Anne

Anne Archibald wrote:

Hi,

First of all, matplotlib is a very nice piece of software - I use it
all the time for interactive data analysis, for posters, and for
published papers.

For visualizing low-resolution vector fields, quiver() is quite nice.
But for high-resolution vector fields there is a technique called
"line integral convolution" that can be quite effective. Matplotlib
doesn't seem to implement it, but I have written a simple line
integral convolution operator I'd be happy to contribute. It's written
in cython, and is quite fast and fairly general. There is no
high-level interface like quiver() yet, in part because of the many
ways you might want to use the operator. In fact I don't know for sure
that it belongs here rather than (say) scipy.

Are you interested in including this in matplotlib?

Possibly. As you note, the question may be whether this belongs in scipy, or in an mpl toolkit, or in mpl itself. I had to google to find out what "line integral convolution" is. It does look like it would be a nice thing to offer. I presume the actual plot would be via image (or maybe pcolorfast if anything but a uniform rectangular grid is used). Is that correct, or are there other ways the output of the operator would be displayed?

So far, mpl has not had pyrex or cython as a build dependency, but cython is used in the basemap toolkit. Personally, I like pyrex/cython very much and would be happy to see one of them--presumably cython--used
in mpl. Others may disagree.

Based on what you say above, my guess is that scipy might be the most natural home for your code, in which case we would want to have an example of its use for generating an image in the set of mpl examples. But if something more elaborate is needed in the way of a high-level interface, then either putting it in mpl directly, or in a toolkit, would be good.

Eric

···

Anne

-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
Build the coolest Linux based applications with Moblin SDK & win great prizes
Grand prize is a trip for two to an Open Source event anywhere in the world
http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Matplotlib-devel mailing list
Matplotlib-devel@lists.sourceforge.net
matplotlib-devel List Signup and Options

Anne-

I would be interested in seeing the code, regardless of where it finds a home.

Would you mind sharing, or would you rather wait?

-Rob

···

On Oct 13, 2008, at 4:59 PM, Anne Archibald wrote:

but I have written a simple line
integral convolution operator I'd be happy to contribute.

----
Rob Hetland, Associate Professor
Dept. of Oceanography, Texas A&M University
http://pong.tamu.edu/~rob
phone: 979-458-0096, fax: 979-845-6331

2008/10/16 Rob Hetland <hetland@...345...>:

but I have written a simple line
integral convolution operator I'd be happy to contribute.

Anne-

I would be interested in seeing the code, regardless of where it finds a
home.

Would you mind sharing, or would you rather wait?

Well, just now I'm absolutely swamped with other work, but I thought
it would be good to make the code I have available:
http://www.scipy.org/Cookbook/LineIntegralConvolution

Good luck with it!
Anne

···

On Oct 13, 2008, at 4:59 PM, Anne Archibald wrote: