source: ether_eccad/trunk/API_EXTRACT/src/grid.h @ 68

#ifndef  _GRID_H_
#define _GRID_H_
#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <values.h>

#include <list.h>

// la structure interval de temps
typedef struct {
   char debut[10];
   char fin[10];
}Temporal, *ptrTemporal;

ptrTemporal allocTemporal(char *datdeb, char *datfin);
void freeTemporal(ptrTemporal t);

// la structure parametre
typedef struct {
   int id;
   int unit;
   char *name;
   double  mini;
   double  maxi;
}Param, *ptrParam;

ptrParam allocParam(int id, int unit, char *name, double mini, double maxi);

void freeParam(ptrParam);
int inParamId( int *id, ptrParam p);

// la structure Grid
typedef struct {
   int width;
   int height;
   unsigned int *datas;
   char date[10];
   unsigned int mini;    // valeurs extremes calculées
   unsigned int maxi;
}Grid, *ptrGrid;

typedef struct {
   char date[10];    // AAAAMMJJ //
   double value;     // valeur en flottant d'un pt de grille
}Pt, *ptrPt;

ptrPt allocPt( char *date, double value);
#define freePt(p) free((p))

// les times series
typedef struct {
   
   double lat;        //position
   double lon;
   LIST  pts;
   ptrTemporal dates; //extremes
}Serial, *ptrSerial;

/* structure d'une zone geographique */
typedef struct{
  double latMin;
  double latMax;
  double lonMin;
  double lonMax;
}Zone, *ptrZone;

#define lat2y(zone, lat)\
 (double)(((double)(lat) -(zone)->latMin ) / ((zone)->latMax - (zone)->latMin))
#define lon2x(zone, lon)\
 (double)(((double)(lon) - (zone)->lonMin) / ((zone)->lonMax - (zone)->lonMin))

#define x2lon(zone, x)\
 (zone)->lonMin + (double)x * ((zone)->lonMax - (zone)->lonMin);
#define y2lat(zone, y)\
 (zone)->latMin + (double)y * ((zone)->latMax - (zone)->latMin);

/* conversions 0-1 -> grille */
#define x2i(x, grid) ((int)( (double)(x) * (double)((grid)->width  ) ))
#define y2j(y, grid) ((int)( (double)(y) * (double)((grid)->height ) ))
#define i2x(i, grid) (((double)(i) + 0.5) / (double)((grid)->width  ))
#define j2y(j, grid) (((double)(j) + 0.5) / (double)((grid)->height ))

/* conversions latlon -> grid */
#define lat2j(lat, g, zone) y2j(lat2y((zone), (lat)), (g))
#define lon2i(lon, g, zone) x2i(lon2x((zone), (lon)), (g))
#define i2lon(i, g, zone) x2lon(zone, i2x(i, (g)) )
#define j2lat(j, g, zone) y2lat(zone, j2y(j, (g)) )


ptrGrid allocGrid(int width, int height);
ptrGrid dupGrid(ptrGrid g);
void freeGrid(ptrGrid);
void freeGrids(LIST grids);


//grille géocaliséO
typedef struct {
   ptrParam param;
   ptrGrid grid;
   ptrZone zone;
}GridGL, *ptrGridGL;

ptrGridGL allocGridGL(ptrGrid g, ptrZone z);
void freeGridGL(ptrGridGL);

// convertisseurs codage grille <->flottant
#define int2v(i, min, max)                              \
   (double)(min) +                                      \
    (                                                   \
      ((double)(max)-(double)(min)) /                   \
      (double)((unsigned)MAXINT +(unsigned)MAXINT)      \
    ) * (double)(i)                                     
 

#define v2int(v, min, max)  \
 (unsigned int)( \
    ((double)(v) - (double)(min)) * \
    ((double)((unsigned)MAXINT +(unsigned)MAXINT) /   \
     (double)((double)(max)-(double)(min))) + \
    (double)0.5 \
                 )
// dénition de la zone globale (uniquement enlecture )
extern const Zone globalZone;

// macros de calcul de pas d'une grille en fonction du domaine geographique
#define pasLat(z, g) (double)((z)->latMax -(z)->latMin)/(g)->height;
#define pasLon(z, g) (double)((z)->lonMax -(z)->lonMin)/(g)->width;

#endif //GRID_H_