I would like to bring back a topic that has been addressed a while ago in this forum but I find that no clear answer is available yet. I have a dataset of S1 and S2 images:
For S2, I am re-sampling and subsetting to my ROI
For S1, I am pre-processing the IW GRD products in SNAP using the following chain: apply orbit file, calibration, speckle filter, terrain correction (including UTM projection) and subset (to the same ROI as before).
My doubt is about stacking both datasets after that. Based on some previous discussion on this forum, I have seen that some people have tried the corregistration tool available for SAR.
But based on those replies:
However, based on this post
It seems in that case there was an offset of around 40 m when doing the terrain correction to S1 images and then stacking with Sentinel-2
So, what should be the procedure or at least the one leading to a higher accuracy? How can we check the accuracy (visualizing the overlapping images, adding some control points with known coordinates, etc.)
Once S1 images are terrain corrected, re-projected in UTM (as S2) and stacked together with S2, I guess the geo-location accuracy of overlapping S1/S2 pixels in the stack is defined by the orthorectification algorithm.
Which level/order of accuracy can we expect in terms of pixel alignment? (sub-pixel?)
Achieving centimeter accuracy presumes that one perform additional corrections for finer perturbations like atmospheric path delay, frame shift, and solid Earth tides. AFAIK, these are not included in the SNAP s1tbx processing.
An alternative way to snap Sentinel-1 GRDH pixel grid to Sentinel-2 pixel grid is to set the ‘standardGridOrigin’ of the output Sentinel-1 product when applying the Range Doppler Terrain Correction.
This could be done by setting <standardGridOriginX> and <standardGridOriginY> parameters in the Range Doppler Terrain Correction operator.
I believe this solution could reduce as less as possible the introduction of artifacts. Later, a stacking operator should be used in order to stack together (multitemporal) Sentinel-1 and Sentinel-2 data.
I have shared a GPF graph that is preprocessing Sentinel-1 GRDH data taking into account this solution to snap output images to Sentinel-2 grids: