sub rgb_hash()  {
	$debug=0;
	$raw="pf-2015.0531-249465.pprgb.7360x4912.raw";
	%e2at=%rgb2c=();  # Init Event -> Accululated Time hash; RGB -> Count hash
	&time_event('init', \%e2at, $debug*0);  # Initialize timing structures
	&time_event('Read_Raw', \%e2at, $debug*0); 
	$rsize = -s $raw;
	open(IN, "<$raw");
	binmode IN;
	$read = read(IN, $buf, $rsize);
	close IN;  # Use it n cluse it!
	die("Read ERROR!: Read $read != file $rsize\n")  unless($read==$rsize);
	printf("Read %1.6f MB in raw file $raw\n", $read/1.0E6);
	&time_event('# Pre-extend array', \%e2at, $debug*0); 
	$ushort[108456960+1] = 0x0000FFFF0000;  # Pre-extend array
	&time_event('Unpack', \%e2at, $debug*0); 
	@ushort=unpack("S*", $buf);  $su = scalar(@ushort)  ||  -1;
	printf("Ushort[$su], length(buf)=%1.6f, Bytes/ushort %1.9f\n",
		($bl=length($buf)) / 1.0E6, $bl/$su);
	for($ii=0; $ii < scalar @ushort; $ii+=3)  {
		# Extract the R, G and B values
		($rr, $gg, $bb) = @ushort[$ii .. $ii+2];  # Array slice
		$quad=pack('Q', $rr, $gg, $bb, 0x0000);
		$mult=$rr*$gg*$bb;  # Titanic 281 Trillion size number
		$white=$rr + $gg*65535 + $bb*65535*65535;  # 281 Trillion size number
		$bs=$rr | ($gg << 16) | ($bb << 32);  # << WORKS!
		$rgb2c{$bs}++;  # Increment count for this color

		if(!($ii%6E6))  {  # Document time on every 100k unit
			$jj++;
			
			printf("R=0x%04hx, g=0x%04hx, b=0x%04hx, bs=0x%016qx, q=0x%016qx, m=0x%016qx, W=0x%016qx\n",
				$rr, $gg, $bb, $bs, $quad, $mult, $white);
		}
	}
	&time_event("Count distinct colors", \%e2at, $debug*0);   
	printf("RGB hash has %d distinct colors.\n", scalar keys %rgb2c);
	&time_event('term', \%e2at, $debug*0);  # Finalize timing structures
	&show_event(\%e2at, $debug*0);  # Automatically stops last event timing
}