Sentinel 2 mosaic

[Olivier]

Hi Javier,
Do you still have the same problem? Was it fixed since then?
Is it also the case with version 4.0?
Please let us know, I inform Julien, Nicolas and Omar about this behaviour, but they surely are already.
Thanks a lot for these elements of observation.
Cheers
Olivier

[Tanemahuta]

Hi Olivier,

There are still problems with the mosaic tool (and re-projection too).
When trying to select a pre-defined CRS (ESPG:2193), you can find it OK but it is not put into the white space.

Cheers,
Mike

[plitkey]

Hello,
Now that the products are in several files, they need to be combined…
I’m trying to mosaic resampled products in .dim files and data folders. I’m curious about the pixel size - how come the default is not the pixel size of the resampled product? The fill in box in the GUI is in degrees and the smallest it accepts if 0.001. My mosaic is very coarse at best and at default value I get huge pixels. Question is: how to get 10m pixels when mosaicing BEAM-DIMAP products that are resampled to 10m?

-Paula

[marpet]

In the Mosaic is WGS_84 as default projection selected. Its CRS is defined in degree. That’s why the pixel size is given in degree by default.
What’s wrong here is that you can’t define a small pixel size correctly. Below 0.001 it only shows 0.0. But it uses the values you have entered. So it is only a display issue. However this is wrong.

If you change in the upper part the projection, e.g. to UTM, you can enter the pixel size in meter.

[plitkey]

Thanks!

It went well untill I clicked run and got an error message about my cloud and cloud shadow condition, when I tried to remove the condition the whole mosaic app froze and I could not select anything. I restarted SNAP and started new mosaic, but the tabs Map Project Definition and Variables & Conditions remained frozen. Maybe I need to restart unity? Or reboot?

Fix: It seems I did not select the source, it was not enough to + it, I needed to select it also… now it seems to work

-Paula

[marpet]

OK, good that it works now.

[plitkey]

Well. It wrote a 104 GB mosaic from two granules in 1h 43 min. The mosaic looks universally grey in SNAP. It seems that something went wrong.

-Paula

[ABraun]

maybe there was a mismatch between the coordinate system of the dataset and the boundary coordinates.
I often see extemely huge files being created when UTM coordinates are applied as boundaries for WGS84 projected files.

Just an idea…

[plitkey]

I didn’t give boundary coordinates, maybe that was the problem. I managed to run the mosaic for a small area in the border of the granules just fine. I have subsetted the same area from older products and then the .data folder size was ~266 MB, with this mosaicked product it is ~683 MB. I don’t understand why there is double bands in the band listing … what is in the _count bands?

[marpet]

I have recently explained what’s in the count bands:

L5 & Sentinel2 - Change Detection...many steps to the finish line

The '*_count' bands are created by the mosaic step. They count how many source pixels contributed to one target pixel. This is often useful when mosaicking multiple products. They are also used in the valid expression of the B# bands. You can remove the count bands but you also have to change the valid expression of the B# bands. I agree that it would be a good option for the Mosaic to disable the generation of the count bands. I've put it into our issue tracker. To ease the workflow a bit …

[mikeV]

Hi,
I have the same problem to create a mosaic. I get the error of Product 'sourceProduct2' contains bands of different sizes. Currently it is not possible to use it for mosaicking.
I tried doing the mosaicing operation with the SNAP tutorial, but I cannot achieve to do it.

But when I do the mosaicing operation with one band (for example band 2 from each image) I have not any error and I will get the mosaic.

[Ardo]

If I don’t get the multiple reader dialog when opening SENTINEL-2 data . How do I enable it?
Need it as preparation for mosaicing, getting the same resolution for all bands as indicated above.
Running W7 and W8.
Regards
/Jonas

[obarrilero]

Hi,

The 10/20/60 m resolution readers are no longer available, but you will get the same result if you use the resampling operator after opening the multi-resolution product.

About the 10m, 20m and 60m readers removed in S2TBX 4.0

Hi, I consider a pity the fact of having removed 10m, 20m and 60m readers removed in S2TBX 4.0. In one side, it is a good idea to make the user think about the interpolation algorithm for up or downscaling. But on the other side, knowing that the removed readers used the nearest neighbor interpolation algorithm, it was very useful while it was possible to let images as native. Javier

[pettehadi]

I have problem with mosaicing the sentinel 2 images. I tried resampling and subseting but it is not fixed. Do you find a way to fix mosaicing problem anybody?

[Poanh]

I have a problem with mosaic 2 sentinel images . I subset them then make a mosaic. However, it said can not display source product, …" as below. what is a matter?

[marpet]

This is a known issue and it is actually already solved. The fix will be released soon.
Meanwhile you can try to set the mosaic boundary manually. This should help

See also:

Error when trying to mosaic Sentinel 2 - weird result stretched over the globe

Hello, I'm trying to mosaic some adjacent Sentinel 2 tiles together over an area of Russia. I've resampled the tiles to 10 metres and put them into the mosaicing operation where I've chosen bands 2, 3 and 4 and set it to use "Auto UTM" as shown below. [image] [image] [image] The result I get seems weird, where the mosaic is stretched over a huge part of northern Russia and is clearly incorrect. [image] For reference, the actual distribution of the tiles are as below: [image] An…

[falahfakhri]

Please take a look at this video it is already released now

[ThomasLankester]

Hi, I followed the video to investigate the blending operation at native resolution for 1 band (8a) at native 20m resolution (tab 2 window below and final results) as the final result screenshot shows, the processing seems to have picked up some aspect of the cloud (mask?) with the rest of the pixel values set to NaN. I am using SNAP v5.0.8 on a 64-bit kubuntu OS.
If I do the same processing but not at 30m resolution instead of the native (20m) resolution then I do get a mosaic but the ‘blend’ option is ignored.

RE the video:
-it does not state the nature of the blend algorithm, is it distance weighted?
-why does the ‘expression’ field in the third control tab seems to have no function, unlike the one in the basic ‘Mosaic’ tool?

[kraftek]

Hello Thomas.

To answer to your questions:

• the blending is just a simple pairwise alpha blending (on two images, the first one has weight 1, the second one 0.5);
• the ‘expression’, in the case of multi-size products, has not too much sense. Expressions were designed to work with pixel values. In the case of, let’s say, B8 and B8a, since a pixel in B8a corresponds to 4 pixels in B8, it doesn’t make sense to use these values in an expression (which value from B8 will it use?);
  If you do the mosaic at 30m, that means you resample the bands at 30m and the resulting mosaic is no more a multi-size one, but a single size. AFAIK, the blending, in this case, doesn’t work if there are NaN values (it’s a JAI - the image processing library used - limitation).
  For the result you’ve obtained for multi-size, could you please specify the products used?

Cosmin.

[ThomasLankester]

Hi Cosmin,

I used S2 tiles S2A_MSIL1C_20170618T105621_N205_R094_T30SVF_20170618T110415 and S2A_MSIL1C_20170628T105651_N205_R094_T30SVF_20170618T110445.

With the blending, if the same blend is applied across the overlap zone, will this not lead to lines (radiometric jumps) along both overlap edges?

For the parameter expression I thought this would allow filtering using flag or classification information to be used (as with band maths). Otherwise, now are bad (invalid, snow, cloud etc.) pixels to be removed?

Tom