isequal

Сравните два projcrs или geocrs объекты для эквивалентности

Синтаксис

Описание

пример

tf = isequal(crs1,crs2) возвращает логический 1 TRUE) если системы координат (CRSs) crs1 и crs2 эквивалентны. В противном случае это возвращает логический 0 ложь).

  • Два projcrs объекты эквивалентны, если у них есть тот же географический CRS, метод проекции, параметры проекции и единица длины.

  • Два geocrs объекты эквивалентны, если у них есть то же имя данной величины, опорный сфероид, нулевой меридиан и угловой модуль.

isequal функция не рассматривает имена CRS.

Примеры

свернуть все

Получите информацию приблизительно два растровых набора данных с помощью georasterinfo функция. Присвойте их системы координат переменным. В этом случае системами координат является projcrs объекты.

info1 = georasterinfo('MtWashington-ft.grd');
info2 = georasterinfo('MountDartmouth-ft.grd');
 
crs1 = info1.CoordinateReferenceSystem;
crs2 = info2.CoordinateReferenceSystem;

Сравните системы координат.

tf = isequal(crs1,crs2)
tf = logical
   1

Две системы координат могут быть эквивалентными, даже если строки, возвращенные wktstring функция не идентична.

Возвратите информацию о файле форм как структура. Найдите систему координат путем запроса CoordinateReferenceSystem поле .

info = shapeinfo('landareas.shp');
g1 = info.CoordinateReferenceSystem
g1 = 
  geocrs with properties:

             Name: "WGS 84"
            Datum: "World Geodetic System 1984"
         Spheroid: [1x1 referenceEllipsoid]
    PrimeMeridian: 0
        AngleUnit: "degree"

Создайте вторую систему координат при помощи geocrs функция.

g2 = geocrs(4326)
g2 = 
  geocrs with properties:

             Name: "WGS 84"
            Datum: "World Geodetic System 1984"
         Spheroid: [1x1 referenceEllipsoid]
    PrimeMeridian: 0
        AngleUnit: "degree"

Системы координат имеют то же имя данной величины, опорный сфероид, нулевой меридиан и угловой модуль. Поэтому системы координат эквивалентны.

isequal(g1,g2)
ans = logical
   1

Обратите внимание на то, что строки WKT не идентичны. Строка для g2 содержит дополнительную информацию об области использования ("USAGE") и код полномочий ("ID").

wktstring(g1,'Format','Formatted')
ans = 
    "GEOGCRS["WGS 84",
         DATUM["World Geodetic System 1984",
             ELLIPSOID["WGS 84",6378137,298.257223563,
                 LENGTHUNIT["metre",1]],
             ID["EPSG",6326]],
         PRIMEM["Greenwich",0,
             ANGLEUNIT["Degree",0.0174532925199433]],
         CS[ellipsoidal,2],
             AXIS["longitude",east,
                 ORDER[1],
                 ANGLEUNIT["Degree",0.0174532925199433]],
             AXIS["latitude",north,
                 ORDER[2],
                 ANGLEUNIT["Degree",0.0174532925199433],
                 ID["EPSG",7030]]]"

wktstring(g2,'Format','Formatted')
ans = 
    "GEOGCRS["WGS 84",
         DATUM["World Geodetic System 1984",
             ELLIPSOID["WGS 84",6378137,298.257223563,
                 LENGTHUNIT["metre",1]]],
         PRIMEM["Greenwich",0,
             ANGLEUNIT["degree",0.0174532925199433]],
         CS[ellipsoidal,2],
             AXIS["geodetic latitude (Lat)",north,
                 ORDER[1],
                 ANGLEUNIT["degree",0.0174532925199433]],
             AXIS["geodetic longitude (Lon)",east,
                 ORDER[2],
                 ANGLEUNIT["degree",0.0174532925199433]],
         USAGE[
             SCOPE["Horizontal component of 3D system."],
             AREA["World."],
             BBOX[-90,-180,90,180]],
         ID["EPSG",4326]]"

Входные параметры

свернуть все

Первая система координат в виде projcrs или geocrs объект.

Вторая система координат в виде projcrs или geocrs объект.

Советы

  • Если вы ожидаете, что две системы координат будут эквивалентны, но isequal функция возвращает 0 ложь), можно смотреть строки WKT для различий с помощью wktstring функция.

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

Объекты

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