kegg的某种svg展示图

#!perl use warnings; use strict; use SVG; use File::Basename qw(basename); die "perl $0 <sample.list> <pathway.sorted> <outprefix> Note: support 24 colors to sample~ " if @ARGV != 3; my (@file,@pwtmp,@pw2sp); my $i = 0; open PW,$ARGV[1] or die $!; while(<PW>) { chomp; my @tmp = split / /; $pwtmp[$i]{$tmp[0]} = 1; $i ++; } my $f = 0; open SP,$ARGV[0] or die $!; while(<SP>) { chomp; my @ft = split; my $name = basename($ft[0]); my ($id) = $ft[0] =~ /^([^.]+)/; push @file,$ft[1]; open FA,$ft[0] or die $!; <FA>; for(my $n = 0; $n < @pwtmp; $n ++) { $pw2sp[$f][$n] = 1; } while(my $line = <FA>) { chomp($line); my @tmp = split / /,$line; for(my $j = 0; $j < @pwtmp; $j ++) { if(exists $pwtmp[$j]{$tmp[0]}) { $pw2sp[$f][$j] = $tmp[4]; last; } } } $f ++; } my $fn = @file; my $pwn = @pwtmp; my $svg = SVG->new(width => ($fn * 50 + 1100),height => ($pwn * 50 + 300)); $svg->rect(x => 55,y => 55,width => (($fn + 1) * 50 + 10),height => (($pwn + 1) * 50 + 10),fill => "white",stroke => "grey",'stroke-width' => 8); for(my $j = 0; $j < @pwtmp; $j ++) { my $n = $j + 1; $svg->line(x1 => 60,y1 => (60 + $n * 50),x2 => (60 + ($fn+1) * 50),y2 => (60 + $n * 50),stroke => "orange",'stroke-width' => 3); foreach my $pc(keys %{$pwtmp[$j]}) { $svg->text(x => (($fn + 1) * 50 + 70),y => 65 + $n * 50,style => {'font' => 'Times new roman','font-size' => 30,'font-weight' => 'bold'})->cdata($pc); } } my @color = ('red','green','blue','yellow','cyan','purple','sandybrown','navy','violet','darksalmon','bisque','fuchsia','lightcoral','goldenrod','silver','brown','lawngreen','skyblue','indigo','forestgreen','lightseagreen','aqua','lime'); #my @color = ('red','red','yellow'); my $column = 0; my $row = 0; for(my $j = 0; $j < @file; $j ++) { my $n = $j + 1; my $rx = 75 + $n * 50 + int(length($file[$j]) / 2) * 14; my $ry = 80 + ($pwn + 1) * 50 + length($file[$j]) * 12; $svg->line(x1 => (60 + $n * 50),y1 => 60,x2 => (60 + $n * 50),y2 => (60 + ($pwn + 1) * 50),'stroke-width' => 3); $svg->text(x => (60 + $n * 50),y => (80 + ($pwn + 1) * 50),'font-size' => 25,'font-weight' => 'bold'},'transform' => "rotate(⑷5,$rx,$ry)")->cdata($file[$j]); for(my $k = 0; $k < @{$pw2sp[$j]}; $k ++) { if($pw2sp[$j][$k] >= 0.05) { $svg->circle(cx => (60 + $n * 50),cy => (60 + ($k+1) * 50),r => 1,fill => $color[$j],stroke => 'black'); }elsif($pw2sp[$j][$k] < 0.05 and $pw2sp[$j][$k] >= 0.01){ $svg->circle(cx => (60 + $n * 50),r => 20,stroke => 'black'); }else{ $svg->circle(cx => (60 + $n * 50),r => 45,stroke => 'black'); } } if(($n / 6) =~ /./) { $svg->circle(cx => (250 + (@file + 1) * 50 + $column * 100),cy => (75 + ($pwn + 1) * 50 + $row * 75),fill => "$color[$j]",stroke => 'black'); $svg->text(x => (250 + (@file + 1) * 50 + $column * 100),y => (110 + ($pwn + 1) * 50 + $row * 75),'font-size'=> 15,'font-weight' => 'bold'})->cdata($file[$j]); $column ++; }else{ $svg->circle(cx => (250 + (@file + 1) * 50 + $column * 100),'font-weight' => 'bold'})->cdata($file[$j]); $row ++; $column = 0; } } $svg->circle(cx => (100 + (@file + 1) * 50),cy => (75 + ($pwn + 1) * 50),fill => 'sandybrown',stroke => 'black'); $svg->circle(cx => (100 + (@file + 1) * 50),cy => (150 + ($pwn + 1) * 50),stroke => 'black'); $svg->text(x => (122 + (@file + 1) * 50),y => (75 + ($pwn + 1) * 50),'font-weight' => 'bold'})->cdata("FDR < 0.05"); $svg->text(x => (147 + (@file + 1) * 50),y => (150 + ($pwn + 1) * 50),'font-weight' => 'bold'})->cdata("FDR < 0.01"); open OUT,"> $ARGV[2].svg" or die $!; print OUT $svg->xmlify;