In [18]:
import rasterio
import numpy as np
In [19]:
path = './dummy.asc'
with open(path) as f:
for l in f:
print l,
In [50]:
ds = rasterio.open(path)
out = np.ma.zeros((5,4), dtype=ds.dtypes[0])
returned = ds.read_band(1, out=out, masked=True)
In [58]:
out
Out[58]:
In [52]:
returned
Out[52]:
In [32]:
a = np.array([1.1,2.2,3.3,4.4,5.5,6.6,7.7,8.8])
print np.around(a)
out = np.zeros((8))
print out
result = np.around(a,out=out)
print result
print out
In [54]:
ma = np.ma.array([1,2,3])
In [57]:
ma.fill_value = -9999
ma
Out[57]:
In [62]:
out.mask = out == ma.fill_value
In [3]:
import numpy as np
a = np.ma.array([1, 2, 3, 4], mask=[True, False, False, True])
b = np.ma.array([1, 2, 3, 4], mask=[False, True, False, True])
# ufunc not in-place
c = np.add(a, b)
# object operator
d = a + b
# ufunc in_place
np.add(a, b, out=a)
print 'not in-place:\n{}\n{}\n'.format(repr(c), c.data)
print 'object operator:\n{}\n{}\n'.format(repr(d), d.data)
print 'in-place:\n{}\n{}\n'.format(repr(a), a.data)
In [2]:
Out[2]:
In [ ]: