# Make a four-view picture of the photometry of a luminaire

# This is inspired by objpict.pl that renders four-view

# Written by Axel Jacobs <axel@jaloxa.eu>

my $td = tempdir( CLEANUP => 1 );

use File::Copy qw(copy);

my $res = 1024;                # Default output image dimensions

my $maxsize = 0.001;           # max luminaire size after scaling

my $lumi = "$td/lumi.rad";     # Fitting as given on cmd line, or as generated by ies2rad

my $lumi2 = "$td/lumi2.rad";   # Fitting scaled to size

my $mat = "$td/lt.mat";

my $room1 = "$td/room1.rad";

my $room2 = "$td/room2.rad";

my $room3 = "$td/room3.rad";

my $room4 = "$td/room4.rad";

my $oct1 = "$td/lt1.oct";

my $oct2 = "$td/lt2.oct";

my $oct3 = "$td/lt3.oct";

my $oct4 = "$td/lt4.oct";

while ($#ARGV >= 0) {

	$_ = $ARGV[0];

	if (m/-i/) {         # File is an IES file, not Radiance luminaire

		$is_ies = 1;

	} elsif (m/-d/) {    # Resolution of the output HDR image

		$res = $ARGV[1];

		shift @ARGV;

    } elsif (m/^-/) { 	 # Oops! Illegal option

        die("ltpict: bad option \"$ARGV[0]\"\n");

    } else {

		last;    # No more options.  What's left is the actual file name.

    shift @ARGV;

# We need one Radiance luminaire or an IES file

if ($is_ies == 0) {

	# Input file is a Radiance luminaire

	copy $ARGV[0], $lumi or

			die("ltpict: Cannot copy file '$ARGV[0]' to $lumi\n");

} else {

	# Input file is IES photometry

	system "ies2rad -p $td -o lumi $ARGV[0]";

}

my $res2 = $res / 2;    # Each rendering is half the size of final composite

# Scale luminaire and center at origin

my $dimstr = `getbbox -h $lumi`;

$dimstr =~ s/^\s+//;                # remove leading spaces

# Move luminaire so centre is at origin, and scale

my $scale = $maxsize / $size;

open(FH, ">$lumi2") or

		die("ltpict: Cannot write to temporary file $lumi");

print FH "!xform -t $xtrans $ytrans $ztrans -s $scale $lumi";

open(FH, ">$mat") or

		die("ltpict: Cannot write to temporary file $mat");

my $offset = 0.1;

open(FH, ">$room1") or die("ltpict: Cannot write to temporary file $room1");

print FH "wall_mat polygon box.4620  0  0  12  -$offset -5 5  -$offset 5 5  -$offset 5 -5  -$offset -5 -5";

open(FH, ">$room2") or die("ltpict: Cannot write to temporary file $room2");

print FH "wall_mat polygon box.1540  0  0  12  5 $offset -5  5 $offset 5  -5 $offset 5  -5 $offset -5";

open(FH, ">$room3") or die("ltpict: Cannot write to temporary file $room3");

open(FH, ">$room4") or die("ltpict: Cannot write to temporary file $room4");

# Get bbox again, for the translated and scaled luminaire.

$dimstr = `getbbox -h $lumi2`;

system "oconv $mat $room1 $lumi2 > $oct1";

system "oconv $mat $room2 $lumi2 > $oct2";

system "oconv $mat $room3 $lumi2 > $oct3";

system "oconv $mat $room4 $lumi2 > $oct4";

my $rpict_cmd = "rpict -ab 0 -ds 0 -dv -av 0 0 0 -x $res2 -y $res2";

system "$rpict_cmd $vw1 $oct1 > $td/right.hdr";

system "$rpict_cmd $vw2 $oct2 > $td/front.hdr";

system "$rpict_cmd $vw3 $oct3 > $td/down.hdr";

system "$rpict_cmd $vw4 $oct4 > $td/up.hdr";

my $vtl = "$td/vtl.hdr";     # The two parallel views

my $vta = "$td/vta.hdr";     # The two fisheye views

my $pcond_cmd = "pcond -l";

system "pcompos -a 2 $td/right.hdr  $td/front.hdr > $td/rf.hdr";

system "pcompos -a 2 $td/down.hdr $td/up.hdr > $td/du.hdr";

# Combine top two images with bottom row.  Output HDR.