{"id":8225,"date":"2019-03-05T01:55:08","date_gmt":"2019-03-05T00:55:08","guid":{"rendered":"http:\/\/www.sigterritoires.fr\/?p=8225"},"modified":"2019-03-05T08:21:51","modified_gmt":"2019-03-05T07:21:51","slug":"starting-with-arcgis-pro-create-your-own-source-of-elevations-3","status":"publish","type":"post","link":"https:\/\/www.sigterritoires.fr\/index.php\/en\/starting-with-arcgis-pro-create-your-own-source-of-elevations-3\/","title":{"rendered":"Starting with ArcGis Pro: create your own source of elevations ( 3)"},"content":{"rendered":"\n

To conclude the creation of our own source of elevation , started in\nprevious articles( Starting with ArcGis\nPro: create<\/a> your<\/a> own\nsource of elevations ( 1)<\/a> and (2)\n) we are\ngoing discuss the use of the raster calculator with ArCGis Pro.
\nWe have two different areas in our work surface: <\/p>\n\n\n\n

a terrestrial area <\/li><\/ul>\n\n\n\n
<\/p>\n\n\n\n
- a marine area<\/li><\/ul>\n\n\n\n\n\n\n\n
  $\"\"$ <\/figure>\n\n\n\n
  \n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\nWe\nhave finished calculating the marine area with the command Topo to raster: \n\n\n\n<\/p>\n\n\n\n
  $\"\"$ <\/figure>\n\n\n\n
  We will create a continuous DEM by following the following steps :
  \n1-cutting the marine DEM with the \u00a0\u00bb sea \u00a0\u00bb polygon
  \n2-cutting of the terrestrial DTM with the \u00a0\u00bb earth \u00a0\u00bb polygons
  \n3-filling of the empty cells eventually created in the boundary of both areas (\nland-sea ) <\/p>\n\n\n\n
  Cutting<\/strong> the<\/strong> marine<\/strong> DEM<\/strong> with the<\/strong> \u201csea<\/strong> \u00ab\u00a0<\/strong> polygon<\/strong> \n\nWe will select the polygon \u00ab\u00a0sea\u00a0\u00bb of\nthe polygon terrestrial-sea layer ( in Map tab of the ribbon \u00ab\u00a0Select\u201d\nbutton)
  \nIn the Geoprocessing window , click on Toolboxes and go to Data Management\nTools -> Raster -> Raster Processing -> Clip \n\n\n\n<\/p>\n\n\n\n
  $\"\"$ <\/figure>\n\n\n\n
  \n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\nYou\nmust indicate that you want to cut the raster with the layer terrestrial_sea AND\ncheck Use Input Features for geometry Clipping. This makes ArcGIS Pro to use,\nonly, the entities selected for cutting , in this specific case, the \u00a0\u00bb sea\n\u00a0\u00bb polygon. \n\n\n\n<\/p>\n\n\n\n
  $\"\"$ <\/figure>\n\n\n\n
  Cutting<\/strong>\nof the<\/strong>\nterrestrial<\/strong>\nDTM<\/strong>\nwith the<\/strong>\n\u201cearth<\/strong>\u201d<\/strong> polygons<\/strong>\n\nWe\nstart by loading the whole terrestrial area DTM in the view window\n\n\n\n<\/p>\n\n\n\n
  $\"\"$ <\/figure>\n\n\n\n
  We will perform the\nsame operation as before, but by using the \u201cterrestrial\u201d polygons.
  \nSince the \u00a0\u00bb sea \u00a0\u00bb polygon is always selected, just click on the terrestrial_sea layer in the legend window -> Select -> select switch to invert\nthe selection and get all \u00ab\u00a0earth\u00a0\u00bb polygons\nselected. \n\nThe\nRaster Clip command is executed again , but this time on the srtm raster . It just\nreplaces the name of the input raster.
  \nThe result is loaded\ninto the view : \n\n\n\n<\/p>\n\n\n\n
  $\"\"$ <\/figure>\n\n\n\n
  Now that we have both sections of the desired DTM, we need to join them in an unique raster. To that end we will use the ArcGis Pro Raster Calculator .\n\n\nIn\nthe Geoprocessing window, return to the Toolboxes and click Spatial Analyst\nTools -> Map Algebra -> Raster Calculator \n\n\n\n<\/p>\n\n\n\n
  $\"\"$ <\/figure>\n\n\n\n
  The formula is the\nfollowing :
  \nCon (<\/strong> IsNull<\/strong> (\u00ab\u00a0<\/strong> MntMer<\/strong>\n\u00ab\u00a0), \u00ab\u00a0srtm_36_031_Clip\u00a0\u00bb, \u00ab\u00a0MntMer\u00a0\u00bb)<\/strong> <\/p>\n\n\n\n
  In French , this\nmeans that, depending on the following condition (Con):\nif the pixel of the MntMer raster is empty ( IsNull (\u00a0\u00bb MntMer \u00ab\u00a0)), we will retain the\nvalue of the terrestrial raster , if not we will retain the value of the sea\nraster. \n\nThe\nonly problem with this formula\nis that some pixels on the terrestrial-sea boundary\nare empty on both input\nraster (the pixels that appear in white):\n \n\n\n\n<\/p>\n\n\n\n
  $\"\"$ <\/figure>\n\n\n\n
  It is enough to replace these null values by 0, which corresponds to the height of the coastline.\n<\/p>\n\n\n\n
  To this end , in the raster calculator you must type the formula
  \nCon (<\/strong> IsNull<\/strong> (\u00ab\u00a0Mntcommun1\u00a0\u00bb),\n0, \u00ab\u00a0Mntcommun1\u00a0\u00bb)<\/strong> \n\nThe\nresult is the following : \n\n\n\n<\/p>\n\n\n\n
  $\"\"$ <\/figure>\n\n\n\n
  Now we can click on\nthe Scene tab and edit the dtm used for the 3D\nview . \n\nBy\ndefault the Aerial image appears in\nthe dtm provided with ArcGis Pro: \n\n\n\n<\/p>\n\n\n\n
  $\"\"$ <\/figure>\n\n\n\n
  \n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\nBy clicking the ribbon, MAP tab, Add Preset button-> Ground and\npointing to dtm that we have just created, we get: \n\n\n\n<\/p>\n\n\n\n
  $\"\"$ <\/figure>\n\n\n\n
  \n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\nEven if it is subtle enough , we can see that the absolute\nflat surface of the marine part of the first picture is replaced by the actual seabed elevation.\n\n\n\n<\/p>\n","protected":false},"excerpt":{"rendered":"
  To conclude the creation of our own source of elevation , started in previous articles( Starting with ArcGis Pro: create your own source of elevations ( 1) and (2) ) we are going discuss the use…<\/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-8225","post","type-post","status-publish","format-standard","hentry","category-non-classe-en"],"aioseo_notices":[],"jetpack_featured_media_url":"","jetpack_shortlink":"https:\/\/wp.me\/p6XU0A-28F","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/8225","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=8225"}],"version-history":[{"count":0,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/8225\/revisions"}],"wp:attachment":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media?parent=8225"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/categories?post=8225"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/tags?post=8225"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}