POF-ML

Examples

In the following sections, you will be able to look at some of the examples of POF-ML we have made over the years

Deforestation in Cameroon Kalissaye Ramsar site Overlay RCSE of Sentinel-1 offshore scenes Sentinel-2 occurrence map Landsat7 pan-sharpening Sentinel-2 and Landsat-8 pan-sharpening Active fires and burnt areas detection DEM isocurve Land / water mask on Lac Naivasha (Kenya) Sentinel-1 Land/Water RCSE MERIS Tristimulus composite (BEAM VISAT) Parc Marin des Mangroves (mv,mh,cv) Sentinel-1 Landsat Collection 2 Level 2 - Surface Temperature Sentinel-1 - Flooding Valencia

Sentinel Vision - SED-734 - Deforestation in Cameroon

25.09.2020 by Zhour Najoui

Look at the result in Sentinel Vision by clicking on this link.

Product list

Sentinel2A/MSI - S2A_MSIL1C_20151209T093402_N0204_R136_T32NQM_20151209T093921
Sentinel2A/MSI - S2A_MSIL1C_20151209T093402_N0204_R136_T32NPM_20151209T093921
Sentinel2A/MSI - S2A_MSIL1C_20191208T093401_N0208_R136_T32NQM_20191208T104325
Sentinel2A/MSI - S2A_MSIL1C_20191208T093401_N0208_R136_T32NPM_20191208T104325

POF-ML script

lut(
	mask(
		mask(
			dif(
				ndi(mean(2:QT_B08,3:QT_B08),mean(2:QT_B04,3:QT_B04),0,0.1),
				ndi(mean(0:QT_B08,1:QT_B08),mean(0:QT_B04,1:QT_B04),0,0.6)
			),
			neg(ndi(mean(0:QT_B11,1:QT_B11),mean(0:QT_B08,1:QT_B08),-0.9,-0.01)),
			80
		),
		neg(ndi(mean(2:QT_B11,3:QT_B11),mean(2:QT_B08,3:QT_B08),-0.3,-0.5)),
		1
	),
	reds,-60,60,0,0,-,10
)

Sentinel Vision - SED-740 - Kalissaye Ramsar site

02.10.2020 by Zhour Najoui

Look at the result in Sentinel Vision by clicking on this link.

Product list

Sentinel2B/MSI - S2B_MSIL1C_20191129T113329_N0208_R080_T28PCV_20191201T132309

POF-ML script

Layer “Bare_sols - MNDWI”

lut(
    mask(ndi(QT_B03,QT_B11,-0.5,0.5),neg(ndi(QT_B03,QT_B11,-0.5,0.5)),230),
    lng_f7f71e_f7f71e,0,255,0
)

Layer “Water - MNDWI”

lut(
    mask(ndi(QT_B03,QT_B11,-0.5,0.5),ndi(QT_B03,QT_B11,-0.5,0.6),255),
    lng_f7866f_40ff00_039c0b_b4edbe_dbffe2_aafaf7_1b0aff,0,255,0,0,-,10
)

Layer “Vegetation - NDWI-SM”

lut(
    mask(ndi(QT_B08,QT_B11,-0.1,0.8),neg(ndi(QT_B08,QT_B11,0.1,0.5)),180),
    lng_276b05_0ffa98,0,255,0,0,-,10
)

Layer “Vegetation2 - NDVI”

lut(
    mask(ndi(QT_B04,QT_B08,-0,-0.7),ndi(QT_B04,QT_B08,-0,-0.7),10),
    lng_282e2d_26ff00,0,255,0,0,-,10
)

Layer “Built-up_areas + bare_sols - dif(NDBI,NDVI)”

lut(
    mask(ndi(QT_B11,QT_B08,0.1,0.8),ndi(QT_B11,QT_B08,-0.05,0.4),40),
    lng_ff2f00/50_ff2f00,0,255,0,0,-,10
)

VtWeb - Overlay RCSE of Sentinel-1 scenes offshore

31.08.2021 by Serge RIAZANOFF

Look at the result in VtWeb by clicking on this link.

Product list

Sentinel1B/CSAR/IW - S1B_IW_GRDH_1SDV_20210821T001458_20210821T001523_028337_03618E_64D4
Sentinel1B/CSAR/IW - S1B_IW_GRDH_1SDV_20210821T001523_20210821T001548_028337_03618E_C2F5

POF-ML script

over(
    sar(
        stc(
            pp(s1nc(band(0:QT_vv))),
            1:RCS_vv,
            140,
            60,
            500
        ),
        0:rasterCrsOrtho,
        null,
        null
    ),
    sar(
        stc(
            pp(s1nc(band(1:QT_vv))),
            1:RCS_vv,
            140,
            60,
            500
        ),
        1:rasterCrsOrtho,
        null,
        null
    )
)

Sentinel Vision - SED-944 Sentinel-2 occurrence map

05.10.2021 by Kévin GROSS

Look at the result in Sentinel Vision by clicking on this link.

Product list

OCC/Sentinel-2 - S2A_OCCURRENCE_MAP_2020
OCC/Sentinel-2 - S2B_OCCURRENCE_MAP_2020

POF-ML script

lut(
    for(QT_occ,sum),
    jet2,
    1,
    2000,
    2
)

lut(
    for(QT_occ,sum),
    jet2,
    1,
    2000,
    log
)

Sentinel Vision - SED-942 - Landsat-7 pan-sharpening

04.10.2021 by Kevin Gross

Look at the result in Sentinel Vision by clicking on this link.

Product list

LANDSAT7/ETM - LE71980302001264NSG00

POF-ML script

pan(
    geo(
        sta(band(QT_ETM7),HT_ETM7,0.02,0.1),
        RasterCRS_ETM1,
        bl,
        0.5,
        0.5
    ),
    geo(
        sta(band(QT_ETM5),HT_ETM5,0.02,0.1),
        RasterCRS_ETM1,
        bl,
        0.5,
        0.5
    ),
    geo(
        sta(band(QT_ETM3),HT_ETM3,0.02,0.1),
        RasterCRS_ETM1,
        bl,
        0.5,
        0.5
    ),
    geo(
        gra(band(QT_ETM8),8C,0.75,128),
        RasterCRS_ETM8,
        bl,
        1,
        1
    )
)

VtWeb - Sentinel-2 and Landsat-8 pan-sharpening

18.10.2021 by Serge Riazanoff

Look at the result in VtWeb by clicking on this link.

Product list

Sentinel2B/MSI - S2B_MSIL1C_20210826T120319_N0301_R023_T28RBS_20210826T153525

POF-ML script

pan(
    sta(QT_B12,HT_B12,0.02,0.005),
    sta(QT_B11,HT_B11,0.02,0.005),
    sta(QT_B02,HT_B02,0.02,0.005),
    geo(
        gra(
            poly(
                band(
                    QT_B08
                ),
                0,
                1
            ),
            8C,
            0.01,
            0
        ),
        RasterCRS_10,
        bl,
        0,
        0
    )
)

CAFWS - Active fires and burnt areas detection

02.11.2021 by Zhour Najoui

Look at the result in CAFWS by clicking on this link.

Product list

Sentinel2A/MSI - S2A_MSIL1C_20161224T090352_N0204_R007_T34PCS_20161224T090547

POF-ML script

Active fires

lut(for(if(mult(poly(pow(div(mult(mult(QT_B06,QT_B07),QT_B8A),QT_B04),0.5),1,1),poly(div(sub(QT_B12,QT_B8A),pow(add(QT_B12,QT_B8A),0.5)),1,1)),le,-162751.477,uint(255)),over),rainbowsoft,0,255,0)

Burnt areas

lut(if(dem(clsm30),eq,1,for(if(s2lai(QT_B03,QT_B04,QT_B05,QT_B06,QT_B07,QT_B8A,QT_B11,QT_B12,QT_BVZ,QT_BSZ,dif(QT_BSA,QT_BVA)),le,0.428,if(QT_B07,le,2875,if(ndvi(QT_B08,QT_B12),le,0.158,if(QT_B11,le,4624,uint(255)))),if(s2lai(QT_B03,QT_B04,QT_B05,QT_B06,QT_B07,QT_B8A,QT_B11,QT_B12,QT_BVZ,QT_BSZ,dif(QT_BSA,QT_BVA)),le,0.428,if(QT_B07,le,2875,if(ndvi(QT_B08,QT_B12),gt,0.158,if(QT_B02,gt,2096.5,if(QT_B08,le,2118,uint(255))))),if(s2lai(QT_B03,QT_B04,QT_B05,QT_B06,QT_B07,QT_B8A,QT_B11,QT_B12,QT_BVZ,QT_BSZ,dif(QT_BSA,QT_BVA)),le,0.426,if(QT_B07,gt,2923,if(QT_B02,gt,2749.5,if(QT_B11,le,3176,uint(255)))),if(s2lai(QT_B03,QT_B04,QT_B05,QT_B06,QT_B07,QT_B8A,QT_B11,QT_B12,QT_BVZ,QT_BSZ,dif(QT_BSA,QT_BVA)),gt,0.426,if(s2lai(QT_B03,QT_B04,QT_B05,QT_B06,QT_B07,QT_B8A,QT_B11,QT_B12,QT_BVZ,QT_BSZ,dif(QT_BSA,QT_BVA)),le,0.47,if(QT_B11,le,3239,if(ndvi(QT_B08,QT_B12),le,0.169,if(mult(poly(pow(div(mult(mult(QT_B06,QT_B07),QT_B8A),QT_B04),0.5),1,-1),poly(div(sub(QT_B12,QT_B8A),pow(add(QT_B12,QT_B8A),0.5)),1,1)),le,54864.904,uint(255))))),if(s2lai(QT_B03,QT_B04,QT_B05,QT_B06,QT_B07,QT_B8A,QT_B11,QT_B12,QT_BVZ,QT_BSZ,dif(QT_BSA,QT_BVA)),gt,0.426,if(s2lai(QT_B03,QT_B04,QT_B05,QT_B06,QT_B07,QT_B8A,QT_B11,QT_B12,QT_BVZ,QT_BSZ,dif(QT_BSA,QT_BVA)),gt,0.47,if(s2lai(QT_B03,QT_B04,QT_B05,QT_B06,QT_B07,QT_B8A,QT_B11,QT_B12,QT_BVZ,QT_BSZ,dif(QT_BSA,QT_BVA)),le,0.546,if(QT_B03,gt,2573,if(QT_B08,le,3317.5,uint(255)))))))))),over)),drywet,255,0,0)

Sentinel Vision - EVT-023 - DEM isocurve

12.01.2022 by Kevin Gross

Look at the result in Sentinel Vision by clicking on this link.

Product list

DEM/COP-DEM/GLO-30 - COP-DEM_GLO-30-DGED_2020-1

POF-ML script

lut(
    iso(
        QT_elevation,
        2,             // width of the iso value lines (here 1 or 2 pixels)
        step,          // step mode to specify an iso value step
        200            // iso-curve step
    ),
    jet,
    8,                 // minimum index set on the min of the lut
    22,                // maximum index set on the max of the lut
    0
)

VtWeb - Land / water mask on Lac Naivasha (Kenya)

14.02.2022 by Kevin Gross

Look at the result in VtWeb by clicking on this link.

Product list

Sentinel2A/MSI - S2A_MSIL1C_20220113T074301_N0301_R092_T36MZE_20220113T084734

POF-ML script

rgb(
    if(
        QT_B07,
        lt,
        1100,
        sta(QT_B04,HT_B04,0.01,0.8),
        str(QT_B11,1678,7333)
    ),
    if(
        QT_B07,
        lt,
        1100,
        sta(QT_B03,HT_B03,0.01,0.8),
        str(QT_B08,2240,7474)
    ),
    if(
        QT_B07,
        lt,
        1100,
        sta(QT_B02,HT_B02,0.01,0.8),
        str(QT_B02,1308,2714)
    )
)

VtWeb - Sentinel-1 Land/Water RCSE (multi scene)

08.04.2022 by Kevin Gross

Look at the result in VtWeb by clicking on this link.

Product list

Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20220330T064326_20220330T064351_042547_051344_9E81
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20220330T064301_20220330T064326_042547_051344_63B0
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20220330T064236_20220330T064301_042547_051344_3CA3
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20220330T064211_20220330T064236_042547_051344_42A4
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20220330T064146_20220330T064211_042547_051344_E678
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20220330T064121_20220330T064146_042547_051344_B292

POF-ML script

over(
    lut(
        lsm(
            for(sar(stc(pp(s1nc(band(QT_vv))),5:RCS_vv,100,50,500),rasterCrsOrtho,null,null),over),
            clsm30,
            land
        ),
        blues,
        255,
        1,
        0
    ),
    rgb(
        for(str(sar(pp(s1cal(s1nc(band(QT_vv)),g0)),rasterCrsOrtho,cdem30,egm2008),-16,-3),over),
        for(str(sar(pp(s1cal(s1nc(band(QT_vh)),g0)),rasterCrsOrtho,cdem30,egm2008),-22,-10),over),
        for(
            ndi(
                sar(pp(s1cal(s1nc(band(QT_vh)),g0)),rasterCrsOrtho,cdem30,egm2008),
                sar(pp(s1cal(s1nc(band(QT_vv)),g0)),rasterCrsOrtho,cdem30,egm2008),
                -1,
                0.5
            ),
            over
        )
    )
)

VtWeb - MERIS Tristimulus composite (BEAM VISAT)

07.12.2022 by Kevin Gross

Look at the result in VtWeb by clicking on this link.

Product list

ENVISAT/MERIS/FRS - MER_FRS_1PPEPA20090321_105155_000004252077_00266_36893_3403

POF-ML script

rgb(
    str(log(poly(sum(
        poly(poly(geo(band(QT_radiance2),rasterCrs,bc),0,0.010603477),0,0.35),
        poly(poly(geo(band(QT_radiance5),rasterCrs,bc),0,0.00932934),0,0.60),
        poly(geo(band(QT_radiance6),rasterCrs,bc),0,0.008188124),
        poly(poly(geo(band(QT_radiance7),rasterCrs,bc),0,0.0068480745),0,0.13)
    ),1,1)),3.75,4.75),
    str(log(poly(sum(
        poly(poly(geo(band(QT_radiance3),rasterCrs,bc),0,0.01158433),0,0.21),
        poly(poly(geo(band(QT_radiance4),rasterCrs,bc),0,0.01070568),0,0.50),
        poly(geo(band(QT_radiance5),rasterCrs,bc),0,0.00932934),
        poly(poly(geo(band(QT_radiance6),rasterCrs,bc),0,0.008188124),0,0.38)
    ),1,1)),3.86,4.83),
    str(log(poly(sum(
        poly(poly(geo(band(QT_radiance1),rasterCrs,bc),0,0.009474839),0,0.21),
        poly(poly(geo(band(QT_radiance2),rasterCrs,bc),0,0.010603477),0,1.75),
        poly(poly(geo(band(QT_radiance3),rasterCrs,bc),0,0.01158433),0,0.47),
        poly(poly(geo(band(QT_radiance4),rasterCrs,bc),0,0.01070568),0,0.16)
    ),1,1)),4.71,5.31)
)

VtWeb - Sentinel-1 - Parc Marin des Mangroves (mv,mh,cv)

05.12.2022 by Serge RIAZANOFF

Look at the result in VtWeb by clicking on this link.

Product list

Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210107T045215_20210107T045240_036027_0438B1_23D8
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210119T045214_20210119T045239_036202_043EDB_1EC7
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210131T045214_20210131T045239_036377_0444E9_7247
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210212T045213_20210212T045238_036552_044B00_532D
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210224T045213_20210224T045238_036727_045118_2386
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210308T045213_20210308T045238_036902_045730_6259
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210320T045213_20210320T045238_037077_045D4F_F322
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210401T045213_20210401T045238_037252_046355_1339
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210413T045214_20210413T045239_037427_046963_5491
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210425T045215_20210425T045240_037602_046F6D_AB92
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210507T045215_20210507T045240_037777_047569_5B29
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210519T045215_20210519T045240_037952_047AB3_B142
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210531T045216_20210531T045241_038127_047FF3_1745
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210612T045217_20210612T045242_038302_048523_B44D
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210624T045218_20210624T045243_038477_048A5A_BDAE
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210706T045219_20210706T045244_038652_048FA1_C735
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210718T045219_20210718T045244_038827_0494E0_8901
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210730T045220_20210730T045245_039002_049A13_68F6
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210823T045221_20210823T045246_039352_04A5D3_FD7E
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210916T045222_20210916T045247_039702_04B1CF_C91D
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20211127T045222_20211127T045247_040752_04D629_7254
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210928T045223_20210928T045248_039877_04B7D2_8F9A
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20210904T045222_20210904T045247_039527_04ABD6_6114
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20211103T045223_20211103T045248_040402_04CA0B_C2E5
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20211010T045223_20211010T045248_040052_04BDD9_0FF0
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20211115T045223_20211115T045248_040577_04D027_E66F
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20211209T045222_20211209T045247_040927_04DC3A_59CB
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20211022T045223_20211022T045248_040227_04C3F4_C613
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20211221T045221_20211221T045246_041102_04E21A_17DA

POF-ML script

rgb(
    str(
        for(
            sar(pp(s1cal(s1nc(band(
               QT_vv
            )),g0))),rasterCrsOrtho,cdem30,egm2008),
            mean
        ),
        -11,
        -3
    ),
    str(
        for(
            sar(pp(s1cal(s1nc(band(
               QT_vh
               )),g0)),rasterCrsOrtho,cdem30,egm2008),
            mean
        ),
        -19,
        -9
    ),
    str(
        for(
            sar(pp(s1cal(s1nc(band(
               QT_vv
               )),g0)),rasterCrsOrtho,cdem30,egm2008),
            cvf
        ),
        -0.2,
        -0
    )
)

VtWeb - Landsat Collection 2 Level 2 - Surface Temperature

13.09.2023 by Kevin Gross

Look at the result in VtWeb by clicking on this link.

Product list

LANDSAT-C2/OLI-TIRS/L2 - LC09_L2SP_199026_20230605_20230607_02_T1

POF-ML script

lut(
    poly(
        poly(
            QT_B10,
            149,
            0.00341802
        ),
        -273.15,
        1
    ),
    temperature,
    20,
    40,
    0
)

VtWeb - Sentinel-1 - Flooding Valencia

04.11.2024 by Kevin Gross

Look at the result in VtWeb by clicking on this link.

Product list

Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20241031T180256_20241031T180321_056350_06E71E_479F
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20241019T180256_20241019T180321_056175_06E02E_2D52
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20240913T180255_20240913T180320_055650_06CB70_373A
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20240901T180254_20240901T180319_055475_06C483_C134
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20240925T180255_20240925T180320_055825_06D25C_830A
Sentinel1A/CSAR/IW - S1A_IW_GRDH_1SDV_20241007T180256_20241007T180321_056000_06D943_D46B

POF-ML script

rgb(
    str(
        sar(pp(s1cal(s1nc(band(0:QT_vh)),g0)),rasterCrsOrtho,cdem30,egm2008),
        -25,
        -4
    ),
    str(
        sar(pp(s1cal(s1nc(band(0:QT_vh)),g0)),rasterCrsOrtho,cdem30,egm2008),
        -25,
        -4
    ),
    str(
        for(sar(pp(s1cal(s1nc(band(QT_vh)),g0)),rasterCrsOrtho,cdem30,egm2008),mean,1,2,3,4,5),
        -25,
        -4
    )
)