petal_plot_v4RG.pl

#!/usr/bin/perl 
eval "use GD";
use FindBin 1.51 qw( $RealBin );
use lib $RealBin;
use GD;
use Math::Bezier;
use strict;
use Math::VecStat qw(min max);
use strict;
use warnings;
use Data::Dumper qw(Dumper);
use List::MoreUtils qw(uniq minmax );
use xytoradial;
use drawimage;  
use GD::Simple;
use GD::SVG;
#use GD::Font;
use SVG;
use POSIX;
use mycolors;
use Color::Rgb;
my @argume=@ARGV;
my $infile; 
my $outfile;

if (@argume < 1){ die "Please provide .txt input file"; }
if (@argume == 2){
($infile, $outfile) = @argume;
}else{
  $infile=$argume[0];
  $outfile= $infile =~ s/.txt/.svg/r;
  #$outfile= $infile =~ s/.txt/.png/r;
}
my $rgb = new Color::Rgb(rgb_txt=>'rgb.txt');
my $filename = $infile;#'CD3_data.txt';
open(my $fh, '<:encoding(UTF-8)', $filename)
  or die "Could not open file '$filename' $!";

my %point_data; 
my @temp_axis_names; 
my %axis_names;
my @allvalues;
my @rowarray; 
my @defaul_im_size=(10000,10000);
#-------------------x(c1),y(c2)----------

my $outer_cir= ($defaul_im_size[0]/2)-($defaul_im_size[0]/2)*0.1;#450; #radius of outer circle upper limit 
my $inner_cir=($defaul_im_size[0])*0.1;#100; #radius of inner circle lower limit
my $outer_axis=($defaul_im_size[0]/2)*0.95;#490; #radius of outer circle upper limit 
my $inner_axis=($defaul_im_size[0]/2)*0.15;#49; #radius of inner circle lower limit
my @maxmin_lim= (0,1);
my $col="blue";
my @colorcode=mycolors::findcolor($col);
#print Dumper \@colorcode;
@allvalues=@maxmin_lim;
my %clone_color_code;
while (my $row = <$fh>) {
  chomp $row;
  @rowarray= split /\s+/,$row;
  push @allvalues,($rowarray[2],$rowarray[3]);
  push @temp_axis_names, $rowarray[1];
  push @{$point_data{$rowarray[0]}{$rowarray[1]}}, @rowarray[2..3];
  #@{$clone_color_code{$rowarray[0]}}=mycolors::findcolor($rowarray[4]);
  my @rghx = map $_ , unpack 'C*', pack 'H*', $rowarray[4];
  @{$clone_color_code{$rowarray[0]}}=$rgb->hex2rgb($rowarray[4]);
  #print Dumper @rghx;
  #$point_data{$rowarray[0]}{$rowarray[1]}= \@rowarray;
  
}

 
#print Dumper \%point_data;
my @chekaxis_names = uniq @temp_axis_names;
my $size=@chekaxis_names;
#print "$size\n";
if ($size<3){ 
   die "Less than 3 defined";
}

#print "-----------------##################------------------\n";
for my $ky (@temp_axis_names){
    push @{$axis_names{$ky}}, $ky;
}

#print Dumper \%axis_names;
my @minmax = (minmax @allvalues);
if($minmax[0]>$maxmin_lim[0]){ die "Minimum limit is higher than data minimum"; }
if($minmax[1]<$maxmin_lim[1]){ die "Maximum limit is higher than data maximum"; }
#print "$minmax[0]\n";
@maxmin_lim=(0,$maxmin_lim[1]);
my %axs_ang = xytoradial::get_axis_angles(%axis_names);



#print Dumper \@test;

#print "$axs_ang{'a1'}[0]\n";

my $img=drawimage::createimage(@defaul_im_size);
my $blue = $img->colorAllocate(0,0,255);
drawimage::circ_out(($defaul_im_size[0],$defaul_im_size[1],$inner_cir,$outer_cir,$img));

#my $poly = new GD::Polygon;
#$img->string(gdSmallFont,2,10,"Peachy Keen",$blue);
#$img->openPolygon($poly,$colx);
my %axis_col= ('a1'=>"black",#"yellow",
                'a2'=>"black",
                'a3'=>"black",);#"cyan");
#print Dumper \%point_data;
foreach my $outkey (keys %axs_ang) {
  
  my ($x1,$y1,$x2,$y2)=drawimage::axisdraw(($defaul_im_size[0],$defaul_im_size[1],$outer_axis,$inner_axis,$axs_ang{$outkey}[1],$img));
  my ($reg,$green,$blue)=mycolors::findcolor($axis_col{$outkey});
  my $alpha=0;
  my $bl = $img->colorAllocateAlpha($reg,$green,$blue,$alpha);
  $img ->setThickness(40);
  $img->line($x1,$y1,$x2,$y2,$bl);
  #print Dumper \($x1,$y1,$x2,$y2);
  foreach my $innerkey (keys %{$point_data{$outkey}}) {
    my $tempvalue=$point_data{$outkey}{$innerkey}[1];
    my $norm_val= $inner_cir+($tempvalue-$minmax[0])/($minmax[1]-$minmax[0])*($outer_cir-$inner_cir);
    my %temhash=('im_size'=> [@defaul_im_size],
                  'value'=> $norm_val,
                  'theta'=> $axs_ang{$outkey});
  
  my @xycord=xytoradial::get_xy_loc(%temhash);
  my $crs=$point_data{$outkey}{$innerkey}[3];
  #drawimage::addPnt(($xycord[0],$xycord[1],$point_data{$outkey}{$innerkey}[2],$crs,$img));
  #print Dumper \%temhash;
  }
}
my %clone_stack_cord;
my $stck_width=($defaul_im_size[0])*0.04;#100;
my $stck_widthBG=($defaul_im_size[0])*0.08;#100;
my $no_curve_pt=100;
my %point_data2 = %point_data;
for my $pkeys (keys %point_data){
my %tem=%{$point_data{$pkeys}};
#print Dumper \%tem;
for my $tmkey (keys %tem){
my %temrecod;
my @temval=@{$tem{$tmkey}};
my $nval1=$inner_cir+($temval[0]-$minmax[0])/($minmax[1]-$minmax[0])*($outer_cir-$inner_cir);
my $nval2=$inner_cir+($temval[1]-$minmax[0])/($minmax[1]-$minmax[0])*($outer_cir-$inner_cir);

%temrecod=('im_size'=> [@defaul_im_size],
            'value' => [$nval1,$nval2],
            'theta' => $axs_ang{$tmkey}[1],
            'width' => $stck_width/2
);
#print Dumper \%temrecod;
my @points=(xytoradial::get_stack_cord(%temrecod));
# @points = ($x1c1,$y1b1,$x1,$y1,$x1c2,$y1b2,$x2c2,$y2b2,$x2,$y2,$x2c1,$y2b1)
#               0    1     2  3    4     5     6     7    8   9   10    11 
push @{$clone_stack_cord{$pkeys}{$tmkey}}, @points;
my @test;
my $alpha=0;
$img ->setThickness(1);
my $bl = $img->colorAllocateAlpha(@{$clone_color_code{$pkeys}},$alpha);
my $poly = new GD::Polygon;
#my $poly = GD::SVG::Polygon->new;
while (@points) {
  my @newpnt= splice( @points, 0, 2 );
  #print "($newpnt[0], $newpnt[1])\n";
  #$poly->addPt($newpnt[0], $newpnt[1]);
}
#$img->filledPolygon($poly,$bl);
}
}
my %clone_ribion;
my $shif_angle=60;
my $angle_steps=5;
for my $clonkey (keys %clone_stack_cord){
  my %temclone = %{$clone_stack_cord{$clonkey}};
  #print Dumper \%temclone;
  my @all_ax = sort(keys %temclone);
   if (@all_ax == 1){
    #continue;
  }
  if (@all_ax == 2){
    my ($ax1, $ax2)= @all_ax;
    my @ax1vari = @{$temclone{$ax1}};
    my @ax2vari = @{$temclone{$ax2}};
    my ($bs1x1,$bs1y1,$bs1x2,$bs1y2);
    my ($bs2x1,$bs2y1,$bs2x2,$bs2y2);
    my ($c1x,$b1x,$c2x,$b2x);
    if (($ax1 eq "a1") && ($ax2 ne "a2")){
            # @points = ($x1c1,$y1b1,$x1,$y1,$x1c2,$y1b2,$x2c2,$y2b2,$x2,$y2,$x2c1,$y2b1)
      #               0    1     2  3    4     5     6     7    8   9   10    11 
       ($c1x,$b1x,$c2x,$b2x) = ($ax1vari[0]-$ax1vari[2],$ax1vari[1]-$ax1vari[3],$ax1vari[4]-$ax1vari[2],$ax1vari[5]-$ax1vari[3]);
       
        my %tem1=('axis' => $ax1,
                  'im_size'=> [@defaul_im_size],
                  'cords'=> [$ax1vari[2],$ax1vari[3],$ax1vari[8],$ax1vari[9]],
                  'bcs' => [-$c1x,-$b1x,-$c2x,-$b2x],
                  'theta'=> $axs_ang{$ax1}[1],
                  'factor' => -1,
                  'width' => $stck_width/2,
                  'steps' => $angle_steps,
                  'tranthe' => $axs_ang{$ax1}[1],#-$shif_angle,
                  'shiftangele'=> $shif_angle);
        #($bs1x1,$bs1y1,$bs1x2,$bs1y2)
        my %bezhash1=xytoradial::get_rib_corda3a1(%tem1);
        
        ($c1x,$b1x,$c2x,$b2x) = ($ax2vari[0]-$ax2vari[2],$ax2vari[1]-$ax2vari[3],$ax2vari[4]-$ax2vari[2],$ax2vari[5]-$ax2vari[3]);
        my %tem2=('axis' => $ax2,
                  'im_size'=> [@defaul_im_size],
                  'cords'=> [$ax2vari[2],$ax2vari[3],$ax2vari[8],$ax2vari[9]],
                  'bcs' => [$c1x,$b1x,$c2x,$b2x],
                  'theta'=> $axs_ang{$ax2}[1],
                  'factor' => 1,
                  'width' => $stck_width/2,
                  'steps' => $angle_steps,
                  'tranthe' => $axs_ang{$ax2}[1],#+$shif_angle,
                  'shiftangele'=> $shif_angle);
        #($bs2x1,$bs2y1,$bs2x2,$bs2y2)
        my %bezhash2=xytoradial::get_rib_corda3a1(%tem2);
        my @points1;
        my @bez_input1=($ax1vari[4],$ax1vari[5]);
        push @bez_input1,@{$bezhash1{'x1y1'}};
        push @bez_input1,@{$bezhash2{'x1y1'}};
        push @bez_input1,($ax2vari[0],$ax2vari[1]);
        #my $bezier1 = Math::Bezier->new($ax1vari[4],$ax1vari[5],$bs1x1,$bs1y1,$bs2x1,$bs2y1,$ax2vari[0],$ax2vari[1]);
        my $bezier1 = Math::Bezier->new(@bez_input1);
        push @points1, ($ax1vari[4],$ax1vari[5]);
        push @points1, $bezier1->curve($no_curve_pt);
        push @points1, ($ax2vari[0],$ax2vari[1]);
        my @bezier_points;
        while (@points1) {
            push @bezier_points, [ splice( @points1, 0, 2 ) ];
        }
        #push @{$clone_ribion{$clonkey}{'rib1'}}, @points1;
        push @bezier_points, [$ax2vari[10],$ax2vari[11]];
        my @points2;
        my @bez_input2=($ax2vari[10],$ax2vari[11]);
        push @bez_input2,@{$bezhash2{'x2y2'}};
        push @bez_input2,@{$bezhash1{'x2y2'}};
        push @bez_input2,($ax1vari[6],$ax1vari[7]);
        #my $bezier2 = Math::Bezier->new($ax2vari[10],$ax2vari[11],$bs2x2,$bs2y2,$bs1x2,$bs1y2,$ax1vari[6],$ax1vari[7]);
        my $bezier2 = Math::Bezier->new(@bez_input2);
        push @points2, ($ax2vari[10],$ax2vari[11]);
        @points2=$bezier2->curve($no_curve_pt);
        push @points2, ($ax1vari[6],$ax1vari[7]);
        my @bezier2_points2;
        while (@points2) {
            push @bezier_points, [ splice( @points2, 0, 2 ) ];
        }
        push @{$clone_ribion{$clonkey}{'rib1'}}, @bezier_points;
        #push @{$clone_ribion{$clonkey}{'rib1'}},( @points1,@points2);
        #print "$ax1\t$tem1{'theta'}[1]\t$ax2\t$tem2{'theta'}[1]\n";
    }else{
      # @points = ($x1c1,$y1b1,$x1,$y1,$x1c2,$y1b2,$x2c2,$y2b2,$x2,$y2,$x2c1,$y2b1)
      #               0    1     2  3    4     5     6     7    8   9   10    11 
       ($c1x,$b1x,$c2x,$b2x) = ($ax1vari[0]-$ax1vari[2],$ax1vari[1]-$ax1vari[3],$ax1vari[4]-$ax1vari[2],$ax1vari[5]-$ax1vari[3]);
        my %tem1=('axis' => $ax1,
                  'im_size'=> [@defaul_im_size],
                  'cords'=> [$ax1vari[2],$ax1vari[3],$ax1vari[8],$ax1vari[9]],
                  'bcs' => [$c1x,$b1x,$c2x,$b2x],
                  'theta'=> $axs_ang{$ax1}[1],
                  'factor' => 1,
                  'width' => $stck_width/2,
                  'steps' => $angle_steps,
                  'tranthe' => $axs_ang{$ax1}[1],#+$shif_angle,
                  'shiftangele'=> $shif_angle);
        
        #($bs1x1,$bs1y1,$bs1x2,$bs1y2)=xytoradial::get_rib_cord(%tem1);
        my %bezhash1=xytoradial::get_rib_cord(%tem1);
        ($c1x,$b1x,$c2x,$b2x) = ($ax2vari[0]-$ax2vari[2],$ax2vari[1]-$ax2vari[3],$ax2vari[4]-$ax2vari[2],$ax2vari[5]-$ax2vari[3]);
        my %tem2=('axis' => $ax2,
                  'im_size'=> [@defaul_im_size],
                  'cords'=> [$ax2vari[2],$ax2vari[3],$ax2vari[8],$ax2vari[9]],
                  'bcs' => [-$c1x,-$b1x,-$c2x,-$b2x],
                  'theta'=> $axs_ang{$ax2}[1],
                  'factor' => -1,
                  'width' => $stck_width/2,
                  'steps' => $angle_steps,
                  'tranthe' => $axs_ang{$ax2}[1],#-$shif_angle,
                  'shiftangele'=> $shif_angle);
        
        #($bs2x1,$bs2y1,$bs2x2,$bs2y2)=xytoradial::get_rib_cord(%tem2);
        my %bezhash2=xytoradial::get_rib_cord(%tem2);
        my @points1;
        my @bez_input1=($ax1vari[0],$ax1vari[1]);
        push @bez_input1,@{$bezhash1{'x1y1'}};
        push @bez_input1,@{$bezhash2{'x1y1'}};
        push @bez_input1,($ax2vari[4],$ax2vari[5]);
        
        my $bezier1 = Math::Bezier->new(@bez_input1);
        #my $bezier1 = Math::Bezier->new($ax1vari[0],$ax1vari[1],$bs1x1,$bs1y1,$bs2x1,$bs2y1,$ax2vari[4],$ax2vari[5]);
        push @points1, ($ax1vari[0],$ax1vari[1]);
        push @points1, $bezier1->curve($no_curve_pt);
        push @points1, ($ax2vari[4],$ax2vari[5]);
        my @bezier_points;
        while (@points1) {
            push @bezier_points, [ splice( @points1, 0, 2 ) ];
        }
        #push @{$clone_ribion{$clonkey}{'rib1'}}, @points1;
        push @bezier_points, [$ax2vari[6],$ax2vari[7]];
        my @points2;
        my @bez_input2=($ax2vari[6],$ax2vari[7]);
        push @bez_input2,@{$bezhash2{'x2y2'}};
        push @bez_input2,@{$bezhash1{'x2y2'}};
        push @bez_input2,($ax1vari[10],$ax1vari[11]);
        
        my $bezier2 = Math::Bezier->new(@bez_input2);
        #my $bezier2 = Math::Bezier->new($ax2vari[6],$ax2vari[7],$bs2x2,$bs2y2,$bs1x2,$bs1y2,$ax1vari[10],$ax1vari[11]);
        
        push @points2, ($ax2vari[6],$ax2vari[7]);
        @points2=$bezier2->curve($no_curve_pt);
        push @points2, ($ax1vari[10],$ax1vari[11]);
        my @bezier2_points2;
        while (@points2) {
            push @bezier_points, [ splice( @points2, 0, 2 ) ];
        }
        push @{$clone_ribion{$clonkey}{'rib1'}}, @bezier_points;
        #push @{$clone_ribion{$clonkey}{'rib1'}},( @points1,@points2);
        #print "$ax1\t$tem1{'theta'}[1]\t$ax2\t$tem2{'theta'}[1]\n";
    }
  }


  if (@all_ax == 3){
    my %axishash=('a1a2'=> ['a1','a2'],
                  'a2a3'=> ['a2','a3'],
                  'a3a1'=> ['a1','a3']);
    for my $ribkey (keys %axishash){
    #print "$ribkey\n";    
    my ($ax1, $ax2)= sort(@{$axishash{$ribkey}});  
    
    my @ax1vari = @{$temclone{$ax1}};
    my @ax2vari = @{$temclone{$ax2}};
    my ($bs1x1,$bs1y1,$bs1x2,$bs1y2);
    my ($bs2x1,$bs2y1,$bs2x2,$bs2y2);
    my ($c1x,$b1x,$c2x,$b2x);
    if (($ax1 eq "a1") && ($ax2 ne "a2")){
            # @points = ($x1c1,$y1b1,$x1,$y1,$x1c2,$y1b2,$x2c2,$y2b2,$x2,$y2,$x2c1,$y2b1)
      #               0    1     2  3    4     5     6     7    8   9   10    11 
       ($c1x,$b1x,$c2x,$b2x) = ($ax1vari[0]-$ax1vari[2],$ax1vari[1]-$ax1vari[3],$ax1vari[4]-$ax1vari[2],$ax1vari[5]-$ax1vari[3]);
       
        my %tem1=('axis' => $ax1,
                  'im_size'=> [@defaul_im_size],
                  'cords'=> [$ax1vari[2],$ax1vari[3],$ax1vari[8],$ax1vari[9]],
                  'bcs' => [-$c1x,-$b1x,-$c2x,-$b2x],
                  'theta'=> $axs_ang{$ax1}[1],
                  'factor' => -1,
                  'width' => $stck_width/2,
                  'steps' => $angle_steps,
                  'tranthe' => $axs_ang{$ax1}[1],#-$shif_angle,
                  'shiftangele'=> $shif_angle);
        #($bs1x1,$bs1y1,$bs1x2,$bs1y2)
        my %bezhash1=xytoradial::get_rib_corda3a1(%tem1);
        
        ($c1x,$b1x,$c2x,$b2x) = ($ax2vari[0]-$ax2vari[2],$ax2vari[1]-$ax2vari[3],$ax2vari[4]-$ax2vari[2],$ax2vari[5]-$ax2vari[3]);
        my %tem2=('axis' => $ax2,
                  'im_size'=> [@defaul_im_size],
                  'cords'=> [$ax2vari[2],$ax2vari[3],$ax2vari[8],$ax2vari[9]],
                  'bcs' => [$c1x,$b1x,$c2x,$b2x],
                  'theta'=> $axs_ang{$ax2}[1],
                  'factor' => 1,
                  'width' => $stck_width/2,
                  'steps' => $angle_steps,
                  'tranthe' => $axs_ang{$ax2}[1],#+$shif_angle,
                  'shiftangele'=> $shif_angle);
        #($bs2x1,$bs2y1,$bs2x2,$bs2y2)
        my %bezhash2=xytoradial::get_rib_corda3a1(%tem2);
        my @points1;
        my @bez_input1=($ax1vari[4],$ax1vari[5]);
        push @bez_input1,@{$bezhash1{'x1y1'}};
        push @bez_input1,@{$bezhash2{'x1y1'}};
        push @bez_input1,($ax2vari[0],$ax2vari[1]);
        #my $bezier1 = Math::Bezier->new($ax1vari[4],$ax1vari[5],$bs1x1,$bs1y1,$bs2x1,$bs2y1,$ax2vari[0],$ax2vari[1]);
        my $bezier1 = Math::Bezier->new(@bez_input1);
        push @points1, ($ax1vari[4],$ax1vari[5]);
        push @points1, $bezier1->curve($no_curve_pt);
        push @points1, ($ax2vari[0],$ax2vari[1]);
        my @bezier_points;
        while (@points1) {
            push @bezier_points, [ splice( @points1, 0, 2 ) ];
        }
        #push @{$clone_ribion{$clonkey}{'rib1'}}, @points1;
        push @bezier_points, [$ax2vari[10],$ax2vari[11]];
        my @points2;
        my @bez_input2=($ax2vari[10],$ax2vari[11]);
        push @bez_input2,@{$bezhash2{'x2y2'}};
        push @bez_input2,@{$bezhash1{'x2y2'}};
        push @bez_input2,($ax1vari[6],$ax1vari[7]);
        #my $bezier2 = Math::Bezier->new($ax2vari[10],$ax2vari[11],$bs2x2,$bs2y2,$bs1x2,$bs1y2,$ax1vari[6],$ax1vari[7]);
        my $bezier2 = Math::Bezier->new(@bez_input2);
        push @points2, ($ax2vari[10],$ax2vari[11]);
        @points2=$bezier2->curve($no_curve_pt);
        push @points2, ($ax1vari[6],$ax1vari[7]);
        my @bezier2_points2;
        while (@points2) {
            push @bezier_points, [ splice( @points2, 0, 2 ) ];
        }
        push @{$clone_ribion{$clonkey}{$ribkey}}, @bezier_points;
        #push @{$clone_ribion{$clonkey}{'rib1'}},( @points1,@points2);
        #print "$ax1\t$tem1{'theta'}[1]\t$ax2\t$tem2{'theta'}[1]\n";
    }else{
      # @points = ($x1c1,$y1b1,$x1,$y1,$x1c2,$y1b2,$x2c2,$y2b2,$x2,$y2,$x2c1,$y2b1)
      #               0    1     2  3    4     5     6     7    8   9   10    11 
       ($c1x,$b1x,$c2x,$b2x) = ($ax1vari[0]-$ax1vari[2],$ax1vari[1]-$ax1vari[3],$ax1vari[4]-$ax1vari[2],$ax1vari[5]-$ax1vari[3]);
        my %tem1=('axis' => $ax1,
                  'im_size'=> [@defaul_im_size],
                  'cords'=> [$ax1vari[2],$ax1vari[3],$ax1vari[8],$ax1vari[9]],
                  'bcs' => [$c1x,$b1x,$c2x,$b2x],
                  'theta'=> $axs_ang{$ax1}[1],
                  'factor' => 1,
                  'width' => $stck_width/2,
                  'steps' => $angle_steps,
                  'tranthe' => $axs_ang{$ax1}[1],#+$shif_angle,
                  'shiftangele'=> $shif_angle);
        
        #($bs1x1,$bs1y1,$bs1x2,$bs1y2)=xytoradial::get_rib_cord(%tem1);
        my %bezhash1=xytoradial::get_rib_cord(%tem1);
        ($c1x,$b1x,$c2x,$b2x) = ($ax2vari[0]-$ax2vari[2],$ax2vari[1]-$ax2vari[3],$ax2vari[4]-$ax2vari[2],$ax2vari[5]-$ax2vari[3]);
        my %tem2=('axis' => $ax2,
                  'im_size'=> [@defaul_im_size],
                  'cords'=> [$ax2vari[2],$ax2vari[3],$ax2vari[8],$ax2vari[9]],
                  'bcs' => [-$c1x,-$b1x,-$c2x,-$b2x],
                  'theta'=> $axs_ang{$ax2}[1],
                  'factor' => -1,
                  'width' => $stck_width/2,
                  'steps' => $angle_steps,
                  'tranthe' => $axs_ang{$ax2}[1],#-$shif_angle,
                  'shiftangele'=> $shif_angle);
        
        #($bs2x1,$bs2y1,$bs2x2,$bs2y2)=xytoradial::get_rib_cord(%tem2);
        my %bezhash2=xytoradial::get_rib_cord(%tem2);
        my @points1;
        my @bez_input1=($ax1vari[0],$ax1vari[1]);
        push @bez_input1,@{$bezhash1{'x1y1'}};
        push @bez_input1,@{$bezhash2{'x1y1'}};
        push @bez_input1,($ax2vari[4],$ax2vari[5]);
        
        my $bezier1 = Math::Bezier->new(@bez_input1);
        #my $bezier1 = Math::Bezier->new($ax1vari[0],$ax1vari[1],$bs1x1,$bs1y1,$bs2x1,$bs2y1,$ax2vari[4],$ax2vari[5]);
        push @points1, ($ax1vari[0],$ax1vari[1]);
        push @points1, $bezier1->curve($no_curve_pt);
        push @points1, ($ax2vari[4],$ax2vari[5]);
        my @bezier_points;
        while (@points1) {
            push @bezier_points, [ splice( @points1, 0, 2 ) ];
        }
        #push @{$clone_ribion{$clonkey}{'rib1'}}, @points1;
        push @bezier_points, [$ax2vari[6],$ax2vari[7]];
        my @points2;
        my @bez_input2=($ax2vari[6],$ax2vari[7]);
        push @bez_input2,@{$bezhash2{'x2y2'}};
        push @bez_input2,@{$bezhash1{'x2y2'}};
        push @bez_input2,($ax1vari[10],$ax1vari[11]);
        
        my $bezier2 = Math::Bezier->new(@bez_input2);
        #my $bezier2 = Math::Bezier->new($ax2vari[6],$ax2vari[7],$bs2x2,$bs2y2,$bs1x2,$bs1y2,$ax1vari[10],$ax1vari[11]);
        
        push @points2, ($ax2vari[6],$ax2vari[7]);
        @points2=$bezier2->curve($no_curve_pt);
        push @points2, ($ax1vari[10],$ax1vari[11]);
        my @bezier2_points2;
        while (@points2) {
            push @bezier_points, [ splice( @points2, 0, 2 ) ];
        }
        push @{$clone_ribion{$clonkey}{$ribkey}}, @bezier_points;
        #push @{$clone_ribion{$clonkey}{'rib1'}},( @points1,@points2);
        #print "$ax1\t$tem1{'theta'}[1]\t$ax2\t$tem2{'theta'}[1]\n";
    }
    }

  }
  
  #print Dumper \@all_ax;
  
}
#print Dumper \%clone_ribion;
my $alpha=40;
my $alphablack=100;
$img ->setThickness(0);
for my $keyclo (keys %clone_ribion){
my %xtemhash = %{$clone_ribion{$keyclo}};
for my $rib (keys %xtemhash){
my @xyz = @{$xtemhash{$rib}};
#my @newtemarray = @{$xtemhash{$rib}};
my @newtemarray;
foreach my $dcop (@xyz){
  push @newtemarray, $dcop;
}
my $bl = $img->colorAllocateAlpha(@{$clone_color_code{$keyclo}},$alpha);


my $bl1 = $img->colorAllocateAlpha(@{$clone_color_code{$keyclo}},$alphablack); 

my $poly = new GD::Polygon;
#print Dumper \@newtemarray;
#while (@newtemarray){
foreach my $txxt (@newtemarray){
my @temxsp = @{$txxt};
#print Dumper \@temxsp;
#my @newploypt = splice(splice(@newtemarray, 0,2), 0,2); 
my @temsplic = splice(@temxsp, 0,2);

my @newploypt = splice(@temsplic, 0,2); 
#print Dumper \@newploypt;
$poly->addPt($newploypt[0], $newploypt[1]);
}
if ($keyclo eq 'BG_H1' | $keyclo eq 'BG_H5'){
  #$img->filledPolygon($poly,$bl1);
}else{
$img->filledPolygon($poly,$bl);
}
}
}
#drawimage::circ_out(($defaul_im_size[0],$defaul_im_size[1],$inner_cir,$outer_cir,$img));
#drawimage::write_img($img);
for my $pkeys (keys %point_data2){
my %tem=%{$point_data2{$pkeys}};
#print Dumper \%tem;
for my $tmkey (keys %tem){
my %temrecod;
my @temval=@{$tem{$tmkey}};
my $nval1=$inner_cir+($temval[0]-$minmax[0])/($minmax[1]-$minmax[0])*($outer_cir-$inner_cir);
my $nval2=$inner_cir+($temval[1]-$minmax[0])/($minmax[1]-$minmax[0])*($outer_cir-$inner_cir);

%temrecod=('im_size'=> [@defaul_im_size],
            'value' => [$nval1,$nval2],
            'theta' => $axs_ang{$tmkey}[1],
            'width' => $stck_width/3
);
#print Dumper \%temrecod;
my @points=(xytoradial::get_stack_cord(%temrecod));
# @points = ($x1c1,$y1b1,$x1,$y1,$x1c2,$y1b2,$x2c2,$y2b2,$x2,$y2,$x2c1,$y2b1)
#               0    1     2  3    4     5     6     7    8   9   10    11 
push @{$clone_stack_cord{$pkeys}{$tmkey}}, @points;
my @test;
my $alpha=0;
my $alphaBG=0;
$img ->setThickness(1);
my $bl = $img->colorAllocateAlpha(@{$clone_color_code{$pkeys}},$alpha);
my $blBG = $img->colorAllocateAlpha(@{$clone_color_code{$pkeys}},$alphaBG);
my $poly = new GD::Polygon;
#my $poly = GD::SVG::Polygon->new;
while (@points) {
  my @newpnt= splice( @points, 0, 2 );
  #print "($newpnt[0], $newpnt[1])\n";
  $poly->addPt($newpnt[0], $newpnt[1]);
}
    if($pkeys eq 'BG_H1' | $pkeys eq 'BG_H5'){
        $img->setThickness(40);
        $img->polygon($poly,$blBG);
    }
}
}

for my $pkeys (keys %point_data2){
my %tem=%{$point_data2{$pkeys}};
#print Dumper \%tem;
for my $tmkey (keys %tem){
my %temrecod;
my @temval=@{$tem{$tmkey}};
my $nval1=$inner_cir+($temval[0]-$minmax[0])/($minmax[1]-$minmax[0])*($outer_cir-$inner_cir);
my $nval2=$inner_cir+($temval[1]-$minmax[0])/($minmax[1]-$minmax[0])*($outer_cir-$inner_cir);

%temrecod=('im_size'=> [@defaul_im_size],
            'value' => [$nval1,$nval2],
            'theta' => $axs_ang{$tmkey}[1],
            'width' => $stck_width/3
);
#print Dumper \%temrecod;
my @points=(xytoradial::get_stack_cord(%temrecod));
# @points = ($x1c1,$y1b1,$x1,$y1,$x1c2,$y1b2,$x2c2,$y2b2,$x2,$y2,$x2c1,$y2b1)
#               0    1     2  3    4     5     6     7    8   9   10    11 
push @{$clone_stack_cord{$pkeys}{$tmkey}}, @points;
my @test;
my $alpha=0;
my $alphaBG=0;
$img ->setThickness(1);
my $bl = $img->colorAllocateAlpha(@{$clone_color_code{$pkeys}},$alpha);
my $blBG = $img->colorAllocateAlpha(@{$clone_color_code{$pkeys}},$alphaBG);
my $poly = new GD::Polygon;
#my $poly = GD::SVG::Polygon->new;
while (@points) {
  my @newpnt= splice( @points, 0, 2 );
  #print "($newpnt[0], $newpnt[1])\n";
  $poly->addPt($newpnt[0], $newpnt[1]);
}
    if($pkeys ne 'BG_H1' & $pkeys ne 'BG_H5'){
        $img ->setThickness(1);
    $img->filledPolygon($poly,$bl);
    }
}
}
#drawimage::circ_out(($defaul_im_size[0],$defaul_im_size[1],$inner_cir,$outer_cir,$img));

my @img_write= ($outfile,$img);
drawimage::write_img(@img_write);
#**********************************************************#
#----------------------------------------------------------#
#------Add ticks and tick values to the plots--------------#

#my $file=$outfile;
my @svgfile;
open my $in, "<:encoding(utf8)", $outfile or die "$outfile: $!";
my $iter=1;
while (my $line = <$in>) {
    push @svgfile,$line;
     $iter=$iter+1;
    
}
close $in;
my @plotinfo = split /[\/]/, $infile;
#print Dumper @plotinfo;
my $infofile = $plotinfo[$#plotinfo] =~ s/.txt/.tsv/r;
$infofile = "maxsum_".$infofile;
$plotinfo[$#plotinfo]=$infofile;
my $readinfo= join"/",@plotinfo;
print "$readinfo\n";
#my @defaul_im_size=(1000,1000);
#my $outer_cir= ($defaul_im_size[0]/2)-($defaul_im_size[0]/2)*0.1;#450; #radius of outer circle upper limit 
#my $inner_cir=($defaul_im_size[0])*0.1;
#print $outer_cir, $inner_cir;
my $myx=50;
my @myy;
#----------Number of intervals-------------##
my $num_ticks=10;
#-------------------------------------------#
my $spce= ($outer_cir-$inner_cir)/$num_ticks;
for (my $ix=0; $ix <= $num_ticks; $ix++){
    push @myy,$inner_cir+($spce*$ix);
}
#print Dumper @myy;
open(my $fhinfo, '<:encoding(UTF-8)', $readinfo)
  or die "Could not open file '$readinfo' $!";
my %infohash;
my @maxval;
while (my $row = <$fhinfo>) {
    chomp $row;
    my @rowarray= split /\s+/,$row;
    push @maxval,$rowarray[$#rowarray];
}
close $fhinfo;
my @mymax = max(@maxval);
print Dumper $mymax[0];


#my @svgkeys = sort(keys %svgfile);
my $des=$outfile;
open(DES,'>',$des) or die $!;
foreach my $iter (0..2){
    print  DES $svgfile[$iter];
}
my $tickint= $mymax[0]/$num_ticks;

my $tickval=0;#$mymax[0];
my $icont=0;
foreach(@myy){
    my $newy=($defaul_im_size[0]/2)-$_;
    
#print  DES qq{\t<line x1="50" x2="450" y1="$newy" y2="$newy" stroke="black" stroke-width= "2"/>\n};
if($icont % 2 == 0 ){
print  DES qq{\t<text x="5150" y="$newy" dy="-10" font-family="arial" font-size="450px" fill="black">$tickval</text>\n};  #font-weight="bold"
}
$tickval=$tickval+$tickint;
$icont++;

}
foreach my $iter (3..$#svgfile){
    print  DES $svgfile[$iter];
}

close $des;
=code
perl petal_plot_v3.pl
##--------------------------------------
  if (@all_ax == 1){
    continue;
  }
  if (@all_ax == 2){
    my ($ax1, $ax2)= @all_ax;
    my @ax1vari = @{$temclone{$ax1}};
    my @ax2vari = @{$temclone{$ax2}};
    my ($bs1x1,$bs1y1,$bs1x2,$bs1y2);
    my ($bs2x1,$bs2y1,$bs2x2,$bs2y2);
    my ($c1x,$b1x,$c2x,$b2x);
    if (($ax1 eq "a1") && ($ax2 ne "a2")){
            # @points = ($x1c1,$y1b1,$x1,$y1,$x1c2,$y1b2,$x2c2,$y2b2,$x2,$y2,$x2c1,$y2b1)
      #               0    1     2  3    4     5     6     7    8   9   10    11 
       ($c1x,$b1x,$c2x,$b2x) = ($ax1vari[0]-$ax1vari[2],$ax1vari[1]-$ax1vari[3],$ax1vari[4]-$ax1vari[2],$ax1vari[5]-$ax1vari[3]);
       
        my %tem1=('axis' => $ax1,
                  'im_size'=> [@defaul_im_size],
                  'cords'=> [$ax1vari[2],$ax1vari[3],$ax1vari[8],$ax1vari[9]],
                  'bcs' => [-$c1x,-$b1x,-$c2x,-$b2x],
                  'theta'=> $axs_ang{$ax1}[1],
                  'factor' => -1,
                  'steps' => $angle_steps,
                  'tranthe' => $axs_ang{$ax1}[1],#-$shif_angle,
                  'shiftangele'=> $shif_angle);
        #($bs1x1,$bs1y1,$bs1x2,$bs1y2)
        my %bezhash1=xytoradial::get_rib_corda3a1(%tem1);
        
        ($c1x,$b1x,$c2x,$b2x) = ($ax2vari[0]-$ax2vari[2],$ax2vari[1]-$ax2vari[3],$ax2vari[4]-$ax2vari[2],$ax2vari[5]-$ax2vari[3]);
        my %tem2=('axis' => $ax2,
                  'im_size'=> [@defaul_im_size],
                  'cords'=> [$ax2vari[2],$ax2vari[3],$ax2vari[8],$ax2vari[9]],
                  'bcs' => [$c1x,$b1x,$c2x,$b2x],
                  'theta'=> $axs_ang{$ax2}[1],
                  'factor' => 1,
                  'steps' => $angle_steps,
                  'tranthe' => $axs_ang{$ax2}[1],#+$shif_angle,
                  'shiftangele'=> $shif_angle);
        #($bs2x1,$bs2y1,$bs2x2,$bs2y2)
        my %bezhash2=xytoradial::get_rib_corda3a1(%tem2);
        my @points1;
        my @bez_input1=($ax1vari[4],$ax1vari[5]);
        push @bez_input1,@{$bezhash1{'x1y1'}};
        push @bez_input1,@{$bezhash2{'x1y1'}};
        push @bez_input1,($ax2vari[0],$ax2vari[1]);
        #my $bezier1 = Math::Bezier->new($ax1vari[4],$ax1vari[5],$bs1x1,$bs1y1,$bs2x1,$bs2y1,$ax2vari[0],$ax2vari[1]);
        my $bezier1 = Math::Bezier->new(@bez_input1);
        push @points1, ($ax1vari[4],$ax1vari[5]);
        push @points1, $bezier1->curve($no_curve_pt);
        push @points1, ($ax2vari[0],$ax2vari[1]);
        my @bezier_points;
        while (@points1) {
            push @bezier_points, [ splice( @points1, 0, 2 ) ];
        }
        #push @{$clone_ribion{$clonkey}{'rib1'}}, @points1;
        push @bezier_points, [$ax2vari[10],$ax2vari[11]];
        my @points2;
        my @bez_input2=($ax2vari[10],$ax2vari[11]);
        push @bez_input2,@{$bezhash2{'x2y2'}};
        push @bez_input2,@{$bezhash1{'x2y2'}};
        push @bez_input2,($ax1vari[6],$ax1vari[7]);
        #my $bezier2 = Math::Bezier->new($ax2vari[10],$ax2vari[11],$bs2x2,$bs2y2,$bs1x2,$bs1y2,$ax1vari[6],$ax1vari[7]);
        my $bezier2 = Math::Bezier->new(@bez_input2);
        push @points2, ($ax2vari[10],$ax2vari[11]);
        @points2=$bezier2->curve($no_curve_pt);
        push @points2, ($ax1vari[6],$ax1vari[7]);
        my @bezier2_points2;
        while (@points2) {
            push @bezier_points, [ splice( @points2, 0, 2 ) ];
        }
        push @{$clone_ribion{$clonkey}{'rib1'}}, @bezier_points;
        #push @{$clone_ribion{$clonkey}{'rib1'}},( @points1,@points2);
        #print "$ax1\t$tem1{'theta'}[1]\t$ax2\t$tem2{'theta'}[1]\n";
    }else{
      # @points = ($x1c1,$y1b1,$x1,$y1,$x1c2,$y1b2,$x2c2,$y2b2,$x2,$y2,$x2c1,$y2b1)
      #               0    1     2  3    4     5     6     7    8   9   10    11 
       ($c1x,$b1x,$c2x,$b2x) = ($ax1vari[0]-$ax1vari[2],$ax1vari[1]-$ax1vari[3],$ax1vari[4]-$ax1vari[2],$ax1vari[5]-$ax1vari[3]);
        my %tem1=('axis' => $ax1,
                  'im_size'=> [@defaul_im_size],
                  'cords'=> [$ax1vari[2],$ax1vari[3],$ax1vari[8],$ax1vari[9]],
                  'bcs' => [$c1x,$b1x,$c2x,$b2x],
                  'theta'=> $axs_ang{$ax1}[1],
                  'factor' => 1,
                  'steps' => $angle_steps,
                  'tranthe' => $axs_ang{$ax1}[1],#+$shif_angle,
                  'shiftangele'=> $shif_angle);
        
        #($bs1x1,$bs1y1,$bs1x2,$bs1y2)=xytoradial::get_rib_cord(%tem1);
        my %bezhash1=xytoradial::get_rib_cord(%tem1);
        ($c1x,$b1x,$c2x,$b2x) = ($ax2vari[0]-$ax2vari[2],$ax2vari[1]-$ax2vari[3],$ax2vari[4]-$ax2vari[2],$ax2vari[5]-$ax2vari[3]);
        my %tem2=('axis' => $ax2,
                  'im_size'=> [@defaul_im_size],
                  'cords'=> [$ax2vari[2],$ax2vari[3],$ax2vari[8],$ax2vari[9]],
                  'bcs' => [-$c1x,-$b1x,-$c2x,-$b2x],
                  'theta'=> $axs_ang{$ax2}[1],
                  'factor' => -1,
                  'steps' => $angle_steps,
                  'tranthe' => $axs_ang{$ax2}[1],#-$shif_angle,
                  'shiftangele'=> $shif_angle);
        
        #($bs2x1,$bs2y1,$bs2x2,$bs2y2)=xytoradial::get_rib_cord(%tem2);
        my %bezhash2=xytoradial::get_rib_cord(%tem2);
        my @points1;
        my @bez_input1=($ax1vari[0],$ax1vari[1]);
        push @bez_input1,@{$bezhash1{'x1y1'}};
        push @bez_input1,@{$bezhash2{'x1y1'}};
        push @bez_input1,($ax2vari[4],$ax2vari[5]);
        
        my $bezier1 = Math::Bezier->new(@bez_input1);
        #my $bezier1 = Math::Bezier->new($ax1vari[0],$ax1vari[1],$bs1x1,$bs1y1,$bs2x1,$bs2y1,$ax2vari[4],$ax2vari[5]);
        push @points1, ($ax1vari[0],$ax1vari[1]);
        push @points1, $bezier1->curve($no_curve_pt);
        push @points1, ($ax2vari[4],$ax2vari[5]);
        my @bezier_points;
        while (@points1) {
            push @bezier_points, [ splice( @points1, 0, 2 ) ];
        }
        #push @{$clone_ribion{$clonkey}{'rib1'}}, @points1;
        push @bezier_points, [$ax2vari[6],$ax2vari[7]];
        my @points2;
        my @bez_input2=($ax2vari[6],$ax2vari[7]);
        push @bez_input2,@{$bezhash2{'x2y2'}};
        push @bez_input2,@{$bezhash1{'x2y2'}};
        push @bez_input2,($ax1vari[10],$ax1vari[11]);
        
        my $bezier2 = Math::Bezier->new(@bez_input2);
        #my $bezier2 = Math::Bezier->new($ax2vari[6],$ax2vari[7],$bs2x2,$bs2y2,$bs1x2,$bs1y2,$ax1vari[10],$ax1vari[11]);
        
        push @points2, ($ax2vari[6],$ax2vari[7]);
        @points2=$bezier2->curve($no_curve_pt);
        push @points2, ($ax1vari[10],$ax1vari[11]);
        my @bezier2_points2;
        while (@points2) {
            push @bezier_points, [ splice( @points2, 0, 2 ) ];
        }
        push @{$clone_ribion{$clonkey}{'rib1'}}, @bezier_points;
        #push @{$clone_ribion{$clonkey}{'rib1'}},( @points1,@points2);
        #print "$ax1\t$tem1{'theta'}[1]\t$ax2\t$tem2{'theta'}[1]\n";
    }
  }

  ##------------------------------


my $filename1 = 'line_data.txt';
open(my $fh1, '<:encoding(UTF-8)', $filename1)
  or die "Could not open file '$filename1' $!";
my %line_hash;
while (my $row = <$fh1>) {
  chomp $row;
  @rowarray= split /\s+/,$row;
  #print Dumper \@rowarray;
  push @{$line_hash{$rowarray[0]}}, @rowarray[1..$#rowarray];
  #$point_data{$rowarray[0]}{$rowarray[1]}= \@rowarray;
  
}
#print Dumper \%line_hash;
my %line_axis_xy;
for my $lkey (keys %line_hash){
  my ($clone,$ax1,$ax2,$val1,$val2,$colr,$crs)=@{$line_hash{$lkey}};
  my %two_ax=($ax1=>$val1,
              $ax2=>$val2);
  my @line_cord;
  foreach my $tkey (keys %two_ax){
    my $nval=$inner_cir+($two_ax{$tkey}-$minmax[0])/($minmax[1]-$minmax[0])*($outer_cir-$inner_cir);
    my %tem1=('im_size'=> [@defaul_im_size],
                  'value'=> $nval,
                  'theta'=> $axs_ang{$tkey});
    
    my ($tempx,$tempy)=xytoradial::get_xy_loc(%tem1);
    push @line_cord, ($tempx,$tempy);
    push @{$line_axis_xy{$lkey}{$tkey}},($defaul_im_size[0],$defaul_im_size[1],$tempx,$tempy,$nval,$axs_ang{$tkey}[1]); 
  }
  my $bezier = Math::Bezier->new(drawimage::draw_curv_lines(%{$line_axis_xy{$lkey}}));
  my @points= $bezier->curve(40);
  my @bezier_points;
  while (@points) {
    push @bezier_points, [ splice( @points, 0, 2 ) ];
  }
  my $bl=$img->colorAllocate(mycolors::findcolor($colr));
  $img ->setThickness(1);
  for my $iz (0 .. @bezier_points-2){
    #print Dumper \$bezier_points[$iz][0];
    $img->line($bezier_points[$iz][0],$bezier_points[$iz][1], $bezier_points[ $iz + 1 ][0],$bezier_points[ $iz + 1 ][1],$bl);
    }
  
}