[37] | 1 | #!/usr/bin/perl |
---|
| 2 | # |
---|
| 3 | # vripsubvollist: This program, given a list of subvol .conf |
---|
| 4 | # files, will generate the commands to vrip each of them. |
---|
| 5 | |
---|
| 6 | |
---|
| 7 | sub printUsage { |
---|
| 8 | print STDERR "\n"; |
---|
| 9 | print STDERR "Usage:\n"; |
---|
| 10 | print STDERR "vripsubvollist [-crunch] [-rampscale s] voxelsize *subvol*.conf > outcommands\n"; |
---|
| 11 | print STDERR "\n"; |
---|
| 12 | print STDERR " -crunch Will generate plycrunched .set files...\n"; |
---|
| 13 | print STDERR " -rampscale <s> Will tell vrip to use rampscale s.\n"; |
---|
| 14 | print STDERR " -passtovrip \"X\" will pass X to vrip.\n"; |
---|
| 15 | print STDERR " -passtovripsurf \"X\" will pass X to vripsurf.\n"; |
---|
| 16 | print STDERR "\n"; |
---|
| 17 | print STDERR "and it will generate a list of the commands necessary\n"; |
---|
| 18 | print STDERR "to run vripnew, vripsurf, and plycull on the volumes.\n"; |
---|
| 19 | print STDERR "\n"; |
---|
| 20 | exit(-1); |
---|
| 21 | } |
---|
| 22 | |
---|
| 23 | if ($#ARGV < 1 || substr($ARGV[0], 0, 2) eq "-h") { |
---|
| 24 | &printUsage(); |
---|
| 25 | } |
---|
| 26 | |
---|
| 27 | # Extract the options, if present |
---|
| 28 | $DOCRUNCH = 0; |
---|
| 29 | while (substr($ARGV[0], 0,1) eq "-") { |
---|
| 30 | if ($ARGV[0] eq "-crunch") { |
---|
| 31 | $DOCRUNCH=1; |
---|
| 32 | splice(@ARGV, 0, 1); |
---|
| 33 | } elsif ($ARGV[0] eq "-rampscale") { |
---|
| 34 | $RAMPSCALE = $ARGV[1]; |
---|
| 35 | if ($RAMPSCALE eq "") { |
---|
| 36 | print STDERR "Error: no rampscale???\n\n"; |
---|
| 37 | &printUsage; |
---|
| 38 | } |
---|
| 39 | splice(@ARGV, 0, 2); |
---|
| 40 | } elsif ($ARGV[0] eq "-passtovrip") { |
---|
| 41 | $PASSTOVRIP = $ARGV[1]; |
---|
| 42 | splice(@ARGV, 0, 2); |
---|
| 43 | } elsif ($ARGV[0] eq "-passtovripsurf") { |
---|
| 44 | $PASSTOVRIPSURF = $ARGV[1]; |
---|
| 45 | splice(@ARGV, 0, 2); |
---|
| 46 | } else { |
---|
| 47 | print STDERR "Error: unhandled arg $ARGV[0]...\n"; |
---|
| 48 | &printUsage(); |
---|
| 49 | exit(-1); |
---|
| 50 | } |
---|
| 51 | } |
---|
| 52 | |
---|
| 53 | |
---|
| 54 | |
---|
| 55 | # Get and check voxel size |
---|
| 56 | $voxelsize = 0.0+$ARGV[0]; |
---|
| 57 | if ($voxelsize <= 0) { |
---|
| 58 | die("Error, voxelsize not greater than zero...\n"); |
---|
| 59 | } |
---|
| 60 | |
---|
| 61 | # Get directory |
---|
| 62 | $cwd = `pwd`; |
---|
| 63 | chop $cwd; |
---|
| 64 | |
---|
| 65 | #Generate the command for each file name... |
---|
| 66 | for ($n=1; $n <= $#ARGV; $n++) { |
---|
| 67 | $conf = $ARGV[$n]; |
---|
| 68 | $bbox = $conf; |
---|
| 69 | $bbox =~ s/.conf/.bbox.ply/g; |
---|
| 70 | $vri = $conf; |
---|
| 71 | $vri =~ s/.conf/.vri/g; |
---|
| 72 | $plyv = $conf; |
---|
| 73 | $plyv =~ s/.conf/_beyondbbox.ply/g; |
---|
| 74 | $ply = $conf; |
---|
| 75 | $ply =~ s/.conf/.ply/g; |
---|
| 76 | # Generate the command. Include "npri -h 250", which sets |
---|
| 77 | # its priority really low... |
---|
| 78 | if ($DOCRUNCH) { |
---|
| 79 | $crunchstr = " ply2crunchset -l 5 $ply;"; |
---|
| 80 | } |
---|
| 81 | if (defined($RAMPSCALE)) { |
---|
| 82 | $rampscalestr = " -rampscale $RAMPSCALE "; |
---|
| 83 | } else { |
---|
| 84 | $rampscalestr = ""; |
---|
| 85 | } |
---|
| 86 | |
---|
| 87 | # Remove the plyv (overlapping ply files?) |
---|
| 88 | $rm_plyv = ""; |
---|
| 89 | # $rm_plyv = "/bin/rm $plyv;"; |
---|
| 90 | |
---|
| 91 | # Amount of safety in plyculling triangles |
---|
| 92 | $epsilon = 2.1 * $voxelsize; |
---|
| 93 | |
---|
| 94 | print STDOUT "cd $cwd; ". |
---|
| 95 | "time vripnew $vri $conf $bbox $voxelsize $rampscalestr $PASSTOVRIP -noui; ". |
---|
| 96 | "time vripsurf -no_remove_slivers $vri $plyv $PASSTOVRIPSURF -noui; ". |
---|
| 97 | "time plycull $bbox $plyv $ply $epsilon; ". |
---|
| 98 | "$rm_plyv $crunchstr\n"; |
---|
| 99 | $plys = $plys . " $ply"; |
---|
| 100 | } |
---|
| 101 | |
---|
| 102 | # print STDOUT "echo plymerge $plys > all_unshared.ply\n"; |
---|
| 103 | # print STDOUT "echo plyshared < all_unshared.ply > all.ply\n"; |
---|