Εισαγωγή στην Επιστήμη των Η/Υ ΙΙ

Μάθημα 5

Δεδομένα – Ψηφιοποίηση με Matlab

http://seismo.geology.upatras.gr/comp/http://seismo.geology.upatras.gr/comp/

Είδη δεδομένων

• Vector – διανυσματικά δεδομένα

• Raster – “πλεγματικά” δεδομένα

Παραδείγματα

Είδη δεδομένων

• Τα διανυσματικά δεδομένα προκύπτουν από την ψηφιοποίηση χαρτών π.χ. όρια γεωλογικών σχηματισμών, ρήγματα, ποτάμια, ακτογραμμές κλπ

• Τα πλεγματικά δεδομένα προκύπτουν απευθείας π.χ. από δορυφόρους, κάμερες κλπ αλλά συνήθως από παρεμβολή (contouring) και δημιουργία κανάβου από τιμές δειγμάτων

Είδη δεδομένων

Είδη δεδομένων

• Δεδομένα ακτογραμμών (διανυσματικά)

– GSHHS (Global Self-consistent Hierarchical High resolution shoreline data)

– World Data Bank II or CIA Data Bank (WDB) – World Vector Shoreline (WVS)

Είδη δεδομένων

• Παράδειγμα

nan nan40.569739 13.91555940.569152 13.92054640.566805 13.92465340.565632 13.93022740.566512 13.93638740.572086 13.93756040.576486 13.93462740.578833 13.93022740.580886 13.92553340.582646 13.92054640.582646 13.91526640.577073 13.91379940.572379 13.91555940.569739 13.915559nan nan40.457677 13.92700040.456797 13.92964040.452984 13.933454

Είδη δεδομένων

Είδη δεδομένων

Πλεγματικά δεδομένα τοπογραφίας

• ETOPO5 – 5 min ανάλυση

Πλεγματικά δεδομένα τοπογραφίας

• ETOPO2 – 2 min ανάλυση

http://www.gfdl.noaa.gov/products/vis/data/datasets/etopo2_topography.html

Πλεγματικά δεδομένα τοπογραφίας

• GTOPO30 – 30 sec ανάλυση ~ 1km

http://edc.usgs.gov/products/elevation/gtopo30/gtopo30.html

Πλεγματικά δεδομένα τοπογραφίας• SRTM – 30 και 90m ανάλυση

http://www2.jpl.nasa.gov/srtm/

Πλεγματικά δεδομένα τοπογραφίας• SRTM – 30 και 90m ανάλυση

Πλεγματικά δεδομένα τοπογραφίας

• Το μέγεθος των αρχείων μεγαλώνει όσο αυξάνει η ανάλυση

• ETOPO5 ~ 20Mb

• SRTM ~ 3Mb 1x1 ~ 200Gb

Προβολές

Προβολές

> 60 προβολές (εντολή axesm)

help mapproj

Matlab Mapping toolbox

• Worldmap

• figure• worldmap(‘greece','patch')• scaleruler• hidem(gca)

Matlab Mapping toolbox

• Plotm φτιάχνει το χάρτη

• Inputm

• textm

worldmap greece[lat,lon] = inputm(4);plotm(lat,lon,'-+r')hidem(gca)

[lat,long] = inputm(npts)

Matlab Mapping toolbox

load geoidworldmap(geoid,geoidlegend)contourcmap(10,'jet','colorbar','on','location','horizontal')

Matlab Mapping toolbox

• meshmload topoaxesm millermeshm(topo,topolegend,[90 180])

Demos

• viewmaps

Ψηφιοποίηση με Matlab• Ψηφιοποίηση στην οθόνη

• Ψηφιοποίηση με digitizer

Ψηφιοποίηση με Matlab

• Εισαγωγή της εικόνας στη Matlab

• Δημιουργία ενός συστήματος συντεταγμένων

• Εισαγωγή δεδομένων

Ψηφιοποίηση με Matlab• Imread

• A = imread(filename,fmt)

Η συνάρτηση ginput

• [x,y] = ginput(n)

• Επιστρέφει συντεταγμένες σημείων από ένα γράφημα

Η συνάρτηση minputfunction data = minput(imagefile)% Specify the limits of the imagexmin = input('Specify xmin! ');xmax = input('Specify xmax! ');ymin = input('Specify ymin! ');ymax = input('Specify ymax! ');

% Read image and displayB = imread(imagefile);a = size(B,2); b = size(B,1);imshow(B);

% Define upper left and lower right corner of imagedisp('Click on lower left and upper right cr, then <return>')[xcr,ycr]= ginput;XMIN=xmin-((xmax-xmin)*xcr(1,1)/(xcr(2,1)-xcr(1,1)));XMAX=xmax+((xmax-xmin)*(a-xcr(2,1))/(xcr(2,1)-xcr(1,1)));YMIN=ymin-((ymax-ymin)*ycr(1,1)/(ycr(2,1)-ycr(1,1)));YMAX=ymax+((ymax-ymin)*(b-ycr(2,1))/(ycr(2,1)-ycr(1,1)));

% Digitize data pointsdisp('Click on data points to digitize, then <return>')[xdata,ydata]= ginput;XDATA = XMIN + ((XMAX-XMIN)*xdata / size(B,2));YDATA = YMIN + ((YMAX-YMIN)*ydata / size(B,1));data(:,1) = XDATA; data(:,2) = YDATA;