avhrrgoode

Считайте результат обработки данных AVHRR, сохраненный в Проекции Гуда

Синтаксис

[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) считывает данные из набора данных Усовершенствованного радиометра очень Высокого разрешения (AVHRR) с номинальным разрешением 1 км, который хранится в проекции Гуда. Данные в этом формате включают безразмерный индекс растительности наборы данных Глобальных характеристик типа поверхности (GLCC) и (NDVI). region задает географический охват файла, с помощью следующих значений:

  • 'g' или 'global'

  • 'af' или 'africa'

  • 'ap' или 'australia/pacific'

  • 'ea' или 'eurasia'

  • 'na' или 'north america'

  • 'sa' или 'south america'

filename скаляр строки или вектор символов, задающий имя файла данных. Выведите Z определенная геолокацию сетка данных с координатами latgrat и longrat в модулях степеней. Z, latgrat, и longrat имеют класс дважды. Спроектированные координаты, которые лежат в прерванных областях проекции, установлены в NaN. Масштабный коэффициент 100 применяется к исходному набору данных, так, чтобы Z содержит каждую 100-ю точку в обоих направлениях X и Y.

[...] = avhrrgoode(region,filename,scalefactor) использует целочисленный scalefactor прореживать данные. Масштабный коэффициент 1 возвращает каждую точку. Масштабный коэффициент 10 возвращает каждую 10-ю точку. Значение по умолчанию равняется 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,...
nrows,ncols)
заменяет размерности для стандартного формата файла для выбранной области. Этот синтаксис полезен для данных, хранимых на CD-ROM, который, возможно, был усеченным, чтобы соответствовать. Некоторые наборы глобальных данных были распределены с 16 347 строками и 40 031 столбцом данных по CD-ROMs. Размер по умолчанию для наборов глобальных данных является 17 347 строками и 40 031 столбцом данных.

[...] = avhrrgoode(region,filename,scalefactor,latlim,lonlim,gsize,...
nrows,ncols,resolution)
читает набор данных с пространственным разрешением, заданным в метрах. Задайте resolution как любой 1000 или 8000 (метры). Если пустой, полное разрешение 1 000 метров принято. Данные также доступны в 8 000-метровом разрешении. Безразмерные данные об индексе растительности при 8-километровом пространственном разрешении имеют 2 168 строк и 5 004 столбца.

[...] = avhrrgoode(region,filename,scalefactor,latlim,lonlim,gsize,...
nrows,ncols,resolution,precision)
читает набор данных, ожидая целочисленный precision заданный. Если пустой, 'uint8' принят. 'uint16' подходит для некоторых файлов. Проверяйте метаданные (.txt или README) файл в ftp GLCC папка для спецификации формата файла и содержимого. В любом случае, Z преобразован в класс дважды.

Фон

Соединенные Штаты поддерживают семейство спутниковых приборов для измерения изменений климата в рамках программы Earth Observing System (EOS). Предшественниками данных EOS являются наборы данных, полученных NOAA и НАСА по программе Pathfinder. Это данные, полученные с прибора Advanced High Resolution Radiometer (AVHRR), устанавливавшегося на полярно-орбитальных спутниках NOAA, NOAA-7,-9, и-11, и имеющие пространственное разрешение приблизительно 1 км. Данные прибора AVHRR обрабатываются в специальные индексы отдельно для суши, моря и атмосферы. Данные по суше обрабатываются в безразмерный индекс растительности (NDVI) или классификацию типов поверхности и хранятся в двоичных файлах в проекциях Плате-Карре, Гуда и Ламберта. Морские данные обрабатываются в температуры поверхности и хранятся в форматах HDF. avhrrgoode считывает данные о суше, сохраненные в проекции Гуда с глобальным и континентальным покрытием с разрешением 1 км. Также могут считываться 8-километровые данные с глобальным покрытием.

Ограничения

Большинство файлов хранит данные в масштабированных целых числах. Хотя эта функция возвращает данные как double, масштабирование при преобразовании из целого числа в число с плавающей точкой не выполняется. Проверяйте файл данных README на соответствующие масштабные коэффициенты.

Примеры

Пример 1 — прореженное классифицированное глобальное покрытие GLCC

Считайте и отобразите каждую 50-ю точку из файла Глобальных характеристик типа поверхности (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 для Калифорнии

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

Советы

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

Проект AVHRR и наборы данных описаны в и предоставлены различными американскими правительственными веб-сайтами. Смотрите запись для Глобальных характеристик типа поверхности (GLCC) в технических замечаниях, упомянутых ниже.

Примечание

Для получения дополнительной информации на располагающихся данных о карте для загрузки по Интернету, см. следующую документацию на веб-сайте MathWorks: Найдите Картографические данные Онлайн.

Смотрите также

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