In [1]:

    

# NDVI

%matplotlib inline
import mpld3; mpld3.enable_notebook()

from datetime import datetime
from datacube.analytics.analytics_engine import AnalyticsEngine
from datacube.execution.execution_engine import ExecutionEngine
from datacube.analytics.utils.analytics_utils import plot


    

    




/home/rc/miniconda/lib/python2.7/site-packages/datacube-1.0.3+70.gc24cf61.dirty-py2.7.egg/datacube/index/postgres/tables/_dataset.py:11: SADeprecationWarning: The SQLAlchemy PostgreSQL dialect has been renamed from 'postgres' to 'postgresql'. The new URL format is postgresql[+driver]://<user>:<pass>@<host>/<dbname>
  from sqlalchemy.dialects import postgres
/home/rc/miniconda/lib/python2.7/site-packages/osgeo/gdal.py:80: DeprecationWarning: gdal.py was placed in a namespace, it is now available as osgeo.gdal
  DeprecationWarning)

In [2]:

    

# initialise AE/EE

a = AnalyticsEngine()
e = ExecutionEngine()


    

    




/home/rc/miniconda/lib/python2.7/site-packages/datacube-1.0.3+70.gc24cf61.dirty-py2.7.egg/datacube/api/_api.py:64: UserWarning: No application name set during API initialization
  warnings.warn('No application name set during API initialization')

In [3]:

    

# Lake Burley Griffin
dimensions = {'x':    {'range': (149.07, 149.18)},
              'y':    {'range': (-35.32, -35.28)},
              'time': {'range': (datetime(1990, 1, 1), datetime(1990, 12, 31))}}

b40 = a.create_array(('LANDSAT_5', 'nbar'), ['nir'], dimensions, 'b40')
b30 = a.create_array(('LANDSAT_5', 'nbar'), ['red'], dimensions, 'b30')

ndvi = a.apply_expression([b40, b30], '((array1 - array2) / (array1 + array2))', 'ndvi')


    

In [4]:

    

e.execute_plan(a.plan)


    

In [5]:

    

plot(e.cache['b40'])


    

    




plot3D






    




/home/rc/miniconda/lib/python2.7/site-packages/IPython/core/formatters.py:92: DeprecationWarning: DisplayFormatter._ipython_display_formatter_default is deprecated: use @default decorator instead.
  def _ipython_display_formatter_default(self):
/home/rc/miniconda/lib/python2.7/site-packages/IPython/core/formatters.py:669: DeprecationWarning: PlainTextFormatter._singleton_printers_default is deprecated: use @default decorator instead.
  def _singleton_printers_default(self):






    

In [6]:

    

plot(e.cache['b30'])


    

    




plot3D






    

In [7]:

    

plot(e.cache['ndvi'])


    

    




plot3D






    

In [ ]: