root/OpenSceneGraph/trunk/examples/osgvolume/osgvolume.cpp @ 5962

Revision 5962, 63.9 kB (checked in by robert, 8 years ago)

Ported following examples to osgViewer:

osgteapot
osgtessellate
osgtext
osgtexture1D
osgtexture2D
osgtexture3D
osgtexturerectangle
osgvertexprogram
osgvolume
osgwindows
osgsimple
osgkeyboardmouse

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
Line 
1#include <osg/Node>
2#include <osg/Geometry>
3#include <osg/Notify>
4#include <osg/Texture3D>
5#include <osg/TexGen>
6#include <osg/Geode>
7#include <osg/Billboard>
8#include <osg/PositionAttitudeTransform>
9#include <osg/ClipNode>
10#include <osg/AlphaFunc>
11#include <osg/TexGenNode>
12#include <osg/TexEnv>
13#include <osg/TexEnvCombine>
14#include <osg/Material>
15#include <osg/PrimitiveSet>
16#include <osg/Endian>
17
18#include <osgDB/Registry>
19#include <osgDB/ReadFile>
20#include <osgDB/WriteFile>
21#include <osgDB/FileUtils>
22#include <osgDB/FileNameUtils>
23
24#include <osgGA/EventVisitor>
25
26#include <osgUtil/CullVisitor>
27
28#include <osgViewer/Viewer>
29
30#include <iostream>
31
32
33typedef std::vector< osg::ref_ptr<osg::Image> > ImageList;
34
35//  example ReadOperator
36// struct ReadOperator
37// {
38//     inline void luminance(float l) const { rgba(l,l,l,1.0f); }
39//     inline void alpha(float a) const { rgba(1.0f,1.0f,1.0f,a); }
40//     inline void luminance_alpha(float l,float a) const { rgba(l,l,l,a); }
41//     inline void rgb(float r,float g,float b) const { rgba(r,g,b,1.0f); }
42//     inline void rgba(float r,float g,float b,float a) const { std::cout<<"pixel("<<r<<", "<<g<<", "<<b<<", "<<a<<")"<<std::endl; }
43// };
44
45
46template <typename T, class O>   
47void _readRow(unsigned int num, GLenum pixelFormat, T* data,float scale, const O& operation)
48{
49    switch(pixelFormat)
50    {
51        case(GL_LUMINANCE):         { for(unsigned int i=0;i<num;++i) { float l = float(*data++)*scale; operation.luminance(l); } }  break;
52        case(GL_ALPHA):             { for(unsigned int i=0;i<num;++i) { float a = float(*data++)*scale; operation.alpha(a); } }  break;
53        case(GL_LUMINANCE_ALPHA):   { for(unsigned int i=0;i<num;++i) { float l = float(*data++)*scale; float a = float(*data++)*scale; operation.luminance_alpha(l,a); } }  break;
54        case(GL_RGB):               { for(unsigned int i=0;i<num;++i) { float r = float(*data++)*scale; float g = float(*data++)*scale; float b = float(*data++)*scale; operation.rgb(r,g,b); } }  break;
55        case(GL_RGBA):              { for(unsigned int i=0;i<num;++i) { float r = float(*data++)*scale; float g = float(*data++)*scale; float b = float(*data++)*scale; float a = float(*data++)*scale; operation.rgba(r,g,b,a); } }  break;
56        case(GL_BGR):               { for(unsigned int i=0;i<num;++i) { float b = float(*data++)*scale; float g = float(*data++)*scale; float r = float(*data++)*scale; operation.rgb(r,g,b); } }  break;
57        case(GL_BGRA):              { for(unsigned int i=0;i<num;++i) { float b = float(*data++)*scale; float g = float(*data++)*scale; float r = float(*data++)*scale; float a = float(*data++)*scale; operation.rgba(r,g,b,a); } }  break;
58    }
59}
60
61template <class O>   
62void readRow(unsigned int num, GLenum pixelFormat, GLenum dataType, unsigned char* data, const O& operation)
63{
64    switch(dataType)
65    {
66        case(GL_BYTE):              _readRow(num,pixelFormat, (char*)data,            1.0f/128.0f,        operation); break;
67        case(GL_UNSIGNED_BYTE):     _readRow(num,pixelFormat, (unsigned char*)data,   1.0f/255.0f,        operation); break;
68        case(GL_SHORT):             _readRow(num,pixelFormat, (short*) data,          1.0f/32768.0f,      operation); break;
69        case(GL_UNSIGNED_SHORT):    _readRow(num,pixelFormat, (unsigned short*)data,  1.0f/65535.0f,      operation); break;
70        case(GL_INT):               _readRow(num,pixelFormat, (int*) data,            1.0f/2147483648.0f, operation); break;
71        case(GL_UNSIGNED_INT):      _readRow(num,pixelFormat, (unsigned int*) data,   1.0f/4294967295.0f, operation); break;
72        case(GL_FLOAT):             _readRow(num,pixelFormat, (float*) data,          1.0f,               operation); break;
73    }
74}
75
76template <class O>   
77void readImage(osg::Image* image, const O& operation)
78{
79    if (!image) return;
80   
81    for(int r=0;r<image->r();++r)
82    {
83        for(int t=0;t<image->t();++t)
84        {
85            readRow(image->s(), image->getPixelFormat(), image->getDataType(), image->data(0,t,r), operation);
86        }
87    }
88}
89
90//  example ModifyOperator
91// struct ModifyOperator
92// {
93//     inline void luminance(float& l) const {}
94//     inline void alpha(float& a) const {}
95//     inline void luminance_alpha(float& l,float& a) const {}
96//     inline void rgb(float& r,float& g,float& b) const {}
97//     inline void rgba(float& r,float& g,float& b,float& a) const {}
98// };
99
100
101template <typename T, class M>   
102void _modifyRow(unsigned int num, GLenum pixelFormat, T* data,float scale, const M& operation)
103{
104    float inv_scale = 1.0f/scale;
105    switch(pixelFormat)
106    {
107        case(GL_LUMINANCE):         { for(unsigned int i=0;i<num;++i) { float l = float(*data)*scale; operation.luminance(l); *data++ = T(l*inv_scale); } }  break;
108        case(GL_ALPHA):             { for(unsigned int i=0;i<num;++i) { float a = float(*data)*scale; operation.alpha(a); *data++ = T(a*inv_scale); } }  break;
109        case(GL_LUMINANCE_ALPHA):   { for(unsigned int i=0;i<num;++i) { float l = float(*data)*scale; float a = float(*(data+1))*scale; operation.luminance_alpha(l,a); *data++ = T(l*inv_scale); *data++ = T(a*inv_scale); } }  break;
110        case(GL_RGB):               { for(unsigned int i=0;i<num;++i) { float r = float(*data)*scale; float g = float(*(data+1))*scale; float b = float(*(data+2))*scale; operation.rgb(r,g,b); *data++ = T(r*inv_scale); *data++ = T(g*inv_scale); *data++ = T(b*inv_scale); } }  break;
111        case(GL_RGBA):              { for(unsigned int i=0;i<num;++i) { float r = float(*data)*scale; float g = float(*(data+1))*scale; float b = float(*(data+2))*scale; float a = float(*(data+3))*scale; operation.rgba(r,g,b,a); *data++ = T(r*inv_scale); *data++ = T(g*inv_scale); *data++ = T(g*inv_scale); *data++ = T(a*inv_scale); } }  break;
112        case(GL_BGR):               { for(unsigned int i=0;i<num;++i) { float b = float(*data)*scale; float g = float(*(data+1))*scale; float r = float(*(data+2))*scale; operation.rgb(r,g,b); *data++ = T(b*inv_scale); *data++ = T(g*inv_scale); *data++ = T(r*inv_scale); } }  break;
113        case(GL_BGRA):              { for(unsigned int i=0;i<num;++i) { float b = float(*data)*scale; float g = float(*(data+1))*scale; float r = float(*(data+2))*scale; float a = float(*(data+3))*scale; operation.rgba(r,g,b,a); *data++ = T(g*inv_scale); *data++ = T(b*inv_scale); *data++ = T(r*inv_scale); *data++ = T(a*inv_scale); } }  break;
114    }
115}
116
117template <class M>   
118void modifyRow(unsigned int num, GLenum pixelFormat, GLenum dataType, unsigned char* data, const M& operation)
119{
120    switch(dataType)
121    {
122        case(GL_BYTE):              _modifyRow(num,pixelFormat, (char*)data,            1.0f/128.0f,        operation); break;
123        case(GL_UNSIGNED_BYTE):     _modifyRow(num,pixelFormat, (unsigned char*)data,   1.0f/255.0f,        operation); break;
124        case(GL_SHORT):             _modifyRow(num,pixelFormat, (short*) data,          1.0f/32768.0f,      operation); break;
125        case(GL_UNSIGNED_SHORT):    _modifyRow(num,pixelFormat, (unsigned short*)data,  1.0f/65535.0f,      operation); break;
126        case(GL_INT):               _modifyRow(num,pixelFormat, (int*) data,            1.0f/2147483648.0f, operation); break;
127        case(GL_UNSIGNED_INT):      _modifyRow(num,pixelFormat, (unsigned int*) data,   1.0f/4294967295.0f, operation); break;
128        case(GL_FLOAT):             _modifyRow(num,pixelFormat, (float*) data,          1.0f,               operation); break;
129    }
130}
131
132template <class M>   
133void modifyImage(osg::Image* image, const M& operation)
134{
135    if (!image) return;
136   
137    for(int r=0;r<image->r();++r)
138    {
139        for(int t=0;t<image->t();++t)
140        {
141            modifyRow(image->s(), image->getPixelFormat(), image->getDataType(), image->data(0,t,r), operation);
142        }
143    }
144}
145
146struct PassThroughTransformFunction
147{
148    unsigned char operator() (unsigned char c) const { return c; }
149};
150
151
152struct ProcessRow
153{
154    virtual ~ProcessRow() {}
155
156    virtual void operator() (unsigned int num,
157                    GLenum source_pixelFormat, unsigned char* source,
158                    GLenum dest_pixelFormat, unsigned char* dest) const 
159    {
160        switch(source_pixelFormat)
161        {
162        case(GL_LUMINANCE):
163        case(GL_ALPHA):
164            switch(dest_pixelFormat)
165            {
166            case(GL_LUMINANCE):
167            case(GL_ALPHA): A_to_A(num, source, dest); break;
168            case(GL_LUMINANCE_ALPHA): A_to_LA(num, source, dest); break;
169            case(GL_RGB): A_to_RGB(num, source, dest); break;
170            case(GL_RGBA): A_to_RGBA(num, source, dest); break;
171            }
172            break;
173        case(GL_LUMINANCE_ALPHA):
174            switch(dest_pixelFormat)
175            {
176            case(GL_LUMINANCE):
177            case(GL_ALPHA): LA_to_A(num, source, dest); break;
178            case(GL_LUMINANCE_ALPHA): LA_to_LA(num, source, dest); break;
179            case(GL_RGB): LA_to_RGB(num, source, dest); break;
180            case(GL_RGBA): LA_to_RGBA(num, source, dest); break;
181            }
182            break;
183        case(GL_RGB):
184            switch(dest_pixelFormat)
185            {
186            case(GL_LUMINANCE):
187            case(GL_ALPHA): RGB_to_A(num, source, dest); break;
188            case(GL_LUMINANCE_ALPHA): RGB_to_LA(num, source, dest); break;
189            case(GL_RGB): RGB_to_RGB(num, source, dest); break;
190            case(GL_RGBA): RGB_to_RGBA(num, source, dest); break;
191            }
192            break;
193        case(GL_RGBA):
194            switch(dest_pixelFormat)
195            {
196            case(GL_LUMINANCE):
197            case(GL_ALPHA): RGBA_to_A(num, source, dest); break;
198            case(GL_LUMINANCE_ALPHA): RGBA_to_LA(num, source, dest); break;
199            case(GL_RGB): RGBA_to_RGB(num, source, dest); break;
200            case(GL_RGBA): RGBA_to_RGBA(num, source, dest); break;
201            }
202            break;
203        }
204    }
205
206    ///////////////////////////////////////////////////////////////////////////////
207    // alpha sources..   
208    virtual void A_to_A(unsigned int num, unsigned char* source, unsigned char* dest) const
209    {
210        for(unsigned int i=0;i<num;++i)
211        {
212            *dest++ = *source++;
213        }
214    }
215
216    virtual void A_to_LA(unsigned int num, unsigned char* source, unsigned char* dest) const
217    {
218        for(unsigned int i=0;i<num;++i)
219        {
220            *dest++ = *source;
221            *dest++ = *source++;
222        }
223    }
224                   
225    virtual void A_to_RGB(unsigned int num, unsigned char* source, unsigned char* dest) const
226    {
227        for(unsigned int i=0;i<num;++i)
228        {
229            *dest++ = *source;
230            *dest++ = *source;
231            *dest++ = *source++;
232        }
233    }
234
235    virtual void A_to_RGBA(unsigned int num, unsigned char* source, unsigned char* dest) const
236    {
237        for(unsigned int i=0;i<num;++i)
238        {
239            *dest++ = *source;
240            *dest++ = *source;
241            *dest++ = *source;
242            *dest++ = *source++;
243        }
244    }
245
246    ///////////////////////////////////////////////////////////////////////////////
247    // alpha luminiance sources..   
248    virtual void LA_to_A(unsigned int num, unsigned char* source, unsigned char* dest) const
249    {
250        for(unsigned int i=0;i<num;++i)
251        {
252            ++source;
253            *dest++ = *source++;
254        }
255    }
256
257    virtual void LA_to_LA(unsigned int num, unsigned char* source, unsigned char* dest) const
258    {
259        for(unsigned int i=0;i<num;++i)
260        {
261            *dest++ = *source++;
262            *dest++ = *source++;
263        }
264    }
265                   
266    virtual void LA_to_RGB(unsigned int num, unsigned char* source, unsigned char* dest) const
267    {
268        for(unsigned int i=0;i<num;++i)
269        {
270            *dest++ = *source;
271            *dest++ = *source;
272            *dest++ = *source;
273            source+=2;
274        }
275    }
276
277    virtual void LA_to_RGBA(unsigned int num, unsigned char* source, unsigned char* dest) const
278    {
279        for(unsigned int i=0;i<num;++i)
280        {
281            *dest++ = *source;
282            *dest++ = *source;
283            *dest++ = *source++;
284            *dest++ = *source++;
285        }
286    }
287
288    ///////////////////////////////////////////////////////////////////////////////
289    // RGB sources..   
290    virtual void RGB_to_A(unsigned int num, unsigned char* source, unsigned char* dest) const
291    {
292        for(unsigned int i=0;i<num;++i)
293        {
294            unsigned char val = *source;
295            *dest++ = val;
296            source += 3;
297        }
298    }
299
300    virtual void RGB_to_LA(unsigned int num, unsigned char* source, unsigned char* dest) const
301    {
302        for(unsigned int i=0;i<num;++i)
303        {
304            unsigned char val = *source;
305            *dest++ = val;
306            *dest++ = val;
307            source += 3;
308        }
309    }
310                   
311    virtual void RGB_to_RGB(unsigned int num, unsigned char* source, unsigned char* dest) const
312    {
313        for(unsigned int i=0;i<num;++i)
314        {
315            *dest++ = *source++;
316            *dest++ = *source++;
317            *dest++ = *source++;
318        }
319    }
320
321    virtual void RGB_to_RGBA(unsigned int num, unsigned char* source, unsigned char* dest) const
322    {
323        for(unsigned int i=0;i<num;++i)
324        {
325            unsigned char val = *source;
326            *dest++ = *source++;
327            *dest++ = *source++;
328            *dest++ = *source++;
329            *dest++ = val;
330        }
331    }
332
333    ///////////////////////////////////////////////////////////////////////////////
334    // RGBA sources..   
335    virtual void RGBA_to_A(unsigned int num, unsigned char* source, unsigned char* dest) const
336    {
337        for(unsigned int i=0;i<num;++i)
338        {
339            source += 3;
340            *dest++ = *source++;
341        }
342    }
343
344    virtual void RGBA_to_LA(unsigned int num, unsigned char* source, unsigned char* dest) const
345    {
346        for(unsigned int i=0;i<num;++i)
347        {
348            unsigned char val = *source;
349            source += 3;
350            *dest++ = val;
351            *dest++ = *source++;
352        }
353    }
354                   
355    virtual void RGBA_to_RGB(unsigned int num, unsigned char* source, unsigned char* dest) const
356    {
357        for(unsigned int i=0;i<num;++i)
358        {
359            *dest++ = *source++;
360            *dest++ = *source++;
361            *dest++ = *source++;
362            ++source;
363        }
364    }
365
366    virtual void RGBA_to_RGBA(unsigned int num, unsigned char* source, unsigned char* dest) const
367    {
368        for(unsigned int i=0;i<num;++i)
369        {
370            *dest++ = *source++;
371            *dest++ = *source++;
372            *dest++ = *source++;
373            *dest++ = *source++;
374        }
375    }
376};
377
378
379void clampToNearestValidPowerOfTwo(int& sizeX, int& sizeY, int& sizeZ, int s_maximumTextureSize, int t_maximumTextureSize, int r_maximumTextureSize)
380{
381    // compute nearest powers of two for each axis.
382    int s_nearestPowerOfTwo = 1;
383    while(s_nearestPowerOfTwo<sizeX && s_nearestPowerOfTwo<s_maximumTextureSize) s_nearestPowerOfTwo*=2;
384
385    int t_nearestPowerOfTwo = 1;
386    while(t_nearestPowerOfTwo<sizeY && t_nearestPowerOfTwo<t_maximumTextureSize) t_nearestPowerOfTwo*=2;
387
388    int r_nearestPowerOfTwo = 1;
389    while(r_nearestPowerOfTwo<sizeZ && r_nearestPowerOfTwo<r_maximumTextureSize) r_nearestPowerOfTwo*=2;
390
391    sizeX = s_nearestPowerOfTwo;
392    sizeY = t_nearestPowerOfTwo;
393    sizeZ = r_nearestPowerOfTwo;
394}
395
396osg::Image* createTexture3D(ImageList& imageList, ProcessRow& processRow,
397            unsigned int numComponentsDesired,
398            int s_maximumTextureSize,
399            int t_maximumTextureSize,
400            int r_maximumTextureSize )
401{
402    int max_s = 0;
403    int max_t = 0;
404    unsigned int max_components = 0;
405    int total_r = 0;
406    ImageList::iterator itr;
407    for(itr=imageList.begin();
408        itr!=imageList.end();
409        ++itr)
410    {
411        osg::Image* image = itr->get();
412        GLenum pixelFormat = image->getPixelFormat();
413        if (pixelFormat==GL_ALPHA ||
414            pixelFormat==GL_LUMINANCE ||
415            pixelFormat==GL_LUMINANCE_ALPHA ||
416            pixelFormat==GL_RGB ||
417            pixelFormat==GL_RGBA)
418        {
419            max_s = osg::maximum(image->s(), max_s);
420            max_t = osg::maximum(image->t(), max_t);
421            max_components = osg::maximum(osg::Image::computeNumComponents(pixelFormat), max_components);
422            total_r += image->r();
423        }
424        else
425        {
426            osg::notify(osg::NOTICE)<<"Image "<<image->getFileName()<<" has unsuitable pixel format"<< std::hex<< pixelFormat << std::dec << std::endl;
427        }
428    }
429   
430    if (numComponentsDesired!=0) max_components = numComponentsDesired;
431   
432    GLenum desiredPixelFormat = 0;
433    switch(max_components)
434    {
435    case(1):
436        osg::notify(osg::NOTICE)<<"desiredPixelFormat = GL_LUMINANCE" << std::endl;
437        desiredPixelFormat = GL_LUMINANCE;
438        break;
439    case(2):
440        osg::notify(osg::NOTICE)<<"desiredPixelFormat = GL_LUMINANCE_ALPHA" << std::endl;
441        desiredPixelFormat = GL_LUMINANCE_ALPHA;
442        break;
443    case(3):
444        osg::notify(osg::NOTICE)<<"desiredPixelFormat = GL_RGB" << std::endl;
445        desiredPixelFormat = GL_RGB;
446        break;
447    case(4):
448        osg::notify(osg::NOTICE)<<"desiredPixelFormat = GL_RGBA" << std::endl;
449        desiredPixelFormat = GL_RGBA;
450        break;
451    }   
452    if (desiredPixelFormat==0) return 0;
453   
454    // compute nearest powers of two for each axis.
455    int s_nearestPowerOfTwo = 1;
456    while(s_nearestPowerOfTwo<max_s && s_nearestPowerOfTwo<s_maximumTextureSize) s_nearestPowerOfTwo*=2;
457
458    int t_nearestPowerOfTwo = 1;
459    while(t_nearestPowerOfTwo<max_t && t_nearestPowerOfTwo<t_maximumTextureSize) t_nearestPowerOfTwo*=2;
460
461    int r_nearestPowerOfTwo = 1;
462    while(r_nearestPowerOfTwo<total_r && r_nearestPowerOfTwo<r_maximumTextureSize) r_nearestPowerOfTwo*=2;
463
464
465    osg::notify(osg::NOTICE)<<"max image width = "<<max_s<<"  nearest power of two = "<<s_nearestPowerOfTwo<<std::endl;
466    osg::notify(osg::NOTICE)<<"max image height = "<<max_t<<"  nearest power of two = "<<t_nearestPowerOfTwo<<std::endl;
467    osg::notify(osg::NOTICE)<<"max image depth = "<<total_r<<"  nearest power of two = "<<r_nearestPowerOfTwo<<std::endl;
468   
469    // now allocate the 3d texture;
470    osg::ref_ptr<osg::Image> image_3d = new osg::Image;
471    image_3d->allocateImage(s_nearestPowerOfTwo,t_nearestPowerOfTwo,r_nearestPowerOfTwo,
472                            desiredPixelFormat,GL_UNSIGNED_BYTE);
473       
474
475    unsigned int r_offset = (total_r<r_nearestPowerOfTwo) ? r_nearestPowerOfTwo/2 - total_r/2 : 0;
476
477    int curr_dest_r = r_offset;
478
479    // copy across the values from the source images into the image_3d.
480    for(itr=imageList.begin();
481        itr!=imageList.end();
482        ++itr)
483    {
484        osg::Image* image = itr->get();
485        GLenum pixelFormat = image->getPixelFormat();
486        if (pixelFormat==GL_ALPHA ||
487            pixelFormat==GL_LUMINANCE ||
488            pixelFormat==GL_LUMINANCE_ALPHA ||
489            pixelFormat==GL_RGB ||
490            pixelFormat==GL_RGBA)
491        {
492       
493            int num_r = osg::minimum(image->r(), (image_3d->r() - curr_dest_r));
494            int num_t = osg::minimum(image->t(), image_3d->t());
495            int num_s = osg::minimum(image->s(), image_3d->s());
496       
497            unsigned int s_offset_dest = (image->s()<s_nearestPowerOfTwo) ? s_nearestPowerOfTwo/2 - image->s()/2 : 0;
498            unsigned int t_offset_dest = (image->t()<t_nearestPowerOfTwo) ? t_nearestPowerOfTwo/2 - image->t()/2 : 0;
499
500            for(int r=0;r<num_r;++r, ++curr_dest_r)
501            {
502                for(int t=0;t<num_t;++t)
503                {
504                    unsigned char* dest = image_3d->data(s_offset_dest,t+t_offset_dest,curr_dest_r);
505                    unsigned char* source = image->data(0,t,r);
506
507                    processRow(num_s, image->getPixelFormat(), source, image_3d->getPixelFormat(), dest);
508                }
509            }
510        }
511    }
512    return image_3d.release();
513}
514
515
516osg::Image* createNormalMapTexture(osg::Image* image_3d)
517{
518    unsigned int sourcePixelIncrement = 1;
519    unsigned int alphaOffset = 0;
520    switch(image_3d->getPixelFormat())
521    {
522    case(GL_ALPHA):
523    case(GL_LUMINANCE):
524        sourcePixelIncrement = 1;
525        alphaOffset = 0;
526        break;
527    case(GL_LUMINANCE_ALPHA):
528        sourcePixelIncrement = 2;
529        alphaOffset = 1;
530        break;
531    case(GL_RGB):
532        sourcePixelIncrement = 3;
533        alphaOffset = 0;
534        break;
535    case(GL_RGBA):
536        sourcePixelIncrement = 4;
537        alphaOffset = 3;
538        break;
539    default:
540        osg::notify(osg::NOTICE)<<"Source pixel format not support for normal map generation."<<std::endl;
541        return 0;
542    }
543   
544    osg::ref_ptr<osg::Image> normalmap_3d = new osg::Image;
545    normalmap_3d->allocateImage(image_3d->s(),image_3d->t(),image_3d->r(),
546                            GL_RGBA,GL_UNSIGNED_BYTE);
547
548    if (osg::getCpuByteOrder()==osg::LittleEndian) alphaOffset = sourcePixelIncrement-alphaOffset-1;
549
550    for(int r=1;r<image_3d->r()-1;++r)
551    {
552        for(int t=1;t<image_3d->t()-1;++t)
553        {
554            unsigned char* ptr = image_3d->data(1,t,r)+alphaOffset;
555            unsigned char* left = image_3d->data(0,t,r)+alphaOffset;
556            unsigned char* right = image_3d->data(2,t,r)+alphaOffset;
557            unsigned char* above = image_3d->data(1,t+1,r)+alphaOffset;
558            unsigned char* below = image_3d->data(1,t-1,r)+alphaOffset;
559            unsigned char* in = image_3d->data(1,t,r+1)+alphaOffset;
560            unsigned char* out = image_3d->data(1,t,r-1)+alphaOffset;
561
562            unsigned char* destination = (unsigned char*) normalmap_3d->data(1,t,r);
563
564            for(int s=1;s<image_3d->s()-1;++s)
565            {
566
567                osg::Vec3 grad((float)(*left)-(float)(*right),
568                               (float)(*below)-(float)(*above),
569                               (float)(*out) -(float)(*in));
570
571                grad.normalize();
572
573                if (grad.x()==0.0f && grad.y()==0.0f && grad.z()==0.0f)
574                {
575                    grad.set(128.0f,128.0f,128.0f);
576                }
577                else
578                {
579                    grad.x() = osg::clampBetween((grad.x()+1.0f)*128.0f,0.0f,255.0f);
580                    grad.y() = osg::clampBetween((grad.y()+1.0f)*128.0f,0.0f,255.0f);
581                    grad.z() = osg::clampBetween((grad.z()+1.0f)*128.0f,0.0f,255.0f);
582                }
583
584                *(destination++) = (unsigned char)(grad.x()); // scale and bias X.
585                *(destination++) = (unsigned char)(grad.y()); // scale and bias Y.
586                *(destination++) = (unsigned char)(grad.z()); // scale and bias Z.
587
588                *destination++ = *ptr;
589
590                ptr += sourcePixelIncrement;
591                left += sourcePixelIncrement;
592                right += sourcePixelIncrement;
593                above += sourcePixelIncrement;
594                below += sourcePixelIncrement;
595                in += sourcePixelIncrement;
596                out += sourcePixelIncrement;
597            }
598        }
599    }
600   
601    return normalmap_3d.release();
602}
603
604
605
606osg::Node* createCube(float size,float alpha, unsigned int numSlices, float sliceEnd=1.0f)
607{
608
609    // set up the Geometry.
610    osg::Geometry* geom = new osg::Geometry;
611
612    float halfSize = size*0.5f;
613    float y = halfSize;
614    float dy =-size*1.4/(float)(numSlices-1)*sliceEnd;
615
616    //y = -halfSize;
617    //dy *= 0.5;
618
619    osg::Vec3Array* coords = new osg::Vec3Array(4*numSlices);
620    geom->setVertexArray(coords);
621    for(unsigned int i=0;i<numSlices;++i, y+=dy)
622    {
623        (*coords)[i*4+0].set(-halfSize,y,halfSize);
624        (*coords)[i*4+1].set(-halfSize,y,-halfSize);
625        (*coords)[i*4+2].set(halfSize,y,-halfSize);
626        (*coords)[i*4+3].set(halfSize,y,halfSize);
627    }
628   
629    osg::Vec3Array* normals = new osg::Vec3Array(1);
630    (*normals)[0].set(0.0f,-1.0f,0.0f);
631    geom->setNormalArray(normals);
632    geom->setNormalBinding(osg::Geometry::BIND_OVERALL);
633
634    osg::Vec4Array* colors = new osg::Vec4Array(1);
635    (*colors)[0].set(1.0f,1.0f,1.0f,alpha);
636    geom->setColorArray(colors);
637    geom->setColorBinding(osg::Geometry::BIND_OVERALL);
638
639    geom->addPrimitiveSet(new osg::DrawArrays(osg::PrimitiveSet::QUADS,0,coords->size()));
640
641    osg::Billboard* billboard = new osg::Billboard;
642    billboard->setMode(osg::Billboard::POINT_ROT_WORLD);
643    billboard->addDrawable(geom);
644    billboard->setPosition(0,osg::Vec3(0.0f,0.0f,0.0f));
645   
646    return billboard;
647}
648
649class FollowMouseCallback : public osgGA::GUIEventHandler, public osg::StateSet::Callback
650{
651    public:
652   
653        FollowMouseCallback()
654        {
655            _updateTransparency = false;
656            _updateAlphaCutOff = false;
657            _updateSampleDensity = false;
658        }
659
660        FollowMouseCallback(const FollowMouseCallback&,const osg::CopyOp&) {}
661
662        META_Object(osg,FollowMouseCallback);
663
664        virtual void operator() (osg::StateSet* stateset, osg::NodeVisitor* nv)
665        {
666            if (nv->getVisitorType()==osg::NodeVisitor::EVENT_VISITOR)
667            {
668                osgGA::EventVisitor* ev = dynamic_cast<osgGA::EventVisitor*>(nv);
669                if (ev)
670                {
671                    osgGA::GUIActionAdapter* aa = ev->getActionAdapter();
672                    osgGA::EventQueue::Events& events = ev->getEvents();
673                    for(osgGA::EventQueue::Events::iterator itr=events.begin();
674                        itr!=events.end();
675                        ++itr)
676                    {
677                        handle(*(*itr), *aa, stateset, ev);
678                    }
679                }
680            }
681        }
682       
683        virtual bool handle(const osgGA::GUIEventAdapter& ea,osgGA::GUIActionAdapter&, osg::Object* object, osg::NodeVisitor*)
684        {
685            osg::StateSet* stateset = dynamic_cast<osg::StateSet*>(object);
686            if (!stateset) return false;
687           
688            switch(ea.getEventType())
689            {
690                case(osgGA::GUIEventAdapter::MOVE):
691                case(osgGA::GUIEventAdapter::DRAG):
692                {
693                    float v = ea.getY()*0.5f+0.5f;
694                    osg::Uniform* uniform = 0;
695                    if (_updateTransparency && (uniform = stateset->getUniform("transparency"))) uniform->set(v);
696                    if (_updateAlphaCutOff && (uniform = stateset->getUniform("alphaCutOff"))) uniform->set(v);
697                    if (_updateSampleDensity && (uniform = stateset->getUniform("sampleDensity"))) uniform->set(powf(v,5));
698                    break;
699                }
700                case(osgGA::GUIEventAdapter::KEYDOWN):
701                {
702                    if (ea.getKey()=='t') _updateTransparency = true;
703                    if (ea.getKey()=='a') _updateAlphaCutOff = true;
704                    if (ea.getKey()=='d') _updateSampleDensity = true;
705                    break;
706                }
707                case(osgGA::GUIEventAdapter::KEYUP):
708                {
709                    if (ea.getKey()=='t') _updateTransparency = false;
710                    if (ea.getKey()=='a') _updateAlphaCutOff = false;
711                    if (ea.getKey()=='d') _updateSampleDensity = false;
712                    break;
713                }
714                default:
715                    break;
716            }
717            return false;
718        }
719       
720        bool _updateTransparency;
721        bool _updateAlphaCutOff;
722        bool _updateSampleDensity;
723
724};
725
726osg::Node* createShaderModel(osg::ref_ptr<osg::Image>& image_3d, osg::ref_ptr<osg::Image>& /*normalmap_3d*/,
727                       osg::Texture::InternalFormatMode internalFormatMode,
728                       float /*xSize*/, float /*ySize*/, float /*zSize*/,
729                       float /*xMultiplier*/, float /*yMultiplier*/, float /*zMultiplier*/,
730                       unsigned int /*numSlices*/=500, float /*sliceEnd*/=1.0f, float alphaFuncValue=0.02f)
731{
732    osg::Geode* geode = new osg::Geode;
733    osg::StateSet* stateset = geode->getOrCreateStateSet();
734   
735    stateset->setEventCallback(new FollowMouseCallback);
736   
737    stateset->setMode(GL_ALPHA_TEST,osg::StateAttribute::ON);
738
739    // set up the 3d texture itself,
740    // note, well set the filtering up so that mip mapping is disabled,
741    // gluBuild3DMipsmaps doesn't do a very good job of handled the
742    // inbalanced dimensions of the 256x256x4 texture.
743    osg::Texture3D* texture3D = new osg::Texture3D;
744    texture3D->setFilter(osg::Texture3D::MIN_FILTER,osg::Texture3D::LINEAR);
745    texture3D->setFilter(osg::Texture3D::MAG_FILTER,osg::Texture3D::LINEAR);
746    texture3D->setWrap(osg::Texture3D::WRAP_R,osg::Texture3D::CLAMP);
747    texture3D->setWrap(osg::Texture3D::WRAP_S,osg::Texture3D::CLAMP);
748    texture3D->setWrap(osg::Texture3D::WRAP_T,osg::Texture3D::CLAMP);
749    if (image_3d->getPixelFormat()==GL_ALPHA ||
750        image_3d->getPixelFormat()==GL_LUMINANCE)
751    {
752        texture3D->setInternalFormatMode(osg::Texture3D::USE_USER_DEFINED_FORMAT);
753        texture3D->setInternalFormat(GL_INTENSITY);
754    }
755    else
756    {
757        texture3D->setInternalFormatMode(internalFormatMode);
758    }
759
760    texture3D->setImage(image_3d.get());
761
762    stateset->setTextureAttributeAndModes(0,texture3D,osg::StateAttribute::ON);
763
764    osg::Program* program = new osg::Program;
765    stateset->setAttribute(program);
766
767    // get shaders from source
768    std::string vertexShaderFile = osgDB::findDataFile("volume.vert");
769    if (!vertexShaderFile.empty())
770    {
771        program->addShader(osg::Shader::readShaderFile(osg::Shader::VERTEX, vertexShaderFile));
772    }
773    else
774    {
775        char vertexShaderSource[] =
776            "varying vec3 texcoord;\n"
777            "varying vec3 deltaTexCoord;\n"
778            "\n"
779            "void main(void)\n"
780            "{\n"
781            "    texcoord = gl_MultiTexCoord0.xyz;\n"
782            "    gl_Position     = ftransform();  \n"
783            "    deltaTexCoord = normalize(gl_ModelViewMatrixInverse * vec4(0,0,0,1) - gl_Vertex);\n"
784            "}\n";
785
786        osg::Shader* vertex_shader = new osg::Shader(osg::Shader::VERTEX, vertexShaderSource);
787        program->addShader(vertex_shader);
788
789    }
790   
791    std::string fragmentShaderFile = osgDB::findDataFile("volume.frag");
792    if (!fragmentShaderFile.empty())
793    {
794        program->addShader(osg::Shader::readShaderFile(osg::Shader::FRAGMENT, fragmentShaderFile));
795    }
796    else
797    {
798        //////////////////////////////////////////////////////////////////
799        // fragment shader
800        //
801        char fragmentShaderSource[] =
802            "uniform sampler3D baseTexture;\n"
803            "uniform float sampleDensity;\n"
804            "uniform float transparency;\n"
805            "uniform float alphaCutOff;\n"
806            "\n"
807            "varying vec3 deltaTexCoord;\n"
808            "varying vec3 texcoord;\n"
809            "void main(void) \n"
810            "{ \n"
811            "    vec3 deltaTexCoord2 = normalize(deltaTexCoord)*sampleDensity; \n"
812            "\n"
813            "    gl_FragColor = vec4(0.0, 0.0, 0.0, 0.0); \n"
814            "    \n"
815            "    while (texcoord.x>=0.0 && texcoord.x<=1.0 &&\n"
816            "           texcoord.y>=0.0 && texcoord.y<=1.0 &&\n"
817            "           texcoord.z>=0.0 && texcoord.z<=1.0)\n"
818            "    {\n"
819            "       vec4 color = texture3D( baseTexture, texcoord);\n"
820            "       float r = color[3]*transparency;\n"
821            "       if (r>alphaCutOff)\n"
822            "       {\n"
823            "         gl_FragColor.xyz = gl_FragColor.xyz*(1.0-r)+color.xyz*r;\n"
824            "         gl_FragColor.w += r;\n"
825            "       }\n"
826            "       texcoord += deltaTexCoord2; \n"
827            "    }\n"
828            "    if (gl_FragColor.w>1.0) gl_FragColor.w = 1.0; \n"
829            "}\n";
830
831        osg::Shader* fragment_shader = new osg::Shader(osg::Shader::FRAGMENT, fragmentShaderSource);
832        program->addShader(fragment_shader);
833    }
834
835    osg::Uniform* baseTextureSampler = new osg::Uniform("baseTexture",0);
836    stateset->addUniform(baseTextureSampler);
837
838    osg::Uniform* sampleDensity = new osg::Uniform("sampleDensity", 0.01f);
839    stateset->addUniform(sampleDensity);
840
841    osg::Uniform* transpancy = new osg::Uniform("transparency",0.5f);
842    stateset->addUniform(transpancy);
843
844    osg::Uniform* alphaCutOff = new osg::Uniform("alphaCutOff",alphaFuncValue);
845    stateset->addUniform(alphaCutOff);
846
847    stateset->setMode(GL_CULL_FACE, osg::StateAttribute::ON);
848
849    {
850        osg::Geometry* geom = new osg::Geometry;
851
852        osg::Vec3Array* coords = new osg::Vec3Array(8);
853        (*coords)[0].set(0,0,0);
854        (*coords)[1].set(1,0,0);
855        (*coords)[2].set(1,1,0);
856        (*coords)[3].set(0,1,0);
857        (*coords)[4].set(0,0,1);
858        (*coords)[5].set(1,0,1);
859        (*coords)[6].set(1,1,1);
860        (*coords)[7].set(0,1,1);
861        geom->setVertexArray(coords);
862
863        osg::Vec3Array* tcoords = new osg::Vec3Array(8);
864        (*tcoords)[0].set(0,0,0);
865        (*tcoords)[1].set(1,0,0);
866        (*tcoords)[2].set(1,1,0);
867        (*tcoords)[3].set(0,1,0);
868        (*tcoords)[4].set(0,0,1);
869        (*tcoords)[5].set(1,0,1);
870        (*tcoords)[6].set(1,1,1);
871        (*tcoords)[7].set(0,1,1);
872        geom->setTexCoordArray(0,tcoords);
873
874        osg::Vec4Array* colours = new osg::Vec4Array(1);
875        (*colours)[0].set(1.0f,1.0f,1.0,1.0f);
876        geom->setColorArray(colours);
877        geom->setColorBinding(osg::Geometry::BIND_OVERALL);
878
879        osg::DrawElementsUShort* drawElements = new osg::DrawElementsUShort(GL_QUADS);
880        // bottom
881        drawElements->push_back(0);
882        drawElements->push_back(1);
883        drawElements->push_back(2);
884        drawElements->push_back(3);
885       
886        // bottom
887        drawElements->push_back(3);
888        drawElements->push_back(2);
889        drawElements->push_back(6);
890        drawElements->push_back(7);
891
892        // left
893        drawElements->push_back(0);
894        drawElements->push_back(3);
895        drawElements->push_back(7);
896        drawElements->push_back(4);
897
898        // right
899        drawElements->push_back(5);
900        drawElements->push_back(6);
901        drawElements->push_back(2);
902        drawElements->push_back(1);
903
904        // front
905        drawElements->push_back(1);
906        drawElements->push_back(0);
907        drawElements->push_back(4);
908        drawElements->push_back(5);
909
910        // top
911        drawElements->push_back(7);
912        drawElements->push_back(6);
913        drawElements->push_back(5);
914        drawElements->push_back(4);
915
916        geom->addPrimitiveSet(drawElements);
917
918        geode->addDrawable(geom);
919
920    }
921
922    return geode;
923}
924
925osg::Node* createModel(osg::ref_ptr<osg::Image>& image_3d, osg::ref_ptr<osg::Image>& normalmap_3d,
926                       osg::Texture::InternalFormatMode internalFormatMode,
927                       float xSize, float ySize, float zSize,
928                       float xMultiplier, float yMultiplier, float zMultiplier,
929                       unsigned int numSlices=500, float sliceEnd=1.0f, float alphaFuncValue=0.02f)
930{
931    bool two_pass = normalmap_3d.valid() && (image_3d->getPixelFormat()==GL_RGB || image_3d->getPixelFormat()==GL_RGBA);
932
933    osg::Group* group = new osg::Group;
934   
935    osg::TexGenNode* texgenNode_0 = new osg::TexGenNode;
936    texgenNode_0->setTextureUnit(0);
937    texgenNode_0->getTexGen()->setMode(osg::TexGen::EYE_LINEAR);
938    texgenNode_0->getTexGen()->setPlane(osg::TexGen::S, osg::Plane(xMultiplier,0.0f,0.0f,0.5f));
939    texgenNode_0->getTexGen()->setPlane(osg::TexGen::T, osg::Plane(0.0f,yMultiplier,0.0f,0.5f));
940    texgenNode_0->getTexGen()->setPlane(osg::TexGen::R, osg::Plane(0.0f,0.0f,zMultiplier,0.5f));
941   
942    if (two_pass)
943    {
944        osg::TexGenNode* texgenNode_1 = new osg::TexGenNode;
945        texgenNode_1->setTextureUnit(1);
946        texgenNode_1->getTexGen()->setMode(osg::TexGen::EYE_LINEAR);
947        texgenNode_1->getTexGen()->setPlane(osg::TexGen::S, texgenNode_0->getTexGen()->getPlane(osg::TexGen::S));
948        texgenNode_1->getTexGen()->setPlane(osg::TexGen::T, texgenNode_0->getTexGen()->getPlane(osg::TexGen::T));
949        texgenNode_1->getTexGen()->setPlane(osg::TexGen::R, texgenNode_0->getTexGen()->getPlane(osg::TexGen::R));
950
951        texgenNode_1->addChild(texgenNode_0);
952
953        group->addChild(texgenNode_1);
954    }
955    else
956    { 
957        group->addChild(texgenNode_0);
958    }
959
960    osg::BoundingBox bb(-xSize*0.5f,-ySize*0.5f,-zSize*0.5f,xSize*0.5f,ySize*0.5f,zSize*0.5f);
961
962    osg::ClipNode* clipnode = new osg::ClipNode;
963    clipnode->addChild(createCube(1.0f,1.0f, numSlices,sliceEnd));
964    clipnode->createClipBox(bb);
965
966    {
967        // set up the Geometry to enclose the clip volume to prevent near/far clipping from affecting billboard
968        osg::Geometry* geom = new osg::Geometry;
969
970        osg::Vec3Array* coords = new osg::Vec3Array();
971        coords->push_back(bb.corner(0));
972        coords->push_back(bb.corner(1));
973        coords->push_back(bb.corner(2));
974        coords->push_back(bb.corner(3));
975        coords->push_back(bb.corner(4));
976        coords->push_back(bb.corner(5));
977        coords->push_back(bb.corner(6));
978        coords->push_back(bb.corner(7));
979
980        geom->setVertexArray(coords);
981
982        osg::Vec4Array* colors = new osg::Vec4Array(1);
983        (*colors)[0].set(1.0f,1.0f,1.0f,1.0f);
984        geom->setColorArray(colors);
985        geom->setColorBinding(osg::Geometry::BIND_OVERALL);
986
987        geom->addPrimitiveSet(new osg::DrawArrays(osg::PrimitiveSet::POINTS,0,coords->size()));
988
989        osg::Geode* geode = new osg::Geode;
990        geode->addDrawable(geom);
991       
992        clipnode->addChild(geode);
993       
994    }
995
996    texgenNode_0->addChild(clipnode);
997
998    osg::StateSet* stateset = texgenNode_0->getOrCreateStateSet();
999
1000    stateset->setMode(GL_LIGHTING,osg::StateAttribute::ON);
1001    stateset->setMode(GL_BLEND,osg::StateAttribute::ON);
1002    stateset->setAttribute(new osg::AlphaFunc(osg::AlphaFunc::GREATER,alphaFuncValue));
1003   
1004    osg::Material* material = new osg::Material;
1005    material->setDiffuse(osg::Material::FRONT_AND_BACK,osg::Vec4(1.0f,1.0f,1.0f,1.0f));
1006    stateset->setAttributeAndModes(material);
1007   
1008    osg::Vec3 lightDirection(1.0f,-1.0f,1.0f);
1009    lightDirection.normalize();
1010
1011    if (normalmap_3d.valid())
1012    {
1013        if (two_pass)
1014        {
1015
1016            // set up normal texture
1017            osg::Texture3D* bump_texture3D = new osg::Texture3D;
1018            bump_texture3D->setFilter(osg::Texture3D::MIN_FILTER,osg::Texture3D::LINEAR);
1019            bump_texture3D->setFilter(osg::Texture3D::MAG_FILTER,osg::Texture3D::LINEAR);
1020            bump_texture3D->setWrap(osg::Texture3D::WRAP_R,osg::Texture3D::CLAMP);
1021            bump_texture3D->setWrap(osg::Texture3D::WRAP_S,osg::Texture3D::CLAMP);
1022            bump_texture3D->setWrap(osg::Texture3D::WRAP_T,osg::Texture3D::CLAMP);
1023            bump_texture3D->setImage(normalmap_3d.get());
1024
1025            bump_texture3D->setInternalFormatMode(internalFormatMode);
1026
1027            stateset->setTextureAttributeAndModes(0,bump_texture3D,osg::StateAttribute::ON);
1028
1029            osg::TexEnvCombine* tec = new osg::TexEnvCombine;
1030            tec->setConstantColorAsLightDirection(lightDirection);
1031
1032            tec->setCombine_RGB(osg::TexEnvCombine::DOT3_RGB);
1033            tec->setSource0_RGB(osg::TexEnvCombine::CONSTANT);
1034            tec->setOperand0_RGB(osg::TexEnvCombine::SRC_COLOR);
1035            tec->setSource1_RGB(osg::TexEnvCombine::TEXTURE);
1036            tec->setOperand1_RGB(osg::TexEnvCombine::SRC_COLOR);
1037
1038            tec->setCombine_Alpha(osg::TexEnvCombine::REPLACE);
1039            tec->setSource0_Alpha(osg::TexEnvCombine::PRIMARY_COLOR);
1040            tec->setOperand0_Alpha(osg::TexEnvCombine::SRC_ALPHA);
1041            tec->setSource1_Alpha(osg::TexEnvCombine::TEXTURE);
1042            tec->setOperand1_Alpha(osg::TexEnvCombine::SRC_ALPHA);
1043
1044            stateset->setTextureAttributeAndModes(0, tec, osg::StateAttribute::OVERRIDE|osg::StateAttribute::ON);
1045
1046            stateset->setTextureMode(0,GL_TEXTURE_GEN_S,osg::StateAttribute::ON);
1047            stateset->setTextureMode(0,GL_TEXTURE_GEN_T,osg::StateAttribute::ON);
1048            stateset->setTextureMode(0,GL_TEXTURE_GEN_R,osg::StateAttribute::ON);
1049
1050
1051            // set up color texture
1052            osg::Texture3D* texture3D = new osg::Texture3D;
1053            texture3D->setFilter(osg::Texture3D::MIN_FILTER,osg::Texture3D::LINEAR);
1054            texture3D->setFilter(osg::Texture3D::MAG_FILTER,osg::Texture3D::LINEAR);
1055            texture3D->setWrap(osg::Texture3D::WRAP_R,osg::Texture3D::CLAMP);
1056            texture3D->setWrap(osg::Texture3D::WRAP_S,osg::Texture3D::CLAMP);
1057            texture3D->setWrap(osg::Texture3D::WRAP_T,osg::Texture3D::CLAMP);
1058            if (image_3d->getPixelFormat()==GL_ALPHA ||
1059                image_3d->getPixelFormat()==GL_LUMINANCE)
1060            {
1061                texture3D->setInternalFormatMode(osg::Texture3D::USE_USER_DEFINED_FORMAT);
1062                texture3D->setInternalFormat(GL_INTENSITY);
1063            }
1064            else
1065            {
1066                texture3D->setInternalFormatMode(internalFormatMode);
1067            }
1068            texture3D->setImage(image_3d.get());
1069
1070            stateset->setTextureAttributeAndModes(1,texture3D,osg::StateAttribute::ON);
1071
1072            stateset->setTextureMode(1,GL_TEXTURE_GEN_S,osg::StateAttribute::ON);
1073            stateset->setTextureMode(1,GL_TEXTURE_GEN_T,osg::StateAttribute::ON);
1074            stateset->setTextureMode(1,GL_TEXTURE_GEN_R,osg::StateAttribute::ON);
1075
1076            stateset->setTextureAttributeAndModes(1,new osg::TexEnv(),osg::StateAttribute::ON);
1077
1078        }
1079        else
1080        {
1081            osg::ref_ptr<osg::Image> normalmap_3d = createNormalMapTexture(image_3d.get());
1082            osg::Texture3D* bump_texture3D = new osg::Texture3D;
1083            bump_texture3D->setFilter(osg::Texture3D::MIN_FILTER,osg::Texture3D::LINEAR);
1084            bump_texture3D->setFilter(osg::Texture3D::MAG_FILTER,osg::Texture3D::LINEAR);
1085            bump_texture3D->setWrap(osg::Texture3D::WRAP_R,osg::Texture3D::CLAMP);
1086            bump_texture3D->setWrap(osg::Texture3D::WRAP_S,osg::Texture3D::CLAMP);
1087            bump_texture3D->setWrap(osg::Texture3D::WRAP_T,osg::Texture3D::CLAMP);
1088            bump_texture3D->setImage(normalmap_3d.get());
1089
1090            bump_texture3D->setInternalFormatMode(internalFormatMode);
1091
1092            stateset->setTextureAttributeAndModes(0,bump_texture3D,osg::StateAttribute::ON);
1093
1094            osg::TexEnvCombine* tec = new osg::TexEnvCombine;
1095            tec->setConstantColorAsLightDirection(lightDirection);
1096
1097            tec->setCombine_RGB(osg::TexEnvCombine::DOT3_RGB);
1098            tec->setSource0_RGB(osg::TexEnvCombine::CONSTANT);
1099            tec->setOperand0_RGB(osg::TexEnvCombine::SRC_COLOR);
1100            tec->setSource1_RGB(osg::TexEnvCombine::TEXTURE);
1101            tec->setOperand1_RGB(osg::TexEnvCombine::SRC_COLOR);
1102
1103            tec->setCombine_Alpha(osg::TexEnvCombine::MODULATE);
1104            tec->setSource0_Alpha(osg::TexEnvCombine::PRIMARY_COLOR);
1105            tec->setOperand0_Alpha(osg::TexEnvCombine::SRC_ALPHA);
1106            tec->setSource1_Alpha(osg::TexEnvCombine::TEXTURE);
1107            tec->setOperand1_Alpha(osg::TexEnvCombine::SRC_ALPHA);
1108
1109            stateset->setTextureAttributeAndModes(0, tec, osg::StateAttribute::OVERRIDE|osg::StateAttribute::ON);
1110
1111            stateset->setTextureMode(0,GL_TEXTURE_GEN_S,osg::StateAttribute::ON);
1112            stateset->setTextureMode(0,GL_TEXTURE_GEN_T,osg::StateAttribute::ON);
1113            stateset->setTextureMode(0,GL_TEXTURE_GEN_R,osg::StateAttribute::ON);
1114
1115            image_3d = normalmap_3d;
1116        }
1117    }
1118    else
1119    {     
1120        // set up the 3d texture itself,
1121        // note, well set the filtering up so that mip mapping is disabled,
1122        // gluBuild3DMipsmaps doesn't do a very good job of handled the
1123        // inbalanced dimensions of the 256x256x4 texture.
1124        osg::Texture3D* texture3D = new osg::Texture3D;
1125        texture3D->setFilter(osg::Texture3D::MIN_FILTER,osg::Texture3D::LINEAR);
1126        texture3D->setFilter(osg::Texture3D::MAG_FILTER,osg::Texture3D::LINEAR);
1127        texture3D->setWrap(osg::Texture3D::WRAP_R,osg::Texture3D::CLAMP);
1128        texture3D->setWrap(osg::Texture3D::WRAP_S,osg::Texture3D::CLAMP);
1129        texture3D->setWrap(osg::Texture3D::WRAP_T,osg::Texture3D::CLAMP);
1130        if (image_3d->getPixelFormat()==GL_ALPHA ||
1131            image_3d->getPixelFormat()==GL_LUMINANCE)
1132        {
1133            texture3D->setInternalFormatMode(osg::Texture3D::USE_USER_DEFINED_FORMAT);
1134            texture3D->setInternalFormat(GL_INTENSITY);
1135        }
1136        else
1137        {
1138            texture3D->setInternalFormatMode(internalFormatMode);
1139        }
1140
1141        texture3D->setImage(image_3d.get());
1142
1143        stateset->setTextureAttributeAndModes(0,texture3D,osg::StateAttribute::ON);
1144
1145        stateset->setTextureMode(0,GL_TEXTURE_GEN_S,osg::StateAttribute::ON);
1146        stateset->setTextureMode(0,GL_TEXTURE_GEN_T,osg::StateAttribute::ON);
1147        stateset->setTextureMode(0,GL_TEXTURE_GEN_R,osg::StateAttribute::ON);
1148
1149        stateset->setTextureAttributeAndModes(0,new osg::TexEnv(),osg::StateAttribute::ON);
1150    }
1151 
1152    return group;
1153}
1154
1155struct FindRangeOperator
1156{
1157    FindRangeOperator():
1158        _rmin(FLT_MAX),
1159        _rmax(-FLT_MAX),
1160        _gmin(FLT_MAX),
1161        _gmax(-FLT_MAX),
1162        _bmin(FLT_MAX),
1163        _bmax(-FLT_MAX),
1164        _amin(FLT_MAX),
1165        _amax(-FLT_MAX) {}
1166       
1167    mutable float _rmin, _rmax, _gmin, _gmax, _bmin, _bmax, _amin, _amax;
1168
1169    inline void luminance(float l) const { rgb(l,l,l); }
1170    inline void alpha(float a) const { _amin = osg::minimum(a,_amin); _amax = osg::maximum(a,_amax); }
1171    inline void luminance_alpha(float l,float a) const { rgb(l,l,l); alpha(a); }
1172    inline void rgb(float r,float g,float b) const { _rmin = osg::minimum(r,_rmin); _rmax = osg::maximum(r,_rmax); _gmin = osg::minimum(g,_gmin); _gmax = osg::maximum(g,_gmax); _bmin = osg::minimum(b,_bmin); _bmax = osg::maximum(b,_bmax);  }
1173    inline void rgba(float r,float g,float b,float a) const { rgb(r,g,b); alpha(a); }
1174};
1175 
1176struct ScaleOperator
1177{
1178    ScaleOperator():_scale(1.0f) {}
1179    ScaleOperator(float scale):_scale(scale) {}
1180    ScaleOperator(const ScaleOperator& so):_scale(so._scale) {}
1181   
1182    ScaleOperator& operator = (const ScaleOperator& so) { _scale = so._scale; return *this; }
1183
1184    float _scale;
1185
1186    inline void luminance(float& l) const { l*= _scale; }
1187    inline void alpha(float& a) const { a*= _scale; }
1188    inline void luminance_alpha(float& l,float& a) const { l*= _scale; a*= _scale;  }
1189    inline void rgb(float& r,float& g,float& b) const { r*= _scale; g*=_scale; b*=_scale; }
1190    inline void rgba(float& r,float& g,float& b,float& a) const { r*= _scale; g*=_scale; b*=_scale; a*=_scale; }
1191};
1192
1193struct RecordRowOperator
1194{
1195    RecordRowOperator(unsigned int num):_colours(num),_pos(0) {}
1196
1197    mutable std::vector<osg::Vec4>  _colours;
1198    mutable unsigned int            _pos;
1199   
1200    inline void luminance(float l) const { rgba(l,l,l,1.0f); }
1201    inline void alpha(float a) const { rgba(1.0f,1.0f,1.0f,a); }
1202    inline void luminance_alpha(float l,float a) const { rgba(l,l,l,a);  }
1203    inline void rgb(float r,float g,float b) const { rgba(r,g,b,1.0f); }
1204    inline void rgba(float r,float g,float b,float a) const { _colours[_pos++].set(r,g,b,a); }
1205};
1206
1207struct WriteRowOperator
1208{
1209    WriteRowOperator():_pos(0) {}
1210    WriteRowOperator(unsigned int num):_colours(num),_pos(0) {}
1211
1212    std::vector<osg::Vec4>  _colours;
1213    mutable unsigned int    _pos;
1214   
1215    inline void luminance(float& l) const { l = _colours[_pos++].r(); }
1216    inline void alpha(float& a) const { a = _colours[_pos++].a(); }
1217    inline void luminance_alpha(float& l,float& a) const { l = _colours[_pos].r(); a = _colours[_pos++].a(); }
1218    inline void rgb(float& r,float& g,float& b) const { r = _colours[_pos].r(); g = _colours[_pos].g(); b = _colours[_pos].b(); }
1219    inline void rgba(float& r,float& g,float& b,float& a) const {  r = _colours[_pos].r(); g = _colours[_pos].g(); b = _colours[_pos].b(); a = _colours[_pos++].a(); }
1220};
1221
1222osg::Image* readRaw(int sizeX, int sizeY, int sizeZ, int numberBytesPerComponent, int numberOfComponents, const std::string& endian, const std::string& raw_filename)
1223{
1224    std::ifstream fin(raw_filename.c_str(), std::ifstream::binary);
1225    if (!fin) return 0;
1226
1227    GLenum pixelFormat;
1228    switch(numberOfComponents)
1229    {
1230        case 1 : pixelFormat = GL_LUMINANCE; break;
1231        case 2 : pixelFormat = GL_LUMINANCE_ALPHA; break;
1232        case 3 : pixelFormat = GL_RGB; break;
1233        case 4 : pixelFormat = GL_RGBA; break;
1234        default :
1235            osg::notify(osg::NOTICE)<<"Error: numberOfComponents="<<numberOfComponents<<" not supported, only 1,2,3 or 4 are supported."<<std::endl;
1236            return 0;
1237    }
1238
1239   
1240    GLenum dataType;
1241    switch(numberBytesPerComponent)
1242    {
1243        case 1 : dataType = GL_UNSIGNED_BYTE; break;
1244        case 2 : dataType = GL_UNSIGNED_SHORT; break;
1245        case 4 : dataType = GL_UNSIGNED_INT; break;
1246        default :
1247            osg::notify(osg::NOTICE)<<"Error: numberBytesPerComponent="<<numberBytesPerComponent<<" not supported, only 1,2 or 4 are supported."<<std::endl;
1248            return 0;
1249    }
1250   
1251    int s_maximumTextureSize=256, t_maximumTextureSize=256, r_maximumTextureSize=256;
1252   
1253    int sizeS = sizeX;
1254    int sizeT = sizeY;
1255    int sizeR = sizeZ;
1256    clampToNearestValidPowerOfTwo(sizeS, sizeT, sizeR, s_maximumTextureSize, t_maximumTextureSize, r_maximumTextureSize);
1257
1258    osg::ref_ptr<osg::Image> image = new osg::Image;
1259    image->allocateImage(sizeS, sizeT, sizeR, pixelFormat, dataType);
1260   
1261   
1262    bool endianSwap = (osg::getCpuByteOrder()==osg::BigEndian) ? (endian!="big") : (endian=="big");
1263   
1264    unsigned int r_offset = (sizeZ<sizeR) ? sizeR/2 - sizeZ/2 : 0;
1265   
1266    int offset = endianSwap ? numberBytesPerComponent : 0;
1267    int delta = endianSwap ? -1 : 1;
1268    for(int r=0;r<sizeZ;++r)
1269    {
1270        for(int t=0;t<sizeY;++t)
1271        {
1272            char* data = (char*) image->data(0,t,r+r_offset);
1273            for(int s=0;s<sizeX;++s)
1274            {
1275                if (!fin) return 0;
1276               
1277                for(int c=0;c<numberOfComponents;++c)
1278                {
1279                    char* ptr = data+offset;
1280                    for(int b=0;b<numberBytesPerComponent;++b)
1281                    {
1282                        fin.read((char*)ptr, 1);
1283                        ptr += delta;
1284                    }
1285                    data += numberBytesPerComponent;
1286                }
1287            }
1288        }
1289    }
1290
1291
1292    // normalise texture
1293    {
1294        // compute range of values
1295        FindRangeOperator rangeOp;   
1296        readImage(image.get(), rangeOp);
1297        modifyImage(image.get(),ScaleOperator(1.0f/rangeOp._rmax));
1298    }
1299   
1300   
1301    fin.close();
1302
1303    if (dataType!=GL_UNSIGNED_BYTE)
1304    {
1305        // need to convert to ubyte
1306       
1307        osg::ref_ptr<osg::Image> new_image = new osg::Image;
1308        new_image->allocateImage(sizeS, sizeT, sizeR, pixelFormat, GL_UNSIGNED_BYTE);
1309       
1310        RecordRowOperator readOp(sizeS);
1311        WriteRowOperator writeOp;
1312
1313        for(int r=0;r<sizeR;++r)
1314        {
1315            for(int t=0;t<sizeT;++t)
1316            {
1317                // reset the indices to begining
1318                readOp._pos = 0;
1319                writeOp._pos = 0;
1320           
1321                // read the pixels into readOp's _colour array
1322                readRow(sizeS, pixelFormat, dataType, image->data(0,t,r), readOp);
1323                               
1324                // pass readOp's _colour array contents over to writeOp (note this is just a pointer swap).
1325                writeOp._colours.swap(readOp._colours);
1326               
1327                modifyRow(sizeS, pixelFormat, GL_UNSIGNED_BYTE, new_image->data(0,t,r), writeOp);
1328
1329                // return readOp's _colour array contents back to its rightful owner.
1330                writeOp._colours.swap(readOp._colours);
1331            }
1332        }
1333       
1334        image = new_image;
1335    }
1336   
1337    return image.release();
1338   
1339   
1340}
1341
1342enum ColourSpaceOperation
1343{
1344    NO_COLOUR_SPACE_OPERATION,
1345    MODULATE_ALPHA_BY_LUMINANCE,
1346    MODULATE_ALPHA_BY_COLOUR,
1347    REPLACE_ALPHA_WITH_LUMINACE
1348};
1349
1350struct ModulateAlphaByLuminanceOperator
1351{
1352    ModulateAlphaByLuminanceOperator() {}
1353
1354    inline void luminance(float&) const {}
1355    inline void alpha(float&) const {}
1356    inline void luminance_alpha(float& l,float& a) const { a*= l; }
1357    inline void rgb(float&,float&,float&) const {}
1358    inline void rgba(float& r,float& g,float& b,float& a) const { float l = (r+g+b)*0.3333333; a *= l;}
1359};
1360
1361struct ModulateAlphaByColourOperator
1362{
1363    ModulateAlphaByColourOperator(const osg::Vec4& colour):_colour(colour) { _lum = _colour.length(); }
1364   
1365    osg::Vec4 _colour;
1366    float _lum;
1367
1368    inline void luminance(float&) const {}
1369    inline void alpha(float&) const {}
1370    inline void luminance_alpha(float& l,float& a) const { a*= l*_lum; }
1371    inline void rgb(float&,float&,float&) const {}
1372    inline void rgba(float& r,float& g,float& b,float& a) const { a = (r*_colour.r()+g*_colour.g()+b*_colour.b()+a*_colour.a()); }
1373};
1374
1375struct ReplaceAlphaWithLuminanceOperator
1376{
1377    ReplaceAlphaWithLuminanceOperator() {}
1378
1379    inline void luminance(float&) const {}
1380    inline void alpha(float&) const {}
1381    inline void luminance_alpha(float& l,float& a) const { a= l; }
1382    inline void rgb(float&,float&,float&) const { }
1383    inline void rgba(float& r,float& g,float& b,float& a) const { float l = (r+g+b)*0.3333333; a = l; }
1384};
1385
1386void doColourSpaceConversion(ColourSpaceOperation op, osg::Image* image, osg::Vec4& colour)
1387{
1388    switch(op)
1389    {
1390        case (MODULATE_ALPHA_BY_LUMINANCE):
1391            std::cout<<"doing conversion MODULATE_ALPHA_BY_LUMINANCE"<<std::endl;
1392            modifyImage(image,ModulateAlphaByLuminanceOperator());
1393            break;
1394        case (MODULATE_ALPHA_BY_COLOUR):
1395            std::cout<<"doing conversion MODULATE_ALPHA_BY_COLOUR"<<std::endl;
1396            modifyImage(image,ModulateAlphaByColourOperator(colour));
1397            break;
1398        case (REPLACE_ALPHA_WITH_LUMINACE):
1399            std::cout<<"doing conversion REPLACE_ALPHA_WITH_LUMINACE"<<std::endl;
1400            modifyImage(image,ReplaceAlphaWithLuminanceOperator());
1401            break;
1402        default:
1403            break;
1404    }
1405}
1406
1407int main( int argc, char **argv )
1408{
1409    // use an ArgumentParser object to manage the program arguments.
1410    osg::ArgumentParser arguments(&argc,argv);
1411   
1412    // set up the usage document, in case we need to print out how to use this program.
1413    arguments.getApplicationUsage()->setDescription(arguments.getApplicationName()+" is the example which demonstrates use of 3D textures.");
1414    arguments.getApplicationUsage()->setCommandLineUsage(arguments.getApplicationName()+" [options] filename ...");
1415    arguments.getApplicationUsage()->addCommandLineOption("-h or --help","Display this information");
1416    arguments.getApplicationUsage()->addCommandLineOption("-n","Create normal map for per voxel lighting.");
1417    arguments.getApplicationUsage()->addCommandLineOption("-s <numSlices>","Number of slices to create.");
1418    arguments.getApplicationUsage()->addCommandLineOption("--images [filenames]","Specify a stack of 2d images to build the 3d volume from.");
1419    arguments.getApplicationUsage()->addCommandLineOption("--shader","Use OpenGL Shading Language.");
1420    arguments.getApplicationUsage()->addCommandLineOption("--xSize <size>","Relative width of rendered brick.");
1421    arguments.getApplicationUsage()->addCommandLineOption("--ySize <size>","Relative length of rendered brick.");
1422    arguments.getApplicationUsage()->addCommandLineOption("--zSize <size>","Relative height of rendered brick.");
1423    arguments.getApplicationUsage()->addCommandLineOption("--xMultiplier <multiplier>","Tex coord x mulitplier.");
1424    arguments.getApplicationUsage()->addCommandLineOption("--yMultiplier <multiplier>","Tex coord y mulitplier.");
1425    arguments.getApplicationUsage()->addCommandLineOption("--zMultiplier <multiplier>","Tex coord z mulitplier.");
1426    arguments.getApplicationUsage()->addCommandLineOption("--clip <ratio>","clip volume as a ratio, 0.0 clip all, 1.0 clip none.");
1427    arguments.getApplicationUsage()->addCommandLineOption("--maxTextureSize <size>","Set the texture maximum resolution in the s,t,r (x,y,z) dimensions.");
1428    arguments.getApplicationUsage()->addCommandLineOption("--s_maxTextureSize <size>","Set the texture maximum resolution in the s (x) dimension.");
1429    arguments.getApplicationUsage()->addCommandLineOption("--t_maxTextureSize <size>","Set the texture maximum resolution in the t (y) dimension.");
1430    arguments.getApplicationUsage()->addCommandLineOption("--r_maxTextureSize <size>","Set the texture maximum resolution in the r (z) dimension.");
1431    arguments.getApplicationUsage()->addCommandLineOption("--compressed","Enable the usage of compressed textures.");
1432    arguments.getApplicationUsage()->addCommandLineOption("--compressed-arb","Enable the usage of OpenGL ARB compressed textures.");
1433    arguments.getApplicationUsage()->addCommandLineOption("--compressed-dxt1","Enable the usage of S3TC DXT1 compressed textures.");
1434    arguments.getApplicationUsage()->addCommandLineOption("--compressed-dxt3","Enable the usage of S3TC DXT3 compressed textures.");
1435    arguments.getApplicationUsage()->addCommandLineOption("--compressed-dxt5","Enable the usage of S3TC DXT5 compressed textures.");
1436    arguments.getApplicationUsage()->addCommandLineOption("--modulate-alpha-by-luminance","For each pixel multiple the alpha value by the luminance.");
1437    arguments.getApplicationUsage()->addCommandLineOption("--replace-alpha-with-luminance","For each pixel mSet the alpha value to the luminance.");
1438    arguments.getApplicationUsage()->addCommandLineOption("--num-components <num>","Set the number of components to in he target image.");
1439//    arguments.getApplicationUsage()->addCommandLineOption("--raw <sizeX> <sizeY> <sizeZ> <numberBytesPerComponent> <numberOfComponents> <endian> <filename>","read a raw image data");
1440
1441    // construct the viewer.
1442    osgViewer::Viewer viewer;
1443
1444    // if user request help write it out to cout.
1445    if (arguments.read("-h") || arguments.read("--help"))
1446    {
1447        arguments.getApplicationUsage()->write(std::cout);
1448        return 1;
1449    }
1450
1451    std::string outputFile;
1452    while (arguments.read("-o",outputFile)) {}
1453
1454
1455    unsigned int numSlices=500;
1456    while (arguments.read("-s",numSlices)) {}
1457   
1458   
1459    float sliceEnd=1.0f;
1460    while (arguments.read("--clip",sliceEnd)) {}
1461
1462    float alphaFunc=0.02f;
1463    while (arguments.read("--alphaFunc",alphaFunc)) {}
1464
1465
1466    bool createNormalMap = false;
1467    while (arguments.read("-n")) createNormalMap=true;
1468
1469    float xSize=1.0f, ySize=1.0f, zSize=1.0f;
1470    while (arguments.read("--xSize",xSize)) {}
1471    while (arguments.read("--ySize",ySize)) {}
1472    while (arguments.read("--zSize",zSize)) {}
1473
1474    float xMultiplier=1.0f, yMultiplier=1.0f, zMultiplier=1.0f;
1475    while (arguments.read("--xMultiplier",xMultiplier)) {}
1476    while (arguments.read("--yMultiplier",yMultiplier)) {}
1477    while (arguments.read("--zMultiplier",zMultiplier)) {}
1478
1479    int s_maximumTextureSize = 256;
1480    int t_maximumTextureSize = 256;
1481    int r_maximumTextureSize = 256;
1482    int maximumTextureSize = 256;
1483    while(arguments.read("--maxTextureSize",maximumTextureSize))
1484    {
1485        s_maximumTextureSize = maximumTextureSize;
1486        t_maximumTextureSize = maximumTextureSize;
1487        r_maximumTextureSize = maximumTextureSize;
1488    }
1489    while(arguments.read("--s_maxTextureSize",s_maximumTextureSize)) {}
1490    while(arguments.read("--t_maxTextureSize",t_maximumTextureSize)) {}
1491    while(arguments.read("--r_maxTextureSize",r_maximumTextureSize)) {}
1492
1493    osg::Texture::InternalFormatMode internalFormatMode = osg::Texture::USE_IMAGE_DATA_FORMAT;
1494    while(arguments.read("--compressed") || arguments.read("--compressed-arb")) { internalFormatMode = osg::Texture::USE_ARB_COMPRESSION; }
1495
1496    while(arguments.read("--compressed-dxt1")) { internalFormatMode = osg::Texture::USE_S3TC_DXT1_COMPRESSION; }
1497    while(arguments.read("--compressed-dxt3")) { internalFormatMode = osg::Texture::USE_S3TC_DXT3_COMPRESSION; }
1498    while(arguments.read("--compressed-dxt5")) { internalFormatMode = osg::Texture::USE_S3TC_DXT5_COMPRESSION; }
1499   
1500   
1501    // set up colour space operation.
1502    ColourSpaceOperation colourSpaceOperation = NO_COLOUR_SPACE_OPERATION;
1503    osg::Vec4 colourModulate(0.25f,0.25f,0.25f,0.25f);
1504    while(arguments.read("--modulate-alpha-by-luminance")) { colourSpaceOperation = MODULATE_ALPHA_BY_LUMINANCE; }
1505    while(arguments.read("--modulate-alpha-by-colour", colourModulate.x(),colourModulate.y(),colourModulate.z(),colourModulate.w() )) { colourSpaceOperation = MODULATE_ALPHA_BY_COLOUR; }
1506    while(arguments.read("--replace-alpha-with-luminance")) { colourSpaceOperation = REPLACE_ALPHA_WITH_LUMINACE; }
1507       
1508   
1509    unsigned int numComponentsDesired = 0;
1510    while(arguments.read("--num-components", numComponentsDesired)) {}
1511
1512    bool useShader = false;
1513    while(arguments.read("--shader")) { useShader = true; }
1514
1515    osg::ref_ptr<osg::Image> image_3d;
1516
1517    int sizeX, sizeY, sizeZ, numberBytesPerComponent, numberOfComponents;
1518    std::string endian, raw_filename;
1519    while (arguments.read("--raw", sizeX, sizeY, sizeZ, numberBytesPerComponent, numberOfComponents, endian, raw_filename))
1520    {
1521        image_3d = readRaw(sizeX, sizeY, sizeZ, numberBytesPerComponent, numberOfComponents, endian, raw_filename);
1522    }
1523
1524    while (arguments.read("--images"))
1525    {
1526        ImageList imageList;
1527        for(int pos=1;pos<arguments.argc() && !arguments.isOption(pos);++pos)
1528        {
1529            // not an option so assume string is a filename.
1530            osg::Image *image = osgDB::readImageFile( arguments[pos]);
1531
1532            if(image)
1533            {
1534                imageList.push_back(image);
1535            }
1536        }
1537       
1538        // pack the textures into a single texture.
1539        ProcessRow processRow;
1540        image_3d = createTexture3D(imageList, processRow, numComponentsDesired, s_maximumTextureSize, t_maximumTextureSize, r_maximumTextureSize);
1541    }
1542
1543
1544    // any option left unread are converted into errors to write out later.
1545    arguments.reportRemainingOptionsAsUnrecognized();
1546
1547    // report any errors if they have occured when parsing the program aguments.
1548    if (arguments.errors())
1549    {
1550        arguments.writeErrorMessages(std::cout);
1551        return 1;
1552    }
1553
1554    // assume remaining argments are file names of textures.
1555    for(int pos=1;pos<arguments.argc() && !image_3d;++pos)
1556    {
1557        if (!arguments.isOption(pos))
1558        {
1559            // not an option so assume string is a filename.
1560            image_3d = osgDB::readImageFile( arguments[pos]);
1561        }
1562    }
1563   
1564    if (!image_3d) return 0;
1565   
1566    if (colourSpaceOperation!=NO_COLOUR_SPACE_OPERATION)
1567    {
1568        doColourSpaceConversion(colourSpaceOperation, image_3d.get(), colourModulate);
1569    }
1570   
1571    osg::ref_ptr<osg::Image> normalmap_3d = createNormalMap ? createNormalMapTexture(image_3d.get()) : 0;
1572
1573
1574
1575    // create a model from the images.
1576    osg::Node* rootNode = 0;
1577   
1578    if (useShader)
1579    {
1580        rootNode = createShaderModel(image_3d, normalmap_3d,
1581                               internalFormatMode,
1582                               xSize, ySize, zSize,
1583                               xMultiplier, yMultiplier, zMultiplier,
1584                               numSlices, sliceEnd, alphaFunc);
1585    }
1586    else
1587    {
1588        rootNode = createModel(image_3d, normalmap_3d,
1589                               internalFormatMode,
1590                               xSize, ySize, zSize,
1591                               xMultiplier, yMultiplier, zMultiplier,
1592                               numSlices, sliceEnd, alphaFunc);
1593    }
1594   
1595    if (!outputFile.empty())
1596    {   
1597        std::string ext = osgDB::getFileExtension(outputFile);
1598        std::string name_no_ext = osgDB::getNameLessExtension(outputFile);
1599        if (ext=="osg")
1600        {
1601            if (image_3d.valid())
1602            {
1603                image_3d->setFileName(name_no_ext + ".dds");           
1604                osgDB::writeImageFile(*image_3d, image_3d->getFileName());
1605            }
1606            if (normalmap_3d.valid())
1607            {
1608                normalmap_3d->setFileName(name_no_ext + "_normalmap.dds");           
1609                osgDB::writeImageFile(*normalmap_3d, normalmap_3d->getFileName());
1610            }
1611           
1612            osgDB::writeNodeFile(*rootNode, outputFile);
1613        }
1614        else if (ext=="ive")
1615        {
1616            osgDB::writeNodeFile(*rootNode, outputFile);       
1617        }
1618        else if (ext=="dds")
1619        {
1620            osgDB::writeImageFile(*image_3d, outputFile);       
1621        }
1622        else
1623        {
1624            std::cout<<"Extension not support for file output, not file written."<<std::endl;
1625        }
1626       
1627        return 0;
1628    }
1629
1630
1631    if (rootNode)
1632    {
1633
1634        // set the scene to render
1635        viewer.setSceneData(rootNode);
1636       
1637        // the the viewers main frame loop
1638        viewer.run();
1639    }   
1640   
1641    return 0;
1642
1643}
Note: See TracBrowser for help on using the browser.