Spatial subset
From AWF-Wiki
(Difference between revisions)
(→Interactive box selection) |
(→Vector polygon masking) |
||
Line 11: | Line 11: | ||
= Vector polygon masking = | = Vector polygon masking = | ||
* Load the vector file ''\geodata\vector\gadm36_DEU_4.shp''. These are the level 4 administrative boundaries of the country Germany downloaded from https://gadm.org/index.html . | * Load the vector file ''\geodata\vector\gadm36_DEU_4.shp''. These are the level 4 administrative boundaries of the country Germany downloaded from https://gadm.org/index.html . | ||
− | * Right click on Layer name in the Layers Panel{{mitem|text=Properties --> General}}. Change the '''Data source encoding''' from {{button|text=System}} to {{button|text=UTF-8}}. Click {{button|text=UTF-8}}. | + | * Right click on Layer name in the Layers Panel {{mitem|text=Properties --> General}}. Change the '''Data source encoding''' from {{button|text=System}} to {{button|text=UTF-8}}. Click {{button|text=UTF-8}}. |
* Right click on Layer name in the Layers Panel and open the {{button|text=Attribute table}}. | * Right click on Layer name in the Layers Panel and open the {{button|text=Attribute table}}. | ||
* Click [[File:Qgis_select_by_expression.png]] and select by the expression: {{typed|text="NAME_4" ='Göttingen'}}. {{button|text=Select}}. | * Click [[File:Qgis_select_by_expression.png]] and select by the expression: {{typed|text="NAME_4" ='Göttingen'}}. {{button|text=Select}}. | ||
[[File:Qgis_select_munic.png|400px]] | [[File:Qgis_select_munic.png|400px]] | ||
− | * Right click on Layer name in the Layers Panel {{button|text=Save as...}}. Fill the form as shown below. | + | * Right click on Layer name in the Layers Panel . {{mitem|text=Export --> Save Features as...}} {{button|text=Save as...}}. Fill the form as shown below. |
[[File:Qgis_vector_save_as.png|400px]] | [[File:Qgis_vector_save_as.png|400px]] | ||
− | Make sure that the CRS is the same like the image file (here EPSG:32632)! | + | Make sure that the CRS is the same like the image file (here WGS84/UTM Zone 32, EPSG:32632)! |
# Click {{button|text=Run}}. | # Click {{button|text=Run}}. | ||
# {{mitem|text=Raster --> Extraction ---> Clip Raster by Mask Layer...}} | # {{mitem|text=Raster --> Extraction ---> Clip Raster by Mask Layer...}} |
Revision as of 12:26, 29 October 2018
Interactive box selection
- Click the Open Data Source Manager and Add raster layer and select the file /geodata/raster/s2/Subset_S2A_MSIL2A_20170619T_MUL.tif.
- Raster --> Extraction ---> Clip Raster by Extent
- Input file (raster): Choose a loaded file from the drop down list or click Select to browse for a file.
- Output file: click Select to choose a new file name for the subset image.
- Define the clipping extent by clicking... and Select extent from canvas The clipping mode is now active. Hold a left click on the map canvas dragging a box of your region of interest.
- Adjust the extent coordinates in the text file rounding the coordinates (upper left and lower right).
- Click OK.
Vector polygon masking
- Load the vector file \geodata\vector\gadm36_DEU_4.shp. These are the level 4 administrative boundaries of the country Germany downloaded from https://gadm.org/index.html .
- Right click on Layer name in the Layers Panel Properties --> General. Change the Data source encoding from System to UTF-8. Click UTF-8.
- Right click on Layer name in the Layers Panel and open the Attribute table.
- Click and select by the expression: "NAME_4" ='Göttingen'. Select.
- Right click on Layer name in the Layers Panel . Export --> Save Features as... Save as.... Fill the form as shown below.
Make sure that the CRS is the same like the image file (here WGS84/UTM Zone 32, EPSG:32632)!
- Click Run.
- Raster --> Extraction ---> Clip Raster by Mask Layer...
- Input file (raster): Choose a loaded file name fron the drop down list or click Select to browse to a file (see above).
- Output file: click Select to choose a new file name for the subset image.
- Change the clipping mode from Extent to Mask. Choose the newly created polygon of the municipality as Mask Layer.
- Switch the checkbox on: Crop the extent of the target dataset to the extent of the cutline. All image data outside the polygon will be masked.
- Switch the No data value checkbox on and specify the value 0.
- Click OK.