Veel

Kummaline efekt gdalwarp väljundis

Kummaline efekt gdalwarp väljundis


Ma kasutan pythonis gdalwarpit, et eraldada ja uuesti projitseerida peegeldusribasid MODIS-i tasemelt 1B (MOD021 km). Minu kasutatava käsu näide on järgmine:

gdalwarp -of GTIFF -tps -t_srs EPSG: 4326 -tr 0.01 0.01 -r bilinear -overwrite HDF4_EOS: EOS_SWATH: "C:  tmp  MYD021KM.A2007094.0935.006.2012075055042.hdf": MODIS_SWATH_Tg_Tg_Tg_Type .tif

See teisendab hdf alamandmekogumi ilma probleemideta, kuid väljundfailil näib olevat mingi artefakt, eriti servade lähedal. (Vaata allpool)

Teisendasin samad andmed Modis Conversion Toolkit (MCTK) abil ja tundub, et see on palju parem. (Allpool)

Proovisin erinevaid resamplingumeetodeid ja erinevaid väljundprojektsioone, kuid tulemused olid samad! Kas kellelgi on seda probleemi gdalwarpiga?

P.S. Andmeid saab alla laadida sellelt aadressilt: ftp://ladsweb.nascom.nasa.gov/allData/6/MYD021KM/2007/094/MYD021KM.A2007094.0935.006.2012075055042.hdf


Nii et pärast seda probleemi arutamist teiste arendajatega Gdal-Dev meililistis selgus, et MODIS-i vaalu andmeid ei saa gdalwarpi abil ideaalselt ümber projitseerida. Selle põhjuseks on MODIS skaneerimisanduri kikilipsu efekt, millega gdal ei saa hakkama. Parim võimalus on alla laadida seotud geograafilise asukoha andmekogumid (MOD03) ja kasutada soovitud alamandmekogumite uuesti projitseerimiseks MRTswathi.


Vaata videot: Introduction to GIS Analysis with GeoPandas using Python