isequal

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

Свернуть все на странице

Синтаксис

tf = isequal(crs1,crs2)

Описание

пример

tf = isequal(crs1,crs2) возвращает логический 1 (true), если системы координат-привязок (CRS) crs1 и crs2 являются эквивалентными. В противном случае возвращается логический 0 (false).

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

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

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

Примеры

свернуть все

Открыть Live Script

Получите информацию о двух наборах растровых данных с помощью 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
Открыть Live Script

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

Возвращает информацию о файле shapefile как структуре. Найдите систему координат-привязок путем запроса 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["unknown"],
             AREA["World"],
             BBOX[-90,-180,90,180]],
         ID["EPSG",4326]]"

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

свернуть все

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

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

Совет

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

См. также

Объекты

Введенный в R2021a

Mapping Toolbox

Поддержка

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