Split stack
From AWF-Wiki
(Difference between revisions)
(→Using GDAL Translate) |
(→Using OTB SplitImage) |
||
| (13 intermediate revisions by one user not shown) | |||
| Line 1: | Line 1: | ||
=Using OTB SplitImage= | =Using OTB SplitImage= | ||
| − | # Click the {{button|text=Open Data Source Manager}} [[Image:QGIS_3.0_datasource.png|30px]] and then {{button|text=Add raster layer}} [[Image:QGIS_2.0_addrast.png|20px]] and select | + | # Click the {{button|text=Open Data Source Manager}} [[Image:QGIS_3.0_datasource.png|30px]] and then {{button|text=Add raster layer}} [[Image:QGIS_2.0_addrast.png|20px]] and select a multiband raster file. |
# In processing toolbox, type {{typed|text=split image}} into the search field to find the {{button|text=OTB --> Image Manipulation --> SplitImage}} tool and open it. | # In processing toolbox, type {{typed|text=split image}} into the search field to find the {{button|text=OTB --> Image Manipulation --> SplitImage}} tool and open it. | ||
| − | [[File:qgis-otb-splitimage.png| | + | [[File:qgis-otb-splitimage.png|500px]] |
| − | + | Choose the same ''Output pixel type'' like the input pixel type: here it is {{typed|text=uint16}}. | |
| − | + | ||
| − | + | '''Important''': choose a physical output file: {{button|text=Save to file}}. Specify a file name with the extension '''.tif''' | |
| + | Uncheck "Open output file after running algorithm". Open the output files manually in QGIS. | ||
=Using GDAL Translate= | =Using GDAL Translate= | ||
| − | + | * In the processing toolbar, type {{typed|text=Translate}} into the search field to find the {{button|text=GDAL\OGR --> Translate (Convert format)}} tool and open it. | |
| − | + | * Change the output data format to unsigned 16bit integer): {{typed|text=-ot UInt16}} | |
| − | + | * Click the button {{button|text=...}} of '''Converted''' to select path and name of output file. | |
| + | |||
[[File:Qgis_split_stack.png|600px]] | [[File:Qgis_split_stack.png|600px]] | ||
| − | |||
| − | |||
| − | |||
| − | |||
| − | ** Repeat the steps for other bands: change | + | * Open the the OSGEO4W Shell. {{mitem|text=Start --> All Programs --> OSGeo4W --> OSGeo4W Shell}} |
| + | * Copy the content of the ''GDL/OGR console call'' and paste it into the OSGeo4W Shell. | ||
| + | * Edit the command in OSGeo4W Shell by adding the option {{typed|text=-b 1}} to output only the first band. | ||
| + | * Press {{button|text=Enter}} to execute the command in OSGeo4W Shell. | ||
| + | |||
| + | * Repeat the steps for other bands: change the output file name and edit the -b option accordingly. | ||
| − | + | * Create a multiband file with a spectral subset of 4 bands. Change the output file name and edit the command in OSGeo4W Shell additional creation options: | |
| − | ** '''Converted''' output file:{{typed|text=Subset_S2A_MSIL2A_20170619T_B02- | + | **'''Converted''' output file:{{typed|text=Subset_S2A_MSIL2A_20170619T_B02- |
-B03-B04-B08.tif}} | -B03-B04-B08.tif}} | ||
** '''Additional creation parameters:''' {{typed|text=-b 2 -b 3 -b 4 -b 8}} | ** '''Additional creation parameters:''' {{typed|text=-b 2 -b 3 -b 4 -b 8}} | ||
| − | ** Press {{button|Enter}} to execute the command OSGeo4W Shell. | + | ** Press {{button|text=Enter}} to execute the command in OSGeo4W Shell. |
[[Category:Working with Raster Data]] | [[Category:Working with Raster Data]] | ||
Latest revision as of 10:56, 23 May 2022
[edit] Using OTB SplitImage
- Click the Open Data Source Manager
and then Add raster layer 
and select a multiband raster file.
- In processing toolbox, type split image into the search field to find the OTB --> Image Manipulation --> SplitImage tool and open it.
Choose the same Output pixel type like the input pixel type: here it is uint16.
Important: choose a physical output file: Save to file. Specify a file name with the extension .tif Uncheck "Open output file after running algorithm". Open the output files manually in QGIS.
[edit] Using GDAL Translate
- In the processing toolbar, type Translate into the search field to find the GDAL\OGR --> Translate (Convert format) tool and open it.
- Change the output data format to unsigned 16bit integer): -ot UInt16
- Click the button ... of Converted to select path and name of output file.
- Open the the OSGEO4W Shell. Start --> All Programs --> OSGeo4W --> OSGeo4W Shell
- Copy the content of the GDL/OGR console call and paste it into the OSGeo4W Shell.
- Edit the command in OSGeo4W Shell by adding the option -b 1 to output only the first band.
- Press Enter to execute the command in OSGeo4W Shell.
- Repeat the steps for other bands: change the output file name and edit the -b option accordingly.
- Create a multiband file with a spectral subset of 4 bands. Change the output file name and edit the command in OSGeo4W Shell additional creation options:
- Converted output file:Subset_S2A_MSIL2A_20170619T_B02-
-B03-B04-B08.tif
- Additional creation parameters: -b 2 -b 3 -b 4 -b 8
- Press Enter to execute the command in OSGeo4W Shell.
