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

Revision 8777, 65.8 kB (checked in by robert, 6 years ago)

Added --mip command line option which enables Maximum Intensity Projection filtering

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