diff --git a/include/osg/ClusterCullingCallback b/include/osg/ClusterCullingCallback
index 9efe9ffcf..da03c7ff1 100644
--- a/include/osg/ClusterCullingCallback
+++ b/include/osg/ClusterCullingCallback
@@ -34,7 +34,7 @@ class SG_EXPORT ClusterCullingCallback : public Drawable::CullCallback
         /** compute the control point, normal and deviation from the contents of the drawable.*/
         void computeFrom(const osg::Drawable* drawable);
 
-        void set(const osg::Vec3& controlPoint, const osg::Vec3& normal, float deviation);
+        void set(const osg::Vec3& controlPoint, const osg::Vec3& normal, float deviation, float radius);
 
         void setControlPoint(const osg::Vec3& controlPoint) { _controlPoint = controlPoint; }
         const osg::Vec3& getControlPoint() const { return _controlPoint; }
diff --git a/src/osg/ClusterCullingCallback.cpp b/src/osg/ClusterCullingCallback.cpp
index 11c33372e..c719d0aef 100644
--- a/src/osg/ClusterCullingCallback.cpp
+++ b/src/osg/ClusterCullingCallback.cpp
@@ -134,17 +134,18 @@ void ClusterCullingCallback::computeFrom(const osg::Drawable* drawable)
     _radius = sqrtf(cdf._radius2);
 }
 
-void ClusterCullingCallback::set(const osg::Vec3& controlPoint, const osg::Vec3& normal, float deviation)
+void ClusterCullingCallback::set(const osg::Vec3& controlPoint, const osg::Vec3& normal, float deviation, float radius)
 {
     _controlPoint = controlPoint;
     _normal = normal;
     _deviation = deviation;
+    _radius = radius;
 }
 
 
 bool ClusterCullingCallback::cull(osg::NodeVisitor* nv, osg::Drawable* , osg::State*) const
 {
-    return false;
+//    return false;
 
     if (_deviation<=-1.0f)
     {
@@ -153,8 +154,12 @@ bool ClusterCullingCallback::cull(osg::NodeVisitor* nv, osg::Drawable* , osg::St
     }
     
     osg::Vec3 eye_cp = nv->getEyePoint() - _controlPoint;
+    float radius = eye_cp.length();
     
-    float deviation = (eye_cp * _normal)/eye_cp.length();
+    if (radius<_radius) return false;
+    
+    
+    float deviation = (eye_cp * _normal)/radius;
 
 //    osg::notify(osg::NOTICE)<<"ClusterCullingCallback::cull() _normal="<<_normal<<" _controlPointtest="<<_controlPoint<<" eye_cp="<<eye_cp<<std::endl;
 //    osg::notify(osg::NOTICE)<<"                             deviation="<<deviation<<" _deviation="<<_deviation<<" test="<<(deviation < _deviation)<<std::endl;
diff --git a/src/osgTerrain/DataSet.cpp b/src/osgTerrain/DataSet.cpp
index a2ff1c8ed..5d66d7a56 100644
--- a/src/osgTerrain/DataSet.cpp
+++ b/src/osgTerrain/DataSet.cpp
@@ -2173,6 +2173,7 @@ osg::Node* DataSet::DestinationTile::createPolygonal()
 
     float min_dot_product = 1.0f;
     float max_cluster_culling_height = 0.0f;
+    float max_cluster_culling_radius = 0.0f;
 
     for(r=0;r<numRows;++r)
     {
@@ -2215,6 +2216,7 @@ osg::Node* DataSet::DestinationTile::createPolygonal()
                     float local_m = globe_radius*( 1.0/ cos(theta+phi) - 1.0);
                     min_dot_product = osg::minimum(min_dot_product, local_dot_product);
                     max_cluster_culling_height = osg::maximum(max_cluster_culling_height,local_m);      
+                    max_cluster_culling_radius = osg::maximum(max_cluster_culling_radius,static_cast<float>(beta*globe_radius));
                 }
                 else
                 {
@@ -2305,7 +2307,8 @@ osg::Node* DataSet::DestinationTile::createPolygonal()
 
         ccc->set(center_position + transformed_center_normal*max_cluster_culling_height ,
                  transformed_center_normal, 
-                 min_dot_product);
+                 min_dot_product,
+                 max_cluster_culling_radius);
         geometry->setCullCallback(ccc);
     }