title: how to plot geotiff data in specific area (lat/lon) with python
body: <p>I have a geotiff raster data sets with elevation data init and i want to plot it in specific area, such as 60°E - 70° E ，70°S - 80°E. </p>
  
  <p>I have a bit of code from <a href="https://stackoverflow.com/questions/20488765/plot-gdal-raster-using-matplotlib-basemap/20499331#20499331">here</a>,but the <code>pcolormesh</code> seem couldn't plot my geotif.it's all red. <a href="https://i.stack.imgur.com/mHuiS.jpg" rel="nofollow noreferrer">picture</a>.  The picture is shown by <code>imshow</code> as <a href="https://i.stack.imgur.com/X0wm0.jpg" rel="nofollow noreferrer">really picture</a></p>
  
  <p>When I try to make a plot with this code below:</p>
  
  <pre><code>path = "F:\\Mosaic_h1112v28_ps.tif"
  dataset = gdal.Open(path)
  data = dataset.ReadAsArray()
  x0, dx, dxdy, y0, dydx, dy = dataset.GetGeoTransform()
  nrows, ncols = data.shape
  londata = np.linspace(x0, x0+dx*ncols)
  latdata = np.linspace(y0, y0+dy*nrows)
  lons, lats = np.meshgrid(lonarray, latarray)  
  fig = plt.figure(figsize=(8, 8))
  m = Basemap(projection='lcc', lon_0=67.5, lat_0=-68.5, height=950000,
              width=580000, resolution='h') 
  m.drawcoastlines() 
  x, y = m(lons, lats) 
  </code></pre>
  
  <p>Then i dont know how to continue it . I just want to use <code>imshow</code>, but the <code>imshow</code> dont specify area(lat/lon).</p>
  
  <p>I will really appreciate your help.</p>