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); }