SAR change detection with SNAP
From AWF-Wiki
(Difference between revisions)
(→Data transfer) |
(→Display in QGIS) |
||
| (27 intermediate revisions by one user not shown) | |||
| Line 1: | Line 1: | ||
| − | In this tutorial we will use Level-1 Ground Range Detected (GRD) Sentinel-1 data downloaded from the Copernicus Open Access Hub. The example are fire events during | + | In this tutorial we will use Level-1 Ground Range Detected (GRD) Sentinel-1 data downloaded from the Copernicus Open Access Hub. The example are fire events during 2019 in a peat forest located in Eastern Jambi, Sumatra, Indonesia. |
=Data transfer= | =Data transfer= | ||
| − | Download | + | Download example data from STud.IP or use data from your personal region of interest: two Sentinel-1 scenes one pre-event acquisition date and one post-event date. Follow and adapt accordingly the instruction on [[Downloading Sentinel-2 images]]. |
Sentinel-1 data can have different characteristics (product type, polarization, sensor mode and ascending/descending orbit). Make sure, that you have the same characteristics for both images. | Sentinel-1 data can have different characteristics (product type, polarization, sensor mode and ascending/descending orbit). Make sure, that you have the same characteristics for both images. | ||
Set the sensing period (e.g. from 01.06.2019 to 30.11.2019). | Set the sensing period (e.g. from 01.06.2019 to 30.11.2019). | ||
| Line 13: | Line 13: | ||
* Open SNAP Desktop. | * Open SNAP Desktop. | ||
* Go to {{mitem|text=File → Import → SAR sensors → Sentinel-1}} and select the downloaded zip files. | * Go to {{mitem|text=File → Import → SAR sensors → Sentinel-1}} and select the downloaded zip files. | ||
| − | * | + | * In '''Product Explorer''' click on + to unfold the product folders. Mark the '''Amplitude_VV''' bands of both products. |
* Right click '''Open 2 Image Window'''. | * Right click '''Open 2 Image Window'''. | ||
| − | * {{mitem|text=Windows → Tile Horizontally}}. Both Amplitude bands are shown side by side geometrically linked to each other. | + | * {{mitem|text=Windows → Tile Horizontally}}. Both Amplitude bands are shown side by side geometrically linked to each other. Make sure the option {{mitem|text=View --> Synchronise Image View}} and {{mitem|text=View --> Synchronise Image Cursors}}is activated). |
| + | * Adjust the viewer size that they are evenly. | ||
| + | * Zoom to an area of interest where you may detect changes (see below). Compare the pre-event (2019-06-06) and post-event image (2019-11-21). | ||
[[File:snap_sidebyside.png|800px]] | [[File:snap_sidebyside.png|800px]] | ||
=Subsetting= | =Subsetting= | ||
| − | + | Subset the Sentinel-1 scenes reducing the data volume. | |
| − | * | + | * Activate the pre-event image clicking on the tab of the viewer. The active image is marked with a yellow frame. |
| − | + | ||
| − | + | ||
[[File:SNAP_Select_Subset_Area.PNG|600px]] | [[File:SNAP_Select_Subset_Area.PNG|600px]] | ||
| − | * | + | * Right click into the active image and select '''Spatial Subset from View...'''. |
| − | + | ||
[[File:SNAP_Subset.PNG|400px]] | [[File:SNAP_Subset.PNG|400px]] | ||
| − | + | * Besides a spatial subset you can also select a band subset with this tool. Go to the tab {{button|text=Band subset}} and select only one band '''Amplitude_VV'''. | |
| − | + | * Click {{button|text=OK}}. | |
| − | ''' | + | * A new Product subset_* can be found in the Product Explorer. |
| + | + Activate the prost-event image and repeat the work steps above. | ||
| + | ''Note, the resulting images are still in acquisition geometry, whereby the images appear mirrored.'' | ||
=Radiometric and geometric calibration= | =Radiometric and geometric calibration= | ||
=== Apply Orbit File === | === Apply Orbit File === | ||
| − | * Select the | + | * Select the subset '''pre-event''' product from the previous step within the '''Product Explorer'''. |
| − | * Go to {{mitem|text=Radar → Apply Orbit File}} | + | * Go to {{mitem|text=Radar → Apply Orbit File}} |
* Default settings are fine. | * Default settings are fine. | ||
* Select target product directory. | * Select target product directory. | ||
| Line 46: | Line 47: | ||
=== Radiometric Calibration === | === Radiometric Calibration === | ||
* Select the product from the previous step within the '''Product Explorer''' on the left side. | * Select the product from the previous step within the '''Product Explorer''' on the left side. | ||
| − | * Go to {{mitem|text=Radar → Radiometric | + | * Go to {{mitem|text=Radar → Radiometric → Calibrate}} |
* Select target product directory. | * Select target product directory. | ||
* Go to the '''Processing Parameters''' tab and check sigma0. (You could also use gamma0 processing.) | * Go to the '''Processing Parameters''' tab and check sigma0. (You could also use gamma0 processing.) | ||
| Line 55: | Line 56: | ||
* Select the product from the previous step within the '''Product Explorer''' on the left side. | * Select the product from the previous step within the '''Product Explorer''' on the left side. | ||
* For speckle filtering, there are many different options of filter algorithms and window size. | * For speckle filtering, there are many different options of filter algorithms and window size. | ||
| − | * Use the Lee filter | + | * Use the Lee filter with a window size of 3x3. |
* Go to {{mitem|text=Radar → Speckle Filtering → Single Product Speckle Filter}} | * Go to {{mitem|text=Radar → Speckle Filtering → Single Product Speckle Filter}} | ||
* Select target product directory. | * Select target product directory. | ||
| Line 72: | Line 73: | ||
* Now the image should be correctly displayed, in terms of geometry, and be smoothed. | * Now the image should be correctly displayed, in terms of geometry, and be smoothed. | ||
| − | Repeat | + | '''Repeat the radiometric and geometric calibration steps for the post-event subset (second acquisiton date)!''' |
=Change composite= | =Change composite= | ||
=== Create Stack === | === Create Stack === | ||
* Close all products by marking them all '''Product Explorer''' on the left side, right click and '''Close All Products'''. | * Close all products by marking them all '''Product Explorer''' on the left side, right click and '''Close All Products'''. | ||
| − | * | + | * {{mitem|text=File --> Open Product}}. Load both radiometric and geometric corrected images from the previous step with the extension '''_Orb_Cal_Spk_TC.dim'''. |
* Go to {{mitem|text=Radar → Coregistration → Stack Tools → Create Stack}} | * Go to {{mitem|text=Radar → Coregistration → Stack Tools → Create Stack}} | ||
| − | * Click on | + | * Click on [[File:snap_add_open_button.png]] button on the right side to select the images. |
* Leave the default settings at the '''2-CreateStack''' tab. | * Leave the default settings at the '''2-CreateStack''' tab. | ||
* Go to the '''3-Write''' tab and select target product directory. | * Go to the '''3-Write''' tab and select target product directory. | ||
| Line 87: | Line 88: | ||
=== Convert Datatype === | === Convert Datatype === | ||
* Select the product from the previous step within the '''Product Explorer''' on the left side. | * Select the product from the previous step within the '''Product Explorer''' on the left side. | ||
| − | * Go to {{mitem|text=Raster | + | * Go to {{mitem|text=Raster → Data Conversion → Convert datatype}} |
* Save as {{button|text=GeoTIFF}}. | * Save as {{button|text=GeoTIFF}}. | ||
* Select target product directory. | * Select target product directory. | ||
| Line 98: | Line 99: | ||
* {{mitem|text=Data Source Manger --> Browser --> XYZ Tiles --> Google Satellite}} | * {{mitem|text=Data Source Manger --> Browser --> XYZ Tiles --> Google Satellite}} | ||
* Load the change compoisite file with extension '''_Cnv.tif''' {{mitem|text=Data Source Manger --> Raster}}. | * Load the change compoisite file with extension '''_Cnv.tif''' {{mitem|text=Data Source Manger --> Raster}}. | ||
| − | * {{mitem|text=Layer properties --> Symbology}}. Adjust the band assignment: Red = Band 2, Green = Band 1, Blue = Band 1. | + | * {{mitem|text=Layer properties --> Symbology}}. Adjust the band assignment: Red = Band 2, Green = Band 1, Blue = Band 1. |
* Adjust opacity of the change layer on top of Google Satellite. {{mitem|text=Layer Properties --> Transparency --> Global Transparency}}. | * Adjust opacity of the change layer on top of Google Satellite. {{mitem|text=Layer Properties --> Transparency --> Global Transparency}}. | ||
In the bi-temporal change composite red indicates change areas. | In the bi-temporal change composite red indicates change areas. | ||
Latest revision as of 14:16, 25 January 2021
In this tutorial we will use Level-1 Ground Range Detected (GRD) Sentinel-1 data downloaded from the Copernicus Open Access Hub. The example are fire events during 2019 in a peat forest located in Eastern Jambi, Sumatra, Indonesia.
Contents |
[edit] Data transfer
Download example data from STud.IP or use data from your personal region of interest: two Sentinel-1 scenes one pre-event acquisition date and one post-event date. Follow and adapt accordingly the instruction on Downloading Sentinel-2 images. Sentinel-1 data can have different characteristics (product type, polarization, sensor mode and ascending/descending orbit). Make sure, that you have the same characteristics for both images. Set the sensing period (e.g. from 01.06.2019 to 30.11.2019). Select Product type GRD and sensor mode IW. Download the images of the 06.06.2019 and 21.11.2019 (both GRD, VV/VH, IW, Descending).
[edit] Visualization of amplitude bands
- Open SNAP Desktop.
- Go to File → Import → SAR sensors → Sentinel-1 and select the downloaded zip files.
- In Product Explorer click on + to unfold the product folders. Mark the Amplitude_VV bands of both products.
- Right click Open 2 Image Window.
- Windows → Tile Horizontally. Both Amplitude bands are shown side by side geometrically linked to each other. Make sure the option View --> Synchronise Image View and View --> Synchronise Image Cursorsis activated).
- Adjust the viewer size that they are evenly.
- Zoom to an area of interest where you may detect changes (see below). Compare the pre-event (2019-06-06) and post-event image (2019-11-21).
[edit] Subsetting
Subset the Sentinel-1 scenes reducing the data volume.
- Activate the pre-event image clicking on the tab of the viewer. The active image is marked with a yellow frame.
- Right click into the active image and select Spatial Subset from View....
- Besides a spatial subset you can also select a band subset with this tool. Go to the tab Band subset and select only one band Amplitude_VV.
- Click OK.
- A new Product subset_* can be found in the Product Explorer.
+ Activate the prost-event image and repeat the work steps above. Note, the resulting images are still in acquisition geometry, whereby the images appear mirrored.
[edit] Radiometric and geometric calibration
[edit] Apply Orbit File
- Select the subset pre-event product from the previous step within the Product Explorer.
- Go to Radar → Apply Orbit File
- Default settings are fine.
- Select target product directory.
- Click Run.
[edit] Radiometric Calibration
- Select the product from the previous step within the Product Explorer on the left side.
- Go to Radar → Radiometric → Calibrate
- Select target product directory.
- Go to the Processing Parameters tab and check sigma0. (You could also use gamma0 processing.)
- Click Run.
[edit] Single Product Speckle Filter
- Select the product from the previous step within the Product Explorer on the left side.
- For speckle filtering, there are many different options of filter algorithms and window size.
- Use the Lee filter with a window size of 3x3.
- Go to Radar → Speckle Filtering → Single Product Speckle Filter
- Select target product directory.
- Go to the Processing Parameters tab and select Lee in the Filter drop down menue.
- Adapt the window size in X and Y direction accordingly.
- Click Run.
[edit] Range-Doppler-Terrain Correction
- Select the product from the previous step within the Product Explorer on the left side.
- Go to Radar → Geometric → Terrain Correction → Range-Doppler-Terrain Correction
- Select target product directory.
- Go to the Processing Parameters tab and select SRTM 1sec HGT (Auto Download) in the Digital Elevation Model drop down menue.
- Click Run.
- Now the image should be correctly displayed, in terms of geometry, and be smoothed.
Repeat the radiometric and geometric calibration steps for the post-event subset (second acquisiton date)!
[edit] Change composite
[edit] Create Stack
- Close all products by marking them all Product Explorer on the left side, right click and Close All Products.
- File --> Open Product. Load both radiometric and geometric corrected images from the previous step with the extension _Orb_Cal_Spk_TC.dim.
- Go to Radar → Coregistration → Stack Tools → Create Stack
- Click on
button on the right side to select the images.
- Leave the default settings at the 2-CreateStack tab.
- Go to the 3-Write tab and select target product directory.
- Click Run.
[edit] Convert Datatype
- Select the product from the previous step within the Product Explorer on the left side.
- Go to Raster → Data Conversion → Convert datatype
- Save as GeoTIFF.
- Select target product directory.
- Default settings are fine.
- Click Run.
[edit] Display in QGIS
- Open QGIS.
- Data Source Manger --> Browser --> XYZ Tiles --> Google Satellite
- Load the change compoisite file with extension _Cnv.tif Data Source Manger --> Raster.
- Layer properties --> Symbology. Adjust the band assignment: Red = Band 2, Green = Band 1, Blue = Band 1.
- Adjust opacity of the change layer on top of Google Satellite. Layer Properties --> Transparency --> Global Transparency.
In the bi-temporal change composite red indicates change areas.
