<div dir="ltr">Hi Stephen,<div><br></div><div>I sounds like you should be asking on <a href="https://groups.google.com/a/hycom.org/forum/#!forum/forum">https://groups.google.com/a/hycom.org/forum/#!forum/forum</a> rather than gdal-dev.</div><div><br></div><div>-kurt</div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Jun 5, 2016 at 2:10 PM, Stephen Woodbridge <span dir="ltr"><<a href="mailto:woodbri@swoodbridge.com" target="_blank">woodbri@swoodbridge.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Even, Kurt,<br>
<br>
Thank you for the suggestions.<br>
<br>
I'm reading through stuff and searching for solutions. I need access to the 2d layer which is only available via the GLBa0.08 dataset and not in the GLBu0.08 dataset which is easier to deal with. I'll start with the GLBu0.08 dataset for the other data following Kurt suggestions.<br>
<br>
This page describes the problem of georeferencing the GLBa0.08 dataset but the solution is based on windows python script and using ArcGIS, but the idea might be adapted to writing some C/C++ or script to handle creating those artifacts. I want to script the download and processing of the data into geotiff's running on Linux.<br>
<br>
<a href="http://code.nicholas.duke.edu/projects/mget/wiki/HYCOM" rel="noreferrer" target="_blank">http://code.nicholas.duke.edu/projects/mget/wiki/HYCOM</a><br>
<br>
One problem is that the HYCOM site does not have *.grid.[ab] files that I can find. They are supposed to be in the topo directory:<br>
<br>
Index of <a href="ftp://ftp.hycom.org/datasets/GLBa0.08/expt_91.2/topo/" rel="noreferrer" target="_blank">ftp://ftp.hycom.org/datasets/GLBa0.08/expt_91.2/topo/</a><br>
<br>
Name Size Last Modified<br>
depth_GLBa0.08_09.a 57984 KB 8/24/2013 12:00:00 AM<br>
depth_GLBa0.08_09.b 1 KB 8/24/2013 12:00:00 AM<br>
<a href="http://depth_GLBa0.08_09.nc" rel="noreferrer" target="_blank">depth_GLBa0.08_09.nc</a> 173950 KB 11/12/2013 12:00:00 AM<br>
<br>
and it is not obvious how to read and incorporate these files as I think they only work with the HYCOM related executables.<br>
<br>
I did notice the the lat/lon info is already in the .nc for the "u" dataset which can be seen with ncdump. In fact it is also there for the "a" dataset also, but the values 1-3298 and 1-4500 which is the size of the file in pixels, so not helpful.<br>
<br>
I've looked on the HYCOM forum but that has not been very helpful yet.<br>
<br>
Thanks,<br>
-Steve<span class=""><br>
<br>
<br>
On 6/5/2016 2:04 PM, Kurt Schwehr wrote:<br>
</span><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">
Stephen,<br>
<br>
Take a look at these two discussions for starters on working with<br>
hycom:<br>
<br>
<a href="http://gis.stackexchange.com/questions/167155/how-can-i-make-geotiffs-of-individual-depth-layers-in-the-hycom-ocean-circulatio" rel="noreferrer" target="_blank">http://gis.stackexchange.com/questions/167155/how-can-i-make-geotiffs-of-individual-depth-layers-in-the-hycom-ocean-circulatio</a><br>
<br>
<br>
<a href="http://gis.stackexchange.com/questions/170882/how-can-a-make-a-netcdf-with-subgroups-smaller" rel="noreferrer" target="_blank">http://gis.stackexchange.com/questions/170882/how-can-a-make-a-netcdf-with-subgroups-smaller</a><br>
<br>
You should start with the regular grids. e.g.<br>
<br>
<a href="ftp://ftp.hycom.org/datasets/GLBu0.08/expt_91.1/data/hindcasts/2015/hycom_glb_911_2015093000_t000_uv3z.nc" rel="noreferrer" target="_blank">ftp://ftp.hycom.org/datasets/GLBu0.08/expt_91.1/data/hindcasts/2015/hycom_glb_911_2015093000_t000_uv3z.nc</a><br>
<br>
Make an 80 layer geotiff from the u/v and then use mapserver as<br>
Even suggested.<br>
<br>
-kurt<br>
<br>
<br>
On Sat, Jun 4, 2016 at 6:29 AM, Even Rouault<br></span>
<<a href="mailto:even.rouault@spatialys.com" target="_blank">even.rouault@spatialys.com</a> <mailto:<a href="mailto:even.rouault@spatialys.com" target="_blank">even.rouault@spatialys.com</a>>><div><div class="h5"><br>
wrote:<br>
<br>
Hi,<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
I want to convert some HYCOM data into GTiff format and I'm a<br>
little lost on how to get started, whether I can do what I need<br>
just using a vrt or if I'm going to have to write code to process<br>
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> and I<br>
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<br>
GTiff and apply color to using color ranges without loosing the<br>
NODATA. I think I can do later using .vrt, correct?<br>
</blockquote>
<br>
gdaldem color-relief can output to a VRT, but given that this netCDF<br>
file is georeferenced through a geolocation array, you'll likely have<br>
first to run gdalwarp -geoloc to have something useful at the end.<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
I'm not sure how to deal with the georeferencing, because the<br>
download site says: "Native hycom .[ab] data converted to NetCDF on<br>
native Mercator-curvilinear HYCOM horizontal grid" but that does<br>
not seem in sync with the Metadata in the file.<br>
<br>
2) The more complicated case is that there are two files: a)<br>
eastward_sea_water_velocity (aka: u, u-velocity) b)<br>
northward_sea_water_velocity (aka: v, v-velocity) that I would like<br>
to use to create a grid of vectors with arrowheads using u and v to<br>
define the vector and where the color of the vector is related to<br>
the magnitude of its length. Ultimately this data will get rendered<br>
via mapserver so I'm wondering if I can do this via a GTiff or if I<br>
will have to resort to creating a massive point shapefile with<br>
attributes of angle and magnitude and use a symbol. It seems like<br>
it would be best if I can keep the data as a GTiff and then sample<br>
the points and render the vectors on the fly for when you zoom<br>
in/out.<br>
</blockquote>
<br>
See <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<br>
same raster.<br>
<br>
Even<br>
<br>
-- Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" rel="noreferrer" target="_blank">http://www.spatialys.com</a><br>
_______________________________________________ gdal-dev mailing<br></div></div>
list <a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a> <mailto:<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a>><span class=""><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
<br>
<br>
<br>
<br>
-- -- <a href="http://schwehr.org" rel="noreferrer" target="_blank">http://schwehr.org</a><br>
</span></blockquote>
<br><div class="HOEnZb"><div class="h5">
<br>
---<br>
This email has been checked for viruses by Avast antivirus software.<br>
<a href="https://www.avast.com/antivirus" rel="noreferrer" target="_blank">https://www.avast.com/antivirus</a><br>
<br>
</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></div>