{"id":60,"date":"2023-12-15T02:31:09","date_gmt":"2023-12-15T02:31:09","guid":{"rendered":"https:\/\/www.ccsgeo.com\/?p=60"},"modified":"2023-12-15T03:01:49","modified_gmt":"2023-12-15T03:01:49","slug":"using-sentinel-imagery-to-analyze-drummond-cos-coal-mining-operation-in-colombia","status":"publish","type":"post","link":"https:\/\/www.ccsgeo.com\/index.php\/2023\/12\/15\/using-sentinel-imagery-to-analyze-drummond-cos-coal-mining-operation-in-colombia\/","title":{"rendered":"Using Sentinel Imagery to Analyze Drummond Co’s Coal Mining Operation in Colombia"},"content":{"rendered":"\n

In a recent trip to Colombia, I had the opportunity to tour Puerto Drummond in Santa Marta. The port receives coal via train from their mines further south and sends out coal around the world. In order to provide a geographic context to this operation<\/a>, we can use satellite imagery to analyze the operation, from mine to port. <\/p>\n\n\n\n

Finding and processing satellite imagery however can be a challenge. In a previous post<\/a>, I discussed the cost of high resolution satellite imagery and how it can be inaccessible for independent analysis. Luckily there are other options, like Sentinel 2 imagery<\/a>. It is only 10m spatial resolution but that is sufficient for our purposes. And using ArcGIS Pro and the Sentinel 2 image service<\/a> from ArcGIS Online Living Atlas, you can query, visualize, and export imagery directly. <\/p>\n\n\n\n

The first step is to ensure you are signed in to ArcGIS Online through ArcGIS Pro. Most users will already be as ArcGIS Online will be their active licensing Portal. Then, through Catalog, click ‘Portal’, then ‘Living Atlas’ and search for ‘Sentinel 2’: <\/p>\n\n\n

\n
\"Catalog
Sentinel-2 Views image service from Living Atlas<\/figcaption><\/figure>\n<\/div>\n\n\n

Right click the layer and add it to a new map. You can then zoom to your area of interest. Because this is an image service coming from ArcGIS Online, there are many different processing templates available which can alter the rendering of the layer. In addition, it is important to remember that this is a mosaic dataset. So there are Overviews involved, which work to speed up rendering at different resolutions. Our job is to get to the source raster and export it so we can use it and manipulate as needed. <\/p>\n\n\n\n

My preferred method is to use the Select tool and use a rectangular selection to select my general area of interest. Then, through the Attribute pane, I can view the attributes for the individual rasters that make up the service for this area: <\/p>\n\n\n\n

\"Sentinel-2
Selection of individual images within the image service. <\/figcaption><\/figure>\n\n\n\n

You can see in the image above that the item with objectid 15494269 is the actual source image. All of the other images are overviews. You can then take this objectid and use it to set a Display Filter:
<\/p>\n\n\n

\n
\"Display
Definition Query specifying the image objectid. <\/figcaption><\/figure>\n<\/div>\n\n\n

The result is the actual underlying image: <\/p>\n\n\n

\n
\"Image
Result of definition query showing source image. <\/figcaption><\/figure>\n<\/div>\n\n\n

Unfortunately we can’t export the whole image as this would put unnecessary demands on the service. But we can export an image 5000 x 5000 pixels. The key part to remember is that we are looking at the default processing template, which is a Dynamic Range Adjustment (DRA). If you export the raster with this template applied you may see strange outputs. You will need to first change the template to Natural Colour before exporting (right-click image layer > Properties > Processing Templates): <\/p>\n\n\n\n

\"Specifying
Specifying the Natural Color Processing Template. <\/figcaption><\/figure>\n\n\n\n

It may not look as nice but don’t worry! We can enhance the symbology afterwards. Now, right-click the layer > Data > Export Raster. By default, the clipping extent will be too large and the cell-size won’t be correct. The red highlights warns you that you can’t export more than 5000 x 5000 pixels: <\/p>\n\n\n\n

\n
\"Export
Incorrect Raster size with small Cell size. <\/figcaption><\/figure>\n<\/div>\n\n
\n
\"Export
Correct Cell Size with updated Raster Size and a specified clipping geometry. <\/figcaption><\/figure>\n<\/div><\/div>\n\n\n\n

<\/p>\n\n\n\n

<\/p>\n\n\n\n

<\/p>\n\n\n\n

<\/p>\n\n\n\n

<\/p>\n\n\n\n

<\/p>\n\n\n\n

<\/p>\n\n\n\n

<\/p>\n\n\n\n

<\/p>\n\n\n\n

<\/p>\n\n\n\n

<\/p>\n\n\n\n

Change the cell size to 10 x 10 and modify your extent to ensure you have a correct image size. Then, click ‘Export’. You can now use your image for image analysis and to create maps.<\/p>\n\n\n\n

<\/p>\n\n\n\n

\"Puerto
Final map layout using Sentinel-2 Imagery exported from ArcGIS Online’s living Atlas<\/figcaption><\/figure>\n","protected":false},"excerpt":{"rendered":"

In a recent trip to Colombia, I had the opportunity to tour Puerto Drummond in Santa Marta. The port receives coal via train from their mines further south and sends out coal around the world. In order to provide a geographic context to this operation, we can use satellite imagery to analyze the operation, from […]<\/p>\n","protected":false},"author":2,"featured_media":85,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"_crdt_document":"","_monsterinsights_skip_tracking":false,"_monsterinsights_sitenote_active":false,"_monsterinsights_sitenote_note":"","_monsterinsights_sitenote_category":0,"footnotes":""},"categories":[10,11,5,12],"tags":[],"class_list":["post-60","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-arcgis-pro","category-esri","category-satellite-imagery","category-workflow"],"_links":{"self":[{"href":"https:\/\/www.ccsgeo.com\/index.php\/wp-json\/wp\/v2\/posts\/60","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.ccsgeo.com\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.ccsgeo.com\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.ccsgeo.com\/index.php\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/www.ccsgeo.com\/index.php\/wp-json\/wp\/v2\/comments?post=60"}],"version-history":[{"count":20,"href":"https:\/\/www.ccsgeo.com\/index.php\/wp-json\/wp\/v2\/posts\/60\/revisions"}],"predecessor-version":[{"id":89,"href":"https:\/\/www.ccsgeo.com\/index.php\/wp-json\/wp\/v2\/posts\/60\/revisions\/89"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.ccsgeo.com\/index.php\/wp-json\/wp\/v2\/media\/85"}],"wp:attachment":[{"href":"https:\/\/www.ccsgeo.com\/index.php\/wp-json\/wp\/v2\/media?parent=60"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.ccsgeo.com\/index.php\/wp-json\/wp\/v2\/categories?post=60"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.ccsgeo.com\/index.php\/wp-json\/wp\/v2\/tags?post=60"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}