def main(): main_dir = Path(r'P:\Synchronize\IWS\Testings\fourtrans_practice\multisite_phs_spec_corr\5min\v7') os.chdir(main_dir) path_to_shp = r'P:\Synchronize\IWS\QGIS_Neckar\raster\taudem_out_spate_rockenau\watersheds_all.shp' label_field = r'DN' var = 'orig' suff = '__RR1D_RTsum' path_to_ras = f'{var}/{var}{suff}.nc' nc_x_crds_label = 'X' nc_y_crds_label = 'Y' nc_variable_labels = ['OK'] nc_time_label = 'time' overwrite_flag = True path_to_output = Path(f'{var}/{var}{suff}.h5') if overwrite_flag and path_to_output.exists(): os.remove(path_to_output) Ext = Extract(True) res = Ext.extract_from_netCDF( path_to_shp, label_field, path_to_ras, path_to_output, nc_x_crds_label, nc_y_crds_label, nc_variable_labels, nc_time_label) print('res:', res) return
def main(): main_dir = Path(r'P:\Synchronize\IWS\Testings\fourtrans_practice\multisite_phs_spec_corr\temperature_kriging') os.chdir(main_dir) path_to_shp = r'P:\Synchronize\IWS\QGIS_Neckar\raster\taudem_out_spate_rockenau\watersheds.shp' label_field = r'DN' path_to_ras = r'kriging_1km_mag.nc' input_ras_type = 'nc' # path_to_ras = r'P:\Synchronize\IWS\Colleagues_Students\Mischa\lulc_geohyd_ratio_rasters\lower_de_gauss_z3_1km_hydrogeol_einheit_nr_hydmod_lulc_ratios.tif' # input_ras_type = 'gtiff' nc_x_crds_label = 'X' nc_y_crds_label = 'Y' nc_variable_labels = ['OK'] nc_time_label = 'freq' src_epsg = None dst_epsg = None # src_epsg = 4326 # dst_epsg = 31467 # main_dir = Path(r'P:\Downloads\spinterp_2d_nc_crds_test') # os.chdir(main_dir) # # path_to_shp = r'01Small.shp' # # label_field = r'Id' # # path_to_ras = r'pr_SAM-44_ICHEC-EC-EARTH_historical_r12i1p1_SMHI-RCA4_v3_day_19810101-19851231.nc' # input_ras_type = 'nc' # # nc_x_crds_label = 'lon' # nc_y_crds_label = 'lat' # nc_variable_labels = ['pr'] # nc_time_label = 'time' path_to_output = Path(r'neckar_6cats_temp_1991_1991.h5') # path_to_output = 'lower_de_gauss_z3_1km_hydrogeol_einheit_nr_hydmod_lulc_ratios.h5' Ext = Extract(True) res = None if input_ras_type == 'gtiff': res = Ext.extract_from_geotiff( path_to_shp, label_field, path_to_ras, path_to_output, src_epsg, dst_epsg) elif input_ras_type == 'nc': res = Ext.extract_from_netCDF( path_to_shp, label_field, path_to_ras, path_to_output, nc_x_crds_label, nc_y_crds_label, nc_variable_labels, nc_time_label, src_epsg, dst_epsg) else: raise NotImplementedError print('\n') print('res:', res) return
def main(): main_dir = Path(r'P:\Downloads\spinterp_nc_gtiff_test') os.chdir(main_dir) path_to_shp = r'watersheds.shp' label_field = r'DN' path_to_ras = r'tem_spinterp.nc' input_ras_type = 'nc' # path_to_ras = r'lower_de_gauss_z3_2km_atkis_19_extended_hydmod_lulc_ratios.tif' # input_ras_type = 'gtiff' nc_x_crds_label = 'X' nc_y_crds_label = 'Y' nc_variable_labels = ['OK', 'SK'] nc_time_label = 'time' # main_dir = Path(r'P:\Downloads\spinterp_2d_nc_crds_test') # os.chdir(main_dir) # # path_to_shp = r'01Small.shp' # # label_field = r'Id' # # path_to_ras = r'pr_SAM-44_ICHEC-EC-EARTH_historical_r12i1p1_SMHI-RCA4_v3_day_19810101-19851231.nc' # input_ras_type = 'nc' # # nc_x_crds_label = 'lon' # nc_y_crds_label = 'lat' # nc_variable_labels = ['pr'] # nc_time_label = 'time' path_to_output = 'test.h5' Ext = Extract(True) res = None if input_ras_type == 'gtiff': res = Ext.extract_from_geotiff( path_to_shp, label_field, path_to_ras, path_to_output) elif input_ras_type == 'nc': res = Ext.extract_from_netCDF( path_to_shp, label_field, path_to_ras, path_to_output, nc_x_crds_label, nc_y_crds_label, nc_variable_labels, nc_time_label) else: raise NotImplementedError print(res) return