Считывание данных AVHRR, хранящихся в Good 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 км, который хранится в проекции Гуда. Данные в этом формате включают в себя наборы данных индекса неперспективной растительности (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 содержит каждую 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,... переопределяет размеры для стандартного формата файла для выбранной области. Этот синтаксис полезен для данных, хранящихся на 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.
В Соединенных Штатах имеется семейство спутниковых датчиков для измерения изменения климата в рамках программы Системы наблюдения Земли (ЭОС). Предшественниками данных ЭОС являются наборы данных, произведенные НОАА и НАСА в рамках программы Pathfinder. Эти данные получены с помощью датчика радиометра высокого разрешения, установленного на спутниках NOAA Polar Orbiter, NOAA-7, -9 и -11, и имеют пространственное разрешение около 1 км. Данные от датчика AVHRR обрабатываются в отдельные показатели суши, моря и атмосферы. Данные о площади суши обрабатываются в соответствии с неразмерным индексом растительности (NDVI) или классификацией земного покрова и хранятся в двоичных файлах в проекциях Плейт-Карре, Гуде и Ламберта. Морские данные обрабатываются до температуры поверхности и хранятся в форматах HDF. avhrrgoode считывает земельные данные, сохраненные в проекции Гуда с глобальным и континентальным охватом на расстоянии 1 км. Он также может считывать данные протяженностью 8 км с глобальным охватом.
Большинство файлов хранят данные в масштабированных целых числах. Хотя эта функция возвращает данные как двойные, масштабирование от целого к плавающему значению не выполняется. Проверьте файл README данных на наличие соответствующих параметров масштабирования.
Считывайте и отображайте каждую 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)
Считывайте и отображайте каждую точку из файла глобальных характеристик земного покрова (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: Find Geospatial Data Online.