Yup, forgot to enable netcdf. Sorry about that.

Otherwise your directions worked perfectly, thank you very much!

In terms of understanding, exactly, these steps for the next time:

0) How would I have known the data were transposed lon, lat?
0b) How would I have known the spacing of the lat variable was not constant? I 
imagine "because I looked at it?" :) Knowing this, I've gone back and I see it 
now when I do a quick analysis with xarray. What was the method you diagnosed 
this with?
1) What is the motivation to tiling the output?
2) Knowing about the mirror command is cool
3) I don't understand this step.  Why is this referencing back  the nc file?

Cheers
Chris


Chris Marsh, PhD
University of Saskatchewan
chrismarsh.ca<http://chrismarsh.ca/>



On Sun, Feb 2, 2020 at 6:15 AM Even Rouault 
<even.roua...@spatialys.com<mailto:even.roua...@spatialys.com>> wrote:
CAUTION: This email originated from outside of the University of Saskatchewan. 
Do not click links or open attachments unless you recognize the sender and know 
the content is safe. If in doubt, please forward suspicious emails to 
phish...@usask.ca<mailto:phish...@usask.ca>


On dimanche 2 février 2020 00:34:46 CET Marsh, Chris wrote:
> HI Even,
>
> Thanks for the detailed response. I have compiled master 3.1dev however
> strangest thing when I try to run it

> $ export PATH=`pwd`/gdal-install/bin:$PATH
> $ which gdal_translate
> /Users/chris/Documents/science/analysis/gdal-install/bin/gdal_translate
>  $ gdal_translate NETCDF:"`pwd`/SD_20160901.nc":snd_upd out.tif -co
> TILED=YES
 ERROR 4:
> NETCDF:/Users/chris/Documents/science/analysis/SD_20160901.nc:snd_upd: No
> such file or directory

Do you have the netCDF driver included ? Check with
"gdalinfo --format netCDF". This file is also a netCDF 4 file, so you need to
make sure you netCDF lib is built against HDF5.
You may use the "osgeo/gdal" Docker image.
See https://github.com/OSGeo/gdal/blob/master/gdal/docker/README.md

> $ ls
> /Users/chris/Documents/science/analysis/SD_20160901.nc
> /Users/chris/Documents/science/analysis/SD_20160901.nc
>
> If I try to use my homebrew installed gdal_translate it works.
>
> I've also tried without the `pwd` but that doesn't work either. This seems
> like a stupid PEBKAC but I'm out of ideas.

> Trying the 2nd command resulted in
>
> $ /gdal-install/bin/gdalmdimtranslate out.tif out2.tif -co TILED=YES -array
> band=1,transpose=[1,0],view=[::-1,::-1]
 zsh: no matches found:

Ah, just trying with zsh, this is just a quoting issue. You need to surround
the value of the array option with double quotes

... -array "band=1,transpose=[1,0],view=[::-1,::-1]"

--
Spatialys - Geospatial professional services
http://www.spatialys.com
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

Reply via email to