source: osgVisual/src/sky_Silverlining/visual_skySilverLining.cpp @ 130

Last change on this file since 130 was 130, checked in by Torben Dannhauer, 14 years ago

Upgraded to Silverlining 2.0

Todo:
Re-enable DEBUg build
correct cloud altitude errors
enable the new ray casting cloud type
File size: 26.2 KB
Line 
1/* -*-c++-*- osgVisual - Copyright (C) 2009-2010 Torben Dannhauer
2 *
3 * This library is based on OpenSceneGraph, open source and may be redistributed and/or modified under
4 * the terms of the OpenSceneGraph Public License (OSGPL) version 0.0 or
5 * (at your option) any later version.  The full license is in LICENSE file
6 * included with this distribution, and on the openscenegraph.org website.
7 *
8 * osgVisual requires for some proprietary modules a license from the correspondig manufacturer.
9 * You have to aquire licenses for all used proprietary modules.
10 *
11 * This library is distributed in the hope that it will be useful,
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14 * OpenSceneGraph Public License for more details.
15*/
16
17#include <visual_skySilverLining.h>
18
19using namespace osgVisual;
20
21visual_skySilverLining::visual_skySilverLining(osgViewer::Viewer* viewer_)
22{
23        OSG_NOTIFY( osg::ALWAYS ) << "Initialize visual_skySilverlining..." << std::endl;
24
25        atmosphereInitialized = false;
26        postInitialized = false;
27        atmosphere = NULL;
28        viewer = viewer_;
29        lat = 0;
30        lon = 0;
31        upVectorLat = 5;        // Just to ensure deltaLatToUpdate is big enough to enforce the sky to update at the first frame
32        upVectorLon = 5;        // Just to ensure deltaLonToUpdate is big enough to enforce the sky to update at the first frame
33        for( int i=0; i<MAX_CLOUDLAYER_SLOTS; i++ )
34        {
35                cloudLayerSlots.push_back( cloudLayerSlot() );
36                cloudLayerSlots.back().slot_id = i;
37        }
38}
39
40visual_skySilverLining::~visual_skySilverLining(void)
41{
42        this->removeUpdateCallback( updateCallback );
43        if ( atmosphere != NULL )
44                delete atmosphere;
45}
46
47void visual_skySilverLining::skyUpdateCallback::operator()(osg::Node* node, osg::NodeVisitor* nv)
48{
49        //std::cout << "Sky silverlining update callback" << std::endl;
50        // Check if atmosphere is initialized.
51        if (!sky->isInitialized())
52                return;
53
54        // on first time: perform Post-Init.
55        sky->postInit();
56
57        // Update sky
58        double lat, lon, height;
59        util::getWGS84ofCamera( sceneCamera, csn, lat, lon, height );
60        //std::cout << "lat: " << osg::RadiansToDegrees(lat) << ", lon: " << osg::RadiansToDegrees(lon) << ", height: " << height << std::endl;
61        sky->setLocation( lat, lon, height );
62}
63
64void visual_skySilverLining::setDateTime( int year_, int month_, int day_, int hour_, int minute_, int second_, bool daylightSaving_, double timezoneOffset_ )
65{
66        // Check if atmosphere is initialized.
67        if (!isInitialized())
68                return;
69
70        SilverLining::LocalTime t = atmosphere->GetConditions()->GetTime();
71        t.SetYear( year_ );
72        t.SetMonth( month_ );
73        t.SetDay( day_ );
74        t.SetHour( hour_ );
75        t.SetMinutes( minute_ );
76        t.SetSeconds( second_ );
77        t.SetObservingDaylightSavingsTime( daylightSaving_ );
78        t.SetTimeZone( timezoneOffset_ );
79}
80
81void visual_skySilverLining::setTime( int hour_, int minute_, int second_ )
82{
83        // Check if atmosphere is initialized.
84        if (!isInitialized())
85                return;
86       
87        SilverLining::LocalTime t = atmosphere->GetConditions()->GetTime();
88        t.SetHour( hour_ );
89        t.SetMinutes( minute_ );
90        t.SetSeconds( second_ );
91        atmosphere->GetConditions()->SetTime( t );
92}
93
94void visual_skySilverLining::setDate( int year_, int month_, int day_ )
95{
96        // Check if atmosphere is initialized.
97        if (!isInitialized())
98                return;
99
100        SilverLining::LocalTime t = atmosphere->GetConditions()->GetTime();
101        t.SetYear( year_ );
102        t.SetMonth( month_ );
103        t.SetDay( day_ );
104}
105
106void visual_skySilverLining::setDateByEpoch( int secondsSince1970_ )
107{
108        // Check if atmosphere is initialized.
109        if (!isInitialized())
110                return;
111       
112        SilverLining::LocalTime t = atmosphere->GetConditions()->GetTime();
113        t.SetFromEpochSeconds( secondsSince1970_ );
114}
115
116void visual_skySilverLining::setLocation(double lat_, double lon_, double alt_)
117{
118        // Check if atmosphere is initialized.
119        if (!isInitialized())
120                return;
121
122        lat = lat_;
123        lon = lon_;
124        height = alt_;
125       
126        SilverLining::Location loc = atmosphere->GetConditions()->GetLocation();
127        loc.SetAltitude(alt_);
128        loc.SetLatitude(lat_);
129        loc.SetLongitude(lon_);         
130        atmosphere->GetConditions()->SetLocation( loc );
131
132        updateUpVector();
133}
134
135
136bool visual_skySilverLining::isInitialized()
137{
138        // Check if atmosphere is initialized. If not: make a deep lookup. If initialized, perform. Otherwise return
139        if (!atmosphereInitialized)
140        {
141                skySilverLining_atmosphereReference *ar = dynamic_cast<skySilverLining_atmosphereReference *>(viewer->getCamera()->getUserData());
142                if (ar != NULL )
143                {
144                        if (ar->atmosphereInitialized)
145                                atmosphereInitialized = true;
146                }
147        }
148        return(atmosphereInitialized);
149}
150
151void visual_skySilverLining::init(osg::Group *distortedRoot, osg::CoordinateSystemNode *sceneGraphRoot)
152{
153        sceneRoot = sceneGraphRoot;
154
155        // Use projection matrix callback oder fixed Cullsettings?
156        bool useProjMatrixCallback = true;
157
158        // add Sky to SceneGraphRoot
159        sceneGraphRoot->addChild( this );
160
161        // Deactivate culling for the sky node (required by the silverlining sky framework)
162        this->setCullingActive(false);
163
164        // Instantiate an Atmosphere and associate it with this camera. If you have multiple cameras
165        // in multiple contexts, be sure to instantiate seperate Atmosphere objects for each.
166    // ***IMPORTANT!**** Check that the path to the resources folder for SilverLining in SkyDrawable.cpp
167    // SkyDrawable::initializeSilverLining matches with where you installed SilverLining.
168        atmosphere = new SilverLining::Atmosphere(SILVERLINING_LICENSEE, SILVERLINING_LICENSE);
169
170    // Add the sky (calls Atmosphere::BeginFrame and handles initialization once you're in
171    // the rendering thread)
172        skyDrawable = new skySilverLining_skyDrawable(viewer, sceneRoot);
173
174        if(distortedRoot)       // if distortion used:
175        {
176                int rootKids = distortedRoot->getNumChildren();
177                for (int i = 0; i < rootKids; i++)
178                {
179                        osg::Node *n = distortedRoot->getChild(i);
180                        osg::Camera *cam = dynamic_cast<osg::Camera*>(n);
181                        if (cam && cam->getRenderOrder() == osg::Camera::PRE_RENDER)
182                                sceneCamera = cam;
183                }
184        }
185        else    // if no distortion used:
186                sceneCamera = viewer->getCamera();
187
188        osg::Camera *mainCamera = viewer->getCamera();
189        if (!useProjMatrixCallback)
190        {
191                mainCamera->setClearMask(0);
192                mainCamera->setComputeNearFarMode(osg::CullSettings::DO_NOT_COMPUTE_NEAR_FAR);
193                double fovy, aspect, zNear, zFar;
194                mainCamera->getProjectionMatrixAsPerspective(fovy, aspect, zNear, zFar);
195                mainCamera->setProjectionMatrixAsPerspective(fovy, aspect, 2, 125000);
196        }
197        else
198        {
199                cb = new skySilverLining_projectionMatrixCallback( atmosphere, viewer->getCamera(), sceneRoot);
200                sceneCamera->setClampProjectionMatrixCallback(cb);
201                cb->setSkyDrawable(skyDrawable);
202        }
203
204        // append atmosphere pointer to the cameras.
205        sceneCamera->setClearMask(0);
206        osg::ref_ptr<skySilverLining_atmosphereReference> ar = new skySilverLining_atmosphereReference;
207        ar->atmosphere = atmosphere;
208        sceneCamera->setUserData(ar);
209        mainCamera->setUserData(ar);
210
211        // Create and install updateCallback (for position etc.)
212        updateCallback = new skyUpdateCallback( sceneGraphRoot, sceneCamera, this );
213        this->setUpdateCallback( updateCallback );
214
215    // Use a RenderBin to enforce that the sky gets drawn first, then the scene, then the clouds
216        skyDrawable->getOrCreateStateSet()->setRenderBinDetails(-1, "RenderBin");
217
218    // Add the models
219    sceneGraphRoot->getOrCreateStateSet()->setRenderBinDetails(1, "RenderBin");
220
221    // Add the clouds (note, you need this even if you don't want clouds - it calls
222    // Atmosphere::EndFrame() )
223        cloudsDrawable = new skySilverLining_cloudsDrawable(viewer);
224        cloudsDrawable->getOrCreateStateSet()->setRenderBinDetails(99, "RenderBin");
225
226        // Add drawable to this geode to get rendered
227        this->addDrawable(skyDrawable);
228        this->addDrawable(cloudsDrawable);
229
230        //SilverLining::Atmosphere::EnableHDR( true );
231}
232
233void visual_skySilverLining::postInit()
234{
235        // Only allow one execution
236        //if(postInitialized)
237        //      return;
238        //else postInitialized = true;
239
240        // Execute Updatecallback once before adding Clouds.
241        //updateCallback->operator ()(this, NULL);
242
243        //atmosphere->GetConditions()->SetFog( 0.8, 1, 1, 1);   // use this for simulation real fog.
244
245        //Todo: secure memory-manager of timer*. oder remove paragraph
246        //MyMillisecondTimer *timer = new MyMillisecondTimer();
247 //   atmosphere->GetConditions()->SetMillisecondTimer(timer);
248        //atmosphere->GetConditions()->EnableTimePassage(true, -1);
249}
250
251void visual_skySilverLining::updateUpVector()
252{
253        // Allowed deltaLat and deltaLon between Updating the upvector
254        double deltaLatToUpdate = osg::DegreesToRadians( 0.1 ); // 0.5 deg
255        double deltaLonToUpdate = osg::DegreesToRadians( 0.1 ); // 0.5 deg
256
257        if ( fabs(lat-upVectorLat) > deltaLatToUpdate || fabs(lon-upVectorLon) > deltaLonToUpdate )
258        {
259                // Get ellipsoid model
260                osg::EllipsoidModel* ellipsoid = sceneRoot->getEllipsoidModel();
261                if ( !ellipsoid )
262                        return;
263
264                //OSG_NOTIFY( osg::ALWAYS ) << "update Upvector.." << std::endl;
265
266                // Calculate up vector
267                double x,y,z;
268                util::getXYZofCamera(sceneCamera, x, y, z);
269                osg::Vec3d upVector = ellipsoid->computeLocalUpVector( x, y, z );
270                upVector.normalize();
271               
272                // Calculate side vector
273                osg::Matrixd localToWorld;
274                ellipsoid->computeLocalToWorldTransformFromLatLongHeight(lat, lon, height, localToWorld);
275
276                osg::Vec3d sideVector(localToWorld(0, 0), localToWorld(0, 1), localToWorld(0, 2));
277                sideVector.normalize();
278
279                // Update silverlining vectors
280                //std::cout << "upVector: X: " << upVector.x() << ", Y: "<<  upVector.y() << ", Z: "<< upVector.z() << std::endl;
281                //std::cout << "sideVector: X: " << sideVector.x() << ", Y: "<<  sideVector.y() << ", Z: "<< sideVector.z() << std::endl;
282                atmosphere->SetUpVector( upVector.x(), upVector.y(), upVector.z() );
283                atmosphere->SetRightVector( sideVector.x(), sideVector.y(), sideVector.z() );
284
285                // Note new upvector lat/lon
286                upVectorLon = lon;
287                upVectorLat = lat;
288        }       // If update requiered END
289}
290
291void visual_skySilverLining::shutdown()
292{
293        if (isInitialized())
294        {
295                // Remove this Node from scenegraph
296                sceneRoot->removeChild( this );
297               
298                // Remove updatecallback
299                this->removeUpdateCallback( updateCallback );
300                updateCallback = NULL;
301
302                // delete drawables
303                skyDrawable->shutdown();
304                this->removeDrawable(skyDrawable);
305                this->removeDrawable(cloudsDrawable);
306        }
307}
308
309void visual_skySilverLining::setVisibility(double visibility_)
310{
311        if (isInitialized())
312        {
313                atmosphere->GetConditions()->SetVisibility( visibility_ );
314        }
315
316}
317
318double visual_skySilverLining::getVisibility()
319{
320        if (isInitialized())
321        {
322                return atmosphere->GetConditions()->GetVisibility();
323        }
324        else
325                return -1;
326}
327
328void visual_skySilverLining::setTurbidity(double turbidity_)
329{
330        if (isInitialized())
331        {
332                atmosphere->GetConditions()->SetTurbidity( turbidity_ );
333        }
334}
335
336double visual_skySilverLining::getTurbidity()
337{
338        if (isInitialized())
339        {
340                return atmosphere->GetConditions()->GetTurbidity();
341        }
342        else
343                return -1;
344}
345
346void visual_skySilverLining::clearAllWindVolumes()
347{
348        if (isInitialized())
349        {
350                atmosphere->GetConditions()->ClearWindVolumes();
351        }
352}
353
354bool visual_skySilverLining::insideWind(double height_, double& bottom_, double& top_, double& speed_, double& direction_)
355{
356        if (isInitialized())
357        {
358                // Calculation earth radius on the wind positionl approximated through the now used position.
359                double radius;
360                if ( util::calculateEarthRadiusAtWGS84Coordinate(lat, lon, sceneRoot, radius) )
361                {
362                        // go through all wind volumes an check them for inside()
363                        //std::map<int, SilverLining::WindVolume> windvolumes = atmosphere->GetConditions()->GetWindVolumes();
364                        SL_MAP(int, SilverLining::WindVolume) windvolumes = atmosphere->GetConditions()->GetWindVolumes();
365                        for(unsigned int i=0; i<windvolumes.size(); i++)
366                        {
367                                if( windvolumes[i].Inside(radius + height_) )
368                                {
369                                        // save wind data
370                                        bottom_ = windvolumes[i].GetMinAltitude();
371                                        top_ = windvolumes[i].GetMaxAltitude() - radius;
372                                        speed_ = windvolumes[i].GetWindSpeed() - radius;
373                                        direction_ = windvolumes[i].GetDirection();
374                                        // return that wind was found
375                                        return true;
376                                }
377                        }       // For END
378                }       // If valid radius END
379        }       // If initialized() END
380        bottom_ = -1;
381        top_ = -1;
382        speed_ = -1;
383        direction_ = -1;
384        return false;
385}
386
387void visual_skySilverLining::addWindVolume(double bottom_, double top_, double speed_, int direction_)
388{
389        if (isInitialized())           
390        {
391                // Calculation earth radius on current lat lon position
392                double radius;
393                if ( util::calculateEarthRadiusAtWGS84Coordinate(lat, lon, sceneRoot, radius) )
394                {
395                        // correct wind value:
396                        if ( direction_ < 180 )
397                                direction_ += 180;
398                        else direction_ -= 180;
399
400                        // Setting up Wind
401                        SilverLining::WindVolume wv;
402                        wv.SetDirection( direction_ );
403                        wv.SetMinAltitude( radius + bottom_ );
404                        wv.SetMaxAltitude( radius + top_ );
405                        wv.SetWindSpeed( speed_ );
406                        atmosphere->GetConditions()->SetWind(wv);
407                }
408        }
409}
410
411void visual_skySilverLining::setLightPollution(double lightPollution_)
412{
413        if (isInitialized())
414        {
415                return atmosphere->GetConditions()->SetLightPollution( lightPollution_ );
416        }
417}
418
419double visual_skySilverLining::getLightPollution()
420{
421        if (isInitialized())
422        {
423                return atmosphere->GetConditions()->GetLightPollution();
424        }
425        else
426                return -1;
427
428}
429
430void visual_skySilverLining::addCloudLayer(int slot_, double baseLength_, double baseWidth_, double thickness_, double baseHeight_, double density_, CloudTypes cloudtype_ )
431{
432        if (isInitialized() && slot_ >= 0 && slot_ < MAX_CLOUDLAYER_SLOTS)
433        {
434                // Calculation earth radius on current lat lon position
435                double radius;
436                if ( util::calculateEarthRadiusAtWGS84Coordinate(lat, lon, sceneRoot, radius) )
437                {
438                        OSG_NOTIFY(osg::ALWAYS) << "radius:" << radius << std::endl;
439                        // generate Cloud Layer
440                        SilverLining::CloudLayer *cloudLayer_;
441                        cloudLayer_ = SilverLining::CloudLayerFactory::Create(cloudtype_);
442                        cloudLayer_->SetBaseAltitude( baseHeight_ + radius); // todo
443                        //cloudLayer_->SetBaseAltitude( radius - 10000); // todo
444                                               
445                        cloudLayer_->SetThickness(thickness_);
446                        cloudLayer_->SetBaseLength(baseLength_);
447                        cloudLayer_->SetBaseWidth(baseWidth_);
448                        cloudLayer_->SetDensity(density_);
449                        cloudLayer_->SetLayerPosition( 0.0, 0.0 );
450                        cloudLayer_->GenerateShadowMaps(false);
451
452                        //Save cloudlayer information into SLOT
453                        cloudLayerSlots[slot_].used = true;
454                        cloudLayerSlots[slot_].enabled = true;
455                        cloudLayerSlots[slot_].cloudLayerPointer = cloudLayer_;
456                        switch(cloudtype_)
457                        {
458                        case 0: cloudLayerSlots[slot_].typeName = "CIRROCUMULUS";                               // High planar cloud puffs             
459                                break;
460                        case 1: cloudLayerSlots[slot_].typeName = "CIRRUS_FIBRATUS";                    // High, thicker and fibrous clouds that signal changing weather
461                        break;
462                        case 2: cloudLayerSlots[slot_].typeName = "STRATUS";                                    // Low clouds represented as a slab
463                        break;
464                        case 3: cloudLayerSlots[slot_].typeName = "CUMULUS_MEDIOCRIS";          // Low, puffy clouds on fair days
465                        break;
466                        case 4: cloudLayerSlots[slot_].typeName = "CUMULUS_CONGESTUS";          // Large cumulus clouds that could turn into a thunderhead
467                        break;
468                        case 5: cloudLayerSlots[slot_].typeName = "CUMULONIMBUS_CAPPILATUS";    // Big storm clouds.
469                        break;
470                        case 6: cloudLayerSlots[slot_].typeName = "STRATOCUMULUS";                              // Low, dense, puffy clouds with some sun breaks between them.
471                                break;
472                        default: OSG_NOTIFY( osg::FATAL ) << "visual_skySilverlining::addCloudLayer - Invalid cloud type." << std::cout;
473                                 break;
474                        };
475               
476                        // Pass cloudlayer for seeding and adding to sky framework to the render thread
477                        skyDrawable->addCloudLayer( &cloudLayerSlots[slot_] );
478                } // If valid radius END
479        }        // If isInitialized() END
480}
481
482void visual_skySilverLining::removeCloudLayer( int slot_ )
483{
484        if (isInitialized() && slot_ >= 0 && slot_ < MAX_CLOUDLAYER_SLOTS)
485        {
486                atmosphere->GetConditions()->RemoveCloudLayer( cloudLayerSlots[slot_].cloudLayerHandle );
487                cloudLayerSlots[slot_].used = false;
488                cloudLayerSlots[slot_].cloudLayerHandle = -1;
489                cloudLayerSlots[slot_].cloudLayerPointer = NULL;
490                cloudLayerSlots[slot_].enabled = false;
491        }
492}
493
494void visual_skySilverLining::clearAllSlots()
495{
496        if (isInitialized())
497        {
498                atmosphere->GetConditions()->RemoveAllCloudLayers();
499                for( int i=0; i<MAX_CLOUDLAYER_SLOTS; i++ )
500                {
501                        cloudLayerSlots[i].used = false;
502                        cloudLayerSlots[i].cloudLayerHandle = -1;
503                        cloudLayerSlots[i].cloudLayerPointer = NULL;
504                        cloudLayerSlots[i].enabled = false;
505                }
506        }
507}
508
509SilverLining::CloudLayer* visual_skySilverLining::getCloudLayer( int slot_ )
510{
511        if (isInitialized() && slot_ >= 0 && slot_ < MAX_CLOUDLAYER_SLOTS)
512        {
513                if ( cloudLayerSlots[slot_].used )
514                        return cloudLayerSlots[slot_].cloudLayerPointer; 
515        }
516
517        return NULL;
518}
519
520void visual_skySilverLining::setEnabled(int slot_, bool enabled_ )
521{
522        if (isInitialized() && slot_ >= 0 && slot_ < MAX_CLOUDLAYER_SLOTS)
523        {
524                cloudLayerSlots[slot_].enabled = enabled_;
525        }
526}
527
528bool visual_skySilverLining::isEnabled( int slot_ )
529{
530        if (isInitialized() && slot_ >= 0 && slot_ < MAX_CLOUDLAYER_SLOTS)
531        {
532                return cloudLayerSlots[slot_].enabled;
533        }
534        return false;
535}
536
537void visual_skySilverLining::fadeVisibility( int slot_, int fadetimeMS_ )
538{
539        if (isInitialized() && slot_ >= 0 && slot_ < MAX_CLOUDLAYER_SLOTS)
540        {
541                if (cloudLayerSlots[slot_].enabled)
542                {
543                        cloudLayerSlots[slot_].enabled = false;
544                        cloudLayerSlots[slot_].cloudLayerPointer->SetEnabled( false, fadetimeMS_ );
545                }
546                else
547                {
548                        cloudLayerSlots[slot_].enabled = true;
549                        cloudLayerSlots[slot_].cloudLayerPointer->SetEnabled( true, fadetimeMS_ );
550                }
551        }
552}
553
554std::string visual_skySilverLining::getCloudLayerTypeName( int slot_ )
555{
556        if (isInitialized() && slot_ >= 0 && slot_ < MAX_CLOUDLAYER_SLOTS)
557        {
558                return cloudLayerSlots[slot_].typeName;
559        }
560        return "";
561}
562
563void visual_skySilverLining::clearGlobalPrecipitation()
564{
565        if (isInitialized())
566        {
567                atmosphere->GetConditions()->SetPrecipitation( SilverLining::CloudLayer::NONE, 0.0 );
568        }
569}
570
571void visual_skySilverLining::setGlobalPrecipitation( double rate_mmPerHour_rain_, double rate_mmPerHour_drySnow_, double rate_mmPerHour_wetSnow_, double rate_mmPerHour_sleet_ )
572{
573        if ( isInitialized() )
574        {
575                // Delete old Precipitation
576                atmosphere->GetConditions()->SetPrecipitation( SilverLining::CloudLayer::NONE, 0.0 );
577               
578                // Set new Precipitation
579                atmosphere->GetConditions()->SetPrecipitation( SilverLining::CloudLayer::RAIN, rate_mmPerHour_rain_ );
580                atmosphere->GetConditions()->SetPrecipitation( SilverLining::CloudLayer::WET_SNOW, rate_mmPerHour_drySnow_ );
581                atmosphere->GetConditions()->SetPrecipitation( SilverLining::CloudLayer::DRY_SNOW, rate_mmPerHour_wetSnow_ );
582                atmosphere->GetConditions()->SetPrecipitation( SilverLining::CloudLayer::SLEET, rate_mmPerHour_sleet_ );
583        }
584}
585       
586bool visual_skySilverLining::getOverallPrecipitationAtLocation( double& rate_mmPerHour_rain, double& rate_mmPerHour_drySnow, double& rate_mmPerHour_wetSnow, double& rate_mmPerHour_sleet, double lat_, double lon_, double height_ )
587{
588        if (isInitialized())
589        {
590                // Init
591                bool hasPrecipitation = false;
592                double x = 0;
593                double y = 0;
594                double z = 0;
595                //// If -1 : Use sky internal values
596                if ( lat_ == -1 )
597                        lat_ = lat;
598                if ( lon_ == -1 )
599                        lon_ = lon;
600                if ( height_ == -1 )
601                        height_ = height;
602
603                // Set precipitation to zero;
604                rate_mmPerHour_rain = 0;
605                rate_mmPerHour_drySnow = 0;
606                rate_mmPerHour_wetSnow = 0;
607                rate_mmPerHour_sleet = 0;
608
609                // Get global position
610                util::calculateXYZAtWGS84Coordinate(lat_, lon_, height_, sceneRoot, x, y, z);
611
612                // Look up every cloud layer for it's precipitation.
613                for( int i=0; i<MAX_CLOUDLAYER_SLOTS; i++ )
614                {
615                        if ( cloudLayerSlots[i].used )  // IF used, Pointer should be valid
616                        {
617                                if( cloudLayerSlots[i].cloudLayerPointer->HasPrecipitationAtPosition(x, y, z) )
618                                {
619                                        hasPrecipitation = true;
620                                        SL_MAP (int, double) precipitationMap = cloudLayerSlots[i].cloudLayerPointer->GetPrecipitation();
621                                        for( SL_MAP (int, double)::iterator it = precipitationMap.begin(); it != precipitationMap.end(); it++ )
622                                        {
623                                                switch(it->first)
624                                                {
625                                                        case SilverLining::CloudLayer::RAIN : rate_mmPerHour_rain += it->second;
626                                                                break;
627                                                        case SilverLining::CloudLayer::DRY_SNOW : rate_mmPerHour_drySnow += it->second;
628                                                                break;
629                                                        case SilverLining::CloudLayer::WET_SNOW : rate_mmPerHour_wetSnow += it->second;
630                                                                break;
631                                                        case SilverLining::CloudLayer::SLEET : rate_mmPerHour_sleet += it->second;
632                                                                break;
633                                                        default: OSG_NOTIFY( osg::FATAL ) << "ERROR: visual_skySilverLining::getOverallPrecipitationAtLocation() : Wrong precipitation type in map!" << std::endl;
634                                                                break;
635                                                };
636                                        }
637                                }       // If slot has Precipitation END
638                        }       // If used END
639                }       // For all slots END
640
641                OSG_NOTIFY( osg::ALWAYS ) << "Rain: " << rate_mmPerHour_rain << ", dry snow: " << rate_mmPerHour_drySnow << ", wet snow: " << rate_mmPerHour_wetSnow << ", sleet: " << rate_mmPerHour_sleet << std::endl;
642                return hasPrecipitation;
643
644        }       // If initialized END
645        return false;
646}
647
648
649bool visual_skySilverLining::getSlotPrecipitationAtLocation( int slot_, double& rate_mmPerHour_rain, double& rate_mmPerHour_drySnow, double& rate_mmPerHour_wetSnow, double& rate_mmPerHour_sleet, double lat_, double lon_, double height_ )
650{
651        if (isInitialized() && slot_ >= 0 && slot_ < MAX_CLOUDLAYER_SLOTS)
652        {
653                if( cloudLayerSlots[slot_].used )       // IF used, Pointer should be valid
654                {
655                        double x = 0;
656                        double y = 0;
657                        double z = 0;
658                        rate_mmPerHour_rain = 0;
659                        rate_mmPerHour_drySnow = 0;
660                        rate_mmPerHour_wetSnow = 0;
661                        rate_mmPerHour_sleet = 0;
662                        //// If -1 : Use sky internal values
663                        if ( lat_ == -1 )
664                                lat_ = lat;
665                        if ( lon_ == -1 )
666                                lon_ = lon;
667                        if ( height_ == -1 )
668                                height_ = height;
669
670                        // Get global position
671                        util::calculateXYZAtWGS84Coordinate(lat_, lon_, height_, sceneRoot, x, y, z);
672
673                        // Check for precipitation
674                        if( cloudLayerSlots[slot_].cloudLayerPointer->HasPrecipitationAtPosition(x, y, z) )
675                        {
676                                SL_MAP (int, double) precipitationMap = cloudLayerSlots[slot_].cloudLayerPointer->GetPrecipitation();
677                                for( SL_MAP (int, double)::iterator it = precipitationMap.begin(); it != precipitationMap.end(); it++ )
678                                {
679                                        switch(it->first)
680                                        {
681                                                case SilverLining::CloudLayer::RAIN : rate_mmPerHour_rain = it->second;
682                                                        break;
683                                                case SilverLining::CloudLayer::DRY_SNOW : rate_mmPerHour_drySnow = it->second;
684                                                        break;
685                                                case SilverLining::CloudLayer::WET_SNOW : rate_mmPerHour_wetSnow = it->second;
686                                                        break;
687                                                case SilverLining::CloudLayer::SLEET : rate_mmPerHour_sleet = it->second;
688                                                        break;
689                                                default: OSG_NOTIFY( osg::FATAL ) << "ERROR: visual_skySilverLining::getSlotPrecipitationAtLocation() : Wrong precipitation type in map!" << std::endl;
690                                                        break;
691                                        };
692                                }       // FOR END
693                                OSG_NOTIFY( osg::ALWAYS ) << "Rain: " << rate_mmPerHour_rain << ", dry snow: " << rate_mmPerHour_drySnow << ", wet snow: " << rate_mmPerHour_wetSnow << ", sleet: " << rate_mmPerHour_sleet << std::endl;
694                                return true;
695                        }       // If slot has Precipitation END
696                        else 
697                                return false;
698                }       // If used END
699                else
700                        return false;
701        }       // If initialized END
702        return false;
703}
704
705bool visual_skySilverLining::getSlotPrecipitation( int slot_, double& rate_mmPerHour_rain, double& rate_mmPerHour_drySnow, double& rate_mmPerHour_wetSnow, double& rate_mmPerHour_sleet )
706{
707        if (isInitialized() && slot_ >= 0 && slot_ < MAX_CLOUDLAYER_SLOTS)
708        {
709                if( cloudLayerSlots[slot_].used )       // IF used, Pointer should be valid
710                {
711                        rate_mmPerHour_rain = 0;
712                        rate_mmPerHour_drySnow = 0;
713                        rate_mmPerHour_wetSnow = 0;
714                        rate_mmPerHour_sleet = 0;
715
716                        // Check for precipitation
717
718                                SL_MAP (int, double) precipitationMap = cloudLayerSlots[slot_].cloudLayerPointer->GetPrecipitation();
719                                for( SL_MAP (int, double)::iterator it = precipitationMap.begin(); it != precipitationMap.end(); it++ )
720                                {
721                                        switch(it->first)
722                                        {
723                                                case SilverLining::CloudLayer::RAIN : rate_mmPerHour_rain = it->second;
724                                                        break;
725                                                case SilverLining::CloudLayer::DRY_SNOW : rate_mmPerHour_drySnow = it->second;
726                                                        break;
727                                                case SilverLining::CloudLayer::WET_SNOW : rate_mmPerHour_wetSnow = it->second;
728                                                        break;
729                                                case SilverLining::CloudLayer::SLEET : rate_mmPerHour_sleet = it->second;
730                                                        break;
731                                                default: OSG_NOTIFY( osg::FATAL ) << "ERROR: visual_skySilverLining::getSlotPrecipitation() : Wrong precipitation type in map!" << std::endl;
732                                                        break;
733                                        };
734                                }       // FOR END
735                               
736                                if ( rate_mmPerHour_rain>0 || rate_mmPerHour_drySnow>0 || rate_mmPerHour_wetSnow>0 || rate_mmPerHour_sleet>0)
737                                {
738                                        OSG_NOTIFY( osg::ALWAYS ) << "Rain: " << rate_mmPerHour_rain << ", dry snow: " << rate_mmPerHour_drySnow << ", wet snow: " << rate_mmPerHour_wetSnow << ", sleet: " << rate_mmPerHour_sleet << std::endl;
739                                        return true;
740                                }
741                                else
742                                        return false;
743                }       // If used END
744                else
745                        return false;
746        }       // If initialized END
747        return false;
748}
749
750void visual_skySilverLining::clearAllPrecipitation( int slot_ )
751{
752        if (isInitialized() && slot_ >= 0 && slot_ < MAX_CLOUDLAYER_SLOTS)
753        {
754                if( cloudLayerSlots[slot_].used )       // IF used, Pointer should be valid
755                {
756                        cloudLayerSlots[slot_].cloudLayerPointer->SetPrecipitation( SilverLining::CloudLayer::NONE, 0.0 );
757                        cloudLayerSlots[slot_].cloudLayerPointer->SetPrecipitation( SilverLining::CloudLayer::NONE, 0.0 ); // Second Call to clear precipitation rate.
758                }
759        }
760}
761
762void visual_skySilverLining::setSlotPrecipitation( int slot_, double rate_mmPerHour_rain_, double rate_mmPerHour_drySnow_, double rate_mmPerHour_wetSnow_, double rate_mmPerHour_sleet_ )
763{
764        if (isInitialized() && slot_ >= 0 && slot_ < MAX_CLOUDLAYER_SLOTS)
765        {
766                if( cloudLayerSlots[slot_].used )       // IF used, Pointer should be valid
767                {
768                        cloudLayerSlots[slot_].cloudLayerPointer->SetPrecipitation( SilverLining::CloudLayer::RAIN, rate_mmPerHour_rain_ );
769                        cloudLayerSlots[slot_].cloudLayerPointer->SetPrecipitation( SilverLining::CloudLayer::DRY_SNOW, rate_mmPerHour_drySnow_ );
770                        cloudLayerSlots[slot_].cloudLayerPointer->SetPrecipitation( SilverLining::CloudLayer::WET_SNOW, rate_mmPerHour_wetSnow_ );
771                        cloudLayerSlots[slot_].cloudLayerPointer->SetPrecipitation( SilverLining::CloudLayer::SLEET, rate_mmPerHour_sleet_ );
772                }
773        }
774}
Note: See TracBrowser for help on using the repository browser.