[Castor-users] Writing .Cdf file

[Konstantinos P. Chatzipapas, MSc]

Hello, 

It worked, thank you very much. I have one more question. 

I created the files I am attaching. There is something wrong again. The
reconstructed image is empty.. All black.... 

Can you help me on this? 

Many thanks, 

Konstantinos 

Στις 15-05-2017 18:24, tmerlin έγραψε: 

> Hi Konstantinos, 
> 
> The out of range projection index errors usually means the code can't correctly read the binary datafile as you said. Judging from the values in your cdf file, I think something went wrong during your conversion. Expected types are float for the data, and unsigned integers (4 bytes) for time & ids. I suggest trying something like this in your program: 
> 
> uint32_t id1, id2,** timeMS; float** data; 
> 
> // recover data 
> 
> for (id1=0; id1<bins_elts; id1++)
> {
> for (id2=0; id2<nCrystalsTot; id2++)
> {
> fwrite((char*)(&timeMS[id1][id2]) , sizeof(uint32_t), 1, pFile);
> fwrite((char*)(&data[id1][id2]) , sizeof(float), 1, pFile);
> fwrite((char*)(&id1) , sizeof(uint32_t), 1, pFile);
> fwrite((char*)(&id2) , sizeof(uint32_t), 1, pFile);
> }
> } 
> 
> Besides, I am not sure why your datafile header states 188811 events, while I would expect 107928 events from your numbers (nb of projections * nb of pixels). 
> 
> Best regards, 
> 
> Thibaut 
> 
> On 15/05/2017 12:55, Konstantinos P. Chatzipapas, MSc wrote: 
> 
>> Dear users/developers, 
>> 
>> I am trying to create a program that will write experimental results to the appropriate .Cdf ascii format. I am using this loop: 
>> 
>> char buffer [20];
>> sprintf(buffer, "CDF_CstrProj.Cdf");
>> pFile = fopen (buffer, "w");
>> for (id1=0; id1<bins_elts; id1++)
>> {
>> for (id2=0; id2<nCrystalsTot; id2++)
>> {
>> fprintf(pFile,"%6d ",timeMS[id1][id2]);
>> fprintf(pFile,"%6d ",data[id1][id2]);
>> fprintf(pFile,"%6d ",id1);
>> fprintf(pFile,"%6d ",id2);
>> }
>> } 
>> 
>> Where id1 is the projectionID (0-35) and id2 is the pixelID (0-2997). Data are the counts. I am using integers in data but I have also try it with floats. Then by using a .Cdh I try to reconstruct. There is a problem there. Probably castor-recon does not recognize the file correct, because I am getting this fault repeatedly (the "Projection index" only change: 
>> 
>> ***** oIterativeAlgorithm::StepInnerLoopInsideSubsetLoop() -> A problem occured while computing the projection line !
>> ***** iScannerSPECTConv::GetPositionsAndOrientations() -> Projection index (538976288) out of range [0:35] !
>> ***** vProjector::Project() -> A problem occured while getting positions and orientations from scanner !
>> ***** oProjectorManager::ComputeProjectionLine() -> A problem occured while forward projecting an event ! 
>> 
>> I am attaching the file in order that you take a look and help me. 
>> 
>> Thank you in advance, 
>> 
>> Konstantinos 
>> 
>> -- 
>> 
>> -------------------------
>> Konstantinos P. Chatzipapas 
>> Physicist, MSc 
>> Faculty of Medical Physics 
>> UNIVERSITY OF PATRAS 
>> 
>> web page: _Physics Magazine [2]_ 
>> 
>> _______________________________________________
>> Castor-users mailing list
>> Castor-users at lists.castor-project.org
>> http://lists.castor-project.org/listinfo/castor-users [1]
> 
> -- 
> Thibaut MERLIN -- PhD
> 
> Docteur en Imagerie Médicale au Laboratoire de Traitement de l'Information Medicale (LaTIM - INSERM UMR 1101)
> Bâtiment 1
> CHRU Morvan - 2, Av. Foch
> 29609 Brest CEDEX - FRANCE
> Tel: 06.75.12.24.90
> 
> _______________________________________________
> Castor-users mailing list
> Castor-users at lists.castor-project.org
> http://lists.castor-project.org/listinfo/castor-users [1]
 

Links:
------
[1] http://lists.castor-project.org/listinfo/castor-users
[2] http://www.physicsmag.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.castor-project.org/pipermail/castor-users/attachments/20170515/14523e57/attachment-0001.html>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: header.Cdh
URL: <http://lists.castor-project.org/pipermail/castor-users/attachments/20170515/14523e57/attachment-0001.ksh>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: data.Cdf
Type: application/octet-stream
Size: 1726272 bytes
Desc: not available
URL: <http://lists.castor-project.org/pipermail/castor-users/attachments/20170515/14523e57/attachment-0001.obj>