12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788899091929394959697989910010110210310410510610710810911011111211311411511611711811912012112212312412512612712812913013113213313413513613713813914014114214314414514614714814915015115215315415515615715815916016116216316416516616716816917017117217317417517617717817918018118218318418518618718818919019119219319419519619719819920020120220320420520620720820921021121221321421521621721821922022122222322422522622722822923023123223323423523623723823924024124224324424524624724824925025125225325425525625725825926026126226326426526626726826927027127227327427527627727827928028128228328428528628728828929029129229329429529629729829930030130230330430530630730830931031131231331431531631731831932032132232332432532632732832933033133233333433533633733833934034134234334434534634734834935035135235335435535635735835936036136236336436536636736836937037137237337437537637737837938038138238338438538638738838939039139239339439539639739839940040140240340440540640740840941041141241341441541641741841942042142242342442542642742842943043143243343443543643743843944044144244344444544644744844945045145245345445545645745845946046146246346446546646746846947047147247347447547647747847948048148248348448548648748848949049149249349449549649749849950050150250350450550650750850951051151251351451551651751851952052152252352452552652752852953053153253353453553653753853954054154254354454554654754854955055155255355455555655755855956056156256356456556656756856957057157257357457557657757857958058158258358458558658758858959059159259359459559659759859960060160260360460560660760860961061161261361461561661761861962062162262362462562662762862963063163263363463563663763863964064164264364464564664764864965065165265365465565665765865966066166266366466566666766866967067167267367467567667767867968068168268368468568668768868969069169269369469569669769869970070170270370470570670770870971071171271371471571671771871972072172272372472572672772872973073173273373473573673773873974074174274374474574674774874975075175275375475575675775875976076176276376476576676776876977077177277377477577677777877978078178278378478578678778878979079179279379479579679779879980080180280380480580680780880981081181281381481581681781881982082182282382482582682782882983083183283383483583683783883984084184284384484584684784884985085185285385485585685785885986086186286386486586686786886987087187287387487587687787887988088188288388488588688788888989089189289389489589689789889990090190290390490590690790890991091191291391491591691791891992092192292392492592692792892993093193293393493593693793893994094194294394494594694794894995095195295395495595695795895996096196296396496596696796896997097197297397497597697797897998098198298398498598698798898999099199299399499599699799899910001001100210031004100510061007100810091010101110121013101410151016101710181019102010211022102310241025102610271028102910301031103210331034103510361037103810391040104110421043104410451046104710481049105010511052105310541055105610571058105910601061106210631064106510661067106810691070107110721073107410751076107710781079108010811082108310841085108610871088108910901091109210931094109510961097109810991100110111021103110411051106110711081109111011111112111311141115111611171118111911201121112211231124112511261127112811291130113111321133113411351136113711381139114011411142114311441145114611471148114911501151115211531154115511561157115811591160116111621163116411651166116711681169117011711172117311741175117611771178117911801181118211831184118511861187118811891190119111921193119411951196119711981199120012011202120312041205120612071208120912101211121212131214121512161217121812191220122112221223122412251226122712281229123012311232123312341235123612371238123912401241124212431244124512461247124812491250125112521253125412551256125712581259126012611262126312641265126612671268126912701271127212731274127512761277127812791280128112821283128412851286128712881289129012911292129312941295129612971298129913001301130213031304130513061307130813091310131113121313131413151316131713181319132013211322132313241325132613271328132913301331133213331334133513361337133813391340134113421343134413451346134713481349135013511352135313541355135613571358135913601361136213631364136513661367136813691370137113721373137413751376137713781379138013811382138313841385138613871388138913901391139213931394139513961397139813991400140114021403140414051406140714081409141014111412141314141415141614171418141914201421142214231424142514261427142814291430143114321433143414351436143714381439144014411442144314441445144614471448144914501451145214531454145514561457145814591460146114621463146414651466146714681469147014711472147314741475147614771478147914801481148214831484148514861487148814891490149114921493149414951496149714981499150015011502150315041505150615071508150915101511151215131514151515161517151815191520152115221523152415251526152715281529153015311532153315341535153615371538153915401541154215431544154515461547154815491550155115521553155415551556155715581559156015611562156315641565156615671568156915701571157215731574157515761577157815791580158115821583158415851586158715881589159015911592159315941595159615971598159916001601160216031604160516061607160816091610161116121613161416151616161716181619162016211622162316241625162616271628162916301631163216331634163516361637163816391640164116421643164416451646164716481649165016511652165316541655165616571658165916601661166216631664166516661667166816691670167116721673167416751676167716781679168016811682168316841685168616871688168916901691169216931694169516961697169816991700170117021703170417051706170717081709171017111712171317141715171617171718171917201721172217231724172517261727172817291730173117321733173417351736173717381739174017411742174317441745174617471748174917501751175217531754175517561757175817591760176117621763176417651766176717681769177017711772177317741775177617771778177917801781178217831784178517861787178817891790179117921793179417951796179717981799180018011802180318041805180618071808180918101811181218131814181518161817181818191820182118221823182418251826182718281829183018311832183318341835183618371838183918401841184218431844184518461847184818491850185118521853185418551856185718581859186018611862186318641865186618671868186918701871187218731874187518761877187818791880188118821883188418851886188718881889189018911892189318941895189618971898189919001901190219031904190519061907190819091910191119121913191419151916191719181919192019211922192319241925192619271928192919301931193219331934193519361937193819391940194119421943194419451946194719481949195019511952195319541955195619571958195919601961196219631964196519661967196819691970197119721973197419751976197719781979198019811982198319841985198619871988198919901991199219931994199519961997199819992000200120022003200420052006200720082009201020112012201320142015201620172018201920202021202220232024202520262027202820292030203120322033203420352036203720382039204020412042204320442045204620472048204920502051205220532054205520562057205820592060206120622063206420652066206720682069207020712072207320742075207620772078207920802081208220832084208520862087208820892090209120922093209420952096209720982099210021012102210321042105210621072108210921102111211221132114211521162117211821192120212121222123212421252126212721282129213021312132213321342135213621372138213921402141214221432144214521462147214821492150215121522153215421552156215721582159216021612162216321642165216621672168216921702171217221732174217521762177217821792180218121822183218421852186218721882189219021912192219321942195219621972198219922002201220222032204220522062207220822092210221122122213221422152216221722182219222022212222222322242225222622272228222922302231223222332234223522362237223822392240224122422243224422452246224722482249225022512252225322542255225622572258225922602261226222632264226522662267226822692270227122722273227422752276227722782279228022812282228322842285228622872288228922902291229222932294229522962297229822992300230123022303230423052306230723082309231023112312231323142315231623172318231923202321232223232324232523262327232823292330233123322333233423352336233723382339234023412342234323442345234623472348234923502351235223532354235523562357235823592360236123622363236423652366236723682369237023712372237323742375237623772378237923802381238223832384238523862387238823892390239123922393239423952396239723982399240024012402240324042405240624072408240924102411241224132414241524162417241824192420242124222423242424252426242724282429243024312432243324342435243624372438243924402441244224432444244524462447244824492450245124522453245424552456245724582459246024612462246324642465246624672468246924702471247224732474247524762477247824792480248124822483248424852486248724882489249024912492249324942495249624972498249925002501250225032504250525062507250825092510251125122513251425152516251725182519252025212522252325242525252625272528252925302531253225332534253525362537253825392540254125422543254425452546254725482549255025512552255325542555255625572558255925602561256225632564256525662567256825692570257125722573257425752576257725782579258025812582258325842585258625872588258925902591259225932594259525962597259825992600260126022603260426052606260726082609261026112612261326142615261626172618261926202621262226232624262526262627262826292630263126322633263426352636263726382639264026412642264326442645264626472648264926502651265226532654265526562657265826592660266126622663266426652666266726682669267026712672267326742675267626772678267926802681268226832684268526862687268826892690269126922693269426952696269726982699270027012702270327042705270627072708270927102711271227132714271527162717271827192720272127222723272427252726272727282729273027312732273327342735273627372738273927402741274227432744274527462747274827492750275127522753275427552756275727582759276027612762276327642765276627672768276927702771277227732774277527762777277827792780278127822783278427852786278727882789279027912792279327942795279627972798279928002801280228032804280528062807280828092810281128122813281428152816281728182819282028212822282328242825282628272828282928302831283228332834283528362837283828392840284128422843284428452846284728482849285028512852285328542855285628572858285928602861286228632864286528662867286828692870287128722873287428752876287728782879288028812882288328842885288628872888288928902891289228932894289528962897289828992900290129022903290429052906290729082909291029112912291329142915291629172918291929202921292229232924292529262927292829292930293129322933293429352936293729382939294029412942294329442945294629472948294929502951295229532954295529562957295829592960296129622963296429652966296729682969297029712972297329742975297629772978297929802981298229832984298529862987298829892990299129922993299429952996299729982999300030013002300330043005300630073008300930103011301230133014301530163017301830193020302130223023302430253026302730283029303030313032303330343035303630373038303930403041304230433044304530463047304830493050305130523053305430553056305730583059306030613062306330643065306630673068306930703071307230733074307530763077307830793080308130823083308430853086308730883089309030913092309330943095309630973098309931003101310231033104310531063107310831093110311131123113311431153116311731183119312031213122312331243125312631273128312931303131313231333134313531363137313831393140314131423143314431453146314731483149315031513152315331543155315631573158315931603161316231633164316531663167316831693170317131723173317431753176317731783179318031813182318331843185318631873188318931903191319231933194319531963197319831993200320132023203320432053206320732083209321032113212321332143215321632173218321932203221322232233224322532263227322832293230323132323233323432353236323732383239324032413242324332443245324632473248324932503251325232533254325532563257325832593260326132623263326432653266326732683269327032713272327332743275327632773278327932803281328232833284328532863287328832893290329132923293329432953296329732983299330033013302330333043305330633073308330933103311331233133314331533163317331833193320332133223323332433253326332733283329333033313332333333343335333633373338333933403341334233433344334533463347334833493350335133523353335433553356335733583359336033613362336333643365336633673368336933703371337233733374337533763377337833793380338133823383338433853386338733883389339033913392339333943395339633973398339934003401340234033404340534063407340834093410341134123413341434153416341734183419342034213422342334243425342634273428342934303431343234333434343534363437343834393440344134423443344434453446344734483449345034513452345334543455345634573458345934603461346234633464346534663467346834693470347134723473347434753476347734783479348034813482348334843485348634873488348934903491349234933494349534963497349834993500350135023503350435053506350735083509351035113512351335143515351635173518351935203521352235233524352535263527352835293530353135323533353435353536353735383539354035413542354335443545354635473548354935503551355235533554355535563557355835593560356135623563356435653566356735683569357035713572357335743575357635773578357935803581358235833584358535863587358835893590359135923593359435953596359735983599360036013602360336043605360636073608360936103611361236133614361536163617361836193620362136223623362436253626362736283629363036313632363336343635363636373638363936403641364236433644364536463647364836493650365136523653365436553656365736583659366036613662366336643665366636673668366936703671367236733674367536763677367836793680368136823683368436853686368736883689369036913692369336943695369636973698369937003701370237033704370537063707370837093710371137123713371437153716371737183719372037213722372337243725372637273728372937303731373237333734373537363737373837393740374137423743374437453746374737483749375037513752375337543755375637573758375937603761376237633764376537663767376837693770377137723773377437753776377737783779378037813782378337843785378637873788378937903791379237933794379537963797379837993800380138023803380438053806380738083809381038113812381338143815381638173818381938203821382238233824382538263827382838293830383138323833383438353836383738383839384038413842384338443845384638473848384938503851385238533854385538563857385838593860386138623863386438653866386738683869387038713872387338743875387638773878387938803881388238833884388538863887388838893890389138923893389438953896389738983899390039013902390339043905390639073908390939103911391239133914391539163917391839193920392139223923392439253926392739283929393039313932393339343935393639373938393939403941394239433944394539463947394839493950395139523953 |
- /*************************************/
- /* Auteur : Rémi Synave */
- /* Date de création : 10/03/10 */
- /* Date de modification : 10/12/14 */
- /* Version : 0.3 */
- /*************************************/
- /***************************************************************************/
- /* This file is part of a2ri. */
- /* */
- /* a2ri is free software: you can redistribute it and/or modify it */
- /* under the terms of the GNU Lesser General Public License as published */
- /* by the Free Software Foundation, either version 3 of the License, or */
- /* (at your option) any later version. */
- /* */
- /* a2ri is distributed in the hope that it will be useful, */
- /* but WITHOUT ANY WARRANTY; without even the implied warranty of */
- /* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the */
- /* GNU Lesser General Public License for more details. */
- /* */
- /* You should have received a copy of the GNU Lesser General Public */
- /* License along with a2ri. */
- /* If not, see <http://www.gnu.org/licenses/>. */
- /***************************************************************************/
- #include "io.h"
- /********** INTERMEDIATE TYPES AND FUNCTIONS **********/
- /* Les fonctions intermédiaires sont préfixées de IF */
- /* et les types intermédiaires de IT */
- void
- IFaffichage_chargement (
- FILE * is,
- long posmax)
- {
- long poscur;
- poscur = ftell (is);
- printf ("%c%c%c%c%c%c%c%c%c%c%c%c%c%lf %c", 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8,
- 8, 8, ((poscur) * 100.0) / (posmax * 1.0), 37);
- }
- //fonction qui ecrit les aretes dans un fichier vef
- void IFsave_vef_write_edge (
- vf_model * m,
- hashtable * table,
- FILE * os)
- {
- int p1,
- p2;
- int indexarete = 0;
-
- //pour toute les faces, on regarde si l'arete existe deja dans la table de hachage
- //si oui, on ne fait rien
- //sinon on entre l'arete dans la table et on l'ecrit dans le fichier.
- //on enregistre également son index pour l'écriture des faces définies par trois aretes
- for (int i = 0; i < m->nbface; i++)
- {
- p1 = m->fa[i].ve1;
- p2 = m->fa[i].ve2;
- vf_edge *edgetemp = hashtable_look_for (table, p1, p2);
- if (edgetemp == NULL)
- {
- edgetemp = (vf_edge *) malloc (sizeof (vf_edge));
- a2ri_erreur_critique_si (edgetemp == NULL,
- "erreur allocation memoire pour edgetemp\nIFsave_vef_write_edge");
- edgetemp->ve1 = p1;
- edgetemp->ve2 = p2;
- edgetemp->att_int = indexarete;
- indexarete++;
- fprintf (os, "%d %d\n", p1, p2);
- hashtable_add (table, edgetemp);
- }
-
- p1 = m->fa[i].ve2;
- p2 = m->fa[i].ve3;
- edgetemp = hashtable_look_for (table, p1, p2);
- if (edgetemp == NULL)
- {
- edgetemp = (vf_edge *) malloc (sizeof (vf_edge));
- a2ri_erreur_critique_si (edgetemp == NULL,
- "erreur allocation memoire pour edgetemp\nIFsave_vef_write_edge");
- edgetemp->ve1 = p1;
- edgetemp->ve2 = p2;
- edgetemp->att_int = indexarete;
- indexarete++;
- fprintf (os, "%d %d\n", p1, p2);
- hashtable_add (table, edgetemp);
- }
-
- p1 = m->fa[i].ve3;
- p2 = m->fa[i].ve1;
- edgetemp = hashtable_look_for (table, p1, p2);
- if (edgetemp == NULL)
- {
- edgetemp = (vf_edge *) malloc (sizeof (vf_edge));
- a2ri_erreur_critique_si (edgetemp == NULL,
- "erreur allocation memoire pour edgetemp\nIFsave_vef_write_edge");
- edgetemp->ve1 = p1;
- edgetemp->ve2 = p2;
- edgetemp->att_int = indexarete;
- indexarete++;
- fprintf (os, "%d %d\n", p1, p2);
- hashtable_add (table, edgetemp);
- }
- }
- }
- //fonction qui ecrit les aretes dans un fichier gts
- void IFsave_gts_write_edge (
- vf_model * m,
- hashtable * table,
- FILE * os)
- {
- int p1,
- p2;
- int indexarete = 0;
-
- //pour toute les faces, on regarde si l'arete existe deja dans la table de hachage
- //si oui, on ne fait rien
- //sinon on entre l'arete dans la table et on l'ecrit dans le fichier.
- //on enregistre également son index pour l'écriture des faces définies par trois aretes
- for (int i = 0; i < m->nbface; i++)
- {
- p1 = m->fa[i].ve1;
- p2 = m->fa[i].ve2;
- vf_edge *edgetemp = hashtable_look_for (table, p1, p2);
- if (edgetemp == NULL)
- {
- edgetemp = (vf_edge *) malloc (sizeof (vf_edge));
- a2ri_erreur_critique_si (edgetemp == NULL,
- "erreur allocation memoire pour edgetemp\nIFsave_gts_write_edge");
- edgetemp->ve1 = p1;
- edgetemp->ve2 = p2;
- edgetemp->att_int = indexarete;
- indexarete++;
- fprintf (os, "%d %d\n", p1 + 1, p2 + 1);
- hashtable_add (table, edgetemp);
- }
-
- p1 = m->fa[i].ve2;
- p2 = m->fa[i].ve3;
- edgetemp = hashtable_look_for (table, p1, p2);
- if (edgetemp == NULL)
- {
- edgetemp = (vf_edge *) malloc (sizeof (vf_edge));
- a2ri_erreur_critique_si (edgetemp == NULL,
- "erreur allocation memoire pour edgetemp\nIFsave_gts_write_edge");
- edgetemp->ve1 = p1;
- edgetemp->ve2 = p2;
- edgetemp->att_int = indexarete;
- indexarete++;
- fprintf (os, "%d %d\n", p1 + 1, p2 + 1);
- hashtable_add (table, edgetemp);
- }
-
- p1 = m->fa[i].ve3;
- p2 = m->fa[i].ve1;
- edgetemp = hashtable_look_for (table, p1, p2);
- if (edgetemp == NULL)
- {
- edgetemp = (vf_edge *) malloc (sizeof (vf_edge));
- a2ri_erreur_critique_si (edgetemp == NULL,
- "erreur allocation memoire pour edgetemp\nIFsave_gts_write_edge");
- edgetemp->ve1 = p1;
- edgetemp->ve2 = p2;
- edgetemp->att_int = indexarete;
- indexarete++;
- fprintf (os, "%d %d\n", p1 + 1, p2 + 1);
- hashtable_add (table, edgetemp);
- }
- }
- }
- /********** MAIN FUNCTIONS **********/
- /**
- Ouverture d'un fichier contenant un modèle
- @param file chemin du fichier à ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_open_file (
- char *file,
- vf_model * m)
- {
- char ext[5];
-
- memcpy (ext, file + (strlen (file) - 4), 4);
- ext[4] = 0;
- //appel de la bonne fonction en fonction de l'extension du fichier
- if (!strcmp (ext, ".wrl"))
- return (a2ri_vf_load_wrl (file, m));
- if (!strcmp (ext, ".asc"))
- return (a2ri_vf_load_asc (file, m));
- if (!strcmp (ext, ".vef"))
- return (a2ri_vf_load_vef (file, m));
- if (!strcmp (ext, ".pgn"))
- return (a2ri_vf_load_pgn (file, m));
- if (!strcmp (ext, ".ply"))
- return (a2ri_vf_load_ply (file, m));
- if (!strcmp (ext, ".off"))
- return (a2ri_vf_load_off (file, m));
- if (!strcmp (ext, ".obj"))
- return (a2ri_vf_load_obj (file, m));
- if (!strcmp (ext, ".gts"))
- return (a2ri_vf_load_gts (file, m));
- if (!strcmp (ext, ".stl"))
- return (a2ri_vf_load_stl (file, m));
- return 0;
- }
- /**
- Ouverture d'un fichier asc contenant un nuage de point
- @param file chemin du fichier à ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_load_asc (
- char *file,
- vf_model * m)
- {
- FILE *is = NULL;
- double coordx,
- coordy,
- coordz;
- int i;
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- i=0;
- fscanf (is, "%lf %lf %lf", &coordx, &coordy, &coordz);
- while(!feof(is))
- {
- //fscanf (is, "%lf %lf %lf", &coordx, &coordy, &coordz);
- a2ri_vf_add_vertex (m, coordx, coordy, coordz);
- if (i == 0)
- {
- m->xmin = coordx;
- m->xmax = coordx;
- m->ymin = coordy;
- m->ymax = coordy;
- m->zmin = coordz;
- m->zmax = coordz;
- }
- else
- {
- //mis à jour de la bounding box du modèle
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- }
- fscanf (is, "%lf %lf %lf", &coordx, &coordy, &coordz);
- fscanf (is, "%lf %lf %lf", &coordx, &coordy, &coordz);
- i=1;
- }
- fclose (is);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Ouverture d'un fichier pgn contenant un modèle
- @param file chemin du fichier à ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_load_pgn (
- char *file,
- vf_model * m)
- {
- FILE *is = NULL;
- int nbvertex,
- nbface,
- ve1,
- ve2,
- ve3;
- double coordx,
- coordy,
- coordz;
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- //lecture du nombre de sommets et de faces
- fscanf (is, "%d %d", &nbvertex, &nbface);
- //lecture des sommets
- for (int i = 0; i < nbvertex; i++)
- {
- fscanf (is, "%lf %lf %lf", &coordx, &coordy, &coordz);
- a2ri_vf_add_vertex (m, coordx, coordy, coordz);
- if (i == 0)
- {
- m->xmin = coordx;
- m->xmax = coordx;
- m->ymin = coordy;
- m->ymax = coordy;
- m->zmin = coordz;
- m->zmax = coordz;
- }
- else
- {
- //mis à jour de la bounding box du modèle
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- }
- }
- //lecture des faces
- for (int i = 0; i < nbface; i++)
- {
- fscanf (is, "%d %d %d", &ve1, &ve2, &ve3);
- a2ri_vf_add_face (m, ve1, ve2, ve3);
- fscanf (is, "%d", &ve1);
- }
- fclose (is);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Ouverture d'un fichier obj contenant un modèle
- @param file chemin du fichier à ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_load_obj (
- char *file,
- vf_model * m)
- {
- FILE *is = NULL;
- int ve1,
- ve2,
- ve3;
- double coordx,
- coordy,
- coordz;
- char str[250],
- str_ve1[20],
- str_ve2[20],
- str_ve3[20];
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- fscanf (is, "%s", str);
- while (strcmp (str, "v") != 0)
- fscanf (is, "%s", str);
- //lecture des sommets
- while (strcmp (str, "v") == 0 && !feof (is))
- {
- fscanf (is, "%lf %lf %lf", &coordx, &coordy, &coordz);
- //ajout d'un sommet
- a2ri_vf_add_vertex (m, coordx, coordy, coordz);
- fscanf (is, "%s", str);
- //mis à jour de la bounding box
- if (m->nbvertex == 1)
- {
- m->xmin = coordx;
- m->xmax = coordx;
- m->ymin = coordy;
- m->ymax = coordy;
- m->zmin = coordz;
- m->zmax = coordz;
- }
- else
- {
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- }
- }
- while (strcmp (str, "f") != 0 && !feof (is))
- {
- fscanf (is, "%s", str);
- }
- //lecture des faces
- while (strcmp (str, "f") == 0 && !feof (is))
- {
- fscanf (is, "%s %s %s", str_ve1, str_ve2, str_ve3);
- ve1 = atoi (str_ve1);
- ve2 = atoi (str_ve2);
- ve3 = atoi (str_ve3);
- ve1--;
- ve2--;
- ve3--;
- //ajout d'une face
- a2ri_vf_add_face (m, ve1, ve2, ve3);
- fscanf (is, "%s", str);
- if (str[0] >= '0' && str[0] <= '9')
- {
- ve2 = atoi (str);
- ve2--;
- a2ri_vf_add_face (m, ve3, ve2, ve1);
- fscanf (is, "%s", str);
- }
- }
- fclose (is);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Ouverture d'un fichier wrl contenant un modèle
- @param file chemin du fichier à ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_load_wrl1 (
- char *file,
- vf_model * m)
- {
- FILE *is = NULL;
- char str[20];
- char coordx[15],
- coordy[15],
- coordz[15];
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- fscanf (is, "%s", str);
- while (strcmp (str, "point") && strcmp (str, "point["))
- fscanf (is, "%s", str);
- if (!strcmp (str, "point"))
- fscanf (is, "%s", str);
- //lecture des sommets
- fscanf (is, "%s", coordx);
- int i = 0;
- while (strcmp (coordx, "}") && strcmp (coordx, "]"))
- {
- fscanf (is, "%s", coordy);
- fscanf (is, "%s", coordz);
- if (coordz[strlen (coordz) - 1] == ',')
- {
- memcpy (coordz, coordz, strlen (coordz) - 1);
- coordz[strlen (coordz)] = '\0';
- }
- //ajout d'un sommet
- a2ri_vf_add_vertex (m, atof (coordx), atof (coordy), atof (coordz));
- //mis à jour de la bounding box
- if (i == 0)
- {
- i++;
- m->xmin = atof (coordx);
- m->xmax = atof (coordx);
- m->ymin = atof (coordy);
- m->ymax = atof (coordy);
- m->zmin = atof (coordz);
- m->zmax = atof (coordz);
- }
- else
- {
- if (atof (coordx) < m->xmin)
- m->xmin = atof (coordx);
- if (atof (coordx) > m->xmax)
- m->xmax = atof (coordx);
- if (atof (coordy) < m->ymin)
- m->ymin = atof (coordy);
- if (atof (coordy) > m->ymax)
- m->ymax = atof (coordy);
- if (atof (coordz) < m->zmin)
- m->zmin = atof (coordz);
- if (atof (coordz) > m->zmax)
- m->zmax = atof (coordz);
- }
- fscanf (is, "%s", coordx);
- }
- fscanf (is, "%s", str);
- while (strcmp (str, "coordIndex") && strcmp (str, "coordIndex["))
- fscanf (is, "%s", str);
- if (!strcmp (str, "coordIndex"))
- fscanf (is, "%s", str);
- //lecture des faces
- fscanf (is, "%s", coordx);
- while (strcmp (coordx, "}") && strcmp (coordx, "]"))
- {
- fscanf (is, "%s", coordy);
- fscanf (is, "%s", coordz);
- if (coordx[strlen (coordx) - 1] == ',')
- {
- memcpy (coordx, coordx, strlen (coordx) - 1);
- coordz[strlen (coordx)] = '\0';
- }
- if (coordy[strlen (coordy) - 1] == ',')
- {
- memcpy (coordy, coordy, strlen (coordy) - 1);
- coordy[strlen (coordy)] = '\0';
- }
- if (coordz[strlen (coordz) - 1] == ',')
- {
- memcpy (coordz, coordz, strlen (coordz) - 1);
- coordz[strlen (coordz)] = '\0';
- }
- //ajout d'une face
- a2ri_vf_add_face (m, atoi (coordx), atoi (coordy), atoi (coordz));
- // reading the -1
- fscanf (is, "%s", coordy);
- if (coordy[strlen (coordy) - 1] == ',')
- {
- memcpy (coordy, coordy, strlen (coordy) - 1);
- coordy[strlen (coordy)] = '\0';
- }
- if (atoi (coordy) != -1)
- {
- a2ri_vf_add_face (m, atoi (coordz), atoi (coordy), atoi (coordx));
- fscanf (is, "%s", coordx);
- }
- fscanf (is, "%s", coordx);
- }
- fclose (is);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Ouverture d'un fichier wrl contenant un modèle
- @param file chemin du fichier à ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_load_wrl2 (
- char *file,
- vf_model * m)
- {
- FILE *is = NULL;
- char str[20];
- char coordx[15],
- coordy[15],
- coordz[15];
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- fscanf (is, "%s", str);
- while (strcmp (str, "point") && strcmp (str, "point["))
- fscanf (is, "%s", str);
- if (!strcmp (str, "point"))
- fscanf (is, "%s", str);
- //lecture des sommets
- fscanf (is, "%s", coordx);
- int i = 0;
- while (strcmp (coordx, "}") && strcmp (coordx, "]"))
- {
- fscanf (is, "%s", coordy);
- fscanf (is, "%s", coordz);
- if (coordz[strlen (coordz) - 1] == ',')
- {
- memcpy (coordz, coordz, strlen (coordz) - 1);
- coordz[strlen (coordz)] = '\0';
- }
- //ajout d'un sommet
- a2ri_vf_add_vertex (m, atof (coordx), atof (coordy), atof (coordz));
- //mis à jour de la bounding box
- if (i == 0)
- {
- i++;
- m->xmin = atof (coordx);
- m->xmax = atof (coordx);
- m->ymin = atof (coordy);
- m->ymax = atof (coordy);
- m->zmin = atof (coordz);
- m->zmax = atof (coordz);
- }
- else
- {
- if (atof (coordx) < m->xmin)
- m->xmin = atof (coordx);
- if (atof (coordx) > m->xmax)
- m->xmax = atof (coordx);
- if (atof (coordy) < m->ymin)
- m->ymin = atof (coordy);
- if (atof (coordy) > m->ymax)
- m->ymax = atof (coordy);
- if (atof (coordz) < m->zmin)
- m->zmin = atof (coordz);
- if (atof (coordz) > m->zmax)
- m->zmax = atof (coordz);
- }
- fscanf (is, "%s", coordx);
- }
- fscanf (is, "%s", str);
- while (strcmp (str, "coordIndex") && strcmp (str, "coordIndex["))
- fscanf (is, "%s", str);
- if (!strcmp (str, "coordIndex"))
- fscanf (is, "%s", str);
- //lecture des faces
- fscanf (is, "%s", coordx);
- while (strcmp (coordx, "}") && strcmp (coordx, "]"))
- {
- fscanf (is, "%s", coordy);
- fscanf (is, "%s", coordz);
- //ajout d'une face
- a2ri_vf_add_face (m, atoi (coordx), atoi (coordy), atoi (coordz));
- // reading the -1
- fscanf (is, "%s", coordx);
- fscanf (is, "%s", coordx);
- }
- fclose (is);
-
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Ouverture d'un fichier wrl contenant un modèle
- @param file chemin du fichier à ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_load_wrl (
- char *file,
- vf_model * m)
- {
- FILE *is = NULL;
- char str[20];
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- fscanf (is, "%s", str);
- fscanf (is, "%s", str);
- if (!strcmp (str, "v1.0") || !strcmp (str, "V1.0"))
- return a2ri_vf_load_wrl1 (file, m);
- if (!strcmp (str, "v2.0") || !strcmp (str, "V2.0"))
- return a2ri_vf_load_wrl2 (file, m);
- return 0;
- }
- /**
- Ouverture d'un fichier vef contenant un modèle
- @param file chemin du fichier à ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_load_vef (
- char *file,
- vf_model * m)
- {
- FILE *is = NULL;
- int nbvertex,
- nbedge,
- nbface,
- ve1,
- ve2,
- ve3,
- ed1,
- ed2,
- ed3,
- classif;
- double coordx,
- coordy,
- coordz;
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- a2ri_vf_init (m);
- //lecture du nombre de sommets
- fscanf (is, "%d", &nbvertex);
- //lecture du nombre d'aretes
- fscanf (is, "%d", &nbedge);
- //lecture du nombre de faces
- fscanf (is, "%d", &nbface);
- //lecture des sommets
- for (int i = 0; i < nbvertex; i++)
- {
- fscanf (is, "%lf", &coordx);
- fscanf (is, "%lf", &coordy);
- fscanf (is, "%lf", &coordz);
- //ajout d'un sommet
- a2ri_vf_add_vertex (m, coordx, coordy, coordz);
- //mis à jour de la bounding box
- if (i == 0)
- {
- m->xmin = coordx;
- m->xmax = coordx;
- m->ymin = coordy;
- m->ymax = coordy;
- m->zmin = coordz;
- m->zmax = coordz;
- }
- else
- {
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- }
- }
- //lecture des aretes et enregistrement dans un tableau temporaire
- vf_edge *temp = (vf_edge *) malloc (nbedge * sizeof (vf_edge));
- a2ri_erreur_critique_si (temp == NULL,
- "erreur allocation memoire pour temp\na2ri_vf_load_vef");
- for (int i = 0; i < nbedge; i++)
- {
- fscanf (is, "%d", &ve1);
- fscanf (is, "%d", &ve2);
- temp[i].ve1 = ve1;
- temp[i].ve2 = ve2;
- }
- //lecture des faces
- for (int i = 0; i < nbface; i++)
- {
- fscanf (is, "%d", &ed1);
- fscanf (is, "%d", &ed2);
- fscanf (is, "%d", &ed3);
- //recherche des trois sommets contenu dans les trois aretes
- if (temp[ed1].ve1 == temp[ed2].ve1 || temp[ed1].ve1 == temp[ed2].ve2)
- {
- ve1 = temp[ed1].ve1;
- ve3 = temp[ed1].ve2;
- }
- else
- {
- ve1 = temp[ed1].ve2;
- ve3 = temp[ed1].ve1;
- }
- if (temp[ed2].ve1 == temp[ed3].ve1 || temp[ed2].ve1 == temp[ed3].ve2)
- ve2 = temp[ed2].ve1;
- else
- ve2 = temp[ed2].ve2;
- //ajout de la face
- a2ri_vf_add_face (m, ve1, ve2, ve3);
- }
- //lecture des classification si disponible
- for (int i = 0; i < nbface; i++)
- {
- if (!feof (is))
- {
- fscanf (is, "%d", &classif);
- m->fa[i].type = classif;
- }
- }
- fclose (is);
-
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Ouverture d'un fichier gts contenant un modèle
- @param file chemin du fichier à ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_load_gts (
- char *file,
- vf_model * m)
- {
- FILE *is = NULL;
- int nbvertex,
- nbedge,
- nbface,
- ve1,
- ve2,
- ve3,
- ed1,
- ed2,
- ed3,
- classif;
- double coordx,
- coordy,
- coordz;
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- a2ri_vf_init (m);
- //lecture du nombre de sommets
- fscanf (is, "%d", &nbvertex);
- //lecture du nombre d'aretes
- fscanf (is, "%d", &nbedge);
- //lecture du nombre de faces
- fscanf (is, "%d", &nbface);
- //lecture des sommets
- for (int i = 0; i < nbvertex; i++)
- {
- fscanf (is, "%lf", &coordx);
- fscanf (is, "%lf", &coordy);
- fscanf (is, "%lf", &coordz);
- //ajout d'un sommet
- a2ri_vf_add_vertex (m, coordx, coordy, coordz);
- //mis à jour de la bounding box
- if (i == 0)
- {
- m->xmin = coordx;
- m->xmax = coordx;
- m->ymin = coordy;
- m->ymax = coordy;
- m->zmin = coordz;
- m->zmax = coordz;
- }
- else
- {
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- }
- }
- //lecture des aretes et enregistrement dans un tableau temporaire
- vf_edge *temp = (vf_edge *) malloc (nbedge * sizeof (vf_edge));
- a2ri_erreur_critique_si (temp == NULL,
- "erreur allocation memoire pour temp\na2ri_vf_load_gts");
- for (int i = 0; i < nbedge; i++)
- {
- fscanf (is, "%d", &ve1);
- fscanf (is, "%d", &ve2);
- ve1--;
- ve2--;
- temp[i].ve1 = ve1;
- temp[i].ve2 = ve2;
- }
- //lecture des faces
- for (int i = 0; i < nbface; i++)
- {
- fscanf (is, "%d", &ed1);
- fscanf (is, "%d", &ed2);
- fscanf (is, "%d", &ed3);
- ed1--;
- ed2--;
- ed3--;
- //recherche des trois sommets contenu dans les trois aretes
- if (temp[ed1].ve1 == temp[ed2].ve1 || temp[ed1].ve1 == temp[ed2].ve2)
- {
- ve1 = temp[ed1].ve1;
- ve3 = temp[ed1].ve2;
- }
- else
- {
- ve1 = temp[ed1].ve2;
- ve3 = temp[ed1].ve1;
- }
- if (temp[ed2].ve1 == temp[ed3].ve1 || temp[ed2].ve1 == temp[ed3].ve2)
- ve2 = temp[ed2].ve1;
- else
- ve2 = temp[ed2].ve2;
- //ajout de la face
- a2ri_vf_add_face (m, ve1, ve2, ve3);
- }
- //lecture des classification si disponible
- for (int i = 0; i < nbface; i++)
- {
- if (!feof (is))
- {
- fscanf (is, "%d", &classif);
- m->fa[i].type = classif;
- }
- }
- fclose (is);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Ouverture d'un fichier ply contenant un modèle
- @param file chemin du fichier à ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_load_ply (
- char *file,
- vf_model * m)
- {
- int prop = 0,
- nbpropvertex = 0,
- nbpropface = 0,
- nbvertex,
- nbface,
- ve1,
- ve2,
- ve3;
- char str[20];
- double coordx,
- coordy,
- coordz;
- FILE *is = NULL;
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- fscanf (is, "%s", str);
- if (strcmp (str, "ply"))
- return 0;
- while (strcmp (str, "end_header"))
- {
- fscanf (is, "%s", str);
- if (prop == 1 && !strcmp (str, "property"))
- nbpropvertex++;
- if (prop == 2 && !strcmp (str, "property"))
- nbpropface++;
- if (!strcmp (str, "element"))
- {
- fscanf (is, "%s", str);
- //lecture du nombre de sommets
- if (!strcmp (str, "vertex"))
- {
- prop = 1;
- fscanf (is, "%d", &nbvertex);
- }
- //lecture du nombre de face
- else
- {
- if (!strcmp (str, "face"))
- {
- prop = 1;
- fscanf (is, "%d", &nbface);
- }
- }
- }
- }
- //lecture des sommets
- for (int i = 0; i < nbvertex; i++)
- {
- fscanf (is, "%lf %lf %lf", &coordx, &coordy, &coordz);
- //ajout d'un sommet
- a2ri_vf_add_vertex (m, coordx, coordy, coordz);
- //mis à jour de la bounding box
- if (i == 0)
- {
- m->xmin = coordx;
- m->xmax = coordx;
- m->ymin = coordy;
- m->ymax = coordy;
- m->zmin = coordz;
- m->zmax = coordz;
- }
- else
- {
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- }
- for (int j = 0; j < nbpropvertex - 4; j++)
- fscanf (is, "%s", str);
- }
- //lecture des faces
- for (int i = 0; i < nbface; i++)
- {
- fscanf (is, "%d", &ve1);
- fscanf (is, "%d %d %d", &ve1, &ve2, &ve3);
- //ajout d'une face
- a2ri_vf_add_face (m, ve1, ve2, ve3);
- }
- fclose (is);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Ouverture d'un fichier off contenant un modèle
- @param file chemin du fichier à ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_load_off (
- char *file,
- vf_model * m)
- {
- FILE *is = NULL;
- int nbvertex,
- nbface,
- ve1,
- ve2,
- ve3;
- double coordx,
- coordy,
- coordz;
- char str[250];
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- fscanf (is, "%s", str);
- if (strcmp (str, "OFF"))
- return 0;
- //lecture du nombre de sommets et de faces
- fscanf (is, "%d %d %d", &nbvertex, &nbface, &ve1);
- //lecture des sommets
- for (int i = 0; i < nbvertex; i++)
- {
- fscanf (is, "%lf %lf %lf", &coordx, &coordy, &coordz);
- //ajout du sommet
- a2ri_vf_add_vertex (m, coordx, coordy, coordz);
- //mis à jour de la bounding box
- if (i == 0)
- {
- m->xmin = coordx;
- m->xmax = coordx;
- m->ymin = coordy;
- m->ymax = coordy;
- m->zmin = coordz;
- m->zmax = coordz;
- }
- else
- {
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- }
- }
- //lecture des faces
- for (int i = 0; i < nbface; i++)
- {
- fscanf (is, "%d", &ve1);
- fscanf (is, "%d %d %d", &ve1, &ve2, &ve3);
- //ajout d'une face
- a2ri_vf_add_face (m, ve1, ve2, ve3);
- }
- fclose (is);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Ouverture d'un fichier stl contenant un modèle
- @param file chemin du fichier ÃÂ ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_load_stl (
- char *file,
- vf_model * m)
- {
- FILE *is = NULL;
- char str[250];
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- fscanf (is, "%s", str);
- fclose (is);
- if (!strcmp (str, "solid"))
- return a2ri_vf_load_ascii_stl (file, m);
- else
- return a2ri_vf_load_binary_stl (file, m);
- }
- /**
- Ouverture d'un fichier stl binaire contenant un modèle
- @param file chemin du fichier ÃÂ ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_load_binary_stl (
- char *file,
- vf_model * m)
- {
- FILE *is = NULL;
- char str[80];
- unsigned int nbface;
- unsigned char finligne;
- float coordx,
- coordy,
- coordz;
- int ve1 = 0,
- ve2 = 0,
- ve3 = 0;
- point3d ptmin,
- ptmax;
- space_partition sp;
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "rb");
- if (is == NULL)
- return 0;
- fread (&str, sizeof (char), 80, is);
- fread (&nbface, sizeof (unsigned int), 1, is);
- //recherche de la bounding box du modele numerique
- for (unsigned int i = 0; i < nbface; i++)
- {
- fread (&coordx, sizeof (float), 1, is);
- fread (&coordy, sizeof (float), 1, is);
- fread (&coordz, sizeof (float), 1, is);
- for (unsigned int j = 0; j < 3; j++)
- {
- //read the vertex
- fread (&coordx, sizeof (float), 1, is);
- fread (&coordy, sizeof (float), 1, is);
- fread (&coordz, sizeof (float), 1, is);
- if (i == 0 && j == 0)
- {
- m->xmin = coordx;
- m->xmax = coordx;
- m->ymin = coordy;
- m->ymax = coordy;
- m->zmin = coordz;
- m->zmax = coordz;
- }
- else
- {
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- }
- }
- fread (&finligne, 1, 1, is);
- fread (&finligne, 1, 1, is);
- }
- point3d_init (&ptmin, m->xmin, m->ymin, m->zmin);
- point3d_init (&ptmax, m->xmax, m->ymax, m->zmax);
- int nb_part = (int) (pow (nbface/22.0, (1.0 / 3.0))) + 1;
- space_partition_new (&sp, ptmin, ptmax, nb_part, nb_part, nb_part);
- //retour au debut du fichier
- fseek (is, 0, SEEK_SET);
- fread (&str, sizeof (char), 80, is);
- fread (&nbface, sizeof (unsigned int), 1, is);
- for (unsigned int i = 0; i < nbface; i++)
- {
- fread (&coordx, sizeof (float), 1, is);
- fread (&coordy, sizeof (float), 1, is);
- fread (&coordz, sizeof (float), 1, is);
- for (unsigned int j = 0; j < 3; j++)
- {
- fread (&coordx, sizeof (float), 1, is);
- fread (&coordy, sizeof (float), 1, is);
- fread (&coordz, sizeof (float), 1, is); //read the vertex
- point3d_init (&ptmin, coordx, coordy, coordz);
- if (i == 0 && j == 0)
- {
- ve1 = 0;
- ptmin.att_int = space_partition_nb_points (&sp);
- space_partition_add_point (&sp, ptmin);
- a2ri_vf_add_vertex (m, ptmin.x, ptmin.y, ptmin.z);
- }
- else
- {
- point3d *ret=space_partition_contains_point (&sp, ptmin);
- if (ret!=NULL)
- {
- switch (j)
- {
- case 0:
- ve1 = ret->att_int;
- break;
- case 1:
- ve2 = ret->att_int;
- break;
- default:
- ve3 = ret->att_int;
- }
- }
- else
- {
- ptmin.att_int = space_partition_nb_points (&sp);
- switch (j)
- {
- case 0:
- ve1 = ptmin.att_int;
- break;
- case 1:
- ve2 = ptmin.att_int;
- break;
- default:
- ve3 = ptmin.att_int;
- }
- space_partition_add_point (&sp, ptmin);
- a2ri_vf_add_vertex (m, ptmin.x, ptmin.y, ptmin.z);
- }
- }
- }
- a2ri_vf_add_face (m, ve1, ve2, ve3);
- fread (&finligne, 1, 1, is);
- fread (&finligne, 1, 1, is);
- }
- fclose (is);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Ouverture d'un fichier stl ascii contenant un modèle
- @param file chemin du fichier ÃÂ ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_load_ascii_stl (
- char *file,
- vf_model * m)
- {
- FILE *is = NULL;
- int ve1,
- ve2,
- ve3,
- i = 0,
- fin = 0;
- double coordx,
- coordy,
- coordz;
- char str[250];
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- fscanf (is, "%s", str);
- while (strcmp (str, "facet"))
- fscanf (is, "%s", str);
- while (!fin)
- {
- fscanf (is, "%s", str); //read normal
- fscanf (is, "%s", str);
- fscanf (is, "%s", str);
- fscanf (is, "%s", str); //read the 3 values of the normal
- fscanf (is, "%s", str);
- fscanf (is, "%s", str); //outer loop
- fscanf (is, "%s", str); //read vertex
- fscanf (is, "%s", str);
- coordx = atof (str); //coordx
- fscanf (is, "%s", str);
- coordy = atof (str); //coordy
- fscanf (is, "%s", str);
- coordz = atof (str); //coordz
- ve1 = a2ri_vf_search_vertex (*m, coordx, coordy, coordz);
- if (ve1 == -1)
- {
- ve1 = m->nbvertex;
- a2ri_vf_add_vertex (m, coordx, coordy, coordz);
- }
- if (i == 0)
- {
- m->xmin = coordx;
- m->xmax = coordx;
- m->ymin = coordy;
- m->ymax = coordy;
- m->zmin = coordz;
- m->zmax = coordz;
- i++;
- }
- else
- {
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- }
- fscanf (is, "%s", str); //read vertex
- fscanf (is, "%s", str);
- coordx = atof (str); //coordx
- fscanf (is, "%s", str);
- coordy = atof (str); //coordy
- fscanf (is, "%s", str);
- coordz = atof (str); //coordz
- ve2 = a2ri_vf_search_vertex (*m, coordx, coordy, coordz);
- if (ve2 == -1)
- {
- ve2 = m->nbvertex;
- a2ri_vf_add_vertex (m, coordx, coordy, coordz);
- }
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- fscanf (is, "%s", str); //read vertex
- fscanf (is, "%s", str);
- coordx = atof (str); //coordx
- fscanf (is, "%s", str);
- coordy = atof (str); //coordy
- fscanf (is, "%s", str);
- coordz = atof (str); //coordz;
- ve3 = a2ri_vf_search_vertex (*m, coordx, coordy, coordz);
- if (ve3 == -1)
- {
- ve3 = m->nbvertex;
- a2ri_vf_add_vertex (m, coordx, coordy, coordz);
- }
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- a2ri_vf_add_face (m, ve1, ve2, ve3);
- fscanf (is, "%s", str); //read endloop
- fscanf (is, "%s", str); //read endfacet
- fscanf (is, "%s", str); //read facet ou endsolid
- fin = !strcmp (str, "endsolid");
- }
- fclose (is);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Enregistrement d'un modèle dans un format défini par l'extension du nom de fichier
- @param file nom du fichier
- @param m le modèle
- @return aucun
- */
- int
- a2ri_vf_save_file (
- char *file,
- vf_model m)
- {
- char ext[4];
- memcpy (ext, file + (strlen (file) - 3), 3);
- ext[3] = 0;
- //enregistrement du modèle en fonction de l'extension
- if (!strcmp (ext, "wrl"))
- {
- return a2ri_vf_save_wrl (file, m);
- }
- if (!strcmp (ext, "vef"))
- {
- return a2ri_vf_save_vef (file, m);
- }
- if (!strcmp (ext, "pgn"))
- {
- return a2ri_vf_save_pgn (file, m);
- }
- if (!strcmp (ext, "ply"))
- {
- return a2ri_vf_save_ply (file, m);
- }
- if (!strcmp (ext, "off"))
- {
- return a2ri_vf_save_off (file, m);
- }
- if (!strcmp (ext, "obj"))
- {
- return a2ri_vf_save_obj (file, m);
- }
- if (!strcmp (ext, "gts"))
- {
- return a2ri_vf_save_gts (file, m);
- }
- if (!strcmp (ext, "stl"))
- {
- return a2ri_vf_save_binary_stl (file, m);
- }
- return 0;
- }
- /**
- Enregistrement d'un modèle dans un fichier wr@param file chemin du fichier
- @param m modèle à sauvegarder
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_save_wrl (
- char *file,
- vf_model m)
- {
- setlocale(LC_NUMERIC,"C");
-
- if (strlen (file) < 4)
- {
- strcat (file, ".wrl");
- }
- char ext[5];
- memcpy (ext, file + (strlen (file) - 4), 4);
- ext[4] = '\0';
- if (strcmp (ext, ".wrl"))
- strcat (file, ".wrl");
- FILE *os = NULL;
- os = fopen (file, "w");
- if (os == NULL)
- return 0;
- //ecriture de l'en-tete
- fprintf (os, "#VRML V2.0 utf8\n");
- fprintf (os, "Shape{\n");
- fprintf (os, "appearance Appearance {\n");
- fprintf (os, "material Material {\n");
- fprintf (os, "}\n");
- fprintf (os, "}\n");
- fprintf (os, "geometry IndexedFaceSet{\n");
- fprintf (os, "coord Coordinate{\n");
- fprintf (os, "point[\n");
- //ecriture des sommets
- for (int i = 0; i < m.nbvertex; i++)
- fprintf (os, "%f %f %f,\n", m.ve[i].x, m.ve[i].y, m.ve[i].z);
- fprintf (os, "]\n");
- fprintf (os, "}\n");
- //ecriture des faces
- fprintf (os, "coordIndex[\n");
- for (int i = 0; i < m.nbface; i++)
- {
- int nb1 = m.fa[i].ve1;
- int nb2 = m.fa[i].ve2;
- int nb3 = m.fa[i].ve3;
- fprintf (os, "%d %d %d -1\n", nb1, nb2, nb3);
- }
- fprintf (os, "]\n");
- fprintf (os, "}\n");
- fprintf (os, "}\n");
- fclose (os);
-
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Enregistrement d'un modèle dans un fichier wrl version 1.0
- @param file chemin du fichier
- @param m modèle à sauvegarder
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_save_wrl_v1 (
- char *file,
- vf_model m)
- {
- setlocale(LC_NUMERIC,"C");
-
- if (strlen (file) < 4)
- {
- strcat (file, ".wrl");
- }
- char ext[5];
- memcpy (ext, file + (strlen (file) - 4), 4);
- ext[4] = '\0';
- if (strcmp (ext, ".wrl"))
- strcat (file, ".wrl");
- FILE *os = NULL;
- os = fopen (file, "w");
- if (os == NULL)
- return 0;
- //ecriture de l'en-tete
- fprintf (os, "#VRML V1.0 ascii\n");
- fprintf (os, "Separator {\n");
- fprintf (os, "Switch {\n");
- fprintf (os, "DEF Pol{\n");
- fprintf (os, "Separator {\n");
- fprintf (os, "Coordinate3 {\n");
- fprintf (os, "point[\n");
- //ecriture des sommets
- for (int i = 0; i < m.nbvertex; i++)
- fprintf (os, "%f %f %f,\n", m.ve[i].x, m.ve[i].y, m.ve[i].z);
- fprintf (os, "]\n");
- fprintf (os, "}\n");
- //ecriture des faces
- fprintf (os, "IndexedFaceSet {\n");
- fprintf (os, "coordIndex [\n");
- for (int i = 0; i < m.nbface; i++)
- {
- int nb1 = m.fa[i].ve1;
- int nb2 = m.fa[i].ve2;
- int nb3 = m.fa[i].ve3;
- fprintf (os, "%d, %d, %d, -1,\n", nb1, nb2, nb3);
- }
- fprintf (os, "]\n");
- fprintf (os, "}\n");
- fprintf (os, "}\n");
- fprintf (os, "}\n");
- fclose (os);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Enregistrement d'un modèle dans un fichier pgn
- @param file chemin du fichier
- @param m modèle à sauvegarder
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_save_pgn (
- char *file,
- vf_model m)
- {
- setlocale(LC_NUMERIC,"C");
- if (strlen (file) < 4)
- {
- strcat (file, ".pgn");
- }
- char ext[5];
- memcpy (ext, file + (strlen (file) - 4), 4);
- ext[4] = '\0';
- if (strcmp (ext, ".pgn"))
- strcat (file, ".pgn");
- FILE *os = NULL;
- os = fopen (file, "w");
- if (os == NULL)
- return 0;
- //ecriture du nombre de sommets et de faces
- fprintf (os, "%d %d\n", m.nbvertex, m.nbface);
- //ecriture des sommets
- for (int i = 0; i < m.nbvertex; i++)
- fprintf (os, "%lf %lf %lf\n", m.ve[i].x, m.ve[i].y, m.ve[i].z);
- fprintf (os, "\n");
- //ecriture des faces
- for (int i = 0; i < m.nbface; i++)
- fprintf (os, "%d %d %d -1\n", m.fa[i].ve1, m.fa[i].ve2, m.fa[i].ve3);
- fclose (os);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Enregistrement d'un modèle dans un fichier vef
- @param file chemin du fichier
- @param m modèle à sauvegarder
- @return 1 si succès, 0 sinon
- **/
- int
- a2ri_vf_save_vef (
- char *file,
- vf_model m)
- {
- hashtable *table = hashtable_new (m.nbvertex + m.nbface);
- hashtable *tabletemp = NULL;
- int p1,
- p2;
- setlocale(LC_NUMERIC,"C");
- if (strlen (file) < 4)
- {
- strcat (file, ".vef");
- }
- char ext[5];
- memcpy (ext, file + (strlen (file) - 4), 4);
- ext[4] = '\0';
- if (strcmp (ext, ".vef"))
- strcat (file, ".vef");
- FILE *os = NULL;
- os = fopen (file, "w");
- if (os == NULL)
- return 0;
- tabletemp = a2ri_vf_construction_edge_table (&m, NULL, 0);
- //ecriture du nombre de sommets, d'aretes et de faces
- fprintf (os, "%d %d %d\n", m.nbvertex, hashtable_size (tabletemp),
- m.nbface);
- hashtable_free (tabletemp);
- //ecriture des sommets
- for (int i = 0; i < m.nbvertex; i++)
- fprintf (os, "%lf %lf %lf\n", m.ve[i].x, m.ve[i].y, m.ve[i].z);
- fprintf (os, "\n");
- //ecriture des aretes
- IFsave_vef_write_edge (&m, table, os);
- fprintf (os, "\n");
- //ecriture des faces
- for (int i = 0; i < m.nbface; i++)
- {
- p1 = m.fa[i].ve1;
- p2 = m.fa[i].ve2;
- //recherche de l'arete dans la table de hachage
- vf_edge *edgetemp = hashtable_look_for (table, p1, p2);
- fprintf (os, "%d ", edgetemp->att_int);
- p1 = m.fa[i].ve2;
- p2 = m.fa[i].ve3;
- edgetemp = hashtable_look_for (table, p1, p2);
- fprintf (os, "%d ", edgetemp->att_int);
- p1 = m.fa[i].ve3;
- p2 = m.fa[i].ve1;
- edgetemp = hashtable_look_for (table, p1, p2);
- fprintf (os, "%d\n", edgetemp->att_int);
- }
- hashtable_free (table);
- fclose (os);
- setlocale(LC_NUMERIC,"");
- return 0;
- }
- /**
- Enregistrement d'un modèle dans un fichier gts
- @param file chemin du fichier
- @param m modèle à sauvegarder
- @return 1 si succès, 0 sinon
- **/
- int
- a2ri_vf_save_gts (
- char *file,
- vf_model m)
- {
- hashtable *table = hashtable_new (m.nbvertex + m.nbface);
- hashtable *tabletemp = NULL;
- int p1,
- p2;
- setlocale(LC_NUMERIC,"C");
- if (strlen (file) < 4)
- {
- strcat (file, ".gts");
- }
- char ext[5];
- memcpy (ext, file + (strlen (file) - 4), 4);
- ext[4] = '\0';
- if (strcmp (ext, ".gts"))
- strcat (file, ".gts");
- FILE *os = NULL;
- os = fopen (file, "w");
- if (os == NULL)
- return 0;
- tabletemp = a2ri_vf_construction_edge_table (&m, NULL, 0);
- //ecriture du nombre de sommets, d'aretes et de faces
- fprintf (os, "%d %d %d\n", m.nbvertex, hashtable_size (tabletemp),
- m.nbface);
- hashtable_free (tabletemp);
- //ecriture des sommets
- for (int i = 0; i < m.nbvertex; i++)
- fprintf (os, "%lf %lf %lf\n", m.ve[i].x, m.ve[i].y, m.ve[i].z);
- fprintf (os, "\n");
- //ecriture des aretes
- IFsave_gts_write_edge (&m, table, os);
- fprintf (os, "\n");
- //ecriture des faces
- for (int i = 0; i < m.nbface; i++)
- {
- p1 = m.fa[i].ve1;
- p2 = m.fa[i].ve2;
- //recherche de l'arete dans la table de hachage
- vf_edge *edgetemp = hashtable_look_for (table, p1, p2);
- fprintf (os, "%d ", edgetemp->att_int + 1);
- p1 = m.fa[i].ve2;
- p2 = m.fa[i].ve3;
- edgetemp = hashtable_look_for (table, p1, 2);
- fprintf (os, "%d ", edgetemp->att_int + 1);
- p1 = m.fa[i].ve3;
- p2 = m.fa[i].ve1;
- edgetemp = hashtable_look_for (table, p1, p2);
- fprintf (os, "%d\n", edgetemp->att_int + 1);
- }
- hashtable_free (table);
- fclose (os);
- setlocale(LC_NUMERIC,"");
- return 0;
- }
- /**
- Enregistrement d'un modèle dans un fichier off
- @param file chemin du fichier
- @param m modèle à sauvegarder
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_save_off (
- char *file,
- vf_model m)
- {
- setlocale(LC_NUMERIC,"C");
- if (strlen (file) < 4)
- {
- strcat (file, ".off");
- }
- char ext[5];
- memcpy (ext, file + (strlen (file) - 4), 4);
- ext[4] = '\0';
- if (strcmp (ext, ".off"))
- strcat (file, ".off");
- FILE *os = NULL;
- os = fopen (file, "w");
- if (os == NULL)
- return 0;
- //ecriture du nombre de sommets et de faces
- fprintf (os, "OFF\n%d %d 0\n", m.nbvertex, m.nbface);
- //ecriture des sommets
- for (int i = 0; i < m.nbvertex; i++)
- fprintf (os, "%lf %lf %lf\n", m.ve[i].x, m.ve[i].y, m.ve[i].z);
- fprintf (os, "\n");
- //ecriture des faces
- for (int i = 0; i < m.nbface; i++)
- fprintf (os, "3 %d %d %d\n", m.fa[i].ve1, m.fa[i].ve2, m.fa[i].ve3);
- fclose (os);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Enregistrement d'un modèle dans un fichier ply
- @param file chemin du fichier
- @param m modèle à sauvegarder
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_save_ply (
- char *file,
- vf_model m)
- {
- int ve1,
- ve2,
- ve3;
- setlocale(LC_NUMERIC,"C");
- if (strlen (file) < 4)
- {
- strcat (file, ".ply");
- }
- char ext[5];
- memcpy (ext, file + (strlen (file) - 4), 4);
- ext[4] = '\0';
- if (strcmp (ext, ".ply"))
- strcat (file, ".ply");
- FILE *os = NULL;
- os = fopen (file, "w");
- if (os == NULL)
- return 0;
- //ecriture des en-tetes
- fprintf (os, "ply\nformat ascii 1.0\ncomment zipper output\n");
- fprintf (os,
- "element vertex %d\nproperty float x\nproperty float y\nproperty float z\nelement face %d\n",
- m.nbvertex, m.nbface);
- fprintf (os, "property list uchar int vertex_indices\nend_header\n");
- //ecriture des sommets
- for (int i = 0; i < m.nbvertex; i++)
- fprintf (os, "%f %f %f\n", m.ve[i].x, m.ve[i].y, m.ve[i].z);
- //ecriture des faces
- for (int i = 0; i < m.nbface; i++)
- {
- ve1 = m.fa[i].ve1;
- ve2 = m.fa[i].ve2;
- ve3 = m.fa[i].ve3;
- fprintf (os, "3 %d %d %d \n", ve1, ve2, ve3);
- }
- fclose (os);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Enregistrement d'un modèle dans un fichier obj
- @param file chemin du fichier
- @param m modèle à sauvegarder
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vf_save_obj (
- char *file,
- vf_model m)
- {
- setlocale(LC_NUMERIC,"C");
- if (strlen (file) < 4)
- {
- strcat (file, ".obj");
- }
- char ext[5];
- memcpy (ext, file + (strlen (file) - 4), 4);
- ext[4] = '\0';
- if (strcmp (ext, ".obj"))
- strcat (file, ".obj");
- FILE *os = NULL;
- os = fopen (file, "w");
- if (os == NULL)
- return 0;
- //ecriture des sommets
- for (int i = 0; i < m.nbvertex; i++)
- fprintf (os, "v %lf %lf %lf\n", m.ve[i].x, m.ve[i].y, m.ve[i].z);
- fprintf (os, "\n");
- //ecriture des faces
- for (int i = 0; i < m.nbface; i++)
- fprintf (os, "f %d %d %d\n", m.fa[i].ve1 + 1, m.fa[i].ve2 + 1,
- m.fa[i].ve3 + 1);
- fclose (os);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Enregistrement d'un modèle dans un fichier stl binaire
- @param file chemin du fichier
- @param m modèle à sauvegarder
- @return 1 si succés, 0 sinon
- */
- int
- a2ri_vf_save_binary_stl (
- char *file,
- vf_model m)
- {
- setlocale(LC_NUMERIC,"C");
- if (strlen (file) < 4)
- {
- strcat (file, ".stl");
- }
- char ext[5];
- memcpy (ext, file + (strlen (file) - 4), 4);
- ext[4] = '\0';
- if (strcmp (ext, ".stl"))
- strcat (file, ".stl");
- FILE *os = NULL;
- os = fopen (file, "wb");
- if (os == NULL)
- return 0;
- fwrite
- ("generated by A2RI ",
- 1, 80, os);
- fwrite (&(m.nbface), sizeof (unsigned int), 1, os);
- for (int i = 0; i < m.nbface; i++)
- {
- vector3d AB,
- AC,
- normale;
- vector3d_init (&AB,
- m.ve[m.fa[i].ve2].x - m.ve[m.fa[i].ve1].x,
- m.ve[m.fa[i].ve2].y - m.ve[m.fa[i].ve1].y,
- m.ve[m.fa[i].ve2].z - m.ve[m.fa[i].ve1].z);
- vector3d_init (&AC,
- m.ve[m.fa[i].ve3].x - m.ve[m.fa[i].ve1].x,
- m.ve[m.fa[i].ve3].y - m.ve[m.fa[i].ve1].y,
- m.ve[m.fa[i].ve3].z - m.ve[m.fa[i].ve1].z);
- normale = vector3d_vectorialproduct (AB, AC);
- float x,
- y,
- z;
- x = normale.dx;
- y = normale.dy;
- z = normale.dz;
- fwrite (&x, sizeof (float), 1, os);
- fwrite (&y, sizeof (float), 1, os);
- fwrite (&z, sizeof (float), 1, os);
- x = m.ve[m.fa[i].ve1].x;
- y = m.ve[m.fa[i].ve1].y;
- z = m.ve[m.fa[i].ve1].z;
- fwrite (&x, sizeof (float), 1, os);
- fwrite (&y, sizeof (float), 1, os);
- fwrite (&z, sizeof (float), 1, os);
- x = m.ve[m.fa[i].ve2].x;
- y = m.ve[m.fa[i].ve2].y;
- z = m.ve[m.fa[i].ve2].z;
- fwrite (&x, sizeof (float), 1, os);
- fwrite (&y, sizeof (float), 1, os);
- fwrite (&z, sizeof (float), 1, os);
- x = m.ve[m.fa[i].ve3].x;
- y = m.ve[m.fa[i].ve3].y;
- z = m.ve[m.fa[i].ve3].z;
- fwrite (&x, sizeof (float), 1, os);
- fwrite (&y, sizeof (float), 1, os);
- fwrite (&z, sizeof (float), 1, os);
- fwrite (" ", 1, 2, os);
- }
- fclose (os);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Enregistrement d'un modèle dans un fichier stl ascii
- @param file chemin du fichier
- @param m modèle à sauvegarder
- @return 1 si succés, 0 sinon
- */
- int
- a2ri_vf_save_ascii_stl (
- char *file,
- vf_model m)
- {
- setlocale(LC_NUMERIC,"C");
- if (strlen (file) < 4)
- {
- strcat (file, ".stl");
- }
- char ext[5];
- memcpy (ext, file + (strlen (file) - 4), 4);
- ext[4] = '\0';
- if (strcmp (ext, ".stl"))
- strcat (file, ".stl");
- FILE *os = NULL;
- os = fopen (file, "w");
- if (os == NULL)
- return 0;
- fprintf(os,"solid generatedByA2RI\n");
- for (int i = 0; i < m.nbface; i++)
- {
- vector3d AB,
- AC,
- normale;
- vector3d_init (&AB,
- m.ve[m.fa[i].ve2].x - m.ve[m.fa[i].ve1].x,
- m.ve[m.fa[i].ve2].y - m.ve[m.fa[i].ve1].y,
- m.ve[m.fa[i].ve2].z - m.ve[m.fa[i].ve1].z);
- vector3d_init (&AC,
- m.ve[m.fa[i].ve3].x - m.ve[m.fa[i].ve1].x,
- m.ve[m.fa[i].ve3].y - m.ve[m.fa[i].ve1].y,
- m.ve[m.fa[i].ve3].z - m.ve[m.fa[i].ve1].z);
- normale = vector3d_vectorialproduct (AB, AC);
- float x,
- y,
- z;
- x = normale.dx;
- y = normale.dy;
- z = normale.dz;
- fprintf (os,"facet normal %f %f %f\n\touter loop\n\t\tvertex ",x,y,z);
- x = m.ve[m.fa[i].ve1].x;
- y = m.ve[m.fa[i].ve1].y;
- z = m.ve[m.fa[i].ve1].z;
- fprintf (os,"%f %f %f\n\t\tvertex ",x,y,z);
- x = m.ve[m.fa[i].ve2].x;
- y = m.ve[m.fa[i].ve2].y;
- z = m.ve[m.fa[i].ve2].z;
- fprintf (os,"%f %f %f\n\t\tvertex ",x,y,z);
- x = m.ve[m.fa[i].ve3].x;
- y = m.ve[m.fa[i].ve3].y;
- z = m.ve[m.fa[i].ve3].z;
- fprintf (os,"%f %f %f\n\tendloop\nendfacet\n",x,y,z);
- }
- fprintf(os,"endsolid generatedByA2RI\n");
- fclose (os);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Ouverture d'un fichier contenant un modèle
- @param file chemin du fichier à ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vef_open_file (
- char *file,
- vef_model * m)
- {
- char ext[5];
- memcpy (ext, file + (strlen (file) - 4), 4);
- ext[4] = 0;
- //appel de la fonction correspondant à l'extension du fichier
- if (!strcmp (ext, ".wrl"))
- return (a2ri_vef_load_wrl (file, m));
- if (!strcmp (ext, ".vef"))
- return (a2ri_vef_load_vef (file, m));
- if (!strcmp (ext, ".pgn"))
- return (a2ri_vef_load_pgn (file, m));
- if (!strcmp (ext, ".ply"))
- return (a2ri_vef_load_ply (file, m));
- if (!strcmp (ext, ".off"))
- return (a2ri_vef_load_off (file, m));
- if (!strcmp (ext, ".obj"))
- return (a2ri_vef_load_obj (file, m));
- if (!strcmp (ext, ".gts"))
- return (a2ri_vef_load_gts (file, m));
- if (!strcmp (ext, ".stl"))
- return (a2ri_vef_load_stl (file, m));
- return 0;
- }
- /**
- Ouverture d'un fichier pgn contenant un modèle
- @param file chemin du fichier à ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vef_load_pgn (
- char *file,
- vef_model * m)
- {
- FILE *is = NULL;
- double coordx,
- coordy,
- coordz;
- int classif,
- ar1,
- ar2,
- ar3,
- nbvertex,
- nbface;
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- a2ri_vef_init (m);
- //lecture du nombre de sommets
- fscanf (is, "%d", &nbvertex);
- //lecture du nombre de faces
- fscanf (is, "%d", &nbface);
- //lecture des sommets
- for (int i = 0; i < nbvertex; i++)
- {
- fscanf (is, "%lf", &coordx);
- fscanf (is, "%lf", &coordy);
- fscanf (is, "%lf", &coordz);
- a2ri_vef_add_vertex (m, coordx, coordy, coordz);
- //mis à jour de la bounding box
- if (i == 0)
- {
- m->xmin = coordx;
- m->xmax = coordx;
- m->ymin = coordy;
- m->ymax = coordy;
- m->zmin = coordz;
- m->zmax = coordz;
- }
- else
- {
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- }
- }
- //lecture des faces
- for (int i = 0; i < nbface; i++)
- {
- //lecture des trois sommets composant la face
- fscanf (is, "%d", &ar1);
- fscanf (is, "%d", &ar2);
- fscanf (is, "%d", &ar3);
- //ajout des aretes avec vérification au préalable de leur existence
- a2ri_vef_add_edge (m, ar1, ar2, 1);
- a2ri_vef_add_edge (m, ar2, ar3, 1);
- a2ri_vef_add_edge (m, ar3, ar1, 1);
- //recherche des aretes
- int ed1 = a2ri_vef_search_edge (*m, ar1, ar2);
- int ed2 = a2ri_vef_search_edge (*m, ar2, ar3);
- int ed3 = a2ri_vef_search_edge (*m, ar3, ar1);
- //ajout de la faces avec les trois aretes
- a2ri_vef_add_face (m, ed1, ed2, ed3);
- // reading the -1
- fscanf (is, "%d", &ar1);
- }
- //lecture des classifications si disponible
- for (int i = 0; i < nbface; i++)
- {
- if (!feof (is))
- {
- fscanf (is, "%d", &classif);
- m->fa[i].type = classif;
- }
- }
- fclose (is);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Ouverture d'un fichier stl binaire contenant un modèle
- @param file chemin du fichier ÃÂ ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vef_load_binary_stl (
- char *file,
- vef_model * m)
- {
- FILE *is = NULL;
- char str[80];
- unsigned int nbface;
- unsigned char finligne;
- float coordx,
- coordy,
- coordz;
- int ve1,
- ve2,
- ve3,
-
- ed1,
- ed2,
- ed3;
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "rb");
- if (is == NULL)
- return 0;
- fread (&str, sizeof (char), 80, is);
- fread (&nbface, sizeof (unsigned int), 1, is);
- for (unsigned int i = 0; i < nbface; i++)
- {
- fread (&coordx, sizeof (float), 1, is);
- fread (&coordy, sizeof (float), 1, is);
- fread (&coordz, sizeof (float), 1, is); //read the normal
- fread (&coordx, sizeof (float), 1, is);
- fread (&coordy, sizeof (float), 1, is);
- fread (&coordz, sizeof (float), 1, is); //read the first vertex
- ve1 = a2ri_vef_search_vertex (*m, coordx, coordy, coordz);
- if (ve1 == -1)
- {
- ve1 = m->nbvertex;
- a2ri_vef_add_vertex (m, coordx, coordy, coordz);
- }
- if (i == 0)
- {
- m->xmin = coordx;
- m->xmax = coordx;
- m->ymin = coordy;
- m->ymax = coordy;
- m->zmin = coordz;
- m->zmax = coordz;
- }
- else
- {
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- }
- fread (&coordx, sizeof (float), 1, is);
- fread (&coordy, sizeof (float), 1, is);
- fread (&coordz, sizeof (float), 1, is); //read the second vertex
- ve2 = a2ri_vef_search_vertex (*m, coordx, coordy, coordz);
- if (ve2 == -1)
- {
- ve2 = m->nbvertex;
- a2ri_vef_add_vertex (m, coordx, coordy, coordz);
- }
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- fread (&coordx, sizeof (float), 1, is);
- fread (&coordy, sizeof (float), 1, is);
- fread (&coordz, sizeof (float), 1, is); //read the third vertex
- ve3 = a2ri_vef_search_vertex (*m, coordx, coordy, coordz);
- if (ve3 == -1)
- {
- ve3 = m->nbvertex;
- a2ri_vef_add_vertex (m, coordx, coordy, coordz);
- }
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- ed1 = a2ri_vef_search_edge (*m, ve1, ve2);
- ed2 = a2ri_vef_search_edge (*m, ve1, ve3);
- ed3 = a2ri_vef_search_edge (*m, ve2, ve3);
- if (ed1 == -1)
- {
- ed1 = m->nbedge;
- a2ri_vef_add_edge (m, ve1, ve2, 0);
- }
- if (ed2 == -1)
- {
- ed2 = m->nbedge;
- a2ri_vef_add_edge (m, ve1, ve3, 0);
- }
- if (ed3 == -1)
- {
- ed3 = m->nbedge;
- a2ri_vef_add_edge (m, ve2, ve3, 0);
- }
- a2ri_vef_add_face (m, ed2, ed1, ed3);
- fread (&finligne, 1, 1, is);
- fread (&finligne, 1, 1, is);
- }
- fclose (is);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Ouverture d'un fichier stl ascii contenant un modèle
- @param file chemin du fichier ÃÂ ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vef_load_ascii_stl (
- char *file,
- vef_model * m)
- {
- FILE *is = NULL;
- int ve1,
- ve2,
- ve3,
- i = 0,
- fin = 0,
- ed1,
- ed2,
- ed3;
- double coordx,
- coordy,
- coordz;
- char str[250];
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- fscanf (is, "%s", str);
- while (strcmp (str, "facet"))
- fscanf (is, "%s", str);
- while (!fin)
- {
- fscanf (is, "%s", str); //read normal
- fscanf (is, "%s", str);
- fscanf (is, "%s", str);
- fscanf (is, "%s", str); //read the 3 values of the normal
- fscanf (is, "%s", str);
- fscanf (is, "%s", str); //outer loop
- fscanf (is, "%s", str); //read vertex
- fscanf (is, "%s", str);
- coordx = atof (str); //coordx
- fscanf (is, "%s", str);
- coordy = atof (str); //coordy
- fscanf (is, "%s", str);
- coordz = atof (str); //coordz
- ve1 = a2ri_vef_search_vertex (*m, coordx, coordy, coordz);
- if (ve1 == -1)
- {
- ve1 = m->nbvertex;
- a2ri_vef_add_vertex (m, coordx, coordy, coordz);
- }
- if (i == 0)
- {
- m->xmin = coordx;
- m->xmax = coordx;
- m->ymin = coordy;
- m->ymax = coordy;
- m->zmin = coordz;
- m->zmax = coordz;
- i++;
- }
- else
- {
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- }
- fscanf (is, "%s", str); //read vertex
- fscanf (is, "%s", str);
- coordx = atof (str); //coordx
- fscanf (is, "%s", str);
- coordy = atof (str); //coordy
- fscanf (is, "%s", str);
- coordz = atof (str); //coordz
- ve2 = a2ri_vef_search_vertex (*m, coordx, coordy, coordz);
- if (ve2 == -1)
- {
- ve2 = m->nbvertex;
- a2ri_vef_add_vertex (m, coordx, coordy, coordz);
- }
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- fscanf (is, "%s", str); //read vertex
- fscanf (is, "%s", str);
- coordx = atof (str); //coordx
- fscanf (is, "%s", str);
- coordy = atof (str); //coordy
- fscanf (is, "%s", str);
- coordz = atof (str); //coordz;
- ve3 = a2ri_vef_search_vertex (*m, coordx, coordy, coordz);
- if (ve3 == -1)
- {
- ve3 = m->nbvertex;
- a2ri_vef_add_vertex (m, coordx, coordy, coordz);
- }
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- ed1 = a2ri_vef_search_edge (*m, ve1, ve2);
- ed2 = a2ri_vef_search_edge (*m, ve1, ve3);
- ed3 = a2ri_vef_search_edge (*m, ve2, ve3);
- if (ed1 == -1)
- {
- ed1 = m->nbedge;
- a2ri_vef_add_edge (m, ve1, ve2, 0);
- }
- if (ed2 == -1)
- {
- ed2 = m->nbedge;
- a2ri_vef_add_edge (m, ve1, ve3, 0);
- }
- if (ed3 == -1)
- {
- ed3 = m->nbedge;
- a2ri_vef_add_edge (m, ve2, ve3, 0);
- }
- a2ri_vef_add_face (m, ed2, ed1, ed3);
- fscanf (is, "%s", str); //read endloop
- fscanf (is, "%s", str); //read endfacet
- fscanf (is, "%s", str); //read facet ou endsolid
- fin = !strcmp (str, "endsolid");
- }
- fclose (is);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Ouverture d'un fichier stl contenant un modèle
- @param file chemin du fichier ÃÂ ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vef_load_stl (
- char *file,
- vef_model * m)
- {
- FILE *is = NULL;
- char str[250];
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- fscanf (is, "%s", str);
- fclose (is);
- if (!strcmp (str, "solid"))
- return a2ri_vef_load_ascii_stl (file, m);
- else
- return a2ri_vef_load_binary_stl (file, m);
- }
- /**
- Ouverture d'un fichier obj contenant un modèle
- @param file chemin du fichier à ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vef_load_obj (
- char *file,
- vef_model * m)
- {
- FILE *is = NULL;
- int ve1,
- ve2,
- ve3;
- double coordx,
- coordy,
- coordz;
- char str[250];
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- fscanf (is, "%s", str);
- while (strcmp (str, "v") != 0)
- fscanf (is, "%s", str);
- //lecture des sommets
- while (strcmp (str, "g") != 0)
- {
- fscanf (is, "%lf %lf %lf", &coordx, &coordy, &coordz);
- //ajout d'un sommet
- a2ri_vef_add_vertex (m, coordx, coordy, coordz);
- fscanf (is, "%s", str);
- //mis à jour de la bounding box
- if (m->nbvertex == 1)
- {
- m->xmin = coordx;
- m->xmax = coordx;
- m->ymin = coordy;
- m->ymax = coordy;
- m->zmin = coordz;
- m->zmax = coordz;
- }
- else
- {
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- }
- }
- fscanf (is, "%s", str);
- fscanf (is, "%s", str);
- //lecture des faces
- while (!feof (is))
- {
- //lecture des trois sommets
- fscanf (is, "%d %d %d", &ve1, &ve2, &ve3);
- ve1--;
- ve2--;
- ve3--;
- //ajout des aretes avec vérification au préalable de leur existence
- a2ri_vef_add_edge (m, ve1, ve2, 1);
- a2ri_vef_add_edge (m, ve2, ve3, 1);
- a2ri_vef_add_edge (m, ve3, ve1, 1);
- //recherche des aretes
- int ed1 = a2ri_vef_search_edge (*m, ve1, ve2);
- int ed2 = a2ri_vef_search_edge (*m, ve2, ve3);
- int ed3 = a2ri_vef_search_edge (*m, ve3, ve1);
- //ajout de la faces avec les trois aretes
- a2ri_vef_add_face (m, ed1, ed2, ed3);
- fscanf (is, "%s", str);
- //si la face compte plus de trois sommets, on relit et on ajoute des aretes et une seconde face
- if (strcmp (str, "f") != 0)
- {
- ve2 = atoi (str);
- ve2--;
- a2ri_vef_add_edge (m, ve1, ve2, 1);
- a2ri_vef_add_edge (m, ve2, ve3, 1);
- a2ri_vef_add_edge (m, ve3, ve1, 1);
- int ed1 = a2ri_vef_search_edge (*m, ve1, ve2);
- int ed2 = a2ri_vef_search_edge (*m, ve2, ve3);
- int ed3 = a2ri_vef_search_edge (*m, ve3, ve1);
- a2ri_vef_add_face (m, ed1, ed3, ed2);
- fscanf (is, "%s", str);
- }
- }
- fclose (is);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Ouverture d'un fichier wrl contenant un modèle
- @param file chemin du fichier à ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vef_load_wrl (
- char *file,
- vef_model * m)
- {
- FILE *is = NULL;
- char str[20];
- char coordx[15],
- coordy[15],
- coordz[15];
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- fscanf (is, "%s", str);
- while (strcmp (str, "point") && strcmp (str, "point["))
- fscanf (is, "%s", str);
- if (!strcmp (str, "point"))
- fscanf (is, "%s", str);
- // reading the file until coordx = } or ]
- fscanf (is, "%s", coordx);
- int i = 0;
- //lecture des sommets
- while (strcmp (coordx, "}") && strcmp (coordx, "]"))
- {
- fscanf (is, "%s", coordy);
- fscanf (is, "%s", coordz);
- if (coordz[strlen (coordz) - 1] == ',')
- {
- memcpy (coordz, coordz, strlen (coordz) - 1);
- coordz[strlen (coordz)] = '\0';
- }
- a2ri_vef_add_vertex (m, atof (coordx), atof (coordy), atof (coordz));
- //mis à jour de la bounding box
- if (i == 0)
- {
- i++;
- m->xmin = atof (coordx);
- m->xmax = atof (coordx);
- m->ymin = atof (coordy);
- m->ymax = atof (coordy);
- m->zmin = atof (coordz);
- m->zmax = atof (coordz);
- }
- else
- {
- if (atof (coordx) < m->xmin)
- m->xmin = atof (coordx);
- if (atof (coordx) > m->xmax)
- m->xmax = atof (coordx);
- if (atof (coordy) < m->ymin)
- m->ymin = atof (coordy);
- if (atof (coordy) > m->ymax)
- m->ymax = atof (coordy);
- if (atof (coordz) < m->zmin)
- m->zmin = atof (coordz);
- if (atof (coordz) > m->zmax)
- m->zmax = atof (coordz);
- }
- fscanf (is, "%s", coordx);
- }
- fscanf (is, "%s", str);
- while (strcmp (str, "coordIndex") && strcmp (str, "coordIndex["))
- fscanf (is, "%s", str);
- if (!strcmp (str, "coordIndex"))
- fscanf (is, "%s", str);
- //lecture des faces
- fscanf (is, "%s", coordx);
- while (strcmp (coordx, "}") && strcmp (coordx, "]"))
- {
- //lecture des trois sommets composant la face
- fscanf (is, "%s", coordy);
- fscanf (is, "%s", coordz);
- //ajout des aretes avec vérification au préalable de leur existence
- a2ri_vef_add_edge (m, atoi (coordx), atoi (coordy), 1);
- a2ri_vef_add_edge (m, atoi (coordy), atoi (coordz), 1);
- a2ri_vef_add_edge (m, atoi (coordz), atoi (coordx), 1);
- //recherche des aretes
- int ed1 = a2ri_vef_search_edge (*m, atoi (coordx), atoi (coordy));
- int ed2 = a2ri_vef_search_edge (*m, atoi (coordy), atoi (coordz));
- int ed3 = a2ri_vef_search_edge (*m, atoi (coordz), atoi (coordx));
- //ajout de la faces avec les trois aretes
- a2ri_vef_add_face (m, ed1, ed2, ed3);
- // reading the -1
- fscanf (is, "%s", coordx);
- fscanf (is, "%s", coordx);
- }
- fclose (is);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Ouverture d'un fichier vef contenant un modèle
- @param file chemin du fichier à ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vef_load_vef (
- char *file,
- vef_model * m)
- {
- FILE *is = NULL;
- int nbvertex,
- nbedge,
- nbface,
- ve1,
- ve2,
- ed1,
- ed2,
- ed3,
- classif;
- double coordx,
- coordy,
- coordz;
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- a2ri_vef_init (m);
- //lecture du nombre de sommets
- fscanf (is, "%d", &nbvertex);
- //lecture du nombre d'aretes
- fscanf (is, "%d", &nbedge);
- //lecture du nombre de face
- fscanf (is, "%d", &nbface);
- //lecture des sommets
- for (int i = 0; i < nbvertex; i++)
- {
- fscanf (is, "%lf", &coordx);
- fscanf (is, "%lf", &coordy);
- fscanf (is, "%lf", &coordz);
- a2ri_vef_add_vertex (m, coordx, coordy, coordz);
- //mis à jour de la bounding box
- if (i == 0)
- {
- m->xmin = coordx;
- m->xmax = coordx;
- m->ymin = coordy;
- m->ymax = coordy;
- m->zmin = coordz;
- m->zmax = coordz;
- }
- else
- {
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- }
- }
- //lecture des aretes
- for (int i = 0; i < nbedge; i++)
- {
- fscanf (is, "%d", &ve1);
- fscanf (is, "%d", &ve2);
- a2ri_vef_add_edge (m, ve1, ve2, 0);
- }
- //lecture des faces
- for (int i = 0; i < nbface; i++)
- {
- fscanf (is, "%d", &ed1);
- fscanf (is, "%d", &ed2);
- fscanf (is, "%d", &ed3);
- a2ri_vef_add_face (m, ed1, ed2, ed3);
- }
- //lecture des classifications si disponible
- for (int i = 0; i < nbface; i++)
- {
- if (!feof (is))
- {
- fscanf (is, "%d", &classif);
- m->fa[i].type = classif;
- }
- }
- fclose (is);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Ouverture d'un fichier gts contenant un modèle
- @param file chemin du fichier à ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vef_load_gts (
- char *file,
- vef_model * m)
- {
- FILE *is = NULL;
- int nbvertex,
- nbedge,
- nbface,
- ve1,
- ve2,
- ed1,
- ed2,
- ed3,
- classif;
- double coordx,
- coordy,
- coordz;
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- a2ri_vef_init (m);
- //lecture du nombre de sommets
- fscanf (is, "%d", &nbvertex);
- //lecture du nombre d'aretes
- fscanf (is, "%d", &nbedge);
- //lecture du nombre de face
- fscanf (is, "%d", &nbface);
- //lecture des sommets
- for (int i = 0; i < nbvertex; i++)
- {
- fscanf (is, "%lf", &coordx);
- fscanf (is, "%lf", &coordy);
- fscanf (is, "%lf", &coordz);
- a2ri_vef_add_vertex (m, coordx, coordy, coordz);
- //mis à jour de la bounding box
- if (i == 0)
- {
- m->xmin = coordx;
- m->xmax = coordx;
- m->ymin = coordy;
- m->ymax = coordy;
- m->zmin = coordz;
- m->zmax = coordz;
- }
- else
- {
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- }
- }
- //lecture des aretes
- for (int i = 0; i < nbedge; i++)
- {
- fscanf (is, "%d", &ve1);
- fscanf (is, "%d", &ve2);
- ve1--;
- ve2--;
- a2ri_vef_add_edge (m, ve1, ve2, 0);
- }
- //lecture des faces
- for (int i = 0; i < nbface; i++)
- {
- fscanf (is, "%d", &ed1);
- fscanf (is, "%d", &ed2);
- fscanf (is, "%d", &ed3);
- ed1--;
- ed2--;
- ed3--;
- a2ri_vef_add_face (m, ed1, ed2, ed3);
- }
- //lecture des classifications si disponible
- for (int i = 0; i < nbface; i++)
- {
- if (!feof (is))
- {
- fscanf (is, "%d", &classif);
- m->fa[i].type = classif;
- }
- }
- fclose (is);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Ouverture d'un fichier ply contenant un modèle
- @param file chemin du fichier à ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vef_load_ply (
- char *file,
- vef_model * m)
- {
- int prop = 0,
- nbpropvertex = 0,
- nbpropface = 0,
- nbvertex,
- nbface,
- ve1,
- ve2,
- ve3;
- char str[20];
- double coordx,
- coordy,
- coordz;
- FILE *is = NULL;
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- fscanf (is, "%s", str);
- if (strcmp (str, "ply"))
- return 0;
- while (strcmp (str, "end_header"))
- {
- fscanf (is, "%s", str);
- if (prop == 1 && !strcmp (str, "property"))
- nbpropvertex++;
- if (prop == 2 && !strcmp (str, "property"))
- nbpropface++;
- if (!strcmp (str, "element"))
- {
- fscanf (is, "%s", str);
- //lecture du nombre de sommets
- if (!strcmp (str, "vertex"))
- {
- prop = 1;
- fscanf (is, "%d", &nbvertex);
- }
- //lecture du nombre de faces
- else
- {
- prop = 1;
- fscanf (is, "%d", &nbface);
- }
- }
- }
- //lecture des sommets
- for (int i = 0; i < nbvertex; i++)
- {
- fscanf (is, "%lf %lf %lf", &coordx, &coordy, &coordz);
- a2ri_vef_add_vertex (m, coordx, coordy, coordz);
- //mis à jour de la bounding box
- if (i == 0)
- {
- m->xmin = coordx;
- m->xmax = coordx;
- m->ymin = coordy;
- m->ymax = coordy;
- m->zmin = coordz;
- m->zmax = coordz;
- }
- else
- {
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- }
- for (int j = 0; j < nbpropvertex - 4; j++)
- fscanf (is, "%s", str);
- }
- //lecture des faces
- for (int i = 0; i < nbface; i++)
- {
- fscanf (is, "%d", &ve1);
- //lecture des trois sommets composant la face
- fscanf (is, "%d %d %d", &ve1, &ve2, &ve3);
- //ajout des aretes avec vérification au préalable de leur existence
- a2ri_vef_add_edge (m, ve1, ve2, 1);
- a2ri_vef_add_edge (m, ve2, ve3, 1);
- a2ri_vef_add_edge (m, ve3, ve1, 1);
- //recherche des aretes
- int ed1 = a2ri_vef_search_edge (*m, ve1, ve2);
- int ed2 = a2ri_vef_search_edge (*m, ve2, ve3);
- int ed3 = a2ri_vef_search_edge (*m, ve3, ve1);
- //ajout de la faces avec les trois aretes
- a2ri_vef_add_face (m, ed1, ed2, ed3);
- }
- fclose (is);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Ouverture d'un fichier off contenant un modèle
- @param file chemin du fichier à ouvrir
- @param m pointeur sur le modèle
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vef_load_off (
- char *file,
- vef_model * m)
- {
- FILE *is = NULL;
- char str[20];
- double coordx,
- coordy,
- coordz;
- int classif,
- ar1,
- ar2,
- ar3,
- nbvertex,
- nbface;
- setlocale(LC_NUMERIC,"C");
- is = fopen (file, "r");
- if (is == NULL)
- return 0;
- a2ri_vef_init (m);
- fscanf (is, "%s", str);
- if (strcmp (str, "OFF"))
- return 0;
- //lecture du nombre de sommets
- fscanf (is, "%d", &nbvertex);
- //lecture du nombre de faces
- fscanf (is, "%d", &nbface);
- fscanf (is, "%d", &ar1);
- //lecture des sommets
- for (int i = 0; i < nbvertex; i++)
- {
- fscanf (is, "%lf", &coordx);
- fscanf (is, "%lf", &coordy);
- fscanf (is, "%lf", &coordz);
- a2ri_vef_add_vertex (m, coordx, coordy, coordz);
- //mis à jour de la bounding box
- if (i == 0)
- {
- m->xmin = coordx;
- m->xmax = coordx;
- m->ymin = coordy;
- m->ymax = coordy;
- m->zmin = coordz;
- m->zmax = coordz;
- }
- else
- {
- if (coordx < m->xmin)
- m->xmin = coordx;
- if (coordx > m->xmax)
- m->xmax = coordx;
- if (coordy < m->ymin)
- m->ymin = coordy;
- if (coordy > m->ymax)
- m->ymax = coordy;
- if (coordz < m->zmin)
- m->zmin = coordz;
- if (coordz > m->zmax)
- m->zmax = coordz;
- }
- }
- //lecture des faces
- for (int i = 0; i < nbface; i++)
- {
- //lecture des trois sommets composant la face
- fscanf (is, "%d", &ar1);
- fscanf (is, "%d", &ar1);
- fscanf (is, "%d", &ar2);
- fscanf (is, "%d", &ar3);
- //ajout des aretes avec vérification au préalable de leur existence
- a2ri_vef_add_edge (m, ar1, ar2, 1);
- a2ri_vef_add_edge (m, ar2, ar3, 1);
- a2ri_vef_add_edge (m, ar3, ar1, 1);
- //recherche des aretes
- int ed1 = a2ri_vef_search_edge (*m, ar1, ar2);
- int ed2 = a2ri_vef_search_edge (*m, ar2, ar3);
- int ed3 = a2ri_vef_search_edge (*m, ar3, ar1);
- //ajout de la faces avec les trois aretes
- a2ri_vef_add_face (m, ed1, ed2, ed3);
- }
- //lecture des classifications si disponible
- for (int i = 0; i < nbface; i++)
- {
- if (!feof (is))
- {
- fscanf (is, "%d", &classif);
- m->fa[i].type = classif;
- }
- }
- fclose (is);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Enregistrement d'un modèle dans un format défini par l'extension du nom de fichier
- @param file nom du fichier
- @param m le modèle
- @return aucun
- */
- int
- a2ri_vef_save_file (
- char *file,
- vef_model m)
- {
- char ext[4];
-
- memcpy (ext, file + (strlen (file) - 3), 3);
- ext[3] = 0;
- //appel de la fonction d'enregistrement en fonction du type de l'extension du fichier
- if (!strcmp (ext, "wrl"))
- return a2ri_vef_save_wrl (file, m);
- if (!strcmp (ext, "vef"))
- return a2ri_vef_save_vef (file, m);
- if (!strcmp (ext, "pgn"))
- return a2ri_vef_save_pgn (file, m);
- if (!strcmp (ext, "ply"))
- return a2ri_vef_save_ply (file, m);
- if (!strcmp (ext, "off"))
- return a2ri_vef_save_off (file, m);
- if (!strcmp (ext, "obj"))
- return a2ri_vef_save_obj (file, m);
- if (!strcmp (ext, "gts"))
- return a2ri_vef_save_gts (file, m);
- return 0;
- }
- /**
- Enregistrement d'un modèle dans un fichier wrl
- @param file chemin du fichier
- @param m modèle à sauvegarder
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vef_save_wrl (
- char *file,
- vef_model m)
- {
- setlocale(LC_NUMERIC,"C");
- if (strlen (file) < 4)
- {
- strcat (file, ".wrl");
- }
- char ext[5];
- memcpy (ext, file + (strlen (file) - 4), 4);
- ext[4] = '\0';
- if (strcmp (ext, ".wrl"))
- strcat (file, ".wrl");
- FILE *os = NULL;
- os = fopen (file, "w");
- if (os == NULL)
- return 0;
- //ecriture de l'en-tete
- fprintf (os, "#VRML V2.0 utf8\n");
- fprintf (os, "Shape{\n");
- fprintf (os, "appearance Appearance {\n");
- fprintf (os, "material Material {\n");
- fprintf (os, "}\n");
- fprintf (os, "}\n");
- fprintf (os, "geometry IndexedFaceSet{\n");
- fprintf (os, "coord Coordinate{\n");
- fprintf (os, "point[\n");
- //ecriture des sommets
- for (int i = 0; i < m.nbvertex; i++)
- fprintf (os, "%f %f %f,\n", m.ve[i].x, m.ve[i].y, m.ve[i].z);
- fprintf (os, "]\n");
- fprintf (os, "}\n");
- //ecriture des faces
- fprintf (os, "coordIndex[\n");
- for (int i = 0; i < m.nbface; i++)
- {
- int nb1,
- nb2,
- nb3;
- vef_face_get_vertices (m.fa[i], m.ed, &nb1, &nb2, &nb3);
- fprintf (os, "%d %d %d -1\n", nb1, nb2, nb3);
- }
- fprintf (os, "]\n");
- fprintf (os, "}\n");
- fprintf (os, "}\n");
- fclose (os);
-
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Enregistrement d'un modèle dans un fichier pgn
- @param file chemin du fichier
- @param m modèle à sauvegarder
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vef_save_pgn (
- char *file,
- vef_model m)
- {
- setlocale(LC_NUMERIC,"C");
- if (strlen (file) < 4)
- {
- strcat (file, ".pgn");
- }
- char ext[5];
- memcpy (ext, file + (strlen (file) - 4), 4);
- ext[4] = '\0';
- if (strcmp (ext, ".pgn"))
- strcat (file, ".pgn");
- FILE *os = NULL;
- os = fopen (file, "w");
- if (os == NULL)
- return 0;
- //ecriture du nombre de sommets et faces
- fprintf (os, "%d %d \n\n", m.nbvertex, m.nbface);
- //ecriture des sommets
- for (int i = 0; i < m.nbvertex; i++)
- fprintf (os, "%f %f %f\n", m.ve[i].x, m.ve[i].y, m.ve[i].z);
- //ecriture des faces
- for (int i = 0; i < m.nbface; i++)
- {
- int ve1,
- ve2,
- ve3;
- vef_face_get_vertices (m.fa[i], m.ed, &ve1, &ve2, &ve3);
- fprintf (os, "%d %d %d -1\n", ve1, ve2, ve3);
- }
- fclose (os);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Enregistrement d'un modèle dans un fichier vef
- @param file chemin du fichier
- @param m modèle à sauvegarder
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vef_save_vef (
- char *file,
- vef_model m)
- {
- setlocale(LC_NUMERIC,"C");
- if (strlen (file) < 4)
- {
- strcat (file, ".vef");
- }
- char ext[5];
- memcpy (ext, file + (strlen (file) - 4), 4);
- ext[4] = '\0';
- if (strcmp (ext, ".vef"))
- strcat (file, ".vef");
- FILE *os = NULL;
- os = fopen (file, "w");
- if (os == NULL)
- return 0;
- //ecriture des nombre de sommets, aretes et faces
- fprintf (os, "%d %d %d\n\n", m.nbvertex, m.nbedge, m.nbface);
- //ecriture des sommets
- for (int i = 0; i < m.nbvertex; i++)
- fprintf (os, "%f %f %f\n", m.ve[i].x, m.ve[i].y, m.ve[i].z);
- //ecriture des aretes
- for (int i = 0; i < m.nbedge; i++)
- fprintf (os, "%d %d\n", m.ed[i].ve1, m.ed[i].ve2);
- //ecriture des faces
- for (int i = 0; i < m.nbface; i++)
- fprintf (os, "%d %d %d\n", m.fa[i].ed1, m.fa[i].ed2, m.fa[i].ed3);
- fclose (os);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Enregistrement d'un modèle dans un fichier gts
- @param file chemin du fichier
- @param m modèle à sauvegarder
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vef_save_gts (
- char *file,
- vef_model m)
- {
- setlocale(LC_NUMERIC,"C");
- if (strlen (file) < 4)
- {
- strcat (file, ".gts");
- }
- char ext[5];
- memcpy (ext, file + (strlen (file) - 4), 4);
- ext[4] = '\0';
- if (strcmp (ext, ".gts"))
- strcat (file, ".gts");
- FILE *os = NULL;
- os = fopen (file, "w");
- if (os == NULL)
- return 0;
- //ecriture des nombre de sommets, aretes et faces
- fprintf (os, "%d %d %d\n\n", m.nbvertex, m.nbedge, m.nbface);
- //ecriture des sommets
- for (int i = 0; i < m.nbvertex; i++)
- fprintf (os, "%f %f %f\n", m.ve[i].x, m.ve[i].y, m.ve[i].z);
- //ecriture des aretes
- for (int i = 0; i < m.nbedge; i++)
- fprintf (os, "%d %d\n", m.ed[i].ve1 + 1, m.ed[i].ve2 + 1);
- //ecriture des faces
- for (int i = 0; i < m.nbface; i++)
- fprintf (os, "%d %d %d\n", m.fa[i].ed1 + 1, m.fa[i].ed2 + 1,
- m.fa[i].ed3 + 1);
- fclose (os);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Enregistrement d'un modèle dans un fichier off
- @param file chemin du fichier
- @param m modèle à sauvegarder
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vef_save_off (
- char *file,
- vef_model m)
- {
- setlocale(LC_NUMERIC,"C");
- if (strlen (file) < 4)
- {
- strcat (file, ".off");
- }
- char ext[5];
- memcpy (ext, file + (strlen (file) - 4), 4);
- ext[4] = '\0';
- if (strcmp (ext, ".off"))
- strcat (file, ".off");
- FILE *os = NULL;
- os = fopen (file, "w");
- if (os == NULL)
- return 0;
- //ecriture de l'en-tete et des nombres de sommets et faces
- fprintf (os, "OFF\n%d %d 0\n\n", m.nbvertex, m.nbface);
- //ecriture des sommets
- for (int i = 0; i < m.nbvertex; i++)
- fprintf (os, "%f %f %f\n", m.ve[i].x, m.ve[i].y, m.ve[i].z);
- //ecriture des faces
- for (int i = 0; i < m.nbface; i++)
- {
- int ve1,
- ve2,
- ve3;
- vef_face_get_vertices (m.fa[i], m.ed, &ve1, &ve2, &ve3);
- fprintf (os, "3 %d %d %d \n", ve1, ve2, ve3);
- }
- fclose (os);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Enregistrement d'un modèle dans un fichier ply
- @param file chemin du fichier
- @param m modèle à sauvegarder
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vef_save_ply (
- char *file,
- vef_model m)
- {
- setlocale(LC_NUMERIC,"C");
- if (strlen (file) < 4)
- {
- strcat (file, ".ply");
- }
- char ext[5];
- memcpy (ext, file + (strlen (file) - 4), 4);
- ext[4] = '\0';
- if (strcmp (ext, ".ply"))
- strcat (file, ".ply");
- FILE *os = NULL;
- os = fopen (file, "w");
- if (os == NULL)
- return 0;
- //ecriture de l'en-tete avec le nombre de sommets et de faces
- fprintf (os, "ply\nformat ascii 1.0\ncomment zipper output\n");
- fprintf (os,
- "element vertex %d\nproperty float x\nproperty float y\nproperty float z\nelement face %d\n",
- m.nbvertex, m.nbface);
- fprintf (os, "property list uchar int vertex_indices\nend_header\n");
- //ecriture des sommets
- for (int i = 0; i < m.nbvertex; i++)
- fprintf (os, "%f %f %f\n", m.ve[i].x, m.ve[i].y, m.ve[i].z);
- //ecriture des faces
- for (int i = 0; i < m.nbface; i++)
- {
- int ve1,
- ve2,
- ve3;
- vef_face_get_vertices (m.fa[i], m.ed, &ve1, &ve2, &ve3);
- fprintf (os, "3 %d %d %d \n", ve1, ve2, ve3);
- }
- fclose (os);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
- /**
- Enregistrement d'un modèle dans un fichier obj
- @param file chemin du fichier
- @param m modèle à sauvegarder
- @return 1 si succès, 0 sinon
- */
- int
- a2ri_vef_save_obj (
- char *file,
- vef_model m)
- {
- setlocale(LC_NUMERIC,"C");
- if (strlen (file) < 4)
- {
- strcat (file, ".obj");
- }
- char ext[5];
- memcpy (ext, file + (strlen (file) - 4), 4);
- ext[4] = '\0';
- if (strcmp (ext, ".obj"))
- strcat (file, ".obj");
- FILE *os = NULL;
- os = fopen (file, "w");
- if (os == NULL)
- return 0;
- //ecriture des sommets
- for (int i = 0; i < m.nbvertex; i++)
- fprintf (os, "v %f %f %f\n", m.ve[i].x, m.ve[i].y, m.ve[i].z);
- fprintf (os, "\n");
- //ecriture des faces
- for (int i = 0; i < m.nbface; i++)
- {
- int ve1,
- ve2,
- ve3;
- vef_face_get_vertices (m.fa[i], m.ed, &ve1, &ve2, &ve3);
- fprintf (os, "f %d %d %d\n", ve1 + 1, ve2 + 1, ve3 + 1);
- }
- fclose (os);
- setlocale(LC_NUMERIC,"");
- return 1;
- }
|