Matlab - Pasar de Matlab a Arcgis

   
Vista:

Pasar de Matlab a Arcgis

Publicado por mariel (4 intervenciones) el 20/09/2012 20:48:33
Hola a todos

Tengo una duda y espero que alguno pueda ayudarme.

Tengo una serie de archivos de gps que son rutas, a estas les extraje la información que necesitaba y la guarde a parte. Ahora necesito gráficar todas mis rutas juntas (en una sola imagen) hasta aquí todo bien, el problema que tengo es que necesito pasar esto a ArcGis pero no me sirve que sea solo una imagen por que necesito usar la información de estos archivos para crear shapefiles en arcgis.

No se si me explique bien, hay algun modo de hacer esto?? me habian dicho de usar 'geotiffwrite' estuve investigando un poco pero la verdad no entiendo muy bien como usarlo.

Si alguien supiera como hacer esto se los agradeceria mucho, la verdad es que no tengo mucha experiencia con matlab y esto me esta volviendo un poco loca.

gracias
Valora esta pregunta
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder
Imágen de perfil de JOSE JEREMIAS CABALLERO

Pasar de Matlab a Arcgis

Publicado por JOSE JEREMIAS CABALLERO jjcc94@hotmail.com (3440 intervenciones) el 20/09/2012 21:23:17
Mis conocimientos de arcgis son basicos.
Aqui tienes un codigo en matlab sobre mapas y el uso del comando de matlab que mencionas
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
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'))



Saludos.
JOSE JEREMIAS CABALLERO
Asesor de Proyectos con Matlab
programador en matlab
Servicios de programacion Matlab
jjcc94@hotmail.com

http://matlabcaballero.blogspot.com
Valora esta respuesta
Me gusta: Está respuesta es útil y esta claraNo me gusta: Está respuesta no esta clara o no es útil
0
Comentar

Pasar de Matlab a Arcgis

Publicado por mariel (4 intervenciones) el 21/09/2012 16:50:30
Hola Jose

Muchisimas gracias, voy a revisarlo y ver si logro solucionar mi problema
Valora esta respuesta
Me gusta: Está respuesta es útil y esta claraNo me gusta: Está respuesta no esta clara o no es útil
0
Comentar