等值面代码.docxVIP

DateTimedtbg = DateTime.Now; //得到站点要插值的要素类DataTabledt = GetLastDataTable(pName);IFeatureClasspInPointFClass = ModifySationFields(dt, pName); //工作空间路径stringrasterPath = outPath + \\TempData; //要素类图层IFeatureLayerpFeatLayer = UtilityFuntion.ReturnFeatureLayer(pInPointFClass); //这个值可以通过配置文件设置stringpZValueField = string.Empty; if (this.pName == 自动站雨量) {pZValueField = YULANG; } else if (this.pName == 自动站温度) {pZValueField = WENDU; } else if (this.pName == 自动站气压) {pZValueField = QIYA; } //pZValueField = elevation;IFeatureClassDescriptorpFeatClsDes = ReturnFeatureClassDescriptor(pInPointFClass, pZValueField); //克里金插值IRasterpOutRaster = MapSpatialAnalyst.ExcuteRasterKrigie(pFeatClsDes, rasterPath, pFeatLayer); //IFeatureLayerfeatContour = null;if (pOutRaster == null) {return; }IRasterLayerpRasterLayer = new RasterLayerClass();pRasterLayer.CreateFromRaster(pOutRaster);pRasterLayer.Name = yulang;stringinputPath = pRasterLayer.FilePath;stringoutRasterPath = rasterPath + @\yulang;IWorkspaceFactorypFWSF = new ShapefileWorkspaceFactoryClass();IFeatureWorkspacepFWS = pFWSF.OpenFromFile(rasterPath, 0) as IFeatureWorkspace; // 栅格重分类 //ConversionToolspConversionTools = new ConversionTools(axMap.SpatialReference, pFWS);ConversionToolspConversionTools = new ConversionTools(axMap.SpatialReference, pFWS, rasterPath);stringmaskPolygon = outPath +@\Data\CQ.shp; //得到等值面要素类IFeatureClasspFeatureClass = pConversionTools.RasterToContourSurface(pOutRaster, maskPolygon, outRasterPath);stringfieldName = GRIDCODE; //新要素层IFeatureLayerpFeatureLayer = new FeatureLayerClass();pFeatureLayer.Name = 计算结果;pFeatureLayer.FeatureClass = pFeatureClass;MapThememapTheme = new MapTheme();UniqueValueModel model = new UniqueValueModel();model.FeatureLayer = pFeatureLayer;model.Fi