PGIS PRACTICALS

MANUAL

N.G. ACHARYA D.K Marathe college

CHEMBUR TROMBAY EDUCATION SOCIETY
 Index
Sr Practical Aim Pg.No
No. No.
1 1 Creating and Managing Vector Data
1A) Adding vector layer 2

1B) Setting properties

7
1C) Vector Layer Formatting 9
2 2 Calculating line lengths and statistics 10

3 3 3A) Adding raster layers 13

3B) Raster Styling and Analysis 14

4 4 Raster Mosaicking and Clipping 16
5 5 Making a Map 20
6 6 6A) Importing Spreadsheets or CSV files 24

6B) Using Plugin 25

7 7 Working with attributes 26
8 8 Terrain Data and Hill shade analysis 28
9 9 9A) Table Join 33
37
9B) Spatial Join
10 10 10A)Points In Polygon 39

10B)Performing spatial queries 42

11 11 Nearest Neighbor Analysis 46

12 12 Sampling Raster Data using Points or Polygons 56
13 13 Interpolating Point Data 62

1
 PRACTICAL NO: 1

AIM : - Creating and Managing Vector Data:

a) Adding vector layer
b) Setting properties
c) Vector Layer Formatting
Procedure:

a. Adding vector layers (Polygon, Line, Points)

➢ Polygon layers (We have taken 2 layers Matunga, Garden)
➢ Line layers (We have taken 3 layers Small_Roads, Road, Flyover)
➢ Point layers (We have taken 4 layers bank,college,Restaurants,ATM)
b. Setting properties (Labeling, Symbolism)

Our aim is to create map representing a location and its surrounding as follows:

2
1A) Creating Polygon vector layer
➢ Select Project→New
➢ Select Layer→Create Layer→New Shapefile Layer
➢ Select Polygon option from Geometry type.

➢ Fill the appropriate information in each text box.

• File name :
▪ By default the file will be saved in bin folder.
▪ To avoid it click on following button to change the location of file.
➢ Field Panel
➢ Add the Attribute you want to show. (Column Name for Table)
b. Specify Type (DataType:Text Data/Decimal Data/Whole Number/Date) of Attribute
c. Specify the Length of the Attribute. Specify Precision (If Data Type is Decimal)
➢ Click on Add to Field List Button.
➢ You can add as many fields (Column Name) as you want for the layer.
➢ Select Geometry Type as follows
• Click on the following button
➢ The CRS dialog box will appear on screen. Click on the WGS84 option and it will be selected
as follows. click on OK
a) Follow the steps to plot Polygon features.
➢ Select the Polygon Feature( In our case it is chembur for background) from layer panel

3
➢ Click Toggle Editing Button → Click on Add Polygon →Now place the cursor
at the location where you want to place the polygon. for polygon layer minimum 3 points
should be selected

➢ Save the newly added polygon as follows.

➢ Set style for polygon by using property window(Right click on Matunga Layer)
➢ Following screen will appear on the screen. Select pattern as you want and click on OK.

4
➢ Same way we can add one more polygon layer for Gardens.

5
6
1b Creating Line vector layer
➢ Repeat the same steps as we have done for polygon layer.
➢ Select geometry type Line.
➢ Road layer :
➢ To plot road click on Add Line Feature.

➢ Click on the map where you want to draw line.

➢ Once you are done then right click on map (Dotted line turn into solid line)
➢ save your data
➢ set style for Roads in the same way as we have done for polygon

➢ set style for Roads in the same way as we have done for polygon

➢ Road will look as below

7
➢ To label your roads Right click on Road layer .Go to properties window then select label and
set single label property
➢ Following window will appear on the screen

➢ Roads will look like these

➢ To merge roads
• Go to properties of road then select symbology. Click on Advanced button select
Symbol levels.

➢ Check Enable symbol levels option

➢ Click ok & Road will appear as follows

8
1C. Create Point vector layer
➢ Repeat same steps to add point layers as we have done in previous layers.(For
ATM, Restaurants, Banks, Bus Stops etc)

Final output:

9
 PRACTICAL NO.2

AIM: Calculating line lengths and statistics

➢ Go to Layer → Add Layer → Add Vector Layer
➢ Add the following file to project
"\GIS_Workshop\Practicals\Practical_01\D\DATA\IND_rrd\IND_rails.shp"
Press “ADD”
➢ Also add India Administrative Map
“GIS_Workshop\Practicals\Practical_01\D\DATA\IND_adm\IND_adm0.shp”
➢ Double Click on IND_adm0

Select SYMBOLOGY→ Select any outline style from below given options.

Press OK
➢ The display window will appear like

10
➢ In Layer Pane, Right click on IND_rails → Open Attribute Table

➢ Press Toggle Editing button using button , on Attribute table window toolbar.

➢ Press Open Field Calculator using button .

➢ Set the output field as “Track_Len”, field type to “Decimal Number”.

➢ From Function List search $length or go to Geometry → Select $length

➢ Set expression as

Press “OK”
➢ A new column is added to the attribute table with value representing the length of track in
KM.

➢ Press CTRL+S or click on Save Edits option on tool bar

➢ Close the attribute table window.
➢ For calculating the total length of Railway tracks in India.
➢ Select Vector→ Analysis Tools→ Basic Statics for Fields
➢ Select IND_rails layer from input layer. And select Track_Len in “Field to Calculate statistics
on”
➢ Press RUN

11
➢ Open the “output.html” file to get the field statistics.

Analyzed field: Track_Len

Count: 2012
Unique values: 1608
NULL (missing) values: 0
Minimum value: 0.0
Maximum value: 400.48
Range: 400.48
Sum: 60479.320000000014
Mean value: 30.059304174950306
Median value: 14.04
Standard deviation: 39.483220276624444
Coefficient of Variation: 1.313510786770889
Minority (rarest occurring value): 0.03
Majority (most frequently occurring value): 0.0
First quartile: 3.35
Third quartile: 42.855000000000004
Interquartile Range (IQR): 39.505

➢ The above statistics show that the total length of Railway track in India is 60,479.32 KM.

12
 PRACTICAL NO: 3
Exploring and Managing Raster data:
AIM: Adding raster layers
➢ From menu bar select Layer → Add Layer → Add Raster Layer

➢ Select Gridded Population of the World (GPW) v3 dataset from Columbia University,
Population Density Grid for the entire globe in ASCII format and for the year 1990 and 2000.
“\GIS_Workshop\Practicals\Practical_02\A\Data\gl_gpwv3_pdens_90_ascii_one\
glds90ag60.asc”
“\GIS_Workshop\Practicals\Practical_02\A\Data\gl_gpwv3_pdens_90_ascii_one\
glds00ag60.asc”

➢ Go to Project → Properties OR Press the Set CRS option on bottom

right corner.
Select WGS 84 EPSG: 4326 and Press OK

13
3B)AIM:Raster Styling and Analysis
➢ To start with analysis of population data, convert the pixel from grayscale to Color.
➢ Select “glds90ag60.asc” Layer form layer Pane → select property OR double click on it.

➢ Select
➢ Press “APPLY”
➢ Repeat the same for “glds00ag60.asc” Layer

Layer output after applying style.

➢ The objective this experiment is to analyze raster data, as an example we will find areas with
largest population change between 1990 and 2000, by calculating the difference between each
pixel values.
➢ Go to Raster → Raster Calculator

➢ Put the expression "glds00ag60@1" - "glds90ag60@1"

➢ Select the output file location & name and Press OK.
➢ Remove the other two layers i.e. glds00ag60.asc and glds90ag60.asc
➢ Double click on pop_diff layer.

➢ Select

➢ Set Render Type to “Single band Pseudo color”, Interpolation as Discrete, and remove all
classification and add as shown in figure above using button. After all settings press

14
“OK”.
➢ Layer will appear like

➢ Explore an area of your choice and check the raster band value using to verify the
classification rule.
➢ The red pixel shows negative changes and blue shows positive changes.

15
 PRACTICAL N0.4
AIM: Raster Mosaicking and Clipping
A mosaic is a combination or merge of two or more images.
In GIS, a single raster dataset can be created from multiple raster datasets by mosaicking them
together.

In many cases, there will be some overlap of the raster dataset edges that are being mosaicked
together, as shown below.

These overlapping areas can be handled in several ways; for example, you can choose to only
keep

raster data from the first or last dataset, you can blend the overlapping cell values using a
weight-
based algorithm, you can take the mean of the overlapping cell values, or you can take the

minimum or maximum value. When mosaicking discrete data, the First, Minimum, or Maximum
options give the most meaningful results. The Blend and Mean options are best suited for
continuous data. If any of the input rasters are floating point, the output is floating point. If all
the
inputs are integer and First, Minimum, or Maximum is used, the output is integer.
➢ Go to Layer → Add Layer → Add Raster Layer.
➢ Select the following “.tif” raster images for India from data folder.
India-01.tif
India-02.tif
India-03.tif
India-04.tif
India-05.tif

16
➢ Press open
➢ In data source manager | Raster window click Add.

➢ Go to Raster → Miscellaneous → Merge

➢ In the Merge dialog window

➢ Select all layers and Press OK.

17
➢ In Merge dialog window select a file name and location to save merged images.
➢ Save the file to “GIS_Workshop/Practicals/Practical_02/C/” location with the name as
Merge_Files.tif
➢ Press Run and after completion of operation close the Merge window dialog box.

➢ You can now deselect individual layers from layer pane and only keep the merged raster file.

➢ Go to Layer→ Add Vector Layer → Select

\GIS_Workshop\Practicals\Practical_02\C\IndiaAdminBoundry\IND_adm0.shp file.

➢ From layer properties → select → select any one of the following

18
➢ The result will be

➢ Go to Raster → Extraction → Clip Raster by Mask Layer

➢ Select the merge raster image as input and Ind_adm0 as mask layer.
➢ Select a file name and location for clipped raster as /Practical_02/C/Clipped_File.tif.
➢ Press RUN.
AFTER CLIPPING

19
 PRACTICAL NO.5

AIM: Making a Map

➢ Create a new Thematic Map or open and existing one
➢ We will take map of Maharashtra and highlight Mumbai in it
➢ Go to Project → New PrintLayout
➢ Insert a suitable title and press “OK”.

➢ A new Print Layout window will open

➢ Select Add Item → Add Map

20
➢ After adding map go to ItemProperties → Map1→ Layers
Check on Lock Layers and Lock Styles for Layers

This will ensure that if any change in layers or change their styles, the Print Layout view will not
change.
➢ Go to Add Item → Add Picture → Place a picture box at appropriate location.

➢ Also adjust Image Rotation to its appropriate value.

➢ Item Properties → Image Rotation

➢ Add an inset Using Add Item → Add Picture → Select an area to be highlighted on main Map.
➢ Set a frame for Inset by enabling the check box for Frame.

➢ To highlight the area shown in Inset

➢ Select the Picture representing main Map from Items pane.

➢ In Item Properties → Overviews → using icon add an overview.

➢ Select the checkbox Draw Overview
➢ Name the Picture object representing inset (Map1 in our case).

21
➢ The Print Layout will appear like

➢ Add Item → Add Label

➢ Change the Label text To “Mumbai Map”, Set appropriate font size and color using Item
Properties→ Main Properties.

➢ Add Item → Add Legend→ Place the legend indicator at appropriate location.

22
➢ Uncheck auto update and use suitable legend indicator label.
➢ The Print Layout will appear

➢ Add Item → Add Scale Bar

➢ Add Item → Add Label→Add a Label using HTML rendering

➢ A Map can be saved in Image or PDF using Layout → Export as Image / Export as PDF

➢ Save the Map to a location appropriate location as PDF or Image.

➢ Open the PDF or Image from location.

23
 PRACTICAL NO.6
6A)AIM: Importing Spreadsheets or CSV files
➢ Many times the GIS data comes in a table or an Excel spreadsheet or a list lat/long
coordinates,
therefore it has to be imported in a GIS project.
➢ Sample file for Earthquake data will be used in this practical.
➢ Go to Layer → Add Layer → Add Delimited text Layer

➢ Data Source Manager | Delimited Text window will appear

➢ Select the \GIS_Workshop\Practicals\Practical_03\C\Sample.csv file from data folder.

➢ Press ADD and close the window.

➢ Output:

24
6B) AIM: Using Plugins
➢ Core plugins are already part of the standard QGIS installation. To use these, just enable
them.
➢ Open QGIS. Click on Plugins → Manage and Install Plugins....

➢ To enable a plugin, check on the checkbox next to Plugin. This will enable the plugin to use it.
➢ External plugins are available in the QGIS Plugins Repository and need to be installed by the
users
before using them.
➢ Click on Not Installed or Install from ZIP.
➢ Once the plugin is downloaded and installed, you will see a confirmation dialog.
➢ Click on Plugins → <<new Plugin Name>>
➢ The Plugin if marked Experimental plugin can be installed, from Setting→ check on

or

➢ A tab will be added to Plugin Manager Window.

➢ Click on a plugin name and Click Install.

25
 PRACTICAL NO.7

AIM: Working with attributes

➢ Start a new project.
➢ Go to Layer → Add Layer → Add Vector Layer
➢ Select “\GIS_Workshop\Practicals\Practical_04\A\Data\
ne_10m_populated_places_simple.zip”

➢ Right click on Layer in Layer Panel → Open Attribute Table.

➢ Explore various attributes and their values in the Attribute table.
➢ To find the Place with maximum population click on “pop_max” file
➢ On clicking the Select feature using expression button the following window will appear.

➢ Enter pop_max>100 and pop_max<10000 and click button SELECT FEATURES to get all the
places with population between 100 and 10000.

26
➢ Different queries can be performed using the dataset.
➢ Try this

Will give-

PRACTICAL NO.8

27
AIM: Terrain Data and Hill shade analysis
A terrain dataset is a multiresolution, TIN-based surface built from measurements stored as
features
in a geodatabase. Terrain or elevation data is useful for many GIS Analysis like, to generate
various
products from elevation data such as contours, hillshade etc.

➢ Go to Layer → Add Raster Layer → select “10n060e_20101117_gmted_mea300.tif”, from

Data folder

28
➢ The Lower altitude regions are shown using dark color and higher using light shade as seen
on
top region containing Himalaya and Mt Everest.
➢ Mt. Everest - is located at the coordinates 27.9881° N, 86.9253° E.
➢ Enter 86.92, 27.98 in the coordinate field, Scale 900000 and Magnifier 100% at the bottom of
QGIS.

➢ Press enter the view port will be centered on Himalaya Region.

➢ Crop the raster layer only for the region under study.

29
➢ Go to Raster → Extraction→ Clip Raster by Extent
➢ Select the raster layer (if project contains multiple layers).
➢ Select the clipping area by selecting the option Use Canvas Extends if the visible part of map
is to be selected or manually select an area on canvas by using Select Extent on Canvas.
➢ Select the location and file name for storing clipped raster layer.

➢ Press RUN.
➢ Deselect the original layer and keep the clipped one.
➢ The Clipped raster layer is representing altitude are from 103 Meters.

Original Raster Clipped Raster

➢ Counter lines are the lines on a map joining points of equal height above or below sea level.
A
contour interval in surveying is the vertical distance or the difference in the elevation between
the two contour lines in a topographical map.
➢ To derive counter lines from given raster.
➢ Go to Raster → Extraction→ Contour

➢ The Contour configuration window will appear

30
➢ Select the input raster layer name. Set contour interval 100.00 meters, select the output file
name & location and check the option to add output file to project after processing.
➢ Press “RUN”.
➢ The contour layer will appear like this

➢ Label the layer using “ELEV” field and set appropriate symbols for line.
➢ In the Layer panel right click on Contour Raster Layer and select “Open Attribute table”,
➢ Arrange the table in descending order based on the value of “ELEV” column.

31
Compare the above counter line raster layer with the previous Google map image or visit
https://www.google.com/maps/@27.9857765,86.9285378,14.75z/data=!5m1!1e4?hl=en-US
➢ To verify the above contour files using Google Map
➢ Make a copy of Contour Layer, Go to Layer →Save As
➢ Select file format as “Keyhole Markup Language”, set file name, location and Layer Name.
➢ Also set CRS to WGS 84 EPSG:4326

32
 PRACTICAL NO.9
Managing Data Tables and Saptial data Sets:
9A) AIM :Table joins
➢ Start a new project
➢ Go to Layer → Add Layer → Add new Vector Layer
“I:\GIS_Workshop\Practicals\Practical_07\A\Data\tl_2013_06_tract.zip”
➢ We could import this csv file without any further action and it would be imported. But, the
default
type of each column would be a String (text). That is ok except for the D001 field which contains
numbers for the population. Having those imported as text would not allow us to run any
mathematical operations on this column. To tell QGIS to import the field as a number, we need
to create a sidecar file with a .csvt extension.

➢ This file will have only 1 row specifying data types for each column. Save this file as
ca_tracts_pop.csvt in the same directory as the original .csv file.
➢ Go to Layer → Add Layer → Add Delimited Text Layer
And add I:\GIS_Workshop\Practicals\Practical_07\A\Data\ca_tacts_pop.csv”

33
➢ In the layer panel, Right click on “tl_2013_06_tract”, layer and select Properties

➢ Select the option + in Properties, and click on button ADD to add new table join.
➢ In the Add Vector Join window set the following properties and click OK.

34
➢ After performing join

➢ For more clear output, select “tl_2013_06_tact” from Layer Panel, right click and select
properties. Go to Symbology and set the following properties.

35
➢ A detailed and accurate population map of California can be seen as the result. Same
technique
can be used to create maps based on variety of census data.

36
9B) AIM :Spatial Joins
➢ Go to Layer → Add Layer → Add Vector Layer → Select
“I:\GIS_Workshop\Practicals\Practical_07\B\Data\nybb_12c\nybb_13c_av\nybb.shp” and
“I:\GIS_Workshop\Practicals\Practical_07\B\Data\OEM_NursingHomes_001\OEM_NursingHo
mes_001.shp”, from data folder.

➢ Go to attribute table and observe the data.

➢ Table before performing Join
➢ Go to Vector → Data Management Tools → Join Attributes by Location

37
➢ Attribute table after join

➢ Use the Identify Feature Button to select a region to view join data on map Layer.
➢ Output

38
 PRACTICAL NO 10

10A) AIM :Points in polygon analysis

➢ Go to Layer → Add Layer → Add Delimited Text Layer
Select “EarthQuakeDatabase.txt”

Go to Layer → Add Layer → Add vector Layer

“I:\GIS_Workshop\Practicals\Practical_07\C\Data\ne_10m_admin_0_countries.zip”

39
➢ Use the select Feature button to check country wise counting of Earthquakes.

40
41
10B) AIM :Performing spatial queries
➢ Go to Layer → Add Layer → Add Vector Layer and load
“\GIS_Workshop\Practicals\Practical_07\D\Data\ne_10m_populated_places_simple\
ne_10m_popul
ated_places_simple.shp” and
“I:\GIS_Workshop\Practicals\Practical_07\D\Data\ne_10m_rivers_lake_centerlines\
ne_10m_rivers
_lake_centerlines.shp” from project data folder.

➢ Open project Properties → Set CRS “World_Azimuthal_Equidistant EPSG 54032” . The map
will
be re-projected

42
➢ Go to Vector → Geoprocessing Tool → Buffer

➢ Repeat the step to create River Buffer

43
➢ Create a buffer for River

➢ Go to Vector → Research Tool → Select By Location

44
➢ This will highlight only those rivers containing a populated place within 2 KM

45
 PRACTICAL NO.11
Advanced GIS Operations 1:
AIM :Nearest Neighbor Analysis
➢ Go to Layer → add Layer → add Delimited Text Layer and load “EarthQuackeDataBase.txt”
from data file.

➢ Go to Layer → Add Layer → Add vector Layer and from data folder
“\GIS_Workshop\Practicals\Practical_08\A\DATA\ne_10m_populated_places_simple.zip” load
the layer to the project and remove all rows from attribute table other than India.

46
➢ Go to Vector→ Analysis tool → Distance Matrix

➢ Calculate the Distance matrix and perform Nearest Neighbor Analysis

➢ Now you will be able to see the content of our results. The InputID field contains the field
name
from the Earthquake layer. The TargetID field contains the name of the feature from the

47
Populated Places layer that was the closest to the earthquake point. The Distance field is the
distance between the 2 points.

Here select the earthquake layer as the Input point layer and the populated
places ne_10m_populated_places_simpleas the target layer. You also need to select a unique
field from each of these layers which is how your results will be displayed. In this analysis, we
are looking to get only 1 nearest point, so check the Use only the nearest(k) target points, and
enter 1. Name your output file matrix.csv, and click OK. Once the processing finishes,
click Close.

Once the processing finishes, click the Close button in the Distance Matrix dialog. You can now
view the matrix.csv file in Notepad or any text editor. QGIS can import CSV files as well, so
we will add it to QGIS and view it there. Go to Layer ‣Add Layer ‣Add Delimited Text
Layer....

48
9. Browse to the newly created matrix.csv file. Since this file is just text columns, select No
geometry (attribute only table) as theGeometry definition. Click OK.
10. You will see the CSV file loaded as a table. Right-click on the table layer and select Open
Attribute Table.
11. Now you will be able to see the content of our results. The InputID field contains the field
name
from the Earthquake layer. The TargetID field contains the name of the feature from the
Populated Places layer that was the closest to the earthquake point. The Distance field is the
distance between the 2 points.
12. This is very close to the result we were looking for. For some users, this table would be
sufficient. However, we can also integrate this results in our original Earthquake layer using
a Table Join. Right-click on the Earthquake layer, and select Properties.
13. Go to the Joins tab and click on the + button.

14. We want to join the data from our analysis result to this layer. We need to select a field
from each of the layers that has the same values. Select matrix as the Join layer` and InputID as
the Join field. The Target field would be I_D. Leave other options to their default values and
click OK.
15. You will see the join appear in the Joins tab. Click OK.

49
16. Now open the attribute table of the signif layer by right-clicking and selecting Open
Attribute Table.

17. You will see that for every Earthquake feature, we now have an attribute which is the
nearest neighbor (closest populated place) and the distance to the nearest neighbor.
18. We will now explore a way to visualize these results. First, we need to make the table join
permanent by saving it to a new layer. Right-click the signif layer and select Save As....x
19. Click the Browse button next to Save as label and name the output layer
as earthquake_with_places.shp. Make sure the Add saved file to map box is checked and click
ok

50
20. Once the new layer is loaded, you can turn off the visibility of the signif layer. As our dataset
is quite large, we can run our visualization analysis on a subset of the data. QGIS has a neat
feature where you can load a subset of features from a layer without having to export it to a
new layer. Right-click the earthquake_with_places layer and select Properties.\

21. In the source tab, scroll down to the Feature subset section. Click Query Builder.

51
22. For this tutorial, we will visualize the earthquakes and their nearest populated places for
Mexico. Enter the following expression in the Query Builder dialog.

23. You will see that only the points falling within Mexico will be visible in the canvas. Let’s do
the same for the populated places layer. Right-click on

52
the ne_10m_populated_places_simple layer and select Properties.

24. Open the Query Builder dialog from the General tab. Enter the following expression.
"adm0name" = 'Mexico'

25. Now we are ready to create our visualization. We will use a plugin named MMQGIS. Find
and install the plugin. See Using Pluginsfor more details on how to work with plugins. Once you
have the plugin installed, go to MMQGIS ‣Create ‣Hub Lines.

53
26. Select ne_10m_populated_places_simple as the Hub Point Layer and name as the Hub ID
Attribute. Similarly, selectearthquake_with_places as the Spoke Point Layer and matrix_Tar as
the Spoke Hub ID Attribute. The hub lines algorithm will go through each of earthquake points
and create a line that will join it to the populated place which matches the attribute we
specified.
Click Browse and name the Output Shapefile as earthquake_hub_lines.shp. Click OK to start
the processing.

27. The processing may take a few minutes. You can see the progress on the bottom-left corner
of the QGIS window.
28. Once the processing is done, you will see the earthquake_hub_lines layer loaded in QGIS.
You can see that each earthquake point now has a line that connects it to the nearest
populated place.

54
55
 PRACTICAL NO 12
AIM :Sampling Raster Data using Points or Polygons
Many scientific and environmental datasets come as gridded rasters. Elevation data (DEM) is
also distributed as raster files. In these raster files, the parameter that is being represented is
encoded as the pixel values of the raster. Often, one needs to extract the pixel values at certain
locations or aggregate them over some area. This functionality is available in QGIS via two
plugins - Point SamplingTool and Zonal Statistics plugin.

Procedure

1. Go to Layer ‣Add Raster Layer and browse to the

downloaded us.tmax_nohads_ll_{YYYYMMDD}_float.tif file and click Open.
2. Once the layer is loaded, select the Identify tool and click anywhere on the layer. You will see
the temperature value in celsius as the value or Band 1 at that location.

3. Now unzip the downloaded 2013_Gaz_ua_national.zip file and extract

the 2013_Gaz_ua_national.txt file on your disk. Go to Layer ‣Add Delimited Text Layer.

4. In the Create a Layer from Delimited Text File dialog, click Browse and
open 2013_Gaz_ua_national.txt. Choose Tab under Custom delimiters. The point coordinates
are in Latitude and Longitude, so select INTPTLONG as X field and INTPTLAT as Y field.
Check the Use spatial index box and click OK.

56
5. Now we are ready to extract the temperature values from the raster layer. Install
the Point Sampling Tool plugin. See Using Plugins for details on how to install plugins.

6. Open the plugin dialog from Plugins ‣Analyses ‣Point sampling tool.

7. In the Point Sampling Tool dialog, select 2013_Gaz_ua_national as the Layer containing
sampling points. We must explicitely pick the fields from the input layer that we want in the
output layer. Choose GEOID and NAME fields from the2013_Gaz_ua_national layer. We can
sample values from multiple raster band at once, but since our raster has only 1 band, choose
the us.tmax_nohads_ll_{YYYYMMDD}_float: Band 1. Name the output vector layer
as max_temparature_at_urban_locations.shp. Click the OK to start the sampling process.
Click Close once the process finishes.

57
8. You will see a new layer max_temparature_at_urban_locations loaded in QGIS. Use
the Identify tool to click on any point to see the attributes. You will see the us.tmax_no field -
which contains the raster pixel value at the location of the point.

9. First part of our analysis is over. Let’s remove the unnecessary layers. Hold the Shift key and
select max_temparature_at_urban_locations and 2013_Gaz_ua_national layers. Right-click and
select Remove to remove them from QGIS TOC.
10. Go to Layer ‣Add Vector Layer. Browse to the downloaded tl_2013_us_county.zip file and
click Open. Select thetl_2013_us_county.shp as the layer and click OK.
11. The tl_2013_us_county will be added to QGIS. This layer is in EPSG:4269 NAD83 projection.
This doesn’t match the projection of the raster layer. We will re-project this layer
to EPSG:4326 WGS84 projection.

58
12. Right-click the tl_2013_us_county layer and select Save As...

13. In the Save Vector layer as.. dialog, click Browse and name the output file as counties.shp.
Choose Selected CRS from the CRS dropdown menu. Click Browse and select WGS 84 as the
CRS. Check the Add saved file to map and click OK.

15. Enable the Zonal Statistics Plugins. This is a core plugin so it is already installed. See Using
Plugins to know to how enable core plugins.
16. Go to Raster ‣Zonal statistics ‣Zonal statistics.
17. Select us.tmax_nohads_ll_{YYYYMMDD}_float as the Raster layer and counties as

59
the Polygon layer containing the zones. Enter ZS_ as the Output column prefix. Click OK.

18. The analysis may take some time depending on the size of the dataset.

19. Once the processing finishes, select the counties layer. Use the Identify tool and click on any
county polygon. You will see three new attributes added to the
layer: ZS_count, ZS_mean and ZS_sum. These attributes contain the count of raster pixels,
mean of raster pixel values and sum of raster pixel values respectively. Since we are interested
in average temperature, the ZS_meanfield will be the one to use.

20. Let’s style this layer to create a temperature map. Right-click the counties layer and
select Properties.
21. Switch to the Style tab. Choose Graduated style and select ZS_mean as the Column. Choose
a Color Ramp and Mode of your chose. Click Classify to create the classes. Click OK.

60
Output would be:

61
 PRACTICAL NO.13

AIM :Interpolating Point Data

Procedure
1. Open QGIS. Go to Layer ‣Add Layer ‣Add Vector Layer..
2. Browse to the downloaded Shapefiles.zip file and select it. Click Open.
3. In the Select layers to add... dialog, hold the Shift key and
select Arlington_Soundings_2007_stpl83.shp andBoundary2004_550_stpl83.shp layers.
Click OK.

4. You will see the 2 layers loaded in QGIS. The Boundary2004_550_stpl83 layer represents the
boundary of the lake. Un-check the box next to it in the Table of Contents.

5. This will reveal the data from the second layer Arlington_Soundings_2007_stpl83. Though
the data looks like lines, it is a series of points that are very close.

62
6. Click the Zoom icon and select a small area on the screen. As you zoom closer, you will see
the points. Each point represents a reading taken by a Depth Sounder at the location recorded
by a DGPS equipment.

7. Select the Identify tool and click on a point. You will see the Identify Results panel show up
on the left with the attribute value of the point. In this case, the ELEVATION attribute contains
the depth of the lake at the location. As our task is to create a depth profile and elevation
contours, we will use this values as input for the interpolation.

63
8. Make sure you have the Interpolation plugin enabled. See Using Plugins for how to enable
plugins. Once enabled, go toRaster‣Interpolation ‣Interpolation.
9. In the Interpolation dialog, select Arlington_Soundings_2007_stpl83 as the Vector layers in
the Input panel. Select ELEVATION as the Interpolation attribute. Click Add. Change
the Cellsize X and Cellsize Y values to 5. This value is the size of each pixel in the output grid.
Since our source data is in a projected CRS with Feet-US as units, based on our selection, the
grid size will be 5 feet. Click on the ... button next to Output file and name the output file
as elevation_tin.tif. CLick OK.

10. You will see the new later elevation_tin loaded in QGIS. Right-click the layer and select
Zoom to layer.

64
11. Now you will see the full extent of the created surface. Interpolation does not give accurate
results outside the collection area. Let’s clip the resulting surface with the lake boundary. Go
to Raster ‣Extraction ‣Clipper.

12. Name the Output file as elevation_tin_clipped.tif. Select the Cliiped mode as Mask layer.
Select Boundary2004_550_stpl83 as the Mask layer`. Click OK.

65
13. A new raster elevation_clipped will be loaded in QGIS. We will now style this layer to
show the difference in elevations. Note the min and max elevation values from
the elevation_tin layer. Right-click the elevation_tin_clipped layer and select Properties.

14. Go to the Style tab. Select Render type as Singleband pseudocolor. In the Generate new
color map panel, select Spectralcolor ramp. As we want to create a depth-map as opposed to a
height-map, check the Invert box. This will assign blues to deep areas and reds to shallow areas.
Click Classify.

66
15. Switch to the Tranparency tab. We want to remove the black-pixels from our output. Enter
0 as the Additional no data value. Click OK.

67
16. Now you have a elevation relief map for the lake generated from the individual depth
readings.
Let’s generate contours now. Go to Raster ‣Extraction ‣Contours.
17. In the Contour dialog, enter contours as the Output file for contour lines. We will generate
contour lines at 5ft intervals, so enter 5.00 as the Interval between contour lines. Check
the Attribute name box. Click OK.

18. The contour lines will be loaded as contours layer once the processing is finished. Right-click
the layer and select Properties.

68
19. Go to the Labels tab. Check the Label this layer with box and select ELEV as the field.
Select Curved as the Placement type and click OK.

69
Final Output:

70