Binning

In [1]:

from snapista.binning import aggregators
from snapista.binning import BinningVariable, BinningOutputBands, BinningBand, Aggregators, BinningVariables
from snapista import Graph
from snapista import Operator
import lxml.etree as etree
import os

from snapista.binning import aggregators from snapista.binning import BinningVariable, BinningOutputBands, BinningBand, Aggregators, BinningVariables from snapista import Graph from snapista import Operator import lxml.etree as etree import os

In [2]:

os.environ['PATH'] = '/srv/conda/envs/env_snap/snap/bin:' +  os.environ['PATH']

os.environ['PATH'] = '/srv/conda/envs/env_snap/snap/bin:' + os.environ['PATH']

Create an Aggregator¶

In [3]:

ag = aggregators.AggregatorAvg(varName='OTCI', targetName="otci", outputCounts='true', outputSums='true', weightCoeff=1.0)

ag = aggregators.AggregatorAvg(varName='OTCI', targetName="otci", outputCounts='true', outputSums='true', weightCoeff=1.0)

Create a binning operator¶

In [4]:

binning = Operator('Binning')

binning.sourceProductFormat = 'sen3'

binning.aggregatorConfigs = Aggregators([ag,
                                        aggregators.AggregatorAvgOutlier(varName='OTCI', targetName='outlier'),
                                        aggregators.AggregatorMinMax(varName='OTCI', targetName='otci'),
                                        aggregators.AggregatorOnMaxSet(targetName='omvn', onMaxVarName='OTCI', setVarNames='OTCI'),
                                        aggregators.AggregatorPercentile(varName='OTCI', percentage=90, targetName='percent'),
                                        aggregators.AggregatorSum(varName='OTCI', targetName='sum')])

binning.sourceProductPaths = '/data/S3*/*.xml'

binning.outputFile = '/home/vscode/binned_new.dim'

binning = Operator('Binning') binning.sourceProductFormat = 'sen3' binning.aggregatorConfigs = Aggregators([ag, aggregators.AggregatorAvgOutlier(varName='OTCI', targetName='outlier'), aggregators.AggregatorMinMax(varName='OTCI', targetName='otci'), aggregators.AggregatorOnMaxSet(targetName='omvn', onMaxVarName='OTCI', setVarNames='OTCI'), aggregators.AggregatorPercentile(varName='OTCI', percentage=90, targetName='percent'), aggregators.AggregatorSum(varName='OTCI', targetName='sum')]) binning.sourceProductPaths = '/data/S3*/*.xml' binning.outputFile = '/home/vscode/binned_new.dim'

In [5]:

binning

binning

Out[5]:

Operator('Binning', sourceProductPaths='/data/S3*/*.xml', sourceProductFormat='sen3', sourceGraphPaths='None', region='None', startDateTime='None', periodDuration='None', timeFilterMethod='NONE', minDataHour='None', numRows='2160', superSampling='1', maxDistanceOnEarth='-1', maskExpr='None', variableConfigs='None', aggregatorConfigs='<snapista.binning.aggregators.Aggregators object at 0x7f1da3cb10d0>', postProcessorConfig='None', outputType='Product', outputFile='/home/vscode/binned_new.dim', outputFormat='BEAM-DIMAP', bandConfigurations='None', productCustomizerConfig='None', outputBinnedData='false', outputTargetProduct='true', metadataPropertiesFile='./metadata.properties', metadataTemplateDir='.', metadataAggregatorName='NAME', planetaryGridClass='org.esa.snap.binning.support.SEAGrid')

Create the Graph¶

In [6]:

g = Graph()

g.add_node(operator=binning, 
           node_id='binning_1')

g = Graph() g.add_node(operator=binning, node_id='binning_1')

aggregatorConfigs <snapista.binning.aggregators.Aggregators object at 0x7f1da3cb10d0>
bandConfigurations None
postProcessorConfig None
productCustomizerConfig None
variableConfigs None

In [7]:

g.add_node(operator=binning, 
           node_id='binning_1')

g.add_node(operator=binning, node_id='binning_1')

In [8]:

g.view()

g.view()

<graph>
  <version>1.0</version>
  <node id="binning_1">
    <operator>Binning</operator>
    <sources/>
    <parameters class="com.bc.ceres.binding.dom.XppDomElement">
      <aggregators>
        <aggregator>
          <type>AVG</type>
          <varName>OTCI</varName>
          <targetName>otci</targetName>
          <weightCoeff>1.0</weightCoeff>
          <outputCounts>true</outputCounts>
          <outputSums>true</outputSums>
        </aggregator>
        <aggregator>
          <type>AVG_OUTLIER</type>
          <varName>OTCI</varName>
          <targetName>outlier</targetName>
        </aggregator>
        <aggregator>
          <type>MIN_MAX</type>
          <varName>OTCI</varName>
          <targetName>otci</targetName>
        </aggregator>
        <aggregator>
          <type>ON_MAX_SET</type>
          <targetName>omvn</targetName>
          <onMaxVarName>OTCI</onMaxVarName>
          <setVarNames>OTCI</setVarNames>
        </aggregator>
        <aggregator>
          <type>PERCENTILE</type>
          <varName>OTCI</varName>
          <targetName>percent</targetName>
          <percentage>90</percentage>
        </aggregator>
        <aggregator>
          <type>SUM</type>
          <varName>OTCI</varName>
          <targetName>sum</targetName>
        </aggregator>
      </aggregators>
      <maskExpr/>
      <maxDistanceOnEarth>-1</maxDistanceOnEarth>
      <metadataAggregatorName>NAME</metadataAggregatorName>
      <metadataPropertiesFile>./metadata.properties</metadataPropertiesFile>
      <metadataTemplateDir>.</metadataTemplateDir>
      <minDataHour/>
      <numRows>2160</numRows>
      <outputBinnedData>false</outputBinnedData>
      <outputFile>/home/vscode/binned_new.dim</outputFile>
      <outputFormat>BEAM-DIMAP</outputFormat>
      <outputTargetProduct>true</outputTargetProduct>
      <outputType>Product</outputType>
      <periodDuration/>
      <planetaryGridClass>org.esa.snap.binning.support.SEAGrid</planetaryGridClass>
      <region/>
      <sourceGraphPaths/>
      <sourceProductFormat>sen3</sourceProductFormat>
      <sourceProductPaths>/data/S3*/*.xml</sourceProductPaths>
      <startDateTime/>
      <superSampling>1</superSampling>
      <timeFilterMethod>NONE</timeFilterMethod>
    </parameters>
  </node>
</graph>

In [9]:

g.save_graph(filename='s3.xml')

#g.run()

g.save_graph(filename='s3.xml') #g.run()

Processing the graph
Executing processing graph
INFO: org.esa.snap.core.gpf.operators.tooladapter.ToolAdapterIO: Initializing external tool adapters
45%90% done.
INFO: org.esa.s2tbx.dataio.gdal.GDALVersion: GDAL not found on system. Internal GDAL 3.0.0 from distribution will be used. (f1)
INFO: org.esa.s2tbx.dataio.gdal.GDALVersion: Internal GDAL 3.0.0 set to be used by SNAP.
INFO: org.esa.snap.core.util.EngineVersionCheckActivator: Please check regularly for new updates for the best SNAP experience.
INFO: org.esa.s2tbx.dataio.gdal.GDALVersion: Internal GDAL 3.0.0 set to be used by SNAP.

Out[9]:

0

In [ ]:

In [ ]: