DataSourcesRaster

IRasterWorkspace OpenRasterWorkspace(string filePath)
{
    Type factoryType = Type.GetTypeFromProgID(
        "esriDataSourcesRaster.RasterWorkspaceFactory");
    IWorkspaceFactory wsFactory = (IWorkspaceFactory)Activator.CreateInstance
        (factoryType);
    IRasterWorkspace rasterWS = (IRasterWorkspace)wsFactory.OpenFromFile(filePath, 0)
        ;
    return rasterWS;
}

Private Function OpenRasterWorkspace(ByVal filePath As String) As IRasterWorkspace
    Dim factoryType As Type = Type.GetTypeFromProgID("esriDataSourcesRaster.RasterWorkspaceFactory")
    Dim wsFactory As IWorkspaceFactory = (IWorkspaceFactory)Activator.CreateInstance(factoryType)
    Dim rasterWS As IRasterWorkspace = CType(wsFactory.OpenFromFile(filePath, 0), IRasterWorkspace)
    Return rasterWS
End Function

IRasterWorkspaceEx OpenFileGDBWorkspace(string filePath)
{
    Type factoryType = Type.GetTypeFromProgID(
        "esriDataSourcesGDB.FileGDBWorkspaceFactory");
    IWorkspaceFactory wsFactory = (IWorkspaceFactory)Activator.CreateInstance
        (factoryType);
    IRasterWorkspaceEx rasterWS = (IRasterWorkspaceEx)wsFactory.OpenFromFile
        (filePath, 0);
    return rasterWS;
}

Private Function OpenFileGDBWorkspace(ByVal filePath As String) As IRasterWorkspaceEx
    Dim factoryType As Type = Type.GetTypeFromProgID("esriDataSourcesGDB.FileGDBWorkspaceFactory")
    Dim wsFactory As IWorkspaceFactory = CType(Activator.CreateInstance(factoryType), IWorkspaceFactory)
    Dim rasterWS As IRasterWorkspaceEx = CType(wsFactory.OpenFromFile(filePath, 0), IRasterWorkspaceEx)
    Return rasterWS
End Function

IRasterWorkspaceEx OpenSDEWorkspace(string sServer, string sInstance, string sDB,
    string sUser, string sPasswd, string sVersion)
{
    Type factoryType = Type.GetTypeFromProgID(
        "esriDataSourcesGDB.SdeWorkspaceFactory");
    IWorkspaceFactory wsFactory = (IWorkspaceFactory)Activator.CreateInstance
        (factoryType);
    Dim sdeConnection As IPropertySet = New PropertySetClass();
    sdeConnection.SetProperty("Server", sServer);
    sdeConnection.SetProperty("Instance", sInstance);
    sdeConnection.SetProperty("Database", sDB);
    sdeConnection.SetProperty("User", sUser);
    sdeConnection.SetProperty("Password", sPasswd);
    sdeConnection.SetProperty("Version", sVersion);
    IRasterWorkspaceEx ws = CType(wsFactory.Open(sdeConnection, 0),
        IRasterWorkspaceEx);
    Return ws;
}

Private Function OpenSDEWorkspace(ByVal sServer As String, ByVal sInstance As String, ByVal sDB As String, ByVal sUser As String, ByVal sPasswd As String, ByVal sVersion As String) As IRasterWorkspaceEx
    Dim factoryType As Type = Type.GetTypeFromProgID("esriDataSourcesGDB.SdeWorkspaceFactory")
    Dim wsFactory As IWorkspaceFactory = (IWorkspaceFactory)Activator.CreateInstance(factoryType)
    Dim sdeConnection As IPropertySet = New PropertySetClass()
    sdeConnection.SetProperty ("Server", sServer)
    sdeConnection.SetProperty ("Instance", sInstance)
    sdeConnection.SetProperty("Database", sDB)
    sdeConnection.SetProperty ("User", sUser)
    sdeConnection.SetProperty ("Password", sPasswd)
    sdeConnection.SetProperty ("Version", sVersion)
    Dim ws As IRasterWorkspaceEx = CType(wsFactory.Open(sdeConnection, 0), IRasterWorkspaceEx)
    Return ws
End Function

//Open a raster dataset from a file system.
IRasterWorkspace ws = OpenRasterWorkspace(@"d:\data");
IRasterDataset rasterDataset = ws.OpenRasterDataset("airphoto.img");

'Open a raster dataset from a file system.
Dim ws As IRasterWorkspace = OpenRasterWorkspace ("d:\data")
Dim rasterDataset As IRasterDataset = ws.OpenRasterDataset("airphoto.img")

//Open a raster catalog from a file geodatabase workspace.
IRasterWorkspaceEx wsFGDB = OpenFileGDBWorkspace(@"d:\data\fgdb_images.gdb");
IRasterCatalog rasterCatalog = wsFGDB.OpenRasterCatalog("RedlandImages");

//Open a raster dataset from an SDE workspace.
IRasterWorkspaceEx wsSDE = OpenSDEWorkspace("myserver", "5151", "myuser",
    "mydatabase", "mypassws", "SDE.DEFAULT");
IRasterDataset rasterdataset1 = wsSDE.OpenRasterDataset("airphoto");

'Open a raster catalog from a file geodatabase workspace.
Dim wsFGDB As IRasterWorkspaceEx = OpenFileGDBWorkspace("d:\data\fgdb_images.gdb")
Dim rasterCatalog As IRasterCatalog = wsFGDB.OpenRasterCatalog("RedlandImages")

'Open a raster dataset from an SDE workspace.
Dim wsSDE As IRasterWorkspaceEx = OpenSDEWorkspace("myserver", "5151", "myuser", "mydatabase", "mypassws", "SDE.DEFAULT")
Dim rasterdataset1 As IRasterDataset = wsSDE.OpenRasterDataset("airphoto")

public static IMosaicDataset OpenMosaicDataset(IWorkspace ws, String mdName)
{
    IMosaicWorkspaceExtensionHelper MosaicWsHelper = new
        MosaicWorkspaceExtensionHelperClass();
    IMosaicWorkspaceExtension mosaicWsExtension = MosaicWsHelper.FindExtension(ws);
    IMosaicDataset mosaicDataset = mosaicWsExtension.MosaicDataset(mdName);
    return mosaicDataset;
}

Private Function OpenMosaicDataset(ws As IWorkspace, mdName As String)
    Dim MosaicWsHelper As IMosaicWorkspaceExtensionHelper = New MosaicWorkspaceExtensionHelperClass()
    Dim mosaicWsExtension As IMosaicWorkspaceExtension = MosaicWsHelper.FindExtension(ws)
    Dim mosaicDataset = mosaicWsExtension.OpenMosaicDataset(mdName)
End Function

IRasterWorkspace2 ws = (IRasterWorkspace2)OpenRasterWorkspace(@"d:\data");
IPoint origin = new PointClass();
origin.PutCoords(100, 100);
ISpatialReference sr = new UnknownCoordinateSystemClass();
IRasterDataset newRaster = ws.CreateRasterDataset("myimage.tif", "TIFF", origin,
    1024, 1024, 30, 30, 1, rstPixelType.PT_UCHAR, sr, true);

Dim ws As IRasterWorkspace2 = CType(OpenRasterWorkspace("d:\data"), IRasterWorkspace2)
Dim origin As IPoint = New PointClass()
origin.PutCoords(100, 100)
Dim sr As ISpatialReference = New UnknownCoordinateSystemClass()
Dim newRaster As IRasterDataset = ws.CreateRasterDataset("myimage.img", "IMAGINE Image", origin, 1024, 1024, 30, 30, 1, rstPixelType.PT_UCHAR, sr, True)

void GetDatasetXform(IRasterDataset rasterDataset)
{
    //Get the Xform.
    IGeoDataset2 geoDataset = (IGeoDataset2)rasterDataset;
    IGeodataXform xform = geoDataset.GeodataXform;

    //Get the native extent.
    IEnvelope nativeExtent = geoDataset.NativeExtent;

    //Get the native spatial reference.
    ISpatialReference nativeSR = geoDataset.NativeSpatialReference;
}

Private Sub GetDatasetXform(ByVal rasterDataset As IRasterDataset)
    'Get the Xform.
    Dim geoDataset As IGeoDataset2 = CType(rasterDataset, IGeoDataset2)
    Dim xform As IGeodataXform = geoDataset.GeodataXform
    'Get the native extent.
    Dim nativeExtent As IEnvelope = geoDataset.NativeExtent
    'Get the native spatial reference.
    Dim nativeSR As ISpatialReference = geoDataset.NativeSpatialReference
End Sub

void SaveAs(IRasterDataset rasterDataset)
{
    ISaveAs saveAs = (ISaveAs)rasterDataset;
    IWorkspace ws = (IWorkspace)OpenRasterWorkspace(@"d:\data");
    saveAs.SaveAs("mygrid", ws, "GRID");
}

Private Sub SaveAs(ByVal rasterDataset As IRasterDataset)
    Dim saveAs As ISaveAs = CType(rasterDataset, ISaveAs)
    Dim ws As IWorkspace = CType(OpenRasterWorkspace("d:\data"), IWorkspace)
    saveAs.SaveAs("mygrid", ws, "GRID")
End Sub

void SaveAsWithCompression(IRasterDataset rasterDataset)
{
    ISaveAs2 saveAs = (ISaveAs2)rasterDataset;

    //Set the storage properties.
    IRasterStorageDef2 storageDef = (IRasterStorageDef2)new RasterStorageDefClass();
    storageDef.CompressionType =
        esriRasterCompressionType.esriRasterCompressionJPEG2000;
    storageDef.CompressionQuality = 30;
    storageDef.Tiled = true;
    storageDef.TileHeight = 64;
    storageDef.TileWidth = 64;

    //Set the output workspace.
    IWorkspace ws = (IWorkspace)OpenFileGDBWorkspace(@"d:\data\fgdb.gdb");
    saveAs.SaveAsRasterDataset("compressedRaster", ws, "GDB", storageDef);
}

Private Sub SaveAsWithCompression(ByVal rasterDataset As IRasterDataset)
    Dim saveAs As ISaveAs2 = CType(rasterDataset, ISaveAs2)

    'Set the storage properties.
    Dim storageDef As IRasterStorageDef2 = CType(New RasterStorageDefClass(), IRasterStorageDef2)
    storageDef.CompressionType = esriRasterCompressionType.esriRasterCompressionJPEG2000
    storageDef.CompressionQuality = 30
    storageDef.Tiled = True
    storageDef.TileHeight = 64
    storageDef.TileWidth = 64

    'Set the output workspace.
    Dim ws As IWorkspace = CType(OpenFileGDBWorkspace("d:\data\fgdb.gdb"), IWorkspace)
    saveAs.SaveAsRasterDataset("compressedRaster", ws, "GDB", storageDef)
End Sub

void GetRasterDataset(IRaster raster)
{
    IRasterBandCollection bandCol = (IRasterBandCollection)raster;
    IRasterBand rasterBand = bandCol.Item(0);
    IRasterDataset rasterDataset = rasterBand.RasterDataset;
}

Private Sub GetRasterDataset(ByVal raster As IRaster)
    Dim bandCol As IRasterBandCollection = CType(raster, IRasterBandCollection)
    Dim rasterBand As IRasterBand = bandCol.Item(0)
    Dim rasterDataset As IRasterDataset = rasterBand.RasterDataset
End Sub

void GetColormap(IRasterDataset rasterDataset)
{
    IRasterBandCollection bandCol = (IRasterBandCollection)rasterDataset;
    IRasterBand rasterBand = bandCol.Item(0);
    IRasterColormap rasterColormap = rasterBand.Colormap;

    IRaster2 raster = (IRaster2)rasterDataset.CreateDefaultRaster();
    rasterColormap = raster.Colormap;
}

Private Sub GetColormap(ByVal rasterDataset As IRasterDataset)
    Dim bandCol As IRasterBandCollection = CType(rasterDataset, IRasterBandCollection)
    Dim rasterBand As IRasterBand = bandCol.Item(0)
    Dim rasterColormap As IRasterColormap = rasterBand.Colormap

    Dim raster As IRaster2 = CType(rasterDataset.CreateDefaultRaster(), IRaster2)
    rasterColormap = raster.Colormap
End Sub

void RasterAttributeTable(IRaster raster)
{
    IRaster2 raster = (IRaster2)raster;
    ITable rasterTable = raster.AttributeTable;
}

Private Sub RasterAttributeTable(ByVal raster As IRaster)
    Dim raster As IRaster2 = CType(raster, IRaster2)
    Dim rasterTable As ITable = raster.AttributeTable
End Sub

void RasterAttributeTable(IRasterDataset rasterDataset)
{
    IRasterBandCollection bandCol = (IRasterBandCollection)rasterDataset;
    IRasterBand rasterBand = bandCol.Item(0);
    ITable rasterTable = rasterBand.AttributeTable;
}

Dim bandCol As IRasterBandCollection = CType(rasterDataset, IRasterBandCollection)
Dim rasterBand As IRasterBand = bandCol.Item(0)
Dim rasterTable As ITable = rasterBand.AttributeTable

void CalculateStats(IRasterDataset rasterDataset)
{
    IRasterBandCollection bandCollection = (IRasterBandCollection)rasterDataset;
    IRasterBand rasterBand = bandCollection.Item(0);
    IRasterStatistics rasterStats = rasterBand.Statistics;

    rasterStats.SkipFactorX = 100;
    rasterStats.SkipFactorY = 100;
    rasterStats.Recalculate();
}

Private Sub CalculateStats(ByVal rasterDataset As IRasterDataset)
    Dim bandCollection As IRasterBandCollection = CType(rasterDataset, IRasterBandCollection)
    Dim rasterBand As IRasterBand = bandCollection.Item(0)
    Dim rasterStats As IRasterStatistics = rasterBand.Statistics

    rasterStats.SkipFactorX = 100
    rasterStats.SkipFactorY = 100
    rasterStats.Recalculate()
End Sub

void CreateRaster(IRasterDataset rasterDataset)
{
    //Create a default raster.
    IRaster raster = rasterDataset.CreateDefaultRaster();

    //Create a raster with all bands.
    IRasterDataset2 rasterDataset2 = (IRasterDataset2)rasterDataset;
    IRaster raster2 = rasterDataset2.CreateFullRaster();
}

Private Sub CreateRaster(ByVal rasterDataset As IRasterDataset)
    'Create a default raster.
    Dim raster As IRaster = rasterDataset.CreateDefaultRaster()

    'Create a raster with all bands.
    Dim rasterDataset2 As IRasterDataset2 = CType(rasterDataset, IRasterDataset2)
    Dim raster2 As IRaster = rasterDataset2.CreateFullRaster()
End Sub

void CreateRaster(IRasterDataset rasterDataset)
{
    //Create an empty raster and append bands.
    IRasterBandCollection newRaster = new RasterClass();
    newRaster.AppendBands((IRasterBandCollection)rasterDataset);
}

Private Sub CreateRaster(ByVal rasterDataset As IRasterDataset)
    'Create an empty raster and append bands.
    Dim newRaster As IRasterBandCollection = New RasterClass()
    newRaster.AppendBands(CType(rasterDataset, IRasterBandCollection))
End Sub

void SetRasterProperties(IRaster raster, ISpatialReference sr, IWorkspace ws)
{
    IRasterProps rasterProps = (IRasterProps)raster;
    rasterProps.SpatialReference = sr;
    raster.ResampleMethod = rstResamplingTypes.RSP_BilinearInterpolation;

    ISaveAs saveAs = (ISaveAs)raster;
    saveAs.SaveAs("projectedimage.tif", ws, "TIFF");
}

Private Sub SetRasterProperties(ByVal raster As IRaster, ByVal sr As ISpatialReference, ByVal ws As IWorkspace)
    Dim rasterProps As IRasterProps = CType(raster, IRasterProps)
    rasterProps.SpatialReference = sr
    raster.ResampleMethod = rstResamplingTypes.RSP_BilinearInterpolation

    Dim saveAs As ISaveAs = CType(raster, ISaveAs)
    saveAs.SaveAs("projectedimage.tif", ws, "TIFF")
End Sub

void SetAndSaveAsXform(IRaster2 raster, IGeodataXform xform, IWorkspace ws)
{
    //Get the extent and cell size of the raster.
    IRasterProps rasterProps = (IRasterProps)raster;
    IEnvelope extent = rasterProps.Extent;
    IPnt cellSize = rasterProps.MeanCellSize();
    double xCell = cellSize.X;
    double yCell = cellSize.Y;

    //Set the geodataXform.
    raster.GeodataXform = xform;

    //Transform the cell size first, then the extent. The sequence matters.
    xform.TransformCellsize(esriTransformDirection.esriTransformForward, ref xCell,
        ref yCell, extent);
    xform.TransformExtent(esriTransformDirection.esriTransformForward, extent);

    //Put the transformed extent and cell size on the raster.
    rasterProps.Extent = extent;
    rasterProps.Width = (int)(extent.Width / xCell);
    rasterProps.Height = (int)(extent.Width / yCell);

    //Save the transformed image.
    ISaveAs saveAs = (ISaveAs)raster;
    saveAs.SaveAs(@"c:\temp\result.img", null, "IMAGINE Image");
}

Private Sub SetAndSaveAsXform(ByVal raster As IRaster2, ByVal xform As IGeodataXform, ByVal ws As IWorkspace)
    'Get the extent and cell size of the raster.
    Dim rasterProps As IRasterProps = CType(raster, IRasterProps)
    Dim extent As IEnvelope = rasterProps.Extent
    Dim cellSize As IPnt = rasterProps.MeanCellSize()
    Dim xCell As Double = cellSize.X
    Dim yCell As Double = cellSize.Y

    'Set the geodataXform.
    raster.GeodataXform = xform

    'Transform the cell size first, then the extent. The sequence matters.
    xform.TransformCellsize(esriTransformDirection.esriTransformForward, xCell, yCell, extent)
    xform.TransformExtent(esriTransformDirection.esriTransformForward, extent)

    'Put the transformed extent and cell size on the raster.
    rasterProps.Extent = extent
    rasterProps.Width = CInt(Fix(extent.Width / xCell))
    rasterProps.Height = CInt(Fix(extent.Width / yCell))

    'Save the transformed image.
    Dim saveAs As ISaveAs = CType(raster, ISaveAs)
    saveAs.SaveAs("c:\temp\result.img", Nothing, "IMAGINE Image")

void IdentifyOnRaster(IRaster raster, double xmap, double ymap)
{
    IRaster2 raster2 = (IRaster2)raster;
    int col = raster2.ToPixelColumn(xmap);
    int row = raster2.ToPixelRow(ymap);
    double idValue = (double)raster2.GetPixelValue(0, col, row);
}

Private Sub IdentifyOnRaster(ByVal raster As IRaster, ByVal xmap As Double, ByVal ymap As Double)
    Dim raster2 As IRaster2 = CType(raster, IRaster2)
    Dim col As Integer = raster2.ToPixelColumn(xmap)
    Dim row As Integer = raster2.ToPixelRow(ymap)
    Dim idValue As Double = CDbl(raster2.GetPixelValue(0, col, row))
End Sub

void SetFilter(IRaster raster)
{
    //Create the filter.
    IStockConvolutionFilter pixelFilter = new RasterConvolutionFilterClass();
    pixelFilter.Type = esriRasterFilterTypeEnum.esriRasterFilterSmoothing3x3;

    //Set the filter to the raster.
    IPixelOperation pixelOperation = (IPixelOperation)raster;
    pixelOperation.PixelFilter = (IPixelFilter)pixelFilter;
}

Private Sub SetFilter(ByVal raster As IRaster)
    'Create the filter.
    Dim pixelFilter As IStockConvolutionFilter = New RasterConvolutionFilterClass()
    pixelFilter.Type = esriRasterFilterTypeEnum.esriRasterFilterSmoothing3x3

    'Set the filter to the raster.
    Dim pixelOperation As IPixelOperation = CType(raster, IPixelOperation)
    pixelOperation.PixelFilter = CType(pixelFilter, IPixelFilter)
End Sub

void UsingPixelBlock(IRaster raster)
{
    //Create a pixel block using the size provided by the system.
    IRaster2 raster2 = (IRaster2)raster;
    IPixelBlock pixelBlock = raster2.CreateCursorEx(null).PixelBlock;
    int w = pixelBlock.Width;
    int h = pixelBlock.Height;

    //Read the first pixel block.
    IPnt topleftCorner = new PntClass();
    topleftCorner.SetCoords(0, 0);
    raster.Read(topleftCorner, pixelBlock);

    //Modify one pixel value at locations w-1 and h-1 
    //(assume the raster has a pixel type of uchar).
    IPixelBlock3 pixelBlock3 = (IPixelBlock3)pixelBlock;
    System.Array pixels = (System.Array)pixelBlock3.get_PixelData(0);
    pixels.SetValue(Convert.ToByte(10), w - 1, h - 1);
    pixelBlock3.set_PixelData(0, pixels);

    //Write the modified pixel block to the raster dataset.
    IRasterEdit rasterEdit = (IRasterEdit)raster;
    rasterEdit.Write(topleftCorner, pixelBlock);
}

Private Sub UsingPixelBlock(ByVal raster As IRaster)
    'Create a pixel block using the size provided by the system.
    Dim raster2 As IRaster2 = CType(raster, IRaster2)
    Dim pixelBlock As IPixelBlock = raster2.CreateCursorEx(Nothing).PixelBlock
    Dim w As Integer = pixelBlock.Width
    Dim h As Integer = pixelBlock.Height

    'Read the first pixel block.
    Dim topleftCorner As IPnt = New PntClass()
    topleftCorner.SetCoords(0, 0)
    raster.Read(topleftCorner, pixelBlock)

    'Modify one pixel value at locations w-1 and h-1
    '(assume the raster has a pixel type of uchar).
    Dim pixelBlock3 As IPixelBlock3 = CType(pixelBlock, IPixelBlock3)
    Dim pixels As System.Array = CType(pixelBlock3.get_PixelData(0), System.Array)
    pixels.SetValue(Convert.ToByte (10), w - 1, h - 1)
    pixelBlock3.set_PixelData(0, pixels)

    'Write the modified pixel block to the raster dataset.
    Dim rasterEdit As IRasterEdit = CType(raster, IRasterEdit)
    rasterEdit.Write(topleftCorner, pixelBlock)
End Sub

void SaveAsRasterWithNoData(IRaster raster)
{
    //Set NoData.
    IRasterProps rasterProps = (IRasterProps)raster;
    rasterProps.NoDataValue = 255;

    //Or set a different NoData value for each band.
    int[] nodata =
    {
        255, 255, 236
    };
    rasterProps.NoDataValue = nodata;

    //Save.
    ISaveAs saveAs = (ISaveAs)raster;
    saveAs.SaveAs(@"c:\temp\nodata.img", null, "IMAGINE Image");
}

Private Sub SaveAsRasterWithNoData(ByVal raster As IRaster)
    'Set NoData.
    Dim rasterProps As IRasterProps = CType(raster, IRasterProps)
    rasterProps.NoDataValue = 255

    'Or set a different NoData value for each band.
    Dim nodata As Integer() = { 255, 255, 236 }
    rasterProps.NoDataValue = nodata

    'Save.
    Dim saveAs As ISaveAs = CType(raster, ISaveAs)
    saveAs.SaveAs("c:\temp\nodata.img", Nothing, "IMAGINE Image")
End Sub

void AddRasterToRasterCatalog(IRasterDataset rasterDataset, IRasterCatalog
    rasterCatalog)
{
    IFeatureClass fClass = (IFeatureClass)rasterCatalog;
    IFeature feature = fClass.CreateFeature();

    //Create a raster value.
    IRasterValue rasterValue = new RasterValueClass();
    rasterValue.RasterDataset = rasterDataset;

    feature.set_Value(rasterCatalog.NameFieldIndex, rasterValue);
    feature.Store();
}

Private Sub AddRasterToRasterCatalog(ByVal rasterDataset As IRasterDataset, ByVal rasterCatalog As IRasterCatalog)
    Dim fClass As IFeatureClass = CType(rasterCatalog, IFeatureClass)
    Dim feature As IFeature = fClass.CreateFeature()

    'Create a raster value.
    Dim rasterValue As IRasterValue = New RasterValueClass()
    rasterValue.RasterDataset = rasterDataset

    feature.set_Value(rasterCatalog.NameFieldIndex, rasterValue)
    feature.Store()
End Sub

void GetRasterDataset(IRasterCatalog rasterCatalog, int index)
{
    IFeatureClass fClass = (IFeatureClass)rasterCatalog;
    IRasterCatalogItem catalogItem = (IRasterCatalogItem)fClass.GetFeature(index);
    IRasterDataset rasterDataset = catalogItem.RasterDataset;
}

Private Sub GetRasterDataset(ByVal rasterCatalog As IRasterCatalog, ByVal index As Integer)
    Dim fClass As IFeatureClass = CType(rasterCatalog, IFeatureClass)
    Dim catalogItem As IRasterCatalogItem = CType(fClass.GetFeature(index), IRasterCatalogItem)
    Dim rasterDataset As IRasterDataset = catalogItem.RasterDataset
End Sub

public static IRasterFunction CreateHillshadeFunction(IRasterDataset demDS, Double
    altitude, Double azimuth, Double zFactor)
{
    IRasterFunction hillshadeFunction = new HillshadeFunctionClass();
    IHillshadeFunctionArguments arguments = new HillshadeFunctionArgumentsClass();
    arguments.Altitude = altitude;
    arguments.Azimuth = azimuth;
    arguments.ZFactor = zFactor;
    arguments.DEM = demDS;
    hillshadeFunction.Bind(arguments);
    return hillshadeFunction;
}

Public Function CreateHillshadeFunction(ByVal demDS As IRasterDataset, ByVal altitude As Double, ByVal azimuth As Double, ByVal zFactor As Double) As IRasterFunction
    Dim hillshadeFunction As IRasterFunction = New HillshadeFunctionClass()
    Dim arguments As HillshadeFunctionArguments = New HillshadeFunctionArgumentsClass()
    arguments.Altitude = altitude
    arguments.Azimuth = azimuth
    arguments.ZFactor = zFactor
    arguments.DEM = demDS
    hillshadeFunction.Bind(arguments)
    Return hillshadeFunction
End Function

//Create a function raster dataset name.
IFunctionRasterDatasetName functionDSName = new FunctionRasterDatasetNameClass();

//Set the name string and open the function raster dataset.
functionDSName.FullName = "c:\data\image.afr";
IName name = (IName)functionDSName;
IFunctionRasterDataset ds = (IFunctionRasterDataset)name.Open();

'Create a function raster dataset name.
Dim functionDSName As IFunctionRasterDatasetName = New FunctionRasterDatasetNameClass()

'Set the name string and open the function raster dataset.
functionDSName.FullName = "c:\data\image.afr"
Dim Name As IName = CType(functionDSName, IName)
Dim ds As IFunctionRasterDataset = CType(Name.Open(), IFunctionRasterDataset)

public static Void CreateFunctionRasterDataset(String fullName, IRasterFunction
    rasterFunction, IRasterFunctionArguments functionArguments)
{
    //Create a function raster dataset.
    IFunctionRasterDataset functionDS = new FunctionRasterDatasetClass();

    //Define a name for the function raster dataset.
    IFunctionRasterDatasetName functionDsName = new FunctionRasterDatasetNameClass();
    functionDsName.FullName = fullName;
    functionDS.FullName = functionDsName;

    //Initialize the function raster dataset using a given raster function and its 
    //corresponding function arguments.
    functionDS.Init(rasterFunction, functionArguments);

    //Make permanent.
    ITemporaryDataset tempDS = (ITemporaryDataset)functionDS;
    tempDS.MakePermanent();
}

Public Sub CreateFunctionRasterDataset(ByVal fullName As String , ByVal rasterFunction As IRasterFunction, ByVal functionArguments As IRasterFunctionArguments)
    ' Create a function raster dataset.
    Dim functionDS As IFunctionRasterDataset = New FunctionRasterDatasetClass()

    ' Define a name for the function raster dataset.
    Dim functionDsName As IFunctionRasterDatasetName = New FunctionRasterDatasetNameClass()
    functionDsName.FullName = fullName
    functionDS.FullName = functionDsName

    ' Initialize the function raster dataset using a given raster function and its
    ' corresponding function arguments.
    functionDS.Init(rasterFunction, functionArguments)

    ' Make permanent.
    Dim tempDS As ITemporaryDataset = CType(functionDS, ITemporaryDataset)
    tempDS.MakePermanent()
End Sub

public static IMosaicDataset CreateMosaicDataset(IWorkspace fgdbWorkspace, string
    mdName, ISpatialReference mosaicSrs)
{
    // Create the mosaic dataset creation parameters object.
    ICreateMosaicDatasetParameters creationPars = new
        CreateMosaicDatasetParametersClass();
    // Set the number of bands for the mosaic dataset.
    creationPars.BandCount = 3;
    // Set the pixel type of the mosaic dataset.
    creationPars.PixelType = rstPixelType.PT_UCHAR;
    // Create the mosaic workspace extension helper class.
    IMosaicWorkspaceExtensionHelper mosaicExtHelper = new
        MosaicWorkspaceExtensionHelperClass();
    // Find the right extension from the workspace.
    IMosaicWorkspaceExtension mosaicExt = mosaicExtHelper.FindExtension
        (fgdbWorkspace);
    // Use the extension to create a new mosaic dataset, supplying the 
    // spatial reference and the creation parameters object previously created.
    IMosaicDataset theMosaicDataset = mosaicExt.CreateMosaicDataset(mdName,
        mosaicSrs, creationPars);

    return theMosaicDataset;
}

Public Function CreateMosaicDataset(ByVal fgdbWorkspace As IWorkspace , ByVal mdName As String, ByVal mosaicSrs As ISpatialReference ) As IMosaicDataset
    ' Create the mosaic dataset creation parameters object.
    Dim creationPars As ICreateMosaicDatasetParameters = New CreateMosaicDatasetParametersClass()
    ' Set the number of bands for the mosaic dataset.
    creationPars.BandCount = 3
    ' Set the pixel type of the mosaic dataset.
    creationPars.PixelType = rstPixelType.PT_UCHAR
    ' Create the mosaic workspace extension helper class.
    Dim mosaicExtHelper As IMosaicWorkspaceExtensionHelper = New MosaicWorkspaceExtensionHelperClass()
    ' Find the right extension from the workspace.
    Dim mosaicExt As IMosaicWorkspaceExtension = mosaicExtHelper.FindExtension(fgdbWorkspace)
    ' Use the extension to create a new mosaic dataset, supplying the
    ' spatial reference and the creation parameters object previously created.
    Dim theMosaicDataset As IMosaicDataset = mosaicExt.CreateMosaicDataset(mdName, mosaicSrs, creationPars)
    '
    Return theMosaicDataset
End Function

void TransformPoints(IGeodataXform xform, IPointCollection points, IEnvelope extent,
    ref double xcell, ref double ycell)
{
    xform.TransformPoints(esriTransformDirection.esriTransformForward, points);
    xform.TransformExtent(esriTransformDirection.esriTransformForward, extent);
    xform.TransformCellsize(esriTransformDirection.esriTransformForward, ref xcell,
        ref ycell, extent);
}

Private Sub TransformPoints(ByVal xform As IGeodataXform, ByVal points As IPointCollection, ByVal extent As IEnvelope, ByRef xcell As Double, ByRef ycell As Double)
    xform.TransformPoints(esriTransformDirection.esriTransformForward, points)
    xform.TransformExtent(esriTransformDirection.esriTransformForward, extent)
    xform.TransformCellsize(esriTransformDirection.esriTransformForward, xcell, ycell, extent)
End Sub

IPolynomialXform Polynomial(IPointCollection fromPoints, IPointCollection toPoints)
{
    IPolynomialXform polyXform = new PolynomialXformClass();
    polyXform.DefineFromControlPoints(fromPoints, toPoints, 2);
    return polyXform;
}

Private Function Polynomial(ByVal fromPoints As IPointCollection, ByVal toPoints As IPointCollection) As IPolynomialXform
    Dim polyXform As IPolynomialXform = New PolynomialXformClass()
    polyXform.DefineFromControlPoints(fromPoints, toPoints, 2)
    Return polyXform
End Function

void AdjustXform(IPointCollection fromPoints, IPointCollection toPoints)
{
    //Create AdjustXform using control points.
    IAdjustXform xform = new AdjustXformClass();
    xform.DefineFromControlPoints(fromPoints, toPoints);
    xform.PolynomialApproximation = 2;
    //Use interpolation to perform the adjustment.
    xform.NaturalNeighbor = true;
}

Private Sub AdjustXform(ByVal fromPoints As IPointCollection, ByVal toPoints As IPointCollection)
    'Create AdjustXform using control points.
    Dim xform As IAdjustXform = New AdjustXformClass()
    xform.DefineFromControlPoints(fromPoints, toPoints)
    xform.PolynomialApproximation = 2
    'Use interpolation to perform the adjustment.
    xform.NaturalNeighbor = True
End Sub

void CreateApproximationXform(IGeodataXform xform)
{
    IApproximationXform appXform = new ApproximationXformClass();
    appXform.GeodataXform = xform;
    IGeodataXformApproximation gdApproximation = (IGeodataXformApproximation)
        appXform;
    gdApproximation.Tolerance = 0.5;
    gdApproximation.GridSize = 32;
}

Private Sub CreateApproximationXform(ByVal xform As IGeodataXform)
    Dim appXform As IApproximationXform = New ApproximationXformClass()
    appXform.GeodataXform = xform
    Dim gdApproximation As IGeodataXformApproximation = CType(appXform, IGeodataXformApproximation)
    gdApproximation.Tolerance = 0.5
    gdApproximation.GridSize = 32
End Sub

void CoordinateXform(ISpatialReference sourceSR, ISpatialReference targetSR,
    IPointCollection points)
{
    //Create CoordinateXform and set the spatial reference.
    ICoordinateXform2 coorXform = (ICoordinateXform2)new CoordinateXform();
    coorXform.InputSpatialReference = sourceSR;
    coorXform.SpatialReference = targetSR;
    //Define the approximation parameters.
    IGeodataXformApproximation gdApproximation = (IGeodataXformApproximation)
        coorXform;
    gdApproximation.Approximation = true;
    gdApproximation.GridSize = 16;
    //Ensure unit is in the input space.
    gdApproximation.Tolerance = 15;
    //Transform.
    coorXform.TransformPoints(esriTransformDirection.esriTransformForward, points);
}

Private Sub CoordinateXform(ByVal sourceSR As ISpatialReference, ByVal targetSR As ISpatialReference, ByVal points As IPointCollection)
    'Create CoordinateXform and set the spatial reference.
    Dim coorXform As ICoordinateXform2 = CType(New CoordinateXform(), ICoordinateXform2)
    coorXform.InputSpatialReference = sourceSR
    coorXform.SpatialReference = targetSR
    'Define the approximation parameters.
    Dim gdApproximation As IGeodataXformApproximation = CType(coorXform, IGeodataXformApproximation)
    gdApproximation.Approximation = True
    gdApproximation.GridSize = 16
    'Ensure unit is in the input space.
    gdApproximation.Tolerance = 15
    'Transform.
    coorXform.TransformPoints(esriTransformDirection.esriTransformForward, points)
End Sub

void SetCubicResamplingFactor(IRaster raster)
{
    //Get the xformer.
    IRaster2 raster2 = (IRaster2)raster;
    IRasterXformer xformer = raster2.RasterXformer;
    //Set the resampling method and factor (-3, 0).
    raster.ResampleMethod = rstResamplingTypes.RSP_CubicConvolution;
    ISimplePixelResampler resampler = (ISimplePixelResampler)xformer.PixelResampler;
    resampler.CubicConvolutionParameter =  - 0.5;
}

Private Sub SetCubicResamplingFactor(ByVal raster As IRaster)
    'Get the xformer.
    Dim raster2 As IRaster2 = CType(raster, IRaster2)
    Dim xformer As IRasterXformer = raster2.RasterXformer
    'Set the resampling method and factor (-3, 0).
    raster.ResampleMethod = rstResamplingTypes.RSP_CubicConvolution
    Dim resampler As ISimplePixelResampler = CType(xformer.PixelResampler, ISimplePixelResampler)
    resampler.CubicConvolutionParameter = -0.5
End Sub

IPixelFilter FilterKernel()
{
    IRasterConvolutionFilter filter = new RasterConvolutionFilterClass();
    double[] kernel =
    {
         - 1,  - 1,  - 1, 2, 2, 2,  - 1,  - 1,  - 1
    };
    filter.PutCoefficients(ref kernel);
    return (IPixelFilter)filter;
}

Private Function FilterKernel() As IPixelFilter
    Dim Filter As IRasterConvolutionFilter = New RasterConvolutionFilterClass()
    Dim kernel As Double() = { -1, -1, -1, 2, 2, 2, -1, -1, -1 }
    Filter.PutCoefficients(kernel)
    Return CType(Filter, IPixelFilter)
End Function

void SetPixelFilterOnRaster(IRaster raster)
{
    IStockConvolutionFilter filter = (IStockConvolutionFilter)new
        RasterConvolutionFilterClass();
    filter.Type = esriRasterFilterTypeEnum.esriRasterFilterSmoothing3x3;
    IPixelOperation pixelOp = (IPixelOperation)raster;
    pixelOp.PixelFilter = (IPixelFilter)filter;
}

Private Sub SetPixelFilterOnRaster(ByVal raster As IRaster)
    Dim Filter As IStockConvolutionFilter = CType(New RasterConvolutionFilterClass(), IStockConvolutionFilter)
    Filter.Type = esriRasterFilterTypeEnum.esriRasterFilterSmoothing3x3
    Dim pixelOp As IPixelOperation = CType(raster, IPixelOperation)
    pixelOp.PixelFilter = CType(Filter, IPixelFilter)
End Sub

void MosaicRasterCatalog(IRasterWorkspaceEx ws)
{
    //Open a raster catalog.
    IRasterCatalog rasterCatalog = ws.OpenRasterCatalog("images");

    //Mosaic a raster catalog with a Where clause.
    IMosaicRaster mosaicRaster = new MosaicRasterClass();
    mosaicRaster.RasterCatalog = rasterCatalog;
    mosaicRaster.WhereClause = "state = 'CA'";

    //Save the result.
    ISaveAs saveAs = (ISaveAs)mosaicRaster;
    saveAs.SaveAs("mosaicresult", (IWorkspace)ws, "GDB");
}

Private Sub MosaicRasterCatalog(ByVal ws As IRasterWorkspaceEx)
    'Open a raster catalog.
    Dim rasterCatalog As IRasterCatalog = ws.OpenRasterCatalog("images")

    'Mosaic a raster catalog with a Where clause.
    Dim mosaicRaster As IMosaicRaster = New MosaicRasterClass()
    mosaicRaster.RasterCatalog = rasterCatalog
    mosaicRaster.WhereClause = "state = 'CA'"

    'Save the result.
    Dim saveAs As ISaveAs = CType(mosaicRaster, ISaveAs)
    saveAs.SaveAs("mosaicresult", CType(ws, IWorkspace), "GDB")
End Sub

void Mosaic2(IRasterWorkspaceEx ws, IRaster raster1, IRaster raster2)
{
    IRasterCollection rasterCollection = (IRasterCollection)new MosaicRasterClass();
    rasterCollection.Append(raster1);
    rasterCollection.Append(raster2);

    //Save the result.
    ISaveAs saveAs = (ISaveAs)rasterCollection;
    saveAs.SaveAs("mosaicresult", (IWorkspace)ws, "GDB");
}

Private Sub Mosaic2(ByVal ws As IRasterWorkspaceEx, ByVal raster1 As IRaster, ByVal raster2 As IRaster)
    Dim rasterCollection As IRasterCollection = CType(New MosaicRasterClass (), IRasterCollection)
    rasterCollection.Append (raster1)
    rasterCollection.Append (raster2)

    'Save the result.
    Dim saveAs As ISaveAs = CType(rasterCollection, ISaveAs)
    saveAs.SaveAs("mosaicresult", CType(ws, IWorkspace), "GDB")
End Sub

void RasterCatalogLoader(IPropertySet propertySet)
{
    //Set the storage.
    IRasterStorageDef storage = new RasterStorageDefClass();
    storage.CompressionType = esriRasterCompressionType.esriRasterCompressionLZ77;

    //Set the loader.
    IRasterCatalogLoader loader = new RasterCatalogLoaderClass();
    loader.ConnectionProperties = propertySet;
    loader.StorageDef = storage;
    loader.Projected = false;

    loader.Load("state_catalog", @"\c:\data", null);
}

Private Sub RasterCatalogLoader(ByVal propertySet As IPropertySet)
    'Set the storage.
    Dim storage As IRasterStorageDef = New RasterStorageDefClass()
    storage.CompressionType = esriRasterCompressionType.esriRasterCompressionLZ77

    'Set the loader.
    Dim loader As IRasterCatalogLoader = New RasterCatalogLoaderClass()
    loader.ConnectionProperties = propertySet
    loader.StorageDef = storage
    loader.Projected = False

    loader.Load("state_catalog", "\c:\data", Nothing)
End Sub

static void MeasureBuildHeight(IRaster inputRaster, IPoint firstPoint, IPoint
    secondPoint)
{
    //Define the Mensuration class and set the input.
    IMensuration measure = new MensurationClass();
    measure.Raster = inputRaster;
    //Define the mensurement result.
    IHeightMeasurement heightResult = new HeightMeasurementClass();
    //Measure.
    if (measure.CanMeasureHeight == true)
    {
        measure.GetHeightFromBaseAndTop(firstPoint, secondPoint, out heightResult);
        double height = heightResult.HeightMeasurement;
        double uncertainty = heightResult.HeightUncertainty;
    }
}

Private Shared Sub MeasureBuildHeight(inputRaster As IRaster, firstPoint As IPoint, secondPoint As IPoint)
'Define the Mensuration class and set the input.
Dim measure As IMensuration = New MensurationClass()
measure.Raster = inputRaster
'Define the mensurement result.
Dim heightResult As IHeightMeasurement = New HeightMeasurementClass()
'Measure.
If measure.CanMeasureHeight = True Then
    measure.GetHeightFromBaseAndTop(firstPoint, secondPoint, heightResult)
    Dim height As Double = heightResult.HeightMeasurement
    Dim uncertainty As Double = heightResult.HeightUncertainty
End If
End Sub

void Rotate(IRaster raster)
{
    IRasterGeometryProc rasterGP = new RasterGeometryProcClass();
    rasterGP.Rotate(null, 45, raster);
    rasterGP.Register(raster);
}

Private Sub Rotate(ByVal raster As IRaster)
    Dim rasterGP As IRasterGeometryProc = New RasterGeometryProcClass()
    rasterGP.Rotate(Nothing, 45, raster)
    rasterGP.Register(raster)
End Sub

void ColormapToRGB(IRasterDataset rasterDataset)
{
    IRasterColormapToRGBConversion converter = new RasterColormapToRGBConverterClass
        ();
    IRaster raster = converter.CreateRGBRaster(rasterDataset);
    ISaveAs saveAs = (ISaveAs)raster;
    IWorkspace = (IWorkspace)OpenRasterWorkspace(@"c:\temp");
    saveAs.SaveAs("rgbraster.img", ws, "IMAGINE Image");
}

Private Sub ColormapToRGB(ByVal rasterDataset As IRasterDataset)
    Dim converter As IRasterColormapToRGBConversion = New RasterColormapToRGBConverterClass()
    Dim raster As IRaster = converter.CreateRGBRaster(rasterDataset)
    Dim saveAs As ISaveAs = CType(raster, ISaveAs)
    IWorkspace = CType(OpenRasterWorkspace("c:\temp"), IWorkspace)
    saveAs.SaveAs("rgbraster.img", ws, "IMAGINE Image")
End Sub