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

Revision 5636, 64.8 kB (checked in by robert, 8 years ago)

Removed deprecated GUIEventHandler method

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