Return to Question

I need to make a script that processes a NetCDF file that contains 3 days of hourly forecast data from the norwegian meteorological office.

1. The NetCDF file contains various data I need (Precipitation,Tmperature,Wind etc).

2. The NetCDF file is in a lambert projection while I will need to project it into WGS84 UTM 32N.

3. Also I will need to resample from 2.5km (forecast inputs) to 1km(output) grid cells.

4. I need to save it into the IDRISI format .rst

PROBLEM! The original NetCDF is HUGE, covering the whole of scandinavia + neighbouring countries. Thus I will need a system that processes quick.

I managed to do this already with ArcPy, but the process was too slow since for every hour timestep I needed to extract one by one the huge rasters, and only then could I clip them down.

Maybe in GDAL (in Pythn) there is a way to first clip at once the whole netcdf and then continue the processing with a smaller netcdf?

I need to make a script that processes a NetCDF file that contains 3 days of hourly forecast data from the norwegian meteorological office.

1. The NetCDF file contains various data I need (Precipitation,Tmperature,Wind etc).

2. The NetCDF file is in a lambert projection while I will need to project it into WGS84 UTM 32N.

3. Also I will need to resample from 2.5km (forecast inputs) to 1km(output) grid cells.

4. I need to save it into the IDRISI format .rst

PROBLEM! The original NetCDF is HUGE, covering the whole of scandinavia + neighbouring countries. Thus I will need a system that processes quick.

I managed to do this already with ArcPy, but the process was too slow since for every hour timestep I needed to extract one by one the huge rasters, and only then could I clip them down.

Maybe in GDAL (in Python) there is a way to first clip at once the whole netcdf and then continue the processing with a smaller netcdf?

I need to make a script that processes a NetCDF file that contains 3 days of hourly forecast data from the norwegian meteorological office.

1)The NetCDF file contains various data i need (Precipitation,Tmperature,Wind etc).

3)The NetCDF file is in a lambert projection while i will need to project it into WGS84 UTM 32N.

4. Also i will need to resample from 2.5km (forecast inputs) to 1km(output) grid cells.

5. I need to save it into the IDRISI format .rst

PROBLEM! The original NetCDF is HUGE, covering the whole of scandinavia + neighbouring countries. Thus I will need a system that processes quick.

I managed to do this already with ArcPy, but the process was too slow since for every hour timestep i needed to extract one by one the huge rasters, and only then could i clip them down.

Maybe in GDAL (in Python) there is a way to first clip at once the whole netcdf and then continue the processing with a smaller netcdf?

Thanks in advance, Niccolo

I need to make a script that processes a NetCDF file that contains 3 days of hourly forecast data from the norwegian meteorological office.

1. The NetCDF file contains various data I need (Precipitation,Tmperature,Wind etc).

2. The NetCDF file is in a lambert projection while I will need to project it into WGS84 UTM 32N.

3. Also I will need to resample from 2.5km (forecast inputs) to 1km(output) grid cells.

4. I need to save it into the IDRISI format .rst

PROBLEM! The original NetCDF is HUGE, covering the whole of scandinavia + neighbouring countries. Thus I will need a system that processes quick.

I managed to do this already with ArcPy, but the process was too slow since for every hour timestep I needed to extract one by one the huge rasters, and only then could I clip them down.

Maybe in GDAL (in Pythn) there is a way to first clip at once the whole netcdf and then continue the processing with a smaller netcdf?

I need to make a script that processes a NetCDF file that contains 3 days of hourly forecast data from the norwegian meteorological office.

1)The NetCDF file contains various data i need (Precipitation,Tmperature,Wind etc).

2)The Netcdf file covers a huge area (all scandinavia + baltic countries), while I will need just a rectangle that covers a smaller region.

3)The NetCDF file is in a lambert projection while i will need to project it into WGS84 UTM 32N.

4. Also i will need to resample from 2.5km (forecast inputs) to 1km(output) grid cells.

So I would like to ask you how to extract with GDAL in Python in a fast way the rasters from the netcdf already clipping the initial domain to the desired subdomain. I will need to batch prcess this so that i can get for example all of the precipitation raster, hours by hour, out of the three days forecast.

I managed to do this already with ArcPy, but the process was too slow since for every hour timestep i needed to extract one by one the huge rasters, and only then could i clip them down.

Maybe in GDAL there is a way to first clip at once the whole netcdf and then continue the processing with a smaller netcdf?

I need to make a script that processes a NetCDF file that contains 3 days of hourly forecast data from the norwegian meteorological office.

1)The NetCDF file contains various data i need (Precipitation,Tmperature,Wind etc).

3)The NetCDF file is in a lambert projection while i will need to project it into WGS84 UTM 32N.

4. Also i will need to resample from 2.5km (forecast inputs) to 1km(output) grid cells.

5. I need to save it into the IDRISI format .rst

PROBLEM! The original NetCDF is HUGE, covering the whole of scandinavia + neighbouring countries. Thus I will need a system that processes quick.

I managed to do this already with ArcPy, but the process was too slow since for every hour timestep i needed to extract one by one the huge rasters, and only then could i clip them down.

Maybe in GDAL (in Python) there is a way to first clip at once the whole netcdf and then continue the processing with a smaller netcdf?

Thanks in advance, Niccolo