Matlab - Pasar de Matlab a Arcgis

function matlab_arcgis

close all

clear all

basename = 'boston_ovr';

imagefile = [basename '.jpg'];

RGB = imread(imagefile);

worldfile = getworldfilename(imagefile);

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

filename = [basename '.tif'];

geotiffwrite(filename, RGB, R)

figure(1)

usamap(RGB, R)

geoshow(filename)

clear all

nasa = wmsfind('nasa', 'SearchField', 'serverurl');

layerName = 'bluemarbleng';

layer = nasa.refine(layerName,  'SearchField', 'layername', ...

  'MatchType', 'exact');

[A, R] = wmsread(layer(1));

filename = [layerName '.tif'];

geotiffwrite(filename, A, R)

figure(2)

worldmap world

geoshow(filename)

clear all

% Read the two adjacent orthophotos and combine them. 

X_west = imread('concord_ortho_w.tif');

X_east = imread('concord_ortho_e.tif');

X = [X_west X_east];

% Construct referencing objects for the orthophotos and for their

% combination.

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);

% Write the combined image to a GeoTIFF file. Use the code number, 

% 26986, indicating the PCS_NAD83_Massachusetts Projected Coordinate 

% System.

coordRefSysCode = 26986;

filename = 'concord_ortho.tif';

geotiffwrite(filename, X, R, 'CoordRefSysCode', coordRefSysCode);

figure(3)

mapshow(filename)

clear all

[A, R] = geotiffread('boston.tif');

row = [size(A,1)-1024+1 size(A,1)];

col = [1 1024];

subImage = A(row(1):row(2), col(1):col(2), :);

xi = col + [-.5 .5];

yi = row + [-.5 .5];

[xlim, ylim] = R.intrinsicToWorld(xi, yi);

subR = R;

subR.RasterSize = size(subImage);

subR.XLimWorld = sort(xlim);

subR.YLimWorld = sort(ylim);

info = geotiffinfo('boston.tif');

filename = 'boston_subimage.tif';

geotiffwrite(filename, subImage, subR,  ...

'GeoKeyDirectoryTag', info.GeoTIFFTags.GeoKeyDirectoryTag);

figure(4)

mapshow(filename);

clear all

[Z, refmat] = sdtsdemread('9129CATD.ddf');

R = refmatToMapRasterReference(refmat, size(Z));

R.RasterInterpretation = 'postings';

key.GTModelTypeGeoKey  = 1;  % Projected Coordinate System (PCS)

key.GTRasterTypeGeoKey = 2;  % PixelIsPoint

key.ProjectedCSTypeGeoKey = 26719;

filename = '9129.tif';

geotiffwrite(filename, Z, R, 'GeoKeyDirectoryTag', key);

% Plot the outline of the state of New Hampshire in 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(5)

mapshow(x,y)

% Display the GeoTIFF DEM file.

hold on

h = mapshow(filename, 'DisplayType', 'surface');

demcmap(get(h,'ZData'))