Doxygen Source Code Documentation
SUMA_IV_FaceSetsextract.c File Reference
#include "SUMA_suma.h"Go to the source code of this file.
Functions | |
| int * | SUMA_IV_FaceSetsextract (char *IV_filename, int *N_FaceSetList) |
Function Documentation
|
||||||||||||
|
File : Taken from IV_FaceSetsextract.c Author : Ziad Saad Date : Tue Nov 17 19:02:16 CST 1998 Purpose : Extracts the FaceSets coordinates of the nodes in an IV file so that we can manipulate the data to our liking. The program looks for a sequence of characters that seem to indicate the start of the FaceSets list. Once the sequence is found, the series of triplets is read and written out either to a file or to the screen. The sequence of triplets must be terminated by an ending sequence. The starting and ending sequences are hard coded into the program. The program outputs an error message if : If the starting or ending sequences are not found If the number of points read is not a multiple of three If the first number of the first FaceSets triplet and the last character in the starting sequence are not spearated by a space (or tab). You can fix this by manually adding a space. Input paramters : IV_filename (char *) a string specifying the name of the inventor file N_FaceSetList (int *) will give the number of rows in FaceSetList Usage : FaceSetList = SUMA_IV_FaceSetsextract (char *IV_filename, int *N_FaceSetList) Returns : FaceSetList (int *) the facesetlist in the inventor file, an Mx3 integer vector (used to be a matrix before SUMA 1.2) Support : Side effects : Definition at line 59 of file SUMA_IV_FaceSetsextract.c. References i, SUMA_alloc_problem(), SUMA_calloc, SUMA_error_message(), SUMA_free, SUMA_iswordin(), SUMA_malloc, and SUMA_realloc. Referenced by SUMA_Load_Surface_Object_eng().
00060 {/* SUMA_IV_FaceSetsextract */
00061 char s[500],serr[500];
00062 char seq_strt[5][30], seq_end[5][30];
00063 int i, f, ex, si, si_exit, evl, se, se_exit;
00064 int ip, NP, cnt, nospacing, MaxAlloc = 100, *linv, *FaceSetList;
00065 div_t cnt4;
00066 FILE*iv_file;
00067 static char FuncName[]={"SUMA_IV_FaceSetsextract"};
00068
00069
00070 /* intialize the number of points read to 0 */
00071 *N_FaceSetList = 0;
00072
00073 linv = (int *)SUMA_malloc (MaxAlloc*sizeof(int));
00074 if (!linv)
00075 {
00076 SUMA_alloc_problem ("Allocation error in IV-FaceSetExtract");
00077 return (NULL);
00078 }
00079
00080 /*This is the sequence to trigger reading the numbers*/
00081 sprintf (seq_strt[0],"IndexedFaceSet");
00082 sprintf (seq_strt[1],"{");
00083 sprintf (seq_strt[2],"coordIndex");
00084 sprintf (seq_strt[3],"[");
00085
00086 si_exit = 4; /* set equal to the number of strings to be matched */
00087
00088 /*This is a sequence to mark the end of the number list*/
00089 sprintf (seq_end[0],"]");
00090 sprintf (seq_end[1],"}"); /* you do not need that one, it is possible to have other objects following coordIndex */
00091
00092 se_exit = 1; /* set equal to the number of strings to be matched */
00093
00094 iv_file = fopen (IV_filename,"r");
00095 if (iv_file == NULL)
00096 {
00097 SUMA_error_message ("SUMA_IV_FaceSetsextract","Could not open input file ",1);
00098 return (NULL);
00099 }
00100
00101
00102 si = 0;
00103 ex = 1;
00104 cnt = 0;
00105
00106 nospacing = 0; /* this flag is used to when the last number is read and it has the first seq_end*/
00107 /* character attached to it, ie no space ! in between*/
00108 while (ex != EOF && si < si_exit)
00109 {
00110 ex = fscanf (iv_file,"%s",s);
00111
00112 /*evl = equal_strings (s,seq_strt[si]);*/
00113
00114 if (strlen (seq_strt[si]) >= strlen (s))
00115 {
00116 evl = SUMA_iswordin (seq_strt[si],s);
00117 if (evl == 1)
00118 nospacing = 0; /* There is a space between character in starting sequence and first number*/
00119 }
00120 else
00121 {
00122 evl = SUMA_iswordin (s,seq_strt[si]);
00123 if (evl == 1)
00124 nospacing = 1;
00125 }
00126
00127 switch (evl)
00128 {
00129 case 0:
00130 si = 0; /* No match, reset the sequence counter */
00131 break;
00132 case 1:
00133 if (nospacing == 1 && si == si_exit-1) /* it has to be the last character in the sequence*/
00134 {
00135 sprintf (serr,"Must have a space character between first number and last character in start sequence ( %s )",s);
00136 SUMA_error_message ("SUMA_IV_FaceSetsextract",serr,1);
00137 exit (1);
00138
00139 }
00140 si = si +1; /* increment the start sequence counter */
00141 #ifdef DEBUG_3
00142 printf ("found %s --- si = %d\n",s,si);
00143 #endif
00144 break;
00145 default:
00146 break;
00147 }
00148 }
00149
00150 /* Now, read the series of numbers until you encounter the first string of the ending sequence*/
00151 se = 0;
00152 nospacing = 0;
00153
00154 while (ex != EOF && se < se_exit )
00155 {
00156 ex = fscanf (iv_file,"%s",s);
00157 /*evl = equal_strings (s,seq_end[se]);*/
00158
00159 if (strlen (seq_end[se]) >= strlen (s))
00160 {
00161 evl = SUMA_iswordin (seq_end[se],s);
00162 if (evl == 1)
00163 nospacing = 0; /* There is a space between last number and fisrt character in ending sequence */
00164 }
00165 else
00166 {
00167 evl = SUMA_iswordin (s,seq_end[se]);
00168 if (evl == 1)
00169 nospacing = 1;
00170 }
00171
00172 switch (evl)
00173 {
00174 case 0:
00175 f = atoi (s);
00176 linv[cnt] = f;
00177 #ifdef DEBUG_3
00178 printf ("\nNumber (%d): %d is appended to end sequence !\n",cnt,linv[cnt]);
00179 #endif
00180
00181 ++cnt;
00182 if (cnt >= MaxAlloc - 1)
00183 {
00184 MaxAlloc = MaxAlloc + 100;
00185 linv = (int *)SUMA_realloc ((void*) linv, MaxAlloc * sizeof(int));
00186 if (!linv)
00187 {
00188 SUMA_alloc_problem ("Re-Allocation error in IV-FaceSetExtract");
00189 return (NULL);
00190 }
00191
00192 }
00193 se = 0; /* no match for ending sequence, reset the counter */
00194 break;
00195 case 1: /* found the first character in the end sequence */
00196 if (nospacing == 1 && se == 0) /* it has to be the first character in the sequence*/
00197 {
00198 f = atoi (s);
00199 linv[cnt] = f;
00200
00201 ++cnt;
00202 #ifdef DEBUG_3
00203 printf ("\nLast number (%d): %f is appended to end sequence !\n",cnt,f);
00204 #endif
00205 }
00206 #ifdef DEBUG_3
00207 printf ("\nfound %s --- se = %d\n",s,se);
00208 #endif
00209 se = se + 1;
00210 break;
00211 default:
00212 break;
00213 }
00214
00215 }
00216 if (si < si_exit)
00217 {
00218 SUMA_error_message ("IV_FaceSetextract","Could not find specified starting sequence",0);
00219 for (i=0;i<si_exit;++i)
00220 printf ("%s \t",seq_strt[i]);
00221 }
00222 else
00223 {
00224 if (se < se_exit)
00225 {
00226 SUMA_error_message ("IV_FaceSetextract","Could not find specified ending sequence",0);
00227 for (i=0;i<se_exit;++i)
00228 printf ("%s \t",seq_end[i]);
00229 }
00230 else
00231 {
00232 /* check that the number of points read is multiple of 4 */
00233 cnt4 = div (cnt,4);
00234 if (cnt4.rem != 0)
00235 {
00236 SUMA_error_message ("IV_FaceSetextract","number of points read is not multiple of 4 !",0);
00237 #ifdef DEBUG_3
00238 printf ("%f FaceSets sets read !!!\n",(float)cnt/4);
00239 #endif
00240 }
00241
00242 }
00243 }
00244
00245 *N_FaceSetList = cnt4.quot ;
00246
00247 /* Now allocate space for SUMA_IV_FaceSetsextract and fill it up */
00248 NP = 3;
00249 FaceSetList = (int *) SUMA_calloc (*N_FaceSetList * NP, sizeof(int));
00250 if (!FaceSetList)
00251 {
00252 SUMA_alloc_problem("IV_FaceSetextract : Could not allocate");
00253 return(NULL);
00254 }
00255
00256 i = 0;
00257 ip = 0;
00258 while (i < cnt) {
00259 FaceSetList[ip] = linv[i]; ++ip; ++i;
00260 FaceSetList[ip] = linv[i]; ++ip; ++i;
00261 FaceSetList[ip] = linv[i]; ++ip; ++i;
00262 ++i; /* skip the 4th value */
00263 }
00264
00265 fclose (iv_file);
00266
00267 SUMA_free(linv);
00268
00269 return (FaceSetList);
00270
00271 }/* SUMA_IV_FaceSetsextract */
|
