Data & Methods¶
The orthophoto and real estate data for this case study can be downloaded here:
Read netCDF files¶
All data is provided in netcdf-Format and can thus easily be read into python via:
import xarray as xr
data = xr.open_dataset("data/oberschelden.nc")
The orthophoto and real estate data is structured in an xr.DataSet with 4 variables:
data
<xarray.Dataset>
Dimensions: (x: 13224, y: 11108)
Coordinates:
* y (y) float64 5.636e+06 5.636e+06 ... 5.635e+06 5.635e+06
* x (x) float64 4.253e+05 4.253e+05 ... 4.266e+05 4.266e+05
Data variables:
buildings_mask (y, x) uint8 ...
red (y, x) uint8 ...
green (y, x) uint8 ...
blue (y, x) uint8 ...xarray.Dataset
- x: 13224
- y: 11108
- y(y)float645.636e+06 5.636e+06 ... 5.635e+06
array([5636252.35, 5636252.25, 5636252.15, ..., 5635141.85, 5635141.75, 5635141.65]) - x(x)float644.253e+05 4.253e+05 ... 4.266e+05
array([425313.85, 425313.95, 425314.05, ..., 426635.95, 426636.05, 426636.15])
- buildings_mask(y, x)uint8...
- transform :
- [ 1.0000000e-01 0.0000000e+00 4.2531380e+05 0.0000000e+00 -1.0000000e-01 5.6362524e+06]
- crs :
- +init=epsg:25832
- res :
- [0.1 0.1]
- is_tiled :
- 0
- nodatavals :
- [0. 0. 0. 0.]
- scales :
- [1. 1. 1. 1.]
- offsets :
- [0. 0. 0. 0.]
- AREA_OR_POINT :
- Area
[146892192 values with dtype=uint8]
- red(y, x)uint8...
- transform :
- [ 1.0000000e-01 0.0000000e+00 4.2531380e+05 0.0000000e+00 -1.0000000e-01 5.6362524e+06]
- crs :
- +init=epsg:25832
- res :
- [0.1 0.1]
- is_tiled :
- 0
- nodatavals :
- [nan nan nan nan]
- scales :
- [1. 1. 1. 1.]
- offsets :
- [0. 0. 0. 0.]
- AREA_OR_POINT :
- Area
[146892192 values with dtype=uint8]
- green(y, x)uint8...
- transform :
- [ 1.0000000e-01 0.0000000e+00 4.2531380e+05 0.0000000e+00 -1.0000000e-01 5.6362524e+06]
- crs :
- +init=epsg:25832
- res :
- [0.1 0.1]
- is_tiled :
- 0
- nodatavals :
- [nan nan nan nan]
- scales :
- [1. 1. 1. 1.]
- offsets :
- [0. 0. 0. 0.]
- AREA_OR_POINT :
- Area
[146892192 values with dtype=uint8]
- blue(y, x)uint8...
- transform :
- [ 1.0000000e-01 0.0000000e+00 4.2531380e+05 0.0000000e+00 -1.0000000e-01 5.6362524e+06]
- crs :
- +init=epsg:25832
- res :
- [0.1 0.1]
- is_tiled :
- 0
- nodatavals :
- [nan nan nan nan]
- scales :
- [1. 1. 1. 1.]
- offsets :
- [0. 0. 0. 0.]
- AREA_OR_POINT :
- Area
[146892192 values with dtype=uint8]
Task¶
Download all data sets and make yourselves familiar with the data