octave-cvsupdate Mailing List for Octave Forge (Page 18)

SourceForge Headquarters 225 Broadway Suite 1600 San Diego, CA 92101 +1 (858) 422-6466

Revision: 10673
          http://octave.svn.sourceforge.net/octave/?rev=10673&view=rev
Author:   abarth93
Date:     2012-06-22 15:52:36 +0000 (Fri, 22 Jun 2012)
Log Message:
-----------


Added Paths:
-----------
    trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/
    trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/display.m
    trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/ncBaseArray.m
    trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/ncsub.m
    trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/subsasgn.m
    trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/subsref.m

Added: trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/display.m
===================================================================

--- trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/display.m	                        (rev 0)
+++ trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/display.m	2012-06-22 15:52:36 UTC (rev 10673)
@@ -0,0 +1,4 @@
+function display(self)
+disp([self.varname ' from ' self.filename]);
+%self.vinfo
+%self.dims

Added: trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/ncBaseArray.m
===================================================================
--- trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/ncBaseArray.m	                        (rev 0)
+++ trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/ncBaseArray.m	2012-06-22 15:52:36 UTC (rev 10673)
@@ -0,0 +1,62 @@
+% V = ncBaseArray(filename,varname)
+% V = ncBaseArray(filename,varname,'property',value,...)
+% create a ncBaseArray that can be accessed as a normal matlab array.
+% Ths object is a helper object. Users should normally call ncArray.
+%
+% For read access filename can be compressed if it has the extensions
+% ".gz" or ".bz2". It use the function cache_decompress to cache to
+% decompressed files.
+%
+% Data is loaded by ncread and saved by ncwrite. Values equal to _FillValue
+% are thus replaced by NaN and the scaling (add_offset and
+% scale_factor) is applied during loading and saving.
+%
+% Properties:
+%   'tooBigToLoad': if tooBigToLoad is set to true, then only the minimum
+%   data will be loaded. However this can be quite slow.
+%
+% Example:
+%
+% Loading the variable (assuming V is 3 dimensional):
+%
+% x = V(1,1,1);   % load element 1,1,1
+% x = V(:,:,:);   % load the complete array
+% x = V();  x = full(V)  % loads also the complete array
+%
+% Saving data in the netcdf file:
+% V(1,1,1) = x;   % save x in element 1,1,1
+% V(:,:,:) = x;
+%
+% Attributes
+% units = V.units; % get attribute called "units"
+% V.units = 'degree C'; % set attributes;
+%
+% Note: use the '.()' notation if the attribute has a leading underscore 
+% (due to a limitation in the matlab parser):
+%
+% V.('_someStrangeAttribute') = 123;
+%
+% see also cache_decompress 
+
+function retval = ncBaseArray(filename,varname,varargin)
+
+self.tooBigToLoad = false;
+prop = varargin;
+
+for i=1:2:length(prop)
+    if strcmp(prop{i},'tooBigToLoad ')
+        self.tooBigToLoad = prop{i+1};
+    end
+end
+
+self.filename = filename;
+self.varname = varname;
+
+self.vinfo = ncinfo(cached_decompress(filename),varname);
+self.sz = self.vinfo.Size;
+
+self.dims = self.vinfo.Dimensions;
+self.nd = length(self.dims);
+
+retval = class(self,'ncBaseArray',BaseArray(self.sz));
+end

Added: trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/ncsub.m
===================================================================
--- trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/ncsub.m	                        (rev 0)
+++ trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/ncsub.m	2012-06-22 15:52:36 UTC (rev 10673)
@@ -0,0 +1,61 @@
+function [start, count, stride] = ncsub(self,idx)
+
+assert(strcmp(idx.type,'()'))
+
+% number of dimension (including singleton dimensions)
+%n = length(size(self));
+n = self.nd;
+
+% number of subscripts
+ns = length(idx.subs);
+
+if ns == 0
+    % load all
+    start = ones(1,n);
+    count = self.sz;
+    stride = ones(1,n);
+else
+    
+    start = ones(1,ns);
+    count = ones(1,ns);
+    stride = ones(1,ns);
+    
+    % sz is the size padded by 1 if more indices are given than n
+    sz = ones(1,ns);
+    sz(1:length(self.sz)) = self.sz;
+    
+    for i=1:ns
+        if isempty(idx.subs{i})
+            count(i) = 0;
+            
+        elseif strcmp(idx.subs{i},':')
+            count(i) = sz(i);
+            
+        else
+            tmp = idx.subs{i};
+            
+            if length(tmp) == 1
+                start(i) = tmp;
+            else
+                test = tmp(1):tmp(2)-tmp(1):tmp(end);
+                
+                if all(tmp == test)
+                    start(i) = tmp(1);
+                    stride(i) = tmp(2)-tmp(1);
+                    count(i) = (tmp(end)-tmp(1))/stride(i) +1;
+                else
+                    error('indeces');
+                end
+            end
+        end
+    end
+    
+    assert(all(count(n+1:end) == 1 | count(n+1:end) == 0))
+    assert(all(start(n+1:end) == 1))
+    
+    if ~any(count == 0)
+        count = count(1:n);
+        start = start(1:n);
+        stride = stride(1:n);
+    end
+end
\ No newline at end of file

Added: trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/subsasgn.m
===================================================================
--- trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/subsasgn.m	                        (rev 0)
+++ trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/subsasgn.m	2012-06-22 15:52:36 UTC (rev 10673)
@@ -0,0 +1,10 @@
+function self = subsasgn(self,idxs,x)
+%idx
+
+idx = idxs(1);
+assert(strcmp(idx.type,'()'))
+[start, count, stride] = ncsub(self,idx);
+
+if all(count ~= 0)
+    ncwrite(self.filename,self.varname,x,start,stride);
+end

Added: trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/subsref.m
===================================================================
--- trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/subsref.m	                        (rev 0)
+++ trunk/octave-forge/extra/ncArray/inst/@ncBaseArray/subsref.m	2012-06-22 15:52:36 UTC (rev 10673)
@@ -0,0 +1,65 @@
+function x = subsref(self,idx)
+
+assert(length(idx) == 1)
+
+if strcmp(idx.type,'()')
+    % load data
+    
+    if strcmp(idx.type,'()') && length(idx.subs) == 1 ...
+            && length(idx.subs) < self.nd
+        % reference like A([2 3 1 5])
+        
+        if self.tooBigToLoad                                   
+            % number of elements in x
+            ind = idx.subs{1};
+            
+            % transform index to subscript
+            subs = cell(1,self.nd);
+            [subs{:}] = ind2sub(size(self),ind);
+            
+            % make a nice array length(ind) by self.nd
+            subs = cell2mat(subs);
+            
+            % output array
+            x = zeros(size(ind));
+            x(:) = NaN;
+            
+            
+            % get every element
+            % can be quite slow
+            idxe.type = '()';
+            idxe.subs = cell(1,self.nd);
+            
+            for i=1:length(ind)
+                idxe.subs = mat2cell(subs(i,:),1,ones(1,self.nd));
+                x(i) = subsref(self,idxe);
+            end
+        else
+            % load full array
+            tmp = full(self);
+            x = subsref(tmp,idx);
+        end
+    else
+        [start, count, stride] = ncsub(self,idx);
+        
+        if any(count == 0)
+            x = zeros(count);
+        else
+            x = ncread(cached_decompress(self.filename),self.varname,...
+                start,count,stride);
+        end
+    end
+elseif strcmp(idx.type,'.')
+    % load attribute
+    name = idx.subs;
+    index = strmatch(name,{self.vinfo.Attributes(:).Name});
+    
+    if isempty(index)
+        error('variable %s has no attribute called %s',self.varname,name);
+    else
+        x = self.vinfo.Attributes(index).Value;
+    end
+else
+    error('not supported');
+    
+end
\ No newline at end of file

This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.





2001	Jan	Feb	Mar	Apr	May	Jun	Jul (10)	Aug (5)	Sep (3)	Oct (41)	Nov (41)	Dec (33)
2002	Jan (75)	Feb (10)	Mar (170)	Apr (174)	May (66)	Jun (11)	Jul (10)	Aug (44)	Sep (73)	Oct (28)	Nov (139)	Dec (52)
2003	Jan (35)	Feb (93)	Mar (62)	Apr (10)	May (55)	Jun (70)	Jul (37)	Aug (16)	Sep (56)	Oct (31)	Nov (57)	Dec (83)
2004	Jan (85)	Feb (67)	Mar (27)	Apr (37)	May (75)	Jun (85)	Jul (160)	Aug (68)	Sep (104)	Oct (25)	Nov (39)	Dec (23)
2005	Jan (10)	Feb (45)	Mar (43)	Apr (19)	May (108)	Jun (31)	Jul (41)	Aug (23)	Sep (65)	Oct (58)	Nov (44)	Dec (54)
2006	Jan (96)	Feb (27)	Mar (69)	Apr (59)	May (67)	Jun (35)	Jul (13)	Aug (461)	Sep (160)	Oct (399)	Nov (32)	Dec (72)
2007	Jan (316)	Feb (305)	Mar (318)	Apr (54)	May (194)	Jun (173)	Jul (282)	Aug (91)	Sep (227)	Oct (365)	Nov (168)	Dec (18)
2008	Jan (71)	Feb (111)	Mar (155)	Apr (173)	May (70)	Jun (67)	Jul (55)	Aug (83)	Sep (32)	Oct (68)	Nov (80)	Dec (29)
2009	Jan (46)	Feb (18)	Mar (95)	Apr (76)	May (140)	Jun (98)	Jul (84)	Aug (123)	Sep (94)	Oct (131)	Nov (142)	Dec (125)
2010	Jan (128)	Feb (158)	Mar (172)	Apr (134)	May (94)	Jun (84)	Jul (32)	Aug (127)	Sep (167)	Oct (109)	Nov (69)	Dec (78)
2011	Jan (39)	Feb (58)	Mar (52)	Apr (47)	May (56)	Jun (76)	Jul (55)	Aug (54)	Sep (165)	Oct (255)	Nov (328)	Dec (263)
2012	Jan (82)	Feb (147)	Mar (400)	Apr (216)	May (209)	Jun (160)	Jul (86)	Aug (141)	Sep (156)	Oct (6)	Nov	Dec
2015	Jan	Feb	Mar	Apr	May	Jun	Jul (1)	Aug	Sep (1)	Oct	Nov (1)	Dec (2)
2016	Jan	Feb (2)	Mar (2)	Apr (1)	May (1)	Jun (2)	Jul (1)	Aug (1)	Sep	Oct	Nov (1)	Dec
2019	Jan	Feb	Mar (1)	Apr	May	Jun	Jul	Aug	Sep	Oct	Nov	Dec
2021	Jan	Feb	Mar	Apr (3)	May (4)	Jun (8)	Jul (2)	Aug (5)	Sep (9)	Oct	Nov	Dec

octave-cvsupdate Mailing List for Octave Forge (Page 18)

A collection of packages providing extra functionality for GNU Octave

octave-cvsupdate — CVS update notifications