Introduction to Programming with DotSpatial

Raster Data Operations in DotSpatial

Tutorial (2)

Purpose of this tutorial: Become familiar with the following raster operations in DotSpatial:

1. Loading different formats of raster data in DotSpatial.

2. Implementing the hill shade.

3. Change the color of the raster data.

4. Multiply the raster data.

5. Reclassify the raster data.

6. Get the mouse clicked point values on the raster data layer.

This tutorial has 5 important steps.

Step 1: Download the DotSpatial class library

Step 2: Add the DotSpatial reference and change the compile option.

Step 3: Add the DotSpatial Controls into the Visual Studio Toolbox.

Step 4: Copy the Data Extensions folder to the debugs folder of the current project

Step 5: Design the GUI. (GUI - Graphical User Interface)

Step 6: Write the code for implementing the map operations.

Step 1: Download the DotSpatial class library

This step is the same as Tutorial #1 step 1.

Step 2: Add the DotSpatial reference and change the compile option.

1.1) Adding the references.

DotSpatial.Data.Forms.dll, DotSpatial.Symbology.dll, DotSpatial.Controls.dll, DotSpatial.Projections.dll,

DotSpatial.Data.dll, DotSpatial.Topology.dll

Teva ~ VeluppillaiPage 1
 Introduction to Programming with DotSpatial

Fig. 1 Required references.

2.2) Change the compile option.

Change the compile option from .Net Framework 4 Client Profile to .Net Framework4.

This step is the same as the Tutorial # 1 step 2.2.

Step 3: Add the DotSpatial Controls into the Visual Studio Toolbox.

This step is the same as the Tutorial # 1 step 3.

Step 4: Copy the Data Extension folder from your downloaded folder to your current project bin/debug
folder. The .dlls from this folder is necessary for GDAL data provider access.

Fig. 2 Data Extensions folder from DotSpatial downloaded unzip folder

Teva ~ VeluppillaiPage 2
 Introduction to Programming with DotSpatial

Step 5: Design the GUI

Design the GUI as follows:

Fig. 3 Final GUI

Interface design considerations.

1. Add three panel controls. Panel control's properties should be as follows:

Properties Panel1 Panel2 Panel3

Name pnlOperation pnlLegend pnlMap

s

Dock Top Left Fill

2. Add five buttons. Button properties should be as follows:

Properties Button1 Button2 Button3 Button4 Button5

Teva ~ VeluppillaiPage 3
 Introduction to Programming with DotSpatial

Name btnLoadRaste btnHillshade btnChangeColor btnMultiplyRaster btnReclassify

r

Text &Load Raster &Hillshade Change &Color &Multiply Raster &Reclassify

Raster

3. Add two label controls. Labels' properties should be as follows:

Properties Label1 Label2

Name lblElevation lblRasterValue

Text Elevation Row: Column: Value:

4. Add a checkbox and it should have the following properties.

Name: chbRasterValue Text : Raster value at clicked point

5. Add a textbox and it should have the following features.

Name: txtElevation Text: 3000

6. Drag a “Legend” control from the DotSpatial tab under toolbox and drop it on pnlLegend. Legend
properties should be as follows:

Name: Legend1 Dock: Fill

7. Drag a “Map” control from the DotSpatial tab under toolbox and drop it on pnlMap. Map properties
should be as follows:

Name: Map1 Dock: Fill Legend: Legend1

8. Drag an "AppManager" control from DotSpatial tab under toolbox and drop it on the form.

Note: This control is necessary for loading different formats of raster data.

Teva ~ VeluppillaiPage 4
 Introduction to Programming with DotSpatial

Fig. 4 AppManager

9. Set the properties of AppManager1 as follows:

Map: Map1 Legend: Legend1

Step: 6 Loading the different formats of raster data.

6.1) Import the following namespaces in the coding window.

//Required namespaces
using DotSpatial.Symbology;
using DotSpatial.Controls;
using DotSpatial.Data;
using DotSpatial.Topology;

6.2 ) btnLoadRaster click event should be as follows:

private void btnLoadRaster_Click(object sender, EventArgs e)

{
//AddRasterLayer() method is used to load the raster layers on the map
map1.AddRasterLayer();
map1.ZoomToMaxExtent();
}

Implementing the hill shade.

btnHillshade_Click event should have the following code.

private void btnHillshade_Click(object sender, EventArgs e)

{

Teva ~ VeluppillaiPage 5
 Introduction to Programming with DotSpatial

if (map1.Layers.Count > 0)
{

//IMapRasterLayer layer = (IMapRasterLayer)map1.Layers[0];

IMapRasterLayer layer =map1.Layers[0] as IMapRasterLayer;

if (layer == null)
{
MessageBox.Show("Please select a raster layer");
return;
}

layer.Symbolizer.ShadedRelief.ElevationFactor = 1;

layer.Symbolizer.ShadedRelief.IsUsed = true;

layer.WriteBitmap();

}
else
{
MessageBox.Show("Please add a layer to the map.");
}
}

Change the color.

btnChangeColor click event should has the following event.

private void btnChangeColor_Click(object sender, EventArgs e)

{
if (map1.Layers.Count > 0)
{
//change the color of raster

//typecast the first layer to MapRasterLayer

IMapRasterLayer layer = map1.Layers[0] as IMapRasterLayer;

if (layer == null)
{
MessageBox.Show("Please add a raster layer.");
return;
}

//set the color scheme

Teva ~ VeluppillaiPage 6
 Introduction to Programming with DotSpatial

//create an instance for a colorscheme

ColorScheme scheme = new ColorScheme();

//create a new category

ColorCategory category1 = new ColorCategory(2500, 3000,
Color.Red, Color.Yellow);

category1.LegendText = "Elevation 2500 - 3000";

//add the category to the color scheme

scheme.AddCategory(category1);

//create another category

ColorCategory category2 = new ColorCategory(1000, 2500,
Color.Blue, Color.Green);

category2.LegendText = "Elevation 1000 - 2500";

scheme.AddCategory(category2);

//assign new color scheme

layer.Symbolizer.Scheme = scheme;

//refresh the layer display in the map

layer.WriteBitmap();

}
else
{
MessageBox.Show("Please add a layer to the map.");
}

Multiply the raster.

btnMultiplyRaster should has the following event.

private void btnMultiplyRaster_Click(object sender, EventArgs e)

{
if (map1.Layers.Count > 0)
{
IMapRasterLayer layer = map1.Layers[0] as IMapRasterLayer ;

Teva ~ VeluppillaiPage 7
 Introduction to Programming with DotSpatial

if (layer == null)
{
MessageBox.Show("Please select a raster layer");
}

IRaster demRaster = layer.DataSet;

string[] rasterOptions = new string[1];

IRaster newRaster = Raster.CreateRaster("multiply.bgd", null,

demRaster.NumColumns, demRaster.NumRows, 1, demRaster.DataType,
rasterOptions);

//Bounds specify the cellsize and the coordinates of raster corner

newRaster.Bounds = demRaster.Bounds.Copy();

newRaster.NoDataValue = demRaster.NoDataValue;

newRaster.Projection = demRaster.Projection;

for (int i = 0; i <= demRaster.NumRows - 1; i++)

{
for (int j = 0; j <= demRaster.NumColumns - 1; j++)
{
if (demRaster.Value[i, j] != demRaster.NoDataValue)
{
newRaster.Value[i, j] = demRaster.Value[i, j] * 2;
}
}
}

//save the new raster to the file

newRaster.Save();

//add the new raster to the map

map1.Layers.Add(newRaster);

else
{
MessageBox.Show("Please add a layer to the map.");
}

Teva ~ VeluppillaiPage 8
 Introduction to Programming with DotSpatial

Reclassify the raster.

btnReclassify should has the following event.

private void btnReclassify_Click(object sender, EventArgs e)

{
//typecast the selected layer to IMapRasterLayer
IMapRasterLayer layer = map1.Layers.SelectedLayer as
IMapRasterLayer ;

if (layer == null)
{
MessageBox.Show("Please select a raster layer.");
}
else
{
//get the raster dataset
IRaster demRaster = layer.DataSet;

//create a new empty raster with same dimension as original raster

string[] rasterOptions = new string[1];

IRaster newRaster = Raster.CreateRaster("reclassify.bgd", null,

demRaster.NumColumns, demRaster.NumRows, 1, demRaster.DataType,
rasterOptions);
newRaster.Bounds = demRaster.Bounds.Copy();
newRaster.NoDataValue = demRaster.NoDataValue;
newRaster.Projection = demRaster.Projection;

//reclassify raster.
// values >= specified value will have new value 1
// values < specified value will have new value 0

double oldValue = 0;

//get the specified value from the textbox

double specifiedValue = Convert.ToDouble(txtElevation.Text);

for (int i = 0; i <= demRaster.NumRows - 1; i++)

{
for (int j = 0; j <= demRaster.NumColumns - 1; j++)
{
//get the value of original raster
oldValue = demRaster.Value[i, j];

Teva ~ VeluppillaiPage 9
 Introduction to Programming with DotSpatial

if (oldValue >= specifiedValue)

{
newRaster.Value[i, j] = 1;
}
else
{
newRaster.Value[i, j] = 0;
}
}
}

newRaster.Save();

map1.Layers.Add(newRaster);
}
}

Get the mouse clicked point values on the raster data layer.

Select the chbRasterValue's "CheckedChanged" event.

fig.5 ChbRasterValue checkbox's checked changed event.

chbRasterValue_CheckedChanged event should has the following events.

private void chbRasterValue_CheckedChanged(object sender, EventArgs e)

{
if (chbRasterValue.Checked)

Teva ~ VeluppillaiPage 10
 Introduction to Programming with DotSpatial

{
IMapRasterLayer rasterLayer = map1.Layers.SelectedLayer as
IMapRasterLayer ;

if ((rasterLayer != null))
{
//set the map cursor to cross
map1.Cursor = Cursors.Cross;
}
else
{
//if no raster layer is selected, uncheck the checkbox
MessageBox.Show("Please select a raster layer.");
chbRasterValue.Checked = false;
}
}
else
{
//change map cursor back to arrow
map1.Cursor = Cursors.Arrow;
}

Select the Map1's mouse up event.

fig.6 Map1_MouseUp event

map1_MouseUp event should have the following code.

Teva ~ VeluppillaiPage 11
 Introduction to Programming with DotSpatial

private void map1_MouseUp(object sender, MouseEventArgs e)

{
if (chbRasterValue.Checked)
{
//get the layer selected in the legend
IMapRasterLayer rasterLayer = map1.Layers.SelectedLayer as
IMapRasterLayer;

if ((rasterLayer != null))
{

//get the raster data object

IRaster raster = rasterLayer.DataSet;

//convert mouse position to map coordinate

Coordinate coord = map1.PixelToProj(e.Location);

//convert map coordinate to raster row and column

RcIndex rc = raster.Bounds.ProjToCell(coord);

int row = rc.Row;

int column = rc.Column;

//check if clicked point is inside of raster

if ((column > 0 & column < raster.NumColumns & row > 0 &
row < raster.NumRows))
{
//get the raster value at row and column
double value = raster.Value[row, column];

//show the row, column and value in the label

lblRasterValue.Text = string.Format("row: {0} column: {1}
value: {2}", row, column, value);
}
else
{
lblRasterValue.Text = "outside of raster";
}

Teva ~ VeluppillaiPage 12
 Introduction to Programming with DotSpatial

Output screen shot for different operations.

Fig.7 Load a raster data

Teva ~ VeluppillaiPage 13
 Introduction to Programming with DotSpatial

Fig.8 After the hill shade operation

Fig.9 After the color change operation

Teva ~ VeluppillaiPage 14
 Introduction to Programming with DotSpatial

Fig.10 After the raster multiplication

Fig.11 After the reclassify operation

Teva ~ VeluppillaiPage 15
 Introduction to Programming with DotSpatial

Fig.12 Getting the mouse click point value on the raster data file.

Teva ~ VeluppillaiPage 16