<div dir="ltr"><div><div><div>Quick follow-up - Just tried importing one of these problematic layers into stand-alone SAGA (version 2.1.3, x64), which worked. (Only problem was the projection was flipped upside-down) - interesting as it was via the 'GDAL: Import NetCDF' tool.<br><br></div>And the gdal version I was using (from previous e-mail) is GDAL 1.11.2<br><br></div>Any thoughts on this are appreciated - thanks again!<br></div>mike<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jun 19, 2015 at 11:02 AM, Michael Treglia <span dir="ltr"><<a href="mailto:mtreglia@gmail.com" target="_blank">mtreglia@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div><div><div>Hi All,<br><br></div>I'm trying to work with the Global Soil Dataset for Earth System Modeling (available at: <a href="http://globalchange.bnu.edu.cn/research/soilw" target="_blank">http://globalchange.bnu.edu.cn/research/soilw</a>) and having mixed success with the netCDF files - wanted to get any thoughts from folks here as I'm pretty new to working with this format, and running into some issues.<br><br></div>I'm primarily working with the 8-layer files (which are split into two 4-band files each; each band represents a different horizon). I'm able to load netCDFs of some layers (in particular 'gypsum', 'CaCO3', and 'Total Carbon') without issue (using QGIS 2.8 - both import raster and via the netCDF viewer plugin). And running gdalinfo from the command line returns the metadata without issue. <br><br></div>However, I've also been trying to work with Sand, Silt, Clay, and Gravel Content (as % by weight), and not having any success in QGIS - the netCDF files aren't recognized by the netCDF viewer, and can't be pulled in via 'import raster'. Also,  gdalinfo returns the following error: <br>

<p>ERROR 4: `SAND1.nc' not recognised as a supported
file format.</p>

<br></div>I am able to read these data using R - most simply, using the 'raster' package (sample code below), and I can also get the info using ncdf and ncdf4 packages:<br></div>library(raster)<br></div><div>brk <- brick("[file <a href="http://name.nc" target="_blank">name.nc</a>]", varname="[variablename]")<br></div><div>#Can plot using<br></div><div>plot(brk) #plots all 4 layers from netcdf<br></div><div>plot(brk[[1]]) #Plots only first layer from netcdf<br></div><div><br></div>But, I'm trying to do this in QGIS or via gdal utilities for speed  - R is slow for this (may be ways to speed it up, but haven't toyed with it much).<br><br></div>I'd appreciate any ideas on what might be contributing to these issues - thanks in advance for any input!<br><br></div>Mike<br></div>
</blockquote></div><br></div>