← Back to Index
Internal Assessment • Chaitra 9, 2082 • Total: 75 Minutes

GIS Applications & Remote Sensing — All Questions Solved

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.

Key Point: GNSS = the entire family of satellite navigation systems. GPS is just one member of that family.

Major GNSS Constellations

System Country / Region Satellites Status
GPS (Global Positioning System)United States~31Fully Operational
GLONASSRussia~24Fully Operational
GalileoEuropean Union~30Operational
BeiDou (BDS)China~35Fully Operational
NavIC (IRNSS)India7Regional
QZSSJapan4Regional

How GNSS Works

GPS GLONASS Galileo BeiDou Receiver GNSS receiver receives signals from multiple constellations simultaneously
Figure 1.1: GNSS concept — a receiver processes signals from multiple satellite constellations to determine precise location.

Three Segments of GNSS

SegmentDescription
Space SegmentThe satellites orbiting Earth, continuously broadcasting navigation signals.
Control SegmentGround stations that monitor satellite health, update orbits, and synchronize clocks.
User SegmentGNSS receivers (survey-grade, smartphones, vehicles) that compute position from satellite signals.

Advantages of Multi-Constellation GNSS

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

TaskDetails
Define scopeIdentify the municipal boundary, list all infrastructure types: water sources (springs, reservoirs, wells), pipelines (transmission, distribution), and benefited households.
Data model designDesign a geodatabase schema with feature classes for Point (sources, households), Line (pipelines), and Polygon (service zones).
Coordinate systemChoose an appropriate CRS (e.g., MUTM for Nepal-based municipalities) and ensure all data are in the same projection.
Collect existing dataGather any existing maps, engineering drawings, pipe layout plans from the municipality.
Prepare equipmentGNSS receivers (handheld/survey-grade), smartphones with GIS apps, measurement tape, data collection forms.

Step 2: Database Schema Design

Feature ClassGeometryKey Attributes
Water SourcePointSource ID, Name, Type (spring/well/reservoir), Capacity (liters/sec), Elevation, Condition, Photo
PipelineLinePipe ID, Material (HDPE/GI/PVC), Diameter (mm), Length (m), Installation Year, Pressure Zone
Household (Benefited)PointHH ID, Owner Name, Ward No., Connection Type (tap/shared), No. of Members, Status (active/inactive)
Storage TankPointTank ID, Capacity (liters), Type (RCC/ferro-cement), Elevation, Condition
Service ZonePolygonZone ID, Name, Total HH Served, Source Linked, Ward Coverage
Water Supply Infrastructure — GIS Vector Data Model Source Point feature Tank Point feature Transmission Line Service Zone (Polygon) Distribution Household Points Legend Water Source (Point) Storage Tank (Point) Pipeline (Line) Service Zone (Polygon) Household (Point)
Figure 1.2: Schematic of GIS vector data model for water supply infrastructure.

Step 3: Field Data Collection Using GNSS & Mobile GIS Apps

Tool / AppPurposePlatform
QFieldMobile GIS data collection app linked to QGIS projectsAndroid
ArcGIS Field MapsEsri's mobile field collection app with GNSS supportAndroid / iOS
KoBoToolbox / ODK CollectForm-based data collection with GPS, photos, attributesAndroid
SW MapsFree GIS app for collecting points, lines, polygons with external GNSSAndroid
Handheld GNSS ReceiverSurvey-grade positioning (e.g., Trimble, Leica)Hardware
1. Prepare Load forms + base map 2. Navigate Go to each feature 3. Collect GNSS + attributes + photo 4. Trace Lines Walk along pipeline 5. QC in Field Check completeness 6. Sync / Export Upload to GIS 7. Clean & Edit Fix topology in GIS 8. Map Output Thematic maps + reports Complete Field-to-Map Workflow Using GNSS-enabled mobile GIS applications for water supply mapping
Figure 1.3: Step-by-step workflow from field data collection using GNSS/mobile GIS apps to final map output.

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

ProcessToolDescription
Import field dataQGIS / ArcGISLoad collected shapefiles, GeoJSON, or sync from cloud
Topology checkTopology RulesEnsure pipes are connected, no dangles, sources snap to pipe network
Attribute validationField CalculatorCheck for missing values, normalize names, calculate pipe lengths
Projection alignmentReproject toolEnsure all layers are in the same CRS (e.g., MUTM Zone 45)
Service zone creationBuffer / DissolveCreate service area polygons from household clusters

Step 6: Map Production & Outputs

Map / OutputPurpose
Water Source Location MapShows all sources with type classification (spring, well, reservoir)
Pipeline Network MapDisplays the entire pipe network with diameter and material symbolization
Household Coverage MapMaps all benefited households colored by ward or connection status
Service Zone MapShows service area polygons overlaid with population density
Infrastructure Summary ReportTabular report: total sources, pipe length, households served, etc.
Exam Tip: Always mention: (1) the data model/schema, (2) field tools and apps used, (3) GNSS procedures, (4) data processing/cleaning in GIS, and (5) final map outputs. Including a simple diagram of the workflow earns extra marks.

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.

Why do we need interpolation? We can never measure a phenomenon (e.g., elevation, rainfall, temperature) at every single location. Interpolation fills the gaps using mathematical models.

Types of Interpolation

CategoryMethodKey Idea
DeterministicIDWCloser points have more influence (weight decreases with distance)
SplineFits a smooth surface through points like bending a rubber sheet
Natural NeighborUses Voronoi tessellation to assign weights based on area
GeostatisticalKrigingUses a semivariogram model to account for spatial dependency
Co-KrigingUses a correlated secondary variable (e.g., elevation for rainfall)
Interpolation Concept: From Points to Continuous Surface Sample Points 15 22 18 25 12 20 ? What value at the unknown location? IDW Continuous Surface Every cell gets an estimated value
Figure 2.1: Known sample points (left) are used to estimate a continuous raster surface (right) using IDW.

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

IDW: Weight Calculation for Unknown Point ? A Z=20 d₁ = 2 km B Z=30 d₂ = 4 km C Z=10 d₃ = 6 km D Z=25 d₄ = 3 km Closer points (A, D) get HIGHER weights → more influence on result Farther points (B, C) get LOWER weights → less influence on result
Figure 2.2: IDW assigns higher weights to nearby points (A, D) and lower weights to distant points (B, C).

Numerical Example (p = 2)

PointValue ZDistance d (km)1/d²Weight (w)w × Z
A2020.2500.63112.63
B3040.06250.1584.73
C1060.02780.0700.70
D2530.11110.2817.01
Sum0.3961.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)

Effect of Power Parameter (p) on Weight Decay Distance (d) Weight (w) p = 1 p = 2 p = 4 High p → only nearest points matter Low p → distant points still influential
Figure 2.3: As power p increases, weight drops off more steeply with distance. p=2 is the most common choice.
Power (p)Behavior
p = 1Gentle decay — distant points still have moderate influence. Smoother surface.
p = 2 (default)Balanced decay — standard choice in most GIS software.
p > 2Steep 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

FeatureIDWSplineKriging
TypeDeterministicDeterministicGeostatistical
SurfacePasses through pointsSmooth through pointsBest unbiased estimate
Error estimateNoNoYes (variance map)
ComputationFastModerateSlow (needs variogram)
Best forDense, evenly spaced dataSmooth phenomenaSparse data, need uncertainty
Exam Tip: For 10 marks: (1) Define interpolation, (2) State Tobler's Law, (3) Write formula, (4) Explain power parameter, (5) List advantages/limitations. A diagram earns extra marks.

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 / SensorSpatial ResolutionTemporal ResolutionCostSuitability
Landsat 8/9 (OLI)30 m / 15 m pan16 daysFree (USGS)Best — free, 10+ year archive, ideal for municipality-level LULC
Sentinel-2 (MSI)10 m / 20 m5 daysFree (ESA)Excellent — higher resolution, archive from 2015
Landsat 7 (ETM+)30 m16 daysFreeUseful for years before 2013 (has SLC-off issue from 2003)
MODIS250–500 mDailyFree (NASA)Too coarse for municipality-level
Recommended: Use Landsat 8/9 (2013–present) + Landsat 7 (earlier years) for full 10-year coverage. If entirely after 2015, Sentinel-2 at 10 m resolution is preferable.

Image Selection Criteria

CriterionPre-Monsoon ImagePost-Monsoon Image
Season / MonthMarch–May (dry, clear sky)October–November (after monsoon, vegetation peak)
Cloud cover< 10% over study area< 10% (may need compositing)
Number of images1 per year × 10 years = 101 per year × 10 years = 10

Total images needed: 20 scenes (10 pre-monsoon + 10 post-monsoon)

2. Complete Workflow

Complete Workflow: Seasonal Land-Use Change Analysis 1. Data Acquisition & Download 2. Pre-processing & Correction 3. Study Area Extraction (Clip) 4. Land-Use Classification 5. Accuracy Assessment 6. Change Detection Analysis 7. Statistical & Trend Analysis 8. Map Production & Reporting USGS Earth Explorer, Copernicus Hub Atmospheric, radiometric, geometric Municipality boundary shapefile Supervised / unsupervised Confusion matrix, Kappa Post-classification comparison Area statistics, transition matrix Thematic maps, charts, reports Applicable in QGIS, ArcGIS, ERDAS, ENVI, Google Earth Engine
Figure 3.1: Complete 8-step workflow for seasonal land-use change analysis over 10 years.

Pre-processing Details

Correction TypeWhat It DoesWhy Needed
Radiometric CalibrationConverts DN → TOA radiance/reflectanceMakes pixel values comparable across dates
Atmospheric CorrectionRemoves haze → Surface ReflectanceEssential for multi-temporal comparison
Geometric CorrectionEnsures images align spatiallyPixel-to-pixel comparison needs exact match
Cloud MaskingRemoves cloud pixels using QA bandsCritical for post-monsoon images

Classification Scheme

ClassDescriptionExpected Seasonal Change
AgricultureCultivated land, paddy, cropsMost variable — bare pre-monsoon, green post-monsoon
ForestDense/sparse forest, tree coverGreener post-monsoon
Built-upUrban areas, roads, buildingsRelatively stable
Water BodyRivers, ponds, reservoirsExpands post-monsoon, shrinks pre-monsoon
Barren LandRocky, exposed soil, fallowMore visible pre-monsoon
Wetland / MarshSeasonal water-logged areasAppears mainly post-monsoon
Supervised Classification (Recommended): User selects training samples; algorithm (Maximum Likelihood / Random Forest) assigns all pixels. More accurate and consistent across dates than unsupervised methods.

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 ha20 ha150 ha5 ha25 ha
Forest (Pre)30 ha680 ha10 ha0 ha5 ha
Built-up (Pre)5 ha0 ha400 ha0 ha2 ha
Water (Pre)10 ha2 ha0 ha180 ha8 ha
Barren (Pre)40 ha15 ha20 ha5 ha120 ha
Diagonal = no change; off-diagonal = change. Red = significant conversion.

10-Year Trend (Illustrative)

10-Year Land-Use Trend Year Area (ha) 2016201720182019202020212022202320242025 Agri ↓ Built ↑ Forest Water
Figure 3.2: Illustrative 10-year trend — built-up increasing, agricultural land decreasing.

Summary of Key Decisions

DecisionChoiceReason
Satellite dataLandsat 8/9 + Landsat 7Free, 30 m, longest archive
Pre-processingSurface Reflectance (Level-2)Removes atmospheric effects for comparison
ClassificationSupervised (MLC / Random Forest)More accurate, consistent across dates
Change detectionPost-classification comparisonGives from–to transition information
Seasonal timingMar–May vs Oct–NovMaximum seasonal contrast
Exam Tip (30 marks): Structure: (1) Data choice + reason [6], (2) Pre-processing [5], (3) Classification + reason [6], (4) Accuracy assessment [3], (5) Change detection + reason [5], (6) Statistics [3], (7) Map outputs [2]. Include workflow diagram.

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 — Rate of Elevation Change Run (horizontal) Rise Slope angle Slope = arctan(Rise / Run)
Figure 4.1: Slope is the angle of the terrain surface, calculated from rise over run.
Slope RangeTerrain TypeApplication
0°–5°Flat / GentleSuitable for agriculture, construction
5°–15°ModerateRoads feasible, terrace farming
15°–30°SteepErosion risk, needs engineering measures
> 30°Very SteepLandslide 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).

Aspect — Direction of Slope Face N (0°) E (90°) S (180°) W (270°) NE: Morning sun SE: Max sun SW: Afternoon NW: Shade South-facing slopes receive maximum sunlight (N. Hemisphere)
Figure 4.2: Aspect shows slope orientation. South-facing = most sunlight; north-facing = cooler/wetter.

3. Contour Lines

Contour lines are lines connecting points of equal elevation — generated from a DEM for topographic mapping, engineering design, and volume estimation.

Contour Lines from DEM 100 m 200 m 300 m 400 m 500 m Wide spacing = gentle slope Close spacing = steep slope Contour interval = 100 m
Figure 4.3: Closely spaced contours = steep terrain; widely spaced = gentle slopes.

4. Hillshade

Hillshade simulates how terrain looks illuminated from a specific sun angle — creating realistic 3D shaded relief useful for visualization and map backgrounds.

Hillshade — Simulated Terrain Illumination Sun Illuminated Shadow Default: Azimuth 315° (NW), Altitude 45°
Figure 4.4: Hillshade — illuminated (bright) and shadowed (dark) sides of terrain.

Summary Table

ParameterWhat It MeasuresUnitKey Application
SlopeSteepness / gradientDegrees or %Landslide susceptibility, road design, erosion
AspectCompass direction slope facesDegrees (0°–360°)Solar radiation, crop suitability, building orientation
Contour LinesLines of equal elevationMeters (interval)Topographic mapping, volume calculation, engineering
HillshadeSimulated light/shadow0–255 (grayscale)3D visualization, map backgrounds
Exam Tip: For 10 marks, describe exactly 4 parameters — one sentence definition + brief application. Simple figures for slope and aspect strengthen the answer.

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?"

Key Concept: High flow accumulation = stream/river channel. Flow accumulation of 0 = ridgeline/hilltop (no upstream contributors).

2. Hydrological Workflow

Hydrological Workflow 1. Raw DEM Elevation grid 2. Fill Sinks Remove depressions 3. Flow Direction D8 algorithm 4. Flow Accumulation 5. Streams 6. Watershed Each step depends on the previous — never skip.
Figure 5.1: Flow accumulation is Step 4 in the hydrological workflow.

3. How Flow Accumulation Works (D8 Example)

Flow Accumulation — Step-by-Step Example DEM (Elevation) 787269 746756 695344 Flow Direction (D8) Flow Accumulation 000 022 048 0 = Ridge (no upstream) 2 = 2 upstream cells 4 = 4 upstream cells 8 = ALL 8 cells (outlet) High value = stream channel Cell (3,3) = 8 → All water drains here (outlet).
Figure 5.2: Flow accumulation on a 3×3 DEM grid. Bottom-right cell (44m) accumulates all 8 other cells.

4. Stream Extraction via Threshold

Stream = (Flow Accumulation > Threshold)

Common threshold: 500–1000 cells

ThresholdEffect
Low (e.g., 100)Extracts many small tributaries → denser network
High (e.g., 5000)Only major rivers → sparser network

5. Applications

ApplicationHow Flow Accumulation Helps
Stream network extractionIdentifies where streams form based on accumulated drainage
Watershed delineationCombined with flow direction, defines catchment boundaries
Stream orderingStrahler/Shreve ordering on extracted stream links
Flood risk assessmentHigh accumulation = most vulnerable to flooding
Erosion modelingFlow concentration indicates high erosion potential
Infrastructure planningBridges, culverts, drainage design
Exam Tip: Define flow accumulation, explain prerequisites (Fill → Flow Direction → Flow Accumulation), show grid example, mention threshold-based stream extraction.

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.

Low Pass Filter — 3×3 Mean Kernel Kernel: 1/91/91/9 1/91/91/9 1/91/91/9 Input Neighborhood: 101520 129018 141622 Output center pixel: (10+15+20+12+90+18+14+16+22) / 9 = 24.1 Noisy pixel (90) smoothed down to 24.1 — noise removed!
Figure 6.1: Low pass (3×3 mean) averages all 9 pixels, smoothing the anomalous 90 to 24.1.

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.

High Pass Filter — 3×3 Edge Enhancement Kernel Kernel: -1-1-1 -1-1 -1-1-1 +9 Uniform area (no edge): All pixels = 50 → Output = 50 (50×9 + 50×(-8) = 50, no change) At an edge (sharp change): Center=100, neighbors=50 Output = 100×9 + 50×(-8) = 500 (Edge amplified!) High pass produces large values at edges, near-zero in uniform areas.
Figure 6.2: High pass kernel (+9 center, -1 neighbors). Uniform areas unchanged; edges produce large values.

Comparison: Low Pass vs. High Pass

FeatureLow Pass FilterHigh Pass Filter
Also calledSmoothing / AveragingSharpening / Edge enhancement
Passes throughLow spatial frequenciesHigh spatial frequencies
SuppressesHigh frequencies (sharp changes)Low frequencies (gradual changes)
Effect on imageBlurs / smoothsSharpens / enhances edges
Kernel weightsAll positive, equalMixed (+/-), large positive center
NoiseReduces noiseMay amplify noise
Best forNoise removal, generalizationEdge detection, feature extraction
Exam Tip: For 5 marks each: (1) Define, (2) Explain what passes/suppresses, (3) Show kernel, (4) Give 2–3 applications. Drawing the 3×3 kernel matrix earns full marks.