Accuracy Assessment and Comparative Analysis of IDW, Spline and Kriging in Spatial Interpolation of Landform (Topography): An Experimental Study

Size: px
Start display at page:

Download "Accuracy Assessment and Comparative Analysis of IDW, Spline and Kriging in Spatial Interpolation of Landform (Topography): An Experimental Study"

Transcription

1 Journal of Geographc Informaton System, 207, 9, ISSN Onlne: ISSN Prnt: Accuracy Assessment and Comparatve Analyss of IDW, Splne and Krgng n Spatal Interpolaton of Landform (Topography): An Expermental Study Maduako Nnamd Ikechukwu *, Eljah Ebnne, Ufot Idorenyn, Ndukwu Ike Raphael Department of Geonformatcs and Surveyng, Unversty of Ngera, Nsukka, Ngera How to cte ths paper: Ikechukwu, M.N., Ebnne, E., Idorenyn, U. and Raphael, N.I. (207) Accuracy Assessment and Comparatve Analyss of IDW, Splne and Krgng n Spatal Interpolaton of Landform (Topography): An Expermental Study. Journal of Geographc Informaton System, 9, Receved: March 28, 207 Accepted: June 20, 207 Publshed: June 23, 207 Copyrght 207 by authors and Scentfc Research Publshng Inc. Ths work s lcensed under the Creatve Commons Attrbuton Internatonal Lcense (CC BY 4.0). Open Access Abstract It s practcally mpossble and unnecessary to obtan spatal-temporal nformaton of any gven contnuous phenomenon at every pont wthn a gven geographc area. The most practcal approach has always been to obtan nformaton about the phenomenon as n many sample ponts as possble wthn the gven geographc area and estmate the values of the unobserved ponts from the values of the observed ponts through spatal nterpolaton. However, t s mportant that users understand that dfferent nterpolaton methods have ther strength and weaknesses on dfferent datasets. It s not correct to generalze that a gven nterpolaton method (e.g. Krgng, Inverse Dstance Weghtng (IDW), Splne etc.) does better than the other wthout takng nto cognzance, the type and nature of the dataset and phenomenon nvolved. In ths paper, we theoretcally, mathematcally and expermentally evaluate the performance of Krgng, IDW and Splne nterpolaton methods respectvely n estmatng unobserved elevaton values and modelng landform. Ths paper undertakes a comparatve analyss based on the predcton mean error, predcton root mean square error and cross valdaton outputs of these nterpolaton methods. Expermental results for each of the method on both based and normalzed data show that Splne provded a better and more accurate nterpolaton wthn the sample space than the IDW and Krgng methods. The choce of an nterpolaton method should be phenomenon and data set structure dependent. Keywords Spatal Interpolaton, IDW, Krgng, Splne and Modelng Elevaton DOI: /jgs June 23, 207

2 . Introducton Interpolaton ams at fndng the values of a functon f ( x ) for an x between x, x,, x at whch the values of f ( x ) are gven. The dfferent x values 0 gven values n ( ), ( ),, ( ) f = f x f = f x f = f x () 0 0 n n can be obtaned from a mathematcal functon or from an emprcal functon modelled from observatons or experments []. Spatal nterpolaton therefore ams at estmatng values of a spatal phenomenon or functon (temperature, elevaton, etc.) at unobserved/estmated ponts, gven values of the phenomenon at observed/estmated ponts. Spatal nterpolaton has contnued to be an mportant tool for estmatng contnuous spatal envronmental varables for effectve decson makng. Many modelng tools ncludng Geographc Informaton System offer the earth and envronmental scentst the ablty to carry out spatal nterpolaton routnely to generate useful spatal contnuous data for all knds of analyss [2]. Interpolaton becomes very useful and essental n scenaros where, the resoluton, orentaton, or cell sze of a dscretzed surface vares from what s needed. It s also employed when contnuous surface s represented by a data model dfferent from what s desred, and when data spread does not cover an area of nterest totally [3]. Spatal nterpolaton methods gve a means of predctng values of an envronmental parameter at unmeasured locaton usng data from pont measurements wthn the sample space [4]. In an deal stuaton, a fnte set of nputs establshes varatons n an envronmental parameter and they exactly conform to an establshed physcal law. If a relatonshp s establshed, the values of the desred parameter can be correctly estmated. The relatonshp between target varable and the physcal envronment cannot be modeled exactly because of ts complexty [5]. Ths s due to a lack of suffcent knowledge of: (a) the complete lst of nputs nto the model (b) the relatonshp needed to determne the output from these nputs and (c) the mportance of the random part of the system. Estmatng a model usng feld measurement of the parameter of nterest then becomes the only way [6]. For a sample of a target varable Z, denoted as z( s), z( s2),, z( s n ) (where s = ( x, y ) s the locaton context, x, y are the coordnates n geographcal space and n s the number of observatons), the geographcal doman of nterest can be denoted as A. Consderaton s gven to only samples z( s n ) that realzes a process Z = { Z( s), s A } [6]. Accordng to Mtas and Mtasova [7], the choce of sutable nterpolaton methods for dfferent phenomena and dataset presents many dffcultes. The complexty of the modeled felds, dverse spatal data sampled non-optmally, presence of sgnfcant dscontnutes, and noses are common challenges. In addton, datasets obtaned from dverse sources wth dfferent accuraces are usually 3 6 very large ( N 0-0 ). Relable nterpolaton methods sutable for GIS applcatons should meet some necessary requrement such as accuracy and pre- 355

3 dctve power, dmensonal formulaton, drect estmaton of dervatves, applcablty to large datasets, 3 + computatonal effcency, and ease of use. Currently, t s dffcult to fnd a method that fulfls all the above-mentoned requrements for a wde range of georeferenced data. Therefore, the rght choce of the most adequate method wth approprate parameters for applcatons s paramount. Dfferent methods produce dfferent spatal representatons n dfferent datasets; also, n-depth knowledge of the phenomenon n queston s necessary n evaluatng whch of the nterpolaton methods produces results closest to realty. The use of an unsutable method or napproprate parameters can result n a dstorted model of spatal dstrbuton, leadng to potentally wrong decsons based on msleadng spatal nformaton. A wrong nterpolaton results becomes very crtcal when the estmates are nputs for smulatons, as small error or dstorton can cause models to produce false spatal patterns [7]. Whle external factors e.g. data densty, spatal dstrbuton of sample data, surface type, sample sze and samplng desgn, etc. [4] may affect the performance of an nterpolaton method, an understandng of the mathematcal formulaton of these methods wll provde some nformaton on ther sutablty for terran modelng. Ths paper attempts to examne the accuracy of spatal nterpolaton methods n modelng landform (topography) n relaton to ther mathematcal formulaton. The expermental study of ths work employs an area comprsng a slope and a plan as landform-adaptablty test area and focuses on the comparatve analyss of three commonly used nterpolaton methods of Krgng, Splne, and Inverse Dstance Weghtng, IDW. The followng secton summarzes the theoretcal and mathematcal bass of dfferent known nterpolaton methods ncludng the three nterpolaton methods n queston. Secton 3 ntroduces the accuracy analyss methods used n ths paper whle Secton 4 presents the expermental analyss. Secton 5 dscusses the results and Secton 6 concludes. 2. Mathematcal and Theoretcal Concept of IDW, Splne and Krgng Dfferent spatal nterpolaton methods have been developed n dfferent doman for dfferent applcatons. Accordng to [6], many standard lnear models are specal cases of a more general predcton model. Tobler s frst Law of Geography, everythng s related to everythng else, but near thngs are more related than dstant thngs [8] forms the general prncple of many nterpolaton methods. Some of the methods are only sutable for contnuous data felds that assume normal dstrbuton of the dataset. Spatal Interpolaton methods could be classfed nto two major groups: a) Mechancal/determnstc/non-geostatstcal methods; these nclude among other methods, Inverse Dstance Weghtng (IDW) and Splnes. b) Lnear statstcal/stochastc/geostatstcal methods; whch nclude Krgng among others [4], [6]. 356

4 2.. Inverse Dstance Weghtng (IDW) Ths method assumes that the value at an unknown locaton can be approxmated as a weghted average of values at ponts wthn a certan cut-off dstance, or from a gven number of the closest ponts (typcally 0 to 30). Weghts are usually nversely proportonal to a power of dstance [9] whch, at an un-sampled locaton s, leads to an estmator as contaned n Equaton (2) below. n ( ) wz ( s) F s = = = m = m j= ( ) z s p s s s s where p s a parameter (typcally = 2) [7]. IDW s a method that s easy to use and readly avalable; t frequently does not produce the local shape mpled by data and produces local extrema at the data ponts [7]. Some modfcatons have gven rse to a class of multvarate blended IDW surfaces and volumes [9]. The assumpton for IDW s that measured ponts closer to the unknown pont are more lke t than those that are further away n ther values. The weght s gven as: p d λ = (3) n = p d where d s the dstance between x 0 and x, p s a power parameter, and n s the number of measured ponts used for the estmaton. The man factor affectng the accuracy of IDW s the value of the power parameter [0]. Weghts dmnsh as the dstance ncreases, especally when the value of the power parameter ncreases, so nearby samples have a heaver weght and have more nfluence on the estmaton, and the resultant spatal nterpolaton s local [0]. The choce of power parameter and neghborhood sze s arbtrary []. The most popular choce of p s 2 and the resultng method s often called nverse square dstance or nverse dstance squared (IDS). The power parameter can also be chosen based on error measurement (e.g., mnmum mean absolute error), resultng n optmal IDW [2]. The smoothness of the estmated surface vares drectly wth the power parameter, and t s found that the estmated results become less satsfactory when p s and 2 compared wth when ps 4 [2]. IDW s referred to as movng average when ps zero [3], lnear nterpolaton when ps and weghted movng average when ps not equal to [3] Splne Splnes belong to a group of nterpolators called Radal Bass Functons (RBF). Methods n ths group nclude Thn-Plate Splne (TPS), Regularzed Splne wth Tenson, and Inverse Mult-Quadratc Splne [4]. These models use mathematcal functons to connect the sampled data ponts. They produce contnuous elevaton and grade surfaces whle lmtng the bendng of the surface produced to a p j (2) 357

5 mnmum. RBF models are best employed n smooth surfaces for whch the avalable sample data sze s large as ther performance s less than optmum for surfaces wth apprecable varatons spannng short ranges.rbf does not force estmates to mantan the range of the sampled data n these models lke n IDW [4]. Splne functons are the mathematcal equvalents of the flexble ruler cartographers used, called splnes, to ft smooth curves through several fxed ponts. It s a pecewse polynomal consstng of several sectons, each of whch s ftted to a small number of ponts n such a way that each of the sectons jon up at ponts referred to as break ponts. Ths has the advantage of accommodatng local adjustments, f there s a change n the data value at a pont, and s preferable to a smple polynomal nterpolaton because more parameters can be defned, ncludng the amount of smoothng [3]. Splnes are normally ftted usng low order polynomals (.e. second or thrd order) constraned to jon up. They may be two-dmensonal (e.g. when smoothng a contour lne) or three dmensonal (when modelng a surface). The smoothng splne functon also assumes the presence of a measurement error n the data that needs to be smoothed locally [3]. Among the many versons and modfcatons of splne nterpolators, the most wdely used technque s the thn-plate splnes [5] as well as the regularzed splne wth tenson and smoothng [7] Regularzed Splne wth Tenson For regularzed splne wth tenson and smoothng, the predcton s gven by: z( s ) = a + n w R( v ) (4) 0 = where a s a constant and R(v ) s the radal bass functon gven by: R ( v ) = E ( v ) + In ( v ) + C E (5) and v = ϕ h (6) where ( ) [ ] 2 E v s the exponental ntegral functon, C E = s the Euler constant, ϕ s the generalzed tenson parameter and h 0 s the separaton between the new and nterpolaton pont. The coeffcents a and w are obtaned by solvng the system, n = 0 2 w = 0 (7) ( ) δ ϖ ϖ ( ) n = j 0 a + w R v + = z s ; j =,, n. (8) whereϖ0 ϖ are postve weghtng factors for a smoothng parameter at each locaton [7]. The tenson parameter ϕ determnes the dstance over whch the gven ponts nfluence the resultng surface, whle the smoothng parameter controls the vertcal devaton of the surface from the sample locatons. The use of an approprate combnaton of tenson and smoothng produces a surface that correctly fts the emprcal knowledge about the expected varaton [7] Thn Plate Splne Wahba and Wendelberger [6] formulated thn plate splnes (TPS), prevously 358

6 called Laplacan smoothng splnes, for modelng clmatc data [4]. A basc soluton to the b-harmonc equaton, has the form z( r) = r 2 log r (9) where r s the dstance between sample ponts and un-sampled locatons [7]. The relaton below approxmates the surface wth mnmum bend ( ) = ( 0 ) n = f s a a x a y wz s s (0) where the terms a, axay 2, 3 model the lnear porton of the surface defnng a flat plan that best fts all control ponts usng least squares, the last term models the bendng forces due to m sampled ponts, w are control ponts coeffcents and s s 0 s the separaton of sampled pont s and locaton s 0. The unknowns a, axay 2, 3 and w are evaluated usng the relaton L V = w aaa () where L ( ) T K P 2 3 P 0 = T and V s a vector of pont heghts. K s a matrx of the dstance between sampled ponts and P s a matrx of the sampled ponts coordnates. L s obtaned by calculatng the nverse of L or solvng EQ. wth V replaced wth the matrx of the heghts of sampled ponts H padded wth zeros [8]. Once the unknowns are evaluated, one can compute EQ.0 to determne the heghts of unknown ponts.tps computes a smoothng factor by lmtng the Generalzed Cross Valdaton functon, GVC, makng for a comparatvely sturdy model as lmtng the GVC mproves the accuracy of estmatons and s less relant on the accuracy of the model tself. TPS gves a determnaton of spatal accuracy [7]. (2) Inverse Mult-Quadratc Splne The relaton below gves the nverse mult-quadratc splne functon f = + s s ( s) 2 (3) where s s s the Eucldean dstance between control ponts s and the unknown pont s [5]. The surface s modelled by the functon ( ) ( ) z s = n af s (4) = where the weghts such that a are selected to ensure exact estmatons at each data pont n ( ) ( ) z = z s = a f s, =,, n (5) j= j j and s computed by the relaton z = Fa (6) where z s replaced by a vector of sampled data values, F s a square functon matrx gven by, 359

7 FIJ = f ( sj ),, j =,, n (7) The estmaton functon generated wth these weghts s smooth and exact at sampled data ponts [9]. Splnes have been wdely seen as hghly sutable for estmaton of densely sampled heghts and clmatc varables [7], [5]. Among ts dsadvantages, the nablty to ntegrate larger amounts of auxlary maps n modelng the determnstc part of change as well as the arbtrary selecton of the smoothng and tenson parameters have been wdely crtczed [7]. Predctons obtaned from splnes therefore are largely dependent on decsons lke the order of polynomal used, number of break ponts, etc. taken by the user. Splnes may also be modeled not to be exact to avod the generaton of excessvely hgh or low values common wth some exact splnes [3]. Unlke the IDW methods, the values predcted by RBFs are not constraned to the range of measured values,.e., predcted values can be above the maxmum or below the mnmum measured value [4] Krgng Krgng, synonymous to geostatstcal nterpolaton, began n the mnng ndustry as a means of betterng ore reserve estmaton n the early 950 s [6]. Mnng engneer D. G. Krge and statstcan H. S. Schel formulated t. After almost a decade, French mathematcan G. Matheron derved the formulas, establshng the entre feld of lnear geostatstcs []. Krgng s founded on a concept of random functons wth the surface or volume assumed one realzaton of a random functon wth a known spatal covarance [7]. Regonalzed varable theory assumes that the spatal varaton of any varable can be expressed as the sum of the followng three components: a) A structural component havng a constant mean or trend. b) A regonalzed varable, whch s the random but spatally correlated component. c) A random but spatally uncorrelated nose or resdual component. d) Mathematcally, for a random varable z at x, the expresson s Z x = m x + ε x + ε (8) ( ) ( ) ( ) where m( x ) s a structural functon modelng the structural component, ε ( x ) s the spatally auto-correlated stochastc resdual from m( x ) (the regonalzed varable), and ε s random nose whch s normally dstrbuted wth an average of zero and a varance 2 σ [3]. Ordnary Krgng Ordnary Krgng (OK) s a standard verson of Krgng where predctons are based on the model, ( ) = µ + ε( ) Z s s (9) where µ s the fxed statonary functon or global average, and ε ( s ) s the stochastc but spatally correlated part of the varaton. Predctons are made as 360

8 n T ( 0) ( 0) ( ) 0 z s = w s z s = z (20) OK λ = where λ 0 s the vector of krgng weghts ( w ) and z s the vector of n samples at prmary locatons. Krgng, n a way, s an mprovement of the nverse dstance nterpolaton where the key problem of nverse dstance nterpolaton (the determnaton of how much mportance s gven to each neghbor) s addressed n such a way that the estmated weghts account for the true spatal autocorrelaton structure. The novelty Krgng has n the analyss of pont data s the dervaton and plottng of the sem-varance dfferences between the neghborng values [6]. Assumng statonarty, one can estmate a sem-varogram, γ ( h ), for data z( s ), defned as where h s the dstance between pont x and x 0 [2]. Ths relates to the spatal covarance C( h ) by, γ ( h) = var ( ) ( 0 ) 2 z x z x (2) ( ) ( 0) ( ) γ h = C C h (22) where C ( 0) s the sem-varogram value at nfnty(sll) [7]. The Ordnary Krgng (OK) weghts are evaluated by multplyng the covarances, = C C (23) λ 0 0 where C s the covarance matrx derved for an n x n samples matrx wth one addtonal row and column added to ensure the sum of weghts s equal to one, and C 0 s the vector of the covarances at a new locaton [6]. The covarance at a dstance of zero ( C ( 0) ) s by defnton the mean resdual error [6]. Expermental varograms usually have some characterstc features among whch are: ) Low values of h have small varance wth varance ncreasng n drect proporton to h, levelng off at a certan pont to form the sll. 2) At dstances less than the range (the dstance at whch the varance levels off), ponts closer together are more lkely to have smlar values than ponts further apart whle at dstances greater than the range, ponts have no nfluence upon themselves. The range therefore gves an dea of how large the search radus needs to be for a dstance-weghted nterpolaton. 3) The semvarance when h s zero has a postve value referred to as the nugget and ndcates the amount of non-spatally autocorrelated nose [3]. The semvarance dsplayed n an expermental varogram s modeled by a mathematcal functon dependng on the shape of the expermental varogram. A sphercal model s used when the varogram has a classc shape, an exponental model when the approach to the sll s more gradual. A Gaussan model s used when the nugget s small and the varaton s very smooth, and a lnear model when there s no sll. A varogram contanng a trend that has to be modeled separately s ncreasngly steep wth larger values of h. If the nugget varance s large and the varogram shows no tendency to gradually vansh wth smaller values of h, or the dstance between observatons s larger than the range (.e. 36

9 sample ponts are too far apart to nfluence one another), then nterpolaton s not reasonable and the best estmate s the overall mean of the observatons. Anose-flled varogram showng no partcular pattern may mean that the observatons are too few. A varogram that dps at dstances greater than the range to create a hole effect shows the sample space may be too small to reflect some long wave-length varaton n the data [3]. The nterpolated surface s constructed usng statstcal condtons of unbasedness and mnmum varance [7]. Three mportant requrements for Ordnary Krgng are: () The trend functon s fxed () The varogram s nvarant n the entre area of nterest () The target varable s (approxmately) normally dstrbuted. These requrements are often not met and consttute a serous dsadvantage of Ordnary Krgng [6]. However, a major beneft of the varous forms of krgng (and other stochastc nterpolaton methods) s that estmates of the model s predcton errors can be calculated, ncorporated n the analyss, and plotted along wth the predcted surface. Such error nformaton s an mportant tool n the spatal decson makng process [4]. 3. Accuracy Assessment Methodology The accuracy evaluaton ndces commonly used nclude, the Mean Error (ME), Mean Absolute Error (MAE), Mean Squared Error (MSE) and Root Mean Squared Error (RMSE). For n observatons, p predcted value, and o observed value these ndces are evaluated usng the expressons lsted below: ME = n ( p o ) n = MAE = n n p o = MSE = n ( p o ) 2 n = RMSE = n ( p o ) n = 2 2 (24) (25) (26) (27) [4] ME s used for determnng the degree of bas n the estmates often referred to as the bas [0]. Snce postve and negatve estmates counteract each other, the resultant ME tends to be lower than the actual error promptng cauton n ts use as an ndcator of accuracy [4]. RMSE provdes a measure of the error sze, but s senstve to outlers as t places a lot of weght on large errors [4]. MSE suffers the same drawbacks as RMSE. Whereas MAE s less senstve to extreme values [20] and ndcates the extent to whch the estmate can be n error [2]. MAE and RMSE are argued to be smlar measures and they gve estmates of the average error, but do not provde nformaton about the relatve sze of the average dfference and the nature of dfferences comprsng them [20]. Of course, cross-valdaton s used together wth these measurements to assess the performance of the nterpolaton methods [4]. In ths paper the ME and RMSE (aval- 362

10 able on the software used) are used to evaluate the performances of the IDW, Splne and Krgng nterpolaton methods consdered n the expermental analyss part of ths study. Cross valdaton s performed as further test of predcton accuracy. 4. Expermental Analyss The study area as descrbed n Fgure and Fgure 2 s an expanse of land measurng about twenty-one (2) hectares located n Ikot Ukapon lattude 50'0''N and longtude 7 59'0''E n the North-Eastern part of Akwa Ibom state of Ngera. The area s largely hlly wth a plan coverng about twenty (20%) percent of ts total area. Orthometrc heght data, obtaned from the Offce of the State Surveyor General, consst of four hundred and sxty-two (462) randomly sampled ponts spread across the area collected through feld survey usng a Kolda K9-T seres dfferental GPS. Data so obtaned from the feld samples was n Mcrosoft Excel.xls format. The lowest elevaton of the area s 27.08m whle the hghest elevaton s 98.79m above mean sea level. Prelmnary data exploraton usng ESRI s Arc- GIS (verson 9.3) shows the data s not normally dstrbuted as shown by dstrbuton parameters n Table. Fgure. Map of the studyng area showng the dstrbuton of data ponts wthn the study area. 363

11 Fgure 2. The contour map of the studyng area. Table. Dstrbuton parameters of the based elevaton data. Count 462 Mnmum value Maxmum value Mean Standard devaton Skewness Kurtoss Medan st quartle rd quartle IDW Method The followng contour fll surface shown n Fgure 3, s generated for IDW nterpolaton wth power, p of 2, smoothng factor of 0.5and neghborhood sze of 5 for the based data Splne Method Regularzed Splne nterpolaton, mplemented as Radal Bass Functon(RBF), wth order 2 gves the contour fll map n Fgure 4 for power = 2, smoothng factor = 0.5 and neghborhood sze = Krgng Krgng works on the assumpton that the data set s normalzed, therefore, we carred out Box-Cox normalzaton on the data before mplementng Krgng nterpolaton. 364

12 Fgure 3. IDW contour fll map. Fgure 4. The regularzed Splne contour fll map. The data was dvded nto a tranng and test subsets n a rato of 80:20 usng the Geostatstcal Analyst tool n ArcGIS and optmal parameter values (generated by ArcGIS and sets the best possble value for each parameter) used for predctons on the tranng subset. The test subset was then used wth these optmal parameters for valdaton. The data dstrbuton parameters after Box-Cox normalzaton s as shown n Table 2. The contour fll maps generated usng Krgng wth the Gausan model and 365

13 Table 2. Dstrbuton parameter of data after normalzaton usng Box-Cox wth transformaton. count 462 Mnmum value Maxmum value mean Standard devaton skewness kurtoss medan st quartle rd quartle Sphercal model respectvely for auto-calculated values for nugget, sll, mdrange, a lag sze of 54.40, and lag number = 2 are as shown Fgure 5(a) and Fgure 5(b). The maps for IDW, Splne and Krgng after the optmal valdaton of the data are shown n Fgures 6(a)-(c). 5. Result and Dscusson From the predcton errors tabulated n Table 3 and Table 4, the level of bas n estmaton s lowest for Krgng and hghest for IDW as ndcated by the respectve MEs. Ths presence of bas s expected, as the heght data was not random as shown n Table. Box-Cox normalzaton of the data dd not result n a normal dstrbuton ether (see Table 2), ths means the data was very based. Ths scenaro s frequently encountered n practce as feld collecton of elevaton data usually focuses on capturng perceved changes n elevaton rather than on randomness. The bas for the valdated data however, s hghest for Krgng and lowest for Splne lkely because of the reducton n the number of samples used for valdaton. A better measure of the error n predcton, the RMSE, s lowest for Splne and hghest for Krgng for both predcton and valdaton shown n Table 3 and Table 4. Ths ndcates an elevaton model that s closer to what s on ground for Splnes. However, outsde the areas where nterpolaton data where obtaned, Splne produces unrelable predctons. It s therefore not sutable for cases where data outsde the captured area s desred (extrapolaton). IDW produces a model that s better than that of Krgng but not as good as the Splne model. Its predctons outsde the captured area are also better than that of Krgng and Splne. Ths does not mply Krgng s not sutable for terran modelng or wll not perform better than both IDW and Splne. Krgng assumes normal dstrbuton of data and models the spatal dstrbuton of a geographcal event as a realzaton of a functon that s random. Its predctons therefore are dependent on the data satsfyng the statstcal crtera of unbasedness and mnmumvarance. Its mathematcal formulaton makes t unsutable for data that 366

14 (a) (b) Fgure 5. (a) and (b) depcts the Krgng contour fll map after Box-Cox normalzaton created wth Gaussan model and wth sphercal model respectvely. 367

15 (a) (b) (c) Fgure 6. (a) Inverse Dstance Weghted (IDW) Predcton Map; (b) Splne; (c) Ordnary Krgng Predcton Map. 368

16 Table 3. Predcton errors for the three nterpolaton methods at optmal parameters. METHOD IDW SPLINE KRIGING ME RMSE ME s the mean predcton error and RMSE s the root mean square standardzed error of predcton. Table 4. Predcton errors after valdaton. METHOD IDW SPLINE KRIGING ME RMSE ME s the mean predcton error and RMSE s the root mean square standardzed error of predcton. s not normally dstrbuted or dffcult to normalze. Splnes on the other hand use a physcal model varyng n accordance to the varaton n the elastc propertes of the estmaton functon. It tends to do well wth modelng physcal phenomena such as terran. IDW uses a lnear combnaton of values at captured event locatons, assgns weghts by an nverse functon of the separaton between the event locaton to be estmated and ponts captured to estmate values of the unknown locaton. Though weghts are specfed arbtrarly, ArcGIS software provdes an optmal weght management functon that assgns a weght that s most sutable for ponts wthn the captured data set. Predctons are nfluenced by ths weght assgnment but are more relable n terms of error than what s obtaned usng Krgng. It s acknowledged that the Krgng does very well wth covarate data such as temperature data, but the data has to be captured as randomly as possble. Ths s often not acheved. A good knowledge of the data used as well as the strengths and weaknesses of the avalable nterpolaton methods s necessary n decdng on a method to use for nterpolaton for a gven purpose. 6. Concluson In ths study, Splne provdes a more accurate model and result for the elevaton data obtaned drectly from feld survey that was not homogenously randomzed and not normalzed. From the nterpolaton result we obtaned, Splne method outsde the data area also reaffrms that predctons by RBFs are not constraned to the range of measured values,.e., predcted values can be above the maxmum or below the mnmum measured value. Tan and Xu [22] concluded from ther experment on terran modelng usng data from a dgtzed map that IDW gave a better model n terms of accuracy than Splne or Krgng. Ths s most lkely due to means from whch the test data were acqured. Ther test data were dgtzed from a contour map and were homogenously dstrbuted. The knowledge of the source of data may therefore be of mportance n the choce of nterpolaton method. References [] Kreyszg, E. (2003) Advanced Engneerng Mathematcs. 8th Edton, John Wley 369

17 and Sons (Asa) Pte., Ltd., Sngapore. [2] Collns, F.C. and Bolstad, P.V. (996) A Comparson of Spatal Interpolaton Technques n Temperature Estmaton. Proceedngs of the 3rd Internatonal Conference/Workshop on Integratng GIS and Envronmental Modelng, Natonal Center for Geographc Informaton and Analyss, Santa Barbara, Santa Fe, NM; Santa Barbara, CA. [3] Burrough, P.A. and McDonnell, R.A. (998) Prncples of Geographcal Informaton Systems. Oxford Unversty Press, Oxford. [4] L, J. and Heap, A.D. (2008) A Revew of Spatal Interpolaton Methods for Envronmental Scentsts. Geoscence Australa, Record 2008/23, 37. [5] Heuvelnk, G.B.M. and Webster, R. (200) Modelng Sol Varaton: Past, Present, and Future. Geoderma, 00, [6] Hengl, T. (2009) A Practcal Gude to Geostatstcal Mappng. [7] Mtas, L. and Mtasova, H. (2005) Spatal Interpolaton. In: Longley, P.A., Goodchld, M.F., Magure, D.J. and Rhnd, D.W., Eds., Geographc Informaton Systems: Prncples, Technques, Management and Applcatons, 2nd Edton, Vol., Part 2, Chapter [8] Tobler, W.R. (970) A Computer Model Smulaton of Urban Growth n the Detrot Regon. Economc Geography, 46, [9] Watson, D.F. (992) Contourng: A Gude to the Analyss and Dsplay of Spatal Data. Pergamon, Oxford. [0] Isaaks, E.H. and Srvastava, R.M. (989) Appled Geostatstcs. Oxford Unversty Press, New York. [] Webster, R. and Olver, M. (200) Geostatstcs for Envronmental Scentsts. John Wley & Sons, Ltd., Chchester, 27. [2] Rpley, B.D. (98) Spatal Statstcs. John Wley & Sons, New York, [3] Brus, D.J., et al. (996) The Performance of Spatal Interpolaton Methods and Choropleth Maps to Estmate Propertes at Ponts: A Sol Survey Case Study. Envronmetrcs, 7, [4] Unted States Envronmental Protecton Agency (EPA) (2004) Developng Spatally Interpolated Surfaces and Estmatng Uncertanty. EPA-454/R , EPA Publcaton, NC. [5] Hutchnson, M.F. (995) Interpolatng Mean Ranfall Usng Thn Plate Smoothng Splnes. Internatonal Journal of Geographcal Informaton Systems, 9, [6] Wahba, G. and Wendelberger, J. (980) Some New Mathematcal Methods for Varatonal Objectve Analyss Usng Splnes and Cross-Valdaton. Monthly Weather Revew, 08, [7] Belonge, S. (206) Thn Plate Splne. Retreved onlne from MathWorld A Wolfram Web Resource, created by Erc W. Wessten. [8] Lombaert, H. (2006) Manual Regstraton wth Thn Plates

18 [9] Smth, T.E. (206) Notebook on Spatal Data Analyss. [20] Wllmott, C.J. (982) Some Comments on the Evaluaton of Model Performance. Bulletn Amercan Meteorologcal Socety, 63, [2] Nalder, I.A. and Wen, R.W. (998) Spatal Interpolaton of Clmatc Normals: Test of a New Method n the Canadan Boreal Forest. Agrcultural and Forest Meteorology, 92, [22] Tan, Q. and Xu, X. (204) Comparatve Analyss of Spatal Interpolaton Methods. Sensors and Transducers, 65, Submt or recommend next manuscrpt to SCIRP and we wll provde best servce for you: Acceptng pre-submsson nqures through Emal, Facebook, LnkedIn, Twtter, etc. A wde selecton of journals (nclusve of 9 subjects, more than 200 journals) Provdng 24-hour hgh-qualty servce User-frendly onlne submsson system Far and swft peer-revew system Effcent typesettng and proofreadng procedure Dsplay of the result of downloads and vsts, as well as the number of cted artcles Maxmum dssemnaton of your research work Submt your manuscrpt at: Or contact jgs@scrp.org 37

S1 Note. Basis functions.

S1 Note. Basis functions. S1 Note. Bass functons. Contents Types of bass functons...1 The Fourer bass...2 B-splne bass...3 Power and type I error rates wth dfferent numbers of bass functons...4 Table S1. Smulaton results of type

More information

Smoothing Spline ANOVA for variable screening

Smoothing Spline ANOVA for variable screening Smoothng Splne ANOVA for varable screenng a useful tool for metamodels tranng and mult-objectve optmzaton L. Rcco, E. Rgon, A. Turco Outlne RSM Introducton Possble couplng Test case MOO MOO wth Game Theory

More information

Feature Reduction and Selection

Feature Reduction and Selection Feature Reducton and Selecton Dr. Shuang LIANG School of Software Engneerng TongJ Unversty Fall, 2012 Today s Topcs Introducton Problems of Dmensonalty Feature Reducton Statstc methods Prncpal Components

More information

APPLICATION OF A COMPUTATIONALLY EFFICIENT GEOSTATISTICAL APPROACH TO CHARACTERIZING VARIABLY SPACED WATER-TABLE DATA

APPLICATION OF A COMPUTATIONALLY EFFICIENT GEOSTATISTICAL APPROACH TO CHARACTERIZING VARIABLY SPACED WATER-TABLE DATA RFr"W/FZD JAN 2 4 1995 OST control # 1385 John J Q U ~ M Argonne Natonal Laboratory Argonne, L 60439 Tel: 708-252-5357, Fax: 708-252-3 611 APPLCATON OF A COMPUTATONALLY EFFCENT GEOSTATSTCAL APPROACH TO

More information

Determining the Optimal Bandwidth Based on Multi-criterion Fusion

Determining the Optimal Bandwidth Based on Multi-criterion Fusion Proceedngs of 01 4th Internatonal Conference on Machne Learnng and Computng IPCSIT vol. 5 (01) (01) IACSIT Press, Sngapore Determnng the Optmal Bandwdth Based on Mult-crteron Fuson Ha-L Lang 1+, Xan-Mn

More information

Cluster Analysis of Electrical Behavior

Cluster Analysis of Electrical Behavior Journal of Computer and Communcatons, 205, 3, 88-93 Publshed Onlne May 205 n ScRes. http://www.scrp.org/ournal/cc http://dx.do.org/0.4236/cc.205.350 Cluster Analyss of Electrcal Behavor Ln Lu Ln Lu, School

More information

Page 0 of 0 SPATIAL INTERPOLATION METHODS

Page 0 of 0 SPATIAL INTERPOLATION METHODS Page 0 of 0 SPATIAL INTERPOLATION METHODS 2018 1. Introducton Spatal nterpolaton s the procedure to predct the value of attrbutes at unobserved ponts wthn a study regon usng exstng observatons (Waters,

More information

Subspace clustering. Clustering. Fundamental to all clustering techniques is the choice of distance measure between data points;

Subspace clustering. Clustering. Fundamental to all clustering techniques is the choice of distance measure between data points; Subspace clusterng Clusterng Fundamental to all clusterng technques s the choce of dstance measure between data ponts; D q ( ) ( ) 2 x x = x x, j k = 1 k jk Squared Eucldean dstance Assumpton: All features

More information

SLAM Summer School 2006 Practical 2: SLAM using Monocular Vision

SLAM Summer School 2006 Practical 2: SLAM using Monocular Vision SLAM Summer School 2006 Practcal 2: SLAM usng Monocular Vson Javer Cvera, Unversty of Zaragoza Andrew J. Davson, Imperal College London J.M.M Montel, Unversty of Zaragoza. josemar@unzar.es, jcvera@unzar.es,

More information

Proper Choice of Data Used for the Estimation of Datum Transformation Parameters

Proper Choice of Data Used for the Estimation of Datum Transformation Parameters Proper Choce of Data Used for the Estmaton of Datum Transformaton Parameters Hakan S. KUTOGLU, Turkey Key words: Coordnate systems; transformaton; estmaton, relablty. SUMMARY Advances n technologes and

More information

NUMERICAL SOLVING OPTIMAL CONTROL PROBLEMS BY THE METHOD OF VARIATIONS

NUMERICAL SOLVING OPTIMAL CONTROL PROBLEMS BY THE METHOD OF VARIATIONS ARPN Journal of Engneerng and Appled Scences 006-017 Asan Research Publshng Network (ARPN). All rghts reserved. NUMERICAL SOLVING OPTIMAL CONTROL PROBLEMS BY THE METHOD OF VARIATIONS Igor Grgoryev, Svetlana

More information

Simulation: Solving Dynamic Models ABE 5646 Week 11 Chapter 2, Spring 2010

Simulation: Solving Dynamic Models ABE 5646 Week 11 Chapter 2, Spring 2010 Smulaton: Solvng Dynamc Models ABE 5646 Week Chapter 2, Sprng 200 Week Descrpton Readng Materal Mar 5- Mar 9 Evaluatng [Crop] Models Comparng a model wth data - Graphcal, errors - Measures of agreement

More information

X- Chart Using ANOM Approach

X- Chart Using ANOM Approach ISSN 1684-8403 Journal of Statstcs Volume 17, 010, pp. 3-3 Abstract X- Chart Usng ANOM Approach Gullapall Chakravarth 1 and Chaluvad Venkateswara Rao Control lmts for ndvdual measurements (X) chart are

More information

y and the total sum of

y and the total sum of Lnear regresson Testng for non-lnearty In analytcal chemstry, lnear regresson s commonly used n the constructon of calbraton functons requred for analytcal technques such as gas chromatography, atomc absorpton

More information

Classifier Selection Based on Data Complexity Measures *

Classifier Selection Based on Data Complexity Measures * Classfer Selecton Based on Data Complexty Measures * Edth Hernández-Reyes, J.A. Carrasco-Ochoa, and J.Fco. Martínez-Trndad Natonal Insttute for Astrophyscs, Optcs and Electroncs, Lus Enrque Erro No.1 Sta.

More information

Mathematics 256 a course in differential equations for engineering students

Mathematics 256 a course in differential equations for engineering students Mathematcs 56 a course n dfferental equatons for engneerng students Chapter 5. More effcent methods of numercal soluton Euler s method s qute neffcent. Because the error s essentally proportonal to the

More information

Wishing you all a Total Quality New Year!

Wishing you all a Total Quality New Year! Total Qualty Management and Sx Sgma Post Graduate Program 214-15 Sesson 4 Vnay Kumar Kalakband Assstant Professor Operatons & Systems Area 1 Wshng you all a Total Qualty New Year! Hope you acheve Sx sgma

More information

NAG Fortran Library Chapter Introduction. G10 Smoothing in Statistics

NAG Fortran Library Chapter Introduction. G10 Smoothing in Statistics Introducton G10 NAG Fortran Lbrary Chapter Introducton G10 Smoothng n Statstcs Contents 1 Scope of the Chapter... 2 2 Background to the Problems... 2 2.1 Smoothng Methods... 2 2.2 Smoothng Splnes and Regresson

More information

Lecture #15 Lecture Notes

Lecture #15 Lecture Notes Lecture #15 Lecture Notes The ocean water column s very much a 3-D spatal entt and we need to represent that structure n an economcal way to deal wth t n calculatons. We wll dscuss one way to do so, emprcal

More information

An Accurate Evaluation of Integrals in Convex and Non convex Polygonal Domain by Twelve Node Quadrilateral Finite Element Method

An Accurate Evaluation of Integrals in Convex and Non convex Polygonal Domain by Twelve Node Quadrilateral Finite Element Method Internatonal Journal of Computatonal and Appled Mathematcs. ISSN 89-4966 Volume, Number (07), pp. 33-4 Research Inda Publcatons http://www.rpublcaton.com An Accurate Evaluaton of Integrals n Convex and

More information

THE THEORY OF REGIONALIZED VARIABLES

THE THEORY OF REGIONALIZED VARIABLES CHAPTER 4 THE THEORY OF REGIONALIZED VARIABLES 4.1 Introducton It s ponted out by Armstrong (1998 : 16) that Matheron (1963b), realzng the sgnfcance of the spatal aspect of geostatstcal data, coned the

More information

Some Advanced SPC Tools 1. Cumulative Sum Control (Cusum) Chart For the data shown in Table 9-1, the x chart can be generated.

Some Advanced SPC Tools 1. Cumulative Sum Control (Cusum) Chart For the data shown in Table 9-1, the x chart can be generated. Some Advanced SP Tools 1. umulatve Sum ontrol (usum) hart For the data shown n Table 9-1, the x chart can be generated. However, the shft taken place at sample #21 s not apparent. 92 For ths set samples,

More information

R s s f. m y s. SPH3UW Unit 7.3 Spherical Concave Mirrors Page 1 of 12. Notes

R s s f. m y s. SPH3UW Unit 7.3 Spherical Concave Mirrors Page 1 of 12. Notes SPH3UW Unt 7.3 Sphercal Concave Mrrors Page 1 of 1 Notes Physcs Tool box Concave Mrror If the reflectng surface takes place on the nner surface of the sphercal shape so that the centre of the mrror bulges

More information

FEATURE EXTRACTION. Dr. K.Vijayarekha. Associate Dean School of Electrical and Electronics Engineering SASTRA University, Thanjavur

FEATURE EXTRACTION. Dr. K.Vijayarekha. Associate Dean School of Electrical and Electronics Engineering SASTRA University, Thanjavur FEATURE EXTRACTION Dr. K.Vjayarekha Assocate Dean School of Electrcal and Electroncs Engneerng SASTRA Unversty, Thanjavur613 41 Jont Intatve of IITs and IISc Funded by MHRD Page 1 of 8 Table of Contents

More information

BioTechnology. An Indian Journal FULL PAPER. Trade Science Inc.

BioTechnology. An Indian Journal FULL PAPER. Trade Science Inc. [Type text] [Type text] [Type text] ISSN : 0974-74 Volume 0 Issue BoTechnology 04 An Indan Journal FULL PAPER BTAIJ 0() 04 [684-689] Revew on Chna s sports ndustry fnancng market based on market -orented

More information

Hermite Splines in Lie Groups as Products of Geodesics

Hermite Splines in Lie Groups as Products of Geodesics Hermte Splnes n Le Groups as Products of Geodescs Ethan Eade Updated May 28, 2017 1 Introducton 1.1 Goal Ths document defnes a curve n the Le group G parametrzed by tme and by structural parameters n the

More information

Learning the Kernel Parameters in Kernel Minimum Distance Classifier

Learning the Kernel Parameters in Kernel Minimum Distance Classifier Learnng the Kernel Parameters n Kernel Mnmum Dstance Classfer Daoqang Zhang 1,, Songcan Chen and Zh-Hua Zhou 1* 1 Natonal Laboratory for Novel Software Technology Nanjng Unversty, Nanjng 193, Chna Department

More information

A Semi-parametric Regression Model to Estimate Variability of NO 2

A Semi-parametric Regression Model to Estimate Variability of NO 2 Envronment and Polluton; Vol. 2, No. 1; 2013 ISSN 1927-0909 E-ISSN 1927-0917 Publshed by Canadan Center of Scence and Educaton A Sem-parametrc Regresson Model to Estmate Varablty of NO 2 Meczysław Szyszkowcz

More information

CS 534: Computer Vision Model Fitting

CS 534: Computer Vision Model Fitting CS 534: Computer Vson Model Fttng Sprng 004 Ahmed Elgammal Dept of Computer Scence CS 534 Model Fttng - 1 Outlnes Model fttng s mportant Least-squares fttng Maxmum lkelhood estmaton MAP estmaton Robust

More information

Empirical Distributions of Parameter Estimates. in Binary Logistic Regression Using Bootstrap

Empirical Distributions of Parameter Estimates. in Binary Logistic Regression Using Bootstrap Int. Journal of Math. Analyss, Vol. 8, 4, no. 5, 7-7 HIKARI Ltd, www.m-hkar.com http://dx.do.org/.988/jma.4.494 Emprcal Dstrbutons of Parameter Estmates n Bnary Logstc Regresson Usng Bootstrap Anwar Ftranto*

More information

Content Based Image Retrieval Using 2-D Discrete Wavelet with Texture Feature with Different Classifiers

Content Based Image Retrieval Using 2-D Discrete Wavelet with Texture Feature with Different Classifiers IOSR Journal of Electroncs and Communcaton Engneerng (IOSR-JECE) e-issn: 78-834,p- ISSN: 78-8735.Volume 9, Issue, Ver. IV (Mar - Apr. 04), PP 0-07 Content Based Image Retreval Usng -D Dscrete Wavelet wth

More information

A New Approach For the Ranking of Fuzzy Sets With Different Heights

A New Approach For the Ranking of Fuzzy Sets With Different Heights New pproach For the ankng of Fuzzy Sets Wth Dfferent Heghts Pushpnder Sngh School of Mathematcs Computer pplcatons Thapar Unversty, Patala-7 00 Inda pushpndersnl@gmalcom STCT ankng of fuzzy sets plays

More information

Problem Definitions and Evaluation Criteria for Computational Expensive Optimization

Problem Definitions and Evaluation Criteria for Computational Expensive Optimization Problem efntons and Evaluaton Crtera for Computatonal Expensve Optmzaton B. Lu 1, Q. Chen and Q. Zhang 3, J. J. Lang 4, P. N. Suganthan, B. Y. Qu 6 1 epartment of Computng, Glyndwr Unversty, UK Faclty

More information

Modeling Local Uncertainty accounting for Uncertainty in the Data

Modeling Local Uncertainty accounting for Uncertainty in the Data Modelng Local Uncertanty accontng for Uncertanty n the Data Olena Babak and Clayton V Detsch Consder the problem of estmaton at an nsampled locaton sng srrondng samples The standard approach to ths problem

More information

Available online at ScienceDirect. Procedia Environmental Sciences 26 (2015 )

Available online at   ScienceDirect. Procedia Environmental Sciences 26 (2015 ) Avalable onlne at www.scencedrect.com ScenceDrect Proceda Envronmental Scences 26 (2015 ) 109 114 Spatal Statstcs 2015: Emergng Patterns Calbratng a Geographcally Weghted Regresson Model wth Parameter-Specfc

More information

Image Representation & Visualization Basic Imaging Algorithms Shape Representation and Analysis. outline

Image Representation & Visualization Basic Imaging Algorithms Shape Representation and Analysis. outline mage Vsualzaton mage Vsualzaton mage Representaton & Vsualzaton Basc magng Algorthms Shape Representaton and Analyss outlne mage Representaton & Vsualzaton Basc magng Algorthms Shape Representaton and

More information

A mathematical programming approach to the analysis, design and scheduling of offshore oilfields

A mathematical programming approach to the analysis, design and scheduling of offshore oilfields 17 th European Symposum on Computer Aded Process Engneerng ESCAPE17 V. Plesu and P.S. Agach (Edtors) 2007 Elsever B.V. All rghts reserved. 1 A mathematcal programmng approach to the analyss, desgn and

More information

The Research of Support Vector Machine in Agricultural Data Classification

The Research of Support Vector Machine in Agricultural Data Classification The Research of Support Vector Machne n Agrcultural Data Classfcaton Le Sh, Qguo Duan, Xnmng Ma, Me Weng College of Informaton and Management Scence, HeNan Agrcultural Unversty, Zhengzhou 45000 Chna Zhengzhou

More information

A MOVING MESH APPROACH FOR SIMULATION BUDGET ALLOCATION ON CONTINUOUS DOMAINS

A MOVING MESH APPROACH FOR SIMULATION BUDGET ALLOCATION ON CONTINUOUS DOMAINS Proceedngs of the Wnter Smulaton Conference M E Kuhl, N M Steger, F B Armstrong, and J A Jones, eds A MOVING MESH APPROACH FOR SIMULATION BUDGET ALLOCATION ON CONTINUOUS DOMAINS Mark W Brantley Chun-Hung

More information

Synthesizer 1.0. User s Guide. A Varying Coefficient Meta. nalytic Tool. Z. Krizan Employing Microsoft Excel 2007

Synthesizer 1.0. User s Guide. A Varying Coefficient Meta. nalytic Tool. Z. Krizan Employing Microsoft Excel 2007 Syntheszer 1.0 A Varyng Coeffcent Meta Meta-Analytc nalytc Tool Employng Mcrosoft Excel 007.38.17.5 User s Gude Z. Krzan 009 Table of Contents 1. Introducton and Acknowledgments 3. Operatonal Functons

More information

Quality Improvement Algorithm for Tetrahedral Mesh Based on Optimal Delaunay Triangulation

Quality Improvement Algorithm for Tetrahedral Mesh Based on Optimal Delaunay Triangulation Intellgent Informaton Management, 013, 5, 191-195 Publshed Onlne November 013 (http://www.scrp.org/journal/m) http://dx.do.org/10.36/m.013.5601 Qualty Improvement Algorthm for Tetrahedral Mesh Based on

More information

Support Vector Machines

Support Vector Machines /9/207 MIST.6060 Busness Intellgence and Data Mnng What are Support Vector Machnes? Support Vector Machnes Support Vector Machnes (SVMs) are supervsed learnng technques that analyze data and recognze patterns.

More information

A Binarization Algorithm specialized on Document Images and Photos

A Binarization Algorithm specialized on Document Images and Photos A Bnarzaton Algorthm specalzed on Document mages and Photos Ergna Kavalleratou Dept. of nformaton and Communcaton Systems Engneerng Unversty of the Aegean kavalleratou@aegean.gr Abstract n ths paper, a

More information

Cell Count Method on a Network with SANET

Cell Count Method on a Network with SANET CSIS Dscusson Paper No.59 Cell Count Method on a Network wth SANET Atsuyuk Okabe* and Shno Shode** Center for Spatal Informaton Scence, Unversty of Tokyo 7-3-1, Hongo, Bunkyo-ku, Tokyo 113-8656, Japan

More information

TN348: Openlab Module - Colocalization

TN348: Openlab Module - Colocalization TN348: Openlab Module - Colocalzaton Topc The Colocalzaton module provdes the faclty to vsualze and quantfy colocalzaton between pars of mages. The Colocalzaton wndow contans a prevew of the two mages

More information

Improvement of Spatial Resolution Using BlockMatching Based Motion Estimation and Frame. Integration

Improvement of Spatial Resolution Using BlockMatching Based Motion Estimation and Frame. Integration Improvement of Spatal Resoluton Usng BlockMatchng Based Moton Estmaton and Frame Integraton Danya Suga and Takayuk Hamamoto Graduate School of Engneerng, Tokyo Unversty of Scence, 6-3-1, Nuku, Katsuska-ku,

More information

Helsinki University Of Technology, Systems Analysis Laboratory Mat Independent research projects in applied mathematics (3 cr)

Helsinki University Of Technology, Systems Analysis Laboratory Mat Independent research projects in applied mathematics (3 cr) Helsnk Unversty Of Technology, Systems Analyss Laboratory Mat-2.08 Independent research projects n appled mathematcs (3 cr) "! #$&% Antt Laukkanen 506 R ajlaukka@cc.hut.f 2 Introducton...3 2 Multattrbute

More information

Performance Evaluation of Information Retrieval Systems

Performance Evaluation of Information Retrieval Systems Why System Evaluaton? Performance Evaluaton of Informaton Retreval Systems Many sldes n ths secton are adapted from Prof. Joydeep Ghosh (UT ECE) who n turn adapted them from Prof. Dk Lee (Unv. of Scence

More information

Machine Learning 9. week

Machine Learning 9. week Machne Learnng 9. week Mappng Concept Radal Bass Functons (RBF) RBF Networks 1 Mappng It s probably the best scenaro for the classfcaton of two dataset s to separate them lnearly. As you see n the below

More information

An Iterative Solution Approach to Process Plant Layout using Mixed Integer Optimisation

An Iterative Solution Approach to Process Plant Layout using Mixed Integer Optimisation 17 th European Symposum on Computer Aded Process Engneerng ESCAPE17 V. Plesu and P.S. Agach (Edtors) 2007 Elsever B.V. All rghts reserved. 1 An Iteratve Soluton Approach to Process Plant Layout usng Mxed

More information

REFRACTIVE INDEX SELECTION FOR POWDER MIXTURES

REFRACTIVE INDEX SELECTION FOR POWDER MIXTURES REFRACTIVE INDEX SELECTION FOR POWDER MIXTURES Laser dffracton s one of the most wdely used methods for partcle sze analyss of mcron and submcron sze powders and dspersons. It s quck and easy and provdes

More information

For instance, ; the five basic number-sets are increasingly more n A B & B A A = B (1)

For instance, ; the five basic number-sets are increasingly more n A B & B A A = B (1) Secton 1.2 Subsets and the Boolean operatons on sets If every element of the set A s an element of the set B, we say that A s a subset of B, or that A s contaned n B, or that B contans A, and we wrte A

More information

Tsinghua University at TAC 2009: Summarizing Multi-documents by Information Distance

Tsinghua University at TAC 2009: Summarizing Multi-documents by Information Distance Tsnghua Unversty at TAC 2009: Summarzng Mult-documents by Informaton Dstance Chong Long, Mnle Huang, Xaoyan Zhu State Key Laboratory of Intellgent Technology and Systems, Tsnghua Natonal Laboratory for

More information

Review of approximation techniques

Review of approximation techniques CHAPTER 2 Revew of appromaton technques 2. Introducton Optmzaton problems n engneerng desgn are characterzed by the followng assocated features: the objectve functon and constrants are mplct functons evaluated

More information

Load Balancing for Hex-Cell Interconnection Network

Load Balancing for Hex-Cell Interconnection Network Int. J. Communcatons, Network and System Scences,,, - Publshed Onlne Aprl n ScRes. http://www.scrp.org/journal/jcns http://dx.do.org/./jcns.. Load Balancng for Hex-Cell Interconnecton Network Saher Manaseer,

More information

Lobachevsky State University of Nizhni Novgorod. Polyhedron. Quick Start Guide

Lobachevsky State University of Nizhni Novgorod. Polyhedron. Quick Start Guide Lobachevsky State Unversty of Nzhn Novgorod Polyhedron Quck Start Gude Nzhn Novgorod 2016 Contents Specfcaton of Polyhedron software... 3 Theoretcal background... 4 1. Interface of Polyhedron... 6 1.1.

More information

Programming in Fortran 90 : 2017/2018

Programming in Fortran 90 : 2017/2018 Programmng n Fortran 90 : 2017/2018 Programmng n Fortran 90 : 2017/2018 Exercse 1 : Evaluaton of functon dependng on nput Wrte a program who evaluate the functon f (x,y) for any two user specfed values

More information

Wavefront Reconstructor

Wavefront Reconstructor A Dstrbuted Smplex B-Splne Based Wavefront Reconstructor Coen de Vsser and Mchel Verhaegen 14-12-201212 2012 Delft Unversty of Technology Contents Introducton Wavefront reconstructon usng Smplex B-Splnes

More information

2x x l. Module 3: Element Properties Lecture 4: Lagrange and Serendipity Elements

2x x l. Module 3: Element Properties Lecture 4: Lagrange and Serendipity Elements Module 3: Element Propertes Lecture : Lagrange and Serendpty Elements 5 In last lecture note, the nterpolaton functons are derved on the bass of assumed polynomal from Pascal s trangle for the fled varable.

More information

APPLICATION OF MULTIVARIATE LOSS FUNCTION FOR ASSESSMENT OF THE QUALITY OF TECHNOLOGICAL PROCESS MANAGEMENT

APPLICATION OF MULTIVARIATE LOSS FUNCTION FOR ASSESSMENT OF THE QUALITY OF TECHNOLOGICAL PROCESS MANAGEMENT 3. - 5. 5., Brno, Czech Republc, EU APPLICATION OF MULTIVARIATE LOSS FUNCTION FOR ASSESSMENT OF THE QUALITY OF TECHNOLOGICAL PROCESS MANAGEMENT Abstract Josef TOŠENOVSKÝ ) Lenka MONSPORTOVÁ ) Flp TOŠENOVSKÝ

More information

Outline. Type of Machine Learning. Examples of Application. Unsupervised Learning

Outline. Type of Machine Learning. Examples of Application. Unsupervised Learning Outlne Artfcal Intellgence and ts applcatons Lecture 8 Unsupervsed Learnng Professor Danel Yeung danyeung@eee.org Dr. Patrck Chan patrckchan@eee.org South Chna Unversty of Technology, Chna Introducton

More information

Parallelism for Nested Loops with Non-uniform and Flow Dependences

Parallelism for Nested Loops with Non-uniform and Flow Dependences Parallelsm for Nested Loops wth Non-unform and Flow Dependences Sam-Jn Jeong Dept. of Informaton & Communcaton Engneerng, Cheonan Unversty, 5, Anseo-dong, Cheonan, Chungnam, 330-80, Korea. seong@cheonan.ac.kr

More information

Classifying Acoustic Transient Signals Using Artificial Intelligence

Classifying Acoustic Transient Signals Using Artificial Intelligence Classfyng Acoustc Transent Sgnals Usng Artfcal Intellgence Steve Sutton, Unversty of North Carolna At Wlmngton (suttons@charter.net) Greg Huff, Unversty of North Carolna At Wlmngton (jgh7476@uncwl.edu)

More information

A CLASS OF TRANSFORMED EFFICIENT RATIO ESTIMATORS OF FINITE POPULATION MEAN. Department of Statistics, Islamia College, Peshawar, Pakistan 2

A CLASS OF TRANSFORMED EFFICIENT RATIO ESTIMATORS OF FINITE POPULATION MEAN. Department of Statistics, Islamia College, Peshawar, Pakistan 2 Pa. J. Statst. 5 Vol. 3(4), 353-36 A CLASS OF TRANSFORMED EFFICIENT RATIO ESTIMATORS OF FINITE POPULATION MEAN Sajjad Ahmad Khan, Hameed Al, Sadaf Manzoor and Alamgr Department of Statstcs, Islama College,

More information

Course Introduction. Algorithm 8/31/2017. COSC 320 Advanced Data Structures and Algorithms. COSC 320 Advanced Data Structures and Algorithms

Course Introduction. Algorithm 8/31/2017. COSC 320 Advanced Data Structures and Algorithms. COSC 320 Advanced Data Structures and Algorithms Course Introducton Course Topcs Exams, abs, Proects A quc loo at a few algorthms 1 Advanced Data Structures and Algorthms Descrpton: We are gong to dscuss algorthm complexty analyss, algorthm desgn technques

More information

An Optimal Algorithm for Prufer Codes *

An Optimal Algorithm for Prufer Codes * J. Software Engneerng & Applcatons, 2009, 2: 111-115 do:10.4236/jsea.2009.22016 Publshed Onlne July 2009 (www.scrp.org/journal/jsea) An Optmal Algorthm for Prufer Codes * Xaodong Wang 1, 2, Le Wang 3,

More information

Edge Detection in Noisy Images Using the Support Vector Machines

Edge Detection in Noisy Images Using the Support Vector Machines Edge Detecton n Nosy Images Usng the Support Vector Machnes Hlaro Gómez-Moreno, Saturnno Maldonado-Bascón, Francsco López-Ferreras Sgnal Theory and Communcatons Department. Unversty of Alcalá Crta. Madrd-Barcelona

More information

Steps for Computing the Dissimilarity, Entropy, Herfindahl-Hirschman and. Accessibility (Gravity with Competition) Indices

Steps for Computing the Dissimilarity, Entropy, Herfindahl-Hirschman and. Accessibility (Gravity with Competition) Indices Steps for Computng the Dssmlarty, Entropy, Herfndahl-Hrschman and Accessblty (Gravty wth Competton) Indces I. Dssmlarty Index Measurement: The followng formula can be used to measure the evenness between

More information

Type-2 Fuzzy Non-uniform Rational B-spline Model with Type-2 Fuzzy Data

Type-2 Fuzzy Non-uniform Rational B-spline Model with Type-2 Fuzzy Data Malaysan Journal of Mathematcal Scences 11(S) Aprl : 35 46 (2017) Specal Issue: The 2nd Internatonal Conference and Workshop on Mathematcal Analyss (ICWOMA 2016) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES

More information

Learning-Based Top-N Selection Query Evaluation over Relational Databases

Learning-Based Top-N Selection Query Evaluation over Relational Databases Learnng-Based Top-N Selecton Query Evaluaton over Relatonal Databases Lang Zhu *, Wey Meng ** * School of Mathematcs and Computer Scence, Hebe Unversty, Baodng, Hebe 071002, Chna, zhu@mal.hbu.edu.cn **

More information

Term Weighting Classification System Using the Chi-square Statistic for the Classification Subtask at NTCIR-6 Patent Retrieval Task

Term Weighting Classification System Using the Chi-square Statistic for the Classification Subtask at NTCIR-6 Patent Retrieval Task Proceedngs of NTCIR-6 Workshop Meetng, May 15-18, 2007, Tokyo, Japan Term Weghtng Classfcaton System Usng the Ch-square Statstc for the Classfcaton Subtask at NTCIR-6 Patent Retreval Task Kotaro Hashmoto

More information

Solutions to Programming Assignment Five Interpolation and Numerical Differentiation

Solutions to Programming Assignment Five Interpolation and Numerical Differentiation College of Engneerng and Coputer Scence Mechancal Engneerng Departent Mechancal Engneerng 309 Nuercal Analyss of Engneerng Systes Sprng 04 Nuber: 537 Instructor: Larry Caretto Solutons to Prograng Assgnent

More information

Lecture 5: Multilayer Perceptrons

Lecture 5: Multilayer Perceptrons Lecture 5: Multlayer Perceptrons Roger Grosse 1 Introducton So far, we ve only talked about lnear models: lnear regresson and lnear bnary classfers. We noted that there are functons that can t be represented

More information

Principal Component Inversion

Principal Component Inversion Prncpal Component Inverson Dr. A. Neumann, H. Krawczyk German Aerospace Centre DLR Remote Sensng Technology Insttute Marne Remote Sensng Prncpal Components - Propertes The Lnear Inverson Algorthm Optmsaton

More information

Positive Semi-definite Programming Localization in Wireless Sensor Networks

Positive Semi-definite Programming Localization in Wireless Sensor Networks Postve Sem-defnte Programmng Localzaton n Wreless Sensor etworks Shengdong Xe 1,, Jn Wang, Aqun Hu 1, Yunl Gu, Jang Xu, 1 School of Informaton Scence and Engneerng, Southeast Unversty, 10096, anjng Computer

More information

Radial Basis Functions

Radial Basis Functions Radal Bass Functons Mesh Reconstructon Input: pont cloud Output: water-tght manfold mesh Explct Connectvty estmaton Implct Sgned dstance functon estmaton Image from: Reconstructon and Representaton of

More information

High resolution 3D Tau-p transform by matching pursuit Weiping Cao* and Warren S. Ross, Shearwater GeoServices

High resolution 3D Tau-p transform by matching pursuit Weiping Cao* and Warren S. Ross, Shearwater GeoServices Hgh resoluton 3D Tau-p transform by matchng pursut Wepng Cao* and Warren S. Ross, Shearwater GeoServces Summary The 3D Tau-p transform s of vtal sgnfcance for processng sesmc data acqured wth modern wde

More information

UB at GeoCLEF Department of Geography Abstract

UB at GeoCLEF Department of Geography   Abstract UB at GeoCLEF 2006 Mguel E. Ruz (1), Stuart Shapro (2), June Abbas (1), Slva B. Southwck (1) and Davd Mark (3) State Unversty of New York at Buffalo (1) Department of Lbrary and Informaton Studes (2) Department

More information

S.P.H. : A SOLUTION TO AVOID USING EROSION CRITERION?

S.P.H. : A SOLUTION TO AVOID USING EROSION CRITERION? S.P.H. : A SOLUTION TO AVOID USING EROSION CRITERION? Célne GALLET ENSICA 1 place Emle Bloun 31056 TOULOUSE CEDEX e-mal :cgallet@ensca.fr Jean Luc LACOME DYNALIS Immeuble AEROPOLE - Bat 1 5, Avenue Albert

More information

Support Vector Machines

Support Vector Machines Support Vector Machnes Decson surface s a hyperplane (lne n 2D) n feature space (smlar to the Perceptron) Arguably, the most mportant recent dscovery n machne learnng In a nutshell: map the data to a predetermned

More information

Improved Methods for Lithography Model Calibration

Improved Methods for Lithography Model Calibration Improved Methods for Lthography Model Calbraton Chrs Mack www.lthoguru.com, Austn, Texas Abstract Lthography models, ncludng rgorous frst prncple models and fast approxmate models used for OPC, requre

More information

Parameter estimation for incomplete bivariate longitudinal data in clinical trials

Parameter estimation for incomplete bivariate longitudinal data in clinical trials Parameter estmaton for ncomplete bvarate longtudnal data n clncal trals Naum M. Khutoryansky Novo Nordsk Pharmaceutcals, Inc., Prnceton, NJ ABSTRACT Bvarate models are useful when analyzng longtudnal data

More information

Unsupervised Learning

Unsupervised Learning Pattern Recognton Lecture 8 Outlne Introducton Unsupervsed Learnng Parametrc VS Non-Parametrc Approach Mxture of Denstes Maxmum-Lkelhood Estmates Clusterng Prof. Danel Yeung School of Computer Scence and

More information

Hierarchical clustering for gene expression data analysis

Hierarchical clustering for gene expression data analysis Herarchcal clusterng for gene expresson data analyss Gorgo Valentn e-mal: valentn@ds.unm.t Clusterng of Mcroarray Data. Clusterng of gene expresson profles (rows) => dscovery of co-regulated and functonally

More information

The Codesign Challenge

The Codesign Challenge ECE 4530 Codesgn Challenge Fall 2007 Hardware/Software Codesgn The Codesgn Challenge Objectves In the codesgn challenge, your task s to accelerate a gven software reference mplementaton as fast as possble.

More information

Load-Balanced Anycast Routing

Load-Balanced Anycast Routing Load-Balanced Anycast Routng Chng-Yu Ln, Jung-Hua Lo, and Sy-Yen Kuo Department of Electrcal Engneerng atonal Tawan Unversty, Tape, Tawan sykuo@cc.ee.ntu.edu.tw Abstract For fault-tolerance and load-balance

More information

Skew Angle Estimation and Correction of Hand Written, Textual and Large areas of Non-Textual Document Images: A Novel Approach

Skew Angle Estimation and Correction of Hand Written, Textual and Large areas of Non-Textual Document Images: A Novel Approach Angle Estmaton and Correcton of Hand Wrtten, Textual and Large areas of Non-Textual Document Images: A Novel Approach D.R.Ramesh Babu Pyush M Kumat Mahesh D Dhannawat PES Insttute of Technology Research

More information

Outline. Discriminative classifiers for image recognition. Where in the World? A nearest neighbor recognition example 4/14/2011. CS 376 Lecture 22 1

Outline. Discriminative classifiers for image recognition. Where in the World? A nearest neighbor recognition example 4/14/2011. CS 376 Lecture 22 1 4/14/011 Outlne Dscrmnatve classfers for mage recognton Wednesday, Aprl 13 Krsten Grauman UT-Austn Last tme: wndow-based generc obect detecton basc ppelne face detecton wth boostng as case study Today:

More information

The Research of Ellipse Parameter Fitting Algorithm of Ultrasonic Imaging Logging in the Casing Hole

The Research of Ellipse Parameter Fitting Algorithm of Ultrasonic Imaging Logging in the Casing Hole Appled Mathematcs, 04, 5, 37-3 Publshed Onlne May 04 n ScRes. http://www.scrp.org/journal/am http://dx.do.org/0.436/am.04.584 The Research of Ellpse Parameter Fttng Algorthm of Ultrasonc Imagng Loggng

More information

6.854 Advanced Algorithms Petar Maymounkov Problem Set 11 (November 23, 2005) With: Benjamin Rossman, Oren Weimann, and Pouya Kheradpour

6.854 Advanced Algorithms Petar Maymounkov Problem Set 11 (November 23, 2005) With: Benjamin Rossman, Oren Weimann, and Pouya Kheradpour 6.854 Advanced Algorthms Petar Maymounkov Problem Set 11 (November 23, 2005) Wth: Benjamn Rossman, Oren Wemann, and Pouya Kheradpour Problem 1. We reduce vertex cover to MAX-SAT wth weghts, such that the

More information

12/2/2009. Announcements. Parametric / Non-parametric. Case-Based Reasoning. Nearest-Neighbor on Images. Nearest-Neighbor Classification

12/2/2009. Announcements. Parametric / Non-parametric. Case-Based Reasoning. Nearest-Neighbor on Images. Nearest-Neighbor Classification Introducton to Artfcal Intellgence V22.0472-001 Fall 2009 Lecture 24: Nearest-Neghbors & Support Vector Machnes Rob Fergus Dept of Computer Scence, Courant Insttute, NYU Sldes from Danel Yeung, John DeNero

More information

Probability Base Classification Technique: A Preliminary Study for Two Groups

Probability Base Classification Technique: A Preliminary Study for Two Groups Mathematcal Theory and Modelng ISSN 4-5804 (Paper) ISSN 5-05 (Onlne) www.ste.org Probablty Base Classfcaton Technque: A Prelmnary Study for Two Groups Frday Znzendoff Okwonu,* Abdul Rahman Othman. Department

More information

A Five-Point Subdivision Scheme with Two Parameters and a Four-Point Shape-Preserving Scheme

A Five-Point Subdivision Scheme with Two Parameters and a Four-Point Shape-Preserving Scheme Mathematcal and Computatonal Applcatons Artcle A Fve-Pont Subdvson Scheme wth Two Parameters and a Four-Pont Shape-Preservng Scheme Jeqng Tan,2, Bo Wang, * and Jun Sh School of Mathematcs, Hefe Unversty

More information

Optimizing Document Scoring for Query Retrieval

Optimizing Document Scoring for Query Retrieval Optmzng Document Scorng for Query Retreval Brent Ellwen baellwe@cs.stanford.edu Abstract The goal of ths project was to automate the process of tunng a document query engne. Specfcally, I used machne learnng

More information

Analysis of Continuous Beams in General

Analysis of Continuous Beams in General Analyss of Contnuous Beams n General Contnuous beams consdered here are prsmatc, rgdly connected to each beam segment and supported at varous ponts along the beam. onts are selected at ponts of support,

More information

Data Mining: Model Evaluation

Data Mining: Model Evaluation Data Mnng: Model Evaluaton Aprl 16, 2013 1 Issues: Evaluatng Classfcaton Methods Accurac classfer accurac: predctng class label predctor accurac: guessng value of predcted attrbutes Speed tme to construct

More information

LECTURE : MANIFOLD LEARNING

LECTURE : MANIFOLD LEARNING LECTURE : MANIFOLD LEARNING Rta Osadchy Some sldes are due to L.Saul, V. C. Raykar, N. Verma Topcs PCA MDS IsoMap LLE EgenMaps Done! Dmensonalty Reducton Data representaton Inputs are real-valued vectors

More information

Solving two-person zero-sum game by Matlab

Solving two-person zero-sum game by Matlab Appled Mechancs and Materals Onlne: 2011-02-02 ISSN: 1662-7482, Vols. 50-51, pp 262-265 do:10.4028/www.scentfc.net/amm.50-51.262 2011 Trans Tech Publcatons, Swtzerland Solvng two-person zero-sum game by

More information

Sum of Linear and Fractional Multiobjective Programming Problem under Fuzzy Rules Constraints

Sum of Linear and Fractional Multiobjective Programming Problem under Fuzzy Rules Constraints Australan Journal of Basc and Appled Scences, 2(4): 1204-1208, 2008 ISSN 1991-8178 Sum of Lnear and Fractonal Multobjectve Programmng Problem under Fuzzy Rules Constrants 1 2 Sanjay Jan and Kalash Lachhwan

More information

3D vector computer graphics

3D vector computer graphics 3D vector computer graphcs Paolo Varagnolo: freelance engneer Padova Aprl 2016 Prvate Practce ----------------------------------- 1. Introducton Vector 3D model representaton n computer graphcs requres

More information