Zoom to Selected Globe Features

public void ZoomToSelectedGlobeFeatures(ESRI.ArcGIS.GlobeCore.IGlobe globe)
{
        ESRI.ArcGIS.GlobeCore.IGlobeDisplay globeDisplay = globe.GlobeDisplay;
        ESRI.ArcGIS.Analyst3D.ISceneViewer sceneViewer = globeDisplay.ActiveViewer;
        ESRI.ArcGIS.Analyst3D.ICamera camera = sceneViewer.Camera;
        ESRI.ArcGIS.GlobeCore.IGlobeCamera globeCamera = (ESRI.ArcGIS.GlobeCore.IGlobeCamera)camera; // Explicit Cast
        ESRI.ArcGIS.Analyst3D.IScene scene = globeDisplay.Scene;
        ESRI.ArcGIS.Carto.IEnumLayer enumLayer = scene.get_Layers(null, true);

        ESRI.ArcGIS.Geometry.IEnvelope envelope = new ESRI.ArcGIS.Geometry.EnvelopeClass();
        envelope.SetEmpty();
        ESRI.ArcGIS.Geometry.IEnvelope layersExtentEnvelope = new ESRI.ArcGIS.Geometry.EnvelopeClass();
        layersExtentEnvelope.SetEmpty();
        ESRI.ArcGIS.Geometry.IZAware ZAware = (ESRI.ArcGIS.Geometry.IZAware)envelope; // Explicit Cast
        ZAware.ZAware = (true);

        ESRI.ArcGIS.Geodatabase.ISpatialFilter spatialFilter = new ESRI.ArcGIS.Geodatabase.SpatialFilterClass();
        ESRI.ArcGIS.Geometry.ISpatialReference spatialReference = scene.SpatialReference;
        System.Boolean haveFeatures = false;
        enumLayer.Reset();

        ESRI.ArcGIS.Carto.ILayer layer;
        while ((layer = enumLayer.Next()) != null)
        {
                if (layer == null)
                        break;


                if (layer is ESRI.ArcGIS.Carto.IFeatureLayer)
                {
                        ESRI.ArcGIS.Carto.IFeatureLayer featureLayer = (ESRI.ArcGIS.Carto.IFeatureLayer)layer; // Explicit Cast
                        ESRI.ArcGIS.Carto.IFeatureSelection featureSelection = (ESRI.ArcGIS.Carto.IFeatureSelection)layer; // Explicit Cast
                        ESRI.ArcGIS.Geodatabase.ISelectionSet selectionSet = featureSelection.SelectionSet;
                        ESRI.ArcGIS.Geodatabase.IFeatureClass featureClass = featureLayer.FeatureClass;
                        System.String shapeField = featureClass.ShapeFieldName;
                        spatialFilter.GeometryField = shapeField;
                        spatialFilter.set_OutputSpatialReference(shapeField, spatialReference);

                        // The next 2 lines of code are different from many other ArcObjects programming techniques in that the
                        // ICursor Interface variable 'cursor' is initialized to a Null value. It is set by reference with the
                        // call to the Search method; hence the need for the 'out' argument (see MSDN for more information).
                        ESRI.ArcGIS.Geodatabase.ICursor cursor;
                        selectionSet.Search(spatialFilter, true, out cursor);

                        ESRI.ArcGIS.Geodatabase.IFeatureCursor featureCursor = (ESRI.ArcGIS.Geodatabase.IFeatureCursor)cursor; // Explicit Cast

                        System.Boolean getLayerExtent = true;
                        ESRI.ArcGIS.Geodatabase.IFeature feature;  // Automatically initialized to null. Used to test existence of a feature in the featureCursor
                        while ((feature = featureCursor.NextFeature()) != null)
                        {
                                ESRI.ArcGIS.Geometry.IGeometry geometry = feature.Shape;
                                ESRI.ArcGIS.Geometry.IEnvelope featureExtent = geometry.Envelope;
                                envelope.Union(featureExtent);
                                haveFeatures = true;

                                if (getLayerExtent)
                                {
                                        ESRI.ArcGIS.Geodatabase.IGeoDataset geoDataset = (ESRI.ArcGIS.Geodatabase.IGeoDataset)featureLayer; // Explicit Cast
                                        if (geoDataset != null)
                                        {
                                                ESRI.ArcGIS.Geometry.IEnvelope layerExtent = geoDataset.Extent;
                                                layersExtentEnvelope.Union(layerExtent);
                                        }
                                        getLayerExtent = false;
                                }
                        }
                }
        }

        // Since the size of points is very small, we use a special case to zoom in closer

        System.Double width = envelope.Width;
        System.Double height = envelope.Height;
        if (width == 0.0 && height == 0.0)  // Must be a single point, Zoom to 1 x 1 degree area,
        {                                   // or lets say 1/20th of layer extent, whichever is smallest.
                System.Double dim = 1.0;

                System.Boolean bEmpty = layersExtentEnvelope.IsEmpty;
                if (!bEmpty)
                {
                        System.Double layerWidth = layersExtentEnvelope.Width;
                        System.Double layerHeight = layersExtentEnvelope.Height;
                        System.Double layerDim = System.Math.Max(layerWidth, layerHeight) * 0.05;
                        if (layerDim > 0.0)
                                dim = System.Math.Min(1.0, layerDim);
                }

                System.Double xMin = envelope.XMin;
                System.Double yMin = envelope.YMin;

                ESRI.ArcGIS.Geometry.IPoint point = new ESRI.ArcGIS.Geometry.PointClass();
                point.X = xMin;
                point.Y = yMin;

                envelope.Width = dim;
                envelope.Height = dim;
                envelope.CenterAt(point);
        }
        else if (width == 0.0 || height == 0.0)
        {
                System.Double maxDim = System.Math.Max(width, height);
                envelope.Width = maxDim;
                envelope.Height = maxDim;
        }

        globeCamera.SetToZoomToExtents(envelope, globe, sceneViewer);
        sceneViewer.Redraw(true);
}