mfwdtran

Географические особенности проекта, чтобы сопоставить координаты

Синтаксис

[x,y] = mfwdtran(lat,lon)
[x,y,z] = mfwdtran(lat,lon,alt)
[...] = mfwdtran(mstruct,...)

Описание

[x,y] = mfwdtran(lat,lon) применяет прямое преобразование, заданное проекцией карты в текущих осях карты. Можно использовать эту функцию, чтобы преобразовать местоположения точки и строку и вершины полигона, данные в широтах и долготах к плоской, спроектированной системе координат карты.

[x,y,z] = mfwdtran(lat,lon,alt) применяет прямую проекцию к 3-D входу, приводящему к 3-D выводу. Если вход alt пуст или не использован, то alt = 0 принят.

[...] = mfwdtran(mstruct,...) требует допустимой структуры проекции карты в качестве первого аргумента. В этом случае, никакие оси карты необходим.

Примеры

свернуть все

Получите данные о географическом положении для округа Колумбия.

dc = shaperead('usastatelo', 'UseGeoCoords', true,...
    'Selector',{@(name) strcmpi(name,'District of Columbia'),...
    'Name'});
lat = [dc.Lat]';
lon = [dc.Lon]';
[lat lon]
ans = 10×2

   38.9000  -77.0700
   38.9500  -77.1200
   39.0000  -77.0300
   38.9000  -76.9000
   38.7800  -77.0300
   38.8000  -77.0200
   38.8700  -77.0200
   38.9000  -77.0700
   38.9000  -77.0700
       NaN       NaN

Получите зону UTM для этой точки.

dczone = utmzone(mean(lat,'omitnan'),mean(lon,'omitnan'))
dczone = 
'18S'

Настройте систему координат UTM на основе этой информации.

utmstruct = defaultm('utm'); 
utmstruct.zone = dczone;  
utmstruct.geoid = wgs84Ellipsoid;
utmstruct = defaultm(utmstruct)
utmstruct = struct with fields:
     mapprojection: 'utm'
              zone: '18S'
        angleunits: 'degrees'
            aspect: 'normal'
     falsenorthing: 0
      falseeasting: 500000
       fixedorient: []
             geoid: [1x1 referenceEllipsoid]
       maplatlimit: [32 40]
       maplonlimit: [-78 -72]
      mapparallels: []
        nparallels: 0
            origin: [0 -75 0]
       scalefactor: 0.9996
           trimlat: [-80 84]
           trimlon: [-180 180]
             frame: 'off'
             ffill: 100
        fedgecolor: [0.1500 0.1500 0.1500]
        ffacecolor: 'none'
         flatlimit: [32 40]
        flinewidth: 2
         flonlimit: [-3 3]
              grid: 'off'
         galtitude: Inf
            gcolor: [0.1500 0.1500 0.1500]
        glinestyle: ':'
        glinewidth: 0.5000
    mlineexception: []
         mlinefill: 100
        mlinelimit: []
     mlinelocation: 1
      mlinevisible: 'on'
    plineexception: []
         plinefill: 100
        plinelimit: []
     plinelocation: 1
      plinevisible: 'on'
         fontangle: 'normal'
         fontcolor: [0.1500 0.1500 0.1500]
          fontname: 'Helvetica'
          fontsize: 10
         fontunits: 'points'
        fontweight: 'normal'
       labelformat: 'compass'
     labelrotation: 'off'
        labelunits: 'degrees'
     meridianlabel: 'off'
    mlabellocation: 1
    mlabelparallel: 40
       mlabelround: 0
     parallellabel: 'off'
    plabellocation: 1
    plabelmeridian: -78
       plabelround: 0

Спроектируйте данные округа Колумбия из географических координат в координаты карты для зональных 18 UTM.

[x,y] = mfwdtran(utmstruct,lat,lon)
x = 10×1
105 ×

    3.2049
    3.1629
    3.2421
    3.3524
    3.2367
    3.2459
    3.2476
    3.2049
    3.2049
       NaN

y = 10×1
106 ×

    4.3077
    4.3134
    4.3187
    4.3074
    4.2943
    4.2965
    4.3043
    4.3077
    4.3077
       NaN

Смотрите также

| | | | | |

Представлено до R2006a

Для просмотра документации необходимо авторизоваться на сайте