<div dir="ltr"><br><br><div class="gmail_quote">On Tue, 14 Apr 2015 at 05:46 Dominik Schneider <<a href="mailto:Dominik.Schneider@colorado.edu" target="_blank">Dominik.Schneider@colorado.edu</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><p style="margin:1.2em 0px!important">Hi - </p>
<p style="margin:1.2em 0px!important">I have some .mdl files from IDL with .hdr header files that I’d like to convert to netcdf.<br>The following code produces a netcdf file with a subdataset for each band in the mdl file.  Is there any way to define the bands as the time dimension, in this case 4416 hourly timesteps?</p><p style="margin:1.2em 0px!important">Or is there another tool that can convert this?<br></p>
<pre style="font-size:0.85em;font-family:Consolas,Inconsolata,Courier,monospace;font-size:1em;line-height:1.2em;margin:1.2em 0px"><code style="font-size:0.85em;font-family:Consolas,Inconsolata,Courier,monospace;margin:0px 0.15em;padding:0px 0.3em;white-space:pre-wrap;border:1px solid rgb(234,234,234);background-color:rgb(248,248,248);border-radius:3px;display:inline;white-space:pre-wrap;overflow:auto;border-radius:3px;border:1px solid rgb(204,204,204);padding:0.5em 0.7em;display:block!important;display:block;overflow-x:auto;padding:0.5em;color:rgb(51,51,51);background:none repeat scroll 0% 0% rgb(248,248,248)">gdal_translate -of NetCDF swe.mdl <a href="http://swe.nc" target="_blank">swe.nc</a>
</code></pre></div></div></blockquote><div><br></div><div>Can you list the metadata printed out by gdalinfo swe.mdl?  What driver is used? </div><div><br></div><div>Does that source file have time-metadata inside it each subdataset or do you need to assign it? </div><div><br></div><div>Generally, the subdatasets are not considered as a time series - they are more like multiple variables for the same observation (whereas a 3rd and higher dimension/s are often unrolled into a multi-attributed 2D layer and the time/z step is stored on each band -  it's a bit of a mix/clash of worlds). </div><div><br></div><div>Can you point to and example file somewhere?  I imagine you'll need a programmed solution, but it should be pretty easy with R or python or similar.  If R is of interest you can try this and take your question to the R community: </div><div><br></div><div>library(raster)</div><div>r <- stack("swe.mdl")</div><div>r <- setZ(r, [whatever the vector of dates should be])  ## pseudocode</div><div>writeRaster(r, "<a href="http://swe.nc">swe.nc</a>")</div><div><br></div><div>That all assumes quite a lot, including available NetCDF versions and packages, required structure in the output, interpretation of the data read in,  - so it's just included as an indication how it might be done. </div><div><br></div><div>Cheers, Mike. </div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Thanks<br>Dominik​</div></div>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a></blockquote></div></div>