| 1 | #!/usr/bin/perl -w | 
| 2 | # RCSid $Id: rtpict.pl,v 2.17 2020/04/08 02:39:30 greg Exp $ | 
| 3 | # | 
| 4 | # Run rtrace in parallel mode to simulate rpict -n option | 
| 5 | # May also be used to render layered images with -o* option | 
| 6 | # | 
| 7 | #       G. Ward | 
| 8 | # | 
| 9 | use strict; | 
| 10 | # we'll run rpict if no -n or -o* option | 
| 11 | my $nprocs = 1; | 
| 12 | # rtrace options and the associated number of arguments | 
| 13 | my %rtraceC = ('-dt',1, '-dc',1, '-dj',1, '-ds',1, '-dr',1, '-dp',1, | 
| 14 | '-ss',1, '-st',1, '-e',1, '-am',1, '-P',1, '-PP',1, | 
| 15 | '-ab',1, '-af',1, '-ai',1, '-aI',1, '-ae',1, '-aE',1, | 
| 16 | '-av',3, '-aw',1, '-aa',1, '-ar',1, '-ad',1, '-as',1, | 
| 17 | '-me',3, '-ma',3, '-mg',1, '-ms',1, '-lr',1, '-lw',1, | 
| 18 | '-ap',2, '-am',1, '-ac',1, '-aC',1); | 
| 19 | # boolean rtrace options | 
| 20 | my @boolO = ('-w', '-bv', '-dv', '-i', '-u'); | 
| 21 | # view options and the associated number of arguments | 
| 22 | my %vwraysC = ('-vf',1, '-vtv',0, '-vtl',0, '-vth',0, '-vta',0, '-vts',0, '-vtc',0, | 
| 23 | '-x',1, '-y',1, '-vp',3, '-vd',3, '-vu',3, '-vh',1, '-vv',1, | 
| 24 | '-vo',1, '-va',1, '-vs',1, '-vl',1, '-pa',1, '-pj',1); | 
| 25 | # options we need to silently ignore | 
| 26 | my %ignoreC = ('-t',1, '-ps',1, '-pt',1, '-pm',1, '-pd',1); | 
| 27 | # Starting options for rtrace (rpict values) | 
| 28 | my @rtraceA = split(' ', 'rtrace -u- -dt .05 -dc .5 -ds .25 -dr 1 ' . | 
| 29 | '-aa .2 -ar 64 -ad 512 -as 128 -lr 7 -lw 1e-03'); | 
| 30 | my @vwraysA = ('vwrays', '-ff', '-pj', '.67'); | 
| 31 | my @vwrightA = ('vwright', '-vtv'); | 
| 32 | my @rpictA = ('rpict', '-ps', '1'); | 
| 33 | my $outpatt = '^-o[vrxlLRXnNsmM]+'; | 
| 34 | my $refDepth = ""; | 
| 35 | my $irrad = 0; | 
| 36 | my $persist = 0; | 
| 37 | my $outlyr; | 
| 38 | my $outdir; | 
| 39 | my $outpic; | 
| 40 | my $outzbf; | 
| 41 | OPTION:                         # sort through options | 
| 42 | while ($#ARGV >= 0 && "$ARGV[0]" =~ /^[-\@]/) { | 
| 43 | # Check for file inclusion | 
| 44 | if ("$ARGV[0]" =~ /^\@/) { | 
| 45 | open my $handle, '<', substr($ARGV[0], 1) or die "No file for $ARGV[0]\n"; | 
| 46 | shift @ARGV; | 
| 47 | chomp(my @args = <$handle>); | 
| 48 | close $handle; | 
| 49 | unshift @ARGV, split(/\s+/, "@args"); | 
| 50 | next OPTION; | 
| 51 | } | 
| 52 | # Check booleans | 
| 53 | for my $boopt (@boolO) { | 
| 54 | if ("$ARGV[0]" =~ ('^' . $boopt . '[-+01tfynTFYN]?$')) { | 
| 55 | if ("$ARGV[0]" eq '-i') { | 
| 56 | $irrad = ! $irrad; | 
| 57 | } elsif ("$ARGV[0]" =~ /^-i[-+01tfynTFYN]/) { | 
| 58 | $irrad = ("$ARGV[0]" =~ /^-i[+1tyTY]/); | 
| 59 | } | 
| 60 | push @rtraceA, $ARGV[0]; | 
| 61 | push @rpictA, shift(@ARGV); | 
| 62 | next OPTION; | 
| 63 | } | 
| 64 | } | 
| 65 | # Check view options | 
| 66 | if (defined $vwraysC{$ARGV[0]}) { | 
| 67 | push @vwraysA, $ARGV[0]; | 
| 68 | my $isvopt = ("$ARGV[0]" =~ /^-v/); | 
| 69 | push(@vwrightA, $ARGV[0]) if ($isvopt); | 
| 70 | push @rpictA, shift(@ARGV); | 
| 71 | for (my $i = $vwraysC{$rpictA[-1]}; $i-- > 0; ) { | 
| 72 | push @vwraysA, $ARGV[0]; | 
| 73 | push(@vwrightA, $ARGV[0]) if ($isvopt); | 
| 74 | push @rpictA, shift(@ARGV); | 
| 75 | } | 
| 76 | next OPTION; | 
| 77 | } | 
| 78 | # Check rtrace options | 
| 79 | if (defined $rtraceC{$ARGV[0]}) { | 
| 80 | $persist ||= ("$ARGV[0]" =~ /^-PP?$/); | 
| 81 | push @rtraceA, $ARGV[0]; | 
| 82 | push @rpictA, shift(@ARGV); | 
| 83 | for (my $i = $rtraceC{$rpictA[-1]}; $i-- > 0; ) { | 
| 84 | push @rtraceA, $ARGV[0]; | 
| 85 | push @rpictA, shift(@ARGV); | 
| 86 | } | 
| 87 | next OPTION; | 
| 88 | } | 
| 89 | # Check options to ignore | 
| 90 | if (defined $ignoreC{$ARGV[0]}) { | 
| 91 | push @rpictA, shift(@ARGV); | 
| 92 | for (my $i = $ignoreC{$rpictA[-1]}; $i-- > 0; ) { | 
| 93 | push @rpictA, shift(@ARGV); | 
| 94 | } | 
| 95 | next OPTION; | 
| 96 | } | 
| 97 | # Check for output file or number of processes | 
| 98 | if ("$ARGV[0]" eq '-o') { | 
| 99 | shift @ARGV; | 
| 100 | $outpic = shift(@ARGV); | 
| 101 | } elsif ("$ARGV[0]" eq '-z') { | 
| 102 | push @rpictA, shift(@ARGV); | 
| 103 | $outzbf = $ARGV[0]; | 
| 104 | push @rpictA, shift(@ARGV); | 
| 105 | } elsif ("$ARGV[0]" eq '-n') { | 
| 106 | shift @ARGV; | 
| 107 | $nprocs = shift(@ARGV); | 
| 108 | } elsif ("$ARGV[0]" eq '-d') { | 
| 109 | shift @ARGV; | 
| 110 | $refDepth = ' -d ' . shift(@ARGV); | 
| 111 | } elsif ("$ARGV[0]" =~ "$outpatt") { | 
| 112 | push @rtraceA, $ARGV[0]; | 
| 113 | $outlyr = substr(shift(@ARGV), 2); | 
| 114 | $outdir = shift(@ARGV); | 
| 115 | } else { | 
| 116 | die "Unsupported option: " . $ARGV[0] . "\n"; | 
| 117 | } | 
| 118 | } | 
| 119 | die "Number of processes must be positive" if ($nprocs <= 0); | 
| 120 | if (defined $outdir) {          # check conflicting options | 
| 121 | die "Options -o and -o* are mutually exclusive\n" if (defined $outpic); | 
| 122 | die "Options -z and -o* are mutually exclusive\n" if (defined $outzbf); | 
| 123 | } | 
| 124 | if (defined $outpic) {          # redirect output? | 
| 125 | die "File '$outpic' already exists\n" if (-e $outpic); | 
| 126 | open STDOUT, '>', "$outpic"; | 
| 127 | } | 
| 128 | ##################################################################### | 
| 129 | ##### May as well run rpict? | 
| 130 | if ($nprocs == 1 && $persist == 0 && !defined($outdir)) { | 
| 131 | push(@rpictA, $ARGV[0]) if ($#ARGV == 0); | 
| 132 | exec @rpictA ; | 
| 133 | } | 
| 134 | # OK, we're running rtrace in some capacity... | 
| 135 | push @rtraceA, ('-n', "$nprocs"); | 
| 136 | die "Need single octree argument\n" if ($#ARGV != 0); | 
| 137 | my $oct = $ARGV[0]; | 
| 138 | my $view = `@vwrightA 0`; | 
| 139 | chomp $view; | 
| 140 | my @res = split(/\s/, `@vwraysA -d`); | 
| 141 | ##################################################################### | 
| 142 | ##### Generating picture with depth buffer? | 
| 143 | if (defined $outzbf) { | 
| 144 | exec "@vwraysA | @rtraceA -fff -olv @res '$oct' | " . | 
| 145 | "rsplit -ih -iH -f -of '$outzbf' -oh -oH -of3 - | " . | 
| 146 | "pvalue -r -df | getinfo -a 'VIEW=$view'"; | 
| 147 | } | 
| 148 | ##################################################################### | 
| 149 | ##### Base case with output picture only? | 
| 150 | if (! defined $outdir) { | 
| 151 | exec "@vwraysA | @rtraceA -ffc @res '$oct' | getinfo -a 'VIEW=$view'"; | 
| 152 | } | 
| 153 | ##################################################################### | 
| 154 | ##### Layered image output case | 
| 155 | my @rsplitA = ("rsplit", "'-t   '", "-ih", "-iH", "-oh", "-oH"); | 
| 156 | # Supported rtrace -o* options and output file name.typ | 
| 157 | my %rtoutC = ( | 
| 158 | v =>    'radiance.hdr', | 
| 159 | r =>    'r_refl.hdr', | 
| 160 | x =>    'r_unrefl.hdr', | 
| 161 | l =>    'd_effective.dpt', | 
| 162 | L =>    'd_firstsurf.dpt', | 
| 163 | R =>    'd_refl.dpt', | 
| 164 | X =>    'd_unrefl.dpt', | 
| 165 | n =>    'perturbed.nrm', | 
| 166 | N =>    'unperturbed.nrm', | 
| 167 | s =>    'surface.idx', | 
| 168 | m =>    'modifier.idx', | 
| 169 | M =>    'material.idx' | 
| 170 | ); | 
| 171 | # Arguments for rsplit based on output file type | 
| 172 | my %rcodeC = ( | 
| 173 | '.hdr', ['-of3', '!pvalue -r -df -u'], | 
| 174 | '.dpt', ['-of', "!rcode_depth$refDepth -ff"], | 
| 175 | '.nrm', ['-of3', '!rcode_norm -ff'], | 
| 176 | '.idx', ['-oa', '!rcode_ident "-t       "'] | 
| 177 | ); | 
| 178 | if ($irrad) {           # adjust actions for -i+ option | 
| 179 | $rtoutC{'v'} = 'irradiance.hdr'; | 
| 180 | delete $rtoutC{'r'}; | 
| 181 | delete $rtoutC{'x'}; | 
| 182 | delete $rtoutC{'R'}; | 
| 183 | delete $rtoutC{'X'}; | 
| 184 | } | 
| 185 | if (! -d $outdir) { | 
| 186 | mkdir $outdir or die("Cannot create output directory '$outdir'\n"); | 
| 187 | } | 
| 188 | # construct rsplit command | 
| 189 | foreach my $oval (split //, $outlyr) { | 
| 190 | die "Duplicate or unsupported type '$oval' in -o$outlyr\n" | 
| 191 | if (!defined $rtoutC{$oval}); | 
| 192 | my $outfile = "$outdir/$rtoutC{$oval}"; | 
| 193 | die "File '$outfile' already exists\n" if (-e $outfile); | 
| 194 | my ($otyp) = ($rtoutC{$oval} =~ /(\.[^.]+)$/); | 
| 195 | push @rsplitA, $rcodeC{$otyp}[0]; | 
| 196 | push @rsplitA, qq{'$rcodeC{$otyp}[1] > "$outfile"'}; | 
| 197 | delete $rtoutC{$oval}; | 
| 198 | } | 
| 199 | # call rtrace + rsplit | 
| 200 | exec "@vwraysA | @rtraceA -fff @res '$oct' | getinfo -a 'VIEW=$view' | @rsplitA"; |