2021年11月28日星期日
qgis:raster-surface-volume
https://ocw.un-ihe.org/mod/book/view.php?id=6788&chapterid=578
https://github.com/jvdkwast/PyQGIS_Hydro/blob/master/scripts/StageVolume.py
import numpy as np
levels = np.arange(12.9,20.0,0.5).tolist()
for i_level in levels:
r = processing.run("native:rastersurfacevolume", {'INPUT':'OUTPUT.tif','BAND':1,'LEVEL':i_level,'METHOD':1,'OUTPUT_HTML_FILE':'Skip Output'})
print('{}:{}:{}'.format(i_level,round(abs(r['VOLUME']/10**8),2),round(r['AREA']/10**6,2)))
How do I read the hypsometric curve?
Hypsographic Curve of Earth's Surface from ETOPO1
Elevation and Volume relationship curve in QGIS
2021年11月23日星期二
2021年11月22日星期一
2021年11月11日星期四
mikeio笔记
self.items中每个item只有一列数据,
输出的dataframe的表头,是 self.items中每个item的name。
2、Datase的长度,是self.items中item的个数。
3、 2021年11月16日,通过pip 安装的 mikeio1d,不支持res11_file.reaches属性
4、mikeio.dfs0中,DataFrame的to_dfs0方法,是通过import dfs0时,给pandas打补丁,实现的
# Monkey patching onto Pandas classes
pd.DataFrame.to_dfs0 = dataframe_to_dfs0
pd.Series.to_dfs0 = series_to_dfs0