S1 preprocessing Terrain Correction map projection

Hi,
I am trying to process S1 data using gpt. I used gpt for the read, apply orbit File, thermal noise removal, calibration. However when I use it for terrain correction I receive results with wrong projection. I define projection EPSG:32635, but Snap shows WGS84.
The MetaData looks

The used command are as follow:
gpt Terrain-Correction -Ssource=‘[D:/pysnap/04_calib/target.dim](file:///D:/pysnap/04_calib/target.dim)’ -PsourceBands=‘Sigma0_VV,Sigma0_VH’ -PpixelSpacingInMeter=10.0 -PmapProjection=‘GEOGCS[“WGS 84”,DATUM[“WGS_1984”,SPHEROID[“WGS 84”,6378137,298.257223563,AUTHORITY[“EPSG”,“7030”]],AUTHORITY[“EPSG”,“6326”]],PRIMEM[“Greenwich”,0,AUTHORITY[“EPSG”,“8901”]],UNIT[“degree”,0.0174532925199433,AUTHORITY[“EPSG”,“9122”]],AXIS[“Latitude”,NORTH],AXIS[“Longitude”,EAST],AUTHORITY[“EPSG”,“4326”]]’ -PalignToStandardGrid=false -PapplyRadiometricNormalization=false -PauxFile=‘Latest Auxiliary File’ -PdemName=‘SRTM 3Sec’ -PdemResamplingMethod=‘BILINEAR_INTERPOLATION’ -PexternalDEMApplyEGM=true -PexternalDEMNoDataValue=0 -PimgResamplingMethod=‘BILINEAR_INTERPOLATION’ -PincidenceAngleForSigma0=‘Use projected local incidence angle from DEM’ -PnodataValueAtSea=true -PoutputComplex=false -PsaveBetaNought=false -PsaveDEM=false -PsaveGammaNought=false -PsaveIncidenceAngleFromEllipsoid=false -PsaveLatLon=false -PsaveLocalIncidenceAngle=false -PsaveProjectedLocalIncidenceAngle=false -PsaveSelectedSourceBand=true -PsaveSigmaNought=false -PstandardGridOriginX=0 -PstandardGridOriginY=0

and the picture of results are as follow:

There are bug in map_projection, because it must be WGS 84 / UTM zone 34N.

In the command you used it looks like it’s EPSG:4326, not the one you want (EPSG:32634). Edit the command to match this info, WGS 84 / UTM zone 34N - EPSG:32634

1 Like

Thanks, it works!
Now I have UTM-34N projection.
But in next step I have problem with subset.
For subset I define polygon using coordinates in UTM-34N

gpt Subset -Ssource=‘D:/pysnap/06_terCor/target.dim’ -PgeoRegion=‘POLYGON ((699547.87111 6286886.878148, 700008.216467 6287234.832451, 700071.525254 6287021.212184, 699606.912751 6286688.063824, 699547.87111 6286886.878148))’ -PsourceBands=‘Sigma0_VV,Sigma0_VH’ -PcopyMetadata=false -PfullSwath=false -PsubSamplingX=1 -PsubSamplingY=1
INFO: org.esa.snap.core.gpf.operators.tooladapter.ToolAdapterIO: Initializing external tool adapters
INFO: org.esa.snap.core.util.EngineVersionCheckActivator: Please check regularly for new updates for the best SNAP experience.
INFO: org.hsqldb.persist.Logger: dataFileCache open start
WARNING: org.esa.snap.core.gpf.common.SubsetOp: No intersection with source product boundary target
Executing operator…
20%100% done.
INFO: org.esa.snap.core.gpf.common.WriteOp: Start writing product Empty_target to target.dim
Writing…
20%100% done.
INFO: org.esa.snap.core.gpf.common.WriteOp: End writing product Empty_target to target.dim
INFO: org.esa.snap.core.gpf.common.WriteOp: Time: 0.004 s total, Infinity ms per line, Infinity ms per pixel

I found that coordinates must be defined in WGS84 (DD) (EPSG:4326)

1 Like