bin/objline.pl
author Axel Jacobs <axel@jaloxa.eu>
Sun, 26 Jun 2011 16:59:38 +0100
changeset 13 ccfb8a3410da
parent 2 d59fbab896ad
child 32 ea86d26c63ed
permissions -rwxr-xr-x
Made falsecolor.pl executable.
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
2
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
     1
#!/usr/bin/perl -w
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
     2
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
     3
use strict;
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
     4
use warnings;
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
     5
use File::Temp qw/ tempdir /;
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
     6
use File::Basename;
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
     7
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
     8
#my $td = tempdir( CLEANUP => 0 );
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
     9
my $td = "/tmp/objline";
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    10
mkdir $td;
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    11
#print "Temdir: $td\n";
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    12
my @rad = ();
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    13
my @prllmgf = ();
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    14
my @oblqmgf = ();
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    15
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    16
#my $name;
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    17
#my $path;
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    18
#my $extn;
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    19
my $fn = "";
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    20
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    21
my $oblqxf = "-rz 45 -ry 45";
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    22
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    23
if ($#ARGV < 0) {
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    24
	# Input is from STDIN: Capture to file
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    25
	open(FH, ">$td/stdin.rad");
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    26
	while (<>) {
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    27
		print FH;
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    28
	}
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    29
	close FH;
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    30
	# Pretend stdin.rad was passed as a filename
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    31
	@ARGV = ("$td/stdin.rad");
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    32
}
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    33
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    34
foreach (@ARGV) {
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    35
	my ($name, $path, $extn) = fileparse($_);
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    36
	$fn = $name.$extn;
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    37
	push(@rad, "$td/$fn.orig");
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    38
	push(@prllmgf, "$td/$fn.orig.mgf");
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    39
	push(@oblqmgf, "$td/$fn.oblq.mgf");
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    40
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    41
	# Expand in-line commands such as !genbox
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    42
	system "xform $_ > $td/$fn.orig";
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    43
	system "rad2mgf $td/$fn.orig > $td/$fn.orig.mgf";
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    44
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    45
	open (FH, ">$td/$fn.oblq.mgf");
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    46
	print FH "i $fn.orig.mgf $oblqxf > $td/$fn.oblq.mgf\n";
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    47
	close (FH);
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    48
}
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    49
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    50
my $rce = '\'xm=($1+$2)/2;ym=($3+$4)/2;zm=($5+$6)/2';
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    51
$rce .= ';max(a,b):if(a-b,a,b);r=max(max($2-$1,$4-$3),$6-$5)*.52';
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    52
$rce .= ';$1=xm-r;$2=xm+r;$3=ym-r;$4=ym+r;$5=zm-r;$6=zm+r\'';
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    53
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    54
my $rads = join(' ', @rad);
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    55
my $prlldim = `getbbox -h $rads |rcalc -t' ' -e $rce`;
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    56
chomp($prlldim);
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    57
print "prlldim: $prlldim\n";
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    58
my $oblqdim = `xform $oblqxf $rads |getbbox -h|rcalc -t' ' -e $rce`;
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    59
chomp($oblqdim);
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    60
print "oblqdim: $oblqdim\n";
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    61
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    62
my $prlls = join(' ', @prllmgf); #!!!
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    63
print "prlls: $prlls\n";
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    64
my $oblqs = join(' ', @oblqmgf);
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    65
my $cmd;
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    66
foreach ('x', 'y', 'z') {
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    67
	$cmd = "mgf2meta -t .005 $_ " . $prlldim . ' ' . $prlls . " >$td/$_.mta";
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    68
	#system "mgf2meta -t .005 $_ $prlldim $prlls >$td/$_.mta";
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    69
	system "$cmd";
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    70
}
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    71
$cmd = "mgf2meta -t .005 x " . $oblqdim . ' ' . $oblqs . " >$td/o.mta";
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    72
system "$cmd";
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    73
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    74
system("plot4", "$td/x.mta", "$td/y.mta", "$td/z.mta", "$td/o.mta");
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    75
d59fbab896ad Added objline.pl
Axel Jacobs <axel@jaloxa.eu>
parents:
diff changeset
    76
#EOF