avhrrgoode
Чтение продукта данных AVHRR, хранящегося в Goode Projection
Синтаксис
[latgrat,longrat,z] = avhrrgoode(region,filename)
[...] = avhrrgoode(region,filename,scalefactor)
[...] = avhrrgoode(region,filename,scalefactor,latlim,lonlim)
[...] = avhrrgoode(region,filename,scalefactor,latlim,lonlim,gsize)
[...] = avhrrgoode(region,filename,scalefactor,latlim,lonlim,gsize,...
nrows,ncols)
[...] = avhrrgoode(region,filename,scalefactor,latlim,lonlim,gsize,...
nrows,ncols,resolution)
[...] = avhrrgoode(region,filename,scalefactor,latlim,lonlim,gsize,...
nrows,ncols,resolution,precision)
Описание
[latgrat,longrat,z] = avhrrgoode(region,filename) считывает данные из набора данных Advanced Very High Resolution Radiometer (AVHRR) с номинальным разрешением 1 км, которое хранится в проекции Goode. Данные в этом формате включают в себя наборы данных по недиметрическому индексу растительности (NDVI) и глобальным характеристикам типа поверхности (GLCC). region определяет географический охват файла с помощью следующих значений:
'g'или'global''af'или'africa''ap'или'australia/pacific''ea'или'eurasia''na'или'north america''sa'или'south america'
filename - строковый скаляр или вектор символов, задающий имя файла данных. Выходные Z - геолокированная сетка данных с координатами latgrat и longrat в единицах степеней. Z, latgrat, и longrat имеют класс double. Проецируемые координаты, которые находятся в пределах прерывистых областей проекции, устанавливаются в NaN. Коэффициент шкалы 100 применяется к исходному набору данных, так что Z содержит каждые 100th указать в обоих направлениях X и Y.
[...] = avhrrgoode(region,filename,scalefactor) использует целое число scalefactor чтобы уменьшить значение данных. Коэффициент шкалы 1 возвратов каждую точку. Коэффициент шкалы 10 возвратов каждые 10th точка. Значение по умолчанию является 100.
[...] = avhrrgoode(region,filename,scalefactor,latlim,lonlim) возвращает данные для заданной области. Возвращенные данные могут несколько выходить за пределы запрашиваемой области. Пределы являются двухэлементными векторами в единицах степеней, с latlim в области значений [-90 90] и lonlim в области значений [-180 180]. latlim и lonlim должен быть возрастающим. Если latlim и lonlim пусты, возвращается вся область, покрытая файлом данных. Если четырехугольник задан как latlim и lonlim (при проецировании на образование многоугольника в соответствующей проекции Гуда) не пересекает ограничивающий прямоугольник данных в проективных координатах, затем Z, latgrat, и longrat возвращаются как пустые.
[...] = avhrrgoode(region,filename,scalefactor,latlim,lonlim,gsize) управляет размером матриц гратикулы. gsize - двухэлементный вектор, содержащий необходимое количество строк и столбцов. По умолчанию latgrat, и longrat имеют тот же размер, что и Z.
[...] = avhrrgoode(region,filename,scalefactor,latlim,lonlim,gsize,... переопределяет размерности для стандартного формата файла для выбранной области. Этот синтаксис полезен для данных, хранящихся на CD-ROM, которые, возможно, были усечены для соответствия. Некоторые наборы глобальных данных были распространены с 16347 строками и 40031 столбцами данных на компакт-дисках. Размер по умолчанию для глобальных наборов данных - 17347 строк и 40031 столбец данных.
nrows,ncols)
[...] = avhrrgoode(region,filename,scalefactor,latlim,lonlim,gsize,...считывает набор данных с пространственным разрешением, заданным в метрах. Задайте
nrows,ncols,resolution)resolution либо как 1000 или 8000 (м). Если он пуст, принимается полное разрешение 1000 метров. Данные также доступны при разрешении 8000 метров. Данные индекса необмерной растительности при 8-км пространственном разрешении имеют 2168 строк и 5004 столбца.
[...] = avhrrgoode(region,filename,scalefactor,latlim,lonlim,gsize,... считывает набор данных, ожидающий целое число
nrows,ncols,resolution,precision)precision указано. Если пусто, 'uint8' принято. 'uint16' подходит для некоторых файлов. Проверьте метаданные (.txt или README) файл в GLCC ftp папка для спецификации формата файла и содержимого. В любом случае Z преобразуется в класс double.
Фон
Соединенные Штаты поддерживают семейство спутниковых приборов для измерения изменений климата в рамках программы Earth Observing System (EOS). Предшественниками данных EOS являются наборы данных, полученных NOAA и НАСА по программе Pathfinder. Это данные, полученные с прибора Advanced High Resolution Radiometer, летавшего на полярно-орбитальных спутниках NOAA, NOAA-7, -9 и -11, и имеющие пространственные разрешения около 1 км. Этот данный прибора AVHRR обрабатывается в специальные индексы отдельно для суши, моря и атмосферы. Данные о земельных площадях обрабатываются в недиметрический индекс растительности (NDVI) или классификацию типа поверхности и хранятся в двоичных файлах в проекциях Плате-Карре, Гуда и Ламберта. Данные о море обрабатываются до температур поверхности и хранятся в форматах HDF. avhrrgoode считывает данные о земле, сохраненные в прогнозе Гуда с глобальным и континентальным покрытием с разрешением 1 км. Также могут считываться 8-километровые данные с глобальным покрытием.
Ограничения
Большинство файлов сохраняют данные в масштабированных целых числах. Хотя эта функция возвращает данные как double, масштабирование из целого числа в число с плавающей точкой не выполняется. Проверьте файл README данных на соответствие параметрам масштабирования.
Примеры
Пример 1 - Классифицированное покрытие GLCC с пониженной дискретизацией
Считывайте и отображайте каждую 50-ю точку из файла Global Land Cover Characteristics (GLCC), охватывающего весь земной шар, с помощью схемы классификации USGS, названной gusgs2_0g.img. (Чтобы запустить пример, необходимо сначала загрузить файл.)
[latgrat, longrat, Z] = avhrrgoode('global', ...
'gusgs2_0g.img',50);
% Convert the geolocated data grid to an geolocated image.
uniqueClasses = unique(Z);
RGB = ind2rgb8(uint8(Z), jet(numel(uniqueClasses)));
% Display the data as an image using the Goode projection.
origin = [0 0 0];
ellipsoid = [6370997 0];
figure
axesm('MapProjection', 'goode', 'Origin', origin, ...
'Geoid', ellipsoid)
geoshow(latgrat, longrat, RGB, 'DisplayType', 'image');
axis image off
% Plot the coastlines.
hold on
load coastlines
plotm(coastlat,coastlon)
Пример 2 - Классифицированные данные GLCC для Калифорнии
Прочтите и отобразите каждую точку из файла Global Land Cover Characteristics (GLCC), охватывающего Калифорнию с помощью схемы классификации USGS, названной nausgs1_2g.img. Сначала необходимо загрузить файл, чтобы запустить этот пример.
figure
usamap california
mstruct = gcm;
latlim = mstruct.maplatlimit;
lonlim = mstruct.maplonlimit;
scalefactor = 1;
[latgrat, longrat, Z] = ...
avhrrgoode('na', 'nausgs1_2g.img', scalefactor, latlim, lonlim);
geoshow(latgrat, longrat, Z, 'DisplayType', 'texturemap');
% Overlay vector data from usastatehi.shp.
california = shaperead('usastatehi', 'UseGeoCoords', true,...
'BoundingBox', [lonlim;latlim]);
geoshow([california.Lat], [california.Lon], 'Color', 'black');
Совет
Эта функция читает двоичные файлы как есть. Вы не должны использовать программное обеспечение для замены байтов в этих файлах.
Проект АВХРР и наборы данных описываются и предоставляются различными веб-сайтами правительства США. См. запись «Глобальные характеристики типа поверхности» (GLCC) в технической записке, приведенной ниже.
Примечание
Для получения дополнительной информации о поиске данных карты для загрузки через Интернет см. следующую документацию на веб-сайте MathWorks: Find Geospatial Data Online.