LOK-Viewer : Graphical Interface Application for Geomagnetic Data in Geomagnetic Lombok Observatory

Size: px

Start display at page:

Download "LOK-Viewer : Graphical Interface Application for Geomagnetic Data in Geomagnetic Lombok Observatory"

Trevor Rich
5 years ago
Views:

1 LOK-Viewe : Gaphical Iteface Applicatio fo Geogetic Data i Geogetic Lombok Obsevatoy Gii W. Wiiasto (1) Electical Egieeig Dept. Mataam Uivesity Lombok, Idoesia Teti Zubaidah (2) Electical Egieeig Dept. Mataam Uivesity Lombok, Idoesia Bulkis Kaata (3) Electical Egieeig Dept. Mataam Uivesity Lombok, Idoesia Idia P. Kiasih (4) Mathetics Educatio Dept. IKIP Mataam Lombok, Idoesia Abstact This pape is aimed to discuss the developmet pocess of a iteactive applicatio which able to display umeical, gaphical ad statistical pofile of geogetic data fo Lombok Obsevatoy, amed LOK-Viewe. LOK which is a codeame fo Geogetic Lombok Obsevatoy, is beig focused to obtai ad measue data o the eath s getic field. LOK has take that measuemet sice Mach It is impotat to have a dyamic vesio of this data by visualize ad display it cotiuously as a itegated system. This softwae is specific, because it is adapted fom aw data files. It is mea that it is ot associated with ay paticula pogam o applicatio o opeatig systems.. Mea value, sllest value ad lagest value ad stadad deviatio ae showed i statistical sumy pael. This softwae iteface is developed usig GUI-Matlab ad compiled build code with executable file usig Matlab Compile Rutime 2012b. LOK-Viewe Ve.2015a has bee developed ad able to be u diectly stadaloe. This applicatio allows uses to coduct visual obsevatio by oticig eath s getic field gaphical patte show o the iteface with diffeet daily data. It is also ceate a possibility of easiess ad simplicity fo uses, sice they do ot eed to collect a vey lage data ually i ode to obseve its chaacteistics. getism i ay place, we must measue the diectio ad itesity of the field. The Eath's getic field is descibed by seve paametes. These ae decliatio (D), icliatio (I), hoizotal itesity (H), the oth (X) ad east (Y) compoets of the hoizotal itesity, vetical itesity (Z), ad total itesity (F). Paametes that descibe diectio of the getic field ae decliatio ad icliatio. D ad I ae measued i degees, positive east fo D ad positive dow fo I. The itesity of the total field is descibed by H, Z,, X ad Y compoets of the hoizotal itesity [3]. LOK Obsevatoy use Thee-axial fluxgate getomete FGE susped device to measue the vaiatio of getic field i thee vecto compoets.. This tool geeates log data as follows; Keywods LOK-Viewe; Lombok Geogetic Obsevatoy; Geogetic data I. INTRODUCTION Geogetic Obsevatoy Lombok (LOK*) located o Lombok Islad at Latitude: S; Logitude: W that efe to IGRF-11 fo The LOK Obsevatoy is u i coopeatio of the Faculty of Egieeig, Mataam Uivesity ad the GFZ-Postdam sice Mach The geogetic field itesity ad diectios ae cotiously ecoded at the obsevatoy [1]. Lombok Islad has the highest itesity geogetic aolies i Idoesia [2]. Theefoe, its existece will give a cucial cotibutio i tems of fillig the data void fom the global geogetic equato egio, especially Southeast Asia egio. Cotiuous ad log-tem obsevatio to geogetic pheomeo is eeded to be able to udestad some chages i Eath's getic field, which is a life shields agaist sola stoms ad cosmic ays. Geogetic data ae vey useful fo y pactical puposes, icludig ai ad sea avigatio, the seach fo atual esouces, especially mieals ad geothel esouces, as well as fo the mitigatio of atual disastes, icludig eathquake.to measue the Eath's Figue 1: Show fist colum ked with oage dash box thee is six vaiabel fom aw data icludig Time, X_value, Y_value, Z_value, F_value ad delta-f_value. 2d colum which is ked by box k ed dash peset the aw data value. The ame file of data has followed the amig stadat, i.e lo0_ dat. lo0 is Lombok Obsevatoy ad meas local time data poduce with fot yy-mm-dd. Sice we have to deal with a vey lage data size ad we also have to obseve the data i days, we foud that it will be impotat to desig a itegated system which allow us to coduct a futhe data exploig without beig botheed by data iputtig activities. The system should be iteactive, dyamic ad pactical. It should be able to big out statistical sumy ad show a visual display such as gaphic o cetai patte that ifom uses about data chaacteistics ad facilitate them to ife ad detemie the subsequet decisio elated to geogetic obsevatio. At this stage, we develop a applicatio called LOK-Viewe Ve.2015b. It is employed as a gaphic use iteface ad able

to ead aw data by displayig multiple pespectives such as daily time seies of aw data, the aveage of data daily, weekly, mothly util the data ove the whole oe yea.

Figue is compoet o object is called Hadle Gaphics. Each compoet o object o the scee has a uique idetifie, called a hadle. That Allows us to go back ad modify the object at ay time.

2 to ead aw data by displayig multiple pespectives such as daily time seies of aw data, the aveage of data daily, weekly, mothly util the data ove the whole oe yea. This applicatio is developed usig a base Gaphical Use Iteface (GUI) Matlab 2012b. II. APPLICATION DEVELOPMENT APPROACH A. GUI Matlab i a GUI Matlab oe that easily ecogizable as figue. Figue is compoet o object is called Hadle Gaphics. Each compoet o object o the scee has a uique idetifie, called a hadle. That Allows us to go back ad modify the object at ay time. Figue 3: sample GUI Layout B. Flowchat Model I geeal, the algoithm applied to the GUI / iteface Mai widow (page 1) applicatio was desiged usig a flowchat diagam show i figue 2. Based o this algoithm,we wite a souce code o pogam that fill i Callback Fuctio. I this case, Callback Fuctio is a fuctio of object compoet 'Butto 1b (see figue 3) amed by push_ambildata_callback. We iclude The complete pogam late i aothe sectio of this pape. Butto 1b Click Geeally thee s o eed to woy about the hadles, but they e always thee. Show below is a diagam of the Hadle Gaphics object hieachy. These ae all the categoies of object types i MATLAB. The Popety Edito, oe of the Guide tools, is extemely useful fo avigatig the Hadle Gaphics hieachy. If you eed to get o set ay Hadle Gaphics popety o ay object, chaces ae you ca do it most easily with the Popety Edito [4]. Figue 2 show GUI desig pocess appoach. Tedapat dua hal yag uta yaki desig ad wite code. Get file fom existig folde with ( *.dat) file extetio Pasig with space delimete fom ( *.dat) file value Impot ( *.dat) file value o to Vaiable amed data Fid dimesio tix (mx) size vaiabel of data N If values of data tix (mx) >= ? Y Replace value data(i,j) = NaN Check value (i,j) util (m-1, -1) = tue N Y Filte outlies value of data Pasig value of data based o tix coloum ad stoe to ew vaiabel Figue 2: GUI Desig pocess As show i the illustatio figue 3, eve afte you ve used Guide to get the layout ight, you still eed to wite the MATLAB laguage callbacks. Techical oveview about the elatio betwee GUI figue ad callback fuctio o LOK-Viewe ca be see o figue 6 New vaiabel ( Data_X, Data_Y, Data_Z, Data_F, data_deltaf ) & Plottig Figue 4. Flowchat diagam fo Mai Widow iteface fom LOK- Viewe Figue 4 shows flowchat that daw a oveall idea of the i widow. The secod page s algoithm is ot too diffeet fom the oe i i widow, it just eeds to be adjusted slightly. Such as adjustig the comd to specify data hadlig (see ectagle secod box i figue 4).

A. LOK-Viewe Desig Itefaces Figue 5 shows LOK-Viewe s i widow iteface. It is cotai load data butto, exit butto ad diamyc text box that idicate the existig file ame.

3 A. LOK-Viewe Desig Itefaces Figue 5 shows LOK-Viewe s i widow iteface. It is cotai load data butto, exit butto ad diamyc text box that idicate the existig file ame. Thee ae also fou axes plot which able to display gaphics daw fom oigial data, III. DESIGN INTERFACES assisted by adio butto that allow uses to adjust what gaphic will be displayed. Table display o the left side of i widow will show umeical data afte it is loaded with load data butto. Mai widow also had a statistical sumy table o the left-bottom side. It will ease use to gai a statistical ifeece elated to geogetic data 2 1b 10b 9 3b 7b 4b 5b 6b 8b Figue 5: The i widow of Display Applicatio Use Iteface LOK-Viewe. Thee ae 10 boxes labeled by 1b,2,3b,4b,5b,6b,7b,8b,9,10b. 1b) Butto fo load *.DAT data; 2) Text compoet to show ame of data; 3b)Axes1 to show vaiabel X_value ; 4b) Axes2 to show vaiabel Y_value 5b) Axes3 to show vaiabel Z_value ; 6b) Axes4 to show vaiabel delta-f_value ; 7b) Tabledata1 compoet to show *.DAT umeical aw data; 8b) Tabledata2 compoet to show *.DAT statistical fom aw data B. LOK-Viewe GUI Compositio Figue 6 shows the compositio block betwee objects/compoets i i widow iteface as show o figue 5. Mai Fuctio : Object/Compoet 1b Object/Compoet 3b,4b,5b,6b Object/Compoet 7b,8b Callback fuctio 1b Callback fuctio 1b Callback fuctio 1b Object/Compoet 2 Callback fuctio 2 Figue 6 : Diagam Block LOK-Viewe GUI Compositio This GUI Compositio shows that thee ae two i callback fuctio iside the LOK-Viewe i fuctio. These two i callback fuctio will take cotol all of the compoets ad objects C. LOK-Viewe GUI Matlab code This sessio will show us the i widow LOK-Viewe s souce code as displayed o figue 5. These codes allow us to fomulate data pocessig. We also eclose some detail commets that help eade to udestad this pogam GUI-CallBack Fuctio Scipt Mai Widow (Page-1) % Last Modified by GUIDE v Ju :50:31 % Mai_Fuctio % fuctio vaagout = LOKViewe_V03(vaagi) % LOKVIEWER_V03 M-file fo LOKViewe_V03.fig % Begi iitializatio code - DO NOT EDIT gui_sigleto = 1; gui_state = stuct('gui_name', mfileame,... 'gui_sigleto', gui_sigleto,... 'gui_layoutfc', [],... 'gui_callback', []); if agi && ischa(vaagi{1}) gui_state.gui_callback = st2fuc(vaagi{1}); if agout [vaagout{1:agout}] = gui_ifc(gui_state, vaagi{:}); else

4 gui_ifc(gui_state, vaagi{:}); % Ed iitializatio code - DO NOT EDIT OF MAIN FUNCTION % %-----Iitially Sub-Fuctio fom Compoet Object % % --- Executes just befoe LOKViewe_V03 is de visible. fuctio LOKViewe_V03_OpeigFc(hObject, evetdata, hadles, vaagi) % Choose default comd lie output fo LOKViewe_V03 hadles.output = hobject; % Update hadles stuctue guidata(hobject, hadles); % -- Outputs fom this fuctio ae etued to the comd lie. fuctio vaagout = LOKViewe_V03_OutputFc(hObject, evetdata, hadles) % Get default comd lie output fom hadles stuctue vaagout{1} = hadles.output; % -- Executes o butto pess i push_ambildata. (Get file fom exisstig folde see block diagam) fuctio push_ambildata_callback(hobject, evetdata, hadles) %---This is block pat of Read / impot File fileame1=0; %Reset value file %-- Opeig dialog box to gettig file as vaiabel fileame1 fileame1 = uigetfile('*.dat'); %--- Used fo exeptio hadle dialog box if isequal(fileame1,0) etu %--- Used fo gettig fileame1 to text1 box set(hadles.text1,'stig',fileame1); %--- Used fo ead *.DAT file ad stoe to data vaiabel ' delimitei = ' '; %Pasig text with Delimete space headeliesi =0; A = impotdata(fileame1,delimitei,headeliesi); data = A; %Used fo eame ew vaiabel %gettig Matix Dimesio fo data vaiabel SizeData=size(data); Nows=SizeData(1); Ncols=SizeData(2); %filteig Outlie data fo i=1:nows fo j=1:ncols if data(i,j)>= ; data (i,j)=nan; %eplace data with NaN % % used fo Pasig data vaiabel based o Colum Elemet ad stoe to ew vaiabel datax=data(:,1); % Vecto Vaiabel fo time / X-coodiate datay1=data(:,2); % Vecto Vaiabel fo Y1 (data_x / Y1-coodiate) datay2=data(:,3); % Vecto Vaiabel fo Y2 (data_y / Y2-coodiate) datay3=data(:,4); % Vecto Vaiabel fo Y3 (data_z / Y3-coodiate) datay4=data(:,5); % Vecto Vaiabel fo Y4 (data_f/ Y4-coodiate) datay5=data(:,6); % Vecto Vaiabel fo Y5 (data_deltaf / Y5-cood) % STATISTICS BLOCK [Ndata] = [size(datay1,1) size(datay2,1) size(datay3,1) size(datay4,1) size(datay5,1)] % Numbe of data ows [MINdata] = [mi(datay1) mi(datay2) mi(datay3) mi(datay4) mi(datay5)]; % mi data value [MAXdata] = [x(datay1) x(datay2) x(datay3) x(datay4) x(datay5)]; % x data value [MEANdata] = [mea(datay1) mea(datay2) mea(datay3) mea(datay4) mea(datay5)]; )] % mea data value [STDdata] = [std(datay1) std(datay2) std(datay3) std(datay4) std(datay5)]; )]; % std.dev data value stat = cat(1, Ndata, MINdata, MAXdata, MEANdata, STDdata); % Plot data x =datax; % X coodiate y1=datay1; % Y1 o X_value coodiate y2=datay2; % Y2 o Y_value coodiate y3=datay3; % Y3 o Z_value coodiate y5=datay5; % Y5 o delta-f_value coodiate % Plottig X_value o axes1 plot(hadles.axes1,x,y1,'-s','liewidth',2,... 'MakeEdgeColo','k',... 'MakeFaceColo','g',... 'MakeSize',1) %Settig gaphical layout set(hadles.axes1,'xgid','o','ygid','o','xtick',0:23) xlabel(hadles.axes1,'time(h)'); ylabel(hadles.axes1,'(x/dt)t'); % Plottig Y_value o axes2 plot(hadles.axes2,x,y2,'-s','liewidth',2,... 'MakeEdgeColo','k',... 'MakeFaceColo','g',... 'MakeSize',1) set(hadles.axes2,'xgid','o','ygid','o','xtick',0:23) xlabel(hadles.axes2,'time(h)'); ylabel(hadles.axes2,'(y/dt)t'); % Plottig Z_value o axes3 plot(hadles.axes3,x,y3,'-s','liewidth',2,... 'MakeEdgeColo','k',... 'MakeFaceColo','g',... 'MakeSize',1) set(hadles.axes3,'xgid','o','ygid','o','xtick',0:23) xlabel(hadles.axes3,'time(h)'); ylabel(hadles.axes3,'(z/dt)t'); % Plottig delta-f_value o axes3 plot(hadles.axes4,x,y5,'-s','liewidth',1,... 'MakeEdgeColo','k',... 'MakeFaceColo','g',... 'MakeSize',0.5) set(hadles.axes4,'xgid','o','ygid','o','xtick',0:23) xlabel(hadles.axes4,'time(h)'); ylabel(hadles.axes4,'(delta-f/dt)t'); %Callig data vaiabel to display o Table Compoet object myfom = guidata(gcbo); datatabel = data; set(myfom.tbldata,'data',datatabel); % Callig data vaiabel to display o Table Statistic Comp.object myfom2 = guidata(gcbo); datatabel2 = stat; set(myfom2.tblstat,'data',datatabel2); % --- Executes o butto pess i Exit. fuctio Kelua_Callback(hObject, evetdata, hadles) % ---- Customized as follows ---- close(acesto(hobject,'figue')) % END_of_Mai&Sub_Fuctio %

5 IV. RESULTS A. LOK-Viewe GUI Iteface ad Data Repesetatio We y otice fom the uig esult o figue 7, that LOK- Viewe is able to pefom its ole as aw data loade ad also display each gaphic o i widow gaphical iteface (see box A), umeical table (see box B ad box C). The displayed data wee take fom oe sample with attibute lo0_ dat with Time as x-axes fom to B A C Figue 7 : The i widow : display applicatio Iteface of i widow (Page-1) cosist of geogetic daily data. Colou leg compoet object data elatio betwee data Plot (box A ) ad data Table Compoet (box B, C ) object: (ed dash lie) : vaiabel (X/dt)T (blue dash lie) : vaiabel (Z/dt)T (gee dash lie) : vaiabel (Y/dt)T (black dash lie) : vaiabel (delta-f/dt)t B. Readig Geogetic Sigal It appeas that the data o the day of Oct ' is complete ad thee is o loge a outlie value o missig value. The gaph shows geogetic sigal fom to 23:59 local time. Based o aothe elated studies, sigal samplig ca be doe based o the best data (data at cetai times) with miimum oise. We ca take a look this sigal close to be able to obseve the Ulta Low Fequecy (ULF) data pheomeo. C. Statistical Result I the left-bottom side of figue 7 we ca obseve a statistical sumy fo each data seies. It is impotat fo uses to have this ifotio so that they ca gathe some statistical ifeece which suppot data aalysis. Still associated with the aw data o figue 1, to kow thee ae some outlies the data ecoded by the seso device due to exteal factos. We ecouteed a data outlie woth ad hecefoth we filte the data of the outlie. 'N' is the total time seies of data; 'Nmi' is the Miimum Value of evey elemet geogetic vaiatio; 'Nx' is a ximum value of evey elemet geogetic vaiatio; 'Mea' is the mea value of evey elemet geogetic vaiatio; 'Std-Dev' is the value of evey elemet Std.Deviatio geogetic vaiatio; Fom these data the followig values i a ow : each N = 1440; Nmi fo X(T) = e+04; Nmi fo Y(T) = ; Nmi fo Z(T) = e+04; Nmi fo F(T) = e+04; delta-f(t) = ; Nx fo X(T) = e+04; Nx fo Y(T) = ; Nx fo Z(T) = e+04; Nx fo F(T) = e+04; Nx fo delta- F(T) = ; Mea fo X(T) = e+04; Mea fo Y(T) = ; Mea fo Z(T) = e+04; Mea fo F(T) = e+04; Mea fo delta-f(t) = ; Std.Dev fo X(T) = ; Std.Dev fo Y(T) = ; Std.Dev fo Z(T) = ; Std.Dev fo F(T) = ; Std.Dev fo delta-f(t) = ;

6 D. Data Repesetatio (LOK-Viewe Page 2) Figue 8 shows gaphical display fo LOK Obsevatoy data thoughout the Jauay to Decembe Thee ae ed sigal i Apil, August, Septembe, Octobe ad Novembe. Thee is also a empty sigal, which meas that the seso has ot bee able ecod the daily data caused by exteal factos. We still cotiue to wok so that the sigal ca be ecoded popely so that the missig data ca be miimized. 1a 1b 1c Figue 8 : Expose of geogetic sigal vaiatio of vaiabel X ( T), base lie gaphical fom jauay to decembe i yea 2014 with diffeet view (1a),(1b),(1c) whee (1a) is oigi (aw) data ; (1b) is daily mea fom oigi data ; (1c) is daily mea with scatte view. Fo complete sigal vaiatio of Vaiabel Y(T), Z(T) ad deltaf(t) ca be see o Appix 1. V. CONCLUSION AND DISCUSSION This softwae iteface is developed usig GUI-Matlab ad compiled build code with executable (*.exe) file usig Matlab Compile Rutime 2012b. LOK-Viewe Ve.2015a has bee developed ad able to be u diectly stadaloe. This applicatio allows uses to coduct visual obsevatio by oticig eath s getic field gaphical patte show o the iteface with diffeet daily data. It is also ceate a possibility of easiess ad simplicity fo uses, sice they do ot eed to collect a vey lage data ually i ode to obseve its chaacteistics. This applicatio will be developed futhe by addig some featues fo digital sigal pocessig.

7 Ackowledgemet This eseach was suppoted by The Applied Electogetic Techology Reseach Goup, Uivesity of Mataam. It is pat of Uggula Peguua Tiggi Reseach fuded by DIKTI o [1] Last Access o Sept [2] Zubaidah, Teti et.al. [3] Access o Sept Last [4] The Laguage Decembe of techical Computig,MATLAB Refeeces APPENDIX 1: A. Gaphical Iteface to shows oe yeas (jauay decembe 2014) Geogetic data at LOK with diffeet view 2a ja peb ap y ju jul aug se p ok t ov de c 2b ja peb ap y ju jul aug se p ok t ov de c 2c ja peb ap y ju jul aug se p ok t ov de c Figue 9 : Expose of geogetic sigal vaiatio Y(T) with difeet view (2a,2b,2c)

8 APPENDIX 1 (cot d) : B. Gaphical Iteface to shows oe yeas (jauay decembe 2014) Geogetic data at LOK with diffeet view 3a 3b 3c Figue 10 : Expose of geogetic sigal vaiatio Z(T) with difeet view (3a,3b,3c)

9 APPENDIX 1 (cot d) : C. Gaphical Iteface to shows oe yeas (jauay decembe 2014) Geogetic data at LOK with diffeet view 4a 4b 4c Figue 11 : Expose of geogetic sigal vaiatio delta-f(t) with difeet view (4a,4b,4c)

LOK-Viewer : Graphical Interface Application for Geomagnetic Data in Observatory Geomagnetic Lombok

LOK-Viewer : Graphical Interface Application for Geomagnetic Data in Observatory Geomagnetic Lombok Giri W. Wiriasto, Teti Zubaidah, Bulkis Kanata Electrical Engineering Dept. Mataram University Lombok,