exponenta event banner

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(region,filename) считывает данные с усовершенствованного радиометра очень высокого разрешения (AVHRR) с номинальным разрешением 1 км, который хранится в проекции Lambert Equal Area Azimuthal. Данные этого типа включают глобальные характеристики земного покрова (GLCC). region указывает охват файла. Допустимые области перечислены в следующей таблице. filename - строка, задающая имя файла данных. Z - геоориентированная сетка данных с координатами latgrat и longrat в единицах степеней. Масштабный коэффициент 100 применяется к исходному набору данных так, что Z содержит каждую сотую точку как в 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 (при проецировании с образованием многоугольника в соответствующей проекции азимутальной площади Ламберта) не пересекает ограничивающую рамку данных в спроецированных координатах, 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 для указания формата и содержимого файла.

Фон

Соединенные Штаты планируют построить семейство спутниковых датчиков для измерения изменения климата в рамках программы Системы наблюдения Земли (EOS). Ранними предшественниками данных ЭОС являются наборы данных, производимые НОАА и НАСА в рамках программы Pathfinder. Это данные, полученные с помощью датчика радиометра высокого разрешения, выполненного на спутниках NOAA Polar Orbiter, NOAA-7, -9 и -11 с пространственным разрешением около 1 км. Данные от датчика AVHRR обрабатываются в отдельные показатели суши, моря и атмосферы. Данные о площади суши обрабатываются в соответствии с неразмерным индексом растительности или классификацией растительного покрова и хранятся в двоичных файлах в проекциях Плита Карре, Гуде и Ламберт Равная площадь Азимуталь. Морские данные обрабатываются до температуры поверхности и хранятся в форматах HDF. Эта функция считывает данные о земном покрове для континентов, сохраненные в проекции Азимутальского равноденственного района Ламберта на расстоянии 1 км.

Примеры

Пример 1

Считывайте и отображайте каждую сотую точку из файла глобальных характеристик земного покрова (GLCC), охватывающего Северную Америку со схемой классификации USGS, названной nausgs1_2l.img. Для выполнения этого примера необходимо сначала загрузить файл.

[latgrat, longrat, Z] = avhrrlambert('na','nausgs1_2l.img');

Просмотрите данные с помощью проекции Lambert Equal Azimuthal.

origin = [50 -100 0];
ellipsoid = [6370997 0];
figure
axesm('MapProjection', 'eqaazim', 'Origin', ...
   origin, 'Geoid', ellipsoid)
geoshow(latgrat, longrat, Z, 'DisplayType', 'texturemap');  

Пример 2

Считывание и отображение каждой другой точки из файла глобальных характеристик земного покрова (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')

Совет

Эта функция считывает двоичные файлы как есть. В этих файлах не следует использовать программное обеспечение для замены байтов.

Проект AVHRR и наборы данных описаны и предоставлены различными веб-сайтами правительства США.

Примечание

Подробные сведения о поиске данных карты для загрузки через Интернет см. в следующей документации на веб-сайте MathWorks: Find Geospatial Data Online.

См. также

Представлен до R2006a