Earthpy

# Get data for example
data = et.data.get_data("vignette-landsat")

# Set working directory
os.chdir(os.path.join(et.io.HOME, "earth-analytics"))

# Stack the Landsat 8 bands
# This creates a numpy array with each "layer" representing a single band
# You can use the nodata= parameter to mask nodata values
landsat_path = glob(
    os.path.join(
        "data",
        "vignette-landsat",
        "LC08_L1TP_034032_20160621_20170221_01_T1_sr_band*_crop.tif",
    )
)
landsat_path.sort()
array_stack, meta_data = es.stack(landsat_path, nodata=-9999)

Source: earthpy.readthedocs.io

Similar codes

0 vote179 views1 answer
How to set text in QComboBox pyqt5
Python
0 vote171 views1 answer
Input code for python
Python
0 vote177 views1 answer
Fforeveer loop python
Python
0 vote168 views1 answer
Python replace list of ips from yaml file with new list
Python
0 vote164 views1 answer
Randint without repitition
Python
0 vote175 views1 answer
Pandas read csv encoding thai
Python
11 votes181 views3 answers
Django q filter
Python Django
0 vote185 views2 answers
Doormat pattern python
Python
38 votes199 views6 answers
How to create virtual environment in anaconda
Python
0 vote176 views1 answer
90/360
Python