Scribd
Upload
135 views

Generating Potential Irrigable Area Using QGIS Software

This document provides instructions for using QGIS software to generate potential irrigable areas. It involves adding data files as layers, selecting an area boundary, reclassifying a slope raster, extracting land use/cover polygons, and selecting cropland polygons over 0.5 hectares on slopes of 0-3% to identify potential irrigable areas. The process is repeated for other slope classifications and land cover types.

Uploaded by

Rhea Celzo
Copyright
© © All Rights Reserved
Available Formats
Download as PPTX, PDF, TXT or read online on Scribd
135 views

Generating Potential Irrigable Area Using QGIS Software

This document provides instructions for using QGIS software to generate potential irrigable areas. It involves adding data files as layers, selecting an area boundary, reclassifying a slope raster, extracting land use/cover polygons, and selecting cropland polygons over 0.5 hectares on slopes of 0-3% to identify potential irrigable areas. The process is repeated for other slope classifications and land cover types.

Uploaded by

Rhea Celzo
Copyright
© © All Rights Reserved
Available Formats
Download as PPTX, PDF, TXT or read online on Scribd
You are on page 1/ 75

Generating Potential

Irrigable Area Using


QGIS Software
SOFTWAR
E
∙ QGIS 3.14.15-Pi version (can also use higher version)

FILES

∙ IFSAR/LIDAR File
∙ Slope File (can be extracted also from IFSAR/LIDAR)
∙ Sentinel-2 ESRI 2020 Land Use/Land Cover File
∙ https://www.arcgis.com/apps/instant/media/index.html?
appid=fc92d38533d440078f17678ebc20e8e2
∙ Municipal Boundary Shapefile/Soil Map Shapefile
PART Adding files to Layers panel

1
PART Area Boundary Selection

2
PART Slope Reclassification

3
PART LULC Extraction

4
PART
1
1. Open QGIS. Click ‘New Empty Project’ or click ‘Project’ menu and click ‘New’.
*Set Current CRS: WGS 84 / UTM zone 51N
2. On ‘Web’ menu, click ‘QuickMapServices’ (If there’s none, Install it first
in the ‘Plugins’ menu) and click ‘Search QMS’.
3. On the ‘Search QMS’ panel, type in the search box ‘Google Satellite
Hybrid’ and click ‘Add’ (You can close the panel after you add the layer).
4. On ‘Layer’ menu, click ‘Add Layer’ and click ‘Add Raster Layer’.
5. On ‘Data Source Manager | Raster’ dialog box, click on the three dots button
beside the ‘Raster Dataset(s)’ box to browse the Raster Files
(IFSAR/LIDAR/LULC/Slope) we are going to use in the project.
6. On ‘Open GDAL Supported Raster Dataset(s)’ dialog box, browse to the folder
where you saved the Raster File(s). Click the file(s) and click ‘Open’.
7. Going back on ‘Data Source Manager | Raster’ dialog box, the ‘Raster
Dataset(s)’ box has already filled with the Raster File(s) we are going to use in the
project. Click ‘Add’ (Once another dialog box pops up, just click ‘OK’) and click
‘Close’.
8. On ‘Layer’ menu, click ‘Add Layer’ and click ‘Add Vector Layer’.
9. On ‘Data Source Manager | Vector’ dialog box, click on the three dots button
beside the ‘Vector Dataset(s)’ box to browse the Vector Files (Municipal Boundary
Shapefile/Soil Map Shapefile) we are going to use in the project.
10. On ‘Open OGR Supported Vector Dataset(s)’ dialog box, browse to the folder
where you saved the Vector File(s). Click the file(s) and click ‘Open’.
11. Going back on ‘Data Source Manager | Vector’ dialog box, the ‘Vector Dataset(s)’
box has already filled with the Vector File(s) we are going to use in the project. Click
‘Add’ and click ‘Close’.
PART
2
12. On ‘Layers’ panel, click ‘Municipal Boundary’ layer (make sure the ‘Municipal
Boundary’ layer is above the other layer for later selection) then click ‘Select
Features’ in the selection toolbar and select a municipality on the map view/canvass.
13. On ‘Vector’ menu, click ‘Geometry Tools’ and click ‘Collect Geometries’.
14. On ‘Collect Geometries’ dialog box, click the drop down arrow below ‘Input
layer’ and select ‘Municipal Boundary’ vector layer as its input layer. Check the box
for ‘Selected features only’ to extract only the municipality selected earlier on the
map. Click ‘Run’ and click ‘Close’.
15. On ‘Layers’ panel, right click ‘Collected’ layer and click ‘Rename Layer’.
Rename to ‘Amlan Collected’.
16. On ‘Raster’ menu, click ‘Extraction’ and click ‘Clip Raster by Mask Layer’.
17. On ‘Clip Raster by Mask Layer’ dialog box, click the drop down arrow below ‘Input
layer’ and select the ‘Slope’ raster layer as its input layer. Then, click the drop down
arrow below ‘Mask layer’ and select the ‘Amlan Collected’ vector layer as its mask layer
(set ‘Target CRS [optional]’ to ‘EPSG:32651 – WGS 84 / UTM zone 51N’). Click ‘Run’ and
click ‘Close’.
18. Rename the ‘Clipped (mask)’ layer to ‘Amlan Slope’.
(Apply this step if the QGIS version has no Target CRS option when clipping raster file)

19. On ‘Raster’ menu, click ‘Projections’ and click ‘Warp (Reproject)’.

20. On ‘Warp (Reproject)’ dialog box, click the drop down arrow below ‘Input
layer’ and select the ‘Amlan Slope’ raster layer as its input layer. Then, click the
drop down arrow below ‘Target CRS’ and select the ‘Project CRS: EPSG:32651 –
WGS 84 / UTM zone 51N’ as its target CRS. Click ‘Run’ and click ‘Close’.

21. Rename the ‘Reprojected’ layer to ‘Amlan Slope Reprojected’


22. On ‘Processing’ menu, click ‘Toolbox’ (If there’s none, Install it first in the
‘Plugins’ menu).
23. On ‘Processing Toolbox’ panel, click the drop down arrow of ‘Raster
Analysis’ and double click ‘Reclassify by table’.
24. On ‘Reclassify by Table’ dialog box, click the drop down arrow below ‘Raster
layer’ and select ‘Amlan Slope (Reprojected)’ raster layer. Then, click the three
dots button below ‘Reclassification table’.
25. In the Reclassification table, click ‘Add Row’ three times and input the value
as shown in the picture below and click ‘OK’.
26. Going back on ‘Reclassify by Table’ dialog box, in ‘Advanced parameters’, click
the drop down arrow below ‘Range boundaries’ and select ‘min <= value < max’.
Click ‘Run’ and click ‘Close’.
27. On ‘Raster’ menu, click ‘Conversion’ and click ‘Polygonize (Raster to Vector)’.
28. On ‘Polygonize (Raster to Vector)’ dialog box, click the drop down arrow below
‘Input layer’ and select the ‘Reclassified raster’ raster layer as its input layer. Click
‘Run’ and click ‘Close’.
29. Rename the ‘Vectorized’ layer to ‘Amlan Slope Vectorized’.
30. On ‘Processing Toolbox’ panel, click the drop down arrow of ‘Vector
geometry’ and double click ‘Fix geometries’.
31. On ‘Fix Geometries’ dialog box, click the drop down arrow below ‘Input layer’
and select ‘Amlan Slope Vectorized’ vector layer as its input layer. Click ‘Run’ and
click ‘Close’.
32. Rename the ‘Fixed geometries’ layer to ‘Amlan
Slope Fixed geometries’.
PART
3
33. On ‘Layers’ panel, click ‘Amlan Slope Fixed geometries’ layer then click ‘Select
Features by Value’ in the selection toolbar.
34. On ‘Select Features’ dialog box, type ‘1’ (DN value assigned for 0 to 3% slope)
for the DN value. Click ‘Select features’ and click ‘Close’.
35. On ‘Vector’ menu, click ‘Geometry Tools’ and click ‘Collect Geometries’.
36. On ‘Collect Geometries’ dialog box, click the drop down arrow below ‘Input layer’
and select ‘Amlan Slope Fixed geometries’ vector layer as its input layer. Check the
box for ‘Selected features only’. Click ‘Run’ and click ‘Close’.
37. Rename the new ‘Collected’ layer to ‘Amlan 0 to 3 Slope Collected’.
38. On ‘Raster’ menu, click ‘Extraction’ and click ‘Clip Raster by Mask Layer’.
39. On ‘Clip Raster by Mask Layer’ dialog box, click the drop down arrow below ‘Input
layer’ and select the ‘LULC’ raster layer as its input layer. Then, click the drop down
arrow below ‘Mask layer’ and select the ‘Amlan 0 to 3 Slope Collected’ vector layer as
its mask layer (set ‘Target CRS [optional]’ to ‘EPSG:32651 – WGS 84 / UTM zone 51N’).
Click ‘Run’ and click ‘Close’.
40. Rename the new ‘Clipped (mask)’ layer to ‘Amlan 0 to 3 Slope LULC’.
41. On ‘Raster’ menu, click ‘Conversion’ and click ‘Polygonize (Raster to Vector)’.
42. On ‘Polygonize (Raster to Vector)’ dialog box, click the drop down arrow below
‘Input layer’ and select the ‘Amlan 0 to 3 Slope LULC’ raster layer as its input layer.
Click ‘Run’ and click ‘Close’.
43. Rename the new ‘Vectorized’ layer to ‘Amlan 0 to 3 Slope LULC Vectorized’.
44. On ‘Layers’ panel, click ‘Amlan 0 to 3 Slope LULC Vectorized’ then click ‘Open
Field Calculator’ in the attributes toolbar.
45. On ‘Field Calculator’ dialog box, type ‘Area’ for the ‘Output field name’ and
type ‘$area’ in the expression box below ‘Expression’ tab. Click ‘OK’.
46. Click ‘Save Layer Edits’ in the digitizing toolbar and click ‘Toggle Editing’.
PART
4
47. On ‘Layers’ panel, click ‘Amlan 0 to 3 Slope LULC Vectorized’ then click ‘Select
Features by Value’ in the selection toolbar.
48. On ‘Select Features by Value’ dialog box, type ‘5’ (DN value equivalent to
cropland) for the ‘DN’ value and type ‘5000’ (sq. m), equivalent to a half hectare, for
the ‘Area’ then click the drop down arrow beside the area box and select ‘Greater than
(>)’. Click ‘Select features’ and click ‘Close’.
49. On ‘Vector’ menu, click ‘Geometry Tools’ and click ‘Collect Geometries’.
50. On ‘Collect Geometries’ dialog box, click the drop down arrow below ‘Input
layer’ and select ‘Amlan 0 to 3 Slope LULC’ vector layer as its input layer. Check the
box for ‘Selected features only’ then click the three dots below ‘Collected’ and
select ‘Save to File’.
51. On ‘Save file’ dialog box, type ‘Amlan Cropland 0 to 3 Slope’ for the file name
and click ‘Save’.
52. Going back to ‘Collect Geometries’ dialog box, Click ‘Run’ and click ‘Close’.
∙ Repeat Part 4 with different DN value to get Rangeland
∙ Repeat Part 3 with different DN values to get another Slope Classification and
proceed again with Part 4
∙ Repeat Part 2 with different Area Boundary when already done with Part 3
and 4

***END***

You might also like