<div dir="ltr">Stephen,<div><br></div><div>Take a look at these two discussions for starters on working with hycom:</div><div><br></div><div><a href="http://gis.stackexchange.com/questions/167155/how-can-i-make-geotiffs-of-individual-depth-layers-in-the-hycom-ocean-circulatio">http://gis.stackexchange.com/questions/167155/how-can-i-make-geotiffs-of-individual-depth-layers-in-the-hycom-ocean-circulatio</a><br></div><div><br></div><div><a href="http://gis.stackexchange.com/questions/170882/how-can-a-make-a-netcdf-with-subgroups-smaller">http://gis.stackexchange.com/questions/170882/how-can-a-make-a-netcdf-with-subgroups-smaller</a><br></div><div><br></div><div>You should start with the regular grids.  e.g.</div><div><br></div><div><a href="ftp://ftp.hycom.org/datasets/GLBu0.08/expt_91.1/data/hindcasts/2015/hycom_glb_911_2015093000_t000_uv3z.nc">ftp://ftp.hycom.org/datasets/GLBu0.08/expt_91.1/data/hindcasts/2015/hycom_glb_911_2015093000_t000_uv3z.nc</a><br></div><div><br></div><div>Make an 80 layer geotiff from the u/v and then use mapserver as Even suggested.</div><div><br></div><div>-kurt</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Jun 4, 2016 at 6:29 AM, Even Rouault <span dir="ltr"><<a href="mailto:even.rouault@spatialys.com" target="_blank">even.rouault@spatialys.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<span class="">><br>
> I want to convert some HYCOM data into GTiff format and I'm a little<br>
> lost on how to get started, whether I can do what I need just using a<br>
> vrt or if I'm going to have to write code to process the data.<br>
><br>
> To start with I have placed gdalinfo on one file here:<br>
> <a href="http://imaptools.com:8080/dl/hycom-gdalinfo.txt" rel="noreferrer" target="_blank">http://imaptools.com:8080/dl/hycom-gdalinfo.txt</a><br>
><br>
> The data comes from <a href="http://hycom.org/dataserver/glb-analysis" rel="noreferrer" target="_blank">http://hycom.org/dataserver/glb-analysis</a><br>
> and I will need to be using GLBa0.08 dataset.<br>
><br>
> I have two use cases:<br>
><br>
> 1) The simple one is to extract say Band 1 into a georeferenced GTiff<br>
> and apply color to using color ranges without loosing the NODATA. I<br>
> think I can do later using .vrt, correct?<br>
<br>
</span>gdaldem color-relief can output to a VRT, but given that this netCDF file is<br>
georeferenced through a geolocation array, you'll likely have first to run<br>
gdalwarp -geoloc to have something useful at the end.<br>
<span class=""><br>
><br>
> I'm not sure how to deal with the georeferencing, because the download<br>
> site says: "Native hycom .[ab] data converted to NetCDF on native<br>
> Mercator-curvilinear HYCOM horizontal grid" but that does not seem in<br>
> sync with the Metadata in the file.<br>
><br>
> 2) The more complicated case is that there are two files:<br>
>     a)  eastward_sea_water_velocity (aka: u, u-velocity)<br>
>     b)  northward_sea_water_velocity (aka: v, v-velocity)<br>
> that I would like to use to create a grid of vectors with arrowheads<br>
> using u and v to define the vector and where the color of the vector is<br>
> related to the magnitude of its length. Ultimately this data will get<br>
> rendered via mapserver so I'm wondering if I can do this via a GTiff or<br>
> if I will have to resort to creating a massive point shapefile with<br>
> attributes of angle and magnitude and use a symbol. It seems like it<br>
> would be best if I can keep the data as a GTiff and then sample the<br>
> points and render the vectors on the fly for when you zoom in/out.<br>
<br>
</span>See<br>
<a href="http://mapserver.org/input/vector/vector_field.html" rel="noreferrer" target="_blank">http://mapserver.org/input/vector/vector_field.html</a><br>
<br>
You don't need to create a vector file.<br>
<br>
You'll have to preprocess your file to get the u and v bands in the same<br>
raster.<br>
<span class="HOEnZb"><font color="#888888"><br>
Even<br>
<br>
--<br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" rel="noreferrer" target="_blank">http://www.spatialys.com</a><br>
</font></span><div class="HOEnZb"><div class="h5">_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">--<div><a href="http://schwehr.org" target="_blank">http://schwehr.org</a></div></div>
</div>