Use satsense for Sentinel 2 tif about satsense HOT 3 CLOSED

import matplotlib.pyplot as plt

from satsense import Image
from satsense.generators import FullGenerator
from satsense.features import NirNDVI
from satsense import extract_features

# URI to the image
imagefile = 'path/to/your/sentinel2_image.tif'

# Loading the file
image = Image(imagefile, 'quickbird')
image.precompute_normalization()

# Define the step from one window to the next
step_size = [10, 10]
# Define one or more window sizes 
windows=[(10, 10)]

# Define the features to calculate
features = [
  NirNDVI(windows)
]
generator = FullGenerator(image, step_size)

# Extract the features (using 1 cpu)
# extract_features is a generator, therefore use the next() function
feature_vector = next(extract_features(features, generator, n_jobs=1))

plt.imshow(feature_vector.vector[:, :, 0, 0])