ArcGIS栅格处理总结教程.docx

ArcGIS栅格处理总结 RealStyle-赵晓鲁 2012年9月28日 处理栅格数据用到的新的接口有:iconditionalop,Ipnt，IRasterDescriptor, IextractionOp，IrasterLayerExport，IRasterAnalysisEnvironment，ILogicalOp，IReclassOp，IMathOp，INumberRemap，IRemap，ISurfaceop，Irawpixels（很有用） 坐标值位置： 地图坐标最小X、Y值在左下角，IEnvelope最小值是左下角，最大值是右上角，栅格数据最小行列值是在左上角，最大行列值是在右下角。 IRasterAnalysisEnvironment:一般和IextractionOp一起用，用来设置提取出来的栅格数据的一些属性信息。 遍历栅格数据所有数据，将像素值存储在二维数组中 IWorkspaceFactory pWksfact = new RasterWorkspaceFactoryClass(); IRasterWorkspace pRwks = pWksfact.OpenFromFile(@F:\栅格数据处理集\Temp\data, 0) as IRasterWorkspace ; IRaster pRaster = pRwks.OpenRasterDataset(003647.tif).CreateDefaultRaster(); IRasterProps pRprops = pRaster as IRasterProps ; IRasterBandCollection pRbcoll = pRaster as IRasterBandCollection ; int pRbcount = pRbcoll.Count ; IPixelBlock3 pPblock3 = null; GC.Collect();//释放内存，强制回收垃圾内存 byte[,] pByte = new byte[pRprops.Width, pRprops.Height * pRbcount];//内存占有量为byte类型所占的字节数*pRprops.Width*pRprops.Height * pRbcount System.Array pArray = null;//定义一个数组,用来存储pPblock3.get_PixelData(pTindex)获得的像素块的值 for (int pTindex = 0; pTindex pRbcount; pTindex++)//遍历每个波段 { IRasterCursor pRcursor = pRaster.CreateCursor();//128行，全部列是一个pRcursor.Next()， do { pPblock3 = pRcursor.PixelBlock as IPixelBlock3;//将每个128行的值放在pPblock3中 pArray = pPblock3.get_PixelData(pTindex) as System.Array;//将相应波段的数值存储在数组中 for (int pW = 0; pW pArray.GetLength(0); pW++)//列数 { for (int pH = 0; pH pArray.GetLength(1); pH++)//行数 { pByte[pW, pH * (pTindex + 1)] = Convert.ToByte (pArray.GetValue(pW, pH));//给二维数组赋值 } } } while (pRcursor.Next());//获得下一个128行栅格数据 ESRI.ArcGIS.ADF.ComReleaser.ReleaseCOMObject(pRcursor );//对对象进行内存释放 }//遍历每个波段 for (int pT = 0; pT pByte.GetLength(0); pT++)//遍历二维数组 { for (int pTT = 0; pTT pByte.GetLength(1); pTT++) { int pp = Convert.ToInt32(pByte.GetValue(pT, pTT)); } }//遍历二维数组 如果想将Tif转换为bmp，仅仅需要用到的AO的接口是Isaveas。 IWorkspaceFactory pWksfact=newRasterWorkspaceFactoryClass(); IRasterWorkspace pRwks = pWksfact.OpenFromFile(@F:\Temp,0 ) as IRasterWorkspace