Считайте изображение 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');

Обрежьте данные к пределам, заданным xlim и ylim использование 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 поле как 4 326.

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: 10672
     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]