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 }
时间： 2024-10-23 13:31:03