﻿{"id":8025,"date":"2019-02-05T05:48:14","date_gmt":"2019-02-05T04:48:14","guid":{"rendered":"http:\/\/www.sigterritoires.fr\/?p=8025"},"modified":"2019-02-05T08:40:09","modified_gmt":"2019-02-05T07:40:09","slug":"comparative-analysis-of-interpolation-methods-to-generate-the-dem-digital-terrain-model","status":"publish","type":"post","link":"https:\/\/www.sigterritoires.fr\/index.php\/en\/comparative-analysis-of-interpolation-methods-to-generate-the-dem-digital-terrain-model\/","title":{"rendered":"Comparative analysis of interpolation methods to generate the DEM (Digital Terrain Model)"},"content":{"rendered":"\n<p>As you must have already reckoned, we are\nparticularly interested in the transition of the territorial data to 3D.\nThe ArcGis Pro output with two side-by-side windows (2D\nwindow ArcMap type and 3D window ArcScene type) is a novelty that must be\nstudied. In the long run it will become either the standard for most\nGIS, or a vague memory of a failed attempt. This\nwill depend largely on you, whether you find that it adds a plus to your job,\nor that it is, just, another gadget you can do without. <br>\nTherefore,\nyou have to try it, but for that, of course, you must have 3D displayable data.\nThe first being the DEM (Digital Terrain Model), where all\nyour other data will be displayed. <br>\nTo\ngenerate a DEM you have to start with vector data (points, lines, surfaces)\ncontaining information about height (elevation), which one interpolates to have\na continuous XYZ surface. <br>\nThe\naccuracy of the generated terrain model depends on the interpolation mechanism\nused. It is, therefore, necessary to study the comparative\nperformance of different methods in this context. We have\ncompared the general interpolation techniques, namely the inverse of the\nweighted distance (IDW), the kriging, the ANUDEM method (topo to raster in\nArcGis), the natural neighbour, and the Spline method. <\/p>\n\n\n\n<!--more-->\n\n\n\n<p><strong>The different types of\ninterpolation methods.<\/strong> <\/p>\n\n\n\n<p>Different interpolation methods applied to the same data can\nproduce different results. Therefore, it is necessary to\nevaluate the comparative relevance of these methods. <br>\nThe\ninterpolation methods are based on the principle of spatial autocorrelation,\nwhich assumes that the closer the points, the more similar they are.\nIn the literature you will find many methods of interpolation,\nthey\nare generally classified into two categories: local and global methods.\n<br>\nLocal\nmethods predict the value of a point based on the values \u200b\u200bof points in the\nneighborhood. The most used local methods\nare: <\/p>\n\n\n\n<p>\u0095\nthe inverse of the weighted distance (IDW) method, <\/p>\n\n\n\n<p>\u0095\nthe local polynomial method, <\/p>\n\n\n\n<p>\u0095\nthe natural neighbour method (NN), and <\/p>\n\n\n\n<p>\u0095\nradial basic methods (Spline). <\/p>\n\n\n\n<p>On the other hand, global interpolation methods, such as\npolynomial interpolation functions, use all available sampling points to\ngenerate forecasts for a particular point. These\nmethods facilitate the evaluation and elimination of global phenomena (such as\ntrend) in physical data. <br>\nAll\nthese methods are also called deterministic as opposed to geostatistical\nmethods. <\/p>\n\n\n\n<p><strong>Kriging<\/strong>\n<\/p>\n\n\n\n<p>It is a geostatistical interpolation method that uses a\nvariogram ( analysis of the variability of the data according to the distance\nbetween them). The variogram depends on the\nspatial distribution of the data rather than the actual values. When\napplying the kriging method one can see results for input points different from\nthe input value. <\/p>\n\n\n\n<p><strong>The IDW method<\/strong>\n<\/p>\n\n\n\n<p>It is a local deterministic interpolation technique that\ncalculates the value of a point by averaging the values \u200b\u200bof the neighbouring points\nweighted by the inverse of the distance at the calculated point: the closer the\npoints, the more the affected weighting is strong. <br>\nIt\nconsiders that the points closer to the location to be calculated will have\nmore influence. <\/p>\n\n\n\n<p><strong>The natural neighbour method<\/strong>\n<\/p>\n\n\n\n<p>This method searches for the subset of samples closest to a\npoint and applies a weighting based on the area where they are located.\nIt is a local deterministic method and the interpolated\nheights are necessarily within the range of values \u200b\u200bused. It\ndoes not produce peaks, pits, ridges or valleys that are not already present in\nthe input samples and adapts locally to the input data structure.\nIt does not require any configuration by the user and works\nequally well for data distributed regularly as well as irregularly.\n<\/p>\n\n\n\n<p><strong>The Spline interpolation\nmethod<\/strong> <\/p>\n\n\n\n<p>This method uses a mathematical function to minimize the\ncurvature of the surface and produces a smooth surface that matches exactly the\nentry points. <\/p>\n\n\n\n<p><strong>The ANUDEM method (topo to\nraster of ArcGis 3D Analyst)<\/strong> <\/p>\n\n\n\n<p>This method uses an interpolation technique specifically\ndesigned to create a surface that best represents a natural drainage surface\nand preserves both ridge lines and stream networks. <\/p>\n\n\n\n<p><strong>Results of the comparison\nof interpolation methods<\/strong> \n\nWe compared the different\nmethods with the IGN repository for several projects, calculating the mean\nsquare error. <br>\nIn\norder to study the sensitivity of the interpolation methods according to the\nnature of the terrain, we grouped the results according to three types of\nslope: zones with little slope, zones of steep slope, and zones having a\nmixture of the two. The results are presented in\nqualitative form (stars) but it can be said that the best results have an EQM\nof the order of 50 cm and the least good an EQM of the order of 2m.\n\n\n\n\n<\/p>\n\n\n\n<figure class=\"wp-block-image\"><img data-recalc-dims=\"1\" loading=\"lazy\" decoding=\"async\" width=\"525\" height=\"77\" data-attachment-id=\"8026\" data-permalink=\"https:\/\/www.sigterritoires.fr\/index.php\/en\/comparative-analysis-of-interpolation-methods-to-generate-the-dem-digital-terrain-model\/51-3\/\" data-orig-file=\"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2019\/02\/51.png?fit=525%2C77&amp;ssl=1\" data-orig-size=\"525,77\" data-comments-opened=\"1\" data-image-meta=\"{&quot;aperture&quot;:&quot;0&quot;,&quot;credit&quot;:&quot;&quot;,&quot;camera&quot;:&quot;&quot;,&quot;caption&quot;:&quot;&quot;,&quot;created_timestamp&quot;:&quot;0&quot;,&quot;copyright&quot;:&quot;&quot;,&quot;focal_length&quot;:&quot;0&quot;,&quot;iso&quot;:&quot;0&quot;,&quot;shutter_speed&quot;:&quot;0&quot;,&quot;title&quot;:&quot;&quot;,&quot;orientation&quot;:&quot;0&quot;}\" data-image-title=\"51\" data-image-description=\"\" data-image-caption=\"\" data-large-file=\"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2019\/02\/51.png?fit=525%2C77&amp;ssl=1\" src=\"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2019\/02\/51.png?resize=525%2C77&#038;ssl=1\" alt=\"\" class=\"wp-image-8026\" srcset=\"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2019\/02\/51.png?w=525&amp;ssl=1 525w, https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2019\/02\/51.png?resize=300%2C44&amp;ssl=1 300w\" sizes=\"auto, (max-width: 525px) 100vw, 525px\" \/><\/figure>\n\n\n\n<p>In general, we can say that <strong>IDW<\/strong> methods\nand\n<strong>Kriging<\/strong> adapt quite well, whatever the variations of\nterrain. Other methods are generally more sensitive to variations.\nThe <strong>ANUDEM<\/strong> method has excellent performance when it\nwas a question of the calculation of ridges and areas of flow stream ..\n<br>\nIn\nlow slope areas, <strong>Kriging<\/strong> and the <strong>Natural Neighbor<\/strong>\ngive very good results, and we advise to adopt them in this case.\nIn areas of steep slope, the average differences are less\nmarked and must analyze them more on a case-by-case basis. The <strong>Natural\nNeighbor is<\/strong> best applied on small areas of study. For\nthe calculation of the flow zones the <strong>ANUDEM<\/strong> method gives\nthe best results. <br>\nThe\n<strong>natural neighbour<\/strong> method showed almost optimal values \u200b\u200bon smooth\nsurfaces. <br>\n<strong>Spline-<\/strong> based methods fit a minimum curved surface\nthrough the entry points. It preserves trends in the\nsample data and adapts to rapid changes in gradient or slope. <br>\nIn\nall cases , <strong>kriging<\/strong> gives good results, even for steep declivity areas\nas well as for areas with both steep and steep slopes. This\nmethod takes into account the auto-correlation structures of the heights of the\nzone, in order to define the optimal weights. But,\nin return, the method needs a qualified user, with good knowledge of\ngeostatistics. <\/p>\n","protected":false},"excerpt":{"rendered":"<p>As you must have already reckoned, we are particularly interested in the transition of the territorial data to 3D. The ArcGis Pro output with two side-by-side windows (2D window ArcMap type and 3D window ArcScene type)&hellip;<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"give_campaign_id":0,"_bbp_topic_count":0,"_bbp_reply_count":0,"_bbp_total_topic_count":0,"_bbp_total_reply_count":0,"_bbp_voice_count":0,"_bbp_anonymous_reply_count":0,"_bbp_topic_count_hidden":0,"_bbp_reply_count_hidden":0,"_bbp_forum_subforum_count":0,"sfsi_plus_gutenberg_text_before_share":"","sfsi_plus_gutenberg_show_text_before_share":"","sfsi_plus_gutenberg_icon_type":"","sfsi_plus_gutenberg_icon_alignemt":"","sfsi_plus_gutenburg_max_per_row":"","_monsterinsights_skip_tracking":false,"_monsterinsights_sitenote_active":false,"_monsterinsights_sitenote_note":"","_monsterinsights_sitenote_category":0,"jetpack_post_was_ever_published":false,"_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_jetpack_memberships_contains_paid_content":false,"footnotes":""},"categories":[1260],"tags":[],"class_list":["post-8025","post","type-post","status-publish","format-standard","hentry","category-non-classe-en"],"aioseo_notices":[],"jetpack_featured_media_url":"","jetpack_shortlink":"https:\/\/wp.me\/p6XU0A-25r","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/8025","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/comments?post=8025"}],"version-history":[{"count":0,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/8025\/revisions"}],"wp:attachment":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media?parent=8025"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/categories?post=8025"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/tags?post=8025"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}