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

Geologic and tectonic framework of the Indian shield

  Major Terms and Regions Explained 1. Indian Shield The Indian Shield refers to the ancient, stable core of the Indian Plate made of hard crystalline rocks. It comprises Archean to Proterozoic rocks that have remained tectonically stable over billions of years. Important Geological Features and Regions ▪️ Ch – Chhattisgarh Basin A sedimentary basin part of the Bastar Craton . Contains rocks of Proterozoic age , mainly sedimentary. Important for understanding the evolution of central India. ▪️ CIS – Central Indian Shear Zone A major tectonic shear zone , separating the Bundelkhand and Bastar cratons . It records intense deformation and metamorphism . Acts as a suture zone , marking ancient tectonic collisions. ▪️ GR – Godavari Rift A rift valley formed due to stretching and thinning of the Earth's crust. Associated with sedimentary basins and hydrocarbon resources . ▪️ M – Madras Block An Archean crustal block in...

Evaluation and Characteristics of Himalayas

Time Period Event / Process Geological Evidence Key Terms & Concepts Late Precambrian – Palaeozoic (>541 Ma – ~250 Ma) India part of Gondwana , north bordered by Cimmerian Superterranes, separated from Eurasia by Paleo-Tethys Ocean . Pan-African granitic intrusions (~500 Ma), unconformity between Ordovician conglomerates & Cambrian sediments. Gondwana, Paleo-Tethys Ocean, Pan-African orogeny, unconformity, granitic intrusions, Cimmerian Superterranes. Early Carboniferous – Early Permian (~359 – 272 Ma) Rifting between India & Cimmerian Superterranes → Neotethys Ocean formation. Rift-related sediments, passive margin sequences. Rifting, Neotethys Ocean, passive continental margin. Norian (210 Ma) – Callovian (160–155 Ma) Gondwana split into East & West; India part of East Gondwana with Australia & Antarctica. Rift basins, oceanic crust formation. Continental breakup, East Gondwana, West Gondwana, oceanic crust. Early Cretaceous (130–125 Ma) India broke fr...

Seismicity and Earthquakes, Isostasy and Gravity

1. Seismicity and Earthquakes in the Indian Subcontinent Key Concept: Seismicity Definition : The occurrence, frequency, and magnitude of earthquakes in a region. In India, seismicity is high due to active tectonic processes . Plate Tectonics 🌏 Indian Plate : Moves northward at about 5 cm/year. Collision with Eurasian Plate : Causes intense crustal deformation , mountain building (Himalayas), and earthquakes. This is an example of a continental-continental collision zone . Seismic Zones of India Classified into Zone II, III, IV, V (Bureau of Indian Standards, BIS). Zone V = highest hazard (e.g., Himalayas, Northeast India). Zone II = lowest hazard (e.g., parts of peninsular India). Earthquake Hazards ⚠️ Himalayas: prone to large shallow-focus earthquakes due to active thrust faulting. Northeast India: complex subduction and strike-slip faults . Examples: 1897 Shillong Earthquake (Magnitude ~8.1) 1950 Assam–Tib...

Vector geoprocessing - Clipping, Erase, identify, Union & Intersection

Think of your vector data (points, lines, polygons) like shapes drawn on a transparent sheet. Geoprocessing is just cutting, joining, or comparing those shapes to get new shapes or information. 1. Clipping ✂️ Imagine you have a big map and you only want to keep a part of it (like cutting a photo into a smaller rectangle). You use another shape (like the boundary of a district) to "clip" and keep only what is inside. Result: Only the data inside the clipping shape remains. 2. Erase 🚫 Opposite of clipping. You remove (erase) the area of one shape from another shape. Example: You have a city map and want to remove all the park areas from it. 3. Identify 🔍 This checks which features from one layer fall inside (or touch) another layer. Example: Identify all the schools inside a flood zone. 4. Union 🤝 Combines two shapes together and keeps everything from both. Works like stacking two transparent sheets and redrawing t...

vector data analysis in GIS Surface Analysis – Interpolation – IDW

1. Surface Analysis 🗺️ This is when we try to understand and visualize how a value changes across a surface (like land). The values might be temperature, rainfall, elevation, pollution levels, etc. We often start with only some points where we know the value, but we want to guess the values everywhere in between. 2. Interpolation 📍➡️📍 Interpolation is a way of estimating unknown values between known points. Imagine you know the temperature at a few weather stations, but you want to know the temperature everywhere in between. GIS uses math to "fill in the blanks" between the points. 3. IDW (Inverse Distance Weighted) 🎯 One popular interpolation method. The idea: Points that are closer to you have more influence than points farther away. Example: If you're standing between two rain gauges, the closer one's reading will affect your estimated rainfall more than the farther one. "Inverse Distance" means: The ...