<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Hi Anna,<br>
<br>
Thank you for the reply.<br>
<br>
The output from running gdalinfo NETCDF:"<a moz-do-not-send="true"
href="http://global_30s_2009_153.nc/" target="_blank">global_30s_2009_153.nc</a>":lai
is below. I do not see anything in the metadata that suggests a
problem. Perhaps a more experienced eye can find the issue.<br>
<br>
Thanks!!<br>
<br>
-Damian<br>
<br>
Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute<br>
Driver: netCDF/Network Common Data Format<br>
Files: global_30s_2009_153.nc<br>
Size is 43200, 21600<br>
Coordinate System is `'<br>
Origin = (-180.000000000000000,90.000000000000000)<br>
Pixel Size = (0.008333333767951,-0.008333333767951)<br>
Metadata:<br>
lai#long_name=leaf area index<br>
lai#scale_factor=0.1<br>
lai#units=m2/m2<br>
lat#long_name=coordinate latitude<br>
lat#units=degrees_north<br>
lon#long_name=coordinate longitude<br>
lon#units=degrees_east<br>
NC_GLOBAL#citation=Yuan, H., Dai, Y., Xiao, Z., Ji, D., Shangguan,
W., 2011. Reprocessing the MODIS Leaf Area Index Products for Land
Surface and Climate Modelling. Remote Sensing of Environment,
115(5), 1171-1187. doi:10.1016/j.rse.2011.01.001<br>
NC_GLOBAL#creation_date=2011/09/27<br>
NC_GLOBAL#institution=Land-Atmosphere Interaction Research Group
at Beijing Normal University (<a class="moz-txt-link-freetext" href="http://globalchange.bnu.edu.cn">http://globalchange.bnu.edu.cn</a>)<br>
NC_GLOBAL#title=LAI data sets for land surface and climate
modelling<br>
NETCDF_DIM_EXTRA={time}<br>
NETCDF_DIM_time_DEF={1,4}<br>
NETCDF_DIM_time_VALUES=153<br>
time#long_name=day of year 2009<br>
time#units=day<br>
Corner Coordinates:<br>
Upper Left (-180.0000000, 90.0000000) <br>
Lower Left (-180.0000000, -90.0000094) <br>
Upper Right ( 180.0000188, 90.0000000) <br>
Lower Right ( 180.0000188, -90.0000094) <br>
Center ( 0.0000094, -0.0000047) <br>
Band 1 Block=43200x1 Type=Int16, ColorInterp=Undefined<br>
NoData Value=-32767<br>
Offset: 0, Scale:0.100000001490116<br>
Metadata:<br>
long_name=leaf area index<br>
NETCDF_DIM_time=153<br>
NETCDF_VARNAME=lai<br>
scale_factor=0.1<br>
units=m2/m2<br>
<br>
<br>
<div class="moz-cite-prefix">On 03/20/2014 11:12 PM, Anna Petráová
wrote:<br>
</div>
<blockquote
cite="mid:CAE0EDEqLiornYKdutVWSDdmtmhD56M3UxEoHpdZkSiM2GXeyWw@mail.gmail.com"
type="cite">
<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
moz-do-not-send="true"
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 moz-do-not-send="true"
href="http://globalchange.bnu.edu.cn/research/lai/"
target="_blank">http://globalchange.bnu.edu.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 moz-do-not-send="true"
href="http://global_30s_2009_153.nc" target="_blank">global_30s_2009_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
moz-do-not-send="true" href="http://r.info"
target="_blank">r.info</a> reports min of 0 and max of
70. <a moz-do-not-send="true" 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 moz-do-not-send="true"
href="http://global_30s_2009_153.nc/" target="_blank">global_30s_2009_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>
_______________________________________________<br>
grass-user mailing list<br>
<a moz-do-not-send="true"
href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a moz-do-not-send="true"
href="http://lists.osgeo.org/mailman/listinfo/grass-user"
target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
</blockquote>
</div>
<br>
</div>
</div>
</blockquote>
<br>
</body>
</html>