Skip to main content

Kriging in GIS and variogram

Kriging is an advanced spatial interpolation technique used in GIS (Geographic Information System) that estimates values for unknown locations based on the values observed at nearby known locations. It is a geostatistical method that takes into account not only the distances between points but also the spatial correlation or variability in the data.

Unlike simpler interpolation methods like IDW, which assume a constant variation across the study area, kriging incorporates the spatial autocorrelation of the data to produce more accurate and precise estimates. Kriging considers the spatial arrangement and patterns of the data points to generate a surface that honors the underlying spatial structure.

The key principle behind kriging is the variogram, which quantifies the spatial correlation between pairs of points at different distances. The variogram measures how the values of nearby points vary from each other as a function of distance. It provides information about the spatial dependence or variability in the dataset.

The kriging process involves three main steps:

1. Variogram modeling: The first step in kriging is to construct a variogram, which is a plot of the semivariance (a measure of dissimilarity or variability) against distance or lag between pairs of points. The variogram helps to understand the spatial structure of the data and determine the range, sill, and nugget effect. Based on the variogram, a mathematical model is fitted to describe the spatial correlation.

2. Interpolation: Once the variogram is modeled, kriging calculates the weights or coefficients for the known points based on their spatial relationship to the target location. The weights are determined through a process known as kriging equations, which consider the variogram and covariance between points. These equations generate the optimal weights that minimize the prediction error.

   - Ordinary Kriging (OK): Assumes a constant mean value across the study area.
   - Simple Kriging (SK): Accounts for an unknown mean value, estimating it from the data.
   - Universal Kriging (UK): Incorporates additional spatially correlated variables (covariates) in addition to the location coordinates.

3. Prediction: The final step is the estimation of values at the unknown locations using the calculated weights. Kriging provides not only the predicted values but also the estimation error or uncertainty associated with each prediction. This information can be valuable in decision-making processes.

Kriging is particularly useful when dealing with spatial datasets that exhibit spatial autocorrelation, anisotropy (directional dependence), or irregularly spaced points. It provides a more sophisticated approach to spatial interpolation by considering the inherent spatial relationships in the data.

GIS software typically provides various kriging algorithms and tools that allow users to model the variogram, perform the interpolation, and generate kriging predictions and associated error maps.

Comments

Popular posts from this blog

Quantitative Revolution in Geography

The Quantitative Revolution was a major change in Geography during the 1950s and 1960s. It introduced mathematical, statistical and scientific methods to understand geographical phenomena. Before this revolution, geographers mainly described places, regions, landscapes and differences between areas. After the revolution, they increasingly used numerical data, statistical analysis, mathematical models and theories to explain why geographical patterns occur and how they may change. For example: Traditional approach: Describe the distribution of population in different regions. Quantitative approach: Use statistical methods to identify population patterns, explain the factors influencing population distribution and predict future population growth. The revolution helped establish Geography as a spatial science, focusing on the location, distribution, relati...

Satellite Remote Sensing: Concepts and Imaging Systems

  Satellite remote sensing relies on detectors (sensors) that measure reflected/emitted electromagnetic radiation from the Earth. How sensors collect data depends on: Spectral coverage → Multispectral (few bands), Hyperspectral (hundreds), Thermal, Microwave. Detector type → Discrete detectors, linear arrays, or area arrays. Scanning mechanism → Scanning mirrors (whiskbroom) vs. linear arrays (pushbroom). Multispectral Imaging Using Discrete Detectors and Scanning Mirrors (Whiskbroom scanners) Principle : A single detector (or a few detectors) measures radiation one pixel at a time. A scanning mirror sweeps across-track (perpendicular to the satellite path) to build up the image line by line. The forward motion of the satellite provides the along-track dimension. Known as a whiskbroom scanner . Characteristics : Good calibration stability. Narrow instantaneous field of view (IFOV). Susceptible to mechanical wear (moving ...

Spectral Signature vs. Spectral Reflectance Curve

Spectral Signature  A spectral signature is the unique pattern in which an object: absorbs energy reflects energy emits energy across different wavelengths of the electromagnetic spectrum. ✔ Key Points Every natural and man-made object on Earth interacts with sunlight differently. These interactions produce a distinct pattern , just like a "fingerprint". Sensors on satellites record these patterns as digital numbers (DN values) . These patterns help to identify and differentiate objects such as vegetation, soil, water, snow, buildings, minerals, etc. ✔ Examples of Spectral Signatures Healthy vegetation → High reflectance in NIR , strong absorption in red Water → Strong absorption in NIR and SWIR , low reflectance Dry soil → Gradual increase in reflectance from visible to NIR Snow → High reflectance in visible , low in SWIR ✔ Why Spectral Signature Matters It allows: Land cover classification Chan...

Geometric Correction

When satellite or aerial images are captured, they often contain distortions (errors in shape, scale, or position) caused by many factors — like Earth's curvature, satellite motion, terrain height (relief), or the Earth's rotation . These distortions make the image not properly aligned with real-world coordinates (latitude and longitude). 👉 Geometric correction is the process of removing these distortions so that every pixel in the image correctly represents its location on the Earth's surface. After geometric correction, the image becomes geographically referenced and can be used with maps and GIS data. Types  1. Systematic Correction Systematic errors are predictable and can be modeled mathematically. They occur due to the geometry and movement of the satellite sensor or the Earth. Common systematic distortions: Scan skew – due to the motion of the sensor as it scans the Earth. Mirror velocity variation – scanning mirror moves at a va...

RADIOMETRIC CORRECTION

  Radiometric correction is the process of removing sensor and environmental errors from satellite images so that the measured brightness values (Digital Numbers or DNs) truly represent the Earth's surface reflectance or radiance. In other words, it corrects for sensor defects, illumination differences, and atmospheric effects. 1. Detector Response Calibration Satellite sensors use multiple detectors to scan the Earth's surface. Sometimes, each detector responds slightly differently, causing distortions in the image. Calibration adjusts all detectors to respond uniformly. This includes: (a) De-Striping Problem: Sometimes images show light and dark vertical or horizontal stripes (banding). Caused by one or more detectors drifting away from their normal calibration — they record higher or lower values than others. Common in early Landsat MSS data. Effect: Every few lines (e.g., every 6th line) appear consistently brighter or darker. Soluti...