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[id]=MOD10_L2.A2019047.1110.006.2019047191627.hdf&g=G1592489329-NSIDC_ECS&tl=1549808631!4!!&q=MOD10_L2&ok=MOD10_L2
image_info=gdal.Info(download_path) dataset=re.findall(r'(SUBDATASET_1_NAME)\=([^\=\n]*)\n', image_info) 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 MODIS Reprojection Tool (MRT), which has been the suggested solution in most of the related questions I found online, has been retired this February.