python - Plot .tif GDAL raster using matplotlib Basemap -


i have bit of code here using plotting geotiff image (landsat):

from mpl_toolkits.basemap import basemap osgeo import osr, gdal import matplotlib.pyplot plt import numpy np  # read data , metadata ds = gdal.open(filename)  data = ds.readasarray() gt = ds.getgeotransform() proj = ds.getprojection()  xres = gt[1] yres = gt[5]  # edge coordinates , add half resolution  # go center coordinates xmin = gt[0] + xres * 0.5 xmax = gt[0] + (xres * ds.rasterxsize) - xres * 0.5 ymin = gt[3] + (yres * ds.rasterysize) + yres * 0.5 ymax = gt[3] - yres * 0.5  ds = none  # create grid of xy coordinates in original projection xy_source = np.mgrid[xmin:xmax+xres:xres, ymax+yres:ymin:yres] # create figure , basemap object m = basemap(projection='robin', lon_0=0, resolution='c')  # create projection objects convertion inproj = osr.spatialreference() inproj.importfromwkt(proj)  # target projection basemap object outproj = osr.spatialreference() outproj.importfromproj4(m.proj4string)  size = xy_source[0,:,:].size ct = osr.coordinatetransformation(inproj, outproj) xy_target = np.array(ct.transformpoints(xy_source.reshape(2, size).t))

but, fails @ ct.transformpoints(xy_source.reshape(2, size).t)) , i'm not sure why. error gives me:

typeerror: in method 'coordinatetransformation_transformpoints', argument 1 of type 'osrcoordinatetransformationshadow *'

which not understand. osr guru's out there?

thanks reading.

edit 1 projection of .tiff

>>> print proj out[20]: 'projcs["wgs 84 / utm zone 34n",geogcs["wgs 84",datum["wgs_1984", spheroid["wgs 84",6378137,298.257223563,authority["epsg","7030"]], authority["epsg","6326"]],primem["greenwich",0], unit["degree",0.0174532925199433],authority["epsg","4326"]], projection["transverse_mercator"],parameter["latitude_of_origin",0], parameter["central_meridian",21],parameter["scale_factor",0.9996], parameter["false_easting",500000],parameter["false_northing",0], unit["metre",1,authority["epsg","9001"]],authority["epsg","32634"]]'

also, 

>>> m.proj4string out[43]: '+lon_0=0.0 +y_0=8615499.05007 +r=6370997.0 +proj=robin +x_0=16986796.165 +units=m '

the solution install proj4 package. otherwise, input projection not understood osr. 

conda install -c https://conda.binstar.org/jgomezdans proj4

Comments

Post a Comment

Popular posts from this blog

javascript - Jquery show_hide, what to add in order to make the page scroll to the bottom of the hidden field once button is clicked -

i have button on website when clicked, reveals map. button can seen when page first loads up, when clicked, section holding map slides open , goes out of screen since page loads @ top (normal). want page scroll down bottom of div select when button pushed, if map out of screen. if page scrolled down, , clicks show map button, don't want page jump bottom of map. i tried figure out myself, no luck. seasoned jquery coders able give me guidance? many thanks! map = mapdiv.gmap3("get"); infobox = new infobox({ pixeloffset: new google.maps.size(-220, -310), closeboxurl: '', enableeventpropagation: true }); mapdiv.delegate('.infobox .close','click',function () { infobox.close(); }); jquery(".slidingdiv").hide(); jquery(".show_hide").show(); jquery('.show_hide').click(function(){ lastcenter=map.getcenter(); jquery(".slidingdiv&quo
Read more

javascript - Highcharts multi-color line -

Image
i looking create highchart multi-colored line. if goes down under value, become red, , when comes on value, become green . here picture example: or maybe highcharts have other charts models differentiate high or low level color? you can achieve coloring using lineargradient coloring: color: { lineargradient: { x1: 0, x2: 0, y1: 0, y2: 1 }, stops: [ [0, 'red'], [0.25, 'yellow'], [0.50, 'green'], [0.75, 'yellow'], [1, 'red'] ] } here 0 in stops relates minimum value of series , 1 maximum value, , in between percentages between values. if find percentages according values prior creating chart use correct stops values make relate specific values in data (at least approximately). see this jsfiddle example .
Read more

javascript - Enter key does not work in search box -

tried search guys sorry. input in search box doesnt work (or search) when hit enter when user clicks search button. doing wrong? here code. thanks. <script type="text/javascript"> $(document).ready(function() { $("#submit").click(function() { var url = "http://asla.org/awardssearch.html"; url += "?s=" + $('#googlecse').val(); window.location = url; }); $('#googlecse').keydown(function(event) { if (event.keycode == '13') { $("#submit").click(); } else { //do nothing :) } }); }); </script> <form class="navbar-form navbar-right" role="search"> <div class="search"> <input id="googlecse" type="text" onblur="if(this.value=='')this.value=this.defaultvalue;" onfocus="if(this.value==this.defaultvalue)t
Read more