A.
Slope
1. In QGIS application, load the municipal boundary shapefile of Naval
and add OSM base map.
[Link] the canvas so that the whole AOI is displayed.
3. Use SRTM-Downloader Plugin, click the ‘set the canvas extent’ and
‘download’. You will also need to login your account.
4. SRTM image will automatically be loaded in the canvas.
5. Reproject the raster into EPSG:32651-WGS 84/UTM zone 51N (used to
map locations in the Philippines with accurate, meter-based coordinates)
and clip it into the AOI.
6. Go to Raster > Analysis > Slope.
7. For slope classification:
0-5° - Flat Terrain
5-15° - Gentle Slope
15-30° - Moderate Slope
30-45° - Steep Slope
> 45° - Ver Steep
Classify the raster using ‘Reclassify by Table’ in Raster Analysis.
Result:
8. Area Computation
a. Convert the raster to vector
b. Dissolve the vector to combine the polygons with the same
class.
c. Open the attribute table and click field calculator.
d. Save toggle editing
B. Soil
1. Data Collection
a. Add a new script and run. The available data is from SoilGrids,
specifically the silt, clay, and sand layers. These data will be
processed in QGIS for soil classification
b. In tasks tab, download the data. The data will be saved directly
to a Gdrive.
2. Data Processing
a. Load the three Raster image into the QGIS.
b. Reproject
c. Convert the raster’s data into percentage using raster
calculator.
d. Fill the no-data gaps in the raster images, and repeat the
process for each image.
e. Clip each filled images into the AOI.
f. For classification, use the raster calculator.
C. Rainfall
1. The data is obtained from CHIRPS via Google Earth Engine. Use the
same process as with the soil data: run the script, and the data will
be saved to Google Drive.
Script:
var naval = naval
var chirps =
chirpsDaily
.filterDate('2014-01-01', '2024-12-31')
.filterBounds(naval); var years =
[Link](2014, 2024); var
annualRainfall = [Link](
[Link](function(year){
var start = [Link](year, 1, 1);
var end = [Link](1, 'year');
var yearlySum = [Link](start, end).sum();
return [Link]('year', year);
})
);
var meanAnnualRainfall =
[Link]().clip(naval); var visParams = {
min: 0,
max: 4000,
palette: ['blue', 'limegreen', 'yellow', 'orange',
'red']
};
[Link](naval, 11);
[Link](meanAnnualRainfall, visParams, 'Mean Annual Rainfall 2014-
2024');
[Link]({
collection: points,
description: 'Naval_Mean_Annual_Rainfall_Points',
folder: 'GEE_exports',
fileNamePrefix:
'naval_mean_annual_rainfall_points_2014_2024',
fileFormat: 'SHP'
});
2. The data must be a point data for interpolation.
3. In the Processing Toolbox, search for IDW Interpolation. Run the
tool, and the resulting output will be a raster image.
4. Clip the interpolated layer into the AOI.
5. Reproject
6. Reclassify by table. (Equal interval)
7. Result