Bonjour,
Je ne sais pas progrmmer en C et je dois utiliser un programme déjà fait.
Je dois compiler mon programme avec :
gcc -Wall -lm syntheseW.c -o syntheseW
La compilation se fait mais ensuite (mon but est d'obtenir des images) lorsque j'utilise un script qui utilise l'executable syntheseW des problèmes interviennent dans mes images:
(...)
Creation de l'image
/mnt/porto/sde1/vanessa/fichier_imageW/le2004_1_6_ecartypeHH.raw
/mnt/porto/sde1/vanessa/fichier_imageW/le2004_1_6_ecartypeVV.raw
/mnt/porto/sde1/vanessa/fichier_imageW/le2004_1_6_nbvalueVV.raw
/mnt/porto/sde1/vanessa/fichier_imageW/le2004_1_6_nbvalueHH.raw
il ya 7542 pixel sans valeur dans /mnt/porto/sde1/vanessa/fichier_imageW/le2004_1_6_imageHH.raw
il ya 7550 pixel sans valeur dans /mnt/porto/sde1/vanessa/fichier_imageW/le2004_1_6_imageVV.raw
./MainScript: line 73: 25218 Segmentation fault (core dumped) $dir$sep$synthese $1 $2 $3 $4 $5 $6 $7 $8 $9 $appendfile
Si je compile en remplaçant lm par g , la compilation ne se fait pas et une page d'erreur apparait :
/home/stage/vfoitier/tmp/cczTnV7g.o: dans la fonction « correctS »:
/mnt/porto/sde1/vanessa/codes_walid/syntheseW.c:71: référence indéfinie vers « cos »
/home/stage/vfoitier/tmp/cczTnV7g.o: dans la fonction « Statistiques »:
/mnt/porto/sde1/vanessa/codes_walid/syntheseW.c:211: référence indéfinie vers « pow »
/mnt/porto/sde1/vanessa/codes_walid/syntheseW.c:215: référence indéfinie vers « sqrt »
/home/stage/vfoitier/tmp/cczTnV7g.o: dans la fonction « readFile »:
/mnt/porto/sde1/vanessa/codes_walid/syntheseW.c:303: référence indéfinie vers « pow »
/mnt/porto/sde1/vanessa/codes_walid/syntheseW.c:309: référence indéfinie vers « pow »
/home/stage/vfoitier/tmp/cczTnV7g.o: dans la fonction « writePicture »:
/mnt/porto/sde1/vanessa/codes_walid/syntheseW.c:347: référence indéfinie vers « log10 »
/home/stage/vfoitier/tmp/cczTnV7g.o: dans la fonction « CreateImgEcart »:
/mnt/porto/sde1/vanessa/codes_walid/syntheseW.c:415: référence indéfinie vers « pow »
/mnt/porto/sde1/vanessa/codes_walid/syntheseW.c:415: référence indéfinie vers « pow »
/mnt/porto/sde1/vanessa/codes_walid/syntheseW.c:420: référence indéfinie vers « pow »
/mnt/porto/sde1/vanessa/codes_walid/syntheseW.c:420: référence indéfinie vers « pow »
/mnt/porto/sde1/vanessa/codes_walid/syntheseW.c:443: référence indéfinie vers « log10 »
/mnt/porto/sde1/vanessa/codes_walid/syntheseW.c:458: référence indéfinie vers « log10 »
collect2: ld a retourné 1 code d'état d'exécution
Cela veut il dire que le programme syntheseW.c serait faut?
Quelle est la signification du "-lm"?
J'espère qu ej'ai été assez clair, je débute et je ne connais pas grand chose à la programmation......
Merci d'avance...
Forum Programmation.c erreur de compilation
20
juil.
2006
# Compilation
Posté par JaguarWan . Évalué à 4.
"-g" permet de compiler avec les informations de debug, c'est un bon réflexe pour trouver ce qui ne va pas, mais il faut l'ajouter sans enlever les autres paramètres:
gcc -g -Wall -lm syntheseW.c -o syntheseW
Par contre, pour voir ce qui ne va pas il faut ensuite lancer le programme dans un debugger, gdb en l'occurence.
gdb ./syntheseW
Cependant, même en trouvant où le programme plante, corriger l'erreur ne sera peut être pas trivial.
Bonne chance :)
[^] # Re: Compilation
Posté par vanessa15 . Évalué à 2.
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.....
[^] # Re: Compilation
Posté par JaguarWan . Évalué à 2.
Pour voir où il plante dans gdb, il faut utiliser la commande bt (backtrace), qui affichera la liste des fonctions appelées et la ligne de code fautive.
Si tu postes le code source complet (avec /mnt/porto/sdb1/otig/frison/lib/hdf4/hdf.h) et le fichier raw qui le fait planter, je pourrais jeter un coup d'oeil.
[^] # Re: Compilation
Posté par vanessa15 . Évalué à 2.
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.....
[^] # Re: Compilation
Posté par JaguarWan . Évalué à 2.
jaguarwan@selenyx:~$ gcc -g syntheseW.c -lm -o syntheseW
jaguarwan@selenyx:~$ ./syntheseW
./syntheseW [date_deb][nb_jours][maille][long_min][long_max][lat_min][lat_max][dossier entre][dossier sortie][appendfile]
jaguarwan@selenyx:~$
Pourrais-tu m'indiquer les paramètres avec lequel tu lance le programme ? En renommant un fichier vide quelconque en écarttype.raw, je devrait pouvoir reproduire le crash :)
[^] # Re: Compilation
Posté par vanessa15 . Évalué à 1.
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.....
[^] # Re: Compilation
Posté par JaguarWan . Évalué à 2.
jaguarwan@selenyx:~/input$ ls
le2004_1_6_ecartypeHH.raw le2004_1_6_nbvalueHH.raw
le2004_1_6_ecartypeVV.raw le2004_1_6_nbvalueVV.raw
J'ai aussi créer un dossier output, que j'ai laissé vide, et j'ai ensuite appelé le programme comme cela:
(gdb) run 1 15 0.125 -25 20 0 25 input output sortie
Starting program: /home/jaguarwan/syntheseW 1 15 0.125 -25 20 0 25 input output sortie
********************
Chargement des données
Program exited with code 0377.
(gdb)
Apparemment il faut des fichiers raw qui contiennent quelque chose de correct... L'idéal serait que tu mettes en ligne les fichiers qui produisent le crash, ou que tu me donnes une méthode pour en générer :)
[^] # Re: Compilation
Posté par vanessa15 . Évalué à 1.
J'essayerai demain,
Je suis chez moi et je n'ai pas les fichiers........
[^] # Re: Compilation
Posté par JaguarWan . Évalué à 3.
J'ai supposé peut être abusivement que i dans ce contexte devait faire référence aux coordonnées du pixel courant, il faudrait donc vérifier ça auprès du programmeur original :)
Il y avait d'autres petites erreurs liées à la mémoire que j'ai corrigées, on ne sait jamais. Je n'ai pas touché au reste.
J'ai mis en ligne le source modifié ici:
http://mammouth.tuxfamily.org/syntheseW.c
et un patch indiquant toutes les modifications là:
http://mammouth.tuxfamily.org/syntheseW.diff
[^] # Re: Compilation
Posté par vanessa15 . Évalué à 1.
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...
[^] # Re: Compilation
Posté par vanessa15 . Évalué à 1.
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.
# Pourquoi compiler?
Posté par Sixel . Évalué à 2.
"Il faut" (Ezekiel 18:4) "forniquer" (Corinthiens 6:9, 10) "avec des chiens" (Thessaloniciens 1:6-9) "morts" (Timothée 3:1-10).
[^] # Re: Pourquoi compiler?
Posté par vanessa15 . Évalué à 2.
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....
[^] # Re: Pourquoi compiler?
Posté par Sixel . Évalué à 1.
- une compilation incorrecte
- des paramètres incorrects
En supposant que l'auteur sache compiler correctement son code source (et que celui-ci ne soit pas buggé), le binaire qu'il te fournirait ne devrait planter que si tes paramètres sont incorrects.
Si c'est le cas, change tes paramètres.
Sinon, c'était ta compilation qui posait problème.
"Il faut" (Ezekiel 18:4) "forniquer" (Corinthiens 6:9, 10) "avec des chiens" (Thessaloniciens 1:6-9) "morts" (Timothée 3:1-10).
Suivre le flux des commentaires
Note : les commentaires appartiennent à celles et ceux qui les ont postés. Nous n’en sommes pas responsables.