C3+Ae9.2写的生成多波段栅格数据的函数
参数的实例为：
//savepath = "C:\\tmpdata\\";
//savename = "testfile.img";
//FYIR3FileInfo 为我保存的数据。
具体到其他使用要修改的，这个是生成四波段的，如果用的不是四波段，要修改里面的波段数据。
生成的效果是可以的，既可以在ae的mapcontrol控件中显示，也可以输出为img格式的文件。
struct RasterInfo
{
public int lNumbands;
}
/// <summary>
/// 生成多波段IR3L栅格文件数据功能
/// </summary>
/// <param name="FYIR3FileInfo"></param>
/// <param name="savepath"></param>
/// <param name="savefilename"></param>
/// <returns></returns>
public ILayer CreateThreeBandIR3L(IR3LFile FYIR3FileInfo, string savepath, string savefilename)
{
//The file must not be there. Otherwise, dataset can not be got.
if (File.Exists(savepath + savefilename))
{
MessageBox.Show("File Exist!");
IRasterLayer pRasterLayer = new RasterLayerClass();
pRasterLayer.CreateFromFilePath(savepath + savefilename);
ILayer pLayer;
pLayer = pRasterLayer;
pLayer.Name = "New Raster";
return pLayer;
}
// This function creates a new img file in the given workspace
// and then assigns pixel values
try
{
IRasterDataset rasterDataset = null;
IPoint originPoint = new PointClass();
originPoint.PutCoords(FYIR3FileInfo.dblLeftDownLong, FYIR3FileInfo.dblLeftDownLat);
// Create the dataset
IRasterWorkspace2 rasterWorkspace2 = null;
double xPixelSize;
double yPixelSize;
//像元大小
xPixelSize = (FYIR3FileInfo.dblcenterLong - FYIR3FileInfo.dblLeftDownLong) / (FYIR3FileInfo.intNumX / 2);
yPixelSize = (FYIR3FileInfo.dblcenterLat - FYIR3FileInfo.dblLeftDownLat) / (FYIR3FileInfo.intNumY / 2);
//坐标系
ISpatialReferenceFactory2 spatRefFact = new SpatialReferenceEnvironmentClass();
//m_ProjectedCoordinateSystem = spatRefFact.CreateProjectedCoordinateSystem((int)esriSRProjCSType.esriSRProjCS_WGS1984UTM_49N);
ESRI.ArcGIS.Geometry.IGeographicCoordinateSystem m_ProjectedCoordinateSystem;
m_ProjectedCoordinateSystem = spatRefFact.CreateGeographicCoordinateSystem((int)esriSRGeoCSType.esriSRGeoCS_WGS1984);
rasterWorkspace2 = createRasterWorkspace(savepath);
rasterDataset = rasterWorkspace2.CreateRasterDataset(savefilename, "IMAGINE Image", originPoint, FYIR3FileInfo.intNumX, FYIR3FileInfo.intNumY, xPixelSize, yPixelSize, 4, rstPixelType.PT_UCHAR, m_ProjectedCoordinateSystem, true);//rstPixelType.PT_UCHAR,new UnknownCoordinateSystemClass()
RasterInfo rstInfo = new RasterInfo();
rstInfo.lNumbands = 4;
IRawPixels[] rawPixels = new IRawPixels[rstInfo.lNumbands];
IPixelBlock3[] pixelBlock3 = new IPixelBlock3[rstInfo.lNumbands];
IPnt pixelBlockOrigin = null;
IPnt pixelBlockSize = null;
IRasterBandCollection rasterBandCollection;
IRasterProps rasterProps;
// QI for IRawPixels and IRasterProps
rasterBandCollection = (IRasterBandCollection)rasterDataset;
for (int i = 0; i < rstInfo.lNumbands; i ++ )
{
rawPixels[i] = (IRawPixels)rasterBandCollection.Item(i);
}
//调用调色板
string PalFilePath;
CreateIMG.Class.Clsfengyun.Pal TmppalInfo;
TmppalInfo = new CreateIMG.Class.Clsfengyun.Pal();
PalFilePath = Application.StartupPath + "\\系统文件\\pal\\I-01.pal";
ReadPal(PalFilePath, ref TmppalInfo);
int tmpColorValue;
System.Array[] pixelData = new System.Array[rstInfo.lNumbands];
rasterProps = (IRasterProps)rawPixels[0];
for (int i = 0; i < rstInfo.lNumbands; i++)
{
rasterProps = (IRasterProps)rawPixels[i];
// Create pixelblock
pixelBlockOrigin = new DblPntClass();
pixelBlockOrigin.SetCoords(0, 0);
pixelBlockSize = new DblPntClass();
pixelBlockSize.SetCoords(rasterProps.Width, rasterProps.Height);
pixelBlock3[i] = (IPixelBlock3)rawPixels[i].CreatePixelBlock(pixelBlockSize);
// Read pixelblock
rawPixels[i].Read(pixelBlockOrigin, (IPixelBlock)pixelBlock3[i]);
// Get pixeldata array
pixelData[i] = (System.Array)pixelBlock3[i].get_PixelDataByRef(0);
for (int ii = 0; ii < rasterProps.Width; ii++)
for (int jj = 0; jj < rasterProps.Height; jj++)
{
if (i == 0)
{
tmpColorValue = TmppalInfo.PalColor[FYIR3FileInfo.data[ii, jj]].tmpR;
pixelData[i].SetValue(Convert.ToByte(tmpColorValue), ii, jj);
}
else if (i == 1)
{
tmpColorValue = TmppalInfo.PalColor[FYIR3FileInfo.data[ii, jj]].tmpG;
pixelData[i].SetValue(Convert.ToByte(tmpColorValue), ii, jj);
}
else if (i == 2)
{
tmpColorValue = TmppalInfo.PalColor[FYIR3FileInfo.data[ii, jj]].tmpB;
pixelData[i].SetValue(Convert.ToByte(tmpColorValue), ii, jj);
}
else
{
tmpColorValue = TmppalInfo.PalColor[FYIR3FileInfo.data[ii, jj]].tmpValue;
pixelData[i].SetValue(Convert.ToByte(tmpColorValue), ii, jj);
}
}
pixelBlock3[i].set_PixelData(0, (System.Object)pixelData[i]);
System.Object cachePointer;
cachePointer = rawPixels[i].AcquireCache();
rawPixels[i].Write(pixelBlockOrigin, (IPixelBlock)pixelBlock3[i]);
rawPixels[i].ReturnCache(cachePointer);
}
IRasterLayer pRasterLayer = new RasterLayerClass();
pRasterLayer.CreateFromDataset(rasterDataset);
ILayer pLayer;
pLayer = pRasterLayer;
pLayer.Name = "New Raster";
return pLayer;
}
catch (System.Exception ex)
{
MessageBox.Show(ex.Message);
//System.Diagnostics.Debug.WriteLine(ex.Message);
return null;
}
}