Пакет: matlab.io.hdfeos.gd
Задайте проекцию сетки
defProj(gridID,projCode,zoneCode,sphereCode,projParm)
defProj(gridID,projCode,zoneCode,sphereCode,projParm) задает проекцию GCTP на сетке, заданную как gridID. The projCode аргумент может быть одним из следующих значений.
'geo' | Географический |
'utm' | Универсальный Поперечный Меркатор |
'albers' | Равная площадь Альберса |
'lamcc' | Ламберт Конформный Коник |
'ps' | Полярная стереография |
'polyc' | Поликонический |
'tm' | Поперечный Меркатор |
'lamaz' | Равная площадь Ламберта Азимутала |
'snsoid' | Синусоидальный |
'hom' | Hotine Косой Меркатор |
'som' | Космический Косой Меркатор |
'good' | Прерванный Гуд Гомолозин |
'cea' | Цилиндрическая равная площадь |
'bcea' | Бехрманн, цилиндрическая равная площадь |
'isinus' | Целочисленный Синусоидальный |
Если projCode является 'geo', затем zoneCode, sphereCode, и projParm должно быть задано как []. Любые другие значения для этих параметров игнорируются.
zoneCode - код зоны Универсального Поперечного Меркатора. Он должен быть указан как -1 для других проекций.
sphereCode - имя сфероида GCTP или соответствующий числовой код.
projParm является вектором до 13 элементов, содержащих специфические для проекции параметры. Для получения дополнительной информации о projCode, zoneCode, sphereCode, и projParm, см. главу 6 Руководства пользователя библиотеки HDF-EOS для проекта ECS, том 1: Обзор и примеры.
Эта функция соответствует GDdefproj функция в библиотеке HDF С API.
Создайте сетку UTM, ограниченную долготой от 54 E до 60 E и широтой от 20 N до 30 N (зона 40). Разделите сетку на 120 интервалов по оси X и 200 интервалов по оси Y.
import matlab.io.hdfeos.* gfid = gd.open('myfile.hdf','create'); uplft = [210584.50041 3322395.95445]; lowrgt = [813931.10959 2214162.53278]; gridID = gd.create(gfid,'UTMGrid',120,200,uplft,lowrgt); gd.defProj(gridID,'utm',40,'Clarke 1866',[]); gd.detach(gridID); gd.close(gfid);
Добавьте полярную стереографическую проекцию северного полушария с истинной шкалой при 90 Н, 0 долготы ниже полюса с помощью сфероида WGS 84.
import matlab.io.hdfeos.* gfid = gd.open('myfile.hdf','create'); gridID = gd.create(gfid,'PolarGrid',100,100,[],[]); projparm = zeros(1,13); projparm(6) = 90000000; gd.defProj(gridID,'ps',[],'WGS 84',projparm); gd.detach(gridID); gd.close(gfid);