Using the Time-Enabled LandsatLoo LandsatLook k Service with GIS Clients
1
Table of Contents Document History ......................................................................................................................................... 3 Introduction .................................................................................................................................................. 4 Using ArcMap (OGC WMS-T) .................................. ................. ................................... ................................... ................................... ................................... ................................ ............... 5 Using QGIS (OGC WMS-T) ............................................................................................................................. 8 Using GDAL (OGC WMS-T) .......................................................................................................................... 15 Using Global Mapper (OGC WMS-T) ........................................................................................................... 16 Using a Web Browser (OGC WMS-T) .......................................................................................................... 21 Appendix A – A – Sentinel Sentinel 2 Image Service ....................................................................................................... 22
2
Document History Date 3/9/2016 3/15/2016 3/23/2016 3/23/2016 3/23/2016
Author Layth Layth Layth Layth Layth
3/23/2016 1/04/2017
Layth Layth
Section(s) All Using ArcMap (REST API) Using ArcMap (REST API) Using QGIS (WMS-T API) Using Global Mappers (WMS-T API) Using GDAL (WMS-T API) All
3
Description Document created Step 4 added a note for ArcGIS v10.2.2 Section removed at customer’s request Sections added at customer’s request Sections added at customer’s request Sections added at customer’s request Updated all URL’s to use HTTPS instead of HTTP. Also included some edits from peer review.
Introduction The LandsatLook image service has been updated to allow filtering by time using the native ArcGIS Representational state transfer (REST) and Simple Object Access Protocol (SOAP) services as well as OpenGIS® Web Map Service Interface Standard-Time (OGC WMS-T) standards. This document will provide examples of how to interact with the image service without using the provided map viewer (https://landsatlook.usgs.gov/). This capability will allow a user to use their GIS Desktop tool or create a different map viewer that uses this image service. The Landsat scenes are only available at map scales below 1:2M. In order to orient the map over the user’s desired location it is recommended that a basemap covering the area of interest be loaded first. An example is ESRI’s World Street Map available from https://services.arcgisonline.com/arcgis/rest/services/ESRI_StreetMap_World_2D/MapServer. These steps have been tested with ArcGIS Desktop v10.3.1, Quantum GI S v2.8 (Wien), Global Mapper v17, Geospatial Data Abstraction Library (GDAL) v1.11.0. The web browser example should work with any web browser.
4
Using ArcMap (OGC WMS-T) 1. Open ArcMap 2. Add the ESRI Street Map service: http://services.arcgisonline.com/arcgis/rest/services/ESRI_StreetMap_World_2D/MapServer/W MTS/1.0.0/WMTSCapabilities.xml.
3.
Add the LandsatLook WMS service https://landsatlook.usgs.gov/arcgis/rest/services/LandsatLook/ImageServer. This will enable the Time Slider button on the ArcMap.
4. Click the Time Slider button on the “Tools” menu to display the Time Slider window.
5
5. Using the ESRI Street Map layer, orient yourself to your desired area.
6. Ensure you have zoomed to a scale of 1:2M or less.
6
7. Use the Time Slider control to find the time you are interested in. Note: Landsat data is collected for each path every 16 days, so you will have to adjust your t ime control forward and backward until you find a day where one of the Landsat satellites passed over your area of interest.
7
Using QGIS (OGC WMS-T) 1. Open QGIS Desktop. 2. In the Browser pane, right-click WMS and choose New Connection…
8
3. Add the ESRI Streets WMTS service with the following URL: https://services.arcgisonline.com/arcgis/rest/services/ESRI_StreetMap_World_2D/MapServer/ WMTS/1.0.0/WMTSCapabilities.xml.
4. Click OK. 5. In the Browser pane, right-click WMS and choose New Connection… 6. Add the LandsatLook WMS service with the following URL: https://landsatlook.usgs.gov/arcgis/services/LandsatLook/ImageServer/WMSServer. 7. On the Main Menu, choose Plugins | Manage and Install Plugins… 8. Search All plugins for TimeManager.
9
9. Install TimeManager from Plugins.
10. Choose Add WMS/MWTS Layer…
11. Choose the ESRI Streets layer.
12. Choose Connect. 13. Choose the one available TileSet.
14. Choose Add and then Close.
10
15. Zoom to your area of interest. Note: QGIS and ESRI scale calculations do not seem to match. It is strongly recommended you select a scale BELOW 1:1,000,000 in QGIS.
16. Choose Add WMS/MWTS Layer…
17. Choose the LandsatLook layer.
11
18. Choose the PNG format. Note: Testing has shown that TIFF will not return correctly transparent images.
19. Choose Add | Close. 20. In the Time Manager window, choose Settings. 21. In the Time Manager Settings Window, choose “Add raster”
12
22. Choose the LandsatLook layer. Enter a start and end date in the format “YYYY-MMDDT00:00:00Z”. For example: 2016-03-01T00:00:00Z.
23. Click OK 24. Click OK to close the Time Manager Settings window.
13
25. Set the Time Manager time interval to 1 Days.
26. You can now enter a date directly in the Time frame Start window.
27. Or you can use the arrows to move through time.
Note: Depending on network performance you may want to increase the timeout from the default 6 seconds to something longer, such as 20 seconds. Access this setting by choosing Settings | Options… | Network from the Main Menu.
14
Using GDAL (OGC WMS-T) 1. Use gdal_translate to retrieve images from the WMS map service. A properly formed WMS request is used. For example: gdal_translate -of PNG -outsize 800 600 "WMS:
https://landsatlook.usgs.gov/arcgis/services/LandsatLook/ImageServer/WMSServer?SERVICE=W MS&VERSION=1.3.0&REQUEST=GetMap&CRS=CRS:84&BBOX=83.61338722172151,27.060066392463657,79.430294773471886,29.748857329051877&WIDTH=879&HEIGHT=565&LAYERS=0&STYLES=&F ORMAT=image/png&BGCOLOR=0xFEFFFF&TRANSPARENT=TRUE&TIME=2000-01-01/2000-0130" test.png 2. Ensure the output format specified by the –of parameter matches the format requested via WMS. For example “-of PNG” and “FORMAT=image/png”. 3. Ensure the bounding box falls within the scale of 1:2,000,000 or below. 4. Time is adjusted using the WMS TIME element. In the example, all scenes from January 2000 are returned. A single date can be returned using the format: TIME=2000-01-03 as an example.
15
Using Global Mapper (OGC WMS-T)
1. Open Global Mapper. 2. Click the Download Online Data button.
3. Click the Add New Source button.
4. Select WMTS Source Type
5. Enter the ESRI Street Map WMTS service: https://services.arcgisonline.com/arcgis/rest/services/ESRI_StreetMap_World_2D/MapServer/ WMTS/1.0.0/WMTSCapabilities.xml 6. Click the Get List of Available Data Layers button.
16
7. Choose ESRI_StreetMap_World_2D.
8. Click OK 9. Add to your desired Category. 10. Highlight the Street Map layer and c lick the Connect button to add it to the map.
17
11. Zoom to your area of interest. Be sure the scale is below 1:2,000,000.
12. Click the Download Online Data button.
13. Click the Add New Source button.
14. Select WMS Source Type
18
15. Enter the LandsatLook WMS Service: https://landsatlook.usgs.gov/arcgis/services/LandsatLook/ImageServer/WMSServer?request=G etCapabilities&service=WMS 16. In the Service Name field, enter the desired time interval such as: &TIME=2000-01-01/2000-0130
17. Give the WMS-T connection a meaningful name and choose a c ategory:
19
18. Click the Connect button to add the layer to the map. Note: Load times for GlobalMapper seem to be longer than other clients. It will load in a few seconds.
19. Other time periods can be loaded by c reating new connections (Start with Ste p 12 in this section).
20
Using a Web Browser (OGC WMS-T) 1. Open a Web Browser (i.e. Google Chrome) 2. Use the following URL as an example: https://igskmncnvs557/arcgis/services/LandsatLook/ImageServer/WMSServer?SERVICE=WMS& VERSION=1.3.0&REQUEST=GetMap&CRS=CRS:84&BBOX=83.61338722172151,27.060066392463657,79.430294773471886,29.748857329051877&WIDTH=879&HEIGHT=565&LAYERS=0&STYLES=&F ORMAT=image/png&BGCOLOR=0xFEFFFF&TRANSPARENT=TRUE&TIME=2000-01-01/2000-01-30
Notice the TIME= element. In this case, TIME=2000-01-01/2000-01-30 will display all images that were acquired between January 1st, 2000 and January 30 th, 2000. Also notice that the BBOX= element must also take into account the maximum scale of 1:2M. If the BBOX is too large, you will get an empty image.
21
Appendix A – Sentinel-2 Image Service The same steps can be used to access Sentinel-2 data service. The WMS link for Sentinel 2 is https://landsatlook.usgs.gov/arcgis/services/Sentinel2/ImageServer/WMSServer The ArcGIS REST link is https://landsatlook.usgs.gov/arcgis/rest/services/Sentinel2/ImageServer
22
