To use PDL for this, you'd use qsorti on a slice to get the relevant column, then dice_axis the original structure on the sorted indices.