Source code for tuiview.coordinatemgr

This module contains the CoordManager class, which manages 
the relationship between coordinates for display and coordinates 
in the raster
# This file is part of 'TuiView' - a simple Raster viewer
# Copyright (C) 2012  Sam Gillingham
# 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.
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# GNU General Public License for more details.
# You should have received a copy of the GNU General Public License
# along with this program; if not, write to the Free Software
# Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.

from __future__ import division
import math

[docs]class CoordManager(object): """ base class for the layer's coordmgr instance. Derived separately for vector and raster layers. """ def __init__(self): # The size of the display window, in display coords self.dspWidth = None self.dspHeight = None
[docs] def setDisplaySize(self, width, height): """ Set display size in units of display coordinate system """ self.dspWidth = width self.dspHeight = height
[docs] def getWorldExtent(self): """ Gets the extent in world coords as a 4 element tuple. To be implmented by derived class """ raise NotImplementedError("getWorldExtent needs to be overridden")
[docs] def setWorldExtent(self, extent): """ Sets the extent in world coords as a 4 element tuple. To be implmented by derived class """ raise NotImplementedError("setWorldExtent needs to be overridden")
[docs] def getFullWorldExtent(self): """ Gets the full extent of the dataset in world coords as a 4 element tuple. To be implmented by derived class """ raise NotImplementedError("getFullWorldExtent needs to be overridden")
[docs] def getWorldCenter(self): """ Gets the center of the extent in world coords """ (left, top, right, bottom) = self.getWorldExtent() x = left + (right - left) / 2.0 y = bottom + (top - bottom) / 2.0 return x, y
[docs] def setWorldCenter(self, wldX, wldY): """ Sets the center of the extent in world coords """ currentX, currentY = self.getWorldCenter() diffX = wldX - currentX diffY = wldY - currentY (left, top, right, bottom) = self.getWorldExtent() extents = (left + diffX, top + diffY, right + diffX, bottom + diffY) self.setWorldExtent(extents)
[docs]class VectorCoordManager(CoordManager): """ Manages coords for a vector layer """ def __init__(self): CoordManager.__init__(self) self.extent = None self.fullExtent = None self.metersperpix = None
[docs] def recalc(self): """ Recalculate self.metersperpix Called when extents or display size changes """ if self.extent is not None and self.dspWidth is not None: metersaccross = self.extent[2] - self.extent[0] self.metersperpix = metersaccross / self.dspWidth
[docs] def setDisplaySize(self, width, height): """ derived implementation - calls recalculates extent """ CoordManager.setDisplaySize(self, width, height) if self.extent is not None: (left, top, right, bottom) = self.extent self.extent = (left, top, left + self.metersperpix * width, top - self.metersperpix * height)
[docs] def getWorldExtent(self): "Get extent in world coords" return self.extent
[docs] def setWorldExtent(self, extent): "Set extent in world coords" self.extent = extent self.recalc()
[docs] def getFullWorldExtent(self): "full extent of dataset" return self.fullExtent
[docs] def setFullWorldExtent(self, extent): "sets the full extent of dataset" self.fullExtent = extent
[docs] def world2display(self, wldX, wldY): """ convert world coords to display coords returns None if outside May have C implementation also. Not sure yet. """ display = None if self.extent is not None: xoff = wldX - self.extent[0] yoff = self.extent[1] - wldY if xoff >= 0 and yoff >= 0: dspX = xoff / self.metersperpix dspY = yoff / self.metersperpix if dspX < self.dspWidth and dspY < self.dspHeight: display = (dspX, dspY) return display
[docs] def display2world(self, dspX, dspY): """ convert display coords to world coords """ world = None if self.extent is not None: xoff = dspX * self.metersperpix yoff = dspY * self.metersperpix world = (self.extent[0] + xoff, self.extent[1] - yoff) return world
[docs] def recalcBottomRight(self): "I don't think we need to do anything here" pass
[docs]class RasterCoordManager(CoordManager): """ Manage the relationship between the coordinate system used for display and the other coordinate systems used in the raster file. An instance of this class represents the current relationship for a single raster, for a single display. Methods are provided for updating the transformation(s), and for transforming between the different coordinate systems. Coordinate systems involved are: display coords - this notionally corresponds to the screen pixels, although technically it is the units which Qt exposes as its viewport coordinates pixel coords - this is the pixel row/column coordinates in the raster file, using the GDAL conventions world coords - this is the projected coordinate system of the raster file, using the GDAL coordinates. In all cases, coordinate pairs are given with the horizontal coordinate first, i.e. (x, y), even when referring to row/col pairs. Thus, a row/col pair will be given as (col, row). """ def __init__(self): CoordManager.__init__(self) # The raster row/col which is to live in the top-left # corner of the display self.pixTop = None self.pixLeft = None # And the bottom-right self.pixBottom = None self.pixRight = None # Ratio of raster pixels to display pixels. # This defines the zoom level. self.imgPixPerWinPix = None # GDAL geotransform array, which defines relationship between # pixel and world coords self.geotransform = None # size of the raster self.datasetSizeX = None self.datasetSizeY = None def __str__(self): """ For debugging, so I can see what I am set to """ return ("dw:%s dh:%s pt:%s pl:%s pb:%s pr:%s z:%s gt:%s" % ( self.dspWidth, self.dspHeight, self.pixTop, self.pixLeft, self.pixBottom, self.pixRight, self.imgPixPerWinPix, self.geotransform))
[docs] def setTopLeftPixel(self, leftcol, toprow): """ Set row/col of the top/left pixel to display. Args are pixel row/column numbers """ self.pixTop = toprow self.pixLeft = leftcol
[docs] def setGeoTransformAndSize(self, transform, xsize, ysize): """ Set the GDAL geotransform array and size """ self.geotransform = transform self.datasetSizeX = xsize self.datasetSizeY = ysize
[docs] def calcZoomFactor(self, right, bottom): """ Calculate the zoom factor, given the currently set top/left pixel to display, and the bottom/right pixel desired in display, for the currently set display size. The zoom factor is calculated to come as close as possible to display the given section of the raster, with the current display size, but will correct for any difference in aspect ratio between the display window and the desired region of the raster. This means that the given bottom/right values are not always actually displayed. One of them, either bottom or right, will be maintained, but the other will be adjusted to match the aspect ratio of the display. For this reason, the bottom/right values are not stored on the object, but instead the calculated zoom factor is stored. The whole of the desired region will be fitted into the display window. """ displayAspectRatio = self.dspWidth / self.dspHeight rastWidth = right - self.pixLeft rastHeight = bottom - self.pixTop # workaround for user drawing # horizontal or vertical lines which # causes problems below if rastWidth == 0: rastWidth = 1 if rastHeight == 0: rastHeight = 1 rastAspectRatio = rastWidth / rastHeight if rastAspectRatio < displayAspectRatio: rastWidth = displayAspectRatio * rastHeight right = self.pixLeft + rastWidth elif rastAspectRatio > displayAspectRatio: rastHeight = rastWidth / displayAspectRatio bottom = self.pixTop + rastHeight self.imgPixPerWinPix = (right - self.pixLeft) / self.dspWidth self.pixBottom = bottom self.pixRight = right
[docs] def recalcBottomRight(self): """ Called when the window shape has changed. The pixBottom and pixRight values are recalculated, based on the new window shape and the existing zoom factor """ self.pixRight = self.pixLeft + self.imgPixPerWinPix * self.dspWidth self.pixBottom = self.pixTop + self.imgPixPerWinPix * self.dspHeight
[docs] def setZoomFactor(self, imgPixPerWinPix): """ Set the zoom to the given value of imgPixPerWinPix. Will then recalcBottomRight(). """ self.imgPixPerWinPix = imgPixPerWinPix self.recalcBottomRight()
[docs] def display2pixel(self, x, y): """ Convert from display units to raster row/col. Returns a tuple (col, row), as floats. """ col = self.pixLeft + x * self.imgPixPerWinPix row = self.pixTop + y * self.imgPixPerWinPix return (col, row)
[docs] def pixel2display(self, col, row): """ Convert raster row/col to display units. Returns a tuple of (x, y). These are int values, as that appears to be all Qt will ever deal with. """ x = int((col - self.pixLeft) / self.imgPixPerWinPix) y = int((row - self.pixTop) / self.imgPixPerWinPix) return (x, y)
[docs] def pixel2displayF(self, col, row): """ Convert raster row/col to display units. Returns a tuple of (x, y). This version returns floats """ x = (col - self.pixLeft) / self.imgPixPerWinPix y = (row - self.pixTop) / self.imgPixPerWinPix return (x, y)
[docs] def pixel2world(self, col, row): """ Convert raster row/col to world coordinate system. Returns a tuple of floats (x, y) """ gt = self.geotransform x = gt[0] + col * gt[1] + row * gt[2] y = gt[3] + col * gt[4] + row * gt[5] return (x, y)
[docs] def world2pixel(self, x, y): """ Convert world x,y coordinates to raster row/col. Returns a tuple (col, row), as floats. """ gt = self.geotransform # Classic 2x2 matrix inversion det = gt[1] * gt[5] - gt[2] * gt[4] col = (gt[5] * (x - gt[0]) - gt[2] * (y - gt[3])) / det row = (-gt[4] * (x - gt[0]) + gt[1] * (y - gt[3])) / det return (col, row)
[docs] def display2world(self, dspX, dspY): """ Convert display (x, y) to world coordinates. Returns a tuple of floats (x, y), in the world coordinate system """ (col, row) = self.display2pixel(dspX, dspY) (wldX, wldY) = self.pixel2world(col, row) return (wldX, wldY)
[docs] def world2display(self, wldX, wldY): """ Convert world (x, y) to display coordinates. Returns a tuple of int values (x, y) in display coordinate system """ (col, row) = self.world2pixel(wldX, wldY) (dspX, dspY) = self.pixel2display(col, row) return (dspX, dspY)
[docs] def getWorldExtent(self): """ Get the extent of the displayed area in world coords A 4 element tuple is returned. """ (left, top) = self.display2world(0, 0) (right, bottom) = self.display2world(self.dspWidth, self.dspHeight) return (left, top, right, bottom)
[docs] def setWorldExtent(self, extent): """ Sets the world extent as a 4 element tuple """ (leftWorld, topWorld, rightWorld, bottomWorld) = extent (left, top) = self.world2pixel(leftWorld, topWorld) self.setTopLeftPixel(left, top) (right, bottom) = self.world2pixel(rightWorld, bottomWorld) self.calcZoomFactor(right, bottom)
[docs] def getFullWorldExtent(self): """ Gets the full extent of the dataset in world coords as a 4 element tuple. """ (left, top) = self.pixel2world(0, 0) (right, bottom) = self.pixel2world(self.datasetSizeX-1, self.datasetSizeY-1) return (left, top, right, bottom)