title: Display a georeferenced DEM surface in 3D matplotlib
body: <p>I want to use a DEM file to generate a simulated terrain surface using matplotlib.  But I do not know how to georeference the raster coordinates to a given CRS.  Nor do I know how to express the georeferenced raster in a format suitable for use in a 3D matplotlib plot, for example as a numpy array.  </p>
  
  <p>Here is my python code so far:</p>
  
  <pre><code>import osgeo.gdal
  
  dataset = osgeo.gdal.Open("MergedDEM")
  
  gt = dataset.GetGeoTransform()
  </code></pre>