Чтение изображения JPEG из файла.

basename = 'boston_ovr';
imagefile = [basename '.jpg'];
RGB = imread(imagefile);

Извлеките имя файла мира из имени файла изображения, прочитайте файл мира и создайте объект пространственной ссылки.

worldfile = getworldfilename(imagefile);
R = worldfileread(worldfile, 'geographic', size(RGB));

Запись данных изображения и ссылок на данные в файл GeoTIFF.

filename = [basename '.tif'];
geotiffwrite(filename, RGB, R)

Создайте пустые оси карты и отобразите карту.

figure
usamap(RGB, R)
geoshow(filename)

Преобразование классического файла TIFF с географической привязкой в файл BigTIFF с плиткой путем извлечения информации из классического файла TIFF. Сначала импортируйте классическое изображение TIFF Бостона и объект ссылки на ячейки карты. Получение метаданных из файла с помощью geotiffinfo.

infilename = 'boston.tif';
[A,R] = readgeoraster(infilename);
info = geotiffinfo(infilename);

Укажите теги для включения в файл BigTIFF с мозаикой. Для этого извлеките тег каталога GeoKey из метаданных. Затем создайте теги, задающие длину и ширину плиток.

geoTags = info.GeoTIFFTags.GeoKeyDirectoryTag;
tiffTags = struct('TileLength',1024,'TileWidth',1024);

Запишите данные в новый файл GeoTIFF. Укажите формат файла как BigTIFF с помощью 'TiffType' пара имя-значение. Включить теги, указав 'GeoKeyDirectoryTag' и 'TiffTags' пары имя-значение.

outfilename = 'boston_bigtiff.tif';
geotiffwrite(outfilename,A,R,'TiffType','bigtiff', ...
                             'GeoKeyDirectoryTag',geoTags, ...
                             'TiffTags',tiffTags)

Убедитесь, что вы написали файл BigTIFF, прочитав файл и запросив теги.

biginfo = geotiffinfo(outfilename);
biginfo.GeoTIFFTags.GeoKeyDirectoryTag

ans = struct with fields:
        GTModelTypeGeoKey: 1
       GTRasterTypeGeoKey: 1
    ProjectedCSTypeGeoKey: 26986
        PCSCitationGeoKey: 'State Plane Zone 2001 NAD = 83'
    ProjLinearUnitsGeoKey: 9003

t = Tiff(outfilename);
getTag(t,'TileLength')

ans = 1024

getTag(t,'TileWidth')

ans = 1024

close(t)

Прочитайте два смежных ортофота и объедините их.

X_west = imread('concord_ortho_w.tif');
X_east = imread('concord_ortho_e.tif');
X = [X_west X_east];

Создайте ссылочные объекты для ортофотов и для их комбинации.

R_west = worldfileread('concord_ortho_w.tfw', 'planar', size(X_west));
R_east = worldfileread('concord_ortho_e.tfw', 'planar', size(X_east));
R = R_west;
R.XLimWorld = [R_west.XLimWorld(1) R_east.XLimWorld(2)];
R.RasterSize = size(X);

Запишите комбинированное изображение в файл GeoTIFF. Используйте кодовый номер (26986), указывающий PCS_NAD83_Massachusetts Спроецированную систему координат (Projected Coordinate System).

coordRefSysCode = 26986;
filename = 'concord_ortho.tif';
geotiffwrite(filename, X, R, 'CoordRefSysCode', coordRefSysCode);

Отображение карты.

figure
mapshow(filename)

Импорт изображения GeoTIFF и отображение ссылочного объекта ячеек для области вокруг Бостона с помощью readgeoraster.

[A,RA] = readgeoraster('boston.tif');

Обрезать данные до пределов, указанных в xlimits и ylimits использование mapcrop.

xlimits = [764318 767678];
ylimits = [2951122 2954482];
[B,RB] = mapcrop(A,RA,xlimits,ylimits);

Получение информации об изображении GeoTIFF с помощью geotiffinfo. Извлеките тег каталога GeoKey из информации.

info = geotiffinfo('boston.tif');
key = info.GeoTIFFTags.GeoKeyDirectoryTag;

Запишите обрезанные данные и тег каталога GeoKey в файл. Убедитесь, что обрезанные данные записаны в файл, отображая его.

filename = 'boston_subimage.tif';
geotiffwrite(filename,B,RB,'GeoKeyDirectoryTag',key)
figure
mapshow(filename)

Запишите в файл GeoTIFF данные об отметке для области вокруг Южного Боулдера в Колорадо. Сначала импортируйте данные фасада и объект ссылки на географические проводки.

[Z,R] = readgeoraster('n39_w106_3arc_v2.dt1','OutputType','double');

Укажите информацию о теге каталога GeoKey для файла GeoTIFF в качестве структуры. Укажите, что данные находятся в географической системе координат, указав GTModelTypeGeoKey поле 2. Укажите, что ссылочный объект использует проводки (а не ячейки), указав GTRasterTypeGeoKey поле 2. Укажите, что данные привязаны к географической системе координат, указав GeographicTypeGeoKey поле 4326.

key.GTModelTypeGeoKey  = 2;
key.GTRasterTypeGeoKey = 2;
key.GeographicTypeGeoKey = 4326;

Запишите данные и тег каталога GeoKey в файл.

filename = 'southboulder.tif';
geotiffwrite(filename,Z,R,'GeoKeyDirectoryTag',key)

Убедитесь, что данные записаны в файл, отображая их на карте.

usamap([39 40],[-106 -105])
g = geoshow(filename,'DisplayType','mesh');
demcmap(g.CData)

Данные о высотах, используемые в этом примере, предоставлены Геологической службой США.

Создайте образец файла TIFF с метаданными RPC. Для этого создайте массив нулей и связанный объект ссылки.

A = zeros(180,360);
latlim = [-90 90];
lonlim = [-180 180];
RA = georefcells(latlim,lonlim,size(A));

Затем создайте RPCCoefficientTag и задайте некоторые поля с типовыми значениями. RPCCoefficientTag объект представляет метаданные RPC в удобочитаемой форме.

rpctag = map.geotiff.RPCCoefficientTag;
rpctag.LineOffset = 1;
rpctag.SampleOffset = 1;
rpctag.LineScale = 2;
rpctag.SampleScale = 2;
rpctag.GeodeticHeightScale = 500;

Запишите изображение, связанный ссылочный объект и RPCCoefficientTag объект к файлу.

geotiffwrite('myfile',A,RA,'RPCCoefficientTag',rpctag)

В этом примере показано, как записать метаданные коэффициента RPC в файл TIFF. В реальном рабочем процессе метаданные коэффициента RPC создаются в соответствии со спецификацией расширения TIFF. В этом примере не показаны особенности создания допустимых метаданных RPC. Для моделирования необработанных метаданных RPC в примере создается образец файла TIFF с метаданными RPC, а затем используется imfinfo для чтения метаданных RPC в необработанной необработанной форме из файла. Затем этот пример записывает необработанные метаданные RPC в файл с помощью geotiffwrite функция.

myimage = zeros(180,360);
latlim = [-90 90];
lonlim = [-180 180];
R = georefcells(latlim,lonlim,size(myimage));

rpctag = map.geotiff.RPCCoefficientTag;
rpctag.LineOffset = 1;
rpctag.SampleOffset = 1;
rpctag.LineScale = 2;
rpctag.SampleScale = 2;
rpctag.GeodeticHeightScale = 500;

geotiffwrite('myfile',myimage,R,'RPCCoefficientTag',rpctag)

info = imfinfo('myfile.tif');
info.UnknownTags

ans = struct with fields:
        ID: 50844
    Offset: 10676
     Value: [1x92 double]

value = info.UnknownTags.Value;

rpcdata = map.geotiff.RPCCoefficientTag(value) 

rpcdata = 
  RPCCoefficientTag with properties:

                BiasErrorInMeters: -1
              RandomErrorInMeters: -1
                       LineOffset: 1
                     SampleOffset: 1
           GeodeticLatitudeOffset: 0
          GeodeticLongitudeOffset: 0
             GeodeticHeightOffset: 0
                        LineScale: 2
                      SampleScale: 2
            GeodeticLatitudeScale: 1
           GeodeticLongitudeScale: 1
              GeodeticHeightScale: 500
        LineNumeratorCoefficients: [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]
      LineDenominatorCoefficients: [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]
      SampleNumeratorCoefficients: [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]
    SampleDenominatorCoefficients: [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]

geotiffwrite('myfile2',myimage,R,'RPCCoefficientTag',rpcdata)

ginfo = geotiffinfo('myfile2');
ginfo.GeoTIFFTags.RPCCoefficientTag

ans = 
  RPCCoefficientTag with properties:

                BiasErrorInMeters: -1
              RandomErrorInMeters: -1
                       LineOffset: 1
                     SampleOffset: 1
           GeodeticLatitudeOffset: 0
          GeodeticLongitudeOffset: 0
             GeodeticHeightOffset: 0
                        LineScale: 2
                      SampleScale: 2
            GeodeticLatitudeScale: 1
           GeodeticLongitudeScale: 1
              GeodeticHeightScale: 500
        LineNumeratorCoefficients: [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]
      LineDenominatorCoefficients: [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]
      SampleNumeratorCoefficients: [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]
    SampleDenominatorCoefficients: [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]