ArcGIS Engine二次开发——计算shapefile面图层要素的面积

前几天，有个同事问我怎么计算面图层的面积，我也是好久没做AE的事情了，简单的查了查，告诉他用IArea接口。到了下午，他的这个问题依旧没有解决，继续求助于我。我百度了下，没有找到相应的文章，于是我意识到这个简单的事情，可能对于很多接触AE不深的人，真的是非常难得事情。最难的可能是不知道怎么做，就像我的同事一样。我很快就告诉他用IArea接口，他却惊讶的问我“你怎么知道的？”，并说他也是查到“似乎是用这个接口”，当然这是第一个层次的，属于还没入门，刚刚开始接触AE，所有的事情基本全都靠搜索引擎来解决，从网上找到代码了问题就解决了，找不到就解决不了。另一个层次的选手是知道要用某个接口，可是就是不知道怎么用，这种层次的，就算是入门了，但是还不精通，遇到问题不知道如何下手。其实无论什么层次，全部都应该学会看AE的类图和帮助，尤其是类图，在安装目录下的Diagram目录里，看多了自然就会得心应手。废话不多说了，下面是我写的一段试验代码，公布一下，希望能帮助那些有需求的童鞋，给他们节省点时间。

using System.Runtime;
using System.Runtime.InteropServices;
using ESRI.ArcGIS.Geodatabase;
using ESRI.ArcGIS.DataSourcesFile;
using ESRI.ArcGIS.Geometry;

        private void Form1_Click(object sender, EventArgs e)
        {
            IWorkspaceFactory pWSF = null;
            double dArea = 0;
            try
            {
                pWSF = new ShapefileWorkspaceFactoryClass();
                IFeatureWorkspace pWS = pWSF.OpenFromFile(@"H:\水体淹没面积", 0) as IFeatureWorkspace;
                IFeatureClass pFC = pWS.OpenFeatureClass("水体淹没面积专题.shp");

                IFeatureCursor pFeatureCur = pFC.Search(null, false);
                IFeature pFeature = pFeatureCur.NextFeature();
                while (pFeature != null)
                {
                    if (pFeature.Shape.GeometryType == esriGeometryType.esriGeometryPolygon)
                    {
                        IArea pArea = pFeature.Shape as IArea;
                        dArea = dArea + pArea.Area;
                    }

                    pFeature = pFeatureCur.NextFeature();
                }

                Marshal.ReleaseComObject(pFeatureCur);
            }
            catch (System.Exception ex)
            {
                
            }

            Marshal.ReleaseComObject(pWSF);
        }