DataSourcesRaster

• Engine
• ArcGIS for Desktop Basic
• ArcGIS for Desktop Standard
• ArcGIS for Desktop Advanced
• Server

• ESRI.ArcGIS.Geodatabase
• ESRI.ArcGIS.DataSourcesRaster

• Classes that access the three raster data models—raster datasets, raster catalogs, and mosaic datasets—from various data sources including file systems, personal geodatabases, file geodatabases, and ArcSDE geodatabases
• Classes that perform geometric transformations including polynomial transformation, Rational Polynomial Coefficient (RPC) camera model, frame camera model, and so on, as well as classes that perform image registration
• Classes used for on-the-fly raster processing, such as creating shaded relief, pan-sharpening, and arithmetic operations
• Classes used for creating mosaics from data of various raster formats and multiple satellite sensors
• Classes used for loading raster datasets and raster catalogs in geodatabases, and other miscellaneous classes

• GeodataXform
• RasterFunction
• RasterFunctionArguments
• RasterBuilder
• Crawler
• PixelResampler
• PixelReader
• MosaicOperator

• Raster data access classes
• Function raster datasets and raster functions
• Mosaic datasets and raster types
• Geodata transformation and pixel filtering
• Image registration objects
• Image mensuration objects
• Raster data loading and other miscellaneous objects

Raster data access classes

• Raster dataset
• Raster catalog
• Mosaic dataset

Workspace

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

Raster datasets

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

Raster bands

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

Raster attribute table, raster color map, and raster statistics

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

Raster

Creating a raster

• The IRasterDataset2.CreateFullRaster method creates a raster with all the properties from the raster dataset, such as number of bands and width and height of the raster dataset.
• The IRasterDataset.CreateDefaultRaster method creates a raster that has a square cell size and contains only three raster bands even if the dataset has more than three bands. The three bands are the default bands used in the raster RGB renderer and are determined by the settings for default raster behavior made on the RasterDefaultsEnv class in the Carto library.

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

Setting raster properties

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

Setting transformations on a raster

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

Editing and displaying a raster

Pixel blocks

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

NoData

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

Raster catalogs

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

Function raster datasets and raster functions

Raster function

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

Function raster dataset

//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

Mosaic datasets and raster types

Mosaic dataset

• A mosaic catalog, which is an internal raster catalog for managing rasters, raster functions, and footprints
• A polygon feature class defining the boundary of the mosaic
• A raster type table that stores the raster types used in the mosaic
• A log table containing the operations performed on the mosaic dataset
• A stereo table providing stereo pair information (optional)
• A color correction table storing the derived local mean and local standard deviation images (optional)
• A seamline table that defines seamlines for mosaicking (optional)

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 created above.
    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 created above.
    Dim theMosaicDataset As IMosaicDataset=mosaicExt.CreateMosaicDataset(mdName, mosaicSrs, creationPars)
    '
    Return theMosaicDataset
End Function

Mosaic function

• The selection of rasters that will participate in the mosaic
• The mosaic method (for example, how the selected rasters are ordered)
• The mosaic operation (for example, how overlapping pixels at the same location are resolved)

Raster type

• The raster data type to be added (for example, regular raster datasets or data from sensors with special metadata information, so that the crawler understands what files to search for)
• The processes to be used in creating the rasters (for example, multiband rasters, single band rasters, or pan-sharpened rasters)

Geodata transformation and pixel filtering

A geodata transformation is a mathematical operation that is used to perform geometric transformations for geodatasets such as raster datasets. ArcGIS supports many geodata transformations, such as coordinate transformation (reprojection), polynomial transformation, rubber sheeting (Adjust) transformation, Spline transformation, and RPC transformation.

Geodata transformations are represented as an abstract class, GeodataXform, and as many concrete geodata transformation classes (Xform) as follows:

• AdjustXform
• SplineXform
• RPCXform
• PolynomialXform
• CoordinateXform
• ApproximationXform
• GCSShiftXform
• GeometricXform
• CompositeXform
• FrameXform
• SPOT5XForm
• LSRXform
• ProjectiveXform
• CSMXForm
• Formosat2Xform
• IdentityXform

You can also create a custom Xform by implementing interfaces of the GeodataXform abstract class. The geodata transformation objects are shown in the following diagram:

All GeodataXform objects can have an output spatial reference that defines the spatial reference of the data after transformation (for example, in output space). The domains of GeodataXform objects, if they exist, are also defined in the output space. All GeodataXform objects perform point (IPointCollection or WKSPoint) transformation, extent (IEnvelope) transformation, and cell size transformation in forward and backward directions. To transform raster data, the RasterXformer object (discussed at the end of this section) provides the behind-the-scenes mechanism by calling PixelReader, PixelResampler, and IRasterXform. The IRasterXform interface is used for geodata transformations that need raster-specific information to optimize operations. Therefore, when creating a custom geodata transformation to work with raster data, implement the IRasterXform interface.

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

PolynomialXform

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

SplineXform

AdjustXform

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

ApproximateXform

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

CoordinateXform

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

RPCXform

Other GeodataXform classes

Raster transformation

• Nearest neighbor
• Bilinear
• Cubic convolution
• Majority

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

Pixel filters

• RasterConvolutionFilter
• PansharpeningFilter
• BackgroundFilter
• ColormapFilter
• LutFilter
• MultibandFilter
• RemapFilter

• Brovey
• IHS (intensity, hue, saturation)
• Esri
• Mean

Brovey method

IHS method

Esri method

Mean method

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

Image registration objects

Image-to-image matching

void ImageMatchingExample(IRaster ras, IRaster ras_ref)
{
    //Compute the matching points from the overlapped images.
    IImageRegistration imageRegistration=new ImageRegistrationClass();
    IImageMatchPoints matchPoints=imageRegistration.ExtractMatchPoints(ras,
        ras_ref);

    //Prepare source and target points from the matching points.
    IPointCollection sourcePoints=new MultipointClass();
    IPointCollection targetPoints=new MultipointClass();
    for (int i=0; i < matchPoints.Count; i++)
    {
        IImageMatchPoint matchPoint=matchPoints.get_Element(i);
        IImageFeaturePoint imageFromPoint=matchPoint.From;
        IImageFeaturePoint imageToPoint=matchPoint.To;
        sourcePoints.AddPoint(imageFromPoint.Location);
        targetPoints.AddPoint(imageToPoint.Location);
    }

    //Use RasterGeometryProc to perform georeferencing using source and target control points.
    IRasterGeometryProc rasProc=new RasterGeometryProcClass();
    rasProc.Warp(sourcePoints, targetPoints,
        esriGeoTransTypeEnum.esriGeoTransPolyOrder2, ras);
    rasProc.Register(ras);
}

Private Sub ImageMatchingExample(ras As IRaster, ras_ref As IRaster)
    'Compute the matching points from the overlapped images.
    Dim imageRegistration As IImageRegistration=New ImageRegistrationClass()
    Dim matchPoints As IImageMatchPoints=imageRegistration.ExtractMatchPoints(ras, ras_ref)
    
    'Prepare source and target points from the matching points.
    Dim sourcePoints As IPointCollection=New MultipointClass()
    Dim targetPoints As IPointCollection=New MultipointClass()
    For i As Integer=0 To matchPoints.Count - 1
        Dim matchPoint As IImageMatchPoint=matchPoints.get_Element(i)
        Dim imageFromPoint As IImageFeaturePoint=matchPoint.From
        Dim imageToPoint As IImageFeaturePoint=matchPoint.[To]
        sourcePoints.AddPoint(imageFromPoint.Location)
        targetPoints.AddPoint(imageToPoint.Location)
    Next
    
    'Use RasterGeometryProc to perform georeferencing using source and target control points.
    Dim rasProc As IRasterGeometryProc=New RasterGeometryProcClass()
    rasProc.Warp(sourcePoints, targetPoints, esriGeoTransTypeEnum.esriGeoTransPolyOrder2, ras)
    rasProc.Register(ras)
End Sub

Image registration of a collection of images

1. Create the feature point table from the input mosaic dataset.
2. create tie points from the feature point table.
3. Compute links from the tie points.
4. Compute adjustments.
5. Apply the solution to the mosaic dataset if the estimated error is acceptable.

void ImageRegistrationExample(IWorkspace ws, IMosaicDataset md)
{
    //Get the spatial reference for creating control points and feature points.
    IGeoDataset md_geo=(IGeoDataset)md;
    ISpatialReference sr=md_geo.SpatialReference;

    //Create the ImageRegistration object.
    IImageRegistration imageRegistration=new ImageRegistrationClass();

    //Prepare the feature point table for feature point extraction.
    ITable featurePoint=imageRegistration.CreateFeaturePointTable(ws, 
        "featurepointTable", sr, "default");

    //Prepare the control point table for calculating control points.
    ITable tiePoint=imageRegistration.CreateControlPointTable(ws, 
        "is_controlpoint", sr, "default");

    //Prepare the solution table for cauclating adjustment.
    ITable solutionTable=imageRegistration.CreateSolutionTable(ws, 
        "ir_solutiontable", sr, "default");

    //Compute feature points from the input mosaic dataset.
    imageRegistration.ComputeFeaturePoints(md, null, null, null, featurePoint, null);

    //Compute tie points from the feature point table.
    imageRegistration.ComputeTiePoints(featurePoint, null, null, null, tiePoint,
        null);

    //Compute links from the tie point.
    imageRegistration.ComputeLinks(tiePoint, null, solutionTable, null);

    //Compute adjustment from the mosaic dataset using the computed links.
    imageRegistration.ComputeAdjustments(md, null, solutionTable, 2, null);

    //If the RMS is acceptable, apply the adjustment to the mosaic dataset.
    imageRegistration.Adjust(md, null,
        esriImageAdjustmentType.esriImageAdjustmentAppend, solutionTable, null);
}

void ImageRegistrationExample(IWorkspace ws, IMosaicDataset md)
{

Image mensuration objects

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 measurement result.
    IHeightMeasurement heightResult=new HeightMeasurementClass();
    //Measure.
    if (measure.CanMeasureHeightBaseToTop == 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 measurement result.
Dim heightResult As IHeightMeasurement=New HeightMeasurementClass()
'Measure.
If measure.CanMeasureHeightBaseToTop=True Then
    measure.GetHeightFromBaseAndTop(firstPoint, secondPoint, heightResult)
    Dim height As Double=heightResult.HeightMeasurement
    Dim uncertainty As Double=heightResult.HeightUncertainty
End If
End Sub

Raster data loading and other miscellaneous objects

Raster mosaicking

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

Raster catalog loader

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

RasterGeometryProc

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

Color map to RGB conversion

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

Other miscellaneous classes