ArcGIS Server中判断点是否在多边形内

最新推荐文章于 2024-08-08 15:18:16 发布

weixin_30507481

最新推荐文章于 2024-08-08 15:18:16 发布

阅读量188

点赞数

原文链接：http://www.cnblogs.com/tuncaysanli/archive/2009/04/18/1438717.html

版权

// 判断点pt落在哪个多边形内

ESRI.ArcGIS.ADF.Identity id = new ESRI.ArcGIS.ADF.Identity();
id.UserName = "Administrator" ;

                     id.Password = "123456";
                    id.Domain = "gis2_trq";
                    string domainstr = "10.150.38.34";

                    ESRI.ArcGIS.ADF.Connection.AGS.AGSServerConnection agsconn = new ESRI.ArcGIS.ADF.Connection.AGS.AGSServerConnection(domainstr, id);
                    agsconn.Connect();
                    if (!agsconn.IsConnected)
                    {
                        agsconn.Dispose();
                        return "";
                    }
                    ESRI.ArcGIS.Server.IServerObjectManager som = agsconn.ServerObjectManager;
                    string servertype = "MapServer";

                    string serverobjectname2 = "MapRs";
                    ESRI.ArcGIS.Server.IServerContext sc2 = som.CreateServerContext(serverobjectname2, servertype);
                    IMapServer pMapServer2 = sc2.ServerObject as IMapServer;
                    IMapServerObjects pMapServerObjs2 = pMapServer2 as IMapServerObjects;
                    IMap pMap2 = pMapServerObjs2.get_Map(pMapServer2.DefaultMapName);

                    ICompositeLayer pCLayer2 = pMap2.get_Layer(3) as ICompositeLayer;
                    IFeatureLayer pFLayer2 = pCLayer2.get_Layer(2) as IFeatureLayer;
                    IFeatureClass pFeatureClass2 = pFLayer2.FeatureClass;

ESRI.ArcGIS.Geodatabase.IQueryFilter qFilter2 = sc.CreateObject("esriGeodatabase.QueryFilter") as ESRI.ArcGIS.Geodatabase.QueryFilter;
qFilter2.WhereClause = null;

                    IFeatureCursor fCursor2 = pFeatureClass2.Search(qFilter2, true);
                    IFeature feature2;
                    string CountryName=String.Empty;
                    ESRI.ArcGIS.Geometry.IPolygon polygon = null;
                    while ((feature2 = fCursor2.NextFeature()) != null)
                    {
                        polygon = feature2.Shape as ESRI.ArcGIS.Geometry.IPolygon;
                        ESRI.ArcGIS.Geometry.IRelationalOperator pRelOperator = polygon as ESRI.ArcGIS.Geometry.IRelationalOperator;
                        if (pRelOperator.Contains(pt))
                        {
                            CountryName = feature2.get_Value(1).ToString();
                        }
                    }

//获得的CountryName 即为点pt所在的多边形的name属性

转载于:https://www.cnblogs.com/tuncaysanli/archive/2009/04/18/1438717.html