in reply to Working with PDL, need index range of consecutive values

Here's some building blocks (I'm too lazy for a full answer):

# I assume you are using 2 different pdls, not rows of the same pdl as + you mentioned in your post # If it's just one pdl, 'dog' it my $value = pdl [ 1, 20, 60, 3, 900, 34, 93, 5, 12, 24, 16, 200 ]; my $mask = pdl [ 1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 1, 1 ]; # get the indices of $mask that are 0 my $indices = which $mask == 0; # detect lower and upper edges of sequential regions of indices my $lower_edges = $indices - 1 != $indices->rotate(1); my $upper_edges = $lower_edges->rotate(-1); # also: # my $upper_edges = $indices + 1 != $indices->rotate(-1) # now slice and dice up your piddles knowing the regions you want.

Alternative: detect edges of regions of equal value in the mask and then only select even (or odd, depending on the first value) numbered regions. I chose to filter first to eliminate data. You would avoid indirection through the $indices pdl latter by detecting region edges in the mask in stead of detecting region edges in the indices of the mask that are 0. Try both, keep the better (for some value of better) one.

I feel there should be a more elegant way to detect edges of regions.

Update: There is a more elegant way:

use PDL::ImageND; my $edge_detector_kern = pdl [ -1, 1 ]; my $upper_edges = convolveND( $mask, $edge_detector_kern ); # where it is -1 is the upper edge of a region of zeroes # where it is 1 is the upper edge of a region of ones

Replies are listed 'Best First'.
Re^2: Working with PDL, need index range of consecutive values
by Anonymous Monk on Feb 21, 2009 at 23:58 UTC
    Thank you for your help in the matter. I am going to implement this soon and will keep you posted on the results! Thanks again.
Re^2: Working with PDL, need index range of consecutive values
by Anonymous Monk on Feb 27, 2009 at 22:43 UTC
    Update: I have tried your more elegant solution and this works perfectly for what I needed! Thank you so much for your help, and your explanation of things. I am very new to PDL (I'm a graduate student in Texas), and my bosses are familiar with PDL but not quite as savvy as you folk. Piddles take some getting used to, but I think I'm starting to get the hang of them.