Puedo crear un archivo NetCDF (utilizando la biblioteca netCDF4 de Python) y cargarlo en QGIS. Sin embargo, como el tamaño de estos archivos es a veces bastante grande, intento comprimir mis datos. Esto funciona muy bien usando
createVariable(h, 'float32', ("y", "x"), zlib=True,complevel=3,chunksizes=(1,len(x)),fill_value=np.nan)
Sin embargo, al jugar un poco se ha visto que al fragmentar ambas dimensiones (por ejemplo, chunksizes=(100,100) o chunksizes=(len(y),len(x))) se reduce el tamaño del archivo resultante en más de un 50%. Desgraciadamente, en este caso ya no soy capaz de cargar esta cuadrícula correctamente en QGIS. El único chunksize que funciona es chunking sólo la variable más a la derecha (1,len(x)). Supongo que esto tiene que ver con el GDAL. Sobre QGIS dice:
QGIS versión 2.14.12-Essen
Compilado con GDAL/OGR 1.10.1
Ejecutando contra GDAL/OGR 1.11.2
Estoy trabajando en Ubuntu14.04 e instalé QGIS a través del PPA debian-ltr. ¿Hay alguna manera de hacer que otros chunksizes funcionen también? ¿Ayudaría una actualización a QGIS 2.18 o es un problema de GDAL?
2 votos
Parece que la compilación de UbuntuGIS gdal incluye soporte para netcdf4/hdf5, que es necesario para los archivos fragmentados. Intente añadir la ppa ubuntugis-unstable y cambie su repo de qgis de debian-ltr a la versión con dependencias de ubuntugis qgis.org/ubuntugis-ltr