<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="PT" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US">Hi,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">There are a couple of issues that show when I run this test case.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">First is that ERROR 3<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">$ gdal_grid -ot Float32 -txe 0.000000 10.000000 -tye 0.000000 10.000000 -outsize 201 201 -a nearest lixo.vrt junk.nc<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Warning 1: Several drivers matching nc extension. Using NETCDF<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">0...10...20...30...40...50...60...70...80...90...100 - done.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">ERROR 3: An attempt to read an undefined ID from virtual variable collection was made<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Repeating the command (that is, overwriting the grid) shows more warnings, namely the “dimension #1 (lon) is not a Longitude/X dimension"<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">$ gdal_grid -ot Float32 -txe 0.000000 10.000000 -tye 0.000000 10.000000 -outsize 201 201 -a nearest lixo.vrt junk.nc<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Warning 1: Several drivers matching nc extension. Using NETCDF<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Warning 1: The dataset has several variables that could be identified as vector fields, but not all share the same primary dimension. Consequently they will be ignored.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Warning 1: dimension #1 (lon) is not a Longitude/X dimension.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Warning 1: dimension #0 (lat) is not a Latitude/Y dimension.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">0...10...20...30...40...50...60...70...80...90...ERROR 3: An attempt to read an undefined ID from virtual variable collection was made<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Which also show in gdalinfo<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">$ gdalinfo junk.nc<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Warning 1: dimension #1 (lon) is not a Longitude/X dimension.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Warning 1: dimension #0 (lat) is not a Latitude/Y dimension.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Driver: netCDF/Network Common Data Format<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">…<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">I see that the x,y coordinates are always named “lon,lat” but data is not in geographics those names are not true. In GMT we call them “x,y” what data represent cartesian units or “lon,lat” when they are geographic. As
a consequence of this, gdalinfo report those “not a Longitude/X dimension" warnings for the former case when the nc grid has not an “spatial_ref” attribute.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">The data used in this example is attached<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Joaquim<o:p></o:p></span></p>
</div>
</body>
</html>