<html>
<head>
<meta http-equiv="Content-Type" content="text/html;
charset=windows-1252">
</head>
<body>
<p>See <a class="moz-txt-link-freetext" href="https://github.com/OSGeo/gdal/issues/4581">https://github.com/OSGeo/gdal/issues/4581</a><br>
</p>
<div class="moz-cite-prefix">Le 20/10/2021 à 02:01, Joaquim Manuel
Freire Luís a écrit :<br>
</div>
<blockquote type="cite"
cite="mid:AM6PR04MB3959F174B69DBFED38FECC83A6BE9@AM6PR04MB3959.eurprd04.prod.outlook.com">
<meta http-equiv="Content-Type" content="text/html;
charset=windows-1252">
<meta name="Generator" content="Microsoft Word 15 (filtered
medium)">
<style>@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;}p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
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;}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]-->
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US">That the geotransform is
not detected<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 xxx.grd<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Warning 1: dimension #1
(x) is not a Longitude___/X dimension.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Warning 1: dimension #0
(y) 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">Corner Coordinates:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Upper Left ( 0.0,
0.0)<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Lower Left ( 0.0,
59.0)<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Upper Right ( 37.0,
0.0)<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Lower Right ( 37.0,
59.0)<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Center ( 18.5,
29.5)<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 can’t get why it’s
complaining on “dimension #1 (x) is not a Longitude/X
dimension.”<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Tried to check the full
docs (netcdfdataset.cpp) but not even with a debugger I can
understand the why.<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">Due to it the x,y
coordinates are not fished from file, either from Julia or
from GMT<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">NOTE: this works but
does not use GDAL<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">grdinfo xxx.grd<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">xxx.grd: Title: z<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">xxx.grd: Command:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">xxx.grd: Remark:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">xxx.grd: Gridline node
registration used [Cartesian grid]<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">xxx.grd: Grid file
format: nf = GMT netCDF format (32-bit float), CF-1.7<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">xxx.grd: x_min: 580155
x_max: 581235 x_inc: 30 name: x n_columns: 37<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="FR">xxx.grd: y_min: 4416795
y_max: 4418535 y_inc: 30 name: y n_rows: 59<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">xxx.grd: v_min:
0.205308780074 v_max: 0.56782335043 name: z<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">This DOES NOT WORK (it
uses GDAL)<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">grdinfo xxx.grd=gd<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Warning 1: dimension #1
(x) is not a Longitude___/X dimension.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Warning 1: dimension #0
(y) is not a Latitude/Y dimension.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">xxx.grd: Title: Grid
imported via GDAL<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">xxx.grd: Command:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">xxx.grd: Remark:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">xxx.grd: Gridline node
registration used [Cartesian grid]<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">xxx.grd: Grid file
format: gd = Import/export through GDAL<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">xxx.grd: x_min: 0.5
x_max: 36.5 x_inc: 1 name: x n_columns: 37<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="FR">xxx.grd: y_min: 0.5 y_max:
58.5 y_inc: 1 name: y n_rows: 59<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">xxx.grd: v_min:
0.205308780074 v_max: 0.56782335043 name: z<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">xxx.grd: scale_factor: 1
add_offset: 0<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">+proj=utm +zone=10
+ellps=WGS84 +units=m +no_defs<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">What’s so wrong with
this nc file? I have others that look similar in terms of
attributes and things work fine, but for some other cases I
get the terrible<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">Warning 1: dimension #1
(x) is not a Longitude/X dimension.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Warning 1: dimension #0
(y) is not a Latitude/Y 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">and I’m screwed.<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>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<pre class="moz-quote-pre" wrap="">_______________________________________________
gdal-dev mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/gdal-dev">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a>
</pre>
</blockquote>
<pre class="moz-signature" cols="72">--
<a class="moz-txt-link-freetext" href="http://www.spatialys.com">http://www.spatialys.com</a>
My software is free, but my time generally not.</pre>
</body>
</html>