Python netcdf для hec hms dss ⇐ Python
Python netcdf для hec hms dss
С помощью xarray я создал netcdf осадков каждые 5 минут. С помощью Hec HMS 4.11 я могу преобразовать его в формат DSS, проблема в том, что при преобразовании импортируется только начальный столбец, но также и конечный столбец каждого интервала. Есть идеи, как при конвертации в DSS заполнить последний столбец?
Файлы (netcdf и DSS): https://drive.google.com/file/d/17KCpg2 ... sp=sharing
lista=[] для столбца в gdf.columns[:]: пытаться: fecha=pd.to_datetime(column, format='%Y-%m-%d %H:%M:%S') #si la columnsa no es una fecha sale del TRY geo_grid_cubic = make_geocube( гдф, измерения=[столбец], разрешение=(-25, 25), rasterize_function=partial(rasterize_points_griddata, метод="cubic"),) geo_grid_cubic = geo_grid_cubic.rio.clip(contorno.geometry.values, contorno.crs) #añado измерение с фешой da=geo_grid_cubic.assign_coords(время = феча) да = da.expand_dims(dim="время") da=da.rename({столбец: 'Осадки'}) da['Осадки'].attrs['имя'] = 'Осадки' da['Осадки'].attrs['long_name'] = 'Осадки' список.append(да) печать (столбец) кроме: пройти geotiff_da =xr.combine_by_coords(список) geotiff_da = geotiff_da .rename({'x': 'долгота','y': 'широта'}) geotiff_da.rio.write_crs("epsg:25830", inplace=True) Конвертировать netcdf в dss
С помощью xarray я создал netcdf осадков каждые 5 минут. С помощью Hec HMS 4.11 я могу преобразовать его в формат DSS, проблема в том, что при преобразовании импортируется только начальный столбец, но также и конечный столбец каждого интервала. Есть идеи, как при конвертации в DSS заполнить последний столбец?
Файлы (netcdf и DSS): https://drive.google.com/file/d/17KCpg2 ... sp=sharing
lista=[] для столбца в gdf.columns[:]: пытаться: fecha=pd.to_datetime(column, format='%Y-%m-%d %H:%M:%S') #si la columnsa no es una fecha sale del TRY geo_grid_cubic = make_geocube( гдф, измерения=[столбец], разрешение=(-25, 25), rasterize_function=partial(rasterize_points_griddata, метод="cubic"),) geo_grid_cubic = geo_grid_cubic.rio.clip(contorno.geometry.values, contorno.crs) #añado измерение с фешой da=geo_grid_cubic.assign_coords(время = феча) да = da.expand_dims(dim="время") da=da.rename({столбец: 'Осадки'}) da['Осадки'].attrs['имя'] = 'Осадки' da['Осадки'].attrs['long_name'] = 'Осадки' список.append(да) печать (столбец) кроме: пройти geotiff_da =xr.combine_by_coords(список) geotiff_da = geotiff_da .rename({'x': 'долгота','y': 'широта'}) geotiff_da.rio.write_crs("epsg:25830", inplace=True) Конвертировать netcdf в dss
-
- Похожие темы
- Ответы
- Просмотры
- Последнее сообщение
-
-
Извлечение направления потока/стрелки потока/вектора скорости из HEC-RAS
Anonymous » » в форуме Python - 0 Ответы
- 22 Просмотры
-
Последнее сообщение Anonymous
-