def test_median_reduction_over_latlong_old_version(mock_api):
# Test median reduction over lat/long - old version for backwards compatibility
a = AnalyticsEngine(api=mock_api)
e = ExecutionEngine(api=mock_api)
# 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))
}
}
arrays = a.create_array(('LANDSAT_5', 'NBAR'), ['band_40'], dimensions,
'get_data')
median_xy = a.apply_generic_reduction(arrays, ['x', 'y'], 'median(array1)',
'medianXY')
result = e.execute_plan(a.plan)
def main():
a = AnalyticsEngine()
e = ExecutionEngine()
# 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))
}
}
arrays = a.create_array(('LANDSAT_5', 'nbar'), ['nir'], dimensions,
'get_data')
median_xy = a.apply_generic_reduction(arrays, ['y', 'x'], 'median(array1)',
'medianXY')
result = e.execute_plan(a.plan)
plot(e.cache['medianXY'])
def main():
a = AnalyticsEngine()
e = ExecutionEngine()
# 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))}}
arrays = a.create_array(('LANDSAT_5', 'nbar'), ['nir'], dimensions, 'get_data')
median_t = a.apply_generic_reduction(arrays, ['time'], 'median(array1)', 'medianT')
result = e.execute_plan(a.plan)
plot(e.cache['medianT'])
b40_result = e.cache['get_data']['array_result']['nir']
median_result = e.cache['medianT']['array_result']['medianT']
b40_data = Data(x=b40_result[:, ::-1, :], label='B40')
median_data = Data(x=median_result[::-1, :], label='medianT')
long_data = Data(x=b40_result.coords['x'], label='long')
lat_data = Data(x=b40_result.coords['y'], label='lat')
time_data = Data(x=b40_result.coords['time'], label='time')
collection = DataCollection([median_data, b40_data, long_data, lat_data, time_data, ])
app = GlueApplication(collection)
app.start()
def test_median_reduction_over_latlong_old_version(mock_api):
# Test median reduction over lat/long - old version for backwards compatibility
a = AnalyticsEngine(api=mock_api)
e = ExecutionEngine(api=mock_api)
# 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))}}
arrays = a.create_array(('LANDSAT_5', 'NBAR'), ['band_40'], dimensions, 'get_data')
median_xy = a.apply_generic_reduction(arrays, ['x', 'y'], 'median(array1)', 'medianXY')
result = e.execute_plan(a.plan)
def main():
a = AnalyticsEngine()
e = ExecutionEngine()
# 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))}}
arrays = a.create_array(('LANDSAT_5', 'nbar'), ['nir'], dimensions, 'get_data')
median_xy = a.apply_generic_reduction(arrays, ['y', 'x'], 'median(array1)', 'medianXY')
result = e.execute_plan(a.plan)
plot(e.cache['medianXY'])
def main():
a = AnalyticsEngine()
e = ExecutionEngine()
# 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))
}
}
arrays = a.create_array(('LANDSAT_5', 'nbar'), ['nir'], dimensions,
'get_data')
median_t = a.apply_generic_reduction(arrays, ['time'], 'median(array1)',
'medianT')
result = e.execute_plan(a.plan)
plot(e.cache['medianT'])
b40_result = e.cache['get_data']['array_result']['nir']
median_result = e.cache['medianT']['array_result']['medianT']
b40_data = Data(x=b40_result[:, ::-1, :], label='B40')
median_data = Data(x=median_result[::-1, :], label='medianT')
long_data = Data(x=b40_result.coords['x'], label='long')
lat_data = Data(x=b40_result.coords['y'], label='lat')
time_data = Data(x=b40_result.coords['time'], label='time')
collection = DataCollection([
median_data,
b40_data,
long_data,
lat_data,
time_data,
])
app = GlueApplication(collection)
app.start()