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
#include "output.h"
00037
00045 void WriteVolRF(
PPPVOLUME_COMPLEX_DOUBLE volrf)
00046 {
00047
PPPVOLUME_DOUBLE re,im;
00048
int i,j,k;
00049
int X,Y,Z;
00050
00051 X=
IdVolDimX(volrf);
00052 Y=
IdVolDimY(volrf);
00053 Z=
IdVolDimZ(volrf);
00054
00055 re = (
PPPVOLUME_DOUBLE)
IdVolAlloc(X,Y,Z,
VOL_DOUBLE);
00056 im = (
PPPVOLUME_DOUBLE)
IdVolAlloc(X,Y,Z,
VOL_DOUBLE);
00057
for (k=0;k<Z;k++)
for (i=0;i<Y;i++)
for (j=0;j<X;j++)
00058 {
00059 re[k][i][j] = volrf[k][i][j].
re;
00060 im[k][i][j] = volrf[k][i][j].im;
00061 }
00062
IdVolWriteRawFile(
"rf_re.vol",(
PPPVOLUME)re);
00063
IdVolWriteRawFile(
"rf_im.vol",(
PPPVOLUME)im);
00064
00065
IdVolFree(re);
00066
IdVolFree(im);
00067
00068 }
00069
00078 void WriteVolRecUchar(PPPVOLUME_FLOAT volrec,
char * volname)
00079 {
00080
int X,Y,Z;
00081
PPPVOLUME_UCHAR voluchar;
00082
PPIMAGE_UCHAR imauchar;
00083
double min=0.0;
00084
double max=0.0;
00085
double dyn;
00086
int i,j,k;
00087
00088
IdVolWriteRawFile(
"recfloat.raw",(
PPPVOLUME)volrec);
00089 X=
IdVolDimX(volrec);
00090 Y=
IdVolDimY(volrec);
00091 Z=
IdVolDimZ(volrec);
00092
00093
IdVolMinMax((
PPPVOLUME)volrec,&
min,&
max);
00094 dyn=(
max-
min);
00095
if (Z==1)
00096 {
00097 imauchar = (
PPIMAGE_UCHAR)
IdImaAlloc(X,Y,
IMA_UCHAR);
00098
for (i=0;i<Y;i++)
for (j=0;j<X;j++)
00099 imauchar[i][j]=(
unsigned char)(((volrec[0][i][j]-
min)/dyn)*255.0+0.5);
00100
IdImaWriteRawFile(volname,(
PPIMAGE)imauchar);
00101
IdImaFree(imauchar);
00102 }
00103
else
00104 {
00105 voluchar = (
PPPVOLUME_UCHAR)
IdVolAlloc(X,Y,Z,
VOL_UCHAR);
00106
for (k=0;k<Z;k++)
for (i=0;i<Y;i++)
for (j=0;j<X;j++)
00107 voluchar[k][i][j]=(
unsigned char)(((volrec[k][i][j]-
min)/dyn)*255.0+0.5);
00108
IdVolWriteRawFile(volname,(
PPPVOLUME)voluchar);
00109
IdVolFree(voluchar);
00110 }
00111 }
00112
00120 void WriteSignalRFModule(
EXPERIENCE3D * expr,
char * volname)
00121 {
00122
PPPVOLUME_FLOAT rfmod;
00123
00124 rfmod = (
PPPVOLUME_FLOAT)
GetSignalRFModuleFromExperience(expr);
00125
IdVolWriteRawFile(volname,(
PPPVOLUME_FLOAT)rfmod);
00126
IdVolFree(rfmod);
00127 }
00128
00136 void WriteRecIFFTVolumeModule(
EXPERIENCE3D * expr,
char * volname)
00137 {
00138
PPPVOLUME_FLOAT module;
00139
00140 module= (
PPPVOLUME_FLOAT)
RecVolIFFTModule(expr);
00141
IdVolWriteRawFile(volname,module);
00142
00143
IdVolFree(module);
00144 }
00145
00153 void WriteRecIFFTVolumePhase(
EXPERIENCE3D * expr,
char * volname)
00154 {
00155
PPPVOLUME_FLOAT phase;
00156
00157 phase= (
PPPVOLUME_FLOAT)
RecVolIFFTPhase(expr);
00158
IdVolWriteRawFile(volname,phase);
00159
00160
IdVolFree(phase);
00161 }
00162
00171 void MagstateWatch(
EXPERIENCE3D* expr,
char * filename)
00172 {
00173
int g,ii;
00174 FILE *fic;
00175
int x,y;
00176
00177
x=expr->
nfovx/2;
00178 y=expr->
nfovy/2;
00179 fic=fopen(filename,
"w");
00180
for(g=0;g<expr->
nfovz;g++)
00181 {
00182
for(ii=0;ii<=2;ii++)
00183 {
00184 fprintf(fic,
"%lf ",expr->
magstate[ii][
x][y][g]);
00185 }
00186 fprintf(fic,
"%lf ",sqrt(expr->
magstate[0][
x][y][g]*expr->
magstate[0][
x][y][g]+expr->
magstate[1][
x][y][g]*expr->
magstate[1][
x][y][g]));
00187 fprintf(fic,
"%lf ",atan2(expr->
magstate[1][
x][y][g],expr->
magstate[0][
x][y][g])*180/M_PI);
00188 fprintf(fic,
"\n");
00189 }
00190 fclose(fic);
00191 }