Wsrtjtyk

I am trying to plot a raster image using Rasterio, Matplotlib, and Numpy. Although what I am doing is probably irrelevant and general to the actual python code and error.

Here is the code:

import rasterio

import matplotlib.pyplot as plt

import numpy as np



band5 = rasterio.open(r"C:UsersnewDesktopLC08_L1TP_015033_20170822_20170912_01_T1_B5.TIF")



band4 = rasterio.open(r"LC08_L1TP_015033_20170822_20170912_01_T1_B4.TIF")



#rasterio.windows.Window(col_off, row_off, width, height)

window = rasterio.windows.Window(1024, 1024, 1280, 2560)



with rasterio.open(band5) as src:

    subset = src.read(1, window=window)



plt.figure(figsize=(6,8.5))

plt.imshow(subset)

plt.colorbar(shrink=0.5)

plt.title(f'Band 5 Subsetn{window}')

plt.xlabel('Column #')

plt.ylabel('Row #')

Complete error message with traceback:

runfile('C:/Users/new/Desktop/RasterSubsetNDVI.py', wdir='C:/Users/new/Desktop')

Traceback (most recent call last):



  File "<ipython-input-131-5d04fa0ce75f>", line 1, in <module>

    runfile('C:/Users/new/Desktop/RasterSubsetNDVI.py', wdir='C:/Users/new/Desktop')



  File "C:UsersnewMiniconda3libsite-packagesspyder_kernelscustomizespydercustomize.py", line 668, in runfile

    execfile(filename, namespace)



  File "C:UsersnewMiniconda3libsite-packagesspyder_kernelscustomizespydercustomize.py", line 108, in execfile

    exec(compile(f.read(), filename, 'exec'), namespace)



  File "C:/Users/new/Desktop/RasterSubsetNDVI.py", line 12, in <module>

    with rasterio.open(band5) as src:



  File "C:UsersnewMiniconda3libcontextlib.py", line 81, in __enter__

    return next(self.gen)



  File "C:UsersnewMiniconda3libsite-packagesrasterio__init__.py", line 177, in fp_reader

    memfile = MemoryFile(fp.read())



  File "C:UsersnewMiniconda3libsite-packagesrasterioio.py", line 105, in __init__

    file_or_bytes=file_or_bytes, filename=filename, ext=ext)



  File "rasterio/_io.pyx", line 745, in rasterio._io.MemoryFileBase.__init__



ValueError: The truth value of an array with more than one element is ambiguous. Use a.any() or a.all()



    ValueError: The truth value of an array with more than one element is ambiguous. Use a.any() or a.all()

I am using code found here and trying to reproduce it on my machine with my own TIF files. Where it says "Pull an image subset at full resolution". LINK:
https://github.com/geohackweek/tutorial_contents/blob/ba5e9443137a9aca87cdcdcd70e9e6a237cc64ba/raster/notebooks/rasterio-landsat-aws.ipynb

I am running on the IDE Sypder and using Python 3.6

Please, any help is appreciated.

Thank you

Take a careful look at the following lines of your code (others omitted):

band5 = rasterio.open(r"C:UsersnewDesktopLC08_L1TP_015033_20170822_20170912_01_T1_B5.TIF")



# ...



with rasterio.open(band5) as src:

    # ...

From the Rasterio documentation:

Rasterio’s rasterio.open() takes a path and returns a dataset object.

In the second line above, you're passing a dataset to rasterio.open() instead of a path. Perhaps you wanted to pass the same path you used to create band5 from?

It would have been nice if rasterio.open() checked the type of the argument you gave it and gave you an exception that said 'file path is not a string' or suchlike, but sadly it appears not to do this.