avhrrlambert
Чтение продукта данных AVHRR, хранящегося в eqaazim проектирование
Синтаксис
[latgrat,longrat,Z] = avhrrlambert(region,filename)
[...] = avhrrlambert(region,filename,
scalefactor)
[...] = avhrrlambert(region,filename,
scalefactor, latlim, lonlim)
[...] = avhrrlambert(region,filename,
scalefactor, latlim, lonlim, gsize)
[...] = avhrrlambert(region,filename,
scalefactor, latlim, lonlim, gsize,precision)
Описание
[latgrat,longrat,Z] = avhrrlambert( считывает данные из набора данных Advanced Very High Resolution Radiometer (AVHRR) с номинальным разрешением 1 км, который хранится в проекции Lambert Equal Azimuthal. Данные этого типа включают глобальные характеристики типа поверхности (GLCC). region,filename)region задает охват файла. Допустимые области перечислены в следующей таблице. filename - строка, задающая имя файла данных. Z - геолокированная сетка данных с координатами latgrat и longrat в единицах степеней. Коэффициент шкалы 100 применяется к исходному набору данных таким образом, что Z содержит каждую 100-ю точку в обоих X и Y.
| Спецификаторы областей |
|---|
'a' или 'asia' |
'af' или 'africa' |
'ap' или 'australia/pacific' |
'e' или 'europe' |
'na' или 'north america' |
'sa' или 'south america' |
[...] = avhrrlambert( использует целое число region,filename,
scalefactor)scalefactor чтобы уменьшить значение данных. Коэффициент шкалы 1 возвратов каждую точку. Коэффициент шкалы 10 возвратов каждый 10-й пункт. Значение по умолчанию является 100.
[...] = avhrrlambert( возвращает данные для заданной области. Результат может несколько выйти за пределы запрашиваемой области. Пределы являются двухэлементными векторами в единицах степеней, с region,filename,
scalefactor, latlim, lonlim)latlim в области значений [-90 90] и lonlim в области значений [-180 180]. Если latlim и lonlim пусты, возвращается вся область, покрытая файлом данных. Если четырехугольник задан как latlim и lonlim (при проецировании с образованием многоугольника в соответствующей азимутальной проекции Lambert Equal) не пересекает ограничивающий прямоугольник данных в проективных координатах, затем latgrat, longrat, и Z пусты.
[...] = avhrrlambert( управляет размером матриц гратикулы. region,filename,
scalefactor, latlim, lonlim, gsize)gsize - двухэлементный вектор, содержащий необходимое количество строк и столбцов. Если он опущен или пуст, возвращается графическое изображение размера сетки.
[...] = avhrrlambert( считывает набор данных с целым числом region,filename,
scalefactor, latlim, lonlim, gsize,precision)precision указано. Если опущено, 'uint8' принято. 'uint16' подходит для некоторых файлов. Проверьте метаданные (.txt или README) файл в папке ftp для спецификации формата файла и содержимого.
Фон
Соединенные Штаты планируют создать семейство спутниковых приборов для измерения изменений климата в рамках программы Earth Observing System (EOS). Ранними предшественниками данных EOS являются наборы данных, полученных NOAA и НАСА по программе Pathfinder. Это данные, полученные с прибора Advanced High Resolution Radiometer, установленного на полярно-орбитальных спутниках NOAA, NOAA-7, -9 и -11 с пространственным разрешением около 1 км. Этот данный прибора AVHRR обрабатывается в специальные индексы отдельно для суши, моря и атмосферы. Данные о земельных участках обрабатываются в соответствии с недиметрическим индексом растительности или классификацией типа поверхности и хранятся в двоичных файлах в проекциях Plate Carrée, Goode и Lambert Equal Azimuthal. Данные о море обрабатываются до температур поверхности и хранятся в форматах HDF. Эта функция считывает данные типы поверхности континентов, сохраненные в проекции Lambert Equal Area Azimuthal на 1 км.
Примеры
Пример 1
Считывайте и отображайте каждую 100-ю точку из файла Global Land Cover Characteristics (GLCC), охватывающего Северную Америку с помощью схемы классификации USGS, названной nausgs1_2l.img. Чтобы запустить этот пример, необходимо сначала загрузить файл.
[latgrat, longrat, Z] = avhrrlambert('na','nausgs1_2l.img');Отобразите данные с помощью проекции Lambert Equal Area Azimuthal.
origin = [50 -100 0];
ellipsoid = [6370997 0];
figure
axesm('MapProjection', 'eqaazim', 'Origin', ...
origin, 'Geoid', ellipsoid)
geoshow(latgrat, longrat, Z, 'DisplayType', 'texturemap'); 
Пример 2
Прочтите и отобразите каждую другую точку из файла Global Land Cover Characteristics (GLCC), охватывающего Европу с помощью схемы классификации USGS, названной eausgs1_2le.img. Чтобы запустить этот пример, необходимо сначала загрузить файл.
figure
worldmap france
mstruct = gcm;
latlim = mstruct.maplatlimit;
lonlim = mstruct.maplonlimit;
scalefactor = 2;
[latgrat, longrat, Z] = avhrrlambert('e', 'eausgs1_2le.img', ...
scalefactor, latlim, lonlim);
geoshow(latgrat, longrat, Z, 'DisplayType', 'texturemap');
geoshow('landareas.shp','FaceColor','none','EdgeColor','black')
Совет
Эта функция читает двоичные файлы как есть. Вы не должны использовать программное обеспечение для замены байтов в этих файлах.
Проект АВХРР и наборы данных описываются и предоставляются различными веб-сайтами правительства США.
Примечание
Для получения дополнительной информации о поиске данных карты для загрузки через Интернет см. следующую документацию на веб-сайте MathWorks: Find Geospatial Data Online.