Пакет: matlab.io.hdfeos.gd
Задайте проекцию сетки
defProj(gridID,projCode,zoneCode,sphereCode,projParm)
defProj(gridID,projCode,zoneCode,sphereCode,projParm) задает проекцию GCTP на сетке, заданной gridID. Аргумент projCode может быть одним из этих значений.
'geo' | Географический |
'utm' | Universal, поперечная меркаторский |
'albers' | Алберс коническая равная область |
'lamcc' | Ламберт конформное коническое сечение |
'ps' | Полярный стереографический |
'polyc' | Поликонический |
'tm' | Поперечный меркаторский |
'lamaz' | Ламберт азимутальная равная область |
'snsoid' | Синусоидальный |
'hom' | Hotine, наклонный меркаторский |
'som' | Пробел, наклонный меркаторский |
'good' | Прерванный Гуд Хомолозин |
'cea' | Цилиндрическая равная область |
'bcea' | Берманн цилиндрическая равная область |
'isinus' | Синусоидальный Integerized |
Если projCode является 'geo', то zoneCode, sphereCode и projParm должны быть заданы как []. Любые другие значения для этих параметров проигнорированы.
zoneCode является Универсальным Поперечным Меркаторским зональным кодом. Это должно быть задано как-1 для других проекций.
sphereCode является именем сфероида GCTP или соответствующего цифрового кода.
projParm является вектором до 13 элементов, содержащих специфичные для проекции параметры. Для получения дополнительной информации о projCode, zoneCode, sphereCode и projParm, видят Главу 6 Пользовательского Руководства Библиотеки HDF-EOS для Проекта ECS, Объем 1: Обзор и Примеры.
Эта функция соответствует функции GDdefproj в API библиотеки HDF C.
Создайте сетку UTM, ограниченную 54 E к 60 долготам E и от 20 Н до широты на 30 Н (зона 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);