Skip to main content

Processing satellite data directly from TAR files.

Project description

Python-for-Remote-Sensing-and-GIS

PyRSGIS is a powerful module to read, manipulate and export geo-rasters. The module is built on the GDAL library, and is very efficient for various geospatial analysis.

Please find few example below:

Let's import the module by using the below code
import pyrsgis as rg

We will first start with reading raster and perform some basic operations. Make sure your current working directory is the same where the raster files are located, if not use the following command:
import os
os.chdir("d:/yourDirectory")

Please skip the last two lines of code if the directory is already set to the files location.
Considering that the raster bands are stacked already, reading a multispectral data would be:
yourRaster = rg.readtif("yourFilename.tif")

After the above code, various properties of the raster can be assessed.
print(yourRaster.rows)
print(yourRaster.cols)
will give you the number of rows and columns.

The number of bands can be checked using:
print(yourRaster.nbands)

Any particular band can be extarcted using:
yourBand = yourRaster.getband(bandNumber)

The above code returns the band as array which can be visualised using:
yourRaster.display(yourBand)

The extracted band can be exported using:
yourRaster.export(yourBand, "yourOutputFilename.tif")
This saves the extracted band to the same directory.

If stacked properly, the satellite sensor can also be determined.
print(yourRaster.satellite)

If the above code shows the correct satellite sensor correctly, then this getting this should be easy:
print(yourRaster.bandIndex)
This will show correctly the band number for available bands.

The NDVI (Normalised Difference Vgetaton Index) can be computed easily.
yourndvi = yourRaster.ndvi()
This returns the NDVI array,which can be exported using the same command used for the band above.

yourRaster.export(yourndvi, 'yourNDVI.tif', datatype='float')
One thing to notice is that the NDVI is of float datatype, whereas the raw bands are integer datatype. Float data export uses more space on hard drive, so the default has been set to integer. Therefore, to export any float datatype, the argument should be passed explicitly.

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

pyrsgis-0.1.2.tar.gz (5.6 kB view hashes)

Uploaded Source

Built Distribution

pyrsgis-0.1.2-py3-none-any.whl (7.3 kB view hashes)

Uploaded Python 3

Supported by

AWS AWS Cloud computing and Security Sponsor Datadog Datadog Monitoring Fastly Fastly CDN Google Google Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page