Чтение изображения 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 Проективную систему координат.

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 объект метаданных и установите некоторые поля с типичными значениями. The 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]