ArcGis面要素空间连接，取相交面积最大者 C#

核心代码：

         #region JoinWork_IntersectMax
         private void CreateNewFields_IntersectMax(IFeatureClass destFeatureClass,IFeatureClass sourceFeatureClass)
         {
             List<IField> fieldList = new List<IField>();
             foreach (var item in this.chkField.CheckedItems)
             {
                 IField field=sourceFeatureClass.Fields.get_Field( sourceFeatureClass.Fields.FindField(item.ToString()));
                 fieldList.Add(field);
             }
             string[] fieldArray = Utils.FieldUtils.GetFieldsNameArray(destFeatureClass);
             foreach (IField field in fieldList)
             {
                 if (!fieldArray.Contains(field.Name))
                 {
                     Utils.FieldUtils.CreateField(destFeatureClass, field.Type, field.Name, "Source_"+ field.Name, field.Length, field.Precision, field.Scale);
                 }
             }
         }
         private void GetFieldIndex(IFeatureClass destFeatureClass, IFeatureClass sourceFeatureClass, out List<int> destCLSFieldIDList, out  List<int> sourceCLSFieldIDList)
         {
             List<int> destIdList = new List<int>();
             List<int> sourceIdList = new List<int>();
             foreach (var item in this.chkField.CheckedItems)
             {
                 destIdList.Add( destFeatureClass.Fields.FindField(item.ToString()));
                 sourceIdList.Add( sourceFeatureClass.Fields.FindField(item.ToString()));
             }
             destCLSFieldIDList = destIdList;
             sourceCLSFieldIDList = sourceIdList;
         }

         private void  JoinWork_IntersectMax(IFeatureClass featureClassD, IFeatureClass featureClassS)
         {
             IFeatureCursor featureCursorD = null;
             IFeatureCursor featureCursorD_S = null;
             try
             {
                 int index = ;
                 int count = featureClassD.FeatureCount(null);
                 Utils.UserArcMapProgressbar userArcMapProgressbar = new UserArcMapProgressbar(m_hookHelper);
                 userArcMapProgressbar.SetProgressbar(,count, );

                 //create Fields
                 CreateNewFields_IntersectMax(featureClassD, featureClassS);
                 //get link field id list
                 List<int> destCLSFieldIdList;
                 List<int> sourceCLSFieldIdList;
                 GetFieldIndex(featureClassD, featureClassS, out destCLSFieldIdList, out sourceCLSFieldIdList);

                 Utils.UserWorkspace.FeatureWorkspace featureWorkspace = new Utils.UserWorkspace.FeatureWorkspace(featureClassD);
                 featureWorkspace.StartEdit(false);

                 featureCursorD = featureClassD.Update(null, true);
                 IFeature featureD = featureCursorD.NextFeature();
                 while (featureD != null)
                 {
                     userArcMapProgressbar.Position = index;
                     userArcMapProgressbar.Message = $"空间连接_相交，取重叠面积最大者_已处理：{index}/{count}";
                     userArcMapProgressbar.Step();
                     ISpatialFilter spatialFilterD_S = Utils.FilterUtils.SpatialFilter(featureD.Shape, esriSpatialRelEnum.esriSpatialRelIntersects);
                     featureCursorD_S = featureClassS.Search(spatialFilterD_S, true);
                     IFeature featureS = featureCursorD_S.NextFeature();
                     Dictionary<int, double> feature_Area = new Dictionary<int, double>();
                     while (featureS != null)
                     {
                         ITopologicalOperator topologicalOperatorD_S = featureD.Shape as ITopologicalOperator;
                         IGeometry geometryD_S = topologicalOperatorD_S.Intersect(featureS.Shape, esriGeometryDimension.esriGeometry2Dimension);
                         if (!feature_Area.ContainsKey(featureS.OID))
                             feature_Area.Add(featureS.OID, (geometryD_S as IArea).Area);
                         featureS = featureCursorD_S.NextFeature();
                     }
                     if (feature_Area.Count > )
                     {
                         IFeature featureGoal;
                         var matches = feature_Area.Where(item => item.Value == feature_Area.Values.Max());
                         featureGoal = featureClassS.GetFeature(matches.First().Key);
                         for (int i = ; i < sourceCLSFieldIdList.Count; i++)
                         {
                             featureD.set_Value(destCLSFieldIdList[i], featureGoal.get_Value(sourceCLSFieldIdList[i]));
                         }
                         featureCursorD.UpdateFeature(featureD);
                     }
                     featureD = featureCursorD.NextFeature();
                     index++;
                 }
                 userArcMapProgressbar.Hide();
                 featureWorkspace.StopEdit(true);
             }
             catch (Exception exp)
             {
                 ErrorForm err = new ErrorForm(exp.Message + "\r\n" + exp.StackTrace);
                 err.Show();
             }
             finally
             {
                 System.Runtime.InteropServices.Marshal.FinalReleaseComObject(featureCursorD);
                 System.Runtime.InteropServices.Marshal.FinalReleaseComObject(featureCursorD_S);
             }
         }
         #endregion