vanessa15 a écrit 13 commentaires

Bonjour,

Je n'ai pas eu d'ordinateur pendant quelques jours donc je n'ai pas pu avancer.... As tu essayé de ton côté?
J'ai peur que l'installation de GDAL ne se fasse pas correctement, et du coup le problème viendrai de là ;

quand je compile GDAL je tape ;

./configure --prefix=INSTALL_DIR --with-png=internal --with-libtiff=internal --with-jpeg=internal

il compile et tout est ok. Par contre quand je tape make, il compile pendant presque une heure et des pages de charabia défilent, voici ce que m'affiche l'écran pour ne copier que la dernière ligne:

g++ gdal_rasterize.o -o gdal_rasterize.exe /cygdrive/d/OTB/INSTALL_DIR/gdal-1.4.1/.libs/libgdal.a -lz -ldl : Leaving directory '/cygdrive/d/OTB/INSTALL_DIR/gdal-1.4.1/apps'

J'ai également vu s'afficher un

make [1] : enterring directory ......

Qu'est ce que cela veut dire?

Comment utiliser la même version pour les librairies et mon programme? Ils sont indépendants les uns des autres, je télécharge les codes sources et je prend la dernière version à chaque fois.
Il faut que je choisisse des numéros de version identique?
J'ai également choisi les versions qui sont utilisés dans la "notice d'installation", mais cela ne change rien

Vérifie si ce n'est déjà fait que ta version de libgdal fourni bien les fonctionnalités demandées.

Comment faire cela? Si j'ai bien compris libgdal est une librairie de GDAL, il faut donc que je regarde dans le fichier libgdal? Ce que j'ai fait mais il est illisible (avec word pad) et si c'est écrit en C, je ni comprendrais rien!!


Petite précision ; je ne suis pas très douée en programmation ;)

Merci pour cette aide....

Je faisais effectivement référence aux messages comme "Image cree à partir de: 5774570 valeurs" .

Si ces paramètres n'interviennent pas dans les calculs de syntheseV alors c'est génial, je continue comme ça.

Merci beaucoup pour ton aide (encore une fois).
:)

C'est quoi au juste un apendfile?

Malheureusement je ne sais pas quels résultats je dois obtenir, puisque c'est l'objet de mon travail ; étudier des images satellites et obtimiser celles-ci en trouvant une période et une maille correcte. Je n'ai aucun moyen de comparaison, ces données n'ont jamais été utilisées auparavant.

Effectivement j'obtiens des résultats diffférents pour des périodes en entrée différentes (5j,10j...). Donc à priori le programme prend en compte tout ce qu'il ya dans mon fichier en entrée. Mais suis je bien sure qu'il prend TOUT en compte et que cela ne dépend pas (même peu) de mes paramètres date_deb et nbjour. Car je dois en avoir la certitude. De toute façon cela n'est pas très rigoureux d'avoir 2 paramètres "inutiles"....

En ce qui concerne le nombre de valeur à partir de laquelle l'image est crée le problème venait du fait que le fichier en sortie n'écrase pas ce qu'il y avait avant, je dois créer un nouveau fichier. Et alors quelque soit la valeur des paramètres tappés j'obtiens toujours la même valeur .

Il n'est pas possible de faire les choses à l'envers, c'est à dire de lancer sytheseV "manuellement"ou quelque chose dans le genre ?

Je pense qu'il interprète correctement les noms de fichiers.

J'ai lancé MainScript en tapant des paramètres différents par exemple:
MainScript 1 15 0.125 -25 20 0 25 /mnt/porto/sde1/vanessa/Save/Area/15j/ /mnt/porto/sde1/vanessa/tmp/essai_1_15/

Voici le début et la fin de ce qu'il fait:

"Fichier courant: /mnt/porto/sde1/vanessa/Save/Area/15j/MYQS_S2A23613.20040011144 JOUR=1 ANNEE=2004 date_deb+nbjour=16

FILENAME: /mnt/porto/sde1/vanessa/Save/Area/15j/MYQS_S2A23613.20040011144
sp= 55
Fichier courant: /mnt/porto/sde1/vanessa/Save/Area/15j/MYQS_S2A23614.20040011145 JOUR=1 ANNEE=2004 date_deb+nbjour=16

FILENAME: /mnt/porto/sde1/vanessa/Save/Area/15j/MYQS_S2A23614.20040011145
sp= 127

(.............)

Fichier courant: /mnt/porto/sde1/vanessa/Save/Area/15j/MYQS_S2A23821.20040160349 JOUR=16 ANNEE=2004 date_deb+nbjour=16

FILENAME: /mnt/porto/sde1/vanessa/Save/Area/15j/MYQS_S2A23821.20040160349
sp= 127
Fichier courant: /mnt/porto/sde1/vanessa/Save/Area/15j/MYQS_S2A23822.20040160350 JOUR=16 ANNEE=2004 date_deb+nbjour=16

FILENAME: /mnt/porto/sde1/vanessa/Save/Area/15j/MYQS_S2A23822.20040160350
sp= 30
******PROCESSING********
********************
Creation de l'image
/mnt/porto/sde1/vanessa/tmp/essai_1_15/e2004_1_16_ecartypeHH.raw
/mnt/porto/sde1/vanessa/tmp/essai_1_15/e2004_1_16_ecartypeVV.raw
/mnt/porto/sde1/vanessa/tmp/essai_1_15/e2004_1_16_nbvalueVV.raw
/mnt/porto/sde1/vanessa/tmp/essai_1_15/e2004_1_16_nbvalueHH.raw
il ya 128 pixel sans valeur dans /mnt/porto/sde1/vanessa/tmp/essai_1_15/e2004_1_16_imageHH.raw
il ya 123 pixel sans valeur dans /mnt/porto/sde1/vanessa/tmp/essai_1_15/e2004_1_16_imageVV.raw
col=360 ligne=200

Image cree à partir de: 5774570 valeurs"

Je n'ai pas tout mis , c'est beaucoup trop long.
Quelque soit les paramètres que je tappe j'obtiens toujours la même chose, seul date_deb+nbjour change.
Je viens de remarquer que le nombre de valeur à partir de laquelle l'image est crée peut changer mais ceci ne se voit pas dans les images.....
Par exemple, en tapant 1 et 5 j'obtiens:

"Image cree à partir de: 2887285 valeurs"

Mais je ne sais pas si cela à de l'importance ou non.....




[

Rien de bien sorcier...... :)

En fait je doit générer des images satellites sur des périodes différentes, 5jours, 10 jours .....
Le problème c'est que le programme aurait été fait de tel sorte que le paramètre date_deb serait le jour julien et que "nbjour" serait un pas. Donc quand je rentre des données sur 15 jours par exemple, je ne devrait pas obtenir la même chose selon la date du jour et le pas. Or même si je tappe 0 0 j'obtient la même chose, en fait il ne tient compte que de la taille du fichiers en entrée.

Je voudrais être sur que se soit bien ça, cela ne me pose pas de problème qu'il n'y ai pas de pas mais j'ai quand même deux paramètres qui ne servent à rien!

Dans le script , on dit que si la valeur jour qu'on a trouvé dans le nom de fichier est superieure ou egale
# a date_deb et strictement inferieure a MAX:
# on appelle concat dans le dossier original avec le nom du fichier donc cela devrait changé selon les paramètres tapés?
Ou alors cela n'a d'impact que sur le nom du fichier en sortie? Car effectivement mes fichiers sont de la forme "le2004_1_6_imageHH.raw " avec 1 date_deb et 6 date_deb + nbjour.

Dans syntheseV on ne parle pas de ses deux paramètres donc je ne pense pas que ça vienne de là.....

Salut

Le problème venait effectivement de la variable i.
Maintenant tout se fait sans erreur et j'obtient de belles images!!
Merci beaucoup pour ton aide.

Trop fort!

J'essayerai ca lundi, le problème c' est qu'on a plus de nouvelles du programmateur!!
Mais certaines personnes on suivit le programme de loin....On devrait y arriver enfin j'espère!

En tout cas merci beaucoup...

OK
J'essayerai demain,
Je suis chez moi et je n'ai pas les fichiers........

Génial,

Les paramètres ;
[date_deb] 1
[nb_jours] 15
[maille] 0.125
[long_min]-25
[long_max]20
[lat_min]0
lat_max]25
[dossier entre] fichier contenant des données
[dossier sortie]fichier contenant ts les fichiers images
[appendfile]fichier créé par concat.

Sinon en tappant ce que tu as écrit j'obtient la même chose que toi sauf bien sur quand je lance MainScript, toujours le même résultat....

J'espère qu'il ne te manque rien.....

Avec la commande bt, il m'affiche "no stack"

Une autre remarque, quand j'ouvre mes images celles avec image.raw et nbvalue.raw marchent par contre c'est celles avec ecarttype qui ne fonctionnent pas, le fichier est vide!

Egalement, je sais pas si ca a un rapport mais quand je lance le Mainscript, le chiffre devant segmentation fault n'ai jamais le même, tu sais ce qu'il signifie?

Voici le fichier hdf.h:

/****************************************************************************
* NCSA HDF *
* Software Development Group *
* National Center for Supercomputing Applications *
* University of Illinois at Urbana-Champaign *
* 605 E. Springfield, Champaign IL 61820 *
* *
* For conditions of distribution and use, see the accompanying *
* hdf/COPYING file. *
* *
****************************************************************************/

/* $Id: hdf.h,v 1.61.2.2 2003/11/11 04:35:55 bmribler Exp $ */

#ifndef HDF_H
#define HDF_H

#include "h4config.h"
#include "hdfi.h"
#include "hlimits.h"

/* Internal DF structure */
typedef struct
{
uint16 tag; /* tag of element */
uint16 ref; /* ref of element */
}
DFdi;

/* For annotations */
/* enumerated types of the varous annotation types */
typedef enum
{
AN_UNDEF = -1,
AN_DATA_LABEL = 0, /* Data label */
AN_DATA_DESC, /* Data description */
AN_FILE_LABEL, /* File label */
AN_FILE_DESC /* File description */
} ann_type;

/* internal file access codes */

#define DFACC_READ 1
#define DFACC_WRITE 2
#define DFACC_CREATE 4
#define DFACC_ALL 7

#define DFACC_RDONLY 1
#define DFACC_RDWR 3
#define DFACC_CLOBBER 4

/* New file access codes (for Hstartaccess only, currently) */
#define DFACC_BUFFER 8 /* buffer the access to this AID */
#define DFACC_APPENDABLE 0x10 /* make this AID appendable */
#define DFACC_CURRENT 0x20 /* start looking for a tag/ref from the current */
/* location in the DD list (useful for continued */
/* searching ala findfirst/findnext) */

/* External Element File access mode */
/* #define DFACC_CREATE 4 is for creating new external element file */
#define DFACC_OLD 1 /* for accessing existing ext. element file */

/* The magic cookie for Hcache to cache all files */
#define CACHE_ALL_FILES (-2)

/* File access modes */
/* 001--007 for different serial modes */
/* 011--017 for different parallel modes */

#define DFACC_DEFAULT 000
#define DFACC_SERIAL 001
#define DFACC_PARALLEL 011

/* used by Hnextread to determine where to start searching for the
next tag/ref to read */

#define DF_START 0
#define DF_CURRENT 1
#define DF_END 2

/* Used by Hfind to determine the direction to search for tag/ref's in the */
/* file. */

#define DF_FORWARD 1
#define DF_BACKWARD 2

/* return code - since some unix/c routines use 0 and -1 as their return
code, and some assumption had been made in the code about that, it is
important to keep these constants the same values. For explicitly
boolean functions, use TRUE and FALSE */

#define SUCCEED 0
#define FAIL (-1)

/* boolean values, reminder: NEVER compare with numeric values */

#ifndef FALSE
# define FALSE 0
#endif
#ifndef TRUE
# define TRUE (!FALSE)
#endif

/* macros */
#define STREQ(s, t) (HDstrcmp((s), (t)) == 0)
#define NSTREQ(s, t, n) (HDstrncmp((s), (t), (n)) == 0)

/*
* Macros used for variable and function scoping in code.....
*/
#ifndef EXPORT
#define EXPORT
#endif

#ifndef PRIVATE
#define PRIVATE static
#endif

/* Include the Number-type definitions */
#include "hntdefs.h"

/* Include the Tag definitions */
#include "htags.h"

/*
* interlacing supported by the vset.
*/

#define FULL_INTERLACE 0
#define NO_INTERLACE 1

/* type for File ID to send to Hlevel from Vxx interface */
typedef int32 HFILEID;

typedef intn (*hdf_termfunc_t)(void); /* termination function typedef */

/* .................................................................. */

/* Publically accessible functions declarations. This includes all the
functions that are used by application programs. */

#include "hbitio.h"
#include "hcomp.h"
#include "herr.h"
#include "hproto.h"
#include "vg.h" /* Add the Vgroup/Vdata header so the users don't have to */
#include "mfgr.h" /* Add the GR header so the users don't have to */

/* For Pablo Instrumentation */
#ifdef HAVE_PABLO
#include "ProcIDs.h"
#include "trace.h"
#endif /* PABLO */

/* these may eventaully evolve into real-life functions but not yet */
#define HDFopen(f,a,d) Hopen((f), (a), (d))
#define HDFclose(f) Hclose((f))
#define Vstart(f) Vinitialize((f))
#define Vend(f) Vfinish((f))

/* Misc. macros for backward compability */
#define HDgettagname(tag) HDgettagdesc(tag)

/* This is also defined in fmpio.h */
#define MP_PAGEALL 0x01 /* page the whole file i.e. no limit on 'maxcache'*/

#endif /* HDF_H */



Par contre je sais pas comment t'envoyer le fichier raw , c'est écarttype.raw(il est vide)? Je peux que faire des copier coller et tu pourra pas l'exécuter?

merci de ton aide...ca fait plusieurs jours que je galère.....

Le code compilé on ne peut pas le lire disons que c'est du charabia puisque c'est un executable.....donc je ne vois pas ce que tu peux en faire?
Je suppose que c'est syntheseW.c ? C'est assez long mais bon pourquoi pas!
Le voici:
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include </mnt/porto/sdb1/otig/frison/lib/hdf4/hdf.h>
#include <dirent.h>
#include <limits.h>
#define MIN_SIG -32.0
#define MAX_SIG -5.0

#define MAX_LON 359.99

#define MIN_PIX 5
#define MAX_PIX 255

#define TERRE 1027.0

#define TEMP_FILE_NAME "appendfile"

#define MAX_COLS 810



struct save {
float64 lat;
float64 lon;
float64 sigma0;
float64 incidence;
float64 surface;
float64 sigma0_attn_map;
};
typedef struct save data;


struct param {
int date_deb;
int nbjour;
double maille;
double lon[2];
double lat[2];
char * dossier;
char * sortie;
char * appendfilename;
};
typedef struct param parameter;


int getParameter(int argc,char *argv[],parameter *p){
if(argc==11){
p->date_deb=atoi(argv[1]);
p->nbjour=atoi(argv[2]);
p->maille=atof(argv[3]);
p->lon[0]=atof(argv[4]);
p->lon[1]=atof(argv[5]);
p->lat[0]=atof(argv[6]);
p->lat[1]=atof(argv[7]);
p->dossier=(char*)malloc(sizeof(char)*(strlen(argv[8]+1)));
strcpy(p->dossier,argv[8]);
p->sortie=(char*)malloc(sizeof(char)*(strlen(argv[9]+1)));
strcpy(p->sortie,argv[9]);
p->appendfilename=(char*)malloc(sizeof(char)*(strlen(argv[10]+1)));
strcpy(p->appendfilename,argv[10]);
} else{
fprintf(stderr,"./syntheseW [date_deb][nb_jours][maille][long_min][long_max][lat_min][lat_max][dossier entre][dossier sortie][appendfile]\n");
return -1;
}
return 0;
}

float64 correctS(float64 sigma, float64 sigma_attn, float64 incidence){
return (sigma + 2 * (sigma_attn/cos((double)((2*3.141592654*incidence)/360.))) );
}


double getPixelFromNumber(double val){
return val;
}

double getPixelValue4(double val){
/* pour obtenir des couleurs claires*/
return ((val-MIN_SIG)*(250))/(MAX_SIG-MIN_SIG)+5;
}
double INVgetPixelValue4(double val){
return ((val-5.0)*(MAX_SIG-MIN_SIG)/250.0)+MIN_SIG;
}

double getPixelValue5(double val){
/* pour obtenir des couleurs sombres*/
return (260-(((val-MIN_SIG)*(250))/(MAX_SIG-MIN_SIG)+5));
}


int histoStat(char *picturename,char *statfilename,int linecol[2]){

int i;
FILE *file=fopen(statfilename,"ab");
FILE *img=fopen(picturename,"r");

float tmp;
int p5pourcent;
int p5p=-1;
int p95pourcent;
int p95p=-1;
int histo[256];
int histoC[256];
int max=INT_MIN;
int min=INT_MAX;

for(i=0;i<256;++i){
histo[i]=0;
histoC[i]=0;
}

unsigned char pic[linecol[0]*linecol[1]];
fread(pic,sizeof(char),linecol[0]*linecol[1],img);

int current;
int nbvalue;
for(i=0;i<linecol[0]*linecol[1];++i){
current=(int)pic[i];
if(current>4){
nbvalue++;
if(current>max) max=current;
else if(current<min) min=current;
}
}

for(i=0;i<linecol[0]*linecol[1];++i){

current=(int)pic[i];
if(current>4){
histo[current]++;
}
}


histoC[5]=histo[5];

for(i=6;i<256;++i){
histoC[i]=histoC[i-1]+ histo[i];
}


/*On commence a 5 pour ne pas compter les points de la mer */
for(i=5;i<256;++i){
if(histo[i]>histo[max]) max=i;
tmp=histoC[i]/nbvalue;
tmp*=100;

if(tmp>=5 && p5p==-1){
p5pourcent=i;
p5p=0;
}

if(tmp>=95 && p95p==-1){
p95pourcent=i;
p95p=0;
}

}
printf("%d\n",nbvalue);
for(i=0;i<50;++i){printf("%d: %d %d %f\n",i,histo[i],histoC[i],(histoC[i]*1.0)/(nbvalue*1.0));}

fprintf(file,"*********************\n");
fprintf(file,"Image: %s \n",picturename);
fprintf(file,"min pix: %d \t max pix: %d \n",min,max);
//float nb=histo[p5pourcent]/nbvalue;
float nb=(histo[p5pourcent]*1.0)/(nbvalue*1.0);
fprintf(file,"à 5pourcent:%d mesures(%f pourcent)\n",histo[p5pourcent],nb);
nb=(histo[p95pourcent]*1.0)/(nbvalue*1.0);
fprintf(file,"à 95pourcent:%d mesures(%f pourcent)\n",histo[p95pourcent],nb);
nb=(histo[max]*1.0)/(nbvalue*1.0);
fprintf(file,"max pour %d: %d mesures(%f pourcent)\n",max,histo[max],nb);
fprintf(file,"*********************\n");
return 1;
}

int Statistiques(char *picturename,char *statfilename,int linecol[2],double incidence[2]){

int i;

int current;

int max=INT_MIN;
int min=INT_MAX;
double moyenne=0;
double variance=0;
double ecart_type;
int nbvalue=0;

FILE *file=fopen(statfilename,"ab");
FILE *img=fopen(picturename,"r");

unsigned char pic[linecol[0]*linecol[1]];
fread(pic,sizeof(char),linecol[0]*linecol[1],img);

for(i=0;i<linecol[0]*linecol[1];++i){
current=(int)pic[i];
if(current>4){
nbvalue++;
if(current>max) max=current;
else if(current<min) min=current;
moyenne+=current;
}
}
moyenne=moyenne/(nbvalue);

for(i=0;i<linecol[0]*linecol[1];++i){
current=(int)pic[i];
if(current>4){
variance+=(pow(current-moyenne,2));
}
}
variance=variance/(nbvalue);
ecart_type=sqrt(variance);

fprintf(file,"*********************\n");
fprintf(file,"Image: %s \n",picturename);
fprintf(file,"min pix: %d \t max pix: %d \n",min,max);
fprintf(file,"min val: %f \t max val: %f \n",INVgetPixelValue4(min),INVgetPixelValue4(max));
fprintf(file,"moyenne pix: %f \n",moyenne);
fprintf(file,"moyenne val: %f \n",INVgetPixelValue4(moyenne));
fprintf(file,"ecart-type pix: %f \n",ecart_type);
fprintf(file,"incidence min : %f \t incidence max : %f \n",incidence[0],incidence[1]);


fclose(img);
fclose(file);

return 0;
}


/*
*non alloué lors de l'appel de la fonction*
double **somme: tableau de pixel contenant( apres execution de la methode) pour chaque pixel i la somme des sigma0 naturelles
double **nbval: tableau indiquant (apres execution de la methode) le nombre de valeurs additionnées pour le pixel i
====> Il ne faut donc pas oublié de les liberer apres l'execution

*directement utilisable*
double *lat: [0]= lat min,[1]=lat max
double *lon: [0]= lon min,[1]=lon max
int *lignecol: [0]= nbr de ligne,[1]= nbr de colonne
double maille: taille de la maille
FILE *file: fichier d'entré

Return: nombre d'entré utilisable

Cette fonction récupère les données dans le flux *file et remplit le tableau somme et nbval
Ainsi pour un pixel i donné il suffit de faire somme[i]/nbval[i] pour obtenir la moyenne des sigma0 naturelles associées à ce pixel
*/
int readFile(double **sommeHH,double **nbvalHH,double **sommeVV,double **nbvalVV,parameter p,int *lignecol,char *filename,double **incidence){
data mydata;
int x,y;
int k=0;
FILE *file;

/*printf("coco: %s\n",filename);*/
/*file=fopen(filename,"r");*/

file=fopen(filename,"r");
if(file==NULL) return 1;

/*
NE PAS OUBLIER DE LIBERER CES TABLEAUX
*/
*nbvalHH=(double*)calloc(lignecol[0]*lignecol[1],sizeof(double));
*sommeHH=(double*)calloc(lignecol[0]*lignecol[1],sizeof(double));
*nbvalVV=(double*)calloc(lignecol[0]*lignecol[1],sizeof(double));
*sommeVV=(double*)calloc(lignecol[0]*lignecol[1],sizeof(double));

printf("******PROCESSING********\n");


while(feof(file)==0)
{

/*Pour chaque entrée dans le fichier concatené*/
fread(&mydata, sizeof(data),1,file);
double longi;
/*cette manipulation est necessaire pour passer la longitude dans le meme systeme que la lat*/

if(mydata.lon>p.lon[1]){longi=(mydata.lon-MAX_LON); }
else {longi=mydata.lon;}

/*On verifie que l'entrée est dans la zone souhaité*/

if((mydata.lat<=p.lat[1] && mydata.lat>=p.lat[0])&&(longi<=p.lon[1] && longi>=p.lon[0]) )
{
/* printf("surface %f \n",mydata.surface); */
/* calcule des coordonnées en pixel a partir des long/lat de l'entré*/
x=fabs((p.lon[0]-longi)/p.maille);
y=fabs((p.lat[1]-mydata.lat)/p.maille);

/*printf("********************\n");
printf("lat= %f lon=%f\n",mydata.lat,mydata.lon);
printf("sigma0db= %f \n",mydata.sigma0);
printf("sigma0= %f \n",pow(10.0,mydata.sigma0/10));
printf("%d %d \n",x,y);*/
if (mydata.surface==TERRE){
++k;
if(mydata.incidence<50){
(*sommeHH)[y*lignecol[1]+x]+=pow(10.0,correctS(mydata.sigma0,mydata.sigma0_attn_map,mydata.incidence)/10);
(*nbvalHH)[y*lignecol[1]+x]+=1;
if(mydata.incidence<(*incidence)[0]) (*incidence)[0]=mydata.incidence;
else if (mydata.incidence>(*incidence)[1]) (*incidence)[1]=mydata.incidence;
}
else{
(*sommeVV)[y*lignecol[1]+x]+=pow(10.0,correctS(mydata.sigma0,mydata.sigma0_attn_map,mydata.incidence)/10);
(*nbvalVV)[y*lignecol[1]+x]+=1;
if(mydata.incidence<(*incidence)[2]) (*incidence)[2]=mydata.incidence;
else if (mydata.incidence>(*incidence)[3]) (*incidence)[3]=mydata.incidence;
}

}
/*Si l'entré est dans la mer*/
else{
if(mydata.incidence<50){
(*sommeHH)[y*lignecol[1]+x]=-1;
(*nbvalHH)[y*lignecol[1]+x]=-1;
}
else{
(*sommeVV)[y*lignecol[1]+x]=-1;
(*nbvalVV)[y*lignecol[1]+x]=-1;
}
}


}

}
fclose(file);
return k;
}


int writePicture(double *somme,double *nbval,int lignecol[2],char *filename, char *nbvaluefilename ){

int tmmp=0;

FILE *out=fopen(filename,"w");
FILE *out2=fopen(nbvaluefilename,"w");
int i;
for(i=0;i<lignecol[0]*lignecol[1];++i){

if(nbval[i]>0){
double tmp=getPixelValue4(10.0*log10(somme[i]/nbval[i]));
unsigned char p=(unsigned char)tmp;
fprintf(out,"%c",p);
fprintf(out2,"%c",(unsigned char)getPixelFromNumber(nbval[i]));
}

/*Pas de valeur pour ce pixel*/
else if(nbval[i]==0) {
tmmp++;
fprintf(out,"%c",4);
fprintf(out2,"%c",0);
}

/*Mer*/
else if(nbval[i]<0) {
/* printf("dans l eau: %d\n",i);*/
fprintf(out,"%c",0);
fprintf(out2,"%c",0);
}
}
fclose(out);
fclose(out2);
printf("il ya %d pixel sans valeur dans %s\n",tmmp,filename);
return 1;
}


int CreateImgEcart(double *sommeHH,double *nbvalHH,double *sommeVV,double *nbvalVV,parameter p,int *lignecol,char *filename,char *FecartypeVV,char *FecartypeHH,double **incidenceEcart)
{
double *ecartypeHH;
double *ecartypeVV;
data mydata;
int x,y;
int i;
FILE *outHH=fopen(FecartypeHH,"w");
FILE *outVV=fopen(FecartypeVV,"w");
FILE *file=fopen(filename,"r");

ecartypeHH=(double*)calloc(lignecol[0]*lignecol[1],sizeof(double));
ecartypeVV=(double*)calloc(lignecol[0]*lignecol[1],sizeof(double));

while(feof(file)==0)
{
/*Pour chaque entrée dans le fichier concatené*/
fread(&mydata, sizeof(data),1,file);

double longi;
/*cette manipulation est necessaire pour passer la longitude dans le meme systeme que la lat*/

if(mydata.lon>p.lon[1])longi=(mydata.lon-MAX_LON);
else longi=mydata.lon;

/*On verifie que l'entrée est dans la zone souhaité*/

if((mydata.lat<=p.lat[1] && mydata.lat>=p.lat[0])&&(longi<=p.lon[1] && longi>=p.lon[0]) )
{
/* printf("surface %f \n",mydata.surface); */
/* calcule des coordonnées en pixel a partir des long/lat de l'entré*/
x=fabs((p.lon[0]-longi)/p.maille);
y=fabs((p.lat[1]-mydata.lat)/p.maille);

/*printf("********************\n");
printf("lat= %f lon=%f\n",mydata.lat,mydata.lon);
printf("sigma0db= %f \n",mydata.sigma0);
printf("sigma0= %f \n",pow(10.0,mydata.sigma0/10));
printf("%d %d \n",x,y);*/
if (mydata.surface==TERRE){
if(mydata.incidence<50){
ecartypeHH[y*lignecol[1]+x]+=pow(pow(10.0,correctS(mydata.sigma0,mydata.sigma0_attn_map,mydata.incidence)/10)-sommeHH[i]/nbvalHH[i],2);
if(mydata.incidence<(*incidenceEcart)[0]) (*incidenceEcart)[0]=mydata.incidence;
else if (mydata.incidence>(*incidenceEcart)[1]) (*incidenceEcart)[1]=mydata.incidence;
}
else{
ecartypeVV[y*lignecol[1]+x]+=pow(pow(10.0,correctS(mydata.sigma0,mydata.sigma0_attn_map,mydata.incidence)/10)-sommeVV[i]/nbvalVV[i],2);
if(mydata.incidence<(*incidenceEcart)[2]) (*incidenceEcart)[2]=mydata.incidence;
else if (mydata.incidence>(*incidenceEcart)[3]) (*incidenceEcart)[3]=mydata.incidence;
}

}
/*Si l'entré n'est pas dans la zone souhaité*/
else{
if(mydata.incidence<50){
ecartypeHH[y*lignecol[1]+x]=-1;
}
else{
ecartypeVV[y*lignecol[1]+x]=-1;
}
}
}

}


for(i=0;i<lignecol[0]*lignecol[1];++i){

if(nbvalHH[i]>0){
double tmp=getPixelValue4(10.0*log10(ecartypeHH[i]/nbvalHH[i]));
unsigned char p=(unsigned char)tmp;
fprintf(outHH,"%c",p);
}
/*Pas de valeur pour ce pixel*/
else if(nbvalHH[i]==0) {
fprintf(outHH,"%c",4);
}

else if(nbvalHH[i]<0) {
/* printf("dans l eau: %d\n",i);*/
fprintf(outHH,"%c",0);
}

if(nbvalVV[i]>0){
double tmp=getPixelValue4(10.0*log10(ecartypeVV[i]/nbvalVV[i]));
unsigned char p=(unsigned char)tmp;
fprintf(outVV,"%c",p);
}
/*Pas de valeur pour ce pixel*/
else if(nbvalVV[i]==0) {
fprintf(outVV,"%c",4);
}
else if(nbvalVV[i]<0) {
/* printf("dans l eau: %d\n",i);*/
fprintf(outVV,"%c",0);
}

}


fclose(file);
fclose(outHH);
fclose(outVV);
free(ecartypeVV);
free(ecartypeHH);
return 1;
}


int main(int argc,char *argv[])
{
double *incidence;
double *incidenceEcart;
incidence=(double*)malloc(sizeof(double)*4);
incidenceEcart=(double*)malloc(sizeof(double)*4);
incidence[0]=INT_MAX;
incidence[1]=INT_MIN;
incidence[2]=INT_MAX;
incidence[3]=INT_MIN;
incidenceEcart[0]=INT_MAX;
incidenceEcart[1]=INT_MIN;
incidenceEcart[2]=INT_MAX;
incidenceEcart[3]=INT_MIN;
parameter p;
int m=getParameter(argc,argv,&p);
if(m==-1) return -1;
int lignecol[2];

lignecol[0]=fabs((p.lat[1]-p.lat[0])/p.maille);
lignecol[1]=fabs((p.lon[1]-p.lon[0])/p.maille);

double *nbvalHH;
double *sommeHH;
double *nbvalVV;
double *sommeVV;

int i,j;
int tail;

tail=strlen(p.appendfilename)+strlen(p.sortie)+1;
char appendfile[tail];
for(i=0;i<tail;++i) appendfile[i]='\0';
strcat(appendfile,p.sortie);
strcat(appendfile,p.appendfilename);

i=strlen(p.appendfilename)-10;
char suffix[i];
for(j=0;j<i;++j) suffix[j]='\0';
strcpy(suffix,p.appendfilename+i);

printf("********************\n");
printf("Chargement des données\n");

int k=-1;
k=readFile(&sommeHH,&nbvalHH,&sommeVV,&nbvalVV,p,lignecol,appendfile,&incidence);
if(k==1) return -1;
printf("********************\n");
printf("Creation de l'image\n");

/* fprintf(stderr,"suffixe:%s\n",suffix);*/

tail=strlen("_imageHH.raw")+strlen(p.sortie)+1+strlen(suffix);
char imageHH[tail];
for(i=0;i<tail;++i) imageHH[i]='\0';
strcat(imageHH,p.sortie);
strcat(imageHH,suffix);
strcat(imageHH,"_imageHH.raw");


tail=strlen("_imageVV.raw")+strlen(p.sortie)+1+strlen(suffix);
char imageVV[tail];
for(i=0;i<tail;++i) imageVV[i]='\0';
strcat(imageVV,p.sortie);
strcat(imageVV,suffix);
strcat(imageVV,"_imageVV.raw");


tail=strlen("_nbvalueHH.raw")+strlen(p.sortie)+1+strlen(suffix);
char nbvalueHH[tail];
for(i=0;i<tail;++i) nbvalueHH[i]='\0';
strcat(nbvalueHH,p.sortie);
strcat(nbvalueHH,suffix);
strcat(nbvalueHH,"_nbvalueHH.raw");


tail=strlen("_nbvalueVV.raw")+strlen(p.sortie)+1+strlen(suffix);
char nbvalueVV[tail];
for(i=0;i<tail;++i) nbvalueVV[i]='\0';
strcat(nbvalueVV,p.sortie);
strcat(nbvalueVV,suffix);
strcat(nbvalueVV,"_nbvalueVV.raw");


tail=strlen("_ecartypeHH.raw")+strlen(p.sortie)+1+strlen(suffix);
char ecartypeHH[tail];
for(i=0;i<tail;++i) ecartypeHH[i]='\0';
strcat(ecartypeHH,p.sortie);
strcat(ecartypeHH,suffix);
strcat(ecartypeHH,"_ecartypeHH.raw");

char ecartypeVV[tail];
for(i=0;i<tail;++i) ecartypeVV[i]='\0';
strcat(ecartypeVV,p.sortie);
strcat(ecartypeVV,suffix);
strcat(ecartypeVV,"_ecartypeVV.raw");

fprintf(stderr,"%s\n%s\n%s\n%s\n",ecartypeHH,ecartypeVV,nbvalueVV,nbvalueHH);

writePicture(sommeHH,nbvalHH,lignecol,imageHH,nbvalueHH);
writePicture(sommeVV,nbvalVV,lignecol,imageVV,nbvalueVV);
CreateImgEcart(sommeHH,nbvalHH,sommeVV,nbvalVV,p,lignecol,appendfile,ecartypeVV,ecartypeHH,&incidenceEcart);

printf(" col=%d ligne=%d\n\n",lignecol[1],lignecol[0]);
printf("Image cree à partir de: %d valeurs\n",k);
/*printf("Incidence minimal: %f\nIncidence maximal:%f\n\n\n",incidence[0],incidence[1]);
*/
tail=strlen("_stat.txt")+strlen(p.sortie)+1+strlen(suffix);
char stat[tail];
for(i=0;i<tail;++i) stat[i]='\0';
strcat(stat,p.sortie);
strcat(stat,suffix);
strcat(stat,"_stat.txt");

Statistiques(imageHH,stat,lignecol,incidence);
Statistiques(imageVV,stat,lignecol,incidence+2);
Statistiques(ecartypeHH,stat,lignecol,incidenceEcart);
Statistiques(ecartypeVV,stat,lignecol,incidenceEcart+2);
//histoStat(nbvalueHH,stat,lignecol);
//histoStat(nbvalueVV,stat,lignecol);

free(nbvalHH);
free(sommeHH);
free(nbvalVV);
free(sommeVV);
free(p.dossier);
free(p.sortie);
free(p.appendfilename);
free(incidence);
free(incidenceEcart);
free(suffix);


return 1;
}

Je l'ai déja "épluché" mais comme je ni connais rien cela ne m'avance pas à grand chose!!

Enfin merci beaucoup....

Merci.

C'est la première fois que je viens sur un forum et je ne pensais quand rentrant de manger j'aurai des réponses!

J'avais déjà lancé
gcc -g -Wall -lm syntheseW.c -o syntheseW
il compile mais j'ai toujours des erreurs de segmentation fault quand je lance mon script....

J'avais également essayé la commande gdb puis lancé run, mais voici la page:
GNU gdb 5.3
Copyright 2002 Free Software Foundation, Inc.
GDB is free software, covered by the GNU General Public License, and you are
welcome to change it and/or distribute copies of it under certain conditions.
Type "show copying" to see the conditions.
There is absolutely no warranty for GDB. Type "show warranty" for details.
This GDB was configured as "i686-pc-linux-gnu"...
(gdb) run
Starting program: /mnt/porto/sde1/vanessa/codes_walid/syntheseW
./syntheseW [date_deb][nb_jours][maille][long_min][long_max][lat_min][lat_max][dossier entre][dossier sortie][appendfile]
Segmentation fault

Je sais pas si c'est une bonne idée ce "run" j'ai trouvé ca sur un site mais toujours est il que ca ne marche pas!

voila.....