acrwrite.c
Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
#include <stdio.h>
00043
#include <stdlib.h>
00044
#include <string.h>
00045
#include <ctype.h>
00046
00047
00048
#include "idsig.h"
00049
#include "idcnt.h"
00050
#include "idima.h"
00051
#include "idvol.h"
00052
#include "idseq.h"
00053
#include "idacr.h"
00054
#include "idacr-restricted.h"
00055
#include "iderr.h"
00056
#include "idio.h"
00057
#include "idprint.h"
00058
00059
00060
00061
00062
00063
00064
00065
00066
00067
00068
00069
00070
00071
00072
00073
00074
00075
00076
00077
00078
00079 #define DEBUG 0
00080
00081 char *
_IdRecCode=
"ACRNEMA_LIBIDO_1.1";
00082
00083
00084
static ACR_ELEMENT *header=0;
00085
00086
00087
00088
00089
00090
00091
00092
00093
00094
00095
00096
00097
00098
00099
00100
00101
00102
00103
00104
00105
00106
00107
00108
00109
00110
00111
00112
00113
00114
00115
00116
00117
00118
00119
00131 int IdAcrWriteFile(
char *nf,
void *obj)
00132
00133
00134
00135
00136 {
00137
return IdAcrWriteFileExtended(nf,obj,NULL);
00138
00139
00140
00141
00142
00143
00144
00145
00146 }
00147
00148
00149
00150
00151
00152
00153
00154
00155
00156
00157
00158
00159
00160
00161
00162
00163
00164
00165
00166
00167
00168
00169
00170
00171
00172
00173
00174
00175
00176
00177
00178
00179
00180
00181
00182
00183
00184
00185
00186
00187
00188
00189
00190
00191
00204 int IdAcrWriteFileExtended(
char *nf,
void *obj,
ACR_ELEMENT *listeAAjouter)
00205
00206
00207
00208
00209
00210
00211
00212
00213
00214 {
00215 FILE *fp,*fphisto;
00216
PSIGNAL sip = NULL;
00217
PCONTOUR cop = NULL;
00218
PPIMAGE imp = NULL;
00219
PPPVOLUME vop = NULL;
00220
PPPPSEQUENCE sep = NULL;
00221
00222
int i,j,k,ty;
00223
int ndim;
00224
short int sh;
00225 uint32_t i32;
00226
ACR_ELEMENT *p;
00227
00228
int retval=1;
00229
00230
if(obj==NULL) {
IdErrPrintf(
" writeacr --> Pointeur NULL !!!\n");
00231
IdErrno=
IDERR_POINTER_IS_NULL;
00232
return (0);
00233 }
00234
00235 fphisto=fopen(
"_t_histo.tmp",
"r");
00236
00237
00238
00239 header=
_IdAcrAddElement(header,
_IdAcrMakeElement(0x0008,0x0000,4,
BD,0));
00240 header=
_IdAcrAddElement(header,
_IdAcrMakeElement(0x0008,0x0001,4,
BD,0));
00241 header=
_IdAcrAddElement(header,
00242
_IdAcrMakeElement(0x0008,0x0010,strlen(
_IdRecCode),
AT,
00243 (
long int)
_IdRecCode));
00244 header=
_IdAcrAddElement(header,
_IdAcrMakeElement(0x0028,0x0000,4,
BD,0));
00245
00246
if(
IdLibidoType(obj)==
SIG) { ndim=1;
00247 imp=(
PPIMAGE)((
void*)(sip=(
PSIGNAL)obj));
00248 }
00249
else if(
IdLibidoType(obj)==
IMA) { ndim=2;
00250 imp= (
PPIMAGE)obj;
00251 }
00252
else if(
IdLibidoType(obj)==
CNT) { ndim=1;
00253 cop=(
PCONTOUR)obj;
00254 imp=(
PPIMAGE)obj;
00255 }
00256
else if(
IdLibidoType(obj)==
VOL) { ndim=3;
00257 imp=(
PPIMAGE)(vop=(
PPPVOLUME)obj);
00258 }
00259
else if(
IdLibidoType(obj)==
SEQ) { ndim=4;
00260 sep=(
PPPPSEQUENCE)obj;
00261 imp=(
PPIMAGE)obj;
00262 }
00263
00264
else {
IdErrPrintf(
"Type d'Objet non traite par IdAcrWriteFile (%04x)!!!\n",
IdLibidoType(obj));
00265
IdErrno=
IDERR_WRONG_TYPE;
00266
return (0);}
00267
00268
00269 header=
_IdAcrAddElement(header,
_IdAcrMakeElement(0x0028,0x0005,2,
BI,ndim));
00270
00271
00272
00273
00274
00275
00276
00277
00278
00279
00280
00281
00282
00283
00284
00285
00286
00287
00288
00289
00290 header=
_IdAcrAddElement(header,
00291
_IdAcrMakeElement(0x0028,0x0100,2,
BI,
IdSizeOfType(imp)*8));
00292
00293
00294 header=
_IdAcrAddElement(header,
00295
_IdAcrMakeElement(0x0028,0x0101,2,
BI,
IdSizeOfType(imp)*8));
00296
00297
00298 header=
_IdAcrAddElement(header,
00299
_IdAcrMakeElement(0x0028,0x0102,2,
BI,
IdSizeOfType(imp)*8-1));
00300 ty=0;
00301
00302
if(
IdLibidoDataType(imp)==
TY_CHAR) ty=1;
00303
else if(
IdLibidoDataType(imp)==
TY_SHORT) ty=1;
00304
else if(
IdLibidoDataType(imp)==
TY_LONG) ty=1;
00305
else if(
IdLibidoDataType(imp)==
TY_FLOAT) ty=1;
00306
else if(
IdLibidoDataType(imp)==
TY_DOUBLE) ty=1;
00307
00308
else if(
IdLibidoDataType(imp)==
TY_COMPLEX_CHAR) ty=1;
00309
else if(
IdLibidoDataType(imp)==
TY_COMPLEX_SHORT) ty=1;
00310
else if(
IdLibidoDataType(imp)==
TY_COMPLEX_LONG) ty=1;
00311
else if(
IdLibidoDataType(imp)==
TY_COMPLEX_FLOAT) ty=1;
00312
else if(
IdLibidoDataType(imp)==
TY_COMPLEX_DOUBLE) ty=1;
00313
00314
00315
00316 header=
_IdAcrAddElement(header,
00317
_IdAcrMakeElement(0x0028,0x0103,2,
BI,ty));
00318
00319 header=
_IdAcrAddElement(header,
00320
_IdAcrMakeElement(0x7FE0,0x0000,4,
BD,0));
00321
00322
00323
00324
if(
IdLibidoType(obj)==
SIG)
00325 {
00326 sip=(
PSIGNAL )obj;
00327 header=
_IdAcrAddElement(header,
00328
_IdAcrMakeElement(0x0028,0x0010,2,
BI,
IdSigDimX(sip)));
00329
00330 header=
_IdAcrAddElement(header,
00331
_IdAcrMakeElement(0x0028,0x0015,2,
BI,
IdSigUsedNbX(sip)));
00332
00333 header=
_IdAcrAddElement(header,
00334
_IdAcrMakeElement(0x0028,0x00199,2,
BI,((
int *)sip)[-1]));
00335
00336 header=
_IdAcrAddElement(header,
_IdAcrMakeElement(0x7FE0,0x0010,
00337
IdSigDimX(sip)*
IdSizeOfType(sip),
00338
NO,(
long int)sip));
00339
if (
_IdCntPrivate(sip)->_message)
00340 header=
_IdAcrAddElement(header,
00341
_IdAcrMakeElement(0x0028,0x0198, strlen(
_IdCntPrivate(sip)->_message),
AT, (
long int)(
_IdCntPrivate(sip)->_message)));
00342 }
00343
00344
00345
00346
else if(
IdLibidoType(obj)==
CNT)
00347 {
00348
00349 cop=(
PCONTOUR )obj;
00350 header=
_IdAcrAddElement(header,
00351
_IdAcrMakeElement(0x0028,0x0010,2,
BI,
IdCntDimX(cop)));
00352
00353
00354 header=
_IdAcrAddElement(header,
00355
_IdAcrMakeElement(0x0028,0x0015,2,
BI,
IdCntUsedNbX(cop)));
00356
00357 header=
_IdAcrAddElement(header,
00358
_IdAcrMakeElement(0x0028,0x00199,2,
BI,((
int *)cop)[-1]));
00359
00360 header=
_IdAcrAddElement(header,
_IdAcrMakeElement(0x7FE0,0x0010,
00361
IdCntDimX(cop)*
IdSizeOfType(cop),
00362
NO,(
long int)cop));
00363
if (
_IdCntPrivate(cop)->_message)
00364 header=
_IdAcrAddElement(header,
00365
_IdAcrMakeElement(0x0028,0x0198, strlen(
_IdCntPrivate(cop)->_message),
AT, (
long int)(
_IdCntPrivate(cop)->_message)));
00366
00367 }
00368
00369
00370
00371
00372
else if(
IdLibidoType(obj)==
IMA)
00373 {
00374 imp=(
PPIMAGE )obj;
00375 header=
_IdAcrAddElement(header,
00376
_IdAcrMakeElement(0x0028,0x0010,2,
BI,
IdImaDimX(imp)));
00377 header=
_IdAcrAddElement(header,
00378
_IdAcrMakeElement(0x0028,0x0011,2,
BI,
IdImaDimY(imp)));
00379
00380 header=
_IdAcrAddElement(header,
00381
_IdAcrMakeElement(0x0028,0x0015,2,
BI,
IdImaUsedNbX(imp)));
00382
00383 header=
_IdAcrAddElement(header,
00384
_IdAcrMakeElement(0x0028,0x0016,2,
BI,
IdImaUsedNbY(imp)));
00385
00386 header=
_IdAcrAddElement(header,
00387
_IdAcrMakeElement(0x0028,0x00199,2,
BI,((
int *)imp)[-1]));
00388
00389 header=
_IdAcrAddElement(header,
_IdAcrMakeElement(0x7FE0,0x0010,
00390
IdImaDimX(imp)*
IdImaDimY(imp)*
IdSizeOfType(imp),
00391
NO,(
long int)imp));
00392
00393
if (
_IdCntPrivate(imp)->_message)
00394 header=
_IdAcrAddElement(header,
00395
_IdAcrMakeElement(0x0028,0x0198,
00396 strlen(
_IdCntPrivate(imp)->_message),
00397
AT,
00398 (
long int)(
_IdCntPrivate(imp)->_message)));
00399
00400 }
00401
00402
00403
00404
else if(
IdLibidoType(obj)==
VOL)
00405 {
00406 vop=(
PPPVOLUME )obj;
00407 header=
_IdAcrAddElement(header,
00408
_IdAcrMakeElement(0x0028,0x0010,2,
BI,
IdVolDimX(vop)));
00409 header=
_IdAcrAddElement(header,
00410
_IdAcrMakeElement(0x0028,0x0011,2,
BI,
IdVolDimY(vop)));
00411 header=
_IdAcrAddElement(header,
00412
_IdAcrMakeElement(0x0028,0x0012,2,
BI,
IdVolDimZ(vop)));
00413
00414 header=
_IdAcrAddElement(header,
00415
_IdAcrMakeElement(0x0028,0x0015,2,
BI,
IdVolUsedNbX(vop)));
00416
00417 header=
_IdAcrAddElement(header,
00418
_IdAcrMakeElement(0x0028,0x0016,2,
BI,
IdVolUsedNbY(vop)));
00419
00420 header=
_IdAcrAddElement(header,
00421
_IdAcrMakeElement(0x0028,0x0017,2,
BI,
IdVolUsedNbZ(vop)));
00422
00423 header=
_IdAcrAddElement(header,
00424
_IdAcrMakeElement(0x0028,0x00199,2,
BI,((
int *)vop)[-1]));
00425
00426 header=
_IdAcrAddElement(header,
_IdAcrMakeElement(0x7FE0,0x0010,
00427
IdVolDimZ(vop)*
IdVolDimX(vop)*
IdVolDimY(vop)*
00428
IdSizeOfType(vop),
00429
NO,(
long int)vop));
00430
00431
if (
_IdCntPrivate(vop)->_message)
00432 header=
_IdAcrAddElement(header,
00433
_IdAcrMakeElement(0x0028,0x0198,
00434 strlen(
_IdCntPrivate(vop)->_message),
00435
AT,
00436 (
long int)(
_IdCntPrivate(vop)->_message)));
00437
00438 }
00439
00440
else if(
IdLibidoType(obj)==
SEQ)
00441 {
00442 sep=(
PPPPSEQUENCE )obj;
00443 header=
_IdAcrAddElement(header,
00444
_IdAcrMakeElement(0x0028,0x0010,2,
BI,
IdSeqDimX(sep)));
00445 header=
_IdAcrAddElement(header,
00446
_IdAcrMakeElement(0x0028,0x0011,2,
BI,
IdSeqDimY(sep)));
00447 header=
_IdAcrAddElement(header,
00448
_IdAcrMakeElement(0x0028,0x0012,2,
BI,
IdSeqDimZ(sep)));
00449 header=
_IdAcrAddElement(header,
00450
_IdAcrMakeElement(0x0028,0x0013,2,
BI,
IdSeqDimT(sep)));
00451
00452 header=
_IdAcrAddElement(header,
00453
_IdAcrMakeElement(0x0028,0x0015,2,
BI,
IdSeqUsedNbX(sep)));
00454
00455 header=
_IdAcrAddElement(header,
00456
_IdAcrMakeElement(0x0028,0x0016,2,
BI,
IdSeqUsedNbY(sep)));
00457
00458 header=
_IdAcrAddElement(header,
00459
_IdAcrMakeElement(0x0028,0x0017,2,
BI,
IdSeqUsedNbZ(sep)));
00460
00461 header=
_IdAcrAddElement(header,
00462
_IdAcrMakeElement(0x0028,0x0018,2,
BI,
IdSeqUsedNbT(sep)));
00463
00464 header=
_IdAcrAddElement(header,
00465
_IdAcrMakeElement(0x0028,0x00199,2,
BI,((
int *)sep)[-1]));
00466
00467 header=
_IdAcrAddElement(header,
_IdAcrMakeElement(0x7FE0,0x0010,
00468
IdSeqDimT(sep)*
IdSeqDimZ(sep)*
IdSeqDimX(sep)*
IdSeqDimY(sep)*
00469
IdSizeOfType(sep),
00470
NO,(
long int)sep));
00471
00472
if (
_IdCntPrivate(sep)->_message)
00473 header=
_IdAcrAddElement(header,
00474
_IdAcrMakeElement(0x0028,0x0198, strlen(
_IdCntPrivate(sep)->_message),
AT, (
long int)(
_IdCntPrivate(sep)->_message)));
00475
00476 }
00477
00478
00479
00480
if (listeAAjouter != NULL)
00481 {
00482
if (
DEBUG)
00483 {
00484
IdPrintf (
"header %p listeAAjouter %p\n",header,listeAAjouter);
00485
_IdAcrAffListe(header);
00486
IdPrintf (
"\n\n\n");
00487
_IdAcrAffListe(listeAAjouter);
00488
IdPrintf (
"\n\n\n");
00489 }
00490
00491 header=
_IdAcrMergeListes(header,listeAAjouter);
00492
if(!header)
00493
IdPrintf(
"Gros soucis lors de l'incorporation des ACR_ELEMENT utilisateurs... \n");
00494
00495 }
00496
00497
00498
00499
00500
if (
DEBUG) {
IdPrintf (
"\n\nAvant _IdAcrAdjustLength\n");
00501
_IdAcrAffListe(header);
00502
IdPrintf (
"\n\n\n");
00503 }
00504
00505
_IdAcrAdjustLength(header);
00506
00507
00508
if (
DEBUG) {
IdPrintf (
"\n\nApres _IdAcrAdjustLength\n");
00509
_IdAcrAffListe(header);
00510
IdPrintf (
"\n\n\n");
00511 }
00512
00513
if((fp=fopen(nf,
ID_WFILE_BIN))!=0)
00514 {
00515 p=header;
00516
while(p)
00517 {
00518
switch(p->
type)
00519 {
00520
case BI:
00521 fwrite(p,8,1,fp);
00522 sh=(
short int)p->
value;
00523 fwrite(&sh,2,1,fp);
00524
00525
if(
DEBUG)
IdPrintf(
" BI element (en %p) %4x,%4x long.=%d valeur ",
00526 p, p->
group,p->
number,p->
length);
00527
if (
DEBUG)
IdPrintf(
"%d\n",p->
value);
00528
00529
break;
00530
case BD:
00531 fwrite(p,8,1,fp);
00532 i32=p->
value;
00533 fwrite(&i32,4,1,fp);
00534
00535
if(
DEBUG)
IdPrintf(
" BD element (en %p) %4x,%4x long.=%d valeur ",
00536 p, p->
group,p->
number,p->
length);
00537
if (
DEBUG)
IdPrintf(
"%d\n",p->
value);
00538
break;
00539
case AT:
00540 fwrite(p,8,1,fp);
00541 fwrite((
char *)(p->
value),(size_t)p->
length,1,fp);
00542
00543
if(
DEBUG)
IdPrintf(
" AT element (en %p) %4x,%4x long.=%d valeur ",
00544 p, p->
group,p->
number,p->
length);
00545
if (
DEBUG)
IdPrintf(
"%s\n",(
char *)p->
value);
00546
break;
00547
case AN:
00548 fwrite(p,8,1,fp);
00549 fwrite((
char *)(p->
value),(size_t)p->
length,1,fp);
00550
00551
if(
DEBUG)
IdPrintf(
" AN element (en %p) %4x,%4x long.=%d valeur ",
00552 p, p->
group,p->
number,p->
length);
00553
if (
DEBUG)
IdPrintf(
"%s\n",(
char *)p->
value);
00554
00555
break;
00556
case NO:
00557 fwrite(p,8,1,fp);
00558
00559
if(
DEBUG)
IdPrintf(
" NO element (en %p) %4x,%4x long.=%d valeur ",
00560 p, p->
group,p->
number,p->
length);
00561
if (
DEBUG)
IdPrintf(
"%x\n",p->
value);
00562
00563
00564
switch (
IdLibidoType(obj))
00565 {
00566
case SIG:
00567 fwrite(sip,
00568
IdSigDimX(sip)*
IdSizeOfType(sip),1,fp);
00569
break;
00570
00571
case CNT:
00572 fwrite(cop,
00573
IdCntDimX(cop)*
IdSizeOfType(cop),1,fp);
00574
break;
00575
00576
case IMA:
00577
for(i=0;i<
IdImaDimY(imp);i++)
00578 fwrite(imp[i],
00579
IdImaDimX(imp)*
IdSizeOfType(imp),1,fp);
00580
break;
00581
00582
case VOL:
00583
for(j=0;j<
IdVolDimZ(vop);j++)
00584
for(i=0;i<
IdVolDimY(vop);i++)
00585 fwrite(vop[j][i],
00586
IdVolDimX(vop)*
IdSizeOfType(vop),1,fp);
00587
break;
00588
00589
case SEQ:
00590
for(k=0;k<
IdSeqDimT(sep);k++)
00591
for(j=0;j<
IdSeqDimZ(sep);j++)
00592
for(i=0;i<
IdSeqDimY(sep);i++)
00593 fwrite(sep[k][j][i],
00594
IdVolDimX(sep)*
IdSizeOfType(sep),1,fp);
00595
break;
00596
default:
00597
IdPrintf(
"TypeObj %x non supporte par LibIDO\n",
IdLibidoType(obj));
00598
00599 }
00600 }
00601 p=p->
suiv;
00602 }
00603 fclose(fp);
00604 }
00605
00606
else {
00607
IdErrno =
IDERR_OPEN_FAILED;
00608 retval = 0;
00609 }
00610
00611
00612 header=
_IdAcrFreeListe(header);
00613
00614
00615
00616
00617
00618
00619
00620
00621
return(retval);
00622 }
00623
00624
00625
00626
00627
00628
00629
00630
00631
00632
00633
00634
00635
00636
00637
00638
00639
00640
00641
00642
00643
00644
00645
00646
00647
00648
00649
00650
00651
00652
00653
00654
00655
00656
00657
00658
00669 int IdAcrWriteListeAcr(
char *nf,
ACR_ELEMENT *header)
00670
00671
00672
00673
00674
00675
00676 {
00677 FILE *fp;
00678
00679
short int sh;
00680 uint32_t i32;
00681
ACR_ELEMENT *p;
00682
00683
int retval=1;
00684
00685
_IdAcrAdjustLength(header);
00686
00687
if (
DEBUG) {
IdPrintf (
"\n\nApres _IdAcrAdjustLength\n");
00688
_IdAcrAffListe(header);
00689
IdPrintf (
"\n\n\n");
00690 }
00691
00692
if((fp=fopen(nf,
ID_WFILE_BIN))!=0)
00693 {
00694 p=header;
00695
while(p)
00696 {
00697
switch(p->
type)
00698 {
00699
case BI:
00700 fwrite(p,8,1,fp);
00701 sh=(
short int)p->
value;
00702 fwrite(&sh,2,1,fp);
00703
00704
if(
DEBUG)
IdPrintf(
" BI element (en %p) %4x,%4x long.=%d valeur ",
00705 p, p->
group,p->
number,p->
length);
00706
if (
DEBUG)
IdPrintf(
"%d\n",p->
value);
00707
00708
break;
00709
case BD:
00710 fwrite(p,8,1,fp);
00711 i32=p->
value;
00712 fwrite(&i32,4,1,fp);
00713
00714
if(
DEBUG)
IdPrintf(
" BD element (en %p) %4x,%4x long.=%d valeur ",
00715 p, p->
group,p->
number,p->
length);
00716
if (
DEBUG)
IdPrintf(
"%d\n",p->
value);
00717
break;
00718
case AT:
00719 fwrite(p,8,1,fp);
00720 fwrite((
char *)(p->
value),(size_t)p->
length,1,fp);
00721
00722
if(
DEBUG)
IdPrintf(
" AT element (en %p) %4x,%4x long.=%d valeur ",
00723 p, p->
group,p->
number,p->
length);
00724
if (
DEBUG)
IdPrintf(
"%s\n",(
char *)p->
value);
00725
break;
00726
case AN:
00727 fwrite(p,8,1,fp);
00728 fwrite((
char *)(p->
value),(size_t)p->
length,1,fp);
00729
00730
if(
DEBUG)
IdPrintf(
" AN element (en %p) %4x,%4x long.=%d valeur ",
00731 p, p->
group,p->
number,p->
length);
00732
if (
DEBUG)
IdPrintf(
"%s\n",(
char *)p->
value);
00733
00734
break;
00735
case NO:
00736 fwrite(p,8,1,fp);
00737 fwrite((
char *)(p->
value),(size_t)p->
length,1,fp);
00738
00739
00740
if(
DEBUG)
IdPrintf(
" NO element (en %p) %4x,%4x long.=%d valeur ",
00741 p, p->
group,p->
number,p->
length);
00742
if (
DEBUG)
IdPrintf(
"%x\n",p->
value);
00743
00744
00745
00746
00747
00748
00749
00750
00751
00752
00753
00754
00755
00756
00757
00758
00759
00760
00761
00762
00763
00764
00765
00766
00767
00768
00769
00770
00771
00772
00773
00774
00775
00776
00777
00778
00779
00780
00781 }
00782 p=p->
suiv;
00783 }
00784 fclose(fp);
00785 }
00786
00787
else {
00788
IdErrno =
IDERR_OPEN_FAILED;
00789 retval = 0;
00790 }
00791
00792
return(retval);
00793 }
Generated on Wed Oct 19 09:28:32 2005 for SIMRI3D by
1.3.7