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

Radar Remote Sensing SAR

1. Radar and Radar Remote Sensing RADAR stands for Radio Detection and Ranging . It is an active remote sensing system that transmits microwave energy toward the Earth's surface and records the energy that is returned to the sensor as an echo or backscatter . Unlike passive optical remote sensing, radar does not depend on sunlight. Therefore, it can operate day and night and, at suitable wavelengths, can acquire data through clouds, haze and light rain. Hence, radar is widely described as an all-weather, day-and-night remote sensing technology . Basic principle Microwave pulse → transmission → interaction with surface → backscatter/echo → antenna receives signal → signal processing → radar image 2. Microwave Energy Radar systems use microwave electromagnetic radiation , generally in wavelength ranges from approximately 1 mm to 1 m . Important radar bands include: Band Approx. wavelength Common applications ...

Lidar

LiDAR (Light Detection and Ranging) is an active remote sensing technology that measures distances by illuminating a target with laser pulses and analyzing the time it takes for the reflected light to return. Unlike passive systems (e.g., cameras, multispectral sensors), LiDAR provides its own energy source (laser), allowing it to operate both day and night and even penetrate through vegetation canopies . 🔹 How LiDAR Works (Step-by-Step Process) Laser Pulse Emission The system emits rapid, short pulses of laser light (commonly in the near-infrared wavelength, 1064 nm ). Some systems emit up to hundreds of thousands of pulses per second . Interaction with Target Surface The laser beam strikes objects such as vegetation, buildings, or bare ground. Depending on the object's structure, part of the pulse may scatter or reflect. Return Signal Detection The sensor records multiple returns : First Return → typically vegetation canopy tops. ...

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...

Radar image. Polarization in Remote Sensing

L band radars operate on a wavelength of 15-30 cm and a frequency of 1-2 GHz. L band radars are mostly used for clear air turbulence studies. S band radars operate on a wavelength of 8-15 cm and a frequency of 2-4 GHz. Because of the wavelength and frequency, S band radars are not easily attenuated. . Polarization refers to the direction of travel of an electromagnetic wave vector's tip: vertical (up and down), horizontal (left to right), or  circular (rotating in a constant plane left or right). . a synthetic aperture radar (SAR) for high-resolution imaging. a radar altimeter, to measure the ocean topography. echo amplitude a wind scatterometer to measure wind speed and direction. Other types of radars have been flown for Earth observation missions: precipitation radars such as the  Tropical Rainfall Measuring Mission, or cloud radars like the one used on Cloudsat. . RISAT-1 (SAR, ISRO India, 2012) RORSAT (SAR, Soviet Union, 1967-1988) Seasat (SAR, altimeter, scatterometer, US, 19...

Image Transformation

In remote sensing, image transformation is used to create new images by manipulating existing ones, enhancing the information extracted from each pixel. Multi-image manipulation and other transformations allow analysts to highlight specific features, compare time-series data, or reduce data dimensions for better interpretation. 1. Multi-Image Manipulation This technique involves applying mathematical operations to multiple images, typically from different bands, dates, or sensors, to highlight specific features or changes. Common methods include: a. Image Addition Purpose : Highlights areas of similarity or enhances features when the same objects appear in multiple images. How it Works : By adding pixel values from two images, the output image emphasizes areas where brightness is high in both. Example : Enhancing overall brightness or detecting seasonal changes by adding images from different times. b. Image Subtraction Purpose : Identifies differences between images, useful for chan...