/[thuban]/trunk/thuban/Thuban/Model/layer.py

Diff of /trunk/thuban/Thuban/Model/layer.py

Parent Directory | Revision Log | View Patch Patch

-revision 1983 by bh,
Thu Nov 27 15:57:23 2003 UTC
+revision 2644 by bh,
Tue Jul  5 16:30:51 2005 UTC
 Line 1
- # Copyright (c) 2001, 2002, 2003 by Intevation GmbH
+ # Copyright (c) 2001, 2002, 2003, 2004, 2005 by Intevation GmbH
  # Authors:
  # Bernhard Herzog <[email protected]>
  # Jonathan Coles <[email protected]>
+ # Silke Reimer <[email protected]>
  #
  # This program is free software under the GPL (>=v2)
  # Read the file COPYING coming with Thuban for details.
-Line 24 
 from data import SHAPETYPE_POLYGON, SHAP
+Line 25 
 from data import SHAPETYPE_POLYGON, SHAP
  import resource
+ from color import Color
  shapetype_names = {SHAPETYPE_POINT: "Point",
                     SHAPETYPE_ARC: "Arc",
-Line 54 
 class BaseLayer(TitledObject, Modifiable
+Line 56 
 class BaseLayer(TitledObject, Modifiable
          self.issue(LAYER_VISIBILITY_CHANGED, self)
      def HasClassification(self):
-         """Determine if this layer support classifications."""
+         """Determine if this layer supports classifications."""
          return False
      def HasShapes(self):
-Line 66 
 class BaseLayer(TitledObject, Modifiable
+Line 68 
 class BaseLayer(TitledObject, Modifiable
          return self.projection
      def SetProjection(self, projection):
-         """Set the layer's projection"""
+         """Set the layer's projection."""
          self.projection = projection
          self.changed(LAYER_PROJECTION_CHANGED, self)
+     def Type(self):
+         return "Unknown"
  class Layer(BaseLayer):
      """Represent the information of one geodata file (currently a shapefile)
-Line 170 
 class Layer(BaseLayer):
+Line 175 
 class Layer(BaseLayer):
              bbox = self.projection.InverseBBox(bbox)
          return bbox
+     def Type(self):
+         return self.ShapeType();
      def ShapesBoundingBox(self, shapes):
          """Return a bounding box in lat/long coordinates for the given
          list of shape ids.
-Line 226 
 class Layer(BaseLayer):
+Line 234 
 class Layer(BaseLayer):
          form (minx, miny, maxx, maxy) in unprojected coordinates.
          """
          if self.projection is not None:
-             bbox = self.projection.ForwardBBox(bbox)
+             # Ensure that region lies within the layer's bounding box
+             # Otherwise projection of the region would lead to incorrect
+             # values.
+             clipbbox = self.__mangle_bounding_box(bbox)
+             bbox = self.projection.ForwardBBox(clipbbox)
          return self.store.ShapesInRegion(bbox)
      def GetClassificationColumn(self):
-Line 289 
 class Layer(BaseLayer):
+Line 301 
 class Layer(BaseLayer):
          bbox = self.LatLongBoundingBox()
          if bbox is not None:
-             items.append(_("Extent (lat-lon): (%g, %g, %g, %g)") % bbox)
+             items.append(_("Extent (lat-lon): (%g, %g, %g, %g)") % tuple(bbox))
          else:
              items.append(_("Extent (lat-lon):"))
          items.append(_("Shapetype: %s") % shapetype_names[self.ShapeType()])
-Line 302 
 class Layer(BaseLayer):
+Line 314 
 class Layer(BaseLayer):
          return (_("Layer '%s'") % self.Title(), items)
+     def __mangle_bounding_box(self, bbox):
+         # FIXME: This method doesn't make much sense.
+         # See RT #2845 which effectively says:
+         #
+         # If this method, which was originally called ClipBoundingBox,
+         # is supposed to do clipping it shouldn't return the parameter
+         # unchanged when it lies completely outside of the bounding box.
+         # It would be better to return None and return an empty list in
+         # ShapesInRegion (the only caller) in that case.
+         #
+         # This method was introduced to fix a bug that IIRC had
+         # something todo with projections and bounding boxes containing
+         # NaN or INF when the parameter to ShapesInRegion covered the
+         # entire earth or something similarly large).
+         bminx, bminy, bmaxx, bmaxy = bbox
+         lminx, lminy, lmaxx, lmaxy = self.LatLongBoundingBox()
+         if bminx > lmaxx or bmaxx < lminx:
+             left, right = bminx, bmaxx
+         else:
+             left = max(lminx, bminx)
+             right = min(lmaxx, bmaxx)
+         if bminy > lmaxy or bmaxy < lminy:
+             bottom, top = bminy, bmaxy
+         else:
+             bottom = max(lminy, bminy)
+             top = min(lmaxy, bmaxy)
+         return (left, bottom, right, top)
  if resource.has_gdal_support():
      import gdal
-Line 309 
 if resource.has_gdal_support():
+Line 350 
 if resource.has_gdal_support():
  class RasterLayer(BaseLayer):
-     def __init__(self, title, filename, projection = None, visible = True):
+     MASK_NONE  = 0
+     MASK_BIT   = 1
+     MASK_ALPHA = 2
+     def __init__(self, title, filename, projection = None,
+                  visible = True, opacity = 1, masktype = MASK_BIT):
          """Initialize the Raster Layer.
          title -- title for the layer.
-Line 332 
 class RasterLayer(BaseLayer):
+Line 378 
 class RasterLayer(BaseLayer):
          self.bbox = -1
+         self.mask_type = masktype
+         self.opacity = opacity
+         self.image_info = None
          if resource.has_gdal_support():
              #
              # temporarily open the file so that GDAL can test if it's valid.
-Line 341 
 class RasterLayer(BaseLayer):
+Line 392 
 class RasterLayer(BaseLayer):
              if dataset is None:
                  raise IOError()
+             #
+             # while we have the file, extract some basic information
+             # that we can display later
+             #
+             self.image_info = {}
+             self.image_info["nBands"] = dataset.RasterCount
+             self.image_info["Size"] = (dataset.RasterXSize, dataset.RasterYSize)
+             self.image_info["Driver"] = dataset.GetDriver().ShortName
+             # store some information about the individual bands
+             # [min_value, max_value]
+             a = self.image_info["BandData"] = []
+             for i in range(1, dataset.RasterCount+1):
+                 band = dataset.GetRasterBand(i)
+                 a.append(band.ComputeRasterMinMax())
          self.UnsetModified()
      def BoundingBox(self):
-Line 395 
 class RasterLayer(BaseLayer):
+Line 464 
 class RasterLayer(BaseLayer):
          return bbox
+     def Type(self):
+         return "Image"
      def GetImageFilename(self):
          return self.filename
+     def MaskType(self):
+         """Return True if the mask should be used when rendering the layer."""
+         return self.mask_type
+     def SetMaskType(self, type):
+         """Set the type of mask to use.
+         type can be one of MASK_NONE, MASK_BIT, MASK_ALPHA
+         If the state changes, a LAYER_CHANGED message is sent.
+         """
+         if type not in (self.MASK_NONE, self.MASK_BIT, self.MASK_ALPHA):
+             raise ValueError("type is invalid")
+         if type != self.mask_type:
+             self.mask_type = type
+             self.changed(LAYER_CHANGED, self)
+     def Opacity(self):
+         """Return the level of opacity used in alpha blending.
+         """
+         return self.opacity
+     def SetOpacity(self, op):
+         """Set the level of alpha opacity.
+<= op <= 1.
+         The layer is fully opaque when op = 1.
+         """
+         if not (0 <= op <= 1):
+             raise ValueError("op out of range")
+         if op != self.opacity:
+             self.opacity = op
+             self.changed(LAYER_CHANGED, self)
+     def ImageInfo(self):
+         return self.image_info
      def TreeInfo(self):
          items = []

 Legend:



Removed from v.1983
 


changed lines


 
Added in v.2644
 Legend:



Removed from v.1983
 


changed lines


 
Added in v.2644
-Removed from v.1983
+Added in v.2644

[email protected]	ViewVC Help
Powered by ViewVC 1.1.26