#!/usr/bin/perl
use PDL;
use PDL::NiceSlice;

sub printstats {
       my $image = shift;
       my ($mean, $rms, $median, $min, $max) = stats($image);
       printf "mean: %.3g  rms: %.3g  median: %.3g  min: %.3g  max:
%.3g\n", $mean, $rms, $median, $min, $max;
}

($x,$y)=rcols('ct007.1C_1003C273.1.gpfit.0.2',10,11);

$az16=$x(0:5);
$el16=$y(0:5);
$az79=$x(6:8);
$el79=$y(6:8);

printstats($az16);
printstats($el16);
printstats($az79);
printstats($el79);