https://github.com/nunorc/mysnap

utilities for quickly importing and handling ESA SNAP toolbox products
https://github.com/nunorc/mysnap

esa experimental snap underdevelopment

Last synced: about 6 hours ago
JSON representation

utilities for quickly importing and handling ESA SNAP toolbox products

• Host: GitHub
• URL: https://github.com/nunorc/mysnap
• Owner: nunorc
• License: mit
• Created: 2021-11-02T14:54:44.000Z (about 3 years ago)
• Default Branch: master
• Last Pushed: 2021-11-03T10:40:18.000Z (about 3 years ago)
• Last Synced: 2024-01-27T13:40:27.644Z (9 months ago)
• Topics: esa, experimental, snap, underdevelopment
• Language: Python
• Homepage: https://nunorc.github.io/mysnap/
• Size: 2.72 MB
• Stars: 0
• Watchers: 1
• Forks: 0
• Open Issues: 0
• Metadata Files:
  • Readme: README.rst
  • License: LICENSE.txt

Awesome Lists containing this project

README

        

mysnap

===========================

An experimental package with simple utilities for quickly importing

and handling ESA SNAP toolbox products in BEAM-DIMAP format.

Package available from `GitHub `_,

documentation available `here `_.

Quick Start

---------------------------

Install package from the git repository:

.. code-block:: bash

    $ pip install git+https://github.com/nunorc/mysnap@master

Open a product, build images by stacking bands, and split in tiles.

.. code-block:: python

    import mysnap

    # product filename in beam-dimap format

    filename = '~/Products/subset_2021_resampled.dim'

    # load product from file

    product = mysnap.Product(filename)

    # check rasterized size of image and number of bands

    product.meta.ncols, product.meta.nrows, product.meta.nbands  # 2312, 2312, 43

    # get list of bands

    product.bands  # ['b1', 'b2', 'b3', ...

    # access individual band

    product.b1.name    # 'b1'

    product.b1.desc    # 'B1 (443.0)'

    product.b1.array   # array([[1731., 1731., 1731., ...

    product.b1.array.shape   # (2312, 2312)

    # build a stacked array of selected bands

    bands = ['b2', 'b3', 'b4', 'b5', 'b6', 'b7', 'b8', 'b8a', 'b11', 'b12']

    arr = product.stack(bands, normalize=mysnap.normalize)

    arr.shape   # (2312, 2312, 10)

    # split stack in 64x64 tiles

    tiles = mysnap.tiles(arr, size=64)

    tiles.shape   # (1296, 64, 64, 10)

    # build a rgb image

    bands = ['b4', 'b3', 'b2']

    rgb = product.stack(bands, quantiles=[0.025, 0.975], max=3000, normalize=mysnap.normalize)

Acknowledgments

---------------------------

Raster data imported with `rasterio `_.

Thank you to `ESA `_,

the authors of the `SNAP Toolbox `_,

and upstream packages and products.