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

REMOTE SENSING INDICES

Remote sensing indices are band ratios designed to highlight specific surface features (vegetation, soil, water, urban areas, snow, burned areas, etc.) using the spectral reflectance properties of the Earth's surface. They improve classification accuracy and environmental monitoring. 1. Vegetation Indices NDVI – Normalized Difference Vegetation Index Formula: (NIR – RED) / (NIR + RED) Concept: Vegetation reflects strongly in NIR and absorbs in RED due to chlorophyll. Measures: Vegetation greenness & health Uses: Agriculture, drought monitoring, biomass estimation EVI – Enhanced Vegetation Index Formula: G × (NIR – RED) / (NIR + C1×RED – C2×BLUE + L) Concept: Corrects for soil and atmospheric noise. Measures: Vegetation vigor in dense canopies Uses: Tropical rainforest mapping, high biomass regions GNDVI – Green Normalized Difference Vegetation Index Formula: (NIR – GREEN) / (NIR + GREEN) Concept: Uses Green instead of Red ...

Energy Interaction with Atmosphere and Earth Surface

In Remote Sensing , satellites record electromagnetic radiation (EMR) that is reflected or emitted from the Earth. Before reaching the sensor, radiation interacts with: The Atmosphere The Earth's Surface These interactions control how satellite images look and how we interpret them. I. Interaction of EMR with the Atmosphere When solar radiation travels from the Sun to the Earth, four main processes occur: 1. Absorption Definition: Absorption occurs when atmospheric gases absorb radiation at specific wavelengths and convert it into heat. Main absorbing gases: Ozone (O₃) → absorbs Ultraviolet (UV) Carbon dioxide (CO₂) → absorbs Thermal Infrared Water vapour (H₂O) → absorbs Infrared Concept: Atmospheric Windows These are wavelength regions where absorption is very low, allowing radiation to pass through the atmosphere. Remote sensing depends on these windows. For example, satellites like Landsat 8 use visible, near-infrared, and thermal bands located in atmospheric windows. 2. Trans...

Atmospheric Window

The atmospheric window in remote sensing refers to specific wavelength ranges within the electromagnetic spectrum that can pass through the Earth's atmosphere relatively unimpeded. These windows are crucial for remote sensing applications because they allow us to observe the Earth's surface and atmosphere without significant interference from the atmosphere's constituents. Key facts and concepts about atmospheric windows: Visible and Near-Infrared (VNIR) window: This window encompasses wavelengths from approximately 0. 4 to 1. 0 micrometers. It is ideal for observing vegetation, water bodies, and land cover types. Shortwave Infrared (SWIR) window: This window covers wavelengths from approximately 1. 0 to 3. 0 micrometers. It is particularly useful for detecting minerals, water content, and vegetation health. Mid-Infrared (MIR) window: This window spans wavelengths from approximately 3. 0 to 8. 0 micrometers. It is valuable for identifying various materials, incl...

Landsat band composition

Short-Wave Infrared (7, 6 4) The short-wave infrared band combination uses SWIR-2 (7), SWIR-1 (6), and red (4). This composite displays vegetation in shades of green. While darker shades of green indicate denser vegetation, sparse vegetation has lighter shades. Urban areas are blue and soils have various shades of brown. Agriculture (6, 5, 2) This band combination uses SWIR-1 (6), near-infrared (5), and blue (2). It's commonly used for crop monitoring because of the use of short-wave and near-infrared. Healthy vegetation appears dark green. But bare earth has a magenta hue. Geology (7, 6, 2) The geology band combination uses SWIR-2 (7), SWIR-1 (6), and blue (2). This band combination is particularly useful for identifying geological formations, lithology features, and faults. Bathymetric (4, 3, 1) The bathymetric band combination (4,3,1) uses the red (4), green (3), and coastal bands to peak into water. The coastal band is useful in coastal, bathymetric, and aerosol studies because...

Landsat 8 Band designation and Band Combination.

Landsat 8 Band designation and Band Combination.  Landsat 8-9 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS) Bands Wavelength (micrometers) Resolution (meters) Band 1 - Coastal aerosol 0.43-0.45 30 Band 2 - Blue 0.45-0.51 30 Band 3 - Green 0.53-0.59 30 Band 4 - Red 0.64-0.67 30 Band 5 - Near Infrared (NIR) 0.85-0.88 30 Band 6 - SWIR 1 1.57-1.65 30 Band 7 - SWIR 2 2.11-2.29 30 Band 8 - Panchromatic 0.50-0.68 15 Band 9 - Cirrus 1.36-1.38 30 Band 10 - Thermal Infrared (TIRS) 1 10.6-11.19 100 Band 11 - Thermal Infrared (TIRS) 2 11.50-12.51 100 Vineesh V Assistant Professor of Geography, Directorate of Education, Government of Kerala. https://www.facebook.com/Applied.Geography http://geogisgeo.blogspot.com