Запишите файл GeoTIFF
geotiffwrite(___, пишет сетка изображения или данных с одним или несколькими Name,Value)Name,Value парные аргументы, которые управляют различными характеристиками выходного файла.
Считайте изображение 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 в рабочую область.
infilename = 'boston.tif';
[A,R] = geotiffread(infilename);Получите метаданные из файла. Вы передаете это в новую версию BigTIFF файла.
info = geotiffinfo(infilename);
Извлеките значение GeoKeyDirectoryTag поле от info структура возвращена geotiffinfo. Установите значение двух тегов TIFF, которые задают длину и ширину мозаик.
geoTags = info.GeoTIFFTags.GeoKeyDirectoryTag; tiffTags = struct('TileLength',1024,'TileWidth',1024);
Запишите данные в новый файл GeoTIFF, с помощью параметра TiffType, чтобы задать формат файла как 'bigtiff'. Включайте geoTags информация и tiffTags в метаданных в новом файле.
outfilename = 'boston_bigtiff.tif'; geotiffwrite(outfilename,A,R,'GeoKeyDirectoryTag',geoTags, ... 'TiffType','bigtiff','TiffTags',tiffTags)
Считайте данные из сервера WMS.
nasaLayers = wmsfind('nasa', 'SearchField', 'serverurl'); layerName = 'bluemarbleng'; layer = refine(nasaLayers, layerName, 'SearchField', 'layername', ... 'MatchType', 'exact'); [A, R] = wmsread(layer(1));
Запишите данные к файлу GeoTIFF.
filename = [layerName '.tif'];
geotiffwrite(filename, A, R)
Просмотрите данные в файле.
figure
worldmap world
geoshow(filename)
Считайте две смежных ортофотографии и объедините их.
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.
[A, R] = geotiffread('boston.tif');Задайте подмножество данных, которые вы хотите записать.
m = size(A,1);
n = 1024;
firstrow = m-n+1;
lastrow = m;
firstcol = 1;
lastcol = n;
subImage = A(firstrow:lastrow, firstcol:lastcol, :);
xi = [firstcol - .5, lastcol + .5];
yi = [firstrow - .5, lastrow + .5];
[xlimits, ylimits] = intrinsicToWorld(R, xi, yi);
subR = R;
subR.RasterSize = size(subImage);
subR.XLimWorld = sort(xlimits);
subR.YLimWorld = sort(ylimits);
info = geotiffinfo('boston.tif');Запишите подмножество в файл.
filename = 'boston_subimage.tif'; geotiffwrite(filename, subImage, subR, ... 'GeoKeyDirectoryTag', info.GeoTIFFTags.GeoKeyDirectoryTag); figure mapshow(filename);
В этом примере показано, как записать данные о ландшафте SDTS DEM Горы Вашингтон в GeoTIFF. На данные ссылаются к Universal, поперечной меркаторской (UTM), Зона 19, в североамериканской Данной величине 1 927. Это соответствует номеру кода GeoTIFF PCS_NAD27_UTM_zone_19N 26719. Установите растровую интерпретацию на 'регистрации', потому что данными является USGS DEM. Это соответствует растровому типу GeoTIFF PixelIsPoint.
Считайте данные SDTS DEM.
[Z, refmat] = sdtsdemread('9129CATD.ddf'); R = refmatToMapRasterReference(refmat, size(Z), 'postings'); key.GTModelTypeGeoKey = 1; % Projected Coordinate System (PCS) key.GTRasterTypeGeoKey = 2; % PixelIsPoint key.ProjectedCSTypeGeoKey = 26719; filename = '9129.tif'; geotiffwrite(filename, Z, R, 'GeoKeyDirectoryTag', key);
Постройте схему Нью-Гемпшира в UTM.
S = shaperead('usastatelo', 'UseGeoCoords', true, 'Selector',... {@(name) any(strcmp(name,{'New Hampshire'})), 'Name'}); proj = geotiffinfo(filename); [x, y] = projfwd(proj, [S.Lat], [S.Lon]); figure mapshow(x,y)
Отобразите файл GeoTIFF DEM.
hold on h = mapshow(filename, 'DisplayType', 'surface'); demcmap(get(h,'ZData'))
В этом примере показано, как считать файл, содержащий метаданные RPC, измените файл, и затем запишите файл, включая метаданные RPC в новом файле. Чтобы проиллюстрировать этот процесс, пример сначала создает демонстрационный файл TIFF, содержащий Метаданные RPC.
Создайте демонстрационный файл TIFF с метаданными RPC
Создайте демонстрационный тестовый файл TIFF, содержащий метаданные RPC. Для этого тестового файла создайте игрушечное изображение и объект привязки, сопоставленный с изображением.
myimage = zeros(180,360); latlim = [-90 90]; lonlim = [-180 180]; R = georefcells(latlim,lonlim,size(myimage));
Создайте 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',myimage,R,'RPCCoefficientTag',rpctag)
Считайте метаданные и данные изображения из демонстрационного файла TIFF
Считайте метаданные и данные изображения от тестового файла. geotiffinfo возвращает метаданные RPC в поле GeoTIFFTags возвращенной структуры.
info = geotiffinfo('myfile'); [imagedata, R] = geotiffread('myfile');
Запишите модифицированные данные изображения и метаданные RPC к файлу
Сначала измените данные изображения, считанные из файла. Например, установите все пиксели с нулевым знаком на 1 в изображении.
imagedata(imagedata==0)= 1;
Запишите модифицированное изображение в файл, наряду с метаданными RPC из исходного файла.
geotiffwrite('myfile2',imagedata,R,... 'RPCCoefficientTag',info.GeoTIFFTags.RPCCoefficientTag)
В этом примере показано, как записать содействующие метаданные RPC в файл TIFF. В действительном рабочем процессе вы создали бы содействующие метаданные RPC согласно спецификации расширения TIFF. Этот пример не показывает специфические особенности того, как создать допустимые метаданные RPC. Чтобы симулировать необработанные метаданные RPC, пример создает демонстрационный файл TIFF с метаданными RPC и затем использует imfinfo считать эти метаданные RPC в сырых данных, необработанной форме из файла. Пример затем пишет эти необработанные метаданные RPC в файл с помощью geotiffwrite функция.
Создайте необработанные содействующие метаданные RPC
Чтобы симулировать необработанные метаданные RPC, создайте простой тестовый файл и запишите некоторые метаданные RPC в файл. Для этого тестового файла создайте игрушечное изображение и объект привязки, сопоставленный с изображением.
myimage = zeros(180,360); latlim = [-90 90]; lonlim = [-180 180]; R = georefcells(latlim,lonlim,size(myimage));
Создайте объект метаданных 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',myimage,R,'RPCCoefficientTag',rpctag)
Считайте необработанные содействующие метаданные RPC
Считайте содействующие метаданные RPC из тестового файла с помощью imfinfo функция. Когда это сталкивается с незнакомыми метаданными, imfinfo возвращает данные, необработанные, в UnknownTags поле . Отметьте, как поле UnknownTags содержит массив 92, удваивается. Это - необработанные содействующие метаданные RPC, читайте из файла в необработанной форме.
info = imfinfo('myfile.tif');
info.UnknownTagsans = struct with fields:
ID: 50844
Offset: 10672
Value: [1x92 double]
Запишите необработанные метаданные RPC в файл
Запишите необработанные метаданные RPC в файл. Во-первых, извлеките содействующие метаданные RPC из информационной структуры.
value = info.UnknownTags.Value;
Затем создайте объект RPCCoefficientTag, передавая необработанные метаданные RPC (массив 92 удваивается) в качестве аргумента.
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]
Передайте объект RPCCoefficientTag geotiffwrite функционируйте и запишите метаданные RPC в файл.
geotiffwrite('myfile2',myimage,R,'RPCCoefficientTag',rpcdata)
Чтобы проверить, что данные были записаны в файл, считайте метаданные RPC из файла TIFF с помощью geotiffinfo. Сравните возвращенные метаданные RPC с метаданными, записанными в тестовый файл.
ginfo = geotiffinfo('myfile2');
ginfo.GeoTIFFTags.RPCCoefficientTagans =
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]
Если вы работаете с координатами изображений в спроектированной системе координат и R объект растровой привязки карты или матрица привязки, установите GeoKeyDirectoryTag или CoordRefSysCode аргумент, соответственно.
Проверяйте спецификацию GeoTIFF на значения следующих параметров:
'CoordRefSysCode' значение для географических систем координат
'CoordRefSysCode' значение для спроектированных систем координат
Геополе ключа называет для 'GeoKeyDirectoryTag'
Можно также получить 'CoordRefSysCode' значения из файлов данных EPSG (pcs.csv и gcs.csv) в папке:
matlabroot/toolbox/map/mapproj/projdata/epsg_csvRPCCoefficientTag | Tiff | geotiffinfo | geotiffread | imread | imwrite