[Scilab-users] Getting netCDF files into Scilab 6.0.1
Samuel Gougeon
sgougeon at free.fr
Sun Jul 22 11:45:57 CEST 2018
Hello Lester,
Le 22/07/2018 à 11:12, arctica1963 a écrit :
> Hello,
>
> I am looking to use netCDF files generated from the GMT software in Scilab
> 6.0.1, but there is no easy way to work the problem other than exporting the
> data out from the binary grids to either XYZ or ESRI Ascii raster. The
> latter is a more compact form with the basic structure of a header
> describing the data limits and structure:
>
> ncols 1200
> nrows 1200
> xllcorner -10 Lower-left X (Long-Lat limits -10/10/-10/10
> yllcorner -10 Lower-Left Y
> cellsize 0.0166666666667 - increment = 1 arc minute
> nodata_value -9999
> 451.343170166 436.005554199 443.061035156 443.665924072 465.607574463
> 492.191741943 476.50994873 452.265014648 451.439880371 461.659393311 ......
> ...
> ... all of the Z-values
> ...
>
> -4323.2890625 -4315.12451172 -4307.19384766 -4305.27490234 -4311.49902344
> -4317.90771484 -4324.23046875 -4320.47070313 -4308.60107422 -4292.62011719
> -4280.39697266
>
> This is a moderately large file ~20 mB or so. Is there a way to reliably
> read this structure and vectorise X, Y, and Z? The data I am looking at is a
> grid so want to work on the whole grid as eith DEM data or gravity etc.
>
> As a test, I converted the grid to xyz, and tried csvRead, having replaced
> the spaces (which were not a constant number of spaces between values) with
> commas. The file was in excess of 50 mb, so not that huge really, but
> csvRead always failed to read all the file and drop out at one specific line
> (~7608) with a message of inconsistent column numbers, i.e. it though there
> were only 2 columns when in reality the whole thing was 3 (checked in the
> program Surfer). When one looks at the line specified, it is clearly 3
> columns.
Do you really get this error with the new nightly built Scilab release
http://www.scilab.org/en/development/nightly_builds/branch60
indicated by Antoine E. yesterday?
It is due to the same mgetl() bug (for Windows), that spreads all over
Scilab 6.0.1.
The continuous build
https://build.scilab.org/view/Scilab%206.0/job/scilab-6.0-windows-64/
can also be used without waiting for any publication delay.
Samuel
More information about the users
mailing list