int gridread(char *file, void ***data, int datatype, int *nx, int *ny,
float *dx, float *dy, double bndbox[4], double *csize,
float *nodata, int *filetype)
{
FILE *fp;
int channel1,type,iesri;
int i, j, hdrlines = 0;
float value;
char fline[MAXLN], keyword[21], utmetag, utmntag;
float **farr;
int **iarr;
short **sarr;
double adjbndbox[4],utme,utmn;
CELLTYPE *rowbuf;
float floatNull;
if(iesri = GridIOSetup() >= 0)
/*
Open input cell layer with no evaluation (5)
*/
if ((channel1 =
CellLayerOpen(file,READONLY,ROWIO,
&type,csize)) >= 0 )
{
/* The success of these means the file could be opened as an ESRI grid file */
printf("File could be opened as an ESRI grid file\n");
*filetype=1;
/* File types are
0=ASCII
1= ARCVIEW grid via the ESRI Application Programmers Interface */
/* here is the ESRI grid file read stuff - following copyrow example */
*dx = *dy = (float) *csize;
if(type == CELLINT && datatype == RPFLTDTYPE)
{
printf("%s File contains integer data but was read as float\n",file);
}
if(type == CELLFLOAT && datatype != RPFLTDTYPE)
{
printf("%s File contains Float data but was read as Integer\n",file);
}
/*
Get the bounding box of the input cell layer (7)
*/
if (BndCellRead (file, bndbox) < 0)
{
printf ("ERROR: Could not read bounding box of input grid\n");
CellLyrClose(channel1);
GridIOExit();
return(1);
}
/* printf("%f %f %f %f %g\n",bndbox[0],bndbox[1],bndbox[2],bndbox[3],*nodata) */;
/* Bounding box is xllcorner, yllcorner, xurcorner, yurcorner */
/*
Set the Window to the output bounding box and cellsize (9)
*/
if (AccessWindowSet (bndbox, *csize, adjbndbox) < 0)
{
printf ("ERROR: Could not set Window\n");
CellLyrClose(channel1);
GridIOExit();
return(2);
}
/*
Get the number of rows and columns (10)
in the window
*/
*nx = WindowCols();
*ny = WindowRows();
/* Allocate memory and set all type pointers */
*data = matalloc(*nx, *ny, datatype);
farr = (float **) (*data);
iarr = (int **) (*data);
sarr = (short **) (*data);
GetMissingFloat(&floatNull);
*nodata = (datatype == RPFLTDTYPE) ? floatNull: -9999.;
/*
Allocate row buffer (11)
*/
if ((rowbuf = (CELLTYPE *)
CAllocate1 ( *nx, sizeof(CELLTYPE)))
== NULL )
{
printf ("ERROR: Could not allocate memory\n");
CellLyrClose(channel1);
GridIOExit();
return(3);
}
/*
Now copy row into array (12)
*/
for (i=0; i < *ny; i++)
{
GetWindowRow (channel1, i, rowbuf);
if(type == CELLINT)
{
register int *buf = (int *)rowbuf;
if(datatype == RPSHRDTYPE)
{
for(j=0; j < *nx; j++)
{
sarr[j][i]=(buf[j] == MISSINGINT) ? (short) *nodata : (short) buf[j];
}
}
else if(datatype == RPINTDTYPE)
{
for(j=0; j < *nx; j++)
{
iarr[j][i]=(buf[j] == MISSINGINT) ? (int) *nodata : buf[j];
}
}
else
{
for(j=0; j < *nx; j++)
{
farr[j][i]=(buf[j] == MISSINGINT) ? *nodata: (float) buf[j];
}
}
}
else
{
register float *buf = (float *)rowbuf;
/* This is all repeated to get right the typecasting of data from ESRI file into the format
we want */
if(datatype == RPSHRDTYPE)
{
for(j=0; j < *nx; j++)
{
sarr[j][i]=(buf[j] == floatNull) ? (short) *nodata : (short) buf[j];
}
}
else if(datatype == RPINTDTYPE)
{
for(j=0; j < *nx; j++)
{
iarr[j][i]=(buf[j] == floatNull) ? (int) *nodata : (int) buf[j];
}
}
else
{
for(j=0; j < *nx; j++)
{
farr[j][i]= buf[j];
}
}
}
}
/*
Free row buffer (13)
*/
CFree1 ((char *)rowbuf);
/*
Close handles (14)
*/
CellLyrClose(channel1);
/*
Done with the library (15)
*/
GridIOExit();
return(0);
}
/* Here assume file is ASCII Close ESRI stuff. */
CellLyrClose(channel1);
GridIOExit();
*filetype=0;
fp = fopen(file,"r");
if(fp == NULL)
{
printf("\nERROR: Cannot open input file %s\n\n",file);
return(1);
}
/* read ARC-Info header */
while(1)
{
readline(fp, fline);
if(!isalpha(*fline) || *fline == '-')
break;
hdrlines++;
sscanf(fline,"%s %f",keyword,&value);
if(strcmp(keyword,"ncols") == 0 || strcmp(keyword,"NCOLS") == 0)
*nx = (int)value;
else if(strcmp(keyword,"nrows") == 0 || strcmp(keyword,"NROWS") == 0)
*ny = (int)value;
else if(strcmp(keyword,"xllcenter") == 0 || strcmp(keyword,"XLLCENTER") == 0)
{
utmetag = 'c';
utme = value;
}
else if(strcmp(keyword,"xllcorner") == 0 || strcmp(keyword,"XLLCORNER") == 0)
{
utmetag = 'e';
utme = value;
}
else if(strcmp(keyword,"yllcenter") == 0 || strcmp(keyword,"YLLCENTER") == 0)
{
utmntag = 'c';
utmn = value;
}
else if(strcmp(keyword,"yllcorner") == 0 || strcmp(keyword,"YLLCORNER") == 0)
{
utmntag = 'e';
utmn = value;
}
else if(strcmp(keyword,"cellsize") == 0 || strcmp(keyword,"CELLSIZE") == 0)
{
*dx = *dy = value;
*csize = (double) value;
}
else if(strcmp(keyword,"nodata_value") == 0 || strcmp(keyword,"NODATA_VALUE") == 0 ||
strcmp(keyword,"NODATA_value") == 0)
*nodata = value;
}
/* adjust utme and utmn if necessary (we store center of reference cell) */
if(utmetag == 'e') utme = utme + *dx/2;
if(utmntag == 'e') utmn = utmn + *dy/2;
bndbox[0] = utme - *csize/2.;
bndbox[1] = utmn - *csize/2.;
bndbox[2] = bndbox[0] + *csize * (*nx);
bndbox[3] = bndbox[1] + *csize * (*ny);
/* position file pointer for ARC-Info file to beginning of image data */
rewind(fp);
for(i=0; i<hdrlines; i++) readline(fp, fline);
/* convert depending on datatype */
if(datatype == RPSHRDTYPE)
{
sarr = (short **) matalloc(*nx, *ny, datatype);
/* read in the ARC-Info file */
for(i=0; i< *ny; i++)
{
for(j=0; j< *nx; j++)
fscanf(fp,"%hd",&sarr[j][i]);
}
*data = (void **) sarr;
}
else if(datatype == RPINTDTYPE)
{
iarr = (int **) matalloc(*nx, *ny, datatype);
for(i=0; i< *ny; i++)
{
for(j=0; j< *nx; j++)
fscanf(fp,"%d",&iarr[j][i]);
}
*data = (void **) iarr;
}
else if(datatype == RPFLTDTYPE)
{
farr = (float **) matalloc(*nx, *ny, datatype);
/* read in the ARC-Info file */
for(i=0; i< *ny; i++)
{
for(j=0; j< *nx; j++)
{
fscanf(fp,"%f",&farr[j][i]);
}
}
*data = (void **) farr;
}
else
{
printf("\nERROR: unknown datatype (%s).\n\n",datatype);
}
fclose(fp);
return(0);
}
CLIPPER NWROW( PARAMS ) // hWnd, hDC, nGraphRow, hFont, nRowHeight -> nTextRow
{
_retni( GetWindowRow( ( HWND ) _parnl( 1 ), ( HDC ) _parnl( 2 ),
_parni( 3 ), ( HFONT ) _parnl( 4 ), ( WORD ) _parni( 5 ) ) );
}
