Marble/Sentinel2MapTheme: Difference between revisions

From KDE TechBase
(Page moved to community wiki)
 
(90 intermediate revisions by 3 users not shown)
Line 1: Line 1:
= Satellite Map Theme based on the Sentinel 2 Mission =
{{ Moved To Community }}
 
== Introduction/Motivation ==
 
This page is about the process that allows to create a satellite map theme for Marble based on the Sentinel data.
 
== Sentinels Scientific Data Hub ==
 
=== Registering at the SciHub ===
 
Go to the [https://scihub.copernicus.eu/dhus/#/self-registration Registration website of the SciHub]  and enter your name, username, password, e-Mail and country. As a domain you can choose "Land" and "Research" for the usage.
 
Once you press the '''Register''' button you will receive an e-Mail which allows you to finalize your registration in the SciHub system by pressing the included link.
 
As an alternative to creating your own login you can use "guest/guest" as a username/password when logging in.
 
=== Looking for data ===
 
Go to the [https://scihub.copernicus.eu/dhus/#/home SciHub website] and press the '''Login Button''' to enter it.
On the right side you can press the layers button to choose between the Hybrid, Road (OpenStreetMap) and Satellite map (Blue Marble/LandSat) for orientation.
 
Now we are ready to query some actual Sentinel data. In order to limit our query to useful entries we recommend these adjustments: Press the menu button on the left of the query line edit to enter the Advanced Search
 
[[File:ss1.png]]
 
* Check '''Mission: Sentinel-2'''
* Enter '''[0 TO 10]''' in the cloud coverage Section
 
Press the menu button again. The button will turn into a filter icon. Then enter a '''country name''' (e.g. Germany). As a result of the query you will get a list of entries.
 
[[File:ss2.png]]
 
You can browse through the pages. The '''thumb nails''' give you an idea about the quality of the photographic images. You can press the "eye" button to get an enlarged version of the thumbnail and to get more information.
 
* images that contain an '''area that hasn't been included in the map theme''' yet
* images that ideally have '''no cloud cover'''
* images that were taken in '''bright light'''.
* images that includes a maximum amount of complete tiles (with a minimum of tilted cut off black portions).
 
Once you've made a choice we suggest that you write down:
 
* the '''id''' of the entry (which you can retrieve by pressing the "eye" icon for the entry and copy the id from the title line.
* and the '''download url'''
 
in a text file. Good examples are:
 
id:    S2A_OPER_PRD_MSIL1C_PDMC_20160508T220612_R008_V20160508T104027_20160508T104027
url:  https://scihub.copernicus.eu/dhus/odata/v1/Products('4658007a-76a0-455c-bb9a-d9e789b5e1e3')/$value
 
id:    S2A_OPER_PRD_MSIL1C_PDMC_20160508T221741_R008_V20160508T104027_20160508T104027
url:  https://scihub.copernicus.eu/dhus/odata/v1/Products('0ded3149-d971-4524-ab80-536ffc1cbb55')/$value
 
Now click the download link to download the tile. The download will amount to 4-7GB on average. Even if the download gets aborted in the browser it should be possible to continue the download later on. Depending on the connection the download might take a few hours. At times the SciHub server has temporary maintenance down times.
 
Once the zip file is download '''unzip''' it in a directory of your choice. Make sure that your hard disc still has enough space left (approx. 10 GB).
 
== Installing QGIS and plugins ==
 
=== QGIS installation ===
 
In order to process the data you need to install QGIS.
On Ubuntu/Debian you can do this e.g. by opening a terminal and type:
 
sudo apt-get update
 
Type in a terminal:
 
'''sudo apt-get install qgis'''
 
=== Installing the required plugins ===
 
Launch QGIS and enter the '''Plugins''' menu. There choose '''Manage and Install Plugins'''.
Select '''Not Installed''' and choose the following plugins for installation in the list:
 
* '''Semi-Automatic Classification Plugin'''
* '''QTiles'''
* OpenLayers Plugin
 
 
Once installed please right-click the toolbar area and enable the '''SCP Toolbar''' .
 
=== Opening the Sentinel data in QGIS using the SCP plugin ===
 
Launch QGIS and choose '''Layer -> Add Layer -> Add Raster Layer'''  (or the associated toolbar '''checkboard''' button).
Then navigate into the '''GRANULE''' directory. There choose the directory for the correct sub tile. The order of the tiles follows columns.
Inside the sub tile directory use the Ctrl button to '''choose multiple files at once''' - namly those which end on b2-b8:
 
'''b2.jp2, b3.jp2, b4.jp2, b5.jp2, b6.jp2, b7.jp2, b8.jp2'''
 
Press '''Open''' to open them and wait for them to be loaded which might take about 10 secs.
 
Now press the rainbow colored "Band set" button in the SCP toolbar. All files should appear in the upper section. Now apply the following steps:
* Press '''Select All'''
* Press Add '''rasters to set'''
* In the middle to the right press the '''Sort by Name''' button
* In the '''Quick wavelength settings''' press the drop down list combobox, scroll down and choose '''Sentinel-2 [bands 2, 3, 4, ...]'''
* Press '''Create virtual raster of band set'''
* Enter a '''.vrt''' name of your choice as a band set file name.
* Press '''Close'''
 
You should see a colorized version of your loaded tiles now. We suggest that you tweak the colors by right clicking on the vrt file entry on the right and press "Properties". In the Band rendering section we suggest that you make the following adjustments:
 
Red Band: Band 1 - 779 1581
Green Band: Band 2 - 634 1632
Blue Band: Band 3 - 600 2000
 
In the Resampling section choose
 
* Zoomed: in '''Bilinear''' out '''Average''' Oversampling 2,00
 
As a result you should now see a nice relatively realistic looking rendering of the sentinel data. You can zoom into it using the mouse scroll wheel down to street level (10m per pixel).
 
=== Tile Creation ===
 
With the new sentinel layer loaded and proper band sets in place we can now start to create our own tiles:
 
Press Plugins -> QTiles -> QTiles  then
 
* In the ''Output'' section choose '''Directory''' and browse to your directory of choice where the tiles should be stored
* Enter '''sentinel2''' as a Tileset name
* In the ''Extent'' section choose '''Layer Extent'''
* In the ''Zoom'' section choose '''14''' as a maximum zoom level
* In the ''Parameters'' section enable '''Make lines appear less jagged''' and '''Render tiles outside of layers extents (within combined extent)'''.
 
=== Deployment on the server ===
 
TBD

Latest revision as of 04:55, 26 October 2016

This page is now on the community wiki.