I am currently trying to process a MODIS MOD10_L2 product (Snow Cover 5-Min L2 Swath 500m, Version 6) to extract snow cover information. It's in HDF-EOS format and I process it with GDAL in Python and tps interpolation to a geolocated version in EPSG 4326.

The input I used is available here: https://search.earthdata.nasa.gov/search/granules?p=C1000001181-NSIDC_ECS&pg[0][id]=MOD10_L2.A2019047.1110.006.2019047191627.hdf&g=G1592489329-NSIDC_ECS&tl=1549808631!4!!&q=MOD10_L2&ok=MOD10_L2

    dataset=re.findall(r'(SUBDATASET_1_NAME)\=([^\=\n]*)\n', image_info)[0][1]
    ds = gdal.Warp(output_path, dataset, tps=True ,dstSRS='EPSG:4326', outputType=gdal.GDT_Byte)

It works well for regions in the center of the image, giving me great overlay with any type of basemap.

The processed swath My problem is the stripe effect on the edges. It displaces snow cover information significantly. Is there a way to avoid or correct this effect?

enter image description here

The MODIS Reprojection Tool (MRT), which has been the suggested solution in most of the related questions I found online, has been retired this February.

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.