AE 创建

using System;
using System.Drawing;
using System.Runtime.InteropServices;
using ESRI.ArcGIS.ADF;
using ESRI.ArcGIS.ADF.BaseClasses;
using ESRI.ArcGIS.ADF.CATIDs;
using ESRI.ArcGIS.Controls;
using ESRI.ArcGIS.GlobeCore;
using ESRI.ArcGIS.Carto;
using ESRI.ArcGIS.Display;
using ESRI.ArcGIS.Geometry;
using ESRI.ArcGIS.esriSystem;
using ESRI.ArcGIS.Analyst3D;
using System.Windows.Forms;

using ESRI.ArcGIS.DataSourcesFile;
using ESRI.ArcGIS.DataSourcesGDB;
using ESRI.ArcGIS.Framework;
using ESRI.ArcGIS.Geodatabase;
using ESRI.ArcGIS.SystemUI;

namespace TestDrapeApplication
{
    /// <summary>
    /// Command that works in ArcGlobe or GlobeControl
    /// </summary>
    [Guid("2d4a8248-b8d1-459f-aa4d-c73fdd4cde3b")]
    [ClassInterface(ClassInterfaceType.None)]
    [ProgId("TestDrapeApplication.AddShapeFileMemory")]
    public sealed class AddShapeFileMemory : BaseCommand
    {
        #region COM Registration Function(s)
        [ComRegisterFunction()]
        [ComVisible(false)]
        static void RegisterFunction(Type registerType)
        {
            // Required for ArcGIS Component Category Registrar support
            ArcGISCategoryRegistration(registerType);

            //
            // TODO: Add any COM registration code here
            //
        }

        [ComUnregisterFunction()]
        [ComVisible(false)]
        static void UnregisterFunction(Type registerType)
        {
            // Required for ArcGIS Component Category Registrar support
            ArcGISCategoryUnregistration(registerType);

            //
            // TODO: Add any COM unregistration code here
            //
        }

        #region ArcGIS Component Category Registrar generated code
        /// <summary>
        /// Required method for ArcGIS Component Category registration -
        /// Do not modify the contents of this method with the code editor.
        /// </summary>
        private static void ArcGISCategoryRegistration(Type registerType)
        {
            string regKey = string.Format("HKEY_CLASSES_ROOT\\CLSID\\{{{0}}}", registerType.GUID);
            GMxCommands.Register(regKey);
            ControlsCommands.Register(regKey);
        }
        /// <summary>
        /// Required method for ArcGIS Component Category unregistration -
        /// Do not modify the contents of this method with the code editor.
        /// </summary>
        private static void ArcGISCategoryUnregistration(Type registerType)
        {
            string regKey = string.Format("HKEY_CLASSES_ROOT\\CLSID\\{{{0}}}", registerType.GUID);
            GMxCommands.Unregister(regKey);
            ControlsCommands.Unregister(regKey);
        }

        #endregion
        #endregion

        private IGlobeHookHelper m_globeHookHelper = null;

        public AddShapeFileMemory()
        {
            //
            // TODO: Define values for the public properties
            //
            base.m_category = "Add Shape File to map"; //localizable text
            base.m_caption = "Add Shape File";  //localizable text
            base.m_message = "Add shape File to Map";  //localizable text
            base.m_toolTip = "Add Shape File to Map";  //localizable text
            base.m_name = "AddShapeFileMemory";   //unique id, non-localizable (e.g. "MyCategory_MyCommand")

            try
            {
                //
                // TODO: change bitmap name if necessary
                //
                string bitmapResourceName = GetType().Name + ".bmp";
                base.m_bitmap = new Bitmap(GetType(), bitmapResourceName);
            }
            catch (Exception ex)
            {
                System.Diagnostics.Trace.WriteLine(ex.Message, "Invalid Bitmap");
            }
        }

        #region Overriden Class Methods

        /// <summary>
        /// Occurs when this command is created
        /// </summary>
        /// <param name="hook">Instance of the application</param>
        public override void OnCreate(object hook)
        {
            if (hook == null)
                return;

            try
            {
                m_globeHookHelper = new GlobeHookHelperClass();
                m_globeHookHelper.Hook = hook;
                if (m_globeHookHelper.ActiveViewer == null)
                {
                    m_globeHookHelper = null;
                }
            }
            catch
            {
                m_globeHookHelper = null;
            }

            if (m_globeHookHelper == null)
                base.m_enabled = false;
            else
                base.m_enabled = true;

            // TODO:  Add other initialization code
        }

        /// <summary>
        /// Occurs when this command is clicked
        /// </summary>
        public override void OnClick()
        {
            //This is run after adding document with elevation Raster layer
            IGlobeDisplay globeDisplay = new GlobeDisplayClass();

            globeDisplay = m_globeHookHelper.GlobeDisplay;

            double lon = -86.00;
            double lat = 37.9166;
            double alt = 100.0;
            double azimuth = 175.0;
            double range = 1000; //ConvertMetersToDatasetUnits((IGlobe2)globeDisplay.Globe, 1000.0);
            double fov = 60.0;

            m_globeHookHelper.GlobeDisplay.RefreshViewers();
            //Elevation present to here.
            AddElement(ref globeDisplay, azimuth, fov, range, lat, lon, alt);
            // still have elevation with polygon under elevated layer.
            m_globeHookHelper.GlobeDisplay.RefreshViewers();

        }

        #endregion

        private void AddElement(ref IGlobeDisplay globeDisplay, double azimuth, double fov, double range, double lat, double lon, double alt)
        {
            //Depending on the ASCII file type set a variable which will control the flow of the application
            string shapeType = null;
            string sensorID = "Sim.Trip1";
            //shapeType = "points";
            //    shapeType = "polylines";
            shapeType = "polygons";

            string datasetType = "Shapefile";
            IWorkspace workSpace = null;
            workSpace = this.CreateWorkspace();
            IFields fields = CreateFields(shapeType);
            //Create the Feature Class
            IFeatureClass featureClass = CreateFeatureClass(workSpace, null, shapeType, fields, null, null, "");

            AddShape(ref featureClass, ref workSpace, azimuth, fov, range, lat, lon, alt, sensorID);

            //Create a Feature Layer to display the data in ArcMap
            IFeatureLayer featureLayer = new FeatureLayerClass();
            featureLayer.FeatureClass = featureClass;
            featureLayer.Name = featureClass.AliasName;
            featureLayer.Visible = true;
            IGeoFeatureLayer lyr = (IGeoFeatureLayer)featureLayer;
            IRgbColor lyrColor = new RgbColorClass();
            lyrColor.Red = 255;
            lyrColor.Green = 0;
            lyrColor.Blue = 0;
            ISimpleFillSymbol sfSymbol = new SimpleFillSymbolClass();
            sfSymbol.Color = lyrColor;
            sfSymbol.Outline.Width = 0.4;
            sfSymbol.Style = esriSimpleFillStyle.esriSFSSolid;
            ISimpleRenderer lyrSimRend = new SimpleRendererClass();
            lyrSimRend.Symbol = (ISymbol)sfSymbol;

            lyr.Renderer = (IFeatureRenderer)lyrSimRend;

            ((IScene)globeDisplay.Globe).AddLayer(featureLayer, true);
            ILayer lyrCreated = featureLayer;

            ((IActiveView)globeDisplay.Globe).Extent = lyrCreated.AreaOfInterest;
        }
        private IWorkspace CreateWorkspace()
        {
            try
            {

                IWorkspaceFactory workspaceFactory = null;
                // Instantiate a Shapefile workspace factory
                //workspaceFactory = new ShapefileWorkspaceFactoryClass();
                Type factoryType = Type.GetTypeFromProgID("esriDataSourcesFile.ShapefileWorkspaceFactory");
                workspaceFactory = (IWorkspaceFactory)Activator.CreateInstance(factoryType);

                //Create a directory hierarchy to seperate out datasets created for Points, Polylines, and Polygons

                IWorkspaceName workspaceName = workspaceFactory.Create(Application.StartupPath + @"\ShapeFile", "Shapefile", null, 0);
                IName Name = (IName)workspaceName;
                IWorkspace workspace = (IWorkspace)(Name.Open());
                return workspace;
            }
            catch (Exception ex)
            {
                String mess = ex.Message;
                MessageBox.Show(mess);
                return null;
            }

        }

        private IFields CreateFields(string shapeType)
        {

            // Create a new fields collection.
            IFields fields = new FieldsClass();

            // Cast to IFieldsEdit to modify the properties of the fields collection.
            IFieldsEdit fieldsEdit = (IFieldsEdit)fields;

            // Set the number of fields the collection will contain.
            fieldsEdit.FieldCount_2 = 5;

            // Create the ObjectID field.
            IField oidField = new FieldClass();
            IFieldEdit oidFieldEdit = (IFieldEdit)oidField;
            oidFieldEdit.Name_2 = "ObjectID";
            oidFieldEdit.AliasName_2 = "FID";
            oidFieldEdit.Type_2 = esriFieldType.esriFieldTypeOID;
            fieldsEdit.set_Field(0, oidField);

            // Create the myID field.
            IField myIDField = new FieldClass();
            IFieldEdit myIDFieldEdit = (IFieldEdit)myIDField;
            myIDFieldEdit.Name_2 = "SensorID";
            myIDFieldEdit.Type_2 = esriFieldType.esriFieldTypeString;
            fieldsEdit.set_Field(1, myIDField);

            // Create the text field.
            IField textField = new FieldClass();
            IFieldEdit textFieldEdit = (IFieldEdit)textField;
            textFieldEdit.Length_2 = 30; // Only string fields require that you set the length.
            textFieldEdit.Name_2 = "Str_Att1";
            textFieldEdit.Type_2 = esriFieldType.esriFieldTypeString;
            fieldsEdit.set_Field(2, textField);

            // Create the numeric field.
            IField numericField = new FieldClass();
            IFieldEdit numericFieldEdit = (IFieldEdit)numericField;
            numericFieldEdit.Name_2 = "Num_Att1";
            numericFieldEdit.Type_2 = esriFieldType.esriFieldTypeDouble;
            fieldsEdit.set_Field(3, numericField);

            //Create the SHAPE field.
            IField shapeField = new FieldClass();
            IFieldEdit shapeFieldEdit = (IFieldEdit)shapeField;
            shapeFieldEdit.Name_2 = "SHAPE";
            shapeFieldEdit.IsNullable_2 = true;
            shapeFieldEdit.Required_2 = true;

            IGeometryDef geometryDef = new GeometryDefClass();
            IGeometryDefEdit geometryDefEdit = (IGeometryDefEdit)geometryDef;
            geometryDefEdit.AvgNumPoints_2 = 1;
            geometryDefEdit.GridCount_2 = 0;
            geometryDefEdit.HasM_2 = false;
            geometryDefEdit.HasZ_2 = false;

            //Depending on the feature type (ie. Point, Polyline, or Polygon)
            //set the appropriate .Type_2 and .GeometryType_2 Properties.
            if (shapeType == "points")
            {
                shapeFieldEdit.Type_2 = esriFieldType.esriFieldTypeGeometry;
                geometryDefEdit.GeometryType_2 = esriGeometryType.esriGeometryPoint;
            }
            else if (shapeType == "polylines")
            {
                shapeFieldEdit.Type_2 = esriFieldType.esriFieldTypeGeometry;
                geometryDefEdit.GeometryType_2 = esriGeometryType.esriGeometryPolyline;
            }
            else if (shapeType == "polygons")
            {
                shapeFieldEdit.Type_2 = esriFieldType.esriFieldTypeGeometry;
                geometryDefEdit.GeometryType_2 = esriGeometryType.esriGeometryPolygon;
            }

            //Generate a default Spatial Reference
            ISpatialReference spatialReference = MakeSpatialReference(esriSRProjCSType.esriSRProjCS_Sphere_AzimuthalEquidistant);
            geometryDefEdit.SpatialReference_2 = spatialReference;

            shapeFieldEdit.GeometryDef_2 = (GeometryDef)geometryDefEdit;
            fieldsEdit.set_Field(4, shapeFieldEdit);

            return fields;

        }

        ///<summary>Helper to create a Feature Class.</summary>
        public IFeatureClass CreateFeatureClass(IWorkspace workspace, IFeatureDataset featureDataset, System.String featureClassName, IFields fields, ESRI.ArcGIS.esriSystem.UID CLSID, ESRI.ArcGIS.esriSystem.UID CLSEXT, System.String strConfigKeyword)
        {

            if (featureClassName == "")
            {
                return null; // name was not passed in
            }
            featureClassName = "newShapeFile";
            IFeatureClass featureClass = null;
            IFeatureWorkspace featureWorkspace = (IFeatureWorkspace)workspace; // Explicit Cast

            // assign the class id value if not assigned
            if (CLSID == null)
            {
                CLSID = new ESRI.ArcGIS.esriSystem.UIDClass();
                CLSID.Value = "esriGeoDatabase.Feature";
            }

            System.String strShapeField = "";

            // locate the shape field
            //for (Int32 j = 0; j < fields.FieldCount; j++)
            //{
            //    if (fields.get_Field(j).Type == esriFieldType.esriFieldTypeGeometry)
            //    {
            //        strShapeField = fields.get_Field(j).Name;
            //    }
            //}
            strShapeField = "SHAPE";
            // finally create and return the feature class

            // no feature dataset passed in, create at the workspace level
            try
            {
                featureClass = featureWorkspace.CreateFeatureClass(featureClassName, fields, CLSID, CLSEXT, esriFeatureType.esriFTSimple, strShapeField, strConfigKeyword);
            }
            catch (Exception ex)
            {
                String mess = ex.Message;
                MessageBox.Show(mess);
            }
                return featureClass;

        }

        public ISpatialReference MakeSpatialReference(esriSRProjCSType coordinateSystem)
        {

            ISpatialReferenceFactory spatialReferenceFactory = new SpatialReferenceEnvironmentClass();

            //Create a projected coordinate system and define its domain, resolution, and x,y tolerance.
            ISpatialReferenceResolution spatialReferenceResolution = spatialReferenceFactory.CreateProjectedCoordinateSystem(System.Convert.ToInt32(coordinateSystem)) as ISpatialReferenceResolution;
            spatialReferenceResolution.ConstructFromHorizon();
            ISpatialReferenceTolerance spatialReferenceTolerance = spatialReferenceResolution as ISpatialReferenceTolerance;
            spatialReferenceTolerance.SetDefaultXYTolerance();
            ISpatialReference spatialReference = spatialReferenceResolution as ISpatialReference;

            return spatialReference;

        }

        private void AddShape(ref IFeatureClass FeatureClass, ref IWorkspace workspace, double azimuth, double fov, double range, double lat, double lon, double alt, String sensorID)
        {
            // Start an edit session and edit operation.
            IWorkspaceEdit workspaceEdit = (IWorkspaceEdit)workspace;
            workspaceEdit.StartEditing(true);
            workspaceEdit.StartEditOperation();

            int IDAttributeFieldIndex = FeatureClass.FindField("SensorID");
            int StringAttributeFieldIndex = FeatureClass.FindField("Str_Att1");
            int NumericAttributeFieldIndex = FeatureClass.FindField("Num_Att1");

            ComReleaser comReleaser = new ComReleaser();
            // Create a new insert cursor with buffering.
            IFeatureCursor featureCursor = FeatureClass.Insert(true);
            comReleaser.ManageLifetime(featureCursor);

            // Create a feature buffer. This will store the values common to every
            // feature to be installed.
            IFeatureBuffer featureBuffer = FeatureClass.CreateFeatureBuffer();
            comReleaser.ManageLifetime(featureBuffer);

            IPolygon polygon = makePolygon(azimuth, fov, range, lat, lon, alt);
            featureBuffer.Shape = (IGeometry)polygon;
            featureBuffer.set_Value(IDAttributeFieldIndex, sensorID);
            featureBuffer.set_Value(StringAttributeFieldIndex, "Try this");
            featureBuffer.set_Value(IDAttributeFieldIndex, 2);
            featureCursor.InsertFeature(featureBuffer);

            // Flush the cursor.
            featureCursor.Flush();

            // All of the features were successfully inserted; stop the edit operation
            // and stop the edit session, saving the changes made in edit operations.
            workspaceEdit.StopEditOperation();
            workspaceEdit.StopEditing(true);
        }
        #region "Build Polygon"

        /// <summary>
        /// Takes input of azimuth field of view and range with position and altitude
        /// to make a circular arc polygon to add as a symbol for element in layer.
        /// </summary>
        /// <param name="azimuth"></param>
        /// <param name="fov"></param>
        /// <param name="range"></param>
        /// <param name="lat"></param>
        /// <param name="lon"></param>
        /// <param name="alt"></param>
        /// <returns>IPolygon polygon</returns>
        private IPolygon makePolygon(double azimuth, double fov, double range, double lat, double lon, double alt)
        {
            object missing = System.Reflection.Missing.Value;
            //IZAware zAware;
            IPolygon polygon = new PolygonClass();

            IPoint centerPoint = new PointClass();
            IPoint startPoint = new PointClass();
            IPoint endPoint = new PointClass();
            Boolean isCCW = true;
            Boolean isMinor = false;
            centerPoint.SpatialReference = MakeSpatialReference(esriSRProjCSType.esriSRProjCS_Sphere_AzimuthalEquidistant);
            startPoint.SpatialReference = centerPoint.SpatialReference;
            endPoint.SpatialReference = centerPoint.SpatialReference;
            //zAware = centerPoint as IZAware;
            ////zAware.ZAware = true;

            //zAware = startPoint as IZAware;
            ////zAware.ZAware = true;

            //zAware = endPoint as IZAware;
            ////zAware.ZAware = true;

            centerPoint.PutCoords(lon, lat);
            ////centerPoint.Z = alt;
            double Xray = centerPoint.X;
            double Yankee = centerPoint.Y;

            ICircularArc circularArc = new CircularArcClass();

            ISegmentCollection ring1 = new RingClass();
            ring1.AddSegment(circularArc as ISegment, ref missing, ref missing);
            double fovRadian = 0.0;
            fovRadian = (fov * 2 * Math.PI) / 360;
            circularArc.PutCoordsByAngle(centerPoint, GetFromAngle(azimuth, fov), fovRadian, range);

            circularArc.QueryCoords(centerPoint, startPoint, endPoint, ref isCCW, ref isMinor);

            ILine lineOut = new LineClass();
            ILine lineIn = new LineClass();
            lineOut.PutCoords(centerPoint, startPoint);
            lineIn.PutCoords(endPoint, centerPoint);

            ISegmentCollection segColl = new PolygonClass();
            segColl.AddSegment((ISegment)lineOut, ref missing, ref missing);
            segColl.AddSegment((ISegment)circularArc, ref missing, ref missing);
            segColl.AddSegment((ISegment)lineIn, ref missing, ref missing);
            polygon = segColl as IPolygon;
            //IGeometryCollection geometryCollection = polygon as IGeometryCollection;
            //geometryCollection.AddGeometry(segColl as IGeometry, ref missing, ref missing);

            //geometryCollection.GeometriesChanged();
            return polygon;

        }

        /// <summary>
        /// Converts Azimuth and Field of View width to geometric
        /// start angle of circular arc.
        /// </summary>
        /// <param name="azimuth"></param>
        /// <param name="FOV"></param>
        /// <returns>Double startAngle</returns>
        private double GetFromAngle(double azimuth, double FOV)
        {
            double startAngle = 0.0;
            double fovRadian = 0.0;
            fovRadian = (FOV * 2 * Math.PI) / 360;
            startAngle = AzimuthToAngle(azimuth) - (.5 * fovRadian);
            return startAngle;
        }
        /// <summary>
        /// Converts Azimuth to geometric angle of arc.
        /// </summary>
        /// <param name="azimuth"></param>
        /// <returns>double angle</returns>
        private double AzimuthToAngle(double azimuth)
        {
            double angle = 0.0;
            double angleRadian = 0.0;

            angleRadian = (azimuth * 2 * Math.PI) / 360;
            angle = (2 * Math.PI) - (angleRadian - Math.PI / 2);
            if (angle >= 2 * Math.PI)
            {
                angle -= 2 * Math.PI;
            }
            return angle;
        }

        private double ConvertMetersToDatasetUnits(IGlobe2 globe, double meters)
        {
            try
            {
                esriUnits units = globe.GlobeUnits;
                IUnitConverter cvtr = new UnitConverter();
                double mapUnits = cvtr.ConvertUnits(meters, esriUnits.esriMeters, units);
                return mapUnits;
            }
            catch (Exception ex)
            {
                String mess = ex.Message;
                return -1.0;
            }

        }

        #endregion

    }
}