---

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);

}

---