Εισαγωγή στην Επιστήμη των Η/Υ ΙΙ Μάθημα 5 Δεδομένα – Ψηφιοποίηση με Matlab

Εισαγωγή στην Επιστήμη των Η/Υ ΙΙ Μάθημα 5 Δεδομένα – Ψηφιοποίηση με Matlab 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 nan 40.569739 13.915559 40.569152 13.920546 40.566805 13.924653 40.565632 13.930227 40.566512 13.936387 40.572086 13.937560 40.576486 13.934627 40.578833 13.930227 40.580886 13.925533 40.582646 13.920546 40.582646 13.915266 40.577073 13.913799 40.572379 13.915559 40.569739 13.915559 nan nan 40.457677 13.927000 40.456797 13.929640 40.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 [lat,long] = inputm(npts) worldmap greece [lat,lon] = inputm(4); plotm(lat,lon,'-+r') hidem(gca)

Matlab Mapping toolbox load geoid worldmap(geoid,geoidlegend) contourcmap(10,'jet','colorbar','on','location','horizontal')

Matlab Mapping toolbox • meshm load topo axesm miller meshm(topo,topolegend,[90 180])

Demos • viewmaps

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

Ψηφιοποίηση με Matlab • Εισαγωγή της εικόνας στη Matlab • Δημιουργία ενός συστήματος συντεταγμένων • Εισαγωγή δεδομένων

Ψηφιοποίηση με Matlab • Imread • A = imread(filename,fmt)

Η συνάρτηση ginput • [x,y] = ginput(n) • Επιστρέφει συντεταγμένες σημείων από ένα γράφημα

Η συνάρτηση minput function data = minput(imagefile) % Specify the limits of the image xmin = input('Specify xmin! '); xmax = input('Specify xmax! '); ymin = input('Specify ymin! '); ymax = input('Specify ymax! '); % Read image and display B = imread(imagefile); a = size(B,2); b = size(B,1); imshow(B); % Define upper left and lower right corner of image disp('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 points disp('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;

Εισαγωγή στην Επιστήμη των Η/Υ ΙΙ Μάθημα 5 Δεδομένα – Ψηφιοποίηση με Matlab