Gaussian Smoothing with PDL

gistfile1.pl

#!/usr/bin/env perl

use strict;
use warnings;

use PDL;
use PDL::ImageND;
use PDL::Graphics::Prima::Simple;

my $x = sequence(50);
my $data = sin($x/4) + grandom(50)/20;

my $kernel = exp(-(rvals(5)/3)**2);
$kernel /= $kernel->sumover;

my $smoothed = $data->convolveND($kernel);

plot(
  -orig => ds::Pair( 
    $x, $data,
    colors => pdl(255,0,0)->rgb_to_color,
    plotType => ppair::Lines,
  ),
  -smoothed => ds::Pair( 
    $x, $smoothed,
    colors => pdl(0,255,0)->rgb_to_color,
    plotType => ppair::Lines,
  ),
);