Jeff, great, thank you! -- Anja Rösel
Jeff Whitaker wrote: > Anja Roesel wrote: >> Dear list, >> >> I am a matplotlib newbie and I have some simple problems with the >> coordinate reprojection. >> >> I have a landsat scene in UTM Projection and I would like to plot it in >> a polarstereograhic projection (it is in the Arctic) >> >> I tried it like this: >> m=Basemap(resolution='i',projection='npstere',lon_0=-45,boundinglat=70) >> XP,YP = m(X_utm,Y_utm) >> >> but the outcome is nonsense (an array with 1.0's) >> I am wondering if there is an option to specify the input parameters >> like utm zone etc. >> >> Thanks for any help and ideas! >> Anja >> > > Anja: You need to pass the latitude and longitude values (in degrees) > to the Basemap instance when converting to projection coordinates. So, > to convert from UTM coordinates to polar stereographic coordinates you > will need to do something like this: > > map1 = Basemap(<parameters for transverse mercator projection>) > > lons, lats = map1(x, y, inverse=True) # x and y are projection > coordinates on original UTM grid. > > # lons and lats are now latitudes and longitudes of UTM grid (in degrees) > > map2 = > > Basemap(resolution='i',projection='npstere',lon_0=-45,boundinglat=70) > > x, y = map2(lons, lats) > > # x,y are now polar stereo coordinates of UTM grid. > > Or, if you already have the latitudes and logintudes of the original UTM > grid you can skip the first two lines and just pass those to the > stereographic Basemap instance. > > -Jeff > ------------------------------------------------------------------------------ Let Crystal Reports handle the reporting - Free Crystal Reports 2008 30-Day trial. Simplify your report design, integration and deployment - and focus on what you do best, core application coding. Discover what's new with Crystal Reports now. http://p.sf.net/sfu/bobj-july _______________________________________________ Matplotlib-users mailing list Matplotlib-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/matplotlib-users
