Radiometric calibration S1 – 3rd swath

I’ve coregistered the following products w.r.t. to a single master:

S1A_IW_SLC__1SDV_20150408T181804_20150408T181831_005396_006DB8_64C8.zip
S1A_IW_SLC__1SDV_20150420T181805_20150420T181832_005571_007209_9775.zip
S1A_IW_SLC__1SDV_20150502T181806_20150502T181833_005746_00760A_27F7.zip
S1A_IW_SLC__1SDV_20150514T181806_20150514T181833_005921_007A0B_D78B.zip
S1A_IW_SLC__1SDV_20150526T181807_20150526T181834_006096_007E6D_9349.zip
S1A_IW_SLC__1SDV_20150607T181808_20150607T181835_006271_008373_12E5.zip
S1A_IW_SLC__1SDV_20150607T181808_20150607T181835_006271_008373_467B.zip
S1A_IW_SLC__1SDV_20150619T181808_20150619T181835_006446_008878_C59B.zip
S1A_IW_SLC__1SDV_20150701T181808_20150701T181835_006621_008D4A_5473.zip
S1A_IW_SLC__1SDV_20150713T181809_20150713T181836_006796_009238_C81F.zip
S1A_IW_SLC__1SDV_20150725T181810_20150725T181837_006971_009732_6B92.zip
S1A_IW_SLC__1SDV_20150806T181810_20150806T181837_007146_009C05_3E15.zip
S1A_IW_SLC__1SDV_20150818T181811_20150818T181838_007321_00A0C8_9D9B.zip
S1A_IW_SLC__1SDV_20150830T181811_20150830T181838_007496_00A58F_35B8.zip
S1A_IW_SLC__1SDV_20150911T181812_20150911T181839_007671_00AA55_516A.zip
S1A_IW_SLC__1SDV_20150923T181812_20150923T181839_007846_00AEF0_8058.zip
S1A_IW_SLC__1SDV_20151005T181812_20151005T181839_008021_00B3A8_6A52.zip
S1A_IW_SLC__1SDV_20151017T181812_20151017T181839_008196_00B864_BEE8.zip
S1A_IW_SLC__1SDV_20151029T181812_20151029T181839_008371_00BD18_56A4.zip
S1A_IW_SLC__1SDV_20151110T181812_20151110T181839_008546_00C1B3_8BF0.zip
S1A_IW_SLC__1SDV_20151122T181806_20151122T181833_008721_00C6AA_4C4B.zip
S1A_IW_SLC__1SDV_20151204T181806_20151204T181833_008896_00CB99_1BC3.zip
S1A_IW_SLC__1SDV_20151216T181805_20151216T181832_009071_00D06C_14E3.zip
S1A_IW_SLC__1SDV_20151228T181805_20151228T181832_009246_00D565_0A4A.zip
S1A_IW_SLC__1SDV_20160109T181805_20160109T181832_009421_00DA65_2043.zip

If I then compute the standard deviation of the amplitude (per pixel along the time dimension), it seems the third swath has a radiometric issue. In particular, the amplitudes vary a lot, as compared with swath one and two.

Here’s a plot of the standard deviation. Black means the amplitude varies significantly.

Esteban
SkyGeo

Interesting, can you check whether something changed at one point in time or not for that swath?

Have you used the “stddev()” function from the SNAP toolbox?

I have observed similar issues. I suspect there are some raster count issues in the function.

Here’s a plot of the mean amplitude per image as a function of time for one of the black patches:

There’s a clear slump on 10-Nov-2015.

Hi, I used my own stddev function.

But then that plot shows that the IW3 swath for the 10-Nov-2015 image was corrupt, no?

1 Like

Make sure you avoid using no-values.

Additionally when using own stddev() make sure it has the correct raster count/location, considering there might be data missing it could solve the problem. As for the native stddev() function of the SNAP it has hard time differentiating how many values are at the pixel, it rather takes the amount of all acquisitions in the whole area regardless of no-data-values.

Raster count (iterated over all bands):
if $band <= 0.0005 then 0.0 + $previous_raster_count else 1.0 + $previous_raster_count

stddev function iteration (note that the values will be as good as your average):
if $band <=0.0005 then $previous_stddev else sqrt(( pow($previous_stddev*($previous_raster_count-1), 2)+pow($band-average,2)) / $previous_raster_count )

This one should work properly.

Hi Esteban,

Your problem is unexpected - it is possible that there was an issue when extracting/downloading the product. I suggest you check the size of the 6 measurement tiff files (in the product measurement folder) - they should be the same as given below:

1144762256 Jan 28 13:37 s1a-iw1-slc-vh-20151110t181814-20151110t181839-008546-00c1b3-001.tiff
1144762256 Jan 28 13:37 s1a-iw1-slc-vv-20151110t181814-20151110t181839-008546-00c1b3-004.tiff
1365928760 Jan 28 13:37 s1a-iw2-slc-vh-20151110t181812-20151110t181837-008546-00c1b3-002.tiff
1365928760 Jan 28 13:38 s1a-iw2-slc-vv-20151110t181812-20151110t181837-008546-00c1b3-005.tiff
1324940168 Jan 28 13:38 s1a-iw3-slc-vh-20151110t181813-20151110t181838-008546-00c1b3-003.tiff
1324940168 Jan 28 13:38 s1a-iw3-slc-vv-20151110t181813-20151110t181838-008546-00c1b3-006.tiff

I suspect one of the IW3 tiff files is much smaller than above. If this is the case, I suggest you try to unzip again or download again.

Peter

Hi Peter,

I’ve just checked. I’m getting a bad CRC when unzipping the file. I think I missed this because I’m working with the zipped product.

inflating: S1A_IW_SLC__1SDV_20151110T181812_20151110T181839_008546_00C1B3_8BF0.SAFE/measurement/s1a-iw3-slc-vv-20151110t181813-20151110t181838-008546-00c1b3-006.tiff   bad CRC 617be33e  (should be db7719e8)

I think this problem should get flagged by the SNAP product reader.

Thanks for your help! Good catch. I will re-download.

Regards,

Esteban
SkyGeo

1 Like

Could you try again with the new update 2.0.2 The reader has been updated to try to detect this. Also, the product library scan has an option to check for bad CRCs in Zips.