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

Revision 8961, 84.3 kB (checked in by robert, 6 years ago)

Updated shaders

  • 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/Texture1D>
24#include <osg/ImageSequence>
25#include <osg/TexGen>
26#include <osg/Geode>
27#include <osg/Billboard>
28#include <osg/PositionAttitudeTransform>
29#include <osg/ClipNode>
30#include <osg/AlphaFunc>
31#include <osg/TexGenNode>
32#include <osg/TexEnv>
33#include <osg/TexEnvCombine>
34#include <osg/Material>
35#include <osg/PrimitiveSet>
36#include <osg/Endian>
37#include <osg/BlendFunc>
38#include <osg/BlendEquation>
39#include <osg/TransferFunction>
40
41#include <osgDB/Registry>
42#include <osgDB/ReadFile>
43#include <osgDB/WriteFile>
44#include <osgDB/FileUtils>
45#include <osgDB/FileNameUtils>
46
47#include <osgGA/EventVisitor>
48#include <osgGA/TrackballManipulator>
49#include <osgGA/FlightManipulator>
50#include <osgGA/KeySwitchMatrixManipulator>
51
52#include <osgUtil/CullVisitor>
53
54#include <osgViewer/Viewer>
55#include <osgViewer/ViewerEventHandlers>
56
57#include <osg/io_utils>
58
59#include <algorithm>
60#include <iostream>
61
62#include <osgVolume/ImageUtils>
63
64typedef std::vector< osg::ref_ptr<osg::Image> > ImageList;
65
66enum ShadingModel
67{
68    Standard,
69    Light,
70    Isosurface,
71    MaximumIntensityProjection
72};
73
74//  example ReadOperator
75// struct ReadOperator
76// {
77//     inline void luminance(float l) const { rgba(l,l,l,1.0f); }
78//     inline void alpha(float a) const { rgba(1.0f,1.0f,1.0f,a); }
79//     inline void luminance_alpha(float l,float a) const { rgba(l,l,l,a); }
80//     inline void rgb(float r,float g,float b) const { rgba(r,g,b,1.0f); }
81//     inline void rgba(float r,float g,float b,float a) const { std::cout<<"pixel("<<r<<", "<<g<<", "<<b<<", "<<a<<")"<<std::endl; }
82// };
83
84
85
86struct PassThroughTransformFunction
87{
88    unsigned char operator() (unsigned char c) const { return c; }
89};
90
91
92struct ProcessRow
93{
94    virtual ~ProcessRow() {}
95
96    virtual void operator() (unsigned int num,
97                    GLenum source_pixelFormat, unsigned char* source,
98                    GLenum dest_pixelFormat, unsigned char* dest) const 
99    {
100        switch(source_pixelFormat)
101        {
102        case(GL_LUMINANCE):
103        case(GL_ALPHA):
104            switch(dest_pixelFormat)
105            {
106            case(GL_LUMINANCE):
107            case(GL_ALPHA): A_to_A(num, source, dest); break;
108            case(GL_LUMINANCE_ALPHA): A_to_LA(num, source, dest); break;
109            case(GL_RGB): A_to_RGB(num, source, dest); break;
110            case(GL_RGBA): A_to_RGBA(num, source, dest); break;
111            }
112            break;
113        case(GL_LUMINANCE_ALPHA):
114            switch(dest_pixelFormat)
115            {
116            case(GL_LUMINANCE):
117            case(GL_ALPHA): LA_to_A(num, source, dest); break;
118            case(GL_LUMINANCE_ALPHA): LA_to_LA(num, source, dest); break;
119            case(GL_RGB): LA_to_RGB(num, source, dest); break;
120            case(GL_RGBA): LA_to_RGBA(num, source, dest); break;
121            }
122            break;
123        case(GL_RGB):
124            switch(dest_pixelFormat)
125            {
126            case(GL_LUMINANCE):
127            case(GL_ALPHA): RGB_to_A(num, source, dest); break;
128            case(GL_LUMINANCE_ALPHA): RGB_to_LA(num, source, dest); break;
129            case(GL_RGB): RGB_to_RGB(num, source, dest); break;
130            case(GL_RGBA): RGB_to_RGBA(num, source, dest); break;
131            }
132            break;
133        case(GL_RGBA):
134            switch(dest_pixelFormat)
135            {
136            case(GL_LUMINANCE):
137            case(GL_ALPHA): RGBA_to_A(num, source, dest); break;
138            case(GL_LUMINANCE_ALPHA): RGBA_to_LA(num, source, dest); break;
139            case(GL_RGB): RGBA_to_RGB(num, source, dest); break;
140            case(GL_RGBA): RGBA_to_RGBA(num, source, dest); break;
141            }
142            break;
143        }
144    }
145
146    ///////////////////////////////////////////////////////////////////////////////
147    // alpha sources..   
148    virtual void A_to_A(unsigned int num, unsigned char* source, unsigned char* dest) const
149    {
150        for(unsigned int i=0;i<num;++i)
151        {
152            *dest++ = *source++;
153        }
154    }
155
156    virtual void A_to_LA(unsigned int num, unsigned char* source, unsigned char* dest) const
157    {
158        for(unsigned int i=0;i<num;++i)
159        {
160            *dest++ = *source;
161            *dest++ = *source++;
162        }
163    }
164                   
165    virtual void A_to_RGB(unsigned int num, unsigned char* source, unsigned char* dest) const
166    {
167        for(unsigned int i=0;i<num;++i)
168        {
169            *dest++ = *source;
170            *dest++ = *source;
171            *dest++ = *source++;
172        }
173    }
174
175    virtual void A_to_RGBA(unsigned int num, unsigned char* source, unsigned char* dest) const
176    {
177        for(unsigned int i=0;i<num;++i)
178        {
179            *dest++ = *source;
180            *dest++ = *source;
181            *dest++ = *source;
182            *dest++ = *source++;
183        }
184    }
185
186    ///////////////////////////////////////////////////////////////////////////////
187    // alpha luminance sources..   
188    virtual void LA_to_A(unsigned int num, unsigned char* source, unsigned char* dest) const
189    {
190        for(unsigned int i=0;i<num;++i)
191        {
192            ++source;
193            *dest++ = *source++;
194        }
195    }
196
197    virtual void LA_to_LA(unsigned int num, unsigned char* source, unsigned char* dest) const
198    {
199        for(unsigned int i=0;i<num;++i)
200        {
201            *dest++ = *source++;
202            *dest++ = *source++;
203        }
204    }
205                   
206    virtual void LA_to_RGB(unsigned int num, unsigned char* source, unsigned char* dest) const
207    {
208        for(unsigned int i=0;i<num;++i)
209        {
210            *dest++ = *source;
211            *dest++ = *source;
212            *dest++ = *source;
213            source+=2;
214        }
215    }
216
217    virtual void LA_to_RGBA(unsigned int num, unsigned char* source, unsigned char* dest) const
218    {
219        for(unsigned int i=0;i<num;++i)
220        {
221            *dest++ = *source;
222            *dest++ = *source;
223            *dest++ = *source++;
224            *dest++ = *source++;
225        }
226    }
227
228    ///////////////////////////////////////////////////////////////////////////////
229    // RGB sources..   
230    virtual void RGB_to_A(unsigned int num, unsigned char* source, unsigned char* dest) const
231    {
232        for(unsigned int i=0;i<num;++i)
233        {
234            unsigned char val = *source;
235            *dest++ = val;
236            source += 3;
237        }
238    }
239
240    virtual void RGB_to_LA(unsigned int num, unsigned char* source, unsigned char* dest) const
241    {
242        for(unsigned int i=0;i<num;++i)
243        {
244            unsigned char val = *source;
245            *dest++ = val;
246            *dest++ = val;
247            source += 3;
248        }
249    }
250                   
251    virtual void RGB_to_RGB(unsigned int num, unsigned char* source, unsigned char* dest) const
252    {
253        for(unsigned int i=0;i<num;++i)
254        {
255            *dest++ = *source++;
256            *dest++ = *source++;
257            *dest++ = *source++;
258        }
259    }
260
261    virtual void RGB_to_RGBA(unsigned int num, unsigned char* source, unsigned char* dest) const
262    {
263        for(unsigned int i=0;i<num;++i)
264        {
265            unsigned char val = *source;
266            *dest++ = *source++;
267            *dest++ = *source++;
268            *dest++ = *source++;
269            *dest++ = val;
270        }
271    }
272
273    ///////////////////////////////////////////////////////////////////////////////
274    // RGBA sources..   
275    virtual void RGBA_to_A(unsigned int num, unsigned char* source, unsigned char* dest) const
276    {
277        for(unsigned int i=0;i<num;++i)
278        {
279            source += 3;
280            *dest++ = *source++;
281        }
282    }
283
284    virtual void RGBA_to_LA(unsigned int num, unsigned char* source, unsigned char* dest) const
285    {
286        for(unsigned int i=0;i<num;++i)
287        {
288            unsigned char val = *source;
289            source += 3;
290            *dest++ = val;
291            *dest++ = *source++;
292        }
293    }
294                   
295    virtual void RGBA_to_RGB(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            ++source;
303        }
304    }
305
306    virtual void RGBA_to_RGBA(unsigned int num, unsigned char* source, unsigned char* dest) const
307    {
308        for(unsigned int i=0;i<num;++i)
309        {
310            *dest++ = *source++;
311            *dest++ = *source++;
312            *dest++ = *source++;
313            *dest++ = *source++;
314        }
315    }
316};
317
318
319void clampToNearestValidPowerOfTwo(int& sizeX, int& sizeY, int& sizeZ, int s_maximumTextureSize, int t_maximumTextureSize, int r_maximumTextureSize)
320{
321    // compute nearest powers of two for each axis.
322    int s_nearestPowerOfTwo = 1;
323    while(s_nearestPowerOfTwo<sizeX && s_nearestPowerOfTwo<s_maximumTextureSize) s_nearestPowerOfTwo*=2;
324
325    int t_nearestPowerOfTwo = 1;
326    while(t_nearestPowerOfTwo<sizeY && t_nearestPowerOfTwo<t_maximumTextureSize) t_nearestPowerOfTwo*=2;
327
328    int r_nearestPowerOfTwo = 1;
329    while(r_nearestPowerOfTwo<sizeZ && r_nearestPowerOfTwo<r_maximumTextureSize) r_nearestPowerOfTwo*=2;
330
331    sizeX = s_nearestPowerOfTwo;
332    sizeY = t_nearestPowerOfTwo;
333    sizeZ = r_nearestPowerOfTwo;
334}
335
336osg::Image* createTexture3D(ImageList& imageList, ProcessRow& processRow,
337            unsigned int numComponentsDesired,
338            int s_maximumTextureSize,
339            int t_maximumTextureSize,
340            int r_maximumTextureSize,
341            bool resizeToPowerOfTwo)
342{
343    int max_s = 0;
344    int max_t = 0;
345    unsigned int max_components = 0;
346    int total_r = 0;
347    ImageList::iterator itr;
348    for(itr=imageList.begin();
349        itr!=imageList.end();
350        ++itr)
351    {
352        osg::Image* image = itr->get();
353        GLenum pixelFormat = image->getPixelFormat();
354        if (pixelFormat==GL_ALPHA ||
355            pixelFormat==GL_INTENSITY ||
356            pixelFormat==GL_LUMINANCE ||
357            pixelFormat==GL_LUMINANCE_ALPHA ||
358            pixelFormat==GL_RGB ||
359            pixelFormat==GL_RGBA)
360        {
361            max_s = osg::maximum(image->s(), max_s);
362            max_t = osg::maximum(image->t(), max_t);
363            max_components = osg::maximum(osg::Image::computeNumComponents(pixelFormat), max_components);
364            total_r += image->r();
365        }
366        else
367        {
368            osg::notify(osg::NOTICE)<<"Image "<<image->getFileName()<<" has unsuitable pixel format"<< std::hex<< pixelFormat << std::dec << std::endl;
369        }
370    }
371   
372    if (numComponentsDesired!=0) max_components = numComponentsDesired;
373   
374    GLenum desiredPixelFormat = 0;
375    switch(max_components)
376    {
377    case(1):
378        osg::notify(osg::NOTICE)<<"desiredPixelFormat = GL_LUMINANCE" << std::endl;
379        desiredPixelFormat = GL_LUMINANCE;
380        break;
381    case(2):
382        osg::notify(osg::NOTICE)<<"desiredPixelFormat = GL_LUMINANCE_ALPHA" << std::endl;
383        desiredPixelFormat = GL_LUMINANCE_ALPHA;
384        break;
385    case(3):
386        osg::notify(osg::NOTICE)<<"desiredPixelFormat = GL_RGB" << std::endl;
387        desiredPixelFormat = GL_RGB;
388        break;
389    case(4):
390        osg::notify(osg::NOTICE)<<"desiredPixelFormat = GL_RGBA" << std::endl;
391        desiredPixelFormat = GL_RGBA;
392        break;
393    }   
394    if (desiredPixelFormat==0) return 0;
395   
396    // compute nearest powers of two for each axis.
397   
398    int s_nearestPowerOfTwo = 1;
399    int t_nearestPowerOfTwo = 1;
400    int r_nearestPowerOfTwo = 1;
401
402    if (resizeToPowerOfTwo)
403    {
404        while(s_nearestPowerOfTwo<max_s && s_nearestPowerOfTwo<s_maximumTextureSize) s_nearestPowerOfTwo*=2;
405        while(t_nearestPowerOfTwo<max_t && t_nearestPowerOfTwo<t_maximumTextureSize) t_nearestPowerOfTwo*=2;
406        while(r_nearestPowerOfTwo<total_r && r_nearestPowerOfTwo<r_maximumTextureSize) r_nearestPowerOfTwo*=2;
407
408        osg::notify(osg::NOTICE)<<"max image width = "<<max_s<<"  nearest power of two = "<<s_nearestPowerOfTwo<<std::endl;
409        osg::notify(osg::NOTICE)<<"max image height = "<<max_t<<"  nearest power of two = "<<t_nearestPowerOfTwo<<std::endl;
410        osg::notify(osg::NOTICE)<<"max image depth = "<<total_r<<"  nearest power of two = "<<r_nearestPowerOfTwo<<std::endl;
411    }
412    else
413    {
414        s_nearestPowerOfTwo = max_s;
415        t_nearestPowerOfTwo = max_t;
416        r_nearestPowerOfTwo = total_r;
417    }
418   
419    // now allocate the 3d texture;
420    osg::ref_ptr<osg::Image> image_3d = new osg::Image;
421    image_3d->allocateImage(s_nearestPowerOfTwo,t_nearestPowerOfTwo,r_nearestPowerOfTwo,
422                            desiredPixelFormat,GL_UNSIGNED_BYTE);
423       
424
425    unsigned int r_offset = (total_r<r_nearestPowerOfTwo) ? r_nearestPowerOfTwo/2 - total_r/2 : 0;
426
427    int curr_dest_r = r_offset;
428
429    // copy across the values from the source images into the image_3d.
430    for(itr=imageList.begin();
431        itr!=imageList.end();
432        ++itr)
433    {
434        osg::Image* image = itr->get();
435        GLenum pixelFormat = image->getPixelFormat();
436        if (pixelFormat==GL_ALPHA ||
437            pixelFormat==GL_LUMINANCE ||
438            pixelFormat==GL_INTENSITY ||
439            pixelFormat==GL_LUMINANCE_ALPHA ||
440            pixelFormat==GL_RGB ||
441            pixelFormat==GL_RGBA)
442        {
443       
444            int num_r = osg::minimum(image->r(), (image_3d->r() - curr_dest_r));
445            int num_t = osg::minimum(image->t(), image_3d->t());
446            int num_s = osg::minimum(image->s(), image_3d->s());
447       
448            unsigned int s_offset_dest = (image->s()<s_nearestPowerOfTwo) ? s_nearestPowerOfTwo/2 - image->s()/2 : 0;
449            unsigned int t_offset_dest = (image->t()<t_nearestPowerOfTwo) ? t_nearestPowerOfTwo/2 - image->t()/2 : 0;
450
451            for(int r=0;r<num_r;++r, ++curr_dest_r)
452            {
453                for(int t=0;t<num_t;++t)
454                {
455                    unsigned char* dest = image_3d->data(s_offset_dest,t+t_offset_dest,curr_dest_r);
456                    unsigned char* source = image->data(0,t,r);
457
458                    processRow(num_s, image->getPixelFormat(), source, image_3d->getPixelFormat(), dest);
459                }
460            }
461        }
462    }
463    return image_3d.release();
464}
465
466
467osg::Image* createNormalMapTexture(osg::Image* image_3d)
468{
469    osg::notify(osg::NOTICE)<<"Computing NormalMapTexture"<<std::endl;
470
471    GLenum dataType = image_3d->getDataType();
472
473    unsigned int sourcePixelIncrement = 1;
474    unsigned int alphaOffset = 0;
475    switch(image_3d->getPixelFormat())
476    {
477    case(GL_ALPHA):
478    case(GL_LUMINANCE):
479        sourcePixelIncrement = 1;
480        alphaOffset = 0;
481        break;
482    case(GL_LUMINANCE_ALPHA):
483        sourcePixelIncrement = 2;
484        alphaOffset = 1;
485        break;
486    case(GL_RGB):
487        sourcePixelIncrement = 3;
488        alphaOffset = 0;
489        break;
490    case(GL_RGBA):
491        sourcePixelIncrement = 4;
492        alphaOffset = 3;
493        break;
494    default:
495        osg::notify(osg::NOTICE)<<"Source pixel format not support for normal map generation."<<std::endl;
496        return 0;
497    }
498
499
500    osg::ref_ptr<osg::Image> normalmap_3d = new osg::Image;
501    normalmap_3d->allocateImage(image_3d->s(),image_3d->t(),image_3d->r(),
502                            GL_RGBA,GL_UNSIGNED_BYTE);
503
504    if (osg::getCpuByteOrder()==osg::LittleEndian) alphaOffset = sourcePixelIncrement-alphaOffset-1;
505
506    for(int r=1;r<image_3d->r()-1;++r)
507    {
508        for(int t=1;t<image_3d->t()-1;++t)
509        {
510
511            if (dataType==GL_UNSIGNED_BYTE)
512            {       
513                unsigned char* ptr = image_3d->data(1,t,r)+alphaOffset;
514                unsigned char* left = image_3d->data(0,t,r)+alphaOffset;
515                unsigned char* right = image_3d->data(2,t,r)+alphaOffset;
516                unsigned char* above = image_3d->data(1,t+1,r)+alphaOffset;
517                unsigned char* below = image_3d->data(1,t-1,r)+alphaOffset;
518                unsigned char* in = image_3d->data(1,t,r+1)+alphaOffset;
519                unsigned char* out = image_3d->data(1,t,r-1)+alphaOffset;
520
521                unsigned char* destination = (unsigned char*) normalmap_3d->data(1,t,r);
522
523                for(int s=1;s<image_3d->s()-1;++s)
524                {
525
526                    osg::Vec3 grad((float)(*left)-(float)(*right),
527                                   (float)(*below)-(float)(*above),
528                                   (float)(*out) -(float)(*in));
529
530                    grad.normalize();
531
532                    if (grad.x()==0.0f && grad.y()==0.0f && grad.z()==0.0f)
533                    {
534                        grad.set(128.0f,128.0f,128.0f);
535                    }
536                    else
537                    {
538                        grad.x() = osg::clampBetween((grad.x()+1.0f)*128.0f,0.0f,255.0f);
539                        grad.y() = osg::clampBetween((grad.y()+1.0f)*128.0f,0.0f,255.0f);
540                        grad.z() = osg::clampBetween((grad.z()+1.0f)*128.0f,0.0f,255.0f);
541                    }
542
543                    *(destination++) = (unsigned char)(grad.x()); // scale and bias X.
544                    *(destination++) = (unsigned char)(grad.y()); // scale and bias Y.
545                    *(destination++) = (unsigned char)(grad.z()); // scale and bias Z.
546
547                    *destination++ = *ptr;
548
549                    ptr += sourcePixelIncrement;
550                    left += sourcePixelIncrement;
551                    right += sourcePixelIncrement;
552                    above += sourcePixelIncrement;
553                    below += sourcePixelIncrement;
554                    in += sourcePixelIncrement;
555                    out += sourcePixelIncrement;
556                }
557            }
558            else if (dataType==GL_SHORT)
559            {
560                short* ptr = (short*)(image_3d->data(1,t,r)+alphaOffset);
561                short* left = (short*)(image_3d->data(0,t,r)+alphaOffset);
562                short* right = (short*)(image_3d->data(2,t,r)+alphaOffset);
563                short* above = (short*)(image_3d->data(1,t+1,r)+alphaOffset);
564                short* below = (short*)(image_3d->data(1,t-1,r)+alphaOffset);
565                short* in = (short*)(image_3d->data(1,t,r+1)+alphaOffset);
566                short* out = (short*)(image_3d->data(1,t,r-1)+alphaOffset);
567
568                unsigned char* destination = (unsigned char*) normalmap_3d->data(1,t,r);
569
570                for(int s=1;s<image_3d->s()-1;++s)
571                {
572
573                    osg::Vec3 grad((float)(*left)-(float)(*right),
574                                   (float)(*below)-(float)(*above),
575                                   (float)(*out) -(float)(*in));
576
577                    grad.normalize();
578
579                    //osg::notify(osg::NOTICE)<<"normal "<<grad<<std::endl;
580
581                    if (grad.x()==0.0f && grad.y()==0.0f && grad.z()==0.0f)
582                    {
583                        grad.set(128.0f,128.0f,128.0f);
584                    }
585                    else
586                    {
587                        grad.x() = osg::clampBetween((grad.x()+1.0f)*128.0f,0.0f,255.0f);
588                        grad.y() = osg::clampBetween((grad.y()+1.0f)*128.0f,0.0f,255.0f);
589                        grad.z() = osg::clampBetween((grad.z()+1.0f)*128.0f,0.0f,255.0f);
590                    }
591                   
592
593                    *(destination++) = (unsigned char)(grad.x()); // scale and bias X.
594                    *(destination++) = (unsigned char)(grad.y()); // scale and bias Y.
595                    *(destination++) = (unsigned char)(grad.z()); // scale and bias Z.
596
597                    *destination++ = *ptr/128;
598
599                    ptr += sourcePixelIncrement;
600                    left += sourcePixelIncrement;
601                    right += sourcePixelIncrement;
602                    above += sourcePixelIncrement;
603                    below += sourcePixelIncrement;
604                    in += sourcePixelIncrement;
605                    out += sourcePixelIncrement;
606                }
607            }
608            else if (dataType==GL_UNSIGNED_SHORT)
609            {
610                unsigned short* ptr = (unsigned short*)(image_3d->data(1,t,r)+alphaOffset);
611                unsigned short* left = (unsigned short*)(image_3d->data(0,t,r)+alphaOffset);
612                unsigned short* right = (unsigned short*)(image_3d->data(2,t,r)+alphaOffset);
613                unsigned short* above = (unsigned short*)(image_3d->data(1,t+1,r)+alphaOffset);
614                unsigned short* below = (unsigned short*)(image_3d->data(1,t-1,r)+alphaOffset);
615                unsigned short* in = (unsigned short*)(image_3d->data(1,t,r+1)+alphaOffset);
616                unsigned short* out = (unsigned short*)(image_3d->data(1,t,r-1)+alphaOffset);
617
618                unsigned char* destination = (unsigned char*) normalmap_3d->data(1,t,r);
619
620                for(int s=1;s<image_3d->s()-1;++s)
621                {
622
623                    osg::Vec3 grad((float)(*left)-(float)(*right),
624                                   (float)(*below)-(float)(*above),
625                                   (float)(*out) -(float)(*in));
626
627                    grad.normalize();
628
629                    if (grad.x()==0.0f && grad.y()==0.0f && grad.z()==0.0f)
630                    {
631                        grad.set(128.0f,128.0f,128.0f);
632                    }
633                    else
634                    {
635                        grad.x() = osg::clampBetween((grad.x()+1.0f)*128.0f,0.0f,255.0f);
636                        grad.y() = osg::clampBetween((grad.y()+1.0f)*128.0f,0.0f,255.0f);
637                        grad.z() = osg::clampBetween((grad.z()+1.0f)*128.0f,0.0f,255.0f);
638                    }
639
640                    *(destination++) = (unsigned char)(grad.x()); // scale and bias X.
641                    *(destination++) = (unsigned char)(grad.y()); // scale and bias Y.
642                    *(destination++) = (unsigned char)(grad.z()); // scale and bias Z.
643
644                    *destination++ = *ptr/256;
645
646                    ptr += sourcePixelIncrement;
647                    left += sourcePixelIncrement;
648                    right += sourcePixelIncrement;
649                    above += sourcePixelIncrement;
650                    below += sourcePixelIncrement;
651                    in += sourcePixelIncrement;
652                    out += sourcePixelIncrement;
653                }
654            }
655        }
656    }
657   
658   
659    osg::notify(osg::NOTICE)<<"Created NormalMapTexture"<<std::endl;
660   
661    return normalmap_3d.release();
662}
663
664
665
666osg::Node* createCube(float size,float alpha, unsigned int numSlices, float sliceEnd=1.0f)
667{
668
669    // set up the Geometry.
670    osg::Geometry* geom = new osg::Geometry;
671
672    float halfSize = size*0.5f;
673    float y = halfSize;
674    float dy =-size/(float)(numSlices-1)*sliceEnd;
675
676    //y = -halfSize;
677    //dy *= 0.5;
678
679    osg::Vec3Array* coords = new osg::Vec3Array(4*numSlices);
680    geom->setVertexArray(coords);
681    for(unsigned int i=0;i<numSlices;++i, y+=dy)
682    {
683        (*coords)[i*4+0].set(-halfSize,y,halfSize);
684        (*coords)[i*4+1].set(-halfSize,y,-halfSize);
685        (*coords)[i*4+2].set(halfSize,y,-halfSize);
686        (*coords)[i*4+3].set(halfSize,y,halfSize);
687    }
688   
689    osg::Vec3Array* normals = new osg::Vec3Array(1);
690    (*normals)[0].set(0.0f,-1.0f,0.0f);
691    geom->setNormalArray(normals);
692    geom->setNormalBinding(osg::Geometry::BIND_OVERALL);
693
694    osg::Vec4Array* colors = new osg::Vec4Array(1);
695    (*colors)[0].set(1.0f,1.0f,1.0f,alpha);
696    geom->setColorArray(colors);
697    geom->setColorBinding(osg::Geometry::BIND_OVERALL);
698
699    geom->addPrimitiveSet(new osg::DrawArrays(osg::PrimitiveSet::QUADS,0,coords->size()));
700
701    osg::Billboard* billboard = new osg::Billboard;
702    billboard->setMode(osg::Billboard::POINT_ROT_WORLD);
703    billboard->addDrawable(geom);
704    billboard->setPosition(0,osg::Vec3(0.0f,0.0f,0.0f));
705   
706    return billboard;
707}
708
709class FollowMouseCallback : public osgGA::GUIEventHandler, public osg::StateSet::Callback
710{
711    public:
712   
713        FollowMouseCallback(bool shader = false):
714            _shader(shader)
715        {
716            _updateTransparency = false;
717            _updateAlphaCutOff = false;
718            _updateSampleDensity = false;
719        }
720
721        FollowMouseCallback(const FollowMouseCallback&,const osg::CopyOp&) {}
722
723        META_Object(osg,FollowMouseCallback);
724
725        virtual void operator() (osg::StateSet* stateset, osg::NodeVisitor* nv)
726        {
727            if (nv->getVisitorType()==osg::NodeVisitor::EVENT_VISITOR)
728            {
729                osgGA::EventVisitor* ev = dynamic_cast<osgGA::EventVisitor*>(nv);
730                if (ev)
731                {
732                    osgGA::GUIActionAdapter* aa = ev->getActionAdapter();
733                    osgGA::EventQueue::Events& events = ev->getEvents();
734                    for(osgGA::EventQueue::Events::iterator itr=events.begin();
735                        itr!=events.end();
736                        ++itr)
737                    {
738                        handle(*(*itr), *aa, stateset, ev);
739                    }
740                }
741            }
742        }
743       
744        virtual bool handle(const osgGA::GUIEventAdapter& ea,osgGA::GUIActionAdapter&, osg::Object* object, osg::NodeVisitor*)
745        {
746            osg::StateSet* stateset = dynamic_cast<osg::StateSet*>(object);
747            if (!stateset) return false;
748           
749            switch(ea.getEventType())
750            {
751                case(osgGA::GUIEventAdapter::MOVE):
752                case(osgGA::GUIEventAdapter::DRAG):
753                {
754                    float v = (ea.getY()-ea.getYmin())/(ea.getYmax()-ea.getYmin());
755                    if (_shader)
756                    {
757                        osg::Uniform* uniform = 0;
758                        if (_updateTransparency && (uniform = stateset->getUniform("transparency"))) uniform->set(v);
759                        if (_updateAlphaCutOff && (uniform = stateset->getUniform("alphaCutOff"))) uniform->set(v);
760                        if (_updateSampleDensity && (uniform = stateset->getUniform("sampleDensity")))
761                        {
762                            float value = powf(v,5);
763                            osg::notify(osg::INFO)<<"sampleDensity = "<<value<<std::endl;
764                            uniform->set(value);
765                        }
766                    }
767                    else
768                    {                   
769                        if (_updateAlphaCutOff)
770                        {
771                            osg::AlphaFunc* alphaFunc = dynamic_cast<osg::AlphaFunc*>(stateset->getAttribute(osg::StateAttribute::ALPHAFUNC));
772                            if (alphaFunc)
773                            {
774                                alphaFunc->setReferenceValue(v);
775                            }
776                        }
777                       
778                        if (_updateTransparency)
779                        {
780                            osg::Material* material = dynamic_cast<osg::Material*>(stateset->getAttribute(osg::StateAttribute::MATERIAL));
781                            if (material)
782                            {
783                                material->setAlpha(osg::Material::FRONT_AND_BACK,v);
784                            }
785                        }
786                    }
787
788                    break;
789                }
790                case(osgGA::GUIEventAdapter::KEYDOWN):
791                {
792                    if (ea.getKey()=='t') _updateTransparency = true;
793                    if (ea.getKey()=='a') _updateAlphaCutOff = true;
794                    if (ea.getKey()=='d') _updateSampleDensity = true;
795                    break;
796                }
797                case(osgGA::GUIEventAdapter::KEYUP):
798                {
799                    if (ea.getKey()=='t') _updateTransparency = false;
800                    if (ea.getKey()=='a') _updateAlphaCutOff = false;
801                    if (ea.getKey()=='d') _updateSampleDensity = false;
802                    break;
803                }
804                default:
805                    break;
806            }
807            return false;
808        }
809       
810        bool _shader;
811        bool _updateTransparency;
812        bool _updateAlphaCutOff;
813        bool _updateSampleDensity;
814
815};
816
817osg::Node* createShaderModel(ShadingModel shadingModel,
818                       osg::ref_ptr<osg::Image>& image_3d,
819                       osg::Image* normalmap_3d,
820                       osg::TransferFunction1D* tf,
821                       osg::Texture::InternalFormatMode internalFormatMode,
822                       float xSize, float ySize, float zSize,
823                       float /*xMultiplier*/, float /*yMultiplier*/, float /*zMultiplier*/,
824                       unsigned int /*numSlices*/=500, float /*sliceEnd*/=1.0f, float alphaFuncValue=0.02f)
825{
826    osg::Texture::FilterMode minFilter = osg::Texture::LINEAR;
827    osg::Texture::FilterMode magFilter = osg::Texture::LINEAR;
828
829    osg::Group* root = new osg::Group;
830   
831    osg::Geode* geode = new osg::Geode;
832    root->addChild(geode);
833   
834    osg::StateSet* stateset = geode->getOrCreateStateSet();
835   
836    stateset->setEventCallback(new FollowMouseCallback(true));
837   
838    stateset->setMode(GL_ALPHA_TEST,osg::StateAttribute::ON);
839
840   
841    osg::Program* program = new osg::Program;
842    stateset->setAttribute(program);
843
844    // get shaders from source
845   
846    osg::Shader* vertexShader = osgDB::readShaderFile(osg::Shader::VERTEX, "volume.vert");
847    if (vertexShader)
848    {
849        program->addShader(vertexShader);
850    }
851    else
852    {
853        #include "volume_vert.cpp"
854        program->addShader(new osg::Shader(osg::Shader::VERTEX, volume_vert));
855    }
856
857    if (!(normalmap_3d && tf))
858    {
859        // set up the 3d texture itself,
860        // note, well set the filtering up so that mip mapping is disabled,
861        // gluBuild3DMipsmaps doesn't do a very good job of handled the
862        // imbalanced dimensions of the 256x256x4 texture.
863        osg::Texture3D* texture3D = new osg::Texture3D;
864        texture3D->setResizeNonPowerOfTwoHint(false);
865        texture3D->setFilter(osg::Texture3D::MIN_FILTER,minFilter);
866        texture3D->setFilter(osg::Texture3D::MAG_FILTER, magFilter);
867        texture3D->setWrap(osg::Texture3D::WRAP_R,osg::Texture3D::CLAMP_TO_EDGE);
868        texture3D->setWrap(osg::Texture3D::WRAP_S,osg::Texture3D::CLAMP_TO_EDGE);
869        texture3D->setWrap(osg::Texture3D::WRAP_T,osg::Texture3D::CLAMP_TO_EDGE);
870        if (image_3d->getPixelFormat()==GL_ALPHA ||
871            image_3d->getPixelFormat()==GL_LUMINANCE)
872        {
873            texture3D->setInternalFormatMode(osg::Texture3D::USE_USER_DEFINED_FORMAT);
874            texture3D->setInternalFormat(GL_INTENSITY);
875        }
876        else
877        {
878            texture3D->setInternalFormatMode(internalFormatMode);
879        }
880        texture3D->setImage(image_3d.get());
881
882        stateset->setTextureAttributeAndModes(0,texture3D,osg::StateAttribute::ON);
883
884        osg::Uniform* baseTextureSampler = new osg::Uniform("baseTexture",0);
885        stateset->addUniform(baseTextureSampler);
886    }
887   
888
889    if (shadingModel==MaximumIntensityProjection)
890    {
891        if (tf)
892        {
893            osg::Texture1D* texture1D = new osg::Texture1D;
894            texture1D->setImage(tf->getImage());   
895            stateset->setTextureAttributeAndModes(1,texture1D,osg::StateAttribute::ON);
896
897            osg::Shader* fragmentShader = osgDB::readShaderFile(osg::Shader::FRAGMENT, "volume_tf_mip.frag");
898            if (fragmentShader)
899            {
900                program->addShader(fragmentShader);
901            }
902            else
903            {
904                #include "volume_tf_mip_frag.cpp"
905                program->addShader(new osg::Shader(osg::Shader::FRAGMENT, volume_tf_mip_frag));
906            }
907
908            osg::Uniform* tfTextureSampler = new osg::Uniform("tfTexture",1);
909            stateset->addUniform(tfTextureSampler);
910
911        }
912        else
913        {   
914            osg::Shader* fragmentShader = osgDB::readShaderFile(osg::Shader::FRAGMENT, "volume_mip.frag");
915            if (fragmentShader)
916            {
917                program->addShader(fragmentShader);
918            }
919            else
920            {
921                #include "volume_mip_frag.cpp"
922                program->addShader(new osg::Shader(osg::Shader::FRAGMENT, volume_mip_frag));
923            }
924        }
925    }
926    else if (shadingModel==Isosurface)
927    {
928
929        if (tf)
930        {
931            osg::Texture1D* texture1D = new osg::Texture1D;
932            texture1D->setImage(tf->getImage());   
933            texture1D->setResizeNonPowerOfTwoHint(false);
934            texture1D->setFilter(osg::Texture::MIN_FILTER, osg::Texture::LINEAR);
935            texture1D->setFilter(osg::Texture::MAG_FILTER, osg::Texture::LINEAR);
936            texture1D->setWrap(osg::Texture::WRAP_R,osg::Texture::CLAMP_TO_EDGE);
937            stateset->setTextureAttributeAndModes(1,texture1D,osg::StateAttribute::ON);
938
939            osg::Uniform* tfTextureSampler = new osg::Uniform("tfTexture",1);
940            stateset->addUniform(tfTextureSampler);
941
942            osg::Shader* fragmentShader = osgDB::readShaderFile(osg::Shader::FRAGMENT, "volume_tf_iso.frag");
943            if (fragmentShader)
944            {
945                program->addShader(fragmentShader);
946            }
947            else
948            {
949                #include "volume_tf_iso_frag.cpp"
950                program->addShader(new osg::Shader(osg::Shader::FRAGMENT, volume_tf_iso_frag));
951            }
952        }
953        else
954        {   
955            osg::Shader* fragmentShader = osgDB::readShaderFile(osg::Shader::FRAGMENT, "volume_iso.frag");
956            if (fragmentShader)
957            {
958                program->addShader(fragmentShader);
959            }
960            else
961            {
962                #include "volume_iso_frag.cpp"
963                program->addShader(new osg::Shader(osg::Shader::FRAGMENT, volume_iso_frag));
964            }
965        }
966    }
967    else if (normalmap_3d)
968    {
969        osg::notify(osg::NOTICE)<<"Setting up normalmapping shader"<<std::endl;
970
971        osg::Uniform* normalMapSampler = new osg::Uniform("normalMap",1);
972        stateset->addUniform(normalMapSampler);
973
974        osg::Texture3D* normalMap = new osg::Texture3D;
975        normalMap->setImage(normalmap_3d);   
976        normalMap->setResizeNonPowerOfTwoHint(false);
977        normalMap->setInternalFormatMode(internalFormatMode);
978        normalMap->setFilter(osg::Texture3D::MIN_FILTER, osg::Texture::LINEAR);
979        normalMap->setFilter(osg::Texture3D::MAG_FILTER, osg::Texture::LINEAR);
980        normalMap->setWrap(osg::Texture3D::WRAP_R,osg::Texture3D::CLAMP_TO_EDGE);
981        normalMap->setWrap(osg::Texture3D::WRAP_S,osg::Texture3D::CLAMP_TO_EDGE);
982        normalMap->setWrap(osg::Texture3D::WRAP_T,osg::Texture3D::CLAMP_TO_EDGE);
983
984        stateset->setTextureAttributeAndModes(1,normalMap,osg::StateAttribute::ON);
985
986        if (tf)
987        {
988            osg::Texture1D* texture1D = new osg::Texture1D;
989            texture1D->setImage(tf->getImage());   
990            texture1D->setResizeNonPowerOfTwoHint(false);
991            texture1D->setFilter(osg::Texture::MIN_FILTER, osg::Texture::LINEAR);
992            texture1D->setFilter(osg::Texture::MAG_FILTER, osg::Texture::LINEAR);
993            texture1D->setWrap(osg::Texture::WRAP_R,osg::Texture::CLAMP_TO_EDGE);
994            stateset->setTextureAttributeAndModes(0,texture1D,osg::StateAttribute::ON);
995
996            osg::Shader* fragmentShader = osgDB::readShaderFile(osg::Shader::FRAGMENT, "volume-tf-n.frag");
997            if (fragmentShader)
998            {
999                program->addShader(fragmentShader);
1000            }
1001            else
1002            {
1003                #include "volume_tf_n_frag.cpp"
1004                program->addShader(new osg::Shader(osg::Shader::FRAGMENT, volume_tf_n_frag));
1005            }
1006
1007            osg::Uniform* tfTextureSampler = new osg::Uniform("tfTexture",0);
1008            stateset->addUniform(tfTextureSampler);
1009        }
1010        else
1011        {
1012            osg::Shader* fragmentShader = osgDB::readShaderFile(osg::Shader::FRAGMENT, "volume-n.frag");
1013            if (fragmentShader)
1014            {
1015                program->addShader(fragmentShader);
1016            }
1017            else
1018            {
1019                #include "volume_n_frag.cpp"
1020                program->addShader(new osg::Shader(osg::Shader::FRAGMENT, volume_n_frag));
1021            }
1022        }
1023    }
1024    else if (tf)
1025    {
1026        osg::Texture1D* texture1D = new osg::Texture1D;
1027        texture1D->setImage(tf->getImage());   
1028        texture1D->setResizeNonPowerOfTwoHint(false);
1029        texture1D->setFilter(osg::Texture::MIN_FILTER, osg::Texture::LINEAR);
1030        texture1D->setFilter(osg::Texture::MAG_FILTER, osg::Texture::LINEAR);
1031        texture1D->setWrap(osg::Texture::WRAP_R,osg::Texture::CLAMP_TO_EDGE);
1032        stateset->setTextureAttributeAndModes(1,texture1D,osg::StateAttribute::ON);
1033
1034        osg::Uniform* tfTextureSampler = new osg::Uniform("tfTexture",1);
1035        stateset->addUniform(tfTextureSampler);
1036
1037        osg::Shader* fragmentShader = osgDB::readShaderFile(osg::Shader::FRAGMENT, "volume-tf.frag");
1038        if (fragmentShader)
1039        {
1040            program->addShader(fragmentShader);
1041        }
1042        else
1043        {
1044            #include "volume_tf_frag.cpp"
1045            program->addShader(new osg::Shader(osg::Shader::FRAGMENT, volume_tf_frag));
1046        }
1047
1048    }
1049    else
1050    {   
1051
1052        osg::Shader* fragmentShader = osgDB::readShaderFile(osg::Shader::FRAGMENT, "volume.frag");
1053        if (fragmentShader)
1054        {
1055            program->addShader(fragmentShader);
1056        }
1057        else
1058        {
1059            #include "volume_frag.cpp"
1060            program->addShader(new osg::Shader(osg::Shader::FRAGMENT, volume_frag));
1061        }
1062    }
1063
1064    osg::Uniform* sampleDensity = new osg::Uniform("sampleDensity", 0.005f);
1065    stateset->addUniform(sampleDensity);
1066
1067    osg::Uniform* transpancy = new osg::Uniform("transparency",0.5f);
1068    stateset->addUniform(transpancy);
1069
1070    osg::Uniform* alphaCutOff = new osg::Uniform("alphaCutOff",alphaFuncValue);
1071    stateset->addUniform(alphaCutOff);
1072
1073    stateset->setMode(GL_CULL_FACE, osg::StateAttribute::ON);
1074
1075    osg::TexGen* texgen = new osg::TexGen;
1076    texgen->setMode(osg::TexGen::OBJECT_LINEAR);
1077    texgen->setPlane(osg::TexGen::S, osg::Plane(1.0f/xSize,0.0f,0.0f,0.0f));
1078    texgen->setPlane(osg::TexGen::T, osg::Plane(0.0f,1.0f/ySize,0.0f,0.0f));
1079    texgen->setPlane(osg::TexGen::R, osg::Plane(0.0f,0.0f,1.0f/zSize,0.0f));
1080    texgen->setPlane(osg::TexGen::Q, osg::Plane(0.0f,0.0f,0.0f,1.0f));
1081   
1082    stateset->setTextureAttributeAndModes(0, texgen, osg::StateAttribute::ON);
1083
1084    {
1085        osg::Geometry* geom = new osg::Geometry;
1086
1087        osg::Vec3Array* coords = new osg::Vec3Array(8);
1088        (*coords)[0].set(0,0,0);
1089        (*coords)[1].set(xSize,0,0);
1090        (*coords)[2].set(xSize,ySize,0);
1091        (*coords)[3].set(0,ySize,0);
1092        (*coords)[4].set(0,0,zSize);
1093        (*coords)[5].set(xSize,0,zSize);
1094        (*coords)[6].set(ySize,ySize,zSize);
1095        (*coords)[7].set(0,ySize,zSize);
1096        geom->setVertexArray(coords);
1097
1098        osg::Vec4Array* colours = new osg::Vec4Array(1);
1099        (*colours)[0].set(1.0f,1.0f,1.0,1.0f);
1100        geom->setColorArray(colours);
1101        geom->setColorBinding(osg::Geometry::BIND_OVERALL);
1102
1103        osg::DrawElementsUShort* drawElements = new osg::DrawElementsUShort(GL_QUADS);
1104        // bottom
1105        drawElements->push_back(0);
1106        drawElements->push_back(1);
1107        drawElements->push_back(2);
1108        drawElements->push_back(3);
1109       
1110        // bottom
1111        drawElements->push_back(3);
1112        drawElements->push_back(2);
1113        drawElements->push_back(6);
1114        drawElements->push_back(7);
1115
1116        // left
1117        drawElements->push_back(0);
1118        drawElements->push_back(3);
1119        drawElements->push_back(7);
1120        drawElements->push_back(4);
1121
1122        // right
1123        drawElements->push_back(5);
1124        drawElements->push_back(6);
1125        drawElements->push_back(2);
1126        drawElements->push_back(1);
1127
1128        // front
1129        drawElements->push_back(1);
1130        drawElements->push_back(0);
1131        drawElements->push_back(4);
1132        drawElements->push_back(5);
1133
1134        // top
1135        drawElements->push_back(7);
1136        drawElements->push_back(6);
1137        drawElements->push_back(5);
1138        drawElements->push_back(4);
1139
1140        geom->addPrimitiveSet(drawElements);
1141
1142        geode->addDrawable(geom);
1143
1144    }
1145    return root;
1146}
1147
1148osg::Node* createModel(ShadingModel shadeModel,
1149                       osg::ref_ptr<osg::Image>& image_3d,
1150                       osg::ref_ptr<osg::Image>& normalmap_3d,
1151                       osg::Texture::InternalFormatMode internalFormatMode,
1152                       float xSize, float ySize, float zSize,
1153                       float xMultiplier, float yMultiplier, float zMultiplier,
1154                       unsigned int numSlices=500, float sliceEnd=1.0f, float alphaFuncValue=0.02f, bool maximumIntensityProjection = false)
1155{
1156    bool two_pass = normalmap_3d.valid() && (image_3d->getPixelFormat()==GL_RGB || image_3d->getPixelFormat()==GL_RGBA);
1157
1158    osg::BoundingBox bb(-xSize*0.5f,-ySize*0.5f,-zSize*0.5f,xSize*0.5f,ySize*0.5f,zSize*0.5f);
1159
1160
1161    osg::Texture::FilterMode minFilter = osg::Texture::NEAREST;
1162    osg::Texture::FilterMode magFilter = osg::Texture::NEAREST;
1163
1164    float maxAxis = xSize;
1165    if (ySize > maxAxis) maxAxis = ySize;
1166    if (zSize > maxAxis) maxAxis = zSize;
1167
1168    osg::Group* group = new osg::Group;
1169   
1170    osg::TexGenNode* texgenNode_0 = new osg::TexGenNode;
1171    texgenNode_0->setTextureUnit(0);
1172    texgenNode_0->getTexGen()->setMode(osg::TexGen::EYE_LINEAR);
1173    texgenNode_0->getTexGen()->setPlane(osg::TexGen::S, osg::Plane(xMultiplier/xSize,0.0f,0.0f,0.5f));
1174    texgenNode_0->getTexGen()->setPlane(osg::TexGen::T, osg::Plane(0.0f,yMultiplier/ySize,0.0f,0.5f));
1175    texgenNode_0->getTexGen()->setPlane(osg::TexGen::R, osg::Plane(0.0f,0.0f,zMultiplier/zSize,0.5f));
1176   
1177    if (two_pass)
1178    {
1179        osg::TexGenNode* texgenNode_1 = new osg::TexGenNode;
1180        texgenNode_1->setTextureUnit(1);
1181        texgenNode_1->getTexGen()->setMode(osg::TexGen::EYE_LINEAR);
1182        texgenNode_1->getTexGen()->setPlane(osg::TexGen::S, texgenNode_0->getTexGen()->getPlane(osg::TexGen::S));
1183        texgenNode_1->getTexGen()->setPlane(osg::TexGen::T, texgenNode_0->getTexGen()->getPlane(osg::TexGen::T));
1184        texgenNode_1->getTexGen()->setPlane(osg::TexGen::R, texgenNode_0->getTexGen()->getPlane(osg::TexGen::R));
1185
1186        texgenNode_1->addChild(texgenNode_0);
1187
1188        group->addChild(texgenNode_1);
1189    }
1190    else
1191    { 
1192        group->addChild(texgenNode_0);
1193    }
1194
1195    float cubeSize = sqrtf(xSize*xSize+ySize*ySize+zSize*zSize);
1196
1197    osg::ClipNode* clipnode = new osg::ClipNode;
1198    clipnode->addChild(createCube(cubeSize,1.0f, numSlices,sliceEnd));
1199    clipnode->createClipBox(bb);
1200
1201    {
1202        // set up the Geometry to enclose the clip volume to prevent near/far clipping from affecting billboard
1203        osg::Geometry* geom = new osg::Geometry;
1204
1205        osg::Vec3Array* coords = new osg::Vec3Array();
1206        coords->push_back(bb.corner(0));
1207        coords->push_back(bb.corner(1));
1208        coords->push_back(bb.corner(2));
1209        coords->push_back(bb.corner(3));
1210        coords->push_back(bb.corner(4));
1211        coords->push_back(bb.corner(5));
1212        coords->push_back(bb.corner(6));
1213        coords->push_back(bb.corner(7));
1214
1215        geom->setVertexArray(coords);
1216
1217        osg::Vec4Array* colors = new osg::Vec4Array(1);
1218        (*colors)[0].set(1.0f,1.0f,1.0f,1.0f);
1219        geom->setColorArray(colors);
1220        geom->setColorBinding(osg::Geometry::BIND_OVERALL);
1221
1222        geom->addPrimitiveSet(new osg::DrawArrays(osg::PrimitiveSet::POINTS,0,coords->size()));
1223
1224        osg::Geode* geode = new osg::Geode;
1225        geode->addDrawable(geom);
1226       
1227        clipnode->addChild(geode);
1228       
1229    }
1230
1231    texgenNode_0->addChild(clipnode);
1232
1233    osg::StateSet* stateset = texgenNode_0->getOrCreateStateSet();
1234
1235    stateset->setEventCallback(new FollowMouseCallback(false));
1236 
1237    stateset->setMode(GL_LIGHTING,osg::StateAttribute::ON);
1238    stateset->setMode(GL_BLEND,osg::StateAttribute::ON);
1239    stateset->setAttributeAndModes(new osg::AlphaFunc(osg::AlphaFunc::GREATER,alphaFuncValue), osg::StateAttribute::ON);
1240   
1241    osg::Material* material = new osg::Material;
1242    material->setDiffuse(osg::Material::FRONT_AND_BACK,osg::Vec4(1.0f,1.0f,1.0f,1.0f));
1243    stateset->setAttributeAndModes(material);
1244   
1245    if (shadeModel==MaximumIntensityProjection)
1246    {
1247        stateset->setAttribute(new osg::BlendFunc(osg::BlendFunc::ONE, osg::BlendFunc::ONE));
1248        stateset->setAttribute(new osg::BlendEquation(osg::BlendEquation::RGBA_MAX));
1249    }
1250   
1251    osg::Vec3 lightDirection(1.0f,-1.0f,1.0f);
1252    lightDirection.normalize();
1253
1254    if (normalmap_3d.valid())
1255    {
1256        if (two_pass)
1257        {
1258
1259            // set up normal texture
1260            osg::Texture3D* bump_texture3D = new osg::Texture3D;
1261            bump_texture3D->setFilter(osg::Texture3D::MIN_FILTER,minFilter);
1262            bump_texture3D->setFilter(osg::Texture3D::MAG_FILTER, magFilter);
1263            bump_texture3D->setWrap(osg::Texture3D::WRAP_R,osg::Texture3D::CLAMP_TO_EDGE);
1264            bump_texture3D->setWrap(osg::Texture3D::WRAP_S,osg::Texture3D::CLAMP_TO_EDGE);
1265            bump_texture3D->setWrap(osg::Texture3D::WRAP_T,osg::Texture3D::CLAMP_TO_EDGE);
1266            bump_texture3D->setImage(normalmap_3d.get());
1267
1268            bump_texture3D->setInternalFormatMode(internalFormatMode);
1269
1270            stateset->setTextureAttributeAndModes(0,bump_texture3D,osg::StateAttribute::ON);
1271
1272            osg::TexEnvCombine* tec = new osg::TexEnvCombine;
1273            tec->setConstantColorAsLightDirection(lightDirection);
1274
1275            tec->setCombine_RGB(osg::TexEnvCombine::DOT3_RGB);
1276            tec->setSource0_RGB(osg::TexEnvCombine::CONSTANT);
1277            tec->setOperand0_RGB(osg::TexEnvCombine::SRC_COLOR);
1278            tec->setSource1_RGB(osg::TexEnvCombine::TEXTURE);
1279
1280            tec->setOperand1_RGB(osg::TexEnvCombine::SRC_COLOR);
1281
1282            tec->setCombine_Alpha(osg::TexEnvCombine::REPLACE);
1283            tec->setSource0_Alpha(osg::TexEnvCombine::PRIMARY_COLOR);
1284            tec->setOperand0_Alpha(osg::TexEnvCombine::SRC_ALPHA);
1285            tec->setSource1_Alpha(osg::TexEnvCombine::TEXTURE);
1286            tec->setOperand1_Alpha(osg::TexEnvCombine::SRC_ALPHA);
1287
1288            stateset->setTextureAttributeAndModes(0, tec, osg::StateAttribute::OVERRIDE|osg::StateAttribute::ON);
1289
1290            stateset->setTextureMode(0,GL_TEXTURE_GEN_S,osg::StateAttribute::ON);
1291            stateset->setTextureMode(0,GL_TEXTURE_GEN_T,osg::StateAttribute::ON);
1292            stateset->setTextureMode(0,GL_TEXTURE_GEN_R,osg::StateAttribute::ON);
1293
1294
1295            // set up color texture
1296            osg::Texture3D* texture3D = new osg::Texture3D;
1297            texture3D->setResizeNonPowerOfTwoHint(false);
1298            texture3D->setFilter(osg::Texture3D::MIN_FILTER,minFilter);
1299            texture3D->setFilter(osg::Texture3D::MAG_FILTER, magFilter);
1300            texture3D->setWrap(osg::Texture3D::WRAP_R,osg::Texture3D::CLAMP_TO_EDGE);
1301            texture3D->setWrap(osg::Texture3D::WRAP_S,osg::Texture3D::CLAMP_TO_EDGE);
1302            texture3D->setWrap(osg::Texture3D::WRAP_T,osg::Texture3D::CLAMP_TO_EDGE);
1303            if (image_3d->getPixelFormat()==GL_ALPHA ||
1304                image_3d->getPixelFormat()==GL_LUMINANCE)
1305            {
1306                texture3D->setInternalFormatMode(osg::Texture3D::USE_USER_DEFINED_FORMAT);
1307                texture3D->setInternalFormat(GL_INTENSITY);
1308            }
1309            else
1310            {
1311                texture3D->setInternalFormatMode(internalFormatMode);
1312            }
1313            texture3D->setImage(image_3d.get());
1314
1315            stateset->setTextureAttributeAndModes(1,texture3D,osg::StateAttribute::ON);
1316
1317            stateset->setTextureMode(1,GL_TEXTURE_GEN_S,osg::StateAttribute::ON);
1318            stateset->setTextureMode(1,GL_TEXTURE_GEN_T,osg::StateAttribute::ON);
1319            stateset->setTextureMode(1,GL_TEXTURE_GEN_R,osg::StateAttribute::ON);
1320
1321            stateset->setTextureAttributeAndModes(1,new osg::TexEnv(),osg::StateAttribute::ON);
1322
1323        }
1324        else
1325        {
1326            osg::Texture3D* bump_texture3D = new osg::Texture3D;
1327            bump_texture3D->setResizeNonPowerOfTwoHint(false);
1328            bump_texture3D->setFilter(osg::Texture3D::MIN_FILTER,minFilter);
1329            bump_texture3D->setFilter(osg::Texture3D::MAG_FILTER, magFilter);
1330            bump_texture3D->setWrap(osg::Texture3D::WRAP_R,osg::Texture3D::CLAMP_TO_EDGE);
1331            bump_texture3D->setWrap(osg::Texture3D::WRAP_S,osg::Texture3D::CLAMP_TO_EDGE);
1332            bump_texture3D->setWrap(osg::Texture3D::WRAP_T,osg::Texture3D::CLAMP_TO_EDGE);
1333            bump_texture3D->setImage(normalmap_3d.get());
1334
1335            bump_texture3D->setInternalFormatMode(internalFormatMode);
1336
1337            stateset->setTextureAttributeAndModes(0,bump_texture3D,osg::StateAttribute::ON);
1338
1339            osg::TexEnvCombine* tec = new osg::TexEnvCombine;
1340            tec->setConstantColorAsLightDirection(lightDirection);
1341
1342            tec->setCombine_RGB(osg::TexEnvCombine::DOT3_RGB);
1343            tec->setSource0_RGB(osg::TexEnvCombine::CONSTANT);
1344            tec->setOperand0_RGB(osg::TexEnvCombine::SRC_COLOR);
1345            tec->setSource1_RGB(osg::TexEnvCombine::TEXTURE);
1346            tec->setOperand1_RGB(osg::TexEnvCombine::SRC_COLOR);
1347
1348            tec->setCombine_Alpha(osg::TexEnvCombine::MODULATE);
1349            tec->setSource0_Alpha(osg::TexEnvCombine::PRIMARY_COLOR);
1350            tec->setOperand0_Alpha(osg::TexEnvCombine::SRC_ALPHA);
1351            tec->setSource1_Alpha(osg::TexEnvCombine::TEXTURE);
1352            tec->setOperand1_Alpha(osg::TexEnvCombine::SRC_ALPHA);
1353
1354            stateset->setTextureAttributeAndModes(0, tec, osg::StateAttribute::OVERRIDE|osg::StateAttribute::ON);
1355
1356            stateset->setTextureMode(0,GL_TEXTURE_GEN_S,osg::StateAttribute::ON);
1357            stateset->setTextureMode(0,GL_TEXTURE_GEN_T,osg::StateAttribute::ON);
1358            stateset->setTextureMode(0,GL_TEXTURE_GEN_R,osg::StateAttribute::ON);
1359
1360            image_3d = normalmap_3d;
1361        }
1362    }
1363    else
1364    {     
1365        // set up the 3d texture itself,
1366        // note, well set the filtering up so that mip mapping is disabled,
1367        // gluBuild3DMipsmaps doesn't do a very good job of handled the
1368        // imbalanced dimensions of the 256x256x4 texture.
1369        osg::Texture3D* texture3D = new osg::Texture3D;
1370        texture3D->setResizeNonPowerOfTwoHint(false);
1371        texture3D->setFilter(osg::Texture3D::MIN_FILTER,minFilter);
1372        texture3D->setFilter(osg::Texture3D::MAG_FILTER, magFilter);
1373        texture3D->setWrap(osg::Texture3D::WRAP_R,osg::Texture3D::CLAMP_TO_EDGE);
1374        texture3D->setWrap(osg::Texture3D::WRAP_S,osg::Texture3D::CLAMP_TO_EDGE);
1375        texture3D->setWrap(osg::Texture3D::WRAP_T,osg::Texture3D::CLAMP_TO_EDGE);
1376        if (image_3d->getPixelFormat()==GL_ALPHA ||
1377            image_3d->getPixelFormat()==GL_LUMINANCE)
1378        {
1379            texture3D->setInternalFormatMode(osg::Texture3D::USE_USER_DEFINED_FORMAT);
1380            texture3D->setInternalFormat(GL_INTENSITY);
1381        }
1382        else
1383        {
1384            texture3D->setInternalFormatMode(internalFormatMode);
1385        }
1386
1387        texture3D->setImage(image_3d.get());
1388
1389        stateset->setTextureAttributeAndModes(0,texture3D,osg::StateAttribute::ON);
1390
1391        stateset->setTextureMode(0,GL_TEXTURE_GEN_S,osg::StateAttribute::ON);
1392        stateset->setTextureMode(0,GL_TEXTURE_GEN_T,osg::StateAttribute::ON);
1393        stateset->setTextureMode(0,GL_TEXTURE_GEN_R,osg::StateAttribute::ON);
1394
1395        stateset->setTextureAttributeAndModes(0,new osg::TexEnv(),osg::StateAttribute::ON);
1396    }
1397 
1398    return group;
1399}
1400
1401struct ScaleOperator
1402{
1403    ScaleOperator():_scale(1.0f) {}
1404    ScaleOperator(float scale):_scale(scale) {}
1405    ScaleOperator(const ScaleOperator& so):_scale(so._scale) {}
1406   
1407    ScaleOperator& operator = (const ScaleOperator& so) { _scale = so._scale; return *this; }
1408
1409    float _scale;
1410
1411    inline void luminance(float& l) const { l*= _scale; }
1412    inline void alpha(float& a) const { a*= _scale; }
1413    inline void luminance_alpha(float& l,float& a) const { l*= _scale; a*= _scale;  }
1414    inline void rgb(float& r,float& g,float& b) const { r*= _scale; g*=_scale; b*=_scale; }
1415    inline void rgba(float& r,float& g,float& b,float& a) const { r*= _scale; g*=_scale; b*=_scale; a*=_scale; }
1416};
1417
1418struct RecordRowOperator
1419{
1420    RecordRowOperator(unsigned int num):_colours(num),_pos(0) {}
1421
1422    mutable std::vector<osg::Vec4>  _colours;
1423    mutable unsigned int            _pos;
1424   
1425    inline void luminance(float l) const { rgba(l,l,l,1.0f); }
1426    inline void alpha(float a) const { rgba(1.0f,1.0f,1.0f,a); }
1427    inline void luminance_alpha(float l,float a) const { rgba(l,l,l,a);  }
1428    inline void rgb(float r,float g,float b) const { rgba(r,g,b,1.0f); }
1429    inline void rgba(float r,float g,float b,float a) const { _colours[_pos++].set(r,g,b,a); }
1430};
1431
1432struct WriteRowOperator
1433{
1434    WriteRowOperator():_pos(0) {}
1435    WriteRowOperator(unsigned int num):_colours(num),_pos(0) {}
1436
1437    std::vector<osg::Vec4>  _colours;
1438    mutable unsigned int    _pos;
1439   
1440    inline void luminance(float& l) const { l = _colours[_pos++].r(); }
1441    inline void alpha(float& a) const { a = _colours[_pos++].a(); }
1442    inline void luminance_alpha(float& l,float& a) const { l = _colours[_pos].r(); a = _colours[_pos++].a(); }
1443    inline void rgb(float& r,float& g,float& b) const { r = _colours[_pos].r(); g = _colours[_pos].g(); b = _colours[_pos].b(); }
1444    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(); }
1445};
1446
1447osg::Image* readRaw(int sizeX, int sizeY, int sizeZ, int numberBytesPerComponent, int numberOfComponents, const std::string& endian, const std::string& raw_filename)
1448{
1449    std::ifstream fin(raw_filename.c_str(), std::ifstream::binary);
1450    if (!fin) return 0;
1451
1452    GLenum pixelFormat;
1453    switch(numberOfComponents)
1454    {
1455        case 1 : pixelFormat = GL_LUMINANCE; break;
1456        case 2 : pixelFormat = GL_LUMINANCE_ALPHA; break;
1457        case 3 : pixelFormat = GL_RGB; break;
1458        case 4 : pixelFormat = GL_RGBA; break;
1459        default :
1460            osg::notify(osg::NOTICE)<<"Error: numberOfComponents="<<numberOfComponents<<" not supported, only 1,2,3 or 4 are supported."<<std::endl;
1461            return 0;
1462    }
1463
1464   
1465    GLenum dataType;
1466    switch(numberBytesPerComponent)
1467    {
1468        case 1 : dataType = GL_UNSIGNED_BYTE; break;
1469        case 2 : dataType = GL_UNSIGNED_SHORT; break;
1470        case 4 : dataType = GL_UNSIGNED_INT; break;
1471        default :
1472            osg::notify(osg::NOTICE)<<"Error: numberBytesPerComponent="<<numberBytesPerComponent<<" not supported, only 1,2 or 4 are supported."<<std::endl;
1473            return 0;
1474    }
1475   
1476    int s_maximumTextureSize=256, t_maximumTextureSize=256, r_maximumTextureSize=256;
1477   
1478    int sizeS = sizeX;
1479    int sizeT = sizeY;
1480    int sizeR = sizeZ;
1481    clampToNearestValidPowerOfTwo(sizeS, sizeT, sizeR, s_maximumTextureSize, t_maximumTextureSize, r_maximumTextureSize);
1482
1483    osg::ref_ptr<osg::Image> image = new osg::Image;
1484    image->allocateImage(sizeS, sizeT, sizeR, pixelFormat, dataType);
1485   
1486   
1487    bool endianSwap = (osg::getCpuByteOrder()==osg::BigEndian) ? (endian!="big") : (endian=="big");
1488   
1489    unsigned int r_offset = (sizeZ<sizeR) ? sizeR/2 - sizeZ/2 : 0;
1490   
1491    int offset = endianSwap ? numberBytesPerComponent : 0;
1492    int delta = endianSwap ? -1 : 1;
1493    for(int r=0;r<sizeZ;++r)
1494    {
1495        for(int t=0;t<sizeY;++t)
1496        {
1497            char* data = (char*) image->data(0,t,r+r_offset);
1498            for(int s=0;s<sizeX;++s)
1499            {
1500                if (!fin) return 0;
1501               
1502                for(int c=0;c<numberOfComponents;++c)
1503                {
1504                    char* ptr = data+offset;
1505                    for(int b=0;b<numberBytesPerComponent;++b)
1506                    {
1507                        fin.read((char*)ptr, 1);
1508                        ptr += delta;
1509                    }
1510                    data += numberBytesPerComponent;
1511                }
1512            }
1513        }
1514    }
1515
1516
1517    // normalise texture
1518    {
1519        // compute range of values
1520        osg::Vec4 minValue, maxValue;
1521        osgVolume::computeMinMax(image.get(), minValue, maxValue);
1522        osgVolume::modifyImage(image.get(),ScaleOperator(1.0f/maxValue.r()));
1523    }
1524   
1525   
1526    fin.close();
1527
1528    if (dataType!=GL_UNSIGNED_BYTE)
1529    {
1530        // need to convert to ubyte
1531       
1532        osg::ref_ptr<osg::Image> new_image = new osg::Image;
1533        new_image->allocateImage(sizeS, sizeT, sizeR, pixelFormat, GL_UNSIGNED_BYTE);
1534       
1535        RecordRowOperator readOp(sizeS);
1536        WriteRowOperator writeOp;
1537
1538        for(int r=0;r<sizeR;++r)
1539        {
1540            for(int t=0;t<sizeT;++t)
1541            {
1542                // reset the indices to beginning
1543                readOp._pos = 0;
1544                writeOp._pos = 0;
1545           
1546                // read the pixels into readOp's _colour array
1547                osgVolume::readRow(sizeS, pixelFormat, dataType, image->data(0,t,r), readOp);
1548                               
1549                // pass readOp's _colour array contents over to writeOp (note this is just a pointer swap).
1550                writeOp._colours.swap(readOp._colours);
1551               
1552                osgVolume::modifyRow(sizeS, pixelFormat, GL_UNSIGNED_BYTE, new_image->data(0,t,r), writeOp);
1553
1554                // return readOp's _colour array contents back to its rightful owner.
1555                writeOp._colours.swap(readOp._colours);
1556            }
1557        }
1558       
1559        image = new_image;
1560    }
1561   
1562    return image.release();
1563   
1564   
1565}
1566
1567enum ColourSpaceOperation
1568{
1569    NO_COLOUR_SPACE_OPERATION,
1570    MODULATE_ALPHA_BY_LUMINANCE,
1571    MODULATE_ALPHA_BY_COLOUR,
1572    REPLACE_ALPHA_WITH_LUMINACE
1573};
1574
1575struct ModulateAlphaByLuminanceOperator
1576{
1577    ModulateAlphaByLuminanceOperator() {}
1578
1579    inline void luminance(float&) const {}
1580    inline void alpha(float&) const {}
1581    inline void luminance_alpha(float& l,float& a) const { a*= l; }
1582    inline void rgb(float&,float&,float&) const {}
1583    inline void rgba(float& r,float& g,float& b,float& a) const { float l = (r+g+b)*0.3333333; a *= l;}
1584};
1585
1586struct ModulateAlphaByColourOperator
1587{
1588    ModulateAlphaByColourOperator(const osg::Vec4& colour):_colour(colour) { _lum = _colour.length(); }
1589   
1590    osg::Vec4 _colour;
1591    float _lum;
1592
1593    inline void luminance(float&) const {}
1594    inline void alpha(float&) const {}
1595    inline void luminance_alpha(float& l,float& a) const { a*= l*_lum; }
1596    inline void rgb(float&,float&,float&) const {}
1597    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()); }
1598};
1599
1600struct ReplaceAlphaWithLuminanceOperator
1601{
1602    ReplaceAlphaWithLuminanceOperator() {}
1603
1604    inline void luminance(float&) const {}
1605    inline void alpha(float&) const {}
1606    inline void luminance_alpha(float& l,float& a) const { a= l; }
1607    inline void rgb(float&,float&,float&) const { }
1608    inline void rgba(float& r,float& g,float& b,float& a) const { float l = (r+g+b)*0.3333333; a = l; }
1609};
1610
1611void doColourSpaceConversion(ColourSpaceOperation op, osg::Image* image, osg::Vec4& colour)
1612{
1613    switch(op)
1614    {
1615        case (MODULATE_ALPHA_BY_LUMINANCE):
1616            std::cout<<"doing conversion MODULATE_ALPHA_BY_LUMINANCE"<<std::endl;
1617            osgVolume::modifyImage(image,ModulateAlphaByLuminanceOperator());
1618            break;
1619        case (MODULATE_ALPHA_BY_COLOUR):
1620            std::cout<<"doing conversion MODULATE_ALPHA_BY_COLOUR"<<std::endl;
1621            osgVolume::modifyImage(image,ModulateAlphaByColourOperator(colour));
1622            break;
1623        case (REPLACE_ALPHA_WITH_LUMINACE):
1624            std::cout<<"doing conversion REPLACE_ALPHA_WITH_LUMINACE"<<std::endl;
1625            osgVolume::modifyImage(image,ReplaceAlphaWithLuminanceOperator());
1626            break;
1627        default:
1628            break;
1629    }
1630}
1631
1632
1633struct ApplyTransferFunctionOperator
1634{
1635    ApplyTransferFunctionOperator(osg::TransferFunction1D* tf, unsigned char* data):
1636        _tf(tf),
1637        _data(data) {}
1638       
1639    inline void luminance(float l) const
1640    {
1641        osg::Vec4 c = _tf->getInterpolatedValue(l);
1642        //std::cout<<"l = "<<l<<" c="<<c<<std::endl;
1643        *(_data++) = (unsigned char)(c[0]*255.0f + 0.5f);
1644        *(_data++) = (unsigned char)(c[1]*255.0f + 0.5f);
1645        *(_data++) = (unsigned char)(c[2]*255.0f + 0.5f);
1646        *(_data++) = (unsigned char)(c[3]*255.0f + 0.5f);
1647    }
1648     
1649    inline void alpha(float a) const
1650    {
1651        luminance(a);
1652    }
1653   
1654    inline void luminance_alpha(float l,float a) const
1655    {
1656        luminance(l);
1657    }
1658     
1659    inline void rgb(float r,float g,float b) const
1660    {
1661        luminance((r+g+b)*0.3333333);
1662    }
1663   
1664    inline void rgba(float r,float g,float b,float a) const
1665    {
1666        luminance(a);
1667    }
1668   
1669    mutable osg::ref_ptr<osg::TransferFunction1D> _tf;
1670    mutable unsigned char* _data;
1671};
1672
1673osg::Image* applyTransferFunction(osg::Image* image, osg::TransferFunction1D* transferFunction)
1674{
1675    std::cout<<"Applying transfer function"<<std::endl;
1676    osg::Image* output_image = new osg::Image;
1677    output_image->allocateImage(image->s(),image->t(), image->r(), GL_RGBA, GL_UNSIGNED_BYTE);
1678   
1679    ApplyTransferFunctionOperator op(transferFunction, output_image->data());
1680    osgVolume::readImage(image,op);
1681   
1682    return output_image;
1683}
1684
1685osg::TransferFunction1D* readTransferFunctionFile(const std::string& filename)
1686{
1687    std::string foundFile = osgDB::findDataFile(filename);
1688    if (foundFile.empty())
1689    {
1690        std::cout<<"Error: could not find transfer function file : "<<filename<<std::endl;
1691        return 0;
1692    }
1693   
1694    std::cout<<"Reading transfer function "<<filename<<std::endl;
1695
1696    osg::TransferFunction1D::ValueMap valueMap;
1697    std::ifstream fin(foundFile.c_str());
1698    while(fin)
1699    {
1700        float value, red, green, blue, alpha;
1701        fin >> value >> red >> green >> blue >> alpha;
1702        if (fin)
1703        {
1704            std::cout<<"value = "<<value<<" ("<<red<<", "<<green<<", "<<blue<<", "<<alpha<<")"<<std::endl;
1705            valueMap[value] = osg::Vec4(red,green,blue,alpha);
1706        }
1707    }
1708   
1709    if (valueMap.empty())
1710    {
1711        std::cout<<"Error: No values read from transfer function file: "<<filename<<std::endl;
1712        return 0;
1713    }
1714   
1715    osg::TransferFunction1D* tf = new osg::TransferFunction1D;
1716    tf->assign(valueMap, true);
1717   
1718    return tf;
1719}
1720
1721
1722class TestSupportOperation: public osg::GraphicsOperation
1723{
1724public:
1725
1726    TestSupportOperation():
1727        osg::GraphicsOperation("TestSupportOperation",false),
1728        supported(true),
1729        errorMessage(),
1730        maximumTextureSize(256) {}
1731
1732    virtual void operator () (osg::GraphicsContext* gc)
1733    {
1734        OpenThreads::ScopedLock<OpenThreads::Mutex> lock(mutex);
1735
1736        glGetIntegerv( GL_MAX_3D_TEXTURE_SIZE, &maximumTextureSize );
1737       
1738        osg::notify(osg::NOTICE)<<"Max texture size="<<maximumTextureSize<<std::endl;
1739    }
1740       
1741    OpenThreads::Mutex  mutex;
1742    bool                supported;
1743    std::string         errorMessage;
1744    GLint               maximumTextureSize;
1745};
1746
1747
1748
1749int main( int argc, char **argv )
1750{
1751    // use an ArgumentParser object to manage the program arguments.
1752    osg::ArgumentParser arguments(&argc,argv);
1753   
1754    // set up the usage document, in case we need to print out how to use this program.
1755    arguments.getApplicationUsage()->setDescription(arguments.getApplicationName()+" is the example which demonstrates use of 3D textures.");
1756    arguments.getApplicationUsage()->setCommandLineUsage(arguments.getApplicationName()+" [options] filename ...");
1757    arguments.getApplicationUsage()->addCommandLineOption("-h or --help","Display this information");
1758    arguments.getApplicationUsage()->addCommandLineOption("-n","Create normal map for per voxel lighting.");
1759    arguments.getApplicationUsage()->addCommandLineOption("-s <numSlices>","Number of slices to create.");
1760    arguments.getApplicationUsage()->addCommandLineOption("--images [filenames]","Specify a stack of 2d images to build the 3d volume from.");
1761    arguments.getApplicationUsage()->addCommandLineOption("--shader","Use OpenGL Shading Language. (default)");
1762    arguments.getApplicationUsage()->addCommandLineOption("--no-shader","Disable use of OpenGL Shading Language.");
1763    arguments.getApplicationUsage()->addCommandLineOption("--gpu-tf","Aply the transfer function on the GPU. (default)");
1764    arguments.getApplicationUsage()->addCommandLineOption("--cpu-tf","Apply the transfer function on the CPU.");
1765    arguments.getApplicationUsage()->addCommandLineOption("--mip","Use Maximum Intensity Projection (MIP) filtering.");
1766    arguments.getApplicationUsage()->addCommandLineOption("--xSize <size>","Relative width of rendered brick.");
1767    arguments.getApplicationUsage()->addCommandLineOption("--ySize <size>","Relative length of rendered brick.");
1768    arguments.getApplicationUsage()->addCommandLineOption("--zSize <size>","Relative height of rendered brick.");
1769    arguments.getApplicationUsage()->addCommandLineOption("--xMultiplier <multiplier>","Tex coord x mulitplier.");
1770    arguments.getApplicationUsage()->addCommandLineOption("--yMultiplier <multiplier>","Tex coord y mulitplier.");
1771    arguments.getApplicationUsage()->addCommandLineOption("--zMultiplier <multiplier>","Tex coord z mulitplier.");
1772    arguments.getApplicationUsage()->addCommandLineOption("--clip <ratio>","clip volume as a ratio, 0.0 clip all, 1.0 clip none.");
1773    arguments.getApplicationUsage()->addCommandLineOption("--maxTextureSize <size>","Set the texture maximum resolution in the s,t,r (x,y,z) dimensions.");
1774    arguments.getApplicationUsage()->addCommandLineOption("--s_maxTextureSize <size>","Set the texture maximum resolution in the s (x) dimension.");
1775    arguments.getApplicationUsage()->addCommandLineOption("--t_maxTextureSize <size>","Set the texture maximum resolution in the t (y) dimension.");
1776    arguments.getApplicationUsage()->addCommandLineOption("--r_maxTextureSize <size>","Set the texture maximum resolution in the r (z) dimension.");
1777    arguments.getApplicationUsage()->addCommandLineOption("--compressed","Enable the usage of compressed textures.");
1778    arguments.getApplicationUsage()->addCommandLineOption("--compressed-arb","Enable the usage of OpenGL ARB compressed textures.");
1779    arguments.getApplicationUsage()->addCommandLineOption("--compressed-dxt1","Enable the usage of S3TC DXT1 compressed textures.");
1780    arguments.getApplicationUsage()->addCommandLineOption("--compressed-dxt3","Enable the usage of S3TC DXT3 compressed textures.");
1781    arguments.getApplicationUsage()->addCommandLineOption("--compressed-dxt5","Enable the usage of S3TC DXT5 compressed textures.");
1782    arguments.getApplicationUsage()->addCommandLineOption("--modulate-alpha-by-luminance","For each pixel multiple the alpha value by the luminance.");
1783    arguments.getApplicationUsage()->addCommandLineOption("--replace-alpha-with-luminance","For each pixel mSet the alpha value to the luminance.");
1784    arguments.getApplicationUsage()->addCommandLineOption("--num-components <num>","Set the number of components to in he target image.");
1785//    arguments.getApplicationUsage()->addCommandLineOption("--raw <sizeX> <sizeY> <sizeZ> <numberBytesPerComponent> <numberOfComponents> <endian> <filename>","read a raw image data");
1786
1787    // construct the viewer.
1788    osgViewer::Viewer viewer(arguments);
1789
1790    // add the window size toggle handler
1791    viewer.addEventHandler(new osgViewer::WindowSizeHandler);
1792       
1793    {
1794        osg::ref_ptr<osgGA::KeySwitchMatrixManipulator> keyswitchManipulator = new osgGA::KeySwitchMatrixManipulator;
1795
1796        keyswitchManipulator->addMatrixManipulator( '1', "Trackball", new osgGA::TrackballManipulator() );
1797       
1798        osgGA::FlightManipulator* flightManipulator = new osgGA::FlightManipulator();
1799        flightManipulator->setYawControlMode(osgGA::FlightManipulator::NO_AUTOMATIC_YAW);
1800        keyswitchManipulator->addMatrixManipulator( '2', "Flight", flightManipulator );
1801
1802        viewer.setCameraManipulator( keyswitchManipulator.get() );
1803    }
1804
1805    // add the stats handler
1806    viewer.addEventHandler(new osgViewer::StatsHandler);
1807
1808    viewer.getCamera()->setClearColor(osg::Vec4(0.0f,0.0f,0.0f,0.0f));
1809
1810    // if user request help write it out to cout.
1811    if (arguments.read("-h") || arguments.read("--help"))
1812    {
1813        arguments.getApplicationUsage()->write(std::cout);
1814        return 1;
1815    }
1816
1817    std::string outputFile;
1818    while (arguments.read("-o",outputFile)) {}
1819
1820
1821
1822    osg::ref_ptr<osg::TransferFunction1D> transferFunction;
1823    std::string tranferFunctionFile;
1824    while (arguments.read("--tf",tranferFunctionFile))
1825    {
1826        transferFunction = readTransferFunctionFile(tranferFunctionFile);
1827    }
1828
1829    unsigned int numSlices=500;
1830    while (arguments.read("-s",numSlices)) {}
1831   
1832   
1833    float sliceEnd=1.0f;
1834    while (arguments.read("--clip",sliceEnd)) {}
1835
1836    float alphaFunc=0.02f;
1837    while (arguments.read("--alphaFunc",alphaFunc)) {}
1838
1839
1840   
1841    ShadingModel shadingModel = Standard;
1842   
1843    bool maximumIntensityProjection = false;
1844    while(arguments.read("--mip")) shadingModel =  MaximumIntensityProjection;
1845
1846    bool createNormalMap = false;
1847    while (arguments.read("-n"))
1848    {
1849        shadingModel = Light;
1850        createNormalMap=true;
1851    }
1852
1853    while (arguments.read("--isosurface"))
1854    {
1855        shadingModel = Isosurface;
1856    }
1857
1858    float xSize=1.0f, ySize=1.0f, zSize=1.0f;
1859    while (arguments.read("--xSize",xSize)) {}
1860    while (arguments.read("--ySize",ySize)) {}
1861    while (arguments.read("--zSize",zSize)) {}
1862
1863    float xMultiplier=1.0f, yMultiplier=1.0f, zMultiplier=1.0f;
1864    while (arguments.read("--xMultiplier",xMultiplier)) {}
1865    while (arguments.read("--yMultiplier",yMultiplier)) {}
1866    while (arguments.read("--zMultiplier",zMultiplier)) {}
1867
1868    osg::ref_ptr<TestSupportOperation> testSupportOperation = new TestSupportOperation;
1869    viewer.setRealizeOperation(testSupportOperation.get());
1870   
1871    viewer.realize();
1872
1873    int maximumTextureSize = testSupportOperation->maximumTextureSize;
1874    int s_maximumTextureSize = maximumTextureSize;
1875    int t_maximumTextureSize = maximumTextureSize;
1876    int r_maximumTextureSize = maximumTextureSize;
1877    while(arguments.read("--maxTextureSize",maximumTextureSize))
1878    {
1879        s_maximumTextureSize = maximumTextureSize;
1880        t_maximumTextureSize = maximumTextureSize;
1881        r_maximumTextureSize = maximumTextureSize;
1882    }
1883    while(arguments.read("--s_maxTextureSize",s_maximumTextureSize)) {}
1884    while(arguments.read("--t_maxTextureSize",t_maximumTextureSize)) {}
1885    while(arguments.read("--r_maxTextureSize",r_maximumTextureSize)) {}
1886
1887    osg::Texture::InternalFormatMode internalFormatMode = osg::Texture::USE_IMAGE_DATA_FORMAT;
1888    while(arguments.read("--compressed") || arguments.read("--compressed-arb")) { internalFormatMode = osg::Texture::USE_ARB_COMPRESSION; }
1889
1890    while(arguments.read("--compressed-dxt1")) { internalFormatMode = osg::Texture::USE_S3TC_DXT1_COMPRESSION; }
1891    while(arguments.read("--compressed-dxt3")) { internalFormatMode = osg::Texture::USE_S3TC_DXT3_COMPRESSION; }
1892    while(arguments.read("--compressed-dxt5")) { internalFormatMode = osg::Texture::USE_S3TC_DXT5_COMPRESSION; }
1893   
1894   
1895    // set up colour space operation.
1896    ColourSpaceOperation colourSpaceOperation = NO_COLOUR_SPACE_OPERATION;
1897    osg::Vec4 colourModulate(0.25f,0.25f,0.25f,0.25f);
1898    while(arguments.read("--modulate-alpha-by-luminance")) { colourSpaceOperation = MODULATE_ALPHA_BY_LUMINANCE; }
1899    while(arguments.read("--modulate-alpha-by-colour", colourModulate.x(),colourModulate.y(),colourModulate.z(),colourModulate.w() )) { colourSpaceOperation = MODULATE_ALPHA_BY_COLOUR; }
1900    while(arguments.read("--replace-alpha-with-luminance")) { colourSpaceOperation = REPLACE_ALPHA_WITH_LUMINACE; }
1901       
1902    bool resizeToPowerOfTwo = false;
1903   
1904    unsigned int numComponentsDesired = 0;
1905    while(arguments.read("--num-components", numComponentsDesired)) {}
1906
1907    bool useShader = true;
1908    while(arguments.read("--shader")) { useShader = true; }
1909    while(arguments.read("--no-shader")) { useShader = true; }
1910
1911    bool gpuTransferFunction = true;
1912    while(arguments.read("--gpu-tf")) { gpuTransferFunction = true; }
1913    while(arguments.read("--cpu-tf")) { gpuTransferFunction = false; }
1914
1915    typedef std::list< osg::ref_ptr<osg::Image> > Images;
1916    Images images;
1917
1918
1919    std::string vh_filename;
1920    while (arguments.read("--vh", vh_filename))
1921    {
1922        std::string raw_filename, transfer_filename;
1923        int xdim(0), ydim(0), zdim(0);
1924
1925        std::ifstream header(vh_filename.c_str());
1926        if (header)
1927        {
1928            header >> raw_filename >> transfer_filename >> xdim >> ydim >> zdim >> xSize >> ySize >> zSize;
1929        }
1930       
1931        if (xdim*ydim*zdim==0)
1932        {
1933            std::cout<<"Error in reading volume header "<<vh_filename<<std::endl;
1934            return 1;
1935        }
1936       
1937        if (!raw_filename.empty())
1938        {
1939            images.push_back(readRaw(xdim, ydim, zdim, 1, 1, "little", raw_filename));
1940        }
1941       
1942        if (!transfer_filename.empty())
1943        {
1944            std::ifstream fin(transfer_filename.c_str());
1945            if (fin)
1946            {
1947                osg::TransferFunction1D::ValueMap valueMap;
1948                float value = 0.0;
1949                while(fin && value<=1.0)
1950                {
1951                    float red, green, blue, alpha;
1952                    fin >> red >> green >> blue >> alpha;
1953                    if (fin)
1954                    {
1955                        valueMap[value] = osg::Vec4(red/255.0f,green/255.0f,blue/255.0f,alpha/255.0f);
1956                        std::cout<<"value = "<<value<<" ("<<red<<", "<<green<<", "<<blue<<", "<<alpha<<")";
1957                        std::cout<<"  ("<<valueMap[value]<<")"<<std::endl;
1958                    }
1959                    value += 1/255.0;
1960                }
1961
1962                if (valueMap.empty())
1963                {
1964                    std::cout<<"Error: No values read from transfer function file: "<<transfer_filename<<std::endl;
1965                    return 0;
1966                }
1967
1968                transferFunction = new osg::TransferFunction1D;
1969                transferFunction->assign(valueMap, true);
1970            }
1971        }
1972
1973    }
1974   
1975
1976    int sizeX, sizeY, sizeZ, numberBytesPerComponent, numberOfComponents;
1977    std::string endian, raw_filename;
1978    while (arguments.read("--raw", sizeX, sizeY, sizeZ, numberBytesPerComponent, numberOfComponents, endian, raw_filename))
1979    {
1980        images.push_back(readRaw(sizeX, sizeY, sizeZ, numberBytesPerComponent, numberOfComponents, endian, raw_filename));
1981    }
1982
1983    while (arguments.read("--images"))
1984    {
1985        ImageList imageList;
1986        for(int pos=1;pos<arguments.argc() && !arguments.isOption(pos);++pos)
1987        {
1988            // not an option so assume string is a filename.
1989            osg::Image *image = osgDB::readImageFile( arguments[pos]);
1990
1991            if(image)
1992            {
1993                imageList.push_back(image);
1994            }
1995        }
1996       
1997        // pack the textures into a single texture.
1998        ProcessRow processRow;
1999        images.push_back(createTexture3D(imageList, processRow, numComponentsDesired, s_maximumTextureSize, t_maximumTextureSize, r_maximumTextureSize, resizeToPowerOfTwo));
2000    }
2001
2002
2003    // any option left unread are converted into errors to write out later.
2004    arguments.reportRemainingOptionsAsUnrecognized();
2005
2006    // report any errors if they have occurred when parsing the program arguments.
2007    if (arguments.errors())
2008    {
2009        arguments.writeErrorMessages(std::cout);
2010        return 1;
2011    }
2012   
2013
2014    // assume remaining arguments are file names of textures.
2015    for(int pos=1;pos<arguments.argc();++pos)
2016    {
2017        if (!arguments.isOption(pos))
2018        {
2019            std::string filename = arguments[pos];
2020            if (osgDB::getLowerCaseFileExtension(filename)=="dicom")
2021            {
2022                images.push_back(osgDB::readImageFile( filename ));
2023            }
2024            else
2025            {
2026                osgDB::FileType fileType = osgDB::fileType(filename);
2027                if (fileType == osgDB::FILE_NOT_FOUND)
2028                {
2029                    filename = osgDB::findDataFile(filename);
2030                    fileType = osgDB::fileType(filename);
2031                }
2032
2033                if (fileType == osgDB::DIRECTORY)
2034                {
2035                    osgDB::DirectoryContents contents = osgDB::getDirectoryContents(filename);
2036
2037                    std::sort(contents.begin(), contents.end());
2038
2039                    ImageList imageList;
2040                    for(osgDB::DirectoryContents::iterator itr = contents.begin();
2041                        itr != contents.end();
2042                        ++itr)
2043                    {
2044                        std::string localFile = filename + "/" + *itr;
2045                        std::cout<<"contents = "<<localFile<<std::endl;
2046                        if (osgDB::fileType(localFile) == osgDB::REGULAR_FILE)
2047                        {
2048                            // not an option so assume string is a filename.
2049                            osg::Image *image = osgDB::readImageFile(localFile);
2050                            if(image)
2051                            {
2052                                imageList.push_back(image);
2053                            }
2054                        }
2055                    }
2056
2057                    // pack the textures into a single texture.
2058                    ProcessRow processRow;
2059                    images.push_back(createTexture3D(imageList, processRow, numComponentsDesired, s_maximumTextureSize, t_maximumTextureSize, r_maximumTextureSize, resizeToPowerOfTwo));
2060
2061                }
2062                else if (fileType == osgDB::REGULAR_FILE)
2063                {
2064                    // not an option so assume string is a filename.
2065                    images.push_back(osgDB::readImageFile( filename ));
2066                }
2067                else
2068                {
2069                    osg::notify(osg::NOTICE)<<"Error: could not find file: "<<filename<<std::endl;
2070                    return 1;
2071                }
2072            }           
2073        }
2074    }
2075   
2076    if (images.empty())
2077    {
2078        std::cout<<"No model loaded, please specify and volumetric image file on the command line."<<std::endl;
2079        return 1;
2080    }
2081
2082
2083    Images::iterator sizeItr = images.begin();
2084    xSize = (*sizeItr)->s();
2085    ySize = (*sizeItr)->t();
2086    zSize = (*sizeItr)->r();
2087    ++sizeItr;
2088
2089    for(;sizeItr != images.end(); ++sizeItr)
2090    {
2091        if ((*sizeItr)->s() != xSize ||
2092            (*sizeItr)->t() != ySize ||
2093            (*sizeItr)->r() != zSize)
2094        {
2095            std::cout<<"Images in sequence are not of the same dimensions."<<std::endl;
2096            return 1;
2097        }
2098    }
2099
2100
2101#if 1
2102    osg::RefMatrix* matrix = dynamic_cast<osg::RefMatrix*>(images.front()->getUserData());
2103    if (matrix)
2104    {
2105        osg::notify(osg::NOTICE)<<"Image has Matrix = "<<*matrix<<std::endl;
2106        xSize = xSize * (*matrix)(0,0);
2107        ySize = ySize * (*matrix)(1,1);
2108        zSize = zSize * (*matrix)(2,2);
2109    }
2110#endif
2111
2112    osg::Vec4 minValue, maxValue;
2113    bool computeMinMax = false;
2114    for(Images::iterator itr = images.begin();
2115        itr != images.end();
2116        ++itr)
2117    {
2118#if 0   
2119        osg::RefMatrix* matrix = dynamic_cast<osg::RefMatrix*>((*itr)->getUserData());
2120        if (matrix)
2121        {
2122            std::cout<<"matrix = "<<*matrix<<std::endl;
2123        }
2124#endif
2125        if (osgVolume::computeMinMax(itr->get(), minValue, maxValue)) computeMinMax = true;
2126    }
2127   
2128   
2129    if (computeMinMax)
2130    {
2131        osg::notify(osg::NOTICE)<<"Min value "<<minValue<<std::endl;
2132        osg::notify(osg::NOTICE)<<"Max value "<<maxValue<<std::endl;
2133
2134        float minComponent = minValue[0];
2135        minComponent = osg::minimum(minComponent,minValue[1]);
2136        minComponent = osg::minimum(minComponent,minValue[2]);
2137        minComponent = osg::minimum(minComponent,minValue[3]);
2138
2139        float maxComponent = maxValue[0];
2140        maxComponent = osg::maximum(maxComponent,maxValue[1]);
2141        maxComponent = osg::maximum(maxComponent,maxValue[2]);
2142        maxComponent = osg::maximum(maxComponent,maxValue[3]);
2143
2144        float scale = 0.99f/(maxComponent-minComponent);
2145        float offset = -minComponent * scale;
2146
2147
2148        for(Images::iterator itr = images.begin();
2149            itr != images.end();
2150            ++itr)
2151        {       
2152            osgVolume::offsetAndScaleImage(itr->get(),
2153                osg::Vec4(offset, offset, offset, offset),
2154                osg::Vec4(scale, scale, scale, scale));
2155        }
2156    }
2157
2158   
2159    if (colourSpaceOperation!=NO_COLOUR_SPACE_OPERATION)
2160    {
2161        for(Images::iterator itr = images.begin();
2162            itr != images.end();
2163            ++itr)
2164        {       
2165            doColourSpaceConversion(colourSpaceOperation, itr->get(), colourModulate);
2166        }
2167    }
2168   
2169    if (!gpuTransferFunction && transferFunction.valid())
2170    {
2171        for(Images::iterator itr = images.begin();
2172            itr != images.end();
2173            ++itr)
2174        {       
2175            *itr = applyTransferFunction(itr->get(), transferFunction.get());
2176        }
2177    }
2178   
2179    osg::ref_ptr<osg::Image> image_3d = 0;
2180
2181    if (images.size()==1)
2182    {
2183        osg::notify(osg::NOTICE)<<"Single image "<<images.size()<<" volumes."<<std::endl;
2184        image_3d = images.front();
2185    }
2186    else
2187    {
2188        osg::notify(osg::NOTICE)<<"Creating sequence of "<<images.size()<<" volumes."<<std::endl;
2189   
2190        osg::ref_ptr<osg::ImageSequence> imageSequence = new osg::ImageSequence;
2191        image_3d = imageSequence.get();
2192        for(Images::iterator itr = images.begin();
2193            itr != images.end();
2194            ++itr)
2195        {       
2196            imageSequence->addImage(itr->get());
2197        }
2198        imageSequence->play();
2199    }
2200   
2201    osg::ref_ptr<osg::Image> normalmap_3d = 0;
2202    if (createNormalMap)
2203    {
2204        if (images.size()==1)
2205        {
2206            normalmap_3d = createNormalMapTexture(images.front().get());
2207        }
2208        else
2209        {
2210            osg::ref_ptr<osg::ImageSequence> normalmapSequence = new osg::ImageSequence;
2211            normalmap_3d = normalmapSequence.get();
2212            for(Images::iterator itr = images.begin();
2213                itr != images.end();
2214                ++itr)
2215            {       
2216                normalmapSequence->addImage(createNormalMapTexture(itr->get()));
2217            }
2218            normalmapSequence->play();
2219        }
2220    }
2221
2222    // create a model from the images.
2223    osg::Node* rootNode = 0;
2224   
2225    if (useShader)
2226    {
2227        rootNode = createShaderModel(shadingModel,
2228                               image_3d, normalmap_3d.get(),
2229                               (gpuTransferFunction ? transferFunction.get() : 0),
2230                               internalFormatMode,
2231                               xSize, ySize, zSize,
2232                               xMultiplier, yMultiplier, zMultiplier,
2233                               numSlices, sliceEnd, alphaFunc);
2234    }
2235    else
2236    {
2237        rootNode = createModel(shadingModel,
2238                               image_3d, normalmap_3d,
2239                               internalFormatMode,
2240                               xSize, ySize, zSize,
2241                               xMultiplier, yMultiplier, zMultiplier,
2242                               numSlices, sliceEnd, alphaFunc);
2243    }
2244   
2245    if (!outputFile.empty())
2246    {   
2247        std::string ext = osgDB::getFileExtension(outputFile);
2248        std::string name_no_ext = osgDB::getNameLessExtension(outputFile);
2249        if (ext=="osg")
2250        {
2251            if (image_3d.valid())
2252            {
2253                image_3d->setFileName(name_no_ext + ".dds");           
2254                osgDB::writeImageFile(*image_3d, image_3d->getFileName());
2255            }
2256            if (normalmap_3d.valid())
2257            {
2258                normalmap_3d->setFileName(name_no_ext + "_normalmap.dds");           
2259                osgDB::writeImageFile(*normalmap_3d, normalmap_3d->getFileName());
2260            }
2261           
2262            osgDB::writeNodeFile(*rootNode, outputFile);
2263        }
2264        else if (ext=="ive")
2265        {
2266            osgDB::writeNodeFile(*rootNode, outputFile);       
2267        }
2268        else if (ext=="dds")
2269        {
2270            osgDB::writeImageFile(*image_3d, outputFile);       
2271        }
2272        else
2273        {
2274            std::cout<<"Extension not support for file output, not file written."<<std::endl;
2275        }
2276       
2277        return 0;
2278    }
2279
2280
2281    if (rootNode)
2282    {
2283
2284        // set the scene to render
2285        viewer.setSceneData(rootNode);
2286       
2287        // the the viewers main frame loop
2288        viewer.run();
2289    }   
2290   
2291    return 0;
2292
2293}
Note: See TracBrowser for help on using the browser.