AE代码 积累

  1         /// <summary>
  2         /// 计算两点角度
  3         /// </summary>
  4         /// <param name="p1">original point</param>
  5         /// <param name="p2">end point</param>
  6         /// <returns></returns>
  7         private float Angle(IPoint p1, IPoint p2)
  8         {
  9             double tan = Math.Atan(Math.Abs((p2.Y - p1.Y) / (p2.X - p1.X))) * 180 / Math.PI;
 10             if (p2.X > p1.X && p2.Y > p1.Y) //第一象限
 11                 return -(float.Parse(tan.ToString()));
 12             else if (p2.X > p1.X && p2.Y > p1.Y) //第二象限
 13                 return float.Parse(tan.ToString());
 14             else if (p2.X < p1.X && p2.Y > p1.Y) //第三象限
 15                 return float.Parse(tan.ToString()) - 180;
 16             else //第四象限
 17                 return 180-float.Parse(tan.ToString());                
 18         }
 19         
 20         /// <summary>
 21         /// 获取面图层的质心<IPoint>
 22         /// </summary>
 23         /// <param name="pFeatureLayer">面图层</param>
 24         /// <returns>List<></returns>
 25         private List<pointofpolygon> Get_PointFromPolygon(IFeatureLayer pFeatureLayer)
 26         {
 27             ESRI.ArcGIS.Geodatabase.IQueryFilter pQueryFilter = new ESRI.ArcGIS.Geodatabase.QueryFilterClass();
 28             pQueryFilter.WhereClause = "FID IS NOT null";
 29             IArea pArea;
 30             List<pointofpolygon> pts=new List<pointofpolygon>();
 31             //MessageBox.Show(pFeatureLayer.FeatureClass.FeatureCount(pQueryFilter).ToString());
 32             for (int i = 0; i < pFeatureLayer.FeatureClass.FeatureCount(pQueryFilter); i++)
 33             {
 34                 pArea = pFeatureLayer.FeatureClass.GetFeature(i).Shape as IArea;
 35                 pts.Insert(i, new pointofpolygon(pFeatureLayer.FeatureClass.GetFeature(i).OID, pArea.Centroid));
 36             }
 37             return pts;
 38         }
 39 
 40         /// <summary>
 41         /// 计算p1点相对p2点的角度，整栋方向为0度
 42         /// </summary>
 43         /// <param name="p1">polyline上的点</param>
 44         /// <param name="p2">original point,面的重心</param>
 45         /// <returns></returns>
 46         private float Get_Angle(IPoint p1, IPoint p2)
 47         {
 48             if (p1.X == p2.X && p1.Y == p2.Y)
 49                 return -1;
 50             IPoint p3=new PointClass();
 51             p3.Y=p2.Y;
 52             p3.X=p2.X+10;
 53             float c = (float)Math.Sqrt(Math.Pow((p1.X - p2.X), 2) + Math.Pow(p1.Y - p2.Y, 2)) * 
 54                 (float)Math.Sqrt(Math.Pow(p3.X - p2.X, 2) + Math.Pow(p3.Y - p2.Y, 2));
 55             if (c == 0) return -1;
 56             float angle = (float)Math.Acos(((p1.X - p2.X) * (p3.X - p2.X) + (p1.Y - p2.Y) * (p3.Y - p2.Y)) / c);
 57             return angle;
 58         }
 59 
 60         /// <summary>
 61         /// 清空属性表中需要计算的字段的值
 62         /// </summary>
 63         /// <param name="pFLyr">输入图层</param>
 64         private void ClearAttributeTable(IFeatureLayer pFLyr)
 65         {
 66             IDataset dataset = (IDataset)pFLyr.FeatureClass;
 67             IWorkspace workspace = dataset.Workspace;
 68             IWorkspaceEdit workspaceedit = (IWorkspaceEdit)workspace;
 69             workspaceedit.StartEditing(true);
 70             workspaceedit.StartEditOperation();
 71             for (int i = 0; i < pFLyr.FeatureClass.FeatureCount(new QueryFilter()); i++)
 72             {
 73                 IFeature pFeature = FlyrtoF(pFLyr, i);
 74                 pFeature.set_Value(pFeature.Fields.FindField("inner"),null);
 75                 pFeature.set_Value(pFeature.Fields.FindField("near"), null);
 76                 pFeature.set_Value(pFeature.Fields.FindField("far"), null);
 77                 pFeature.set_Value(pFeature.Fields.FindField("direction"), null);
 78                 pFeature.set_Value(pFeature.Fields.FindField("sum"), null);
 79                 pFeature.Store();
 80             }
 81             workspaceedit.StopEditing(true);
 82             workspaceedit.StopEditOperation();
 83         }
 84 
 85         /// <summary>
 86         /// 创建内存中的FeatureClass
 87         /// </summary>
 88         /// <param name="point"></param>
 89         public void CreateFeatureClassInMemory(IPoint point)
 90         {
 91             IWorkspaceFactory pWorkspaceFactory = new InMemoryWorkspaceFactoryClass();
 92             IWorkspaceName pWorkspaceName = pWorkspaceFactory.Create("", "pWorkspace", null, 0);
 93             IName pName = (IName)pWorkspaceName;
 94             IWorkspace pWorkspace = (IWorkspace)pName.Open();
 95             IFeatureWorkspace pFeatureWorkspace = pWorkspace as IFeatureWorkspace;
 96             IFields pFields = new FieldsClass();
 97             IFieldsEdit pFieldsEdit = pFields as IFieldsEdit;
 98             IField pField = new FieldClass();
 99             IFieldEdit pFieldEdit = pField as IFieldEdit;
100             pFieldEdit.Name_2 = "SHAPE";
101             pFieldEdit.Type_2 = esriFieldType.esriFieldTypeGeometry;
102             IGeometryDef pGeometryDef = new GeometryDefClass();
103             IGeometryDefEdit pGeometryDefEdit = pGeometryDef as IGeometryDefEdit;
104             pGeometryDefEdit.GeometryType_2 = esriGeometryType.esriGeometryPoint;
105             //为FeatureClass赋参考系，不写会出错***************************************
106             ISpatialReferenceFactory pSpatialReferenceFactory = new SpatialReferenceEnvironmentClass();
107             ISpatialReference pSpatialReference = pSpatialReferenceFactory.CreateGeographicCoordinateSystem((int)esriSRGeoCSType.esriSRGeoCS_WGS1984);
108             pGeometryDefEdit.SpatialReference_2 = pSpatialReference;
109             //************************************************************************
110             pFieldEdit.GeometryDef_2 = pGeometryDef;
111             pFieldsEdit.AddField(pField);
112             pField = new FieldClass();//不要省略写！容易出问题
113             //pFieldEdit = pField as IFieldEdit;
114             //pFieldEdit.AliasName_2 = "高程";
115             //pFieldEdit.Name_2 = "elevation";
116             //pFieldEdit.Type_2 = esriFieldType.esriFieldTypeDouble;
117             pFieldsEdit.AddField(pField);
118             IFeatureClass pFeatureClass = pFeatureWorkspace.CreateFeatureClass("1", pFields, null, null, esriFeatureType.esriFTSimple, "Shape", "");
119 
120             //插入到新建的FeatureClass中
121             IWorkspaceEdit pWorkspaceEdit = pWorkspace as IWorkspaceEdit;
122             pWorkspaceEdit.StartEditing(true);
123             pWorkspaceEdit.StartEditOperation();
124             IFeatureBuffer pFeatureBuffer = pFeatureClass.CreateFeatureBuffer();
125             IFeatureCursor pFeatureCursor = pFeatureClass.Insert(true);
126             pFeatureBuffer.Shape = point;
127             pFeatureCursor.InsertFeature(pFeatureBuffer);
128             
129             pFeatureCursor.Flush();
130             pWorkspaceEdit.StopEditOperation();
131             pWorkspaceEdit.StopEditing(true);
132 
133             pFCInMemory = pFeatureClass;
134         }
135 
136         /// <summary>
137         /// 计算pt是否在point的可视范围中
138         /// </summary>
139         /// <param name="pRasterLayer">dem</param>
140         /// <param name="point">polyline转的点</param>
141         /// <param name="pt">polygon转的点</param>
142         /// <returns>返回值为1或者0，1为可视，0为不可视</returns>
143         private bool CalViewshed(IRasterLayer pRasterLayer, IPoint point, IPoint pt)
144         {
145             //获取两个点的高程
146             IRasterSurface pRasterSurface = new RasterSurfaceClass();
147             pRasterSurface.PutRaster(pRasterLayer.Raster, 0);
148             ISurface pSurface = pRasterSurface as ISurface;
149             double elv_point = pSurface.GetElevation(point);
150             point.Z = elv_point;
151             double elv_pt = pSurface.GetElevation(pt);
152             pt.Z = elv_pt;
153             //IPolyline pPolyline = AxMapC.TrackLine() as IPolyline;
154             IPoint pPoint = null;
155             Boolean pBool = true;
156             IPolyline pVPolyline = null;
157             IPolyline pInPolyline = null;
158             object pRef = 0.13;
159             pSurface.GetLineOfSight(point, pt, out pPoint, out pVPolyline, out pInPolyline, out pBool, false, false, ref pRef);
160             //如果pt在pVPolyline上，则可视；否则不可视
161             return pBool;
162         }
163 
164         /// <summary>
165         /// 计算点所在栅格的值
166         /// </summary>
167         /// <param name="point"></param>
168         /// <param name="rasterlayer">栅格</param>
169         /// <returns>返回栅格值</returns>
170         private int GetPointInRasterValue(IPoint point, IRasterLayer rasterlayer)
171         {
172             //第一种方法
173             object obj=null;
174             IRasterProps rasterprops = (IRasterProps)rasterlayer.Raster;
175             long dHeight = rasterprops.Height;
176             long dWidth = rasterprops.Width;
177             double dx = rasterprops.MeanCellSize().X;
178             double dy = rasterprops.MeanCellSize().Y;
179             rstPixelType pixelType = rasterprops.PixelType;
180             IPnt pnt = new PntClass();
181             pnt.SetCoords(dx,dy);
182             IPixelBlock pixelBlock = rasterlayer.Raster.CreatePixelBlock(pnt);
183             IPnt pnt1 = new PntClass();
184             double xx = rasterprops.Extent.XMin;
185             double yy = rasterprops.Extent.YMax;
186             int row = (int)Math.Abs((yy - point.Y) / dx)+1;
187             int col = (int)Math.Abs((xx - point.X) / dy)+1;
188             pnt1.SetCoords(row,col);
189             rasterlayer.Raster.Read(pnt1, pixelBlock); ;
190             if (pixelBlock != null)
191             {
192                 obj = pixelBlock.GetVal(0, 0, 0);
193             }
194             if (obj != null)
195                 return (int)obj;
196             else return -1;
197             //第二种方法
198             IRasterSurface pRasterSurface=new RasterSurfaceClass();
199             pRasterSurface.PutRaster(rasterlayer.Raster,0);
200             ISurface pSurface=pRasterSurface as ISurface;
201             double elv = pSurface.GetElevation(point);
202         }
203 
204         /// <summary>
205         /// IRasterLayer类型转换成IGeoDataset
206         /// </summary>
207         /// <param name="pRasterLyr">输入DEM图层</param>
208         /// <returns>pSurfaceOP.Visibility的第一个参数</returns>
209         private IGeoDataset GetGeoDatasetFromLayer(IRasterLayer pRasterLyr)
210         {
211             IRaster pRaster = pRasterLyr.Raster;
212             IRasterBandCollection pRasterBandCollection = pRaster as IRasterBandCollection;
213             IRasterBand pRasterBand = pRasterBandCollection.Item(0);
214             IRasterDataset pRasterDataset = pRasterBand as IRasterDataset;
215             return pRasterDataset as IGeoDataset;
216         }
217 
218         /// <summary>
219         /// IgeoDataset转换成IRasterLayer
220         /// </summary>
221         /// <param name="geodataset"></param>
222         /// <returns></returns>
223         private IRasterLayer GetRasterLayerFromGeoDataset(IGeoDataset geodataset)
224         {
225             IRasterLayer pRasterLayer = new RasterLayerClass();
226             pRasterLayer.CreateFromDataset(geodataset as IRasterDataset);
227             return pRasterLayer;
228         }
229 
230         /// <summary>
231         /// FeatureLayer To FeatureClass
232         /// </summary>
233         /// <param name="pFeatureLyr"></param>
234         /// <returns></returns>
235         private IFeatureClass FLyrToFC(IFeatureLayer pFeatureLyr)
236         {
237             return pFeatureLyr.FeatureClass;
238         }
239         /// <summary>
240         /// FeatureLayer To Feature
241         /// </summary>
242         /// <param name="pFeatureLyr"></param>
243         /// <param name="index"></param>
244         /// <returns></returns>
245         private IFeature FlyrtoF(IFeatureLayer pFeatureLyr,int index)
246         {
247             IFeatureClass pFeatureClass=FLyrToFC(pFeatureLyr);
248             return pFeatureClass.GetFeature(index);
249         }
250 
251         /// <summary>
252         /// 编辑属性表
253         /// </summary>
254         /// <param name="pt_polyline">多条polyline上的节点</param>
255         /// <param name="pt_polygon">polygon的质心</param>
256         /// <param name="pFeatureLyr">polygo Layer</param>
257         private void Calulatation(List<List<IPoint>> pt_polyline,List<pointofpolygon> pt_polygon, IFeatureLayer pFeatureLyr)
258         {
259             //使图层处于编辑状态
260             IDataset dataset = (IDataset)pFeatureLyr.FeatureClass;
261             IWorkspace workspace = dataset.Workspace;
262             IWorkspaceEdit workspaceedit = (IWorkspaceEdit)workspace;
263             workspaceedit.StartEditing(true);
264             workspaceedit.StartEditOperation();
265             IFeature feature;
266             float dir = 0, buf = 0;
267             int direction=0;
268 
269             //修改某一记录某一字段的
270             feature.set_Value(feature.Fields.FindField("direction"), (float)feature.get_Value(feature.Fields.FindField("far")) + 0); 
271             feature.Store();
272             //关闭要素的编辑状态
273             workspaceedit.StopEditing(true);
274             workspaceedit.StopEditOperation();
275         }