[Scilab-users] surf and %nan

[Rafael Guerra]

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-bounces at lists.scilab.org] On Behalf Of Richard llom
Sent: Wednesday, November 22, 2017 9:33 AM
To: users at 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


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20171122/d19cd3cd/attachment.htm>