---

disk2img.c

Go to the documentation of this file.

#include        "sadie.h"

/*-Copyright Information------------------------------------------------------*/
/* Copyright (c) 1988 by the University of Arizona Digital Image Analysis Lab */
/*-General Information--------------------------------------------------------*/
/*                                                                            */
/*   This function reads a SADIE image from disk into main memory.            */
/*                                                                            */
/*----------------------------------------------------------------------------*/
/*-Interface Information------------------------------------------------------*/
void
DISK2IMG (char *name,           /*  I   String, containing the name of the disk file.         */
          long size,            /*  I   The size of the image file.                                    */
          IMAGE ** out          /*  O   Address of a pointer to the image.                    */
/*----------------------------------------------------------------------------*/
  )
{
  static char array[2] = { 0x00, 0x01 };
  char msg[SLEN];
  short nbnd, nlin, npix, nbit;
  long nbytes;
  FILE *fp = NULL;
  long calcsize;
  long readbytes;

  if (TIMES)
    TIMING (T_DISK2IMG);

  /* open image file */
  if (!(fp = fopen (name, "r")))
    {
      sprintf (msg, " Can't open file %s.", name);
      MESSAGE ('E', msg);
      goto the_end;
    }

  /* read image size parameters from file */
  if (fread (&nbnd, sizeof (nbnd), 1, fp) != 1
      || fread (&nlin, sizeof (nlin), 1, fp) != 1
      || fread (&npix, sizeof (npix), 1, fp) != 1
      || fread (&nbit, sizeof (nbit), 1, fp) != 1)
    {
      sprintf (msg, " Can't read image header from file %s.", name);
      MESSAGE ('E', msg);
      goto the_end;
    }

  /* swap bytes, if necessary */
  if (*(short *) array != 1)
    {
/*        //printf("DISK2IMG: Swapping bytes...\n"); */

      nbnd = SWAP ((unsigned char *) &nbnd);
      nlin = SWAP ((unsigned char *) &nlin);
      npix = SWAP ((unsigned char *) &npix);
      nbit = SWAP ((unsigned char *) &nbit);
    }

  /* Make sure the size matches the number of bytes in the file */
  if (size > 0)
    {
      calcsize =
        512L +
        lrnd ((double) nbnd * (double) nlin * (double) npix *
              ((double) nbit / 8.0));
      if (size < calcsize || calcsize <= 512L)
        {
          printf ("size = %d, calcsize = %d\n", size, calcsize);
          printf ("nbnd = %d, nlin = %d, npix = %d, nbit = %d\n", nbnd, nlin,
                  npix, nbit);
          sprintf (msg, "File %s is not in SADIE format!", name);
          MESSAGE ('E', msg);
          goto the_end;
        }
    }

  /* create image of appropriate size */
  if (!CHECKIMG (*out))
    GETMEM (nbnd, nlin, npix, out);
  if (!CHECKIMG (*out))
    {
      sprintf (msg, "Can't allocate memory for image.");
      MESSAGE ('E', msg);
      goto the_end;
    }

  /* read the rest of the header data from file */
  (*out)->nbit = nbit;
  nbytes =
    (long) (sizeof ((*out)->gmin) + sizeof ((*out)->gmax) +
            sizeof ((*out)->pstr) + sizeof ((*out)->text));
  if (fread (&((*out)->gmin), nbytes, 1, fp) != 1)
    {
      sprintf (msg, " Can't read image header from file %s.", name);
      MESSAGE ('E', msg);
      goto the_end;
    }

  /* check number of bits per pixel */
  if ((*out)->nbit != sizeof (PIXEL) * SBPB
      && (sizeof (PIXEL) * SBPB < (*out)->nbit
          || (*out)->nbit > sizeof (long) * SBPB))
    {
      sprintf (msg, " Can't unpack from %d bits/pixel.", (*out)->nbit);
      MESSAGE ('E', msg);
      goto the_end;
    }

  /* read the image data from file */
  nbytes =
    (long)
    ceil (((double) nbnd * (double) nlin * (double) npix * (*out)->nbit) /
          (double) SBPB);
  readbytes = fread ((*out)->data[0][0], 1, nbytes, fp);
  if (readbytes < nbytes)
    {
      sprintf (msg, " Can't read image data from file %s.", name);
      MESSAGE ('E', msg);
      sprintf (msg, " Read %ld bytes instead of %ld bytes.", readbytes,
               (long) nbytes);
      MESSAGE ('E', msg);
    }

  if ((*out)->nbit < sizeof (PIXEL) * SBPB)
    {
      UNPACK ((*out)->nbit, (long) nbnd * (long) nlin * (long) npix,
              (*out)->data[0][0]);
    }

  /* Set image name to match file name */
  sprintf ((*out)->text, "%s", name);

  if (NAMES)
    {
      MESSAGE ('I', "");
      MESSAGE ('I', "DISK2IMG");
      MESSAGE ('I', "");
      sprintf (msg, " Input file:                                 %s", name);
      MESSAGE ('I', msg);
      MESSAGE ('I', " Input format:                               SADIE");
      sprintf (msg, " Header size:                                %d",
               (int) (sizeof ((*out)->nbnd) + sizeof ((*out)->nlin) +
                      sizeof ((*out)->npix) + sizeof ((*out)->nbit) +
                      sizeof ((*out)->gmin) + sizeof ((*out)->gmax) +
                      sizeof ((*out)->pstr) + sizeof ((*out)->text)));
      MESSAGE ('I', msg);
      sprintf (msg, " Number of bands, lines, pixels:             %d, %d, %d",
               (*out)->nbnd, (*out)->nlin, (*out)->npix);
      MESSAGE ('I', msg);
      sprintf (msg, " Number of bits/pixel/band:                  %d",
               (*out)->nbit);
      MESSAGE ('I', msg);
      sprintf (msg,
               " Graylevel minimum, maximum:                 %10.4e, %10.4e",
               (*out)->gmin, (*out)->gmax);
      MESSAGE ('I', msg);
      MESSAGE ('I', " ...............");
    }

  /* close image file */

the_end:
  if (fp)
    {
      if (fclose (fp))
        {
          sprintf (msg, " Can't close file %s.", name);
          MESSAGE ('W', msg);
        }
    }

  if (TIMES)
    TIMING (T_EXIT);

}

---