AE一些代码

来自：http://www.pudn.com/downloads153/sourcecode/graph/detail672576.html

using System; 

using System.Drawing; 

using System.Collections; 

using System.ComponentModel; 

using System.Windows.Forms; 

using System.Data; 

using System.IO; 

using System.Runtime.InteropServices; 

 

using ESRI.ArcGIS.esriSystem; 

using ESRI.ArcGIS.Carto; 

using ESRI.ArcGIS.Controls; 

using ESRI.ArcGIS.ADF; 

using ESRI.ArcGIS.SystemUI; 

using ESRI.ArcGIS.Geodatabase; 

using ESRI.ArcGIS.Display; 

using ESRI.ArcGIS.GeoAnalyst; 

 

using ESRI.ArcGIS.DataSourcesRaster; 

using ESRI.ArcGIS.Geometry; 

using ESRI.ArcGIS.Output;//输出栅格 

namespace RasterInterpolate 

{ 

    public sealed partial class MainForm : Form 

    { 

        #region class private members 

        private IMapControl3 m_mapControl = null; 

        private string m_mapDocumentName = string.Empty; 

        #endregion 

 

        #region class constructor 

        public MainForm() 

        { 

            InitializeComponent(); 

        } 

        #endregion 

        string[] PtName; 

        //double[] Cx; 

        //double[] Cy; 

        double[] Ptx; 

        double[] Pty; 

        double[] Railfall; 

        double[] Ptz; 

        double[] lat; 

        double[] lon; 

        double[] di;//大圆距离 

        private void MainForm_Load(object sender, EventArgs e) 

        { 

            //get the MapControl 

            m_mapControl = (IMapControl3)axMapControl1.Object; 

            m_mapControl.AddShapeFile(Application.StartupPath + @"\shp", "sheet2.shp"); 

           // m_mapControl.LoadMxFile(Application.StartupPath + @"\data\站点.mxd","",""); 

            //disable the Save menu (since there is no document yet) 

            menuSaveDoc.Enabled = false; 

 

             dt=oledata.operate(); 

             PtName = new string[dt.Rows.Count]; 

             //Cx = new double[dt.Rows.Count]; 

             //Cy = new double[dt.Rows.Count]; 

             Ptx = new double[dt.Rows.Count]; 

             Pty = new double[dt.Rows.Count]; 

             Railfall = new double[dt.Rows.Count]; 

             Ptz = new double[dt.Rows.Count]; 

             lon = new double[dt.Rows.Count]; 

             lat=new double[dt.Rows.Count]; 

             di = new double[dt.Rows.Count]; 

             Init();//数据赋值 

        } 

        //private int m_ColumnCount; 

        //private int m_RowCount; 

        //private int m_size; 

        #region Main Menu event handlers 

        private void menuNewDoc_Click(object sender, EventArgs e) 

        { 

            //execute New Document command 

            ICommand command = new CreateNewDocument(); 

            command.OnCreate(m_mapControl.Object); 

            command.OnClick(); 

        } 

 

        private void menuOpenDoc_Click(object sender, EventArgs e) 

        { 

            //execute Open Document command 

            ICommand command = new ControlsOpenDocCommandClass(); 

            command.OnCreate(m_mapControl.Object); 

            command.OnClick(); 

        } 

 

        private void menuSaveDoc_Click(object sender, EventArgs e) 

        { 

            //execute Save Document command 

            if (m_mapControl.CheckMxFile(m_mapDocumentName)) 

            { 

                //create a new instance of a MapDocument 

                IMapDocument mapDoc = new MapDocumentClass(); 

                mapDoc.Open(m_mapDocumentName, string.Empty); 

 

                //Make sure that the MapDocument is not readonly 

                if (mapDoc.get_IsReadOnly(m_mapDocumentName)) 

                { 

                    MessageBox.Show("Map document is read only!"); 

                    mapDoc.Close(); 

                    return; 

                } 

 

                //Replace its contents with the current map 

                mapDoc.ReplaceContents((IMxdContents)m_mapControl.Map); 

 

                //save the MapDocument in order to persist it 

                mapDoc.Save(mapDoc.UsesRelativePaths, false); 

 

                //close the MapDocument 

                mapDoc.Close(); 

            } 

        } 

 

        private void menuSaveAs_Click(object sender, EventArgs e) 

        { 

            //execute SaveAs Document command 

            ICommand command = new ControlsSaveAsDocCommandClass(); 

            command.OnCreate(m_mapControl.Object); 

            command.OnClick(); 

        } 

 

        private void menuExitApp_Click(object sender, EventArgs e) 

        { 

            //exit the application 

            Application.Exit(); 

        } 

        #endregion 

 

        //listen to MapReplaced evant in order to update the statusbar and the Save menu 

        private void axMapControl1_OnMapReplaced(object sender, IMapControlEvents2_OnMapReplacedEvent e) 

        { 

            //get the current document name from the MapControl 

            m_mapDocumentName = m_mapControl.DocumentFilename; 

 

            //if there is no MapDocument, diable the Save menu and clear the statusbar 

            if (m_mapDocumentName == string.Empty) 

            { 

                menuSaveDoc.Enabled = false; 

                statusBarXY.Text = string.Empty; 

            } 

            else 

            { 

                //enable the Save manu and write the doc name to the statusbar 

                menuSaveDoc.Enabled = true; 

                statusBarXY.Text = System.IO.Path.GetFileName(m_mapDocumentName); 

            } 

        } 

 

        private void axMapControl1_OnMouseMove(object sender, IMapControlEvents2_OnMouseMoveEvent e) 

        { 

            statusBarXY.Text = string.Format("{0}, {1}  {2}", e.mapX.ToString("#######.##"), e.mapY.ToString("#######.##"), axMapControl1.MapUnits.ToString().Substring(4)); 

        } 

        public void CreateContourListFromShapefile()//内插生成高程信息 

        { 

            //用反距离IDW插值生成的栅格图像 

            IInterpolationOp pInterpolationOp; 

            pInterpolationOp = new RasterInterpolationOpClass(); 

            // Create the input point object 

 

            IFeatureClass pFeatureClass; 

            IFeatureLayer pFeaturelayer; 

            pFeaturelayer = this.m_mapControl.Map.get_Layer(0) as IFeatureLayer;//这里我的数据在第八层 

            pFeatureClass = pFeaturelayer.FeatureClass; 

 

            // Define the search radius 

            IRasterRadius pRadius; 

            pRadius = new RasterRadiusClass(); 

            object maxDistance = Type.Missing; 

            pRadius.SetVariable(12, ref maxDistance); 

 

            //Create FeatureClassDescriptor using a value field 

            IFeatureClassDescriptor pFCDescriptor; 

            pFCDescriptor = new FeatureClassDescriptorClass(); 

            pFCDescriptor.Create(pFeatureClass, null, "高程（米）"); 

 

            //set { cellsize for output raster in the environment 

            object dCellSize = cellSize; 

            IRasterAnalysisEnvironment pEnv; 

            pEnv = new RasterAnalysis(); 

            pEnv = pInterpolationOp as IRasterAnalysisEnvironment; 

            pEnv.SetCellSize(esriRasterEnvSettingEnum.esriRasterEnvValue, ref dCellSize); 

 

            object objectbarrier = Type.Missing; 

            //Perform the interpolation 

            IGeoDataset rasDataset; 

            rasDataset = pInterpolationOp.IDW((IGeoDataset)pFCDescriptor, 2, pRadius, ref objectbarrier); 

            //IRasterDataset rsDt;//这里就是调用Esri函数的地方 克里金 双线性 那些都可以的 

            //rsDt.CreateDefaultRaster(); 

            IRaster pOutRaster; 

            pOutRaster = rasDataset as IRaster; 

             

            //Add output into ArcMap as a raster layer 

            IRasterLayer pOutRasLayer; 

            

            pOutRasLayer = new RasterLayerClass(); 

            pOutRasLayer.CreateFromRaster(pOutRaster); 

            pOutRasLayer.Name = "Contour"; 

            //m_RowCount = pOutRasLayer.RowCount; 

            //m_ColumnCount = pOutRasLayer.ColumnCount; 

 

            this.m_mapControl.AddLayer(pOutRasLayer, 0);// 

 

            UsingRasterClassifyColorRampRenderer(); 

 

    

            IRawPixels rawPixels = null; 

            IRasterBandCollection rasterBandCollection; 

            IRasterProps rasterProps; 

            // QI for IRawPixels and IRasterProps     

            rasterBandCollection = (IRasterBandCollection)rasDataset; 

            rawPixels = (IRawPixels)rasterBandCollection.Item(0); 

            rasterProps = (IRasterProps)rawPixels; 

          //  int dHeight = rasterProps.Height;//当前栅格数据集的行数 

          //  int dWidth = rasterProps.Width; //当前栅格数据集的列数 

           // double dX = rasterProps.MeanCellSize().X; //栅格的宽度 

           // double dY = rasterProps.MeanCellSize().Y; //栅格的高度 

            IEnvelope extent = rasterProps.Extent; //当前栅格数据集的范围 

            //IPoint pt=new  PointClass(); 

            IPnt pixelBlockSize = null; 

            IPnt pixelBlockOrigin = null; 

            // Create pixelblock     

            pixelBlockOrigin = new DblPntClass(); 

            pixelBlockOrigin.SetCoords(0, 0); 

            IPixelBlock3 pixelBlock3; 

            pixelBlockSize = new DblPntClass(); 

            pixelBlockSize.SetCoords(rasterProps.Width, rasterProps.Height); 

            pixelBlock3 = rawPixels.CreatePixelBlock(pixelBlockSize) as IPixelBlock3; 

            // Read pixelblock     

            rawPixels.Read(pixelBlockOrigin, (IPixelBlock)pixelBlock3); 

            // Get pixeldata array      

            System.Array pixelData; 

            pixelData = (System.Array)pixelBlock3.get_PixelDataByRef(0); 

           // IRasterDataset RsDateset; 

            createFileRasterDataset("c:\\temp", RainfallValue, extent.LowerLeft, rasterProps.Width, rasterProps.Height, pixelData); 

 

            //IRasterLayer pRasterLayer = new RasterLayerClass(); 

            //pRasterLayer.CreateFromDataset(RsDateset); 

            //ILayer pLayer; 

            //pLayer = pRasterLayer; 

            //pLayer.Name = "Result"; 

            //axMapControl1.Map.AddLayer(pLayer); 

            //axMapControl1.ActiveView.Refresh(); 

             

            UsingRasterClassifyColorRampRenderer();//这里是要创建刚另一个图层（结果） 在c：\temp 

            //IWorkspaceFactory pWSF;//保存栅格 

            //IWorkspace pWS; 

            //pWSF = new RasterWorkspaceFactory(); 

            //pWS = pWSF.OpenFromFile("c:\\temp", 0); 

            //IRaster raster; 

            //raster = RsDateset.CreateDefaultRaster(); 

            //ISaveAs isa; 

            //isa = raster as ISaveAs; 

            //isa.SaveAs("test2.tif",pWS, "TIFF"); 

        } 

        public int cellSize; 

        public int classifyNum; 

        public string RainfallValue; 

        private void gIDS法ToolStripMenuItem_Click(object sender, EventArgs e) 

        { 

            FrmSet frmSet = new FrmSet(); 

             

            if(DialogResult.OK==frmSet.ShowDialog()) 

            { 

                classifyNum = frmSet.classifyNum; 

                cellSize = frmSet.cellSize; 

                RainfallValue=frmSet.RailfallValue; 

                if (RainfallValue=="降雨量50") 

                { 

                    Cx = 0.0006; Cy = 0.0005; Cz = 0.4175; 

                } 

                else if (RainfallValue == "降雨量100") 

                { 

                    Cx = 0.0007; 

                    Cy = 0.0006; 

                    Cz = 0.5199; 

                } 

                else if (RainfallValue == "降雨量500") 

                { 

                    Cx = 0.001; 

                    Cy = 0.0009; 

                    Cz = 0.769; 

                } 

                else if (RainfallValue == "降雨量1000") 

                { 

                    Cx = 0.0011; 

                    Cy = 0.001; 

                    Cz = 0.8789; 

                } 

                ResetRailfall(); 

                CreateContourListFromShapefile(); 

            } 

 

        } 

        

        

 

        public static IColor ConvertColorToIColor(Color color) 

        { 

 

            IColor pColor = new RgbColorClass(); 

            pColor.RGB = (int)color.R | (int)(color.G) < 8 | (int)(color.B) < 16; 

            return pColor; 

        }  

    

        private void UsingRasterClassifyColorRampRenderer() //分级渲染  

        { 

            //get { Map 

             

            IMap pMap; 

            pMap = m_mapControl.Map; 

            //get { raster input from layer 

            IRasterLayer pRLayer; 

            pRLayer = (IRasterLayer)pMap.get_Layer(0); 

            IRaster pRaster; 

            pRaster = pRLayer.Raster; 

 

            //Create classfy renderer and QI RasterRenderer interface 

            IRasterClassifyColorRampRenderer pClassRen; 

            pClassRen = new RasterClassifyColorRampRenderer(); 

            IRasterRenderer pRasRen; 

            pRasRen = pClassRen as IRasterRenderer; 

 

            //set { raster for the render and update 

            pRasRen.Raster = pRaster; 

            pClassRen.ClassCount = classifyNum; 

            pRasRen.Update(); 

 

            //Create a color ramp to use 

            IAlgorithmicColorRamp pRamp; 

            pRamp = new AlgorithmicColorRamp(); 

            Color clr = new Color(); 

            clr = Color.FromArgb(255, 255, 128); 

            pRamp.FromColor = ConvertColorToIColor(clr); 

            clr = Color.FromArgb(135, 42, 8); 

            pRamp.ToColor = ConvertColorToIColor(clr); 

            pRamp.Size = 10; 

            bool flag = true; 

            pRamp.CreateRamp(out flag); 

 

            //Create symbol for the classes 

            IFillSymbol pFSymbol; 

            pFSymbol = new SimpleFillSymbol(); 

 

            //loop through the classes and apply the color and label 

            int I; 

            for (I = 0; I = pClassRen.ClassCount - 1; I++) 

            { 

                pFSymbol.Color = pRamp.get_Color(I);//  .Color(I); 

                pClassRen.set_Symbol(I, pFSymbol as ISymbol); 

                pClassRen.set_Label(I, "等级" + Convert.ToString(I + 1)); 

            } // I 

 

            //Update the renderer and plug into layer 

            pRasRen.Update(); 

            pRLayer.Renderer = (IRasterRenderer)pClassRen; 

            m_mapControl.ActiveView.Refresh(); 

            // m_mapControl; 

            this.axTOCControl1.Update(); 

 

            //Release memeory 

 

            pMap = null; 

            pRLayer = null; 

            pRaster = null; 

            pRasRen = null; 

            pClassRen = null; 

            pRamp = null; 

            pFSymbol = null; 

        } 

 

        private void rasterDataToolStripMenuItem_Click(object sender, EventArgs e) 

        { 

            ExportImage(this.m_mapControl.ActiveView); 

        } 

        private void ExportImage(IActiveView pActiveView) 

        { 

            SaveFileDialog pSaveFileDialog = new SaveFileDialog(); 

            pSaveFileDialog.Filter = "JPEG(*.jpg)|*.jpg|AI(*.ai)|*.ai|BMP(*.BMP)|*.bmp|EMF(*.emf)|*.emf|GIF(*.gif)|*.gif|PDF(*.pdf)|*.pdf|PNG(*.png)|*.png|EPS(*.eps)|*.eps|SVG(*.svg)|*.svg|TIFF(*.tif)|*.tif"; 

            pSaveFileDialog.Title = "输出地图"; 

            pSaveFileDialog.RestoreDirectory = true; 

            pSaveFileDialog.FilterIndex = 1; 

 

            if (pSaveFileDialog.ShowDialog() != System.Windows.Forms.DialogResult.OK) 

            { 

                return; 

            } 

            string sFilePath = pSaveFileDialog.FileName; 

            string sPathName = System.IO.Path.GetDirectoryName(sFilePath); 

            string sFileName = System.IO.Path.GetFileNameWithoutExtension(sFilePath); 

            IExport pExporter = null; 

 

            switch (pSaveFileDialog.FilterIndex) 

            { 

 

                case 1: 

                    pExporter = new ExportJPEGClass(); 

                    break; 

                case 2: 

                    pExporter = new ExportBMPClass(); 

                   break; 

                case 3: 

                    pExporter = new ExportEMFClass(); 

                    break; 

                case 4: 

                    pExporter = new ExportGIFClass(); 

                    break; 

                case 5: 

 

                    pExporter = new ExportAIClass(); 

                    break; 

                case 6: 

                    pExporter = new ExportPDFClass(); 

                    break; 

                case 7: 

                    pExporter = new ExportPNGClass(); 

                    break; 

                case 8: 

                    pExporter = new ExportPSClass(); 

                    break; 

                case 9: 

                    pExporter = new ExportSVGClass(); 

                    break; 

                case 10: 

                    pExporter = new ExportTIFFClass(); 

                    break; 

                default: 

                    MessageBox.Show("输出格式错误"); 

                    return; 

            } 

            IEnvelope pEnvelope = new EnvelopeClass(); 

            ITrackCancel pTrackCancel = new CancelTrackerClass(); 

            tagRECT ptagRECT; 

            ptagRECT.left = 0; 

            ptagRECT.top = 0; 

            ptagRECT.right = (int)pActiveView.Extent.Width; 

            ptagRECT.bottom = (int)pActiveView.Extent.Height; 

            int pResolution = (int)(pActiveView.ScreenDisplay.DisplayTransformation.Resolution); 

            pEnvelope.PutCoords(ptagRECT.left, ptagRECT.bottom, ptagRECT.right, ptagRECT.top); 

 

            pExporter.Resolution = pResolution; 

            pExporter.ExportFileName = sFilePath; 

            pExporter.PixelBounds = pEnvelope; 

            pActiveView.Output(pExporter.StartExporting(), pResolution, ref ptagRECT, pActiveView.Extent, pTrackCancel); 

            pExporter.FinishExporting(); 

        } 

       // IRasterDataset 

        private void createFileRasterDataset(string directoryName, string fileName, IPoint originPoint, int rowCount, int columCount, System.Array CountourpixelVale)//创建输出结果图层 

        { 

            //The file must not be there. Otherwise, dataset can not be got. 

            if (File.Exists(directoryName + fileName)) 

            { 

                MessageBox.Show("File Exist!"); 

                IRasterLayer pRasterLayer = new RasterLayerClass(); 

                pRasterLayer.CreateFromFilePath(directoryName + fileName); 

                ILayer pLayer; 

                pLayer = pRasterLayer; 

                pLayer.Name = "Result"; 

                axMapControl1.Map.AddLayer(pLayer); 

                axMapControl1.ActiveView.Refresh(); 

                return;// null; 

            } 

            System.IO.FileStream fs = File.Create("C:\\temp\\" + RainfallValue + "均方差.txt"); 

            StreamWriter sw = new StreamWriter((System.IO.Stream)fs); 

            sw.WriteLine("站点名，x坐标，y坐标," + RainfallValue + "均方差"); 

            // This function creates a new img file in the given workspace   

            // and then assigns pixel values 

            try 

            { 

                IRasterDataset rasterDataset = null; 

               // double cellSize = 500.0;//栅格大小xugy 

              //  IPoint originPoint = new PointClass(); 

                double minX = originPoint.X - cellSize / 2; 

                double minY = originPoint.Y - cellSize / 2; 

                //double minX=GetMinX(); 

                //double minY=GetMinY(); 

                //double maxX=GetMaxX(); 

                //double maxY=GetMaxY(); 

                double minLon = GetMinLon(); 

                double minLat = GetMinLat(); 

                double maxLon = GetMaxLon(); 

                double maxLat = GetMaxLat(); 

               // originPoint.PutCoords(minX - cellSize / 2, minY - cellSize / 2); 

                // Create the dataset    

                IRasterWorkspace2 rasterWorkspace2 = null; 

                // IGeographicCoordinateSystem m_GeographicCoordinateSystem;  

                ISpatialReference SpRef = new UnknownCoordinateSystem() as ISpatialReference; 

               // IProjectedCoordinateSystem m_ProjectedCoordinateSystem; 

                //ISpatialReferenceFactory2 spatRefFact = new SpatialReferenceEnvironmentClass(); 

                //m_ProjectedCoordinateSystem = spatRefFact.CreateProjectedCoordinateSystem((int)esriSRProjCSType.esriSRProjCS_WGS1984UTM_49N); 

                rasterWorkspace2 = createRasterWorkspace(directoryName); 

 

                 

               // int rowCount = (int)Math.Round((maxX - minX)/cellSize)+1; 

             //   int columCount = (int)Math.Round((maxY - minY) / cellSize)+1; 

                //经纬度计算 

                double dLon = (maxLon-minLon)/rowCount; 

                double dlat=(maxLat-minLat)/columCount; 

                rasterDataset = rasterWorkspace2.CreateRasterDataset(fileName, "GRID", originPoint, rowCount, columCount, cellSize, cellSize, 1, rstPixelType.PT_FLOAT, SpRef, true);//"IMAGINE Image",rstPixelType.PT_UCHAR,new UnknownCoordinateSystemClass() 

                IRawPixels rawPixels = null; 

                IPixelBlock3 pixelBlock3 = null; 

                IPnt pixelBlockOrigin = null; 

                IPnt pixelBlockSize = null; 

                IRasterBandCollection rasterBandCollection; 

                IRasterProps rasterProps; 

                // QI for IRawPixels and IRasterProps     

                rasterBandCollection = (IRasterBandCollection)rasterDataset; 

                rawPixels = (IRawPixels)rasterBandCollection.Item(0); 

                rasterProps = (IRasterProps)rawPixels; 

 

                // Create pixelblock     

                pixelBlockOrigin = new DblPntClass(); 

                pixelBlockOrigin.SetCoords(0, 0); 

                pixelBlockSize = new DblPntClass(); 

                pixelBlockSize.SetCoords(rasterProps.Width, rasterProps.Height); 

                pixelBlock3 = rawPixels.CreatePixelBlock(pixelBlockSize) as IPixelBlock3; 

                // Read pixelblock     

                rawPixels.Read(pixelBlockOrigin, (IPixelBlock)pixelBlock3); 

                // Get pixeldata array      

                System.Array pixelData; 

                pixelData = (System.Array)pixelBlock3.get_PixelDataByRef(0); 

                // Loop through all the pixels and assign value    

                for (int i = 0; i  rasterProps.Width; i++)//这个数组中就是输入的Value值 

                { 

                    for (int j = 0; j  rasterProps.Height; j++) 

                    { 

                        //tempPt.X = minX + i * cellSize; 

                      //  object oValue =CountourpixelVale.GetValue(0, 0); 

                      //  tempPt.Y = minY + j * cellSize; 

                        Filldi(minLon + i * dLon, minLat + j * dlat); 

                        float temp = (float)GIDS(minX + i * cellSize, minY + j * cellSize, Convert.ToDouble(CountourpixelVale.GetValue(i, j))); 

                        ////求均方差 

                        //for (int k=0;k<dt.Rows.Count;k++) 

                        //{ 

                        //    MeanSquareDeviation += (temp - Railfall[k])*(temp - Railfall[k]); 

                        //} 

                        //MeanSquareDeviation = Math.Sqrt(MeanSquareDeviation / dt.Rows.Count); 

                        ////均方差写入文件 

                        //sw.WriteLine(i + "," + j + "," + MeanSquareDeviation); 

                        

                        pixelData.SetValue(temp, i, j);//Convert.ToByte((i * j) % 255) 

                    } 

                } 

 

               for (int k=0;k<dt.Rows.Count;k++) 

               { 

                   int rr = Convert.ToInt32(Math.Ceiling((Ptx[k] - minX) / cellSize)); 

                   int cc=Convert.ToInt32(Math.Ceiling((Pty[k]-minY)/cellSize)); 

                   double pp; 

                   if (rr  rowCount && cc  columCount) 

                   { 

                       pp = Convert.ToDouble(pixelData.GetValue(rr, cc)); 

                       sw.WriteLine(PtName[k]+","+Ptx[k] + "," + Pty[k] + "," + Math.Abs(pp - Railfall[k])); 

                   } 

                } 

                sw.Close(); 

                pixelBlock3.set_PixelData(0, (System.Object)pixelData); 

                // Write the pixeldata back      

                //System.Object cachePointer; 

                //cachePointer = rawPixels.AcquireCache(); 

                rawPixels.Write(pixelBlockOrigin, (IPixelBlock)pixelBlock3); 

                //rawPixels.ReturnCache(cachePointer); 

                 

                //IRasterLayer pRasterLayer = new RasterLayerClass(); 

                //pRasterLayer.CreateFromDataset(rasterDataset); 

                //ILayer pLayer; 

                //pLayer = pRasterLayer; 

                //pLayer.Name = "Result"; 

                //axMapControl1.Map.AddLayer(pLayer); 

                //axMapControl1.ActiveView.Refresh(); 

 

                //解除引用 

                rawPixels = null; 

                rasterDataset= null; 

                pixelBlock3= null; 

                rasterBandCollection= null; 

                rasterProps = null; 

                GC.Collect();//强制回收 

                IRasterLayer m_raster = new RasterLayerClass(); 

                m_raster.CreateFromFilePath("c:\\temp\\" + RainfallValue); 

                axMapControl1.AddLayer(m_raster);    

                //IRasterEdit rasterEdit; //保存 

                //rasterEdit = pRasterLayer.Raster as IRasterEdit; 

                //rasterEdit.Write(pixelBlockOrigin, pixelBlock3 as IPixelBlock); 

                //rasterEdit.Refresh(); 

 

                //return rasterDataset; 

            } 

            catch (Exception ex) 

            { 

                sw.Close(); 

                MessageBox.Show("My God, Error!" + ex.Message); 

                //System.Diagnostics.Debug.WriteLine(ex.Message); 

                return;// null; 

            } 

        } 

        double Cx = 0.0007; 

        double Cy = 0.0006; 

        double Cz = 0.5199; 

        private double GIDS(double x,double y,double z)//反梯度法求值 

        { 

            double diSum=0; 

            double tempSum=0; 

            for (int i = 0; i  dt.Rows.Count; i++) 

            { 

                diSum += 1 / (di[i] * di[i]); 

                tempSum += (Railfall[i] + (x - Ptx[i]) * Cx + (y - Pty[i]) * Cy + Cz*(z - Ptz[i])) / (di[i] * di[i]); 

            } 

            tempSum = tempSum / diSum; 

            return tempSum; 

        } 

        private void Filldi(double Templon,double Templat) 

        { 

            for (int i = 0; i  dt.Rows.Count; i++) 

            { 

                di[i] = GetDistance(Templon, Templat,lon[i],lat[i]); 

            } 

        } 

        public IRasterWorkspace2 createRasterWorkspace(string pathName) 

        { 

            // Create RasterWorkspace   

            IWorkspaceFactory workspaceFactory = new RasterWorkspaceFactoryClass(); 

            return workspaceFactory.OpenFromFile(pathName, 0) as IRasterWorkspace2; 

        } 

        OleData oledata = new OleData(); 

        DataTable dt =new DataTable();  

 

        private void Init()//初始化复制，从数据库 

        { 

            for(int i = 0; i  dt.Rows.Count; i++) 

            { 

                 PtName[i]=dt.Rows[i][1].ToString(); 

                 //Cx[i] = Convert.ToDouble(dt.Rows[i][2]);//转换前X坐标 

                 //Cy[i] = Convert.ToDouble(dt.Rows[i][3]); 

                 Ptx[i] = Convert.ToDouble(dt.Rows[i][4]); 

                 Pty[i] = Convert.ToDouble(dt.Rows[i][5]); 

                 lon[i] = Convert.ToDouble(dt.Rows[i][6]); 

                 lat[i] = Convert.ToDouble(dt.Rows[i][7]); 

                 Ptz[i] = Convert.ToDouble(dt.Rows[i][8]); 

                 Railfall[i] = Convert.ToDouble(dt.Rows[i][9]);//降雨量100 

            } 

        } 

        private void ResetRailfall() 

        { 

            for (int i = 0; i  dt.Rows.Count; i++) 

            { 

                Railfall[i] = Convert.ToDouble(dt.Rows[i][RainfallValue]);//降雨量100 

            } 

        } 

        private double GetMinX() 

        { 

            double min = Ptx[0]; 

            for (int i = 1; i  dt.Rows.Count;i++ ) 

            { 

                if (Ptx[i]  min) 

                    min = Ptx[i]; 

            } 

            return min; 

        } 

        private double GetMinY() 

        { 

            double min = Pty[0]; 

            for (int i = 1; i  dt.Rows.Count; i++) 

            { 

                if (Pty[i]  min) 

                    min = Pty[i]; 

            } 

            return min; 

        } 

        private double GetMaxX() 

        { 

            double max = Ptx[0]; 

            for (int i = 1; i  dt.Rows.Count; i++) 

            { 

                if (Ptx[i] > max) 

                    max = Ptx[i]; 

            } 

            return max; 

        } 

        private double GetMaxY() 

        { 

            double max = Pty[0]; 

            for (int i = 1; i  dt.Rows.Count; i++) 

            { 

                if (Pty[i] > max) 

                    max = Pty[i]; 

            } 

            return max; 

        } 

        private double GetMinLon() 

        { 

            double min = lon[0]; 

            for (int i = 1; i  dt.Rows.Count; i++) 

            { 

                if (lon[i]  min) 

                    min = lon[i]; 

            } 

            return min; 

        } 

        private double GetMinLat() 

        { 

            double min = lat[0]; 

            for (int i = 1; i  dt.Rows.Count; i++) 

            { 

                if (lat[i]  min) 

                    min = lat[i]; 

            } 

            return min; 

        } 

        private double GetMaxLon() 

        { 

            double max = lon[0]; 

            for (int i = 1; i  dt.Rows.Count; i++) 

            { 

                if (lon[i] > max) 

                    max = lon[i]; 

            } 

            return max; 

        } 

        private double GetMaxLat() 

        { 

            double max = lat[0]; 

            for (int i = 1; i  dt.Rows.Count; i++) 

            { 

                if (lat[i] > max) 

                    max = lat[i]; 

            } 

            return max; 

        } 

        const double PAI=3.1415926535; 

        const double EARTH_RADUS = 6371004; 

        double GetDistance(double a_dFirLon, double a_dFirLat, double a_dSecdLon, double a_dSecdLat) 

        { 

            double dRes = 0; 

            double dTemp = 0, dTemp2 = 0; 

 

            dTemp = Math.Cos(a_dSecdLat * PAI / 180) * Math.Cos(a_dFirLat * PAI / 180) * Math.Cos((a_dFirLon - a_dSecdLon) * PAI / 180) + Math.Sin(a_dSecdLat * PAI / 180) * Math.Sin(a_dFirLat * PAI / 180); 

            dTemp2 = EARTH_RADUS * Math.Acos(dTemp); 

	        dRes = dTemp2; 

	        return dRes; 

        } 

    } 

}