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 = newRasterWorkspaceFactoryClass();IRasterWorkspace pRwks = pWksfact.OpenFromFile(@F:\栅格数据处理集\Temp\data, 0) asIRasterWorkspace ;IRaster pRaster = pRwks.OpenRasterDataset(003647.tif).CreateDefaultRaster();IRasterProps pRprops = pRaster asIRasterProps ;IRasterBandCollection pRbcoll = pRaster asIRasterBandCollection ;int pRbcount = pRbcoll.Count ;IPixelBlock3 pPblock3 = null;GC.Collect();//释放内存，强制回收垃圾内存byte[,] pByte = newbyte[pRprops.Width, pRprops.Height * pRbcount];//内存占有量为byte类型所占的字节数*pRprops.Width*pRprops.Height * pRbcountSystem.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 asIPixelBlock3;//将每个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。IWorkspaceFactorypWksfact=newRasterWorkspaceFactoryClass();IRasterWorkspace pRwks = pWksfact.OpenFromFile(@F:\Temp,0 ) asIRasterWorkspace ;IRasterDataset pRdataset = pRwks.OpenRasterDataset(c