Question 1: GNSS Definition & GIS Database for Water Supply Infrastructure [5 + 25 Marks]
Part (a): Define GNSS [5 Marks]
GNSS (Global Navigation Satellite System) is a general term describing any satellite constellation that provides positioning, navigation, and timing (PNT) services on a global or regional basis. Unlike "GPS," which refers only to the U.S. system, GNSS is an umbrella term covering all such systems worldwide.
Major GNSS Constellations
| System | Country / Region | Satellites | Status |
|---|---|---|---|
| GPS (Global Positioning System) | United States | ~31 | Fully Operational |
| GLONASS | Russia | ~24 | Fully Operational |
| Galileo | European Union | ~30 | Operational |
| BeiDou (BDS) | China | ~35 | Fully Operational |
| NavIC (IRNSS) | India | 7 | Regional |
| QZSS | Japan | 4 | Regional |
How GNSS Works
Three Segments of GNSS
| Segment | Description |
|---|---|
| Space Segment | The satellites orbiting Earth, continuously broadcasting navigation signals. |
| Control Segment | Ground stations that monitor satellite health, update orbits, and synchronize clocks. |
| User Segment | GNSS receivers (survey-grade, smartphones, vehicles) that compute position from satellite signals. |
Advantages of Multi-Constellation GNSS
- More satellites visible = better accuracy, especially in urban canyons and hilly terrain
- Faster fix — time-to-first-fix (TTFF) is reduced
- Redundancy — if one system degrades, others maintain positioning
- Better coverage in high-latitude and obstructed environments
Part (b): Preparing a GIS Database of Water Supply Infrastructure Using GNSS [25 Marks]
Scenario: You are a GIS expert tasked with creating a complete GIS vector database of water supply infrastructure (source, pipelines, household benefited) in a municipality. You have manpower for field work.
Step 1: Planning & Preparation
| Task | Details |
|---|---|
| Define scope | Identify the municipal boundary, list all infrastructure types: water sources (springs, reservoirs, wells), pipelines (transmission, distribution), and benefited households. |
| Data model design | Design a geodatabase schema with feature classes for Point (sources, households), Line (pipelines), and Polygon (service zones). |
| Coordinate system | Choose an appropriate CRS (e.g., MUTM for Nepal-based municipalities) and ensure all data are in the same projection. |
| Collect existing data | Gather any existing maps, engineering drawings, pipe layout plans from the municipality. |
| Prepare equipment | GNSS receivers (handheld/survey-grade), smartphones with GIS apps, measurement tape, data collection forms. |
Step 2: Database Schema Design
| Feature Class | Geometry | Key Attributes |
|---|---|---|
| Water Source | Point | Source ID, Name, Type (spring/well/reservoir), Capacity (liters/sec), Elevation, Condition, Photo |
| Pipeline | Line | Pipe ID, Material (HDPE/GI/PVC), Diameter (mm), Length (m), Installation Year, Pressure Zone |
| Household (Benefited) | Point | HH ID, Owner Name, Ward No., Connection Type (tap/shared), No. of Members, Status (active/inactive) |
| Storage Tank | Point | Tank ID, Capacity (liters), Type (RCC/ferro-cement), Elevation, Condition |
| Service Zone | Polygon | Zone ID, Name, Total HH Served, Source Linked, Ward Coverage |
Step 3: Field Data Collection Using GNSS & Mobile GIS Apps
| Tool / App | Purpose | Platform |
|---|---|---|
| QField | Mobile GIS data collection app linked to QGIS projects | Android |
| ArcGIS Field Maps | Esri's mobile field collection app with GNSS support | Android / iOS |
| KoBoToolbox / ODK Collect | Form-based data collection with GPS, photos, attributes | Android |
| SW Maps | Free GIS app for collecting points, lines, polygons with external GNSS | Android |
| Handheld GNSS Receiver | Survey-grade positioning (e.g., Trimble, Leica) | Hardware |
Step 4: Detailed Field Procedures
Collecting Water Sources (Points): Stand at the source with GNSS receiver, wait for satellite lock (4+ satellites, PDOP < 3), record the point, fill attributes (type, capacity, condition), take geotagged photographs.
Tracing Pipelines (Lines): Start at a known point, enable continuous GNSS tracking, walk along the pipe route, pause at junctions to add vertices with attributes (diameter, material, burial depth).
Recording Benefited Households (Points): Visit each household, record GNSS point at the tap/connection, fill in owner name, ward number, connection type, household members, and status.
Step 5: Data Processing in GIS
| Process | Tool | Description |
|---|---|---|
| Import field data | QGIS / ArcGIS | Load collected shapefiles, GeoJSON, or sync from cloud |
| Topology check | Topology Rules | Ensure pipes are connected, no dangles, sources snap to pipe network |
| Attribute validation | Field Calculator | Check for missing values, normalize names, calculate pipe lengths |
| Projection alignment | Reproject tool | Ensure all layers are in the same CRS (e.g., MUTM Zone 45) |
| Service zone creation | Buffer / Dissolve | Create service area polygons from household clusters |
Step 6: Map Production & Outputs
| Map / Output | Purpose |
|---|---|
| Water Source Location Map | Shows all sources with type classification (spring, well, reservoir) |
| Pipeline Network Map | Displays the entire pipe network with diameter and material symbolization |
| Household Coverage Map | Maps all benefited households colored by ward or connection status |
| Service Zone Map | Shows service area polygons overlaid with population density |
| Infrastructure Summary Report | Tabular report: total sources, pipe length, households served, etc. |
Question 2: Interpolation & IDW (Inverse Distance Weighted) [10 Marks]
1. What is Interpolation?
Interpolation is the process of estimating unknown values at unsampled locations based on known values at surrounding sampled points. In GIS, it is the primary technique for creating continuous surfaces (rasters) from discrete point observations.
Types of Interpolation
| Category | Method | Key Idea |
|---|---|---|
| Deterministic | IDW | Closer points have more influence (weight decreases with distance) |
| Spline | Fits a smooth surface through points like bending a rubber sheet | |
| Natural Neighbor | Uses Voronoi tessellation to assign weights based on area | |
| Geostatistical | Kriging | Uses a semivariogram model to account for spatial dependency |
| Co-Kriging | Uses a correlated secondary variable (e.g., elevation for rainfall) |
2. Inverse Distance Weighted (IDW) Method
IDW is a deterministic interpolation technique based on Tobler's First Law of Geography: "Everything is related to everything else, but near things are more related than distant things."
The value at an unknown point is estimated as a weighted average of nearby known values, where the weights decrease as the distance increases.
IDW Formula
$$\hat{Z}(x_0) = \frac{\sum_{i=1}^{n} \frac{Z(x_i)}{d_i^p}}{\sum_{i=1}^{n} \frac{1}{d_i^p}}$$
Where:
- $\hat{Z}(x_0)$ = estimated value at the unknown point
- $Z(x_i)$ = known value at point $i$
- $d_i$ = distance between point $i$ and $x_0$
- $p$ = power parameter (commonly $p = 2$)
- $n$ = number of nearby sample points used
How IDW Works — Visual Explanation
Numerical Example (p = 2)
| Point | Value Z | Distance d (km) | 1/d² | Weight (w) | w × Z |
|---|---|---|---|---|---|
| A | 20 | 2 | 0.250 | 0.631 | 12.63 |
| B | 30 | 4 | 0.0625 | 0.158 | 4.73 |
| C | 10 | 6 | 0.0278 | 0.070 | 0.70 |
| D | 25 | 3 | 0.1111 | 0.281 | 7.01 |
| Sum | 0.396 | 1.000 | ≈ 21.4 | ||
Estimated Value ≈ 21.4
Point A (closest, d=2) dominates the result, pulling it towards 20.
Effect of the Power Parameter (p)
| Power (p) | Behavior |
|---|---|
| p = 1 | Gentle decay — distant points still have moderate influence. Smoother surface. |
| p = 2 (default) | Balanced decay — standard choice in most GIS software. |
| p > 2 | Steep decay — only the very closest points dominate. More "peaky" surface. |
Advantages & Limitations
Advantages
- Simple, intuitive, easy to understand
- Fast computation
- No assumptions about data distribution
- Works well when sample points are evenly distributed
- Available in all GIS software
Limitations
- Bull's eye effect — concentric patterns around sample points
- Cannot extrapolate beyond data range
- Sensitive to clustered data
- Does not consider spatial autocorrelation
- No estimate of prediction uncertainty
IDW vs. Other Methods
| Feature | IDW | Spline | Kriging |
|---|---|---|---|
| Type | Deterministic | Deterministic | Geostatistical |
| Surface | Passes through points | Smooth through points | Best unbiased estimate |
| Error estimate | No | No | Yes (variance map) |
| Computation | Fast | Moderate | Slow (needs variogram) |
| Best for | Dense, evenly spaced data | Smooth phenomena | Sparse data, need uncertainty |
Question 3: Seasonal Land-Use Variation Analysis (10 Years) [30 Marks]
Scenario: Analyze seasonal variation of land-use before and after monsoon in a municipality for the last 10 years using remote sensing and GIS tools. Software-agnostic answer expected.
1. Choice of Remote Sensing Data
| Satellite / Sensor | Spatial Resolution | Temporal Resolution | Cost | Suitability |
|---|---|---|---|---|
| Landsat 8/9 (OLI) | 30 m / 15 m pan | 16 days | Free (USGS) | Best — free, 10+ year archive, ideal for municipality-level LULC |
| Sentinel-2 (MSI) | 10 m / 20 m | 5 days | Free (ESA) | Excellent — higher resolution, archive from 2015 |
| Landsat 7 (ETM+) | 30 m | 16 days | Free | Useful for years before 2013 (has SLC-off issue from 2003) |
| MODIS | 250–500 m | Daily | Free (NASA) | Too coarse for municipality-level |
Image Selection Criteria
| Criterion | Pre-Monsoon Image | Post-Monsoon Image |
|---|---|---|
| Season / Month | March–May (dry, clear sky) | October–November (after monsoon, vegetation peak) |
| Cloud cover | < 10% over study area | < 10% (may need compositing) |
| Number of images | 1 per year × 10 years = 10 | 1 per year × 10 years = 10 |
Total images needed: 20 scenes (10 pre-monsoon + 10 post-monsoon)
2. Complete Workflow
Pre-processing Details
| Correction Type | What It Does | Why Needed |
|---|---|---|
| Radiometric Calibration | Converts DN → TOA radiance/reflectance | Makes pixel values comparable across dates |
| Atmospheric Correction | Removes haze → Surface Reflectance | Essential for multi-temporal comparison |
| Geometric Correction | Ensures images align spatially | Pixel-to-pixel comparison needs exact match |
| Cloud Masking | Removes cloud pixels using QA bands | Critical for post-monsoon images |
Classification Scheme
| Class | Description | Expected Seasonal Change |
|---|---|---|
| Agriculture | Cultivated land, paddy, crops | Most variable — bare pre-monsoon, green post-monsoon |
| Forest | Dense/sparse forest, tree cover | Greener post-monsoon |
| Built-up | Urban areas, roads, buildings | Relatively stable |
| Water Body | Rivers, ponds, reservoirs | Expands post-monsoon, shrinks pre-monsoon |
| Barren Land | Rocky, exposed soil, fallow | More visible pre-monsoon |
| Wetland / Marsh | Seasonal water-logged areas | Appears mainly post-monsoon |
Change Detection: Post-Classification Comparison
Reason: Directly produces a "from–to" change matrix, works across different sensors/dates, and provides clear interpretable change classes.
Transition Matrix Example
| Agri (Post) | Forest (Post) | Built-up (Post) | Water (Post) | Barren (Post) | |
|---|---|---|---|---|---|
| Agri (Pre) | 850 ha | 20 ha | 150 ha | 5 ha | 25 ha |
| Forest (Pre) | 30 ha | 680 ha | 10 ha | 0 ha | 5 ha |
| Built-up (Pre) | 5 ha | 0 ha | 400 ha | 0 ha | 2 ha |
| Water (Pre) | 10 ha | 2 ha | 0 ha | 180 ha | 8 ha |
| Barren (Pre) | 40 ha | 15 ha | 20 ha | 5 ha | 120 ha |
10-Year Trend (Illustrative)
Summary of Key Decisions
| Decision | Choice | Reason |
|---|---|---|
| Satellite data | Landsat 8/9 + Landsat 7 | Free, 30 m, longest archive |
| Pre-processing | Surface Reflectance (Level-2) | Removes atmospheric effects for comparison |
| Classification | Supervised (MLC / Random Forest) | More accurate, consistent across dates |
| Change detection | Post-classification comparison | Gives from–to transition information |
| Seasonal timing | Mar–May vs Oct–Nov | Maximum seasonal contrast |
Question 4: Parameters Extracted from Elevation Data [10 Marks]
Elevation data (DEM/DTM) is one of the most powerful GIS datasets. By analyzing how elevation values change across a surface, we can derive numerous terrain parameters.
1. Slope
Slope measures the steepness or gradient of a terrain surface — the rate of change of elevation over horizontal distance, expressed in degrees (0°–90°) or percent.
| Slope Range | Terrain Type | Application |
|---|---|---|
| 0°–5° | Flat / Gentle | Suitable for agriculture, construction |
| 5°–15° | Moderate | Roads feasible, terrace farming |
| 15°–30° | Steep | Erosion risk, needs engineering measures |
| > 30° | Very Steep | Landslide prone, unsuitable for development |
2. Aspect
Aspect is the compass direction a slope faces — measured 0°–360° (0°/360° = North, 90° = East, 180° = South, 270° = West).
3. Contour Lines
Contour lines are lines connecting points of equal elevation — generated from a DEM for topographic mapping, engineering design, and volume estimation.
4. Hillshade
Hillshade simulates how terrain looks illuminated from a specific sun angle — creating realistic 3D shaded relief useful for visualization and map backgrounds.
Summary Table
| Parameter | What It Measures | Unit | Key Application |
|---|---|---|---|
| Slope | Steepness / gradient | Degrees or % | Landslide susceptibility, road design, erosion |
| Aspect | Compass direction slope faces | Degrees (0°–360°) | Solar radiation, crop suitability, building orientation |
| Contour Lines | Lines of equal elevation | Meters (interval) | Topographic mapping, volume calculation, engineering |
| Hillshade | Simulated light/shadow | 0–255 (grayscale) | 3D visualization, map backgrounds |
Question 5: Flow Accumulation in River Channel Analysis [10 Marks]
1. What is Flow Accumulation?
Flow accumulation calculates, for each cell in a DEM, the total number of upstream cells whose water would drain through it. It answers: "How much water collects at each point on the landscape?"
2. Hydrological Workflow
3. How Flow Accumulation Works (D8 Example)
4. Stream Extraction via Threshold
Stream = (Flow Accumulation > Threshold)
Common threshold: 500–1000 cells
| Threshold | Effect |
|---|---|
| Low (e.g., 100) | Extracts many small tributaries → denser network |
| High (e.g., 5000) | Only major rivers → sparser network |
5. Applications
| Application | How Flow Accumulation Helps |
|---|---|
| Stream network extraction | Identifies where streams form based on accumulated drainage |
| Watershed delineation | Combined with flow direction, defines catchment boundaries |
| Stream ordering | Strahler/Shreve ordering on extracted stream links |
| Flood risk assessment | High accumulation = most vulnerable to flooding |
| Erosion modeling | Flow concentration indicates high erosion potential |
| Infrastructure planning | Bridges, culverts, drainage design |
Question 6: Low Pass & High Pass Filters in Remote Sensing [5 + 5 Marks]
Spatial filters are convolution-based image processing techniques that replace each pixel's value with a weighted combination of its neighbors using a moving kernel window.
Part (a): Low Pass Filter [5 Marks]
A low pass filter (smoothing/averaging filter) allows low spatial frequencies to pass through while suppressing high frequencies. It replaces each pixel with the average of its neighbors.
- Smooths the image, reduces sharp variations
- Reduces noise (salt-and-pepper)
- Blurs edges — boundaries become less distinct
- Generalizes — reveals broad patterns
Part (b): High Pass Filter [5 Marks]
A high pass filter (edge enhancement/sharpening filter) allows high spatial frequencies to pass through while suppressing low frequencies. It highlights edges, boundaries, and fine details.
- Sharpens the image — enhances edges and boundaries
- Highlights rapid changes — roads, field edges, lineaments
- Removes gradual trends — suppresses smooth areas
- May amplify noise (side effect)
Comparison: Low Pass vs. High Pass
| Feature | Low Pass Filter | High Pass Filter |
|---|---|---|
| Also called | Smoothing / Averaging | Sharpening / Edge enhancement |
| Passes through | Low spatial frequencies | High spatial frequencies |
| Suppresses | High frequencies (sharp changes) | Low frequencies (gradual changes) |
| Effect on image | Blurs / smooths | Sharpens / enhances edges |
| Kernel weights | All positive, equal | Mixed (+/-), large positive center |
| Noise | Reduces noise | May amplify noise |
| Best for | Noise removal, generalization | Edge detection, feature extraction |