proj

I have a openstreetmap image. left upper corner is lat: 80 lon: -169 and right lower corner is lat:-56 lon: 190 The image has the width: 2618 and the height: 1513 How do I convert given lat/lon-coordinates to the pictures x/y-coordinates? I tried it by myself, but couldn't figure out how to convert the coordinates. I hope someone can answer my question. The question, as asked, is almost nonsensical. The area requested covers almost the entire globe, from 80N 169W to 56S 190E. That leaves out the band pole to pole from 170W to 169W, a disk around the north pole north of 80N, and a larger disk

I ran the Python code below that is an example of "Plotting Maps: Visualizing Haiti Earthquake Crisis Data" on a book, Python for Data Analysis . Page 242-246 The code is supposed to create a plot map of Haiti but I got an error as below: Traceback (most recent call last): File "Haiti.py", line 74, in <module> x, y = m(cat_data.LONGITUDE, cat_data.LATITUDE) File "/usr/local/lib/python2.7/site-packages/mpl_toolkits/basemap/__init__.py", line 1148, in __call__ xout,yout = self.projtran(x,y,inverse=inverse) File "/usr/local/lib/python2.7/site-packages/mpl_toolkits/basemap/proj.py", line 286, in _

I have Googled myself to death.. I am attempting to write 2 php functions that will return X and Y from Lat and Long, in both Mercator and flat non-projected (grid) maps. Problem being every calculation I have ran across assumes your map(s) has the same lat and lon at the corners, and then the result is in meters. ugh Here is what I have .. maps of different sizes, different lat's, long's at the 4 corners. I downloaded the Proj4 php port, but with zero documentation and more code then I need, I was over whelmed ... HELP !! This presumes your lat/long coords are decimal values, and don't change

I've been using using ggplot2 to plot climatic gridded data for years. These are usually projected NetCDF files. Cells are square in model coordinates, but depending on which projection the model uses, might not be so in the real world. My usual approach is to remap the data first on a suitable regular grid, and then plot. This introduces a small modification to the data, usually this is acceptable. However, I've decided this is not good enough anymore: I want to plot the projected data directly, without remapping, as other programs (e.g. ncl ) can, if I am not mistaken, do, without touching