-
-
Notifications
You must be signed in to change notification settings - Fork 3k
/
qgsmaplayerset.cpp
106 lines (94 loc) · 3.74 KB
/
qgsmaplayerset.cpp
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
/***************************************************************************
qgsmaplayerset.cpp - holds a set of layers
----------------------
begin : January 2006
copyright : (C) 2006 by Martin Dobias
email : wonder.sk at gmail dot com
***************************************************************************
* *
* This program is free software; you can redistribute it and/or modify *
* it under the terms of the GNU General Public License as published by *
* the Free Software Foundation; either version 2 of the License, or *
* (at your option) any later version. *
* *
***************************************************************************/
/* $Id$ */
#include "qgslogger.h"
#include "qgsmaplayerset.h"
#include "qgsmaplayerregistry.h"
#include "qgsproject.h"
#include <string>
void QgsMapLayerSet::setLayerSet(const std::deque<QString>& layers)
{
mLayerSet = layers;
updateFullExtent();
}
void QgsMapLayerSet::updateFullExtent()
{
QgsDebugMsg("QgsMapLayerSet::updateFullExtent() called !");
QgsMapLayerRegistry* registry = QgsMapLayerRegistry::instance();
bool projectionsEnabled = (QgsProject::instance()->readNumEntry("SpatialRefSys","/ProjectionsEnabled",0)!=0);
// reset the map canvas extent since the extent may now be smaller
// We can't use a constructor since QgsRect normalizes the rectangle upon construction
mFullExtent.setMinimal();
// iterate through the map layers and test each layers extent
// against the current min and max values
std::deque<QString>::iterator it = mLayerSet.begin();
while(it != mLayerSet.end())
{
QgsMapLayer * lyr = registry->mapLayer(*it);
if (lyr == NULL)
{
QgsLogger::warning("WARNING: layer '" + (*it) + "' not found in map layer registry!");
}
else
{
QgsDebugMsg("Updating extent using " + lyr->name());
QgsDebugMsg("Input extent: " + lyr->extent().stringRep());
// Layer extents are stored in the coordinate system (CS) of the
// layer. The extent must be projected to the canvas CS prior to passing
// on to the updateFullExtent function
if (projectionsEnabled)
{
try
{
if ( ! lyr->coordinateTransform() )
throw QgsCsException( std::string("NO COORDINATE TRANSFORM FOUND FOR LAYER") );
mFullExtent.unionRect(lyr->coordinateTransform()->transformBoundingBox(lyr->extent()));
}
catch (QgsCsException &cse)
{
QgsLogger::warning("Transform error caught in " + QString(__FILE__) + " line " +\
QString::number(__LINE__) + QString(cse.what()));
}
}
else
{
mFullExtent.unionRect(lyr->extent());
}
}
it++;
}
if (mFullExtent.width() == 0.0 || mFullExtent.height() == 0.0)
{
// If all of the features are at the one point, buffer the
// rectangle a bit. If they are all at zero, do something a bit
// more crude.
if (mFullExtent.xMin() == 0.0 && mFullExtent.xMax() == 0.0 &&
mFullExtent.yMin() == 0.0 && mFullExtent.yMax() == 0.0)
{
mFullExtent.set(-1.0, -1.0, 1.0, 1.0);
}
else
{
const double padFactor = 1e-8;
double widthPad = mFullExtent.xMin() * padFactor;
double heightPad = mFullExtent.yMin() * padFactor;
double xmin = mFullExtent.xMin() - widthPad;
double xmax = mFullExtent.xMax() + widthPad;
double ymin = mFullExtent.yMin() - heightPad;
double ymax = mFullExtent.yMax() + heightPad;
mFullExtent.set(xmin, ymin, xmax, ymax);
}
}
}