I'm running some gpt commands (using XML template) to do some band math over L2 at 250 meter resolution, but processing it takes too much to export data to GeoTiff, as you can see in the code below I'm also doing a reprojection when executing the commands. For example it takes about 40 min to execute
I have another gpt XML template that I have only exports and reprojects a single band, but still processing time it's very slow.
Any recommendations, perhaps using GDAL or any other external tool, or approach to solve this issue.
<graph id="fai_processing">
<version>1.0</version>
<node id="reprojectNode">
<operator>Reproject</operator>
<sources>
<source>${source}</source>
</sources>
<parameters>
<crs>
GEOGCS["WGS84(DD)",
DATUM["WGS84",
SPHEROID["WGS84", 6378137.0, 298.257223563]],
PRIMEM["Greenwich", 0.0],
UNIT["degree", 0.017453292519943295],
AXIS["Geodetic longitude", EAST],
AXIS["Geodetic latitude", NORTH]]
</crs>
<resampling>Nearest</resampling>
<orientation>0.0</orientation>
<pixelSizeX>0.002250022500225</pixelSizeX>
<pixelSizeY>0.002250022500225</pixelSizeY>
<orthorectify>false</orthorectify>
<noDataValue>NaN</noDataValue>
<includeTiePointGrids>true</includeTiePointGrids>
<addDeltaBands>false</addDeltaBands>
<maskExpression>!l2_flags.LAND and !l2_flags.HILT and !l2_flags.CLDICE and !l2_flags.STRAYLIGHT and !l2_flags.CHLFAIL</maskExpression>
<transferValidPixelExpression>true</transferValidPixelExpression>
<applyValidPixelExpression>true</applyValidPixelExpression>
</parameters>
</node>
<node id="bandMathNode">
<operator>BandMaths</operator>
<sources>
<source>reprojectNode</source>
</sources>
<parameters>
<targetBands>
<targetBand>
<name>fai</name>
<expression>rhos_859 - (rhos_645 + (rhos_1240 - rhos_645) * ((859 - 645) / (1240 - 645)))</expression>
<description>Floating Algae Index</description>
<type>float32</type>
<noDataValue>-32767.0</noDataValue>
</targetBand>
</targetBands>
</parameters>
</node>
</graph>
Thanks again!