Hi Richard,
Scilab 5.5.2: had to use stacksize('max'), data loaded but surf would not
display 866 x 654 matrix
Scilab 6: it did not work fine with one of the following syntaxes:
asc_data(asc_data==-999) = %nan
surf(asc_data);
or with:
ix = find(asc_data==-999)
asc_data(ix) = %nan
surf(asc_data);
Note that as shown by Samuel and in the test I produced yesterday, the syntax
in your original code produces wrong results:
[a,b] = find(asc_data==-999);
asc_data(a,b) = %nan;
To use the indices a,b this should be modified as:
[a,b] = find(asc_data==-999);
asc_data([a;b]) = %nan;
However this last syntax seemed to work on yesterday's test on a much smaller
matrix but it is not working on your big one... Not sure if this is a bug or
problem with the syntax, needs further investigation.
Regards,
Rafael
-----Original Message-----
From: users [mailto:users-boun...@lists.scilab.org] On Behalf Of Richard llom
Sent: Wednesday, November 22, 2017 9:33 AM
To: users@lists.scilab.org
Subject: Re: [Scilab-users] surf and %nan
Hey Rafael & Samuel,
thanks for the replies.
However I'm really having a hard time to process this data.
This is the raw data in question:
ftp://ftp-cdc.dwd.de/pub/CDC/grids_germany/monthly/radiation_diffuse/grids_germany_monthly_radiation_diffuse_201601.zip
My minimal script:
asc_list = listfiles('*.asc');
asc_idx = 1;
asc_file = asc_list(asc_idx);
[asc_data, comments] = fscanfMat(asc_file);
[a,b] = find(asc_data==-999);
asc_data(a,b) = %nan;
surf(asc_data);
Right now it even fails at processing the -999 to nans... It's like I'm
hitting bugs in scilabs?
regards
richard
_______________________________________________
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users
