# include "Ximage.h"

int LoadPicture (int sock) {

  Header header;
  char *buff;
  int status, bytes_left;
  Section *section;
  KapaImageWidget *image;
  Graphic *graphic;
  double Xoffset, Yoffset, wx;

  graphic = GetGraphic ();
  section = GetActiveSection();
  if (section->image == NULL) {
    section->image = InitImageWidget ();
    SetSectionSizes (section);
  }
  image = section->image;
  
  KiiSendMessage (sock, "%d", graphic->Npixels);

  // when we load a new picture, use the same orientation as the old picture
  Xoffset = 0.0;
  Yoffset = 0.0;
  if (image[0].image[0].matrix.size) {
    Xoffset = image[0].picture.X - 0.5*image[0].image[0].matrix.Naxis[0];
    Yoffset = image[0].picture.Y - 0.5*image[0].image[0].matrix.Naxis[1];
  }

  gfits_init_header (&header);
  header.Naxes = 2;
  KiiScanMessage (sock, "%d %d %d %d %lf %lf", &header.Naxis[0], &header.Naxis[1]);
  KiiScanMessage (sock, "%d %d %d %d %lf %lf", &header.bitpix, &header.unsign, &header.bzero, &header.bscale);
  KiiScanMessage (sock, "%lf %lf %s %s",  &image[0].image[0].zero, &image[0].image[0].range, image[0].image[0].name, image[0].image[0].file);
  KiiScanMessage (sock, "%lf %lf %d", &image[0].image[0].min,  &image[0].image[0].max, &header.size);
  KiiScanMessage (sock, "%lf %f %f %f %f", &image[0].image[0].coords.crval1, &image[0].image[0].coords.crpix1, &image[0].image[0].coords.cdelt1, &image[0].image[0].coords.pc1_1, &image[0].image[0].coords.pc1_2);
  KiiScanMessage (sock, "%lf %f %f %f %f", &image[0].image[0].coords.crval2, &image[0].image[0].coords.crpix2, &image[0].image[0].coords.cdelt2, &image[0].image[0].coords.pc2_1, &image[0].image[0].coords.pc2_2);
  KiiScanMessage (sock, "%s", image[0].image[0].coords.ctype);

  gfits_free_matrix (&image[0].image[0].matrix);
  gfits_create_matrix (&header, &image[0].image[0].matrix);

  // reference point for image is the center pixel
  image[0].picture.X = 0.5*header.Naxis[0] + Xoffset;
  image[0].picture.Y = 0.5*header.Naxis[1] + Yoffset;

  // choose expand for wide to guarantee we fit:
  wx = MAX ((header.Naxis[0] / (float) image[0].wide.dx), (header.Naxis[1] / (float) image[0].wide.dy));
  if (wx > 1.0) {
    image[0].wide.expand = -wx;
  } else {
    image[0].wide.expand = 1.0 / wx;
  }    
  image[0].wide.X = 0.5*header.Naxis[0];
  image[0].wide.Y = 0.5*header.Naxis[1];

  fcntl (sock, F_SETFL, O_NONBLOCK);  

  status = 1;
  buff = image[0].image[0].matrix.buffer;
  bytes_left = header.size;
  image[0].image[0].matrix.size = 0;
  while (bytes_left > 0) {
    status = read (sock, buff, bytes_left);
    if (status == 0) {  /* No more pipe */
      fprintf (stderr, "error: pipe closed\n");
      return (FALSE);
    }
    if (status != -1) { /* pipe has data */
      image[0].image[0].matrix.size += status;
      bytes_left -= status;
      buff = (char *)(buff + status);
    }
  }

  fcntl (sock, F_SETFL, !O_NONBLOCK);  

  if (DEBUG) fprintf (stderr, "read %d bytes\n", image[0].image[0].matrix.size);
  /* it it not obvious this condition should kill kii, but ... */
  if (image[0].image[0].matrix.size != header.size) {  
    fprintf (stderr, "error: expected %d bytes, but got only %d\n", header.size, image[0].image[0].matrix.size);
    return (FALSE);
  }

  if (!USE_XWINDOW) return (TRUE);

  SetColorScale (graphic, image);
  Remap (graphic, image);
  CreateWide (graphic, image);
  if (DEBUG) fprintf (stderr, "remapped image\n");
  Refresh ();
  if (DEBUG) fprintf (stderr, "refreshed\n");
  XFlush (graphic[0].display);

  return (TRUE);
}
