方法（一）打开自定义栅格数据
if (mapControlMain.FocusMap.LayerCount == )return;
IRasterLayer pRasterLayer = mapControlMain.ActiveView.CurrentLayer as IRasterLayer;
if (pRasterLayer == null)return;
IRasterLayer memLayer = OpenMemRasterLayer(pRasterLayer, , );
if(memLayer!=null)
{
mapControlMain.FocusMap.AddLayer(memLayer as ILayer);
mapControlMain.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll);
}
/// <summary>
/// 自定义栅格数据读取（通过内存RasterData构建RasterLayer）
/// </summary>
/// <param name="rasterLayer">栅格数据路径</param>
/// <param name="startX">读取起始位置X</param>
/// <param name="startY">读取起始位置Y</param>
/// <param name="?"></param>
/// <returns></returns>
private IRasterLayer OpenMemRasterLayer(IRasterLayer rasterLayer, int startX, int startY)
{
IRasterLayer memRasterLayer = null;
object array = null;//存储读取出的数据
//栅格图层参数读取
IRasterDataset rasterDataset = rasterLayer.Dataset;
int bandCount = rasterDataset.GetBandCount();
PixelDataType rasterDataType = rasterDataset.GetRasterBand().GetRasterDataType();
//读取栅格数据块大小
int width = ;
int height = ;
switch (rasterDataType)
{
case PixelDataType.UInt16:
array = new UInt16[width * height * bandCount];
break;
case PixelDataType.Byte:
array = new Byte[width * height * bandCount];
break;
default:
break;
}
int[] bandMap = new int[bandCount];
for (int i = ; i < bandCount; i++)
{
bandMap[i] = i + ;
}
//将pRasterLayer中的数据读取至数组中
bool oK = rasterDataset.Read(startX, startY, width, height, array, width, height, rasterDataType, bandCount, bandMap); //栅格读取
double[] geoTrans = rasterDataset.GetGeoTransform();
//创建内存数据集
IRasterDataset newRDataset = DatasetFactory.CreateRasterDataset("", width, height, bandCount, rasterDataType, "MEM", null);
//数组中的数据写入内存数据集中
oK = newRDataset.Write(, , width, height, array, width, height, rasterDataType, bandCount, bandMap); //栅格写入
//设置参数
newRDataset.SetGeoTransform(geoTrans);//设置仿射变换参数
newRDataset.SpatialReference = rasterDataset.SpatialReference;//设置空间参考
//创建RasterLayer
memRasterLayer = PIE.Carto.LayerFactory.CreateDefaultRasterLayer(newRDataset);
return memRasterLayer;
}
