/usr/bin/genskyvec is in radiance 4R1+20120125-1.1.
This file is owned by root:root, with mode 0o755.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 | #!/usr/bin/perl -w
# RCSid $Id: genskyvec.pl,v 2.7 2011/09/29 22:36:11 greg Exp $
#
# Generate Reinhart vector for a given sky description
#
# G. Ward
#
use strict;
my $windoz = ($^O eq "MSWin32" or $^O eq "MSWin64");
my @skycolor = (0.960, 1.004, 1.118);
my $mf = 4;
while ($#ARGV >= 0) {
if ("$ARGV[0]" eq "-c") {
@skycolor = @ARGV[1..3];
shift @ARGV; shift @ARGV; shift @ARGV;
} elsif ("$ARGV[0]" eq "-m") {
$mf = $ARGV[1];
shift @ARGV;
}
shift @ARGV;
}
# Load sky description into line array, separating sun if one
my @skydesc;
my $lightline;
my @sunval;
my $sunline;
my $skyOK = 0;
my $srcmod; # putting this inside loop breaks code(?!)
while (<>) {
push @skydesc, $_;
if (/^\w+\s+light\s+/) {
s/\s*$//; s/^.*\s//;
$srcmod = $_;
$lightline = $#skydesc;
} elsif (defined($srcmod) && /^($srcmod)\s+source\s/) {
@sunval = split(' ', $skydesc[$lightline + 3]);
shift @sunval;
$sunline = $#skydesc;
} elsif (/\sskyfunc\s*$/) {
$skyOK = 1;
}
}
die "Bad sky description!\n" if (! $skyOK);
# Strip out the solar source if present
my @sundir;
if (defined $sunline) {
@sundir = split(' ', $skydesc[$sunline + 3]);
shift @sundir;
undef @sundir if ($sundir[2] <= 0);
splice(@skydesc, $sunline, 5);
}
# Reinhart sky sample generator
my $rhcal = 'DEGREE : PI/180;' .
'x1 = .5; x2 = .5;' .
'alpha : 90/(MF*7 + .5);' .
'tnaz(r) : select(r, 30, 30, 24, 24, 18, 12, 6);' .
'rnaz(r) : if(r-(7*MF-.5), 1, MF*tnaz(floor((r+.5)/MF) + 1));' .
'raccum(r) : if(r-.5, rnaz(r-1) + raccum(r-1), 0);' .
'RowMax : 7*MF + 1;' .
'Rmax : raccum(RowMax);' .
'Rfindrow(r, rem) : if(rem-rnaz(r)-.5, Rfindrow(r+1, rem-rnaz(r)), r);' .
'Rrow = if(Rbin-(Rmax-.5), RowMax-1, Rfindrow(0, Rbin));' .
'Rcol = Rbin - raccum(Rrow) - 1;' .
'Razi_width = 2*PI / rnaz(Rrow);' .
'RAH : alpha*DEGREE;' .
'Razi = if(Rbin-.5, (Rcol + x2 - .5)*Razi_width, 2*PI*x2);' .
'Ralt = if(Rbin-.5, (Rrow + x1)*RAH, asin(-x1));' .
'Romega = if(.5-Rbin, 2*PI, if(Rmax-.5-Rbin, ' .
' Razi_width*(sin(RAH*(Rrow+1)) - sin(RAH*Rrow)),' .
' 2*PI*(1 - cos(RAH/2)) ) );' .
'cos_ralt = cos(Ralt);' .
'Dx = sin(Razi)*cos_ralt;' .
'Dy = cos(Razi)*cos_ralt;' .
'Dz = sin(Ralt);' ;
my ($nbins, $octree, $tregcommand, $suncmd);
if ($windoz) {
$nbins = `rcalc -n -e MF:$mf -e \"$rhcal\" -e \"\$1=Rmax+1\"`;
chomp $nbins;
$octree = "gtv$$.oct";
$tregcommand = "cnt $nbins 16 | rcalc -e MF:$mf -e \"$rhcal\" " .
q{-e "Rbin=$1;x1=rand(recno*.37-5.3);x2=rand(recno*-1.47+.86)" } .
q{-e "$1=0;$2=0;$3=0;$4=Dx;$5=Dy;$6=Dz" } .
"| rtrace -h -ab 0 -w $octree | total -16 -m";
if (@sundir) {
$suncmd = "cnt " . ($nbins-1) .
" | rcalc -e MF:$mf -e \"$rhcal\" -e Rbin=recno " .
"-e \"dot=Dx*$sundir[0] + Dy*$sundir[1] + Dz*$sundir[2]\" " .
"-e \"cond=dot-.866\" " .
q{-e "$1=if(1-dot,acos(dot),0);$2=Romega;$3=recno" };
}
} else {
$nbins = `rcalc -n -e MF:$mf -e \'$rhcal\' -e \'\$1=Rmax+1\'`;
chomp $nbins;
$octree = "/tmp/gtv$$.oct";
$tregcommand = "cnt $nbins 16 | rcalc -of -e MF:$mf -e '$rhcal' " .
q{-e 'Rbin=$1;x1=rand(recno*.37-5.3);x2=rand(recno*-1.47+.86)' } .
q{-e '$1=0;$2=0;$3=0;$4=Dx;$5=Dy;$6=Dz' } .
"| rtrace -h -ff -ab 0 -w $octree | total -if3 -16 -m";
if (@sundir) {
$suncmd = "cnt " . ($nbins-1) .
" | rcalc -e MF:$mf -e '$rhcal' -e Rbin=recno " .
"-e 'dot=Dx*$sundir[0] + Dy*$sundir[1] + Dz*$sundir[2]' " .
"-e 'cond=dot-.866' " .
q{-e '$1=if(1-dot,acos(dot),0);$2=Romega;$3=recno' };
}
}
# Create octree for rtrace
open OCONV, "| oconv - > $octree";
print OCONV @skydesc;
print OCONV "skyfunc glow skyglow 0 0 4 @skycolor 0\n";
print OCONV "skyglow source sky 0 0 4 0 0 1 360\n";
close OCONV;
# Run rtrace and average output for every 16 samples
my @tregval = `$tregcommand`;
unlink $octree;
# Find closest 3 patches to sun and divvy up direct solar contribution
sub numSort1 {
my @a1 = split("\t", $a);
my @b1 = split("\t", $b);
return ($a1[0] <=> $b1[0]);
}
if (@sundir) {
my $somega = ($sundir[3]/360)**2 * 3.141592654**3;
my @bestdir = `$suncmd`;
@bestdir = sort numSort1 @bestdir;
my (@ang, @dom, @ndx);
my $wtot = 0;
for my $i (0..2) {
($ang[$i], $dom[$i], $ndx[$i]) = split(' ', $bestdir[$i]);
$wtot += 1./($ang[$i]+.02);
}
for my $i (0..2) {
my $wt = 1./($ang[$i]+.02)/$wtot * $somega / $dom[$i];
my @scolor = split(' ', $tregval[$ndx[$i]]);
for my $j (0..2) { $scolor[$j] += $wt * $sunval[$j]; }
$tregval[$ndx[$i]] = "$scolor[0]\t$scolor[1]\t$scolor[2]\n";
}
}
# Output our final vector
print @tregval;
|