]> www.wagner.pp.ru Git - oss/fgis.git/blob - lib/dgt_input.c
Removed erroneously checked in binaries
[oss/fgis.git] / lib / dgt_input.c
1 #include <stdio.h>
2 #include <stdlib.h>
3 #include "dgt.h"
4 int read_item(DGT* dgt)
5 /* read next item of dgt file opened for reading */
6 { int to_skip=0;
7 #ifndef LSB_FIRST
8   DGT_ITEM tmp;
9 #endif   
10 #ifdef LSB_FIRST 
11   if(fread(dgt->buffer,1,14,(FILE *)dgt->F)!=14) return 1;
12 #else
13   if(fread(&tmp,1,14,(FILE *)(dgt->F))!=14) return 1;
14   dgt->buffer->ID=swaplong(tmp.ID);
15   swab(&(tmp.x1),&(dgt->buffer->x1),10);
16 #endif
17   if (dgt->buffer->npoints)
18   { if(dgt->buffer->npoints>MAX_LINE_LEN) {to_skip=(dgt->buffer->npoints-MAX_LINE_LEN)*4;
19                                   dgt->buffer->npoints=MAX_LINE_LEN;}
20 #ifdef LSB_FIRST
21    if (fread(&(dgt->buffer->s),4,dgt->buffer->npoints,
22                    (FILE *)dgt->F)!=dgt->buffer->npoints) 
23    { map_error = ME_READ_ERROR; return 1;}
24 #else       
25    if(fread(&(tmp.s),4,dgt->buffer->npoints,
26     (FILE *)dgt->F)!=dgt->buffer->npoints)
27      {map_error=ME_READ_ERROR;return 1;}
28    swab(&(tmp.s),&(dgt->buffer->s),dgt->buffer->npoints*4);
29 #endif
30    if (to_skip) fseek((FILE*)dgt->F,to_skip,SEEK_CUR);
31   }
32    dgt->item_no++;
33    return dgt->eof=(dgt->buffer->ID==0&&dgt->buffer->npoints==0);
34     
35 }
36 void go_first(DGT *dgt)
37 {
38  fseek((FILE *)dgt->F,128L,SEEK_SET);
39  dgt->item_no=0;
40  dgt->eof=dgt->next_item(dgt);
41 }
42
43 DGT* open_dgt(char *filename)
44 /*opens existing file, given by name*/
45 {FILE *f;
46  if (!(f=fopen(filename,"r")))
47   return NULL;
48  return fopen_dgt(f); 
49 }
50 typedef void (*close_func)(void *f); 
51 DGT* fopen_dgt(FILE *f)
52 /* open existing file given by file variable */
53 { DGT* tmp;
54   DGTHEADER h;
55   if(!(tmp=malloc(sizeof(DGT)))) {map_error=ME_OUT_OF_MEMORY;return NULL;}
56   if(fread(&h,1,sizeof(h),f)!=sizeof(h)) {map_error=ME_READ_ERROR;free(tmp);return NULL;}
57   if(!(tmp->buffer=malloc(sizeof(DGT_ITEM)))) {map_error=ME_OUT_OF_MEMORY;free(tmp);return NULL;}
58   tmp->XLeft=swapdouble(h.xleft);
59   tmp->XRight=swapdouble(h.xright);
60   tmp->YTop=swapdouble(h.ytop);
61   tmp->YBottom=swapdouble(h.ybottom);
62   tmp->projection=h.coord_sys;
63   tmp->mode=MAP_INPUT;
64   tmp->item_no=0;
65   tmp->next_item=read_item;
66   tmp->rewind=go_first;
67   tmp->F=f;
68   tmp->dispose=(close_func)fclose;
69   tmp->cmp_item=NULL; 
70   tmp->next_item(tmp);
71   return tmp;
72 }
73 void close_dgt(DGT* dgt)
74 /* Destroys DGT object */
75 { dgt->rewind(dgt);
76   dgt->dispose(dgt->F);
77   free(dgt->buffer);
78   free(dgt);
79 }
80 DGT_ITEM* dgt_get(DGT* dgt)
81 /* allocates current item in memory and returns pointer to it*/
82 { DGT_ITEM* item;
83   if (!(dgt->mode&MAP_INPUT)) {map_error=ME_INVALID_MODE;return NULL;}
84   if (!dgt->eof) {item= alloc_item(dgt->buffer);dgt->next_item(dgt);return item;}
85   else return NULL;
86 }
87 int dgt_x(DGT* dgt,double x)
88 {
89 return ((x-dgt->XLeft)/(dgt->XRight-dgt->XLeft)*65535)-32767;
90 }
91 int dgt_y(DGT* dgt,double y)
92 {
93 return ((y-dgt->YBottom)/(dgt->YTop-dgt->YBottom)*65535)-32767;
94 }
95 double real_x(DGT* dgt,int x)
96 { return (x+32767)/65535.0*(dgt->XRight-dgt->XLeft)+dgt->XLeft;
97 }
98 double real_y(DGT* dgt,int y)
99 {
100 return (y+32767)/65535.0*(dgt->YTop-dgt->YBottom)+dgt->YBottom;
101 }
102
103 DGT_ITEM* alloc_item(DGT_ITEM *item)
104 {DGT_ITEM* new_item;
105  new_item=malloc(14+item->npoints*4);
106  return memcpy(new_item,item,14+item->npoints*4);
107 }
108 void reset_dgt(DGT* dgt)
109 /* rewinds dgt file for first item, reopens for reading,if nessecary */
110 { dgt->rewind(dgt);
111   dgt->eof=dgt->next_item(dgt);
112 }