#!/usr/bin/perl

open(CT,"<califa_vfield_ct.txt");
while($line=<CT>) {
    chop($line);
    ($I,$r,$g,$b)=split(" ",$line);
    $a_I[$n]=$I/256;
    $a_r[$n]=$r/256;
    $a_g[$n]=$g/256;
    $a_b[$n]=$b/256;
    $n++;
}

$a_r[0]=$a_r[1];
$a_g[0]=$a_g[1];
$a_b[0]=$a_b[1];
$a_r[$n-1]=$a_r[$n-2];
$a_g[$n-1]=$a_g[$n-2];
$a_b[$n-1]=$a_b[$n-2];

close(CT);
open(LUT,">califa_vfield_ct.sao");
print LUT "# SAOimage color table\n";
print LUT "# Tue Feb 12 10:02:23 2013\n";
print LUT "PSEUDOCOLOR\n";
print LUT "RED:\n";
for ($i=0;$i<$n;$i=$i+25) {
    printf LUT ("(%3.1f,%3.1f)",$a_I[$i],$a_r[$i]); 
}
print LUT "\n";
print LUT "GREEN:\n";
for ($i=0;$i<$n;$i=$i+25) {
    printf LUT ("(%3.1f,%3.1f)",$a_I[$i],$a_g[$i]); 
}
print LUT "\n";
print LUT "BLUE:\n";
for ($i=0;$i<$n;$i=$i+25) {
    printf LUT ("(%3.1f,%3.1f)",$a_I[$i],$a_b[$i]); 
}
print LUT "\n";
close(LUT);