Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It only takes a minute to sign up.
Sign up to join this community
My problem goes like this:
I have a gridded data in netCDF file. The cells have values in meters. For example: 0,012. I only need a small area of the grid and after clipping it in QGIS using a shapefile mask the values are multiplied with 10 000, result for the same cell being 120. I have to mention that the clipping is done on all the bands (365/366) without creating an alpha band. GDAL code looking fine.
Any ideas why it goes like this?
Have a look at the file's metadata. Most probably the values are stored as integer, with an internal scale factor of 10000. Reading the untouched file, QGIS applies the scale factor automatically.
The clipping is done by the external GDAL command gdal_translate, and that destroys the correct handling of the scale factor from the metadata.
As a workaround, use the raster calculator to divide all cell values by 10000. You might need to switch the cell value data type from integer to real using -ot as well. See http://www.gdal.org/gdal_translate.html for the full coammand line syntax.
