Hfrhyu

Why is Minecraft giving an OpenGL error?

Can an x86 CPU running in real mode be considered to be basically an 8086 CPU?

Why can't I see bouncing of switch on oscilloscope screen?

When a company launches a new product do they "come out" with a new product or do they "come up" with a new product?

RSA: Danger of using p to create q

Can a monk's single staff be considered dual wielded, as per the Dual Wielder feat?

Is it possible to run Internet Explorer on OS X El Capitan?

Languages that we cannot (dis)prove to be Context-Free

Why is consensus so controversial in Britain?

How much RAM could one put in a typical 80386 setup?

Can a vampire attack twice with their claws using Multiattack?

What's that red-plus icon near a text?

Does detail obscure or enhance action?

Is it possible to do 50 km distance without any previous training?

Important Resources for Dark Age Civilizations?

How can I prevent hyper evolved versions of regular creatures from wiping out their cousins?

Revoked SSL certificate

How do I deal with an unproductive colleague in a small company?

Is it unprofessional to ask if a job posting on GlassDoor is real?

How much of data wrangling is a data scientist's job?

expand `ifthenelse` immediately

Alternative to sending password over mail?

Replacing matching entries in one column of a file by another column from a different file

Filter any system log file by date or date range

Creating NDVI using Python, rasterio producing grey image

Calculating NDVI with Rasteriopython rasterio errorReading S3 Sentinel-2 image files with rasterioNDVI image processingUnable to create NDVI image using Python gdal packagebad results NDVI in rasterio with jp2 sentinel-2 bandsRasterio: tool for creating mosaic?Planet NDVI calculation: ModuleNotFoundError: No module named 'rasterio'Rasterio write_band output has no CRSNDVI Python image showing inverse

.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty margin-bottom:0;

0

Using Sentinel data I am creating an NDVI image using the below Python code. The output is a grey image.

import rasterio
from rasterio import plot
import matplotlib.pyplot as plt
import numpy as np
import subprocess
import os

os.listdir('/Users/abc/S2B_MSIL1C_20190228T050739_N0207_R019_T43PGQ_20190228T085706.SAFE/')

band4 = rasterio.open('/Users/abc/S2B_MSIL1C_20190228T050739_N0207_R019_T43PGQ_20190228T085706.SAFE/GRANULE/L1C_T43PGQ_A010341_20190228T051623/IMG_DATA/T43PGQ_20190228T050739_B04.jp2')
band5 = rasterio.open('/Users/abc/S2B_MSIL1C_20190228T050739_N0207_R019_T43PGQ_20190228T085706.SAFE/GRANULE/L1C_T43PGQ_A010341_20190228T051623/IMG_DATA/T43PGQ_20190228T050739_B08.jp2')

bandRed = band4.read(1)
print(os.getcwd())


resize = np.empty(shape=(1,5490*2,5490*2),dtype=np.uint16)

bandNIR = band5.read(1,out=resize)

bandRed = bandRed.astype('float32')
bandNIR = bandNIR.astype('float32')


ndvi = np.zeros(band4.shape, dtype=rasterio.float32)

ndvi_upper = bandNIR + bandRed
ndvi_lower = bandNIR - bandRed

ndvi = ndvi_lower / ndvi_upper

kwargs = band4.meta
kwargs.update(
 dtype=rasterio.float32,
 count=1,
 compress='lzw')

ndviImage = rasterio.open('/Users/abc/ndviImage.tiff','w',driver='Gtiff',
 width=band4.width,
 height = band4.height,
 count=1, crs=band4.crs,
 transform=band4.transform,
 dtype='float32')


ndviImage.write(ndvi,1)

python ndvi sentinel-2 rasterio

share|improve this question

edited Apr 2 at 12:25

Vince

14.8k32849

asked Apr 2 at 11:46

hck3rhck3r

213

• 
  
  
  

  
  2
  

  
  
  
  
  
  

  
  
  Why not? NDVI is just a single number, in your code ndvi = ndvi_lower / ndvi_upper. It is up to you to apply a color palette when you view the image.
  
  – user30184
  Apr 2 at 12:03
  
  

  
  
  
  

add a comment | 

0

Using Sentinel data I am creating an NDVI image using the below Python code. The output is a grey image.

import rasterio
from rasterio import plot
import matplotlib.pyplot as plt
import numpy as np
import subprocess
import os

os.listdir('/Users/abc/S2B_MSIL1C_20190228T050739_N0207_R019_T43PGQ_20190228T085706.SAFE/')

band4 = rasterio.open('/Users/abc/S2B_MSIL1C_20190228T050739_N0207_R019_T43PGQ_20190228T085706.SAFE/GRANULE/L1C_T43PGQ_A010341_20190228T051623/IMG_DATA/T43PGQ_20190228T050739_B04.jp2')
band5 = rasterio.open('/Users/abc/S2B_MSIL1C_20190228T050739_N0207_R019_T43PGQ_20190228T085706.SAFE/GRANULE/L1C_T43PGQ_A010341_20190228T051623/IMG_DATA/T43PGQ_20190228T050739_B08.jp2')

bandRed = band4.read(1)
print(os.getcwd())


resize = np.empty(shape=(1,5490*2,5490*2),dtype=np.uint16)

bandNIR = band5.read(1,out=resize)

bandRed = bandRed.astype('float32')
bandNIR = bandNIR.astype('float32')


ndvi = np.zeros(band4.shape, dtype=rasterio.float32)

ndvi_upper = bandNIR + bandRed
ndvi_lower = bandNIR - bandRed

ndvi = ndvi_lower / ndvi_upper

kwargs = band4.meta
kwargs.update(
 dtype=rasterio.float32,
 count=1,
 compress='lzw')

ndviImage = rasterio.open('/Users/abc/ndviImage.tiff','w',driver='Gtiff',
 width=band4.width,
 height = band4.height,
 count=1, crs=band4.crs,
 transform=band4.transform,
 dtype='float32')


ndviImage.write(ndvi,1)

python ndvi sentinel-2 rasterio

share|improve this question

edited Apr 2 at 12:25

Vince

14.8k32849

asked Apr 2 at 11:46

hck3rhck3r

213

• 
  
  
  

  
  2
  

  
  
  
  
  
  

  
  
  Why not? NDVI is just a single number, in your code ndvi = ndvi_lower / ndvi_upper. It is up to you to apply a color palette when you view the image.
  
  – user30184
  Apr 2 at 12:03
  
  

  
  
  
  

add a comment | 

Using Sentinel data I am creating an NDVI image using the below Python code. The output is a grey image.

import rasterio
from rasterio import plot
import matplotlib.pyplot as plt
import numpy as np
import subprocess
import os

os.listdir('/Users/abc/S2B_MSIL1C_20190228T050739_N0207_R019_T43PGQ_20190228T085706.SAFE/')

band4 = rasterio.open('/Users/abc/S2B_MSIL1C_20190228T050739_N0207_R019_T43PGQ_20190228T085706.SAFE/GRANULE/L1C_T43PGQ_A010341_20190228T051623/IMG_DATA/T43PGQ_20190228T050739_B04.jp2')
band5 = rasterio.open('/Users/abc/S2B_MSIL1C_20190228T050739_N0207_R019_T43PGQ_20190228T085706.SAFE/GRANULE/L1C_T43PGQ_A010341_20190228T051623/IMG_DATA/T43PGQ_20190228T050739_B08.jp2')

bandRed = band4.read(1)
print(os.getcwd())


resize = np.empty(shape=(1,5490*2,5490*2),dtype=np.uint16)

bandNIR = band5.read(1,out=resize)

bandRed = bandRed.astype('float32')
bandNIR = bandNIR.astype('float32')


ndvi = np.zeros(band4.shape, dtype=rasterio.float32)

ndvi_upper = bandNIR + bandRed
ndvi_lower = bandNIR - bandRed

ndvi = ndvi_lower / ndvi_upper

kwargs = band4.meta
kwargs.update(
 dtype=rasterio.float32,
 count=1,
 compress='lzw')

ndviImage = rasterio.open('/Users/abc/ndviImage.tiff','w',driver='Gtiff',
 width=band4.width,
 height = band4.height,
 count=1, crs=band4.crs,
 transform=band4.transform,
 dtype='float32')


ndviImage.write(ndvi,1)

python ndvi sentinel-2 rasterio

share|improve this question

edited Apr 2 at 12:25

Vince

14.8k32849

asked Apr 2 at 11:46

hck3rhck3r

213

import rasterio
from rasterio import plot
import matplotlib.pyplot as plt
import numpy as np
import subprocess
import os

os.listdir('/Users/abc/S2B_MSIL1C_20190228T050739_N0207_R019_T43PGQ_20190228T085706.SAFE/')

band4 = rasterio.open('/Users/abc/S2B_MSIL1C_20190228T050739_N0207_R019_T43PGQ_20190228T085706.SAFE/GRANULE/L1C_T43PGQ_A010341_20190228T051623/IMG_DATA/T43PGQ_20190228T050739_B04.jp2')
band5 = rasterio.open('/Users/abc/S2B_MSIL1C_20190228T050739_N0207_R019_T43PGQ_20190228T085706.SAFE/GRANULE/L1C_T43PGQ_A010341_20190228T051623/IMG_DATA/T43PGQ_20190228T050739_B08.jp2')

bandRed = band4.read(1)
print(os.getcwd())


resize = np.empty(shape=(1,5490*2,5490*2),dtype=np.uint16)

bandNIR = band5.read(1,out=resize)

bandRed = bandRed.astype('float32')
bandNIR = bandNIR.astype('float32')


ndvi = np.zeros(band4.shape, dtype=rasterio.float32)

ndvi_upper = bandNIR + bandRed
ndvi_lower = bandNIR - bandRed

ndvi = ndvi_lower / ndvi_upper

kwargs = band4.meta
kwargs.update(
 dtype=rasterio.float32,
 count=1,
 compress='lzw')

ndviImage = rasterio.open('/Users/abc/ndviImage.tiff','w',driver='Gtiff',
 width=band4.width,
 height = band4.height,
 count=1, crs=band4.crs,
 transform=band4.transform,
 dtype='float32')


ndviImage.write(ndvi,1)

share|improve this question

edited Apr 2 at 12:25

Vince

14.8k32849

asked Apr 2 at 11:46

hck3rhck3r

213

share|improve this question

edited Apr 2 at 12:25

Vince

14.8k32849

asked Apr 2 at 11:46

hck3rhck3r

213

edited Apr 2 at 12:25

Vince

14.8k32849

edited Apr 2 at 12:25

Vince

14.8k32849

asked Apr 2 at 11:46

hck3rhck3r

213

asked Apr 2 at 11:46

hck3rhck3r

213