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

Revision 6941, 64.9 kB (checked in by robert, 7 years ago)

From Martin Lavery and Robert Osfield, Updated examples to use a variation of the MIT License

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