32 #include "dcmtk/dcmdata/dcfilefo.h"
33 #include "dcmtk/dcmdata/dcdeftag.h"
34 #include "dcmtk/dcmdata/dcpixel.h"
35 #include "dcmtk/dcmdata/dcpxitem.h"
36 #include "dcmtk/dcmdata/dcpixseq.h"
37 #include "dcmtk/dcmrt/drtimage.h"
66 for(
int ir = 0; ir <
fNoVoxelY; ir += fCompress ) {
67 for(
int ic = 0; ic <
fNoVoxelX; ic += fCompress ) {
69 int isumrMax =
std::min(ir+fCompress,fNoVoxelY);
70 int isumcMax =
std::min(ic+fCompress,fNoVoxelX);
71 for(
int isumr = ir; isumr < isumrMax; isumr ++ ) {
72 for(
int isumc = ic; isumc < isumcMax; isumc ++ ) {
76 meanHV /= (isumrMax-ir)*(isumcMax-ic);
91 for(
int ir = 0; ir <
fNoVoxelY/fCompress; ir++ ) {
92 for(
int ic = 0; ic <
fNoVoxelX/fCompress; ic++ ) {
94 if( ic !=
fNoVoxelX/fCompress-1) fout <<
" ";
95 if( copyNo%8 == 7 ) fout <<
G4endl;
98 if( copyNo%8 != 0 ) fout <<
G4endl;