Некоторые данные прибора MODIS: различия между версиями

Содержимое удалено Содержимое добавлено
Новый учебник.
 
Строка 4:
 
Однако, использование этих данных совместно с распространенным программным обеспечением связано с рядом сложностей (отчасти связанных с дефектами метаданных, отчасти — с несколько необычной структурой используемого в программе EOS формата данных [[w:Hierarchical Data Format|HDF-EOS.]]) Цель данного учебника — познакомить с некоторыми из типов результатов обработки данных MODIS, на примере взаимодействия с ними с использованием программного пакета [[w:QGIS|QGIS.]]
 
== Данные на сетке 0.05° ==
 
Чтобы получить общее представление о данных MODIS ''уровня 3'' («физические величины на регулярной сетке»), рассмотрим данные в системе координат «широта, долгота» на сетке 0.05°. (Данная сетка ориентирована на использование в задачах моделирования климата и, поэтому, также известна как CMG — англ. climate modelling grid.)
 
=== «Надирные» коэффициенты спектральной яркости (MCD43C4) ===
 
Загрузим ([[File:Qgis mActionAddRasterLayer.svg]]) в QGIS, к примеру, растр <code>Nadir_Reflectance_Band1</code> (коэффициенты спектральной яркости, приведенные к стандартной геометрии освещения—наблюдения; канал MODIS 1 — «красный») сетки <code>MCD_CMG_BRDF_0.05Deg</code> файла [http://e4ftl01.cr.usgs.gov/MOTA/MCD43C4.005/2013.06.18/MCD43C4.A2013169.005.2013191204312.hdf MCD43C4.A2013169.005.2013191204312.hdf].
 
[[File:MCD43C4 A2013169 005 2013191204312 band-1.jpeg|240px|right]]
Установим текущую область ''по границам слоя,'' а также изменим ''стиль'' на ''одноканальное псевдоцветное'' и загрузим ''цветовую карту'' <code>reflectance-colors-10000.txt</code>. Результат должен получиться похожим на изображение справа.
 
Здесь уместно упомянуть, что средствами QGIS невозможно объединить три отдельных растра для получения изображения «в истинных цветах» (каналы MODIS 1, 4, 3.) Кроме того, в описании системы координат данного набора допущена ошибка: несложно убедиться, что с точки зрения QGIS, охват слоя — <math>0.1° × 0.05°,</math> а вовсе не ожидаемый <math>360° × 180°.</math>
 
Исправить оба недостатка можно объединив отдельные растры в [http://gdal.org/gdal_vrttut.html <code>.vrt</code>-файл,] подобно:
 
<pre
>$ gdalbuildvrt -separate mcd43c4.a2013169.vrt \
HDF4_EOS:EOS_GRID:"MCD43C4.A2013169.005.2013191204312.hdf":\
MCD_CMG_BRDF_0.05Deg:Nadir_Reflectance_Band{1,4,3}
</pre>
 
После чего в результирующем <code>.vrt</code>-файле исправляется поле <code><GeoTransform /></code>:
 
<source lang="xml"
><GeoTransform> -1.8e+02, 5.e-02, 0.e+00, 9.e+01, 0.e+00, -5.e-02</GeoTransform>
</source>
 
С другой стороны, применить цветовую карту <code>reflectance-colors-10000.txt</code> к каждому из результирующих каналов средствами QGIS окажется невозможным. Впрочем, соответствующая информации также может быть внесена в <code>.vrt</code>-файл.