An error occurred while loading the file. Please try again.
-
Laura LINDEPERG authored6bde6a6d
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
# -*- coding: utf-8 -*-
"""
Created on Wed Mar 17 15:35:20 2021
@author: laura.lindeperg
"""
import matplotlib.pyplot as plt
import seaborn as sns
import numpy as np
from shapely.geometry import mapping
import rioxarray as rxr
import xarray as xr
import geopandas as gpd
import pandas as pd
# Carte des Zones Humides (raster)
ZHp_map_path = 'C:/Users/laura.lindeperg/Documents/DonneesLaura/ZonesHumidesPotentielles/mph-fr.tif'
ZHp_map = rxr.open_rasterio(ZHp_map_path).squeeze()
# Contour du bassin versant test (shapefile)
BV_W3315010_path = 'C:/Users/laura.lindeperg/Documents/INRAE_2021/CODE/fhysa/TestData/BH/W3315010_contour.shp'
BV_W3315010 = gpd.read_file(BV_W3315010_path)
# extraction du raster propre au contour du bassin versant et enregistrement du raster obtenu
ZHp_clipped = ZHp_map.rio.clip(BV_W3315010.geometry.apply(mapping), BV_W3315010.crs)
ZHp_clipped.rio.to_raster('C:/Users/laura.lindeperg/Documents/DonneesLaura/ZonesHumidesPotentielles/W3315010_ZHp.tif')
# Raster test
ZHp_test_path = 'C:/Users/laura.lindeperg/Documents/DonneesLaura/ZonesHumidesPotentielles/W3315010_ZHp.tif'
ZHp_test = rxr.open_rasterio(ZHp_test_path).squeeze()
# ZHp_test.values.flatten() -> fonction "unique" s'en occupe si ce n'est pas 1D
classes, counts = np.unique(ZHp_test, return_counts=True)
classesO, countsO = np.unique(ZHp_map, return_counts=True)
# Quelques tests
ZHp_test.plot.hist()
plt.show()
f, ax = plt.subplots(figsize=(10, 4))
ZHp_test.plot(ax=ax,
cmap='Greys')
ax.set(title="Final Clipped ZHp")
ax.set_axis_off()
plt.show()