<div dir="ltr">Hi Damian,<br><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Mar 19, 2014 at 1:47 PM, Damian Maddalena <span dir="ltr"><<a href="mailto:damianmaddalena@gmail.com" target="_blank">damianmaddalena@gmail.com</a>></span> wrote:<br>



<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">I am having some difficulty with the import of a netCDF dataset into GRASS 7. The dataset I am attempting to import is a 30 second global LAI dataset that can be found here: <a href="http://globalchange.bnu.edu.cn/research/lai/" target="_blank">http://globalchange.bnu.edu.<u></u>cn/research/lai/</a><br>




<br>
I have successfully imported other netCDF files into this build (worldCLIM and the soils data from the same source as this LAI data, for example), but this particular data is giving me issues. I have tried several files from several of the years available of the global LAI product. I get the same results.<br>




<br>
I am running the following commands in a global lat/lon location:<br>
<br>
g.region n=90 s=-90 e=-180 w=180 res=0.008333333767951  (I have also tried res=00:00:30.)<br></blockquote><div><br></div><div>r.in.gdal does not respect region</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">



r.in.gdal input=NETCDF:"<a href="http://global_30s_2009_153.nc" target="_blank">global_30s_2009_<u></u>153.nc</a>":lai output=testNETCDF -o<br>
<br>
The file imports with some warnings.<br>
<br>
Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute<br>
WARNING: Over-riding projection check<br>
WARNING: Fixing subtle input data rounding error of south boundary<br>
         (-9.38773e-06>4.62963e-08)<br>
<br>
It appears to be null when plotted, though <a href="http://r.info" target="_blank">r.info</a> reports min of 0 and max of 70.  <a href="http://r.info" target="_blank">r.info</a> hints at an issue with the e-w resolution:<br>




<br>
 |            N:        90N    S:        90S   Res: 0:00:30<br>
 |            E: 179:59:59.932408W    W:       180W   Res: 0:00:00.000002<br>
<br>
r.stats seems to indicate that all data is null?<br>
<br>
r.stats testNETCDF -cp<br>
100%<br>
* 900072000 100.00%<br>
<br>
d.histogram plots nothing.<br>
<br>
Any hints would be greatly appreciated.<br></blockquote><div><br></div><div>I wanted to try it myself, but it's too big to download it now. Could you post output of:</div><div><br></div><div>gdalinfo NETCDF:"<a href="http://global_30s_2009_153.nc/" target="_blank">global_30s_2009_<u></u>153.nc</a>":lai</div>

<div><br></div><div>I am not an expert on netcdf but others might have some idea based on the gdalinfo output.</div><div><br></div><div>Regards,</div><div><br></div><div>Anna</div><div><br></div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">

<br>
Thanks!!<br>
<br>
-Damian<br>
<br>
<br>
<br>
<br>
______________________________<u></u>_________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/<u></u>mailman/listinfo/grass-user</a><br>
</blockquote></div><br></div></div>