Identify raster cell

View on GitHubSample viewer app

Get the cell value of a local raster at the tapped location and display the result in a callout.

screenshot

Use case

You may want to identify a raster layer to get its exact cell value in case the approximate value conveyed by its symbology is not sufficient. The information available for the raster cell depends on the type of raster layer being identified. For example, a 3-band satellite or aerial image might provide 8-bit RGB values, whereas a digital elevation model (DEM) would provide floating point z values. By identifying a raster cell of a DEM, you can retrieve the precise elevation of a location.

How to use the sample

Tap on the raster or press on the raster, hold, and move around the raster to identify it and see the raster cell attributes information displayed in a callout.

How it works

  1. Connect to signals emmited by a mouseClicked, mousePressedAndHeld, and mouseMoved on the MapView.
  2. On tapped or pressed, held, and dragged:
    • Call identifyLayerAsync(...) passing in the raster layer, screen point, tolerance, whether to return popups only, and maximum number of results per layer.
    • Connect to the identifyLayerCompleted to get the result of the identify and then get the GeoElement from the layer result and get any RasterCells from them.
    • Create a callout at the calculated map point and populate the callout content with text from the RasterCell attributes.
    • Show the callout.

Relevant API

  • GeoView::identifyLayerAsync(...)
  • IdentifyLayerResult
  • RasterCell
  • RasterCell::attributes
  • RasterLayer

Offline data

To set up the sample's offline data, see the Use offline data in the samples section of the Qt Samples repository overview.

Link Local Location
South Africa data <userhome>/ArcGIS/Runtime/Data/raster/SA_EVI_8Day_03May20

About the data

The data shown is an NDVI classification derived from MODIS imagery between 27 Apr 2020 and 4 May 2020. It comes from the NASA Worldview application. In a normalized difference vegetation index, or NDVI, values range between -1 and +1 with the positive end of the spectrum showing green vegetation.

Tags

band, cell, cell value, continuous, discrete, identify, pixel, pixel value, raster

Sample Code

IdentifyRasterCell.cppIdentifyRasterCell.cppIdentifyRasterCell.hIdentifyRasterCell.qml
Use dark colors for code blocksCopy
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
// [WriteFile Name=IdentifyRasterCell, Category=Layers]
// [Legal]
// Copyright 2020 Esri.
//
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
// http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
// See the License for the specific language governing permissions and
// limitations under the License.
// [Legal]

#ifdef PCH_BUILD
#include "pch.hpp"
#endif // PCH_BUILD

// sample headers
#include "IdentifyRasterCell.h"

// ArcGIS Maps SDK headers
#include "AttributeListModel.h"
#include "CalloutData.h"
#include "Envelope.h"
#include "Error.h"
#include "IdentifyLayerResult.h"
#include "LayerListModel.h"
#include "Map.h"
#include "MapQuickView.h"
#include "MapTypes.h"
#include "Raster.h"
#include "RasterCell.h"
#include "RasterLayer.h"

// Qt headers
#include <QFuture>
#include <QStandardPaths>
#include <QString>
#include <QUuid>

// STL headers
#include <memory>

using namespace Esri::ArcGISRuntime;

// helper method to get cross platform data path
namespace
{
QString defaultDataPath()
{
  QString dataPath;

#ifdef Q_OS_IOS
  dataPath = QStandardPaths::writableLocation(QStandardPaths::DocumentsLocation);
#else
  dataPath = QStandardPaths::writableLocation(QStandardPaths::HomeLocation);
#endif

  return dataPath;
}
} // namespace

IdentifyRasterCell::IdentifyRasterCell(QObject* parent /* = nullptr */):
  QObject(parent),
  m_map(new Map(BasemapStyle::ArcGISOceans, this))
{
  // initialize the raster layer
  const QString filepath = defaultDataPath() + "/ArcGIS/Runtime/Data/raster/SA_EVI_8Day_03May20/SA_EVI_8Day_03May20.tif";
  Raster* raster = new Raster(filepath, this);
  m_rasterLayer = new RasterLayer(raster, this);

  m_map->operationalLayers()->append(m_rasterLayer);
}

IdentifyRasterCell::~IdentifyRasterCell() = default;

void IdentifyRasterCell::init()
{
  // Register the map view for QML
  qmlRegisterType<MapQuickView>("Esri.Samples", 1, 0, "MapView");
  qmlRegisterType<IdentifyRasterCell>("Esri.Samples", 1, 0, "IdentifyRasterCellSample");
}

MapQuickView* IdentifyRasterCell::mapView() const
{
  return m_mapView;
}

// Set the view (created in QML)
void IdentifyRasterCell::setMapView(MapQuickView* mapView)
{
  if (!mapView || mapView == m_mapView)
    return;

  m_mapView = mapView;
  m_mapView->setMap(m_map);

  connectSignals();

  emit mapViewChanged();
}

void IdentifyRasterCell::connectSignals()
{
  connect(m_rasterLayer, &Layer::doneLoading, this, [this](const Error& error)
  {
    if (!error.isEmpty())
    {
      qDebug() << error.message() << error.additionalMessage();
      return;
    }

    m_mapView->setViewpointGeometryAsync(m_rasterLayer->fullExtent());
  });

  connect(m_mapView, &MapQuickView::mouseClicked, this, [this](const QMouseEvent& e)
  {
    m_mapView->identifyLayerAsync(m_rasterLayer, e.position(), 10, false, 1).then(this, [this](IdentifyLayerResult* rawIdentifyResult)
    {
      onIdentifyLayerCompleted_(rawIdentifyResult);
    });
    m_clickedPoint = m_mapView->screenToLocation(e.position().x(), e.position().y());
  });

  connect(m_mapView, &MapQuickView::mousePressedAndHeld, this, [this](const QMouseEvent& e)
  {
    m_mapView->identifyLayerAsync(m_rasterLayer, e.position(), 10, false, 1).then(this, [this](IdentifyLayerResult* rawIdentifyResult)
    {
      onIdentifyLayerCompleted_(rawIdentifyResult);
    });
    m_clickedPoint = m_mapView->screenToLocation(e.position().x(), e.position().y());
    m_mousePressed = true;
  });

  connect(m_mapView, &MapQuickView::mouseReleased, this, [this](QMouseEvent&)
  {
    m_mousePressed = false;
  });

  // if mouse is moved after press-and-hold, then update callout on-the-fly
  connect(m_mapView, &MapQuickView::mouseMoved, this, [this](const QMouseEvent& e)
  {
    if (m_mousePressed)
    {
      m_mapView->identifyLayerAsync(m_rasterLayer, e.position(), 10, false, 1).then(this, [this](IdentifyLayerResult* rawIdentifyResult)
      {
        onIdentifyLayerCompleted_(rawIdentifyResult);
      });
      m_clickedPoint = m_mapView->screenToLocation(e.position().x(), e.position().y());
    }
  });
}

void IdentifyRasterCell::onIdentifyLayerCompleted_(IdentifyLayerResult* rawIdentifyResult)
{
  const auto identifyResult = std::unique_ptr<IdentifyLayerResult>(rawIdentifyResult);
  const auto elements = identifyResult->geoElements();
  for (GeoElement* geoElement : elements)
  {
    if (RasterCell* rasterCell = dynamic_cast<RasterCell*>(geoElement))
    {
      QString calloutString;
      const AttributeListModel* attributes = rasterCell->attributes();
      const QStringList attributeNames = rasterCell->attributes()->attributeNames();

      for (int i = 0; i < attributeNames.size(); ++i)
      {
        const QString value = QVariant((*attributes)[attributeNames[i]]).toString();
        calloutString.append(attributeNames[i] + ": " + value + "\n");
      }

      const double xPoint = rasterCell->geometry().extent().xMin();
      const double yPoint = rasterCell->geometry().extent().yMin();

      calloutString.append("X: " + QString::number(xPoint, 'f', 2) + " Y: " + QString::number(yPoint, 'f', 2));

      m_mapView->calloutData()->setLocation(m_clickedPoint);
      m_mapView->calloutData()->setDetail(calloutString);
      m_mapView->calloutData()->setVisible(true);
      m_calloutData = m_mapView->calloutData();
      emit calloutDataChanged();
    }
  }
}

CalloutData* IdentifyRasterCell::calloutData() const
{
  return m_calloutData;
}

Your browser is no longer supported. Please upgrade your browser for the best experience. See our browser deprecation post for more details.