1 | #include "eccad.h" |
---|
2 | |
---|
3 | |
---|
4 | EccadDB connectEccadDB(const char* server, const char* dbName, |
---|
5 | const char* user, const char* pass, const char* port) |
---|
6 | { |
---|
7 | EccadDB conn = PQsetdbLogin(server, port, NULL, NULL, dbName, user, pass); |
---|
8 | if(!checkConnectionStatus(conn)) |
---|
9 | { |
---|
10 | if (PQstatus(conn) == CONNECTION_BAD) |
---|
11 | { |
---|
12 | fprintf(stderr, "Connection to database '%s' failed.\n", dbName); |
---|
13 | fprintf(stderr, "%s", PQerrorMessage(conn)); |
---|
14 | PQfinish(conn); |
---|
15 | } |
---|
16 | return NULL; |
---|
17 | } |
---|
18 | return conn; |
---|
19 | } |
---|
20 | |
---|
21 | void closeEccadDB(EccadDB conn) |
---|
22 | { |
---|
23 | PQfinish(conn); |
---|
24 | } |
---|
25 | |
---|
26 | gdImagePtr grid2img(EccadDB db, ptrGrid g, ptrParam p, int idMap) |
---|
27 | { |
---|
28 | gdImagePtr im=NULL; |
---|
29 | register int i,j; |
---|
30 | LIST map; |
---|
31 | /*if( idMap == 48) |
---|
32 | // colormap dynamique les valeurs devrait etre codï¿œs en [O-1] |
---|
33 | map= mapOfParamById(db, idMap, (double)0., (double)1.); |
---|
34 | // map= mapOfParamById(db, idMap, (double)0., (double)683.92); |
---|
35 | else */ |
---|
36 | map= mapOfParamById(db, idMap, p->mini, p->maxi); |
---|
37 | fprintf(stderr,"start colorAllocate with %d colors\n", map->count); |
---|
38 | if(map) |
---|
39 | { |
---|
40 | if (im= gdImageCreate(g->width, g->height)) |
---|
41 | { |
---|
42 | fprintf(stderr," width %d height %d\n",g->width, g->height ); |
---|
43 | colorAllocate(im, map); |
---|
44 | for (j=0 ; j < im->sy; j++) |
---|
45 | for (i=0 ; i < im->sx ; i++) |
---|
46 | { |
---|
47 | int value= val2Color(g->datas[j * im->sx +i], map); |
---|
48 | gdImageSetPixel(im, i, j, value); |
---|
49 | } |
---|
50 | } |
---|
51 | } |
---|
52 | fprintf(stderr," ok\n" ); |
---|
53 | return im; |
---|
54 | } |
---|
55 | |
---|
56 | // insertion d'une colormap dans une image gd |
---|
57 | // aucune couleur ne doit etre alloue dans l'image d'origine |
---|
58 | void colorAllocate(gdImagePtr im, LIST map) |
---|
59 | { |
---|
60 | ptrInterval v; |
---|
61 | int i, j=0; |
---|
62 | loop_through_list(map, v, ptrInterval) |
---|
63 | { |
---|
64 | i=gdImageColorAllocate(im, v->color->red ,v->color->green, v->color->blue); |
---|
65 | //fprintf(stderr, "color Allocatealloue a l-ndex %d :%d \n", j, i); |
---|
66 | j++; |
---|
67 | } |
---|
68 | } |
---|
69 | |
---|
70 | char * oid2filename(EccadDB db, int oidGrille, int width, int height, |
---|
71 | char * idParam, int idMap, char *resultFile) |
---|
72 | { |
---|
73 | char *fileName=NULL; |
---|
74 | int grilleId; |
---|
75 | char query[512]; |
---|
76 | |
---|
77 | PQexec(db,("BEGIN")); |
---|
78 | if( (grilleId = lo_open(db, oidGrille, INV_READ)) != -1) |
---|
79 | { |
---|
80 | ptrGrid g= allocGrid( width, height); |
---|
81 | if(g) |
---|
82 | { |
---|
83 | unsigned long K; |
---|
84 | unsigned long size; |
---|
85 | |
---|
86 | size= (long)g->width * (long)g->height * sizeof(int); |
---|
87 | //fprintf(stderr, "MAXINT %ld 2*MAXINT %ld\n", MAXINT, (unsigned)MAXINT + MAXINT -1); |
---|
88 | if( (K= lo_read(db, grilleId, (char*)g->datas, size)) == size) |
---|
89 | { |
---|
90 | PGresult *res; |
---|
91 | sprintf(query, |
---|
92 | "SELECT min_param,max_param FROM parametre WHERE id_param='%s'", |
---|
93 | idParam); |
---|
94 | res=PQexec(db,query); |
---|
95 | if(checkSQLErrors(res)) |
---|
96 | { |
---|
97 | double mini= atof( PQgetvalue(res, 0, 0)); |
---|
98 | double maxi= atof( PQgetvalue(res, 0, 1)); |
---|
99 | LIST map; |
---|
100 | |
---|
101 | if( map = mapOfParamById(db, idMap, mini, maxi)) |
---|
102 | { |
---|
103 | gdImagePtr im; |
---|
104 | if (im= gdImageCreate(width, height)) |
---|
105 | { |
---|
106 | register int i,j; |
---|
107 | FILE *out; |
---|
108 | fprintf(stderr," oidGrille %d mini %f maxi %f\n", oidGrille, mini,maxi); |
---|
109 | colorAllocate(im, map); |
---|
110 | for (j=0 ; j < im->sy; j++) |
---|
111 | for (i=0 ; i < im->sx ; i++) |
---|
112 | { |
---|
113 | int value= val2Color(g->datas[j * im->sx +i], map); |
---|
114 | // if( g->datas[j * im->sx +i] !=0) printf( " %d\n", g->datas[j * im->sx +i]); |
---|
115 | gdImageSetPixel(im, i, j, value); |
---|
116 | } |
---|
117 | // sprintf(query, "/tmp/%d.%d.png", oidGrille, idMap); |
---|
118 | if( out= fopen(resultFile, "w")) |
---|
119 | { |
---|
120 | gdImagePng(im, out); |
---|
121 | fclose(out); |
---|
122 | fileName=resultFile; |
---|
123 | } |
---|
124 | else fprintf(stderr,"Erreur open File %s\n", fileName); |
---|
125 | gdImageDestroy (im); |
---|
126 | } |
---|
127 | else fprintf(stderr,"Erreur gdImageCreate\n"); |
---|
128 | freeMap(map); |
---|
129 | } |
---|
130 | else fprintf(stderr,"Erreur Check Map idMap %d\n", idMap); |
---|
131 | } |
---|
132 | else fprintf(stderr,"Erreur Check min max from idParam %s\n", idParam); |
---|
133 | PQclear(res); |
---|
134 | } |
---|
135 | else fprintf(stderr," Erreur Read OID : %ld sur %d octets\n", K, size); |
---|
136 | fprintf (stderr, " now free grid\n"); |
---|
137 | freeGrid(g); |
---|
138 | } |
---|
139 | lo_close(db, grilleId); |
---|
140 | } |
---|
141 | else fprintf(stderr," Erreur Open OID : %d\n", oidGrille); |
---|
142 | PQexec(db,("END")); |
---|
143 | return fileName; |
---|
144 | } |
---|
145 | |
---|
146 | char * getRaster(const char *dbHost, const char * dbName, const char *dbUser, int idGrille, int idMap) |
---|
147 | { |
---|
148 | char *resultFile; |
---|
149 | struct stat sts; |
---|
150 | char query[256]; |
---|
151 | |
---|
152 | sprintf(query,"/tmp/%d.%d.png", idGrille, idMap); |
---|
153 | resultFile= strdup(query); |
---|
154 | |
---|
155 | if ( stat (query, &sts) != 0) |
---|
156 | { |
---|
157 | char *fileName=NULL; |
---|
158 | EccadDB db; |
---|
159 | |
---|
160 | if((db= PQsetdbLogin(dbHost,"5432", NULL, NULL, dbName, dbUser, NULL)) && |
---|
161 | (PQstatus(db) != CONNECTION_BAD) ) |
---|
162 | |
---|
163 | { |
---|
164 | PGresult *res; |
---|
165 | sprintf(query, |
---|
166 | "SELECT valeurs_grille, id_param, ncol_grille, nlign_grille FROM grille WHERE id_grille = %d;", |
---|
167 | idGrille); |
---|
168 | res=PQexec(db,query); |
---|
169 | if(checkSQLErrors(res) && PQntuples(res)) |
---|
170 | { |
---|
171 | int oid= atoi(PQgetvalue(res, 0, 0)); |
---|
172 | int width= atoi(PQgetvalue(res, 0, 2)); |
---|
173 | int height= atoi(PQgetvalue(res, 0, 3)); |
---|
174 | fileName= oid2filename( db, oid, width, height, |
---|
175 | PQgetvalue(res, 0, 1), idMap, resultFile); |
---|
176 | } |
---|
177 | else fprintf(stderr, "No found id-grille %d\n", idGrille); |
---|
178 | PQclear(res); |
---|
179 | PQfinish(db); |
---|
180 | } |
---|
181 | else fprintf(stderr, "Erreur connect to db %s on host %s\n", |
---|
182 | dbName, dbHost); |
---|
183 | resultFile= fileName; |
---|
184 | } |
---|
185 | return resultFile; |
---|
186 | } |
---|
187 | // ============================ Fonctions ================================== |
---|
188 | |
---|
189 | |
---|
190 | int checkConnectionStatus(EccadDB conn) |
---|
191 | { |
---|
192 | char* erreur; |
---|
193 | |
---|
194 | if (conn==NULL || PQstatus(conn) == CONNECTION_BAD) |
---|
195 | { |
---|
196 | return (int)NULL; |
---|
197 | } |
---|
198 | else{ |
---|
199 | // la connection est ok |
---|
200 | // mise en place de la date en UTC/GMT |
---|
201 | // db->clock_consumed=clock(); |
---|
202 | PGresult* res=PQexec(conn, "set time zone 'UTC';"); |
---|
203 | if (!(res) || |
---|
204 | ((PQresultStatus(res) != PGRES_TUPLES_OK ) && |
---|
205 | (PQresultStatus(res) != PGRES_COMMAND_OK) ) ){ |
---|
206 | PQclear(res); |
---|
207 | return 0==1; |
---|
208 | } |
---|
209 | res= PQexec(conn, "set enable_seqscan='false';"); |
---|
210 | if (!(res) || |
---|
211 | ((PQresultStatus(res) != PGRES_TUPLES_OK ) && |
---|
212 | (PQresultStatus(res) != PGRES_COMMAND_OK) ) ) |
---|
213 | { |
---|
214 | PQclear(res); |
---|
215 | return 0==1; |
---|
216 | } |
---|
217 | res= PQexec(conn, "set datestyle to 'ISO, YMD';"); |
---|
218 | if (!(res) || |
---|
219 | ((PQresultStatus(res) != PGRES_TUPLES_OK ) && |
---|
220 | (PQresultStatus(res) != PGRES_COMMAND_OK) ) ) |
---|
221 | { |
---|
222 | PQclear(res); |
---|
223 | return 0==1; |
---|
224 | } |
---|
225 | res= PQexec(conn, "set autocommit to 'on';"); |
---|
226 | if (!(res) || |
---|
227 | ((PQresultStatus(res) != PGRES_TUPLES_OK ) && |
---|
228 | (PQresultStatus(res) != PGRES_COMMAND_OK) ) ) |
---|
229 | { |
---|
230 | PQclear(res); |
---|
231 | return 0==1; |
---|
232 | } |
---|
233 | } |
---|
234 | return 0==0; |
---|
235 | } |
---|
236 | |
---|
237 | int checkSQLErrors(PGresult* res){ |
---|
238 | if ( !(res) || |
---|
239 | ((PQresultStatus(res) != PGRES_TUPLES_OK ) && |
---|
240 | (PQresultStatus(res) != PGRES_COMMAND_OK) ) ) |
---|
241 | { |
---|
242 | return 0==1; |
---|
243 | } |
---|
244 | return 0==0; |
---|
245 | } |
---|
246 | |
---|
247 | LIST initColors(EccadDB conn) |
---|
248 | { |
---|
249 | LIST colors= make_list(); |
---|
250 | PGresult *res; |
---|
251 | res=PQexec(conn,"SELECT * FROM color;"); |
---|
252 | if(checkSQLErrors(res)) |
---|
253 | { |
---|
254 | register int i; |
---|
255 | for(i=0; i< PQntuples(res); i++) |
---|
256 | { |
---|
257 | ptrRGB c= alloueColor( |
---|
258 | atoi(PQgetvalue(res, i, 0)), |
---|
259 | atoi(PQgetvalue(res, i, 1)), |
---|
260 | atoi(PQgetvalue(res, i, 2)), |
---|
261 | atoi(PQgetvalue(res, i, 3))); |
---|
262 | add_to_tail(colors, c); |
---|
263 | } |
---|
264 | } |
---|
265 | PQclear(res); |
---|
266 | |
---|
267 | return colors; |
---|
268 | } |
---|
269 | LIST initParams(EccadDB conn) |
---|
270 | { |
---|
271 | LIST params= make_list(); |
---|
272 | PGresult *res; |
---|
273 | res=PQexec(conn,"SELECT * FROM parametre;"); |
---|
274 | if(checkSQLErrors(res)) |
---|
275 | { |
---|
276 | register int i; |
---|
277 | for(i=0; i< PQntuples(res); i++) |
---|
278 | { |
---|
279 | double test = 0.00000001; |
---|
280 | |
---|
281 | fprintf(stderr, " PQgetvalue*--> maxPar %s %s", PQgetvalue(res, i, 6), PQgetvalue(res, i, 3)); |
---|
282 | //fprintf(stderr, " PQgetvalue-> maxPar %f", (double)strtod(PQgetvalue(res, i, 6),NULL)); |
---|
283 | |
---|
284 | ptrParam p= allocParam( |
---|
285 | atoi(PQgetvalue(res, i, 0)), //id |
---|
286 | atoi(PQgetvalue(res, i, 2)), //id unit |
---|
287 | PQgetvalue(res, i, 3), // full name |
---|
288 | atof(PQgetvalue(res, i, 6)), // minPpar |
---|
289 | atof(PQgetvalue(res, i, 7))); // maxPar |
---|
290 | add_to_tail(params, p); |
---|
291 | } |
---|
292 | |
---|
293 | } |
---|
294 | PQclear(res); |
---|
295 | return params; |
---|
296 | } |
---|
297 | |
---|
298 | // generation de la table de couleur d'une map de la base |
---|
299 | // idMap : id de la base |
---|
300 | // mini/maxi : valeurs min max du parametre de la base |
---|
301 | // ces deux valeurs sont nï¿œssaire car les intervalles de couleurs |
---|
302 | // sont codes en unsigned maxhort / au extremes du param |
---|
303 | |
---|
304 | LIST mapOfParamById(EccadDB conn, int idMap, double mini, double maxi) |
---|
305 | { |
---|
306 | LIST map=NULL; |
---|
307 | static LIST colors=NULL; |
---|
308 | PGresult *res, *resc; |
---|
309 | char query[256]; |
---|
310 | |
---|
311 | int isDynamic = 0; |
---|
312 | |
---|
313 | /*****************/ |
---|
314 | mini=-1; |
---|
315 | maxi=100; |
---|
316 | int pas=10; |
---|
317 | /****************/ |
---|
318 | |
---|
319 | // mapping des couleurs de la base |
---|
320 | if(!colors) colors=initColors(conn); |
---|
321 | |
---|
322 | fprintf(stderr, "Generation interval for idMap %d With min %f max %f\n", |
---|
323 | idMap, mini, maxi); |
---|
324 | |
---|
325 | |
---|
326 | //verification de la class map dynamique |
---|
327 | sprintf(query, |
---|
328 | "SELECT * FROM class_map WHERE clm_id=96500"); |
---|
329 | |
---|
330 | //sprintf(query, |
---|
331 | // "SELECT * FROM class_map WHERE clm_id=%d;", idClass); |
---|
332 | |
---|
333 | res=PQexec(conn,query); |
---|
334 | if(checkSQLErrors(res) && PQntuples(res)) |
---|
335 | { |
---|
336 | isDynamic = (!strcmp(PQgetvalue(res, 0, 1),"linear") || !strcmp(PQgetvalue(res, 0, 1),"logarithm") || !strcmp(PQgetvalue(res, 0, 1),"exponential")); |
---|
337 | } |
---|
338 | |
---|
339 | // recherche des interval de la colormap |
---|
340 | sprintf(query,"SELECT * FROM class_index JOIN color USING (col_id) WHERE clm_id=%d order by num_order;", idMap); |
---|
341 | res=PQexec(conn,query); |
---|
342 | if(checkSQLErrors(res) && PQntuples(res)) |
---|
343 | { |
---|
344 | // mise en place des intervalles base dans la liste map |
---|
345 | ptrInterval interv; |
---|
346 | register int i; |
---|
347 | map = make_list(); |
---|
348 | fprintf(stderr, "--------*******--------min %ld max %ld\n", |
---|
349 | v2int(mini,mini,maxi), v2int(maxi,mini,maxi)); |
---|
350 | unsigned int cmin,cmax; |
---|
351 | for(i=0; i< PQntuples(res); i++) |
---|
352 | { |
---|
353 | |
---|
354 | if (i==0) { |
---|
355 | cmin=-1; |
---|
356 | cmax=0; |
---|
357 | } |
---|
358 | else { |
---|
359 | cmin=(i-1)*pas; |
---|
360 | cmax=i*pas; |
---|
361 | } |
---|
362 | |
---|
363 | // conversion en unsigned int des intervalles |
---|
364 | /*int id= atoi(PQgetvalue(res, i, 1)); //idColor |
---|
365 | if ((idMap == 2) || (idMap == 48)) cmin= v2int(atof(PQgetvalue(res, i, 3))* maxi,mini, maxi); |
---|
366 | else cmin= v2int(atof(PQgetvalue(res, i, 3)),mini, maxi); |
---|
367 | if ((idMap == 2) || (idMap == 48)) cmax= v2int(atof(PQgetvalue(res, i, 2))* maxi,mini, maxi); |
---|
368 | else cmax= v2int(atof(PQgetvalue(res, i, 2)),mini, maxi); |
---|
369 | */ |
---|
370 | |
---|
371 | int id=atoi(PQgetvalue(res, i, 1));//idColor |
---|
372 | /*if (isDynamic) cmin= v2int(atof(PQgetvalue(res, i, 3))* maxi,mini, maxi); |
---|
373 | else cmin= v2int(atof(PQgetvalue(res, i, 3)),mini, maxi); |
---|
374 | if (isDynamic) cmax= v2int(atof(PQgetvalue(res, i, 2))* maxi,mini, maxi); |
---|
375 | else cmax= v2int(atof(PQgetvalue(res, i, 2)),mini, maxi); |
---|
376 | */ |
---|
377 | |
---|
378 | if (isDynamic) cmin= v2int(cmin*maxi,mini, maxi); |
---|
379 | else cmin= v2int(cmin,mini, maxi); |
---|
380 | |
---|
381 | if (isDynamic) cmax= v2int(cmax*maxi,mini, maxi); |
---|
382 | else cmax= v2int(cmin,mini, maxi); |
---|
383 | |
---|
384 | //affectation des couleurs aux intervalles |
---|
385 | ptrRGB c; |
---|
386 | if( c= (ptrRGB)search_list(colors, &id, (PFI)inColorId) ) |
---|
387 | { |
---|
388 | interv = alloueInterval( c,cmin, cmax); |
---|
389 | fprintf(stderr, "add interval: cmin %s %ld cmax %s %ld id -> col %d\n", |
---|
390 | PQgetvalue(res, i, 3), cmin, PQgetvalue(res, i, 2), cmax, c->id); |
---|
391 | |
---|
392 | /* |
---|
393 | fprintf(stderr, "add interval: min %f max %f id -> col %d\n", |
---|
394 | int2v(cmin,mini,maxi), int2v(cmax,mini,maxi), c->id); |
---|
395 | */ |
---|
396 | add_to_tail( map , interv); |
---|
397 | } |
---|
398 | else |
---|
399 | { |
---|
400 | fprintf(stderr, "No found color for id %d\n", id); |
---|
401 | exit(1); |
---|
402 | } |
---|
403 | } |
---|
404 | } |
---|
405 | else fprintf(stderr, "No found color for idMap %d\n", idMap); |
---|
406 | |
---|
407 | PQclear(res); |
---|
408 | |
---|
409 | return map; |
---|
410 | } |
---|
411 | |
---|
412 | //generation d'une table de couleur en fobction d'un parametre // |
---|
413 | LIST mapOfParam(EccadDB conn, ptrParam p) |
---|
414 | { |
---|
415 | LIST map=NULL; |
---|
416 | PGresult *res; |
---|
417 | char query[256]; |
---|
418 | |
---|
419 | sprintf(query,"%s \'%s\';", |
---|
420 | "SELECT cm_id FROM param_color NATURAL JOIN parametre WHERE fullName_param = ",p->name); |
---|
421 | res=PQexec(conn,query); |
---|
422 | if(checkSQLErrors(res)) |
---|
423 | map= mapOfParamById(conn, atoi(PQgetvalue(res, 0, 0)), p->mini, p->maxi); |
---|
424 | PQclear(res); |
---|
425 | return map; |
---|
426 | } |
---|
427 | |
---|
428 | // w, h : size of origin grid isur zone globale implicite |
---|
429 | // z : zone a extraire en inputs zone extraite en output |
---|
430 | static ptrGrid loadDatasFromBd(EccadDB db,int oidGrille,int w,int h,ptrZone z) |
---|
431 | { |
---|
432 | // la grille origine sans data alloue |
---|
433 | Grid g; |
---|
434 | g.width=w; |
---|
435 | g.height=h; |
---|
436 | // la grille de sortie |
---|
437 | ptrGrid grid=NULL; |
---|
438 | { |
---|
439 | int id; |
---|
440 | PQexec(db,("BEGIN")); |
---|
441 | if( (id = lo_open(db, oidGrille, INV_READ)) != -1) |
---|
442 | { |
---|
443 | double eps=0.000001; |
---|
444 | |
---|
445 | if ((z != NULL) && memcmp(z, &globalZone, sizeof(Zone)) ) |
---|
446 | { |
---|
447 | int iDeb ,iFin, jDeb, jFin; |
---|
448 | double paslon, paslat; |
---|
449 | float x, y, lat, lon; |
---|
450 | unsigned long K; |
---|
451 | // le pas dorigine |
---|
452 | paslat= pasLat(&globalZone, &g); |
---|
453 | paslon= pasLon(&globalZone, &g); |
---|
454 | // la zone qui doit etre extraite |
---|
455 | // fprintf(stderr, |
---|
456 | // "Extrait Origin Lat %f %f pasLat %.2f lon %f %f pasLon %.2f\n", |
---|
457 | //z->latMin,z->latMax,paslat,z->lonMin,z->lonMax,paslon); |
---|
458 | // la zone ᅵ extraire |
---|
459 | if( z->lonMax != z->lonMin) |
---|
460 | z->lonMax=z->lonMax -eps; |
---|
461 | if( z->latMax != z->latMin) |
---|
462 | z->latMax=z->latMax +eps; |
---|
463 | jDeb= lat2j(z->latMin, &g, &globalZone); |
---|
464 | jFin= lat2j(z->latMax, &g, &globalZone); |
---|
465 | iDeb= lon2i(z->lonMin, &g, &globalZone); |
---|
466 | iFin= lon2i(z->lonMax, &g, &globalZone); |
---|
467 | //la borne en width et height ne peut etre atteinte |
---|
468 | if( iDeb == g.width) iDeb--; |
---|
469 | if( jDeb == g.height) jDeb--; |
---|
470 | if( iFin == g.width) iFin--; |
---|
471 | if( jFin == g.height) jFin--; |
---|
472 | |
---|
473 | // recalcul de la zone extraite |
---|
474 | lon = i2lon(iDeb, &g, &globalZone); |
---|
475 | lat = j2lat(jDeb, &g, &globalZone); |
---|
476 | z->lonMax = i2lon(iFin, &g, &globalZone) ; |
---|
477 | z->latMax = j2lat(jFin, &g, &globalZone) |
---|
478 | // ajout des limites de domaines |
---|
479 | // car les positions sont aux centres pour les macro I2 et j2 |
---|
480 | z->lonMax += paslon/2.; |
---|
481 | z->latMax += + paslat/2.; |
---|
482 | z->lonMin = lon - paslon/2.; |
---|
483 | z->latMin = lat - paslat/2.; |
---|
484 | fprintf(stderr, |
---|
485 | "----------------------------------------------------------------------------Extrait Reel Lat min %f max %f lon min %f max %f width %d height %d \n", |
---|
486 | z->latMin,z->latMax,z->lonMin,z->lonMax, iFin -iDeb +1, jFin -jDeb +1); |
---|
487 | |
---|
488 | // extraction dans une nouvelle grille |
---|
489 | grid= allocGrid(iFin -iDeb +1, jFin -jDeb +1);fprintf(stderr,"ok grid \n"); |
---|
490 | if(iFin < iDeb) |
---|
491 | { |
---|
492 | freeGrid(grid); |
---|
493 | grid=NULL; |
---|
494 | fprintf(stderr," Not now ok on partial grid crossing in longitude \n iDeb %d Ifin %d jDeb %d jFin %d\n", iDeb ,iFin, jDeb, jFin); |
---|
495 | } |
---|
496 | else |
---|
497 | { |
---|
498 | int i,j, pos, skip, width; |
---|
499 | unsigned char *data=(unsigned char *)(grid->datas); |
---|
500 | pos= (g.width*jDeb + iDeb) * sizeof(unsigned int); |
---|
501 | lo_lseek(db, id, pos, SEEK_SET); |
---|
502 | width= (iFin -iDeb +1) * sizeof(unsigned int); |
---|
503 | skip= g.width * sizeof(unsigned int) - width ; |
---|
504 | // fprintf(stderr," ok pos %d width %d skip %d\n", pos/4 , width/4 , skip/4); |
---|
505 | for(i=jDeb; i<=jFin; i++, data += width) |
---|
506 | { |
---|
507 | if((K= lo_read(db, id, (unsigned char*)data, width)) != (unsigned long)width) |
---|
508 | { |
---|
509 | fprintf(stderr, "Err readLine ligne %d read %ld sur %d demandes\n", i, K, width); |
---|
510 | fprintf(stderr, "EXIT trap:\n iDeb %d iFin %d jDeb %d jFin %d\n pos %d width %d skip %d\n", |
---|
511 | iDeb , iFin , jDeb , jFin , pos , width , skip); |
---|
512 | exit (1); |
---|
513 | } |
---|
514 | lo_lseek(db, id, skip, SEEK_CUR); |
---|
515 | } |
---|
516 | } |
---|
517 | } |
---|
518 | else |
---|
519 | { |
---|
520 | unsigned long K, size= g.width * g.height * sizeof(int); |
---|
521 | grid= allocGrid(g.width , g.height); |
---|
522 | if( (K= lo_read(db, id, (char*)grid->datas, size)) != size) |
---|
523 | fprintf(stderr," Erreur Read OID : %ld sur %d octets\n", K, size); |
---|
524 | // else |
---|
525 | // fprintf(stderr," Full Read OID : %ld sur %d octets\n", K, size); |
---|
526 | } |
---|
527 | lo_close(db,id); |
---|
528 | } |
---|
529 | else fprintf(stderr," Erreur Open OID : %d\n", oidGrille); |
---|
530 | PQexec(db,("END")); |
---|
531 | } |
---|
532 | fprintf(stderr,"Open OID \n"); |
---|
533 | return grid; |
---|
534 | } |
---|
535 | |
---|
536 | ptrGridGL extractGridFromDb(EccadDB db, int idGrille, ptrZone z) |
---|
537 | { |
---|
538 | static LIST params=NULL; |
---|
539 | ptrGridGL gl=NULL; |
---|
540 | PGresult *res; |
---|
541 | char query[256]; |
---|
542 | |
---|
543 | sprintf(query, |
---|
544 | "SELECT valeurs_grille, id_param, ncol_grille, nlign_grille, date_grille FROM grille WHERE id_grille = %d;", idGrille); |
---|
545 | res=PQexec(db,query); |
---|
546 | if(checkSQLErrors(res) && PQntuples(res)) |
---|
547 | { |
---|
548 | ptrParam p; |
---|
549 | int oid= atoi(PQgetvalue(res, 0, 0)); |
---|
550 | int idParam= atoi(PQgetvalue(res, 0, 1)); |
---|
551 | int w= atoi(PQgetvalue(res, 0, 2)); |
---|
552 | int h= atoi(PQgetvalue(res, 0, 3)); |
---|
553 | if(!params) params=initParams(db); |
---|
554 | if( p= (ptrParam)search_list(params, &idParam, (PFI)inParamId) ) |
---|
555 | { |
---|
556 | ptrGrid g; |
---|
557 | Zone zoneExtr; |
---|
558 | g= loadDatasFromBd( db, oid, w, h, z); |
---|
559 | strcpy(g->date, PQgetvalue(res, 0, 4) ); |
---|
560 | gl= allocGridGL(g, z); |
---|
561 | gl->param=p; |
---|
562 | } |
---|
563 | else fprintf(stderr, " No found para for id %d\n", idParam); |
---|
564 | |
---|
565 | } |
---|
566 | PQclear(res); |
---|
567 | return gl; |
---|
568 | } |
---|
569 | |
---|
570 | LIST getGrids(EccadDB db, int idProduct, ptrParam p, ptrZone z, |
---|
571 | const char* dateBegin, const char* dateEnd) |
---|
572 | { |
---|
573 | int i; |
---|
574 | LIST grids=NULL; |
---|
575 | PGresult *res; |
---|
576 | char query[256]; |
---|
577 | sprintf(query, |
---|
578 | "SELECT valeurs_grille, ncol_grille, nlign_grille, date_grille FROM grille WHERE id_param = %d AND id_produit = %d AND date_grille BETWEEN '%s' AND '%s' order by date_grille", |
---|
579 | p->id, idProduct, dateBegin, dateEnd); |
---|
580 | res=PQexec(db,query); |
---|
581 | fprintf(stderr, "found %d grilles SQLErrors %d\n", PQntuples(res), checkSQLErrors(res)); |
---|
582 | if(checkSQLErrors(res) && PQntuples(res) ) |
---|
583 | { |
---|
584 | // origin size |
---|
585 | int width= atoi(PQgetvalue(res,0,1)); |
---|
586 | int height= atoi(PQgetvalue(res,0,2)); |
---|
587 | for(i=0; i< PQntuples(res); i++) |
---|
588 | { |
---|
589 | int oid= atoi(PQgetvalue(res,i,0) ); |
---|
590 | /*fprintf(stderr, "Load oid %d from %dx%d grille date %s\n", |
---|
591 | oid, width, height, PQgetvalue(res, i, 3) );*/ |
---|
592 | |
---|
593 | fprintf(stderr,"Avant--------- : Extrait Reel Lat min %f max %f lon min %f max %f \n",z->latMin,z->latMax,z->lonMin,z->lonMax); |
---|
594 | |
---|
595 | |
---|
596 | ptrGrid g= loadDatasFromBd( db, oid, width, height, z); |
---|
597 | |
---|
598 | fprintf(stderr,"Aprï¿œs--------- : Extrait Reel Lat min %f max %f lon min %f max %f \n",z->latMin,z->latMax,z->lonMin,z->lonMax); |
---|
599 | // printf("date grille %d\n",g->date); |
---|
600 | if(g) |
---|
601 | { |
---|
602 | char *date= PQgetvalue(res, i, 3); |
---|
603 | sprintf(g->date, "%.4s%.2s%.2s", date, &date[5], &date[8]); |
---|
604 | if(!grids) grids = make_list(); |
---|
605 | add_to_tail(grids, g); |
---|
606 | } |
---|
607 | else fprintf(stderr, "loadGrid error\n"); |
---|
608 | } |
---|
609 | } |
---|
610 | PQclear(res); |
---|
611 | return grids; |
---|
612 | } |
---|
613 | |
---|
614 | |
---|
615 | |
---|
616 | ptrGrid getGrid(EccadDB db, int idGrille, ptrZone z) |
---|
617 | { |
---|
618 | fprintf(stderr, "--getGrid \n");fflush(stderr); |
---|
619 | PGresult *res; |
---|
620 | char query[256]; |
---|
621 | sprintf(query, |
---|
622 | "SELECT valeurs_grille, ncol_grille, nlign_grille, date_grille FROM grille WHERE id_grille = %d;", idGrille); |
---|
623 | fprintf(stderr, "1)-->requete %s \n",query ); |
---|
624 | res=PQexec(db,query); |
---|
625 | fprintf(stderr, "found %d grilles SQLErrors %d\n", PQntuples(res), checkSQLErrors(res)); |
---|
626 | fprintf(stderr, "2)-->requete %s \n",query ); |
---|
627 | //fprintf(stderr, "found %d grilles SQLErrors %d\n", PQntuples(res), checkSQLErrors(res)); |
---|
628 | ptrGrid g = NULL; |
---|
629 | |
---|
630 | if(checkSQLErrors(res) && PQntuples(res) ) |
---|
631 | { |
---|
632 | fprintf(stderr, "organise \n"); |
---|
633 | // origin size |
---|
634 | int width= atoi(PQgetvalue(res,0,1)); |
---|
635 | int height= atoi(PQgetvalue(res,0,2)); |
---|
636 | |
---|
637 | int oid= atoi(PQgetvalue(res,0,0) ); |
---|
638 | //fprintf(stderr,"Avant--------- : Extrait Reel Lat min %f max %f lon min %f max %f \n",z->latMin,z->latMax,z->lonMin,z->lonMax); |
---|
639 | |
---|
640 | g= loadDatasFromBd( db, oid, width, height, z); |
---|
641 | |
---|
642 | //fprintf(stderr,"Apres--------- : Extrait Reel Lat min %f max %f lon min %f max %f \n",z->latMin,z->latMax,z->lonMin,z->lonMax); |
---|
643 | |
---|
644 | |
---|
645 | if(g) |
---|
646 | { |
---|
647 | char *date= PQgetvalue(res, 0, 3); |
---|
648 | sprintf(g->date, "%.4s%.2s%.2s", date, &date[5], &date[8]); |
---|
649 | |
---|
650 | fprintf(stderr, "date grille %s \n", date); |
---|
651 | } |
---|
652 | } |
---|
653 | |
---|
654 | return g; |
---|
655 | } |
---|
656 | |
---|
657 | |
---|
658 | double* getMaxParam(int idGrid, EccadDB conn) |
---|
659 | { |
---|
660 | double* maxParam = NULL; |
---|
661 | PGresult *res; |
---|
662 | char query[256]; |
---|
663 | sprintf(query,"SELECT id_param FROM grille WHERE id_grille=%d;", idGrid); |
---|
664 | res=PQexec(conn,query); |
---|
665 | |
---|
666 | if(checkSQLErrors(res) && PQntuples(res) ) |
---|
667 | { |
---|
668 | int idParam= atoi(PQgetvalue(res,0,0)); |
---|
669 | |
---|
670 | ptrParam p; |
---|
671 | LIST params=initParams(conn); |
---|
672 | if( p= (ptrParam)search_list(params, (NODE)&idParam, (PFI)inParamId) ) |
---|
673 | maxParam = &p->maxi; |
---|
674 | } |
---|
675 | return maxParam; |
---|
676 | } |
---|
677 | |
---|
678 | double* getMinParam(int idGrid, EccadDB conn) |
---|
679 | { |
---|
680 | double* minParam = NULL; |
---|
681 | PGresult *res; |
---|
682 | char query[256]; |
---|
683 | sprintf(query,"SELECT id_param FROM grille WHERE id_grille=%d;", idGrid); |
---|
684 | res=PQexec(conn,query); |
---|
685 | |
---|
686 | if(checkSQLErrors(res) && PQntuples(res) ) |
---|
687 | { |
---|
688 | int idParam= atoi(PQgetvalue(res,0,0)); |
---|
689 | |
---|
690 | ptrParam p; |
---|
691 | LIST params=initParams(conn); |
---|
692 | if( p= (ptrParam)search_list(params, (NODE)&idParam, (PFI)inParamId) ) |
---|
693 | minParam = &p->mini; |
---|
694 | } |
---|
695 | |
---|
696 | return minParam; |
---|
697 | } |
---|
698 | |
---|