Пакет: matlab.io.hdfeos.gd
Пиксельные строки и столбцы для пар широты/долготы
[row,col] = getPixels(gridID,lat,lon)
[row,col] = getPixels(gridID,lat,lon) преобразует пары широты/долготы в основанные на нуле пиксельные координаты строки и столбца. Источник является верхним левым углом пикселя сетки. Если пары широты/долготы вне сетки, то row и col -1.
Эта функция соответствует GDgetpixels функция в API библиотеки HDF-EOS C.
import matlab.io.hdfeos.* gfid = gd.open('grid.hdf','read'); gridID = gd.attach(gfid,'PolarGrid'); cornerlat = [20 50]; cornerlon = [-90 -60]; [row,col] = gd.getPixels(gridID,cornerlat,cornerlon); gd.detach(gridID); gd.close(gfid);