public void ClipRasterWithPolygonTest()
{
if (DotSpatial.Data.DataManager.DefaultDataManager.PreferredProviders.Count == 0)
{
DotSpatial.Data.Rasters.GdalExtension.GdalRasterProvider lGdalRasterProvider = new DotSpatial.Data.Rasters.GdalExtension.GdalRasterProvider();
}
String path = AppDomain.CurrentDomain.BaseDirectory;
String shapeFilePath = System.IO.Path.Combine(path, "Data", "elbe_watershed1.shp");
String rasterFilePath = System.IO.Path.Combine(path, "Data", "kriging.bgd" );
String resultFilePath = System.IO.Path.Combine(path, "Data", "clipResult.bgd" );
DotSpatial.Data.Shapefile lClipPolygon = DotSpatial.Data.Shapefile.OpenFile(shapeFilePath);
DotSpatial.Data.IRaster lGridToClip = DotSpatial.Data.Raster.OpenFile(rasterFilePath, false);
DotSpatial.Data.Raster lGridAfterClip = new DotSpatial.Data.Raster();
lGridAfterClip.Filename = resultFilePath;
DotSpatial.Analysis.ClipRaster.ClipRasterWithPolygon(lClipPolygon.Features[0], lGridToClip, lGridAfterClip.Filename);
IRaster ras2 = Raster.Open(lGridAfterClip.Filename);
Assert.AreEqual(lGridToClip.NoDataValue, ras2.NoDataValue);
}
public void ClipRasterWithPolygonTest()
{
if (DotSpatial.Data.DataManager.DefaultDataManager.PreferredProviders.Count == 0)
{
DotSpatial.Data.Rasters.GdalExtension.GdalRasterProvider lGdalRasterProvider = new DotSpatial.Data.Rasters.GdalExtension.GdalRasterProvider();
}
DotSpatial.Data.Shapefile lClipPolygon = DotSpatial.Data.Shapefile.OpenFile(@"C:\Users\Jiri\Desktop\berounka.shp");
DotSpatial.Data.IRaster lGridToClip = DotSpatial.Data.Raster.OpenFile(@"C:\Users\Jiri\Desktop\kriging2.bgd", false);
DotSpatial.Data.Raster lGridAfterClip = new DotSpatial.Data.Raster();
lGridAfterClip.Filename = @"C:\Users\Jiri\Desktop\kriging2.bgd";
DotSpatial.Analysis.ClipRaster.ClipRasterWithPolygon(lClipPolygon.Features[0], lGridToClip, lGridAfterClip.Filename);
IRaster ras2 = Raster.Open(lGridAfterClip.Filename);
Assert.AreEqual(lGridAfterClip.NoDataValue, ras2.NoDataValue);
}
public void SaveAsTest()
{
if (DotSpatial.Data.DataManager.DefaultDataManager.PreferredProviders.Count == 0)
{
DotSpatial.Data.Rasters.GdalExtension.GdalRasterProvider lGdalRasterProvider = new DotSpatial.Data.Rasters.GdalExtension.GdalRasterProvider();
}
string GridDataFolder = System.IO.Path.GetDirectoryName(System.Reflection.Assembly.GetExecutingAssembly().Location) + @"\..\..\..\Data\Grids\";
DotSpatial.Data.IRaster sourceGrid = Raster.Open(GridDataFolder + @"elev_cm_ESRI\elev_cm_clip2\hdr.adf");
Double sourceGridMaximum = sourceGrid.Maximum;
string savedGridName = GridDataFolder + @"elev_cm.tif";
sourceGrid.SaveAs(savedGridName);
Assert.AreEqual(sourceGrid.Maximum, sourceGridMaximum, 0.0001);
DotSpatial.Data.IRaster savedSourceGrid = Raster.Open(savedGridName);
Assert.AreEqual(sourceGridMaximum, savedSourceGrid.Maximum, 0.0001);
sourceGrid.Close();
savedSourceGrid.Close();
System.IO.File.Delete(savedGridName);
}
