[gdal-dev] Writing Pixel Functions with VRT

[Brian Case]

Yves,

You need to be write a custom application to make use of
<PixelFunctionType>

If you are looking for a command line utility solution you should use
gdal_calc.py


if your end goal is to serve this data via a wms you could try the
thredds server
http://www.unidata.ucar.edu/projects/THREDDS/tech/TDS.html 

http://atmos.telascience.org:8080/thredds/wms/ncom_fukushima_agg/Fukushima_best.ncd?LAYERS=sea_water_velocity&ELEVATION=0&TIME=2011-12-17T00%3A00%3A00.000Z&TRANSPARENT=true&STYLES=vector%2Frainbow&CRS=EPSG%3A4326&COLORSCALERANGE=0%2C1&NUMCOLORBANDS=64&LOGSCALE=false&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&EXCEPTIONS=XML&FORMAT=image%2Fpng&BBOX=135.74485290527,32.501498222351,147.25334655761,41.492508888245&WIDTH=768&HEIGHT=600

On Wed, 2011-12-21 at 20:41 +0100, Yves Jacolin (free) wrote:
> Hello,
> 
> I have some question about writing pixel functions with vrt files. I have a 
> netcdf file with two subset named U and V.
> 
> How can I process sqrt(U^2 + V^2) using this feature with vrt file?
> 
> * do I need to create a C file and build GADL with it?
> * do I need to create a C file and locate it close to vrt file?
> 
> Thanks,
> 
> Y.