<?php

  // 呼び出され方の想定:
  //   basecolortrack_graphic.php?species=medaka&revision=test&target=scaffold1&start=0&end=8000&width=800

function draw_graph($file, $line1){
  $max = 100;
  $min = 0;
  $req = pre_get_request($_REQUEST);

  if($req['start'] < $req['end']){
    $start = $req['start'];
    $end = $req['end'];
  }else {
    $start = $req['end'];
    $end = $req['start'];
  }

  $basewidth = $end - $start +1;

  // イメージの生成
  $height = 40;
  $image = imagecreatetruecolor( $req['width'], $height );
  $color = imagecolorallocate( $image, 0x00, 0x00, 0x00 );
  $white = imagecolorallocate( $image, 0xFF, 0xFF, 0xFF );
  $black = imagecolorallocate( $image, 0x00, 0x00, 0x00 );  


  // color

  if(($color_pos = strpos($line1,'color=')) !== FALSE){
     list($color_r, $color_g, $color_b, $remain) = split("[ ,\t]", substr($line1, $color_pos + 6), 4);
     $color = imagecolorallocate( $image, $color_r, $color_g, $color_b);
  }


#  imagefill( $image, 0, 0, $white);
  imagefilledrectangle($image, 0, 0, $req['width']-1,$height-1, $white);

  //min
  if(($min_pos = strpos($line1,'min=')) !== FALSE){
     list($min, $remain) = split("[ ,\t]", substr($line1, $min_pos + 4), 2);
  }

  //max
  if(($max_pos = strpos($line1,'max=')) !== FALSE){
     list($max, $remain) = split("[ ,\t]", substr($line1, $max_pos + 4), 2);
  }

  // イメージの作成

  while (!feof ($file)) {
    $line = trim(fgets($file));
    if(strncmp($line, "graph", 5) != 0) continue;

//  for($num=0; $num < $num_genes; $num++){
//    $line = $genes[$num];

    $target_line = substr(( strchr ($line, 'target=') ), 7);
    list($target, $remain) = split("[ \t]", $target_line);
    if(strcmp($target, $req['target']) != 0) continue;

    $range_line = substr(strchr ($line, 'range='), 6);
    list($range_start, $range_end, $remain) = split("[, \t]", $range_line);

    if($range_start > $range_end) {
      $tmp = $range_start;
      $range_start = $range_end;
      $range_end = $tmp;
    }
    if(($range_start > $end) || ($range_end < $start)) continue; 
    
    $unit_line = substr(strchr ($line, 'unit='), 5);
    list($unit, $remain) = split("[, \t]", $unit_line);

    if(($nums_pos = strpos ($line, 'nums=')) === FALSE) return;
    $nums = substr($line, $nums_pos + 5);

    for($i=$range_start;$i<=$range_end; $i+=$unit){
      if($i > $end){break;}
      list($num, $remain) = split("[, \t]", $nums, 2);
      $nums = $remain;
      if($i + $unit <= $start){continue;}
      if($req['start'] == $start){
        $startpos = $req['width'] * ($i - $req['start']) / $basewidth;
        $endpos = $req['width'] * ($i + $unit - $req['start']) / $basewidth -1;
      }else {
        $endpos = $req['width'] * ($req['start'] - $i) / $basewidth;
        $startpos = $req['width'] * ($req['start'] - $i - $unit) / $basewidth -1;
      }
      //echo($num);
      $h = floor($height * ( $max - $num) / ($max - $min));
      imagefilledrectangle($image, $startpos, $h, $endpos, $height, $color);
    }
  }

  // イメージの出力
  header( 'Content-type: image/png' );
  imagepng( $image );

  imagedestroy( $image );
  fclose( $file );
}
?>