GPT creating 75 GB temporary file to calibrate and ellipsoid project data

Hi, I am using the GPT to radiometrically calibrate and geometrically correct a S1 scene by using the Ellipsoid Range Doppler Correction. My graph xml can be found at the end. My problem is it is taking about an hour to correct one scene but when i do it in the GUI it takes about 15 minutes. In the process of using the GPT it creates a massive temporary file which I believe is slowing things down.

Am I doing something wrong? I am passing the file, the specific polarizations i want and the output path as variables in a python script. I know this shouldn’t be taking this long and I’ve run out of ideas…

> <graph id="Graph">
>   <version>1.0</version>
>   <node id="Read">
>     <operator>Read</operator>
>     <sources/>
>     <parameters class="com.bc.ceres.binding.dom.XppDomElement">
>       <file>${inputGRD}</file>
>     </parameters>
>   </node>
>   <node id="Calibration">
>     <operator>Calibration</operator>
>     <sources>
>       <sourceProduct refid="Read"/>
>     </sources>
>     <parameters class="com.bc.ceres.binding.dom.XppDomElement">
>       <sourceBands/>
>       <auxFile>Product Auxiliary File</auxFile>
>       <externalAuxFile/>
>       <outputImageInComplex>false</outputImageInComplex>
>       <outputImageScaleInDb>false</outputImageScaleInDb>
>       <createGammaBand>false</createGammaBand>
>       <createBetaBand>false</createBetaBand>
>       <selectedPolarisations>${polarization}</selectedPolarisations>
>       <outputSigmaBand>true</outputSigmaBand>
>       <outputGammaBand>false</outputGammaBand>
>       <outputBetaBand>false</outputBetaBand>
>     </parameters>
>   </node>
>   <node id="Ellipsoid-Correction-RD">
>     <operator>Ellipsoid-Correction-RD</operator>
>     <sources>
>       <sourceProduct refid="Calibration"/>
>     </sources>
>     <parameters class="com.bc.ceres.binding.dom.XppDomElement">
>       <sourceBands>Sigma0_${polarization}</sourceBands>
>       <demName>SRTM 3Sec</demName>
>       <externalDEMFile/>
>       <externalDEMNoDataValue>0.0</externalDEMNoDataValue>
>       <externalDEMApplyEGM>true</externalDEMApplyEGM>
>       <demResamplingMethod>BILINEAR_INTERPOLATION</demResamplingMethod>
>       <imgResamplingMethod>BILINEAR_INTERPOLATION</imgResamplingMethod>
>       <pixelSpacingInMeter>10.0</pixelSpacingInMeter>
>       <pixelSpacingInDegree>8.983152841195215E-5</pixelSpacingInDegree>
>       <mapProjection>PROJCS[&quot;WGS 84 / Auto UTM&quot;, &#xd;
>   GEOGCS[&quot;WGS84(DD)&quot;, &#xd;
>     DATUM[&quot;WGS84&quot;, &#xd;
>       SPHEROID[&quot;WGS84&quot;, 6378137.0, 298.257223563]], &#xd;
>     PRIMEM[&quot;Greenwich&quot;, 0.0], &#xd;
>     UNIT[&quot;degree&quot;, 0.017453292519943295], &#xd;
>     AXIS[&quot;Geodetic longitude&quot;, EAST], &#xd;
>     AXIS[&quot;Geodetic latitude&quot;, NORTH]], &#xd;
>   PROJECTION[&quot;Transverse_Mercator&quot;], &#xd;
>   PARAMETER[&quot;central_meridian&quot;, 3.0], &#xd;
>   PARAMETER[&quot;latitude_of_origin&quot;, 0.0], &#xd;
>   PARAMETER[&quot;scale_factor&quot;, 0.9996], &#xd;
>   PARAMETER[&quot;false_easting&quot;, 500000.0], &#xd;
>   PARAMETER[&quot;false_northing&quot;, 0.0], &#xd;
>   UNIT[&quot;m&quot;, 1.0], &#xd;
>   AXIS[&quot;Easting&quot;, EAST], &#xd;
>   AXIS[&quot;Northing&quot;, NORTH]]</mapProjection>
>       <alignToStandardGrid>false</alignToStandardGrid>
>       <standardGridOriginX>0.0</standardGridOriginX>
>       <standardGridOriginY>0.0</standardGridOriginY>
>       <nodataValueAtSea>true</nodataValueAtSea>
>       <saveDEM>false</saveDEM>
>       <saveLatLon>false</saveLatLon>
>       <saveIncidenceAngleFromEllipsoid>false</saveIncidenceAngleFromEllipsoid>
>       <saveLocalIncidenceAngle>false</saveLocalIncidenceAngle>
>       <saveProjectedLocalIncidenceAngle>false</saveProjectedLocalIncidenceAngle>
>       <saveSelectedSourceBand>true</saveSelectedSourceBand>
>       <outputComplex>false</outputComplex>
>       <applyRadiometricNormalization>false</applyRadiometricNormalization>
>       <saveSigmaNought>true</saveSigmaNought>
>       <saveGammaNought>false</saveGammaNought>
>       <saveBetaNought>false</saveBetaNought>
>       <incidenceAngleForSigma0>Use projected local incidence angle from DEM</incidenceAngleForSigma0>
>       <incidenceAngleForGamma0>Use projected local incidence angle from DEM</incidenceAngleForGamma0>
>       <auxFile>Latest Auxiliary File</auxFile>
>       <externalAuxFile/>
>     </parameters>
>   </node>
>   <node id="Write">
>     <operator>Write</operator>
>     <sources>
>       <sourceProduct refid="Ellipsoid-Correction-RD"/>
>     </sources>
>     <parameters class="com.bc.ceres.binding.dom.XppDomElement">
>       <file>${output}</file>
>       <formatName>BEAM-DIMAP</formatName>
>     </parameters>
>   </node>
>   <applicationData id="Presentation">
>     <Description/>
>     <node id="Read">
>             <displayPosition x="37.0" y="134.0"/>
>     </node>
>     <node id="Calibration">
>       <displayPosition x="230.0" y="278.0"/>
>     </node>
>     <node id="Ellipsoid-Correction-RD">
>       <displayPosition x="319.0" y="241.0"/>
>     </node>
>     <node id="Write">
>             <displayPosition x="455.0" y="135.0"/>
>     </node>
>   </applicationData>
> </graph>

S1GRD_Cal.xml (4.3 KB)

Hi Zee,

Could you please provide the product ID of the Sentinel-1 scene so I can test on my end?

Cheers!
-Alex

Hi Alex,

I’ve got a product here with filename S1A_EW_GRDM_1SDH_20150624T134029_20150624T134129_006516_008A78_1656.zip

The link to download it is here: https://scihub.copernicus.eu/dhus/odata/v1/Products(‘61ec814b-30cd-460c-9f54-f0c83e7de033’)/$value

Or here from the ASF: https://datapool.asf.alaska.edu/GRD_MD/SA/S1A_EW_GRDM_1SDH_20150624T134029_20150624T134129_006516_008A78_1656.zip

Thanks in advance!
Zee