Пакет: matlab.io.hdfeos.gd
Определение проекции сетки
defProj(gridID,projCode,zoneCode,sphereCode,projParm)
defProj(gridID,projCode,zoneCode,sphereCode,projParm) определяет проекцию GCTP на сетку, заданную gridID. projCode аргумент может быть одним из этих значений.
'geo' | Географический |
'utm' | Универсальный поперечный меркатор |
'albers' | Коническая равная площадь Альберса |
'lamcc' | Конформный коник Ламберта |
'ps' | Полярный стереографический |
'polyc' | Поликонический |
'tm' | Поперечный Меркатор |
'lamaz' | Равная площадь Ламберта Азимуталя |
'snsoid' | Синусоидальный |
'hom' | Горячий косой меркатор |
'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 в API библиотеки HDF.
Создайте сетку UTM, ограниченную от 54 до 60 E долготы и от 20 до 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);