I am very interested in studying this radar. But after I read the article and I tried to imitate the way of processing but I was stopped at the point of processing in the snap alone, maybe because I do not know or do not understand the way of advanced processing using the formula in the article. Maybe from you can help me overcome the difficulty ???
retrieving soil moisture from radar is mostly not doable with a single formula. Please have a look at these considerations: Soil moisture mapping using Sentinel-1-Data and SNAP? and After Sigma0
You cannot derive soil moisture using just a single SAR image. Most of the studies of soil moisture retrieval conducted using time series data (as described in the article I have posted above). This procedure is called
change detection and it is widely used.
To derive soil moisture from SAR, to be more precise,
relative soil wetness you can do the following steps:
Data collection. A series of SAR data collection within a period of time
Image pre-processing: Each SAR image should be pre-processed seperately (speckle noise suppression, radiometric calibration, incidence angle normalization if needed and terrain correction).
Change detection: To perform change detection you need to subtract a SAR image acquired under dry weather conditions from a SAR image acquired under wet weather conditions. You should follow the equations described in the article i have posted.
These equations are straightforward to be applied. SO, you can use the
band math in SNAP to apply the equations on your image.
Things to consider
Relative soil wetness from SAR is not possible if your AOI is covered by a thick vegetation. The microwave signal interacts with the leaves and doe not penetrate through to reach the soil.
NDVI index should be calculated using optical data. Then, this layer can be used to mask out the vegetated area in your image.
I hope that hekps
Usefull information @johngan
Thanks you sir.
I have collected radar images some time in real time as you say. But I can only process the radar image until I get the backscatter value. And for the next stage I do not know how to process from backscatter to get a raster image or soil moisture map that contains soil moisture values???. Maybe mr can help me.
And than, I apologize in advance because I am not an expert in processing radar or optical image data, because I am not among geophysicists or geomatics.
First of all, there are quite a few methods developed for deriving soil moisture. I have used the
change detection method. A time series of SAR images is require (you have already collected the images as you mentioned). In terms of pre-processing the SAR images, Figure 1 describes the steps that need to be followed.
Figure 1. SAR image preprocessing
Once each image is pre-processed, then you can apply the equations described in the papers below to derive the
Those papers are very good and describe exactly what you need to do in order to derive soil wetness from SAR.
A summary of steps to be followed is given below:
SAR Data collection*
Time series data collection(GRD) within a time window
Pre-process each individual SAR image as shown in figure 1
Mask out areas for each SAR image with NDVI index above 0.6. (NDVI index is not the best measure though as it does not tell us anything about how tall trees are or the shape of trees canopy.)
Take the average of all the masked time-series SAR images
Soil mositure retrieval
Dry amd Wet reference determination.For change detection method, we need an image acquired under wet weather conditions and an image acquired under dry weather conditions.
Two-tailed test.To derive wet and dry reference, we perform a two-tailed test using z-score. Values that falls beyond 1.96 (which represents the 97.5 percentile) are the extreme high values, while values below -1.96 are the extreme values. Then, we take the average for both high and low values and we have our wet and dry reference
Soil moisture retrieval.Finally, soil moisture retrieval can be derived using the equations shown at the articles
This is a brief summary of what you need to do.
I hope that helps
thank you sir, for the information and knowledge that you gave me,
and hopefully my curiosity about this soil moisture measurement can be answered
Hello！In the paper you provided, I know how to get the soil moisture index, but I am not quite sure what is the relationship between SMI and soil moisture content? Thank you in advance！
Soil moisture content refers to the quantity of water contained in a material, such as soil (called soil moisture). Water content is expressed as a ratio, which ranges from zero (completely dry) up to the value of material’s saturation. Water content is measured in m^3 (volume). In order to derive soil moisture content, we need to know some parameters (e.g volumetric water content) which can be derived from in situ measurements .
Using earth observation (SAR data) techniques, we cannot derive soil moisture content (unless we incorporate in site measurements in our analysis). We can derive soil moisture index (soil wetness) which tells us how much wet or dry the soil is.
If we have in site measured data, we can get the soil moisture data through regression analysis.
you can import the measurements as points (by CSV or Shapefile) and use the correlative plot view to build a regression.
The results of soil moisture depends on the size of your AOI and the number of soil moisture sensors.
If you have to study an area, let’s say 100km x 100km using satellite data, you are gonna need a network of sensors scattered around your AOI (might be 100, 200, e.t.c). These sensors will provide you with accurate measurements of water content and their geographic location using GPS equipment.
In order though to estimate the the soil moisture content for the whole satellite scene, you will probably need to apply machine learning algorithms (Neural Network, SVM, e.t.c) in order to make predictions of the soil moisture content for the whole satellite scene.
I agree. There is no simple relationship. But applying ML would require several spatial input parameters. One single SAR image is not sufficient. I therefore like the idea discussed here to calculate an index value (based on multi-temporal SAR data) which can be directly related to the measured moisture variations.
yes agree, on SAR image is not sufficient. We would need other auxiliary data (optical data potentially as well)
the simplest way to do it is the one you described.
Means “wet weather conditions“ that it must rain while the acquisition is ongoing or short before the acuisition? I thought precipitation while radar meassurements leads to unwanted noise?
I’d say shortly before the image acquisition would be the ideal case, but certainly not during the measurement.
Not sure you are aware of the Copernicus product from the copernicus services on Soil moisture from Sentinel-1 that will be soon coming.
@johngan Thanks for your detailed replies. I’m implementing the change detection technique, as you described in your post from Oct 18 however on Sentinel 1 multi temporal imagery, not ASAR WS imagery as the paper describes.
A few questions I have are,
- is the VH polarization of S1 more appropriate than VV for change detection?
- is it necessary to compute angular correction for S1? since each pixel has a characteristic local angle that is unchanged throughout the year over my study area.
- lastly, is it possible to produce a high resolution (10 m) soil moisture index map using Sentinel 1 imagery? the paper you cited only mentions spatial resolution once, i.e. 75 m ASAR pixels are averaged to a 1 km grid to allow low sensitivity pixels to be processed. In my case multi-temporal imagery (3 yrs) gives me a high sensitivity > 3 dB over the arable and barren land of my study area, hence is it necessary to average the resolution of S1 images or can I retain the 10m resolution?
The problem with VV polarization is that the vertically polarized microwaves are affected by vegetation (HH pol is less sensitive to vegetation) and the penetration through soil is minimized. VV polarization can detect the differenceswithin the vertical vegetation structure from various growth stages. As M.R.Saradjian and M.Hosseini, 2011 mentioned, after experiments, HV cross-polarization proved to be more accurate than HH or VV polarization for areas with some vegetation.
The backscatter coefficient as a function of incidence angle declines. The larger the incidence angle the lower the sensitivity of the microwave. Hence, we can identify larger errors in our measurements at large incidence angles due to the low sensitivity. It is recommended to perform an incidence angle normalization at approximately 30 degrees.
You can perform soil moisture analysis using Sentinel-1 data. there is no need to down-sample your 20m resolution to a coarser one
Thanks @johngan for that reply
to your second point, incidence angle normalization is necessary to correctly analyse spatial distribution of soil moisture in areas where topography varies? Is it still important if I am only concerned with per pixel changes in soil moisture over time.
Further I looked into the forums regarding incidence angle normalization, specifically these posts,  and  and this paper my interpretation is that the normalization process gives us gamma-naught which is sigma-naught/cos(local incidence angle). is this correct?