[Proj] pj_transform function

Hugo Tyson Hugo.Tyson at cpslocates.com
Fri Oct 6 05:15:20 EDT 2006


Sandeep Kumar Jakkaraju [sandeepkumar.jakkaraju at gmail.com] wrote:
> So what should be the  projection definition then .. to make
> it general..???
> 
> I am trying to convert the geographical co-ordinates to
> screen co-ordinates ..
> so i require to convert the lat/lon to some projection and
> the convert that X,Y to pixels...
> hope this is correct ...

Ah.  So you just want like a map of the whole world, as 
a rectangle, sort of thing?  That's not a "projection"
in the proj4 sense, as I understand it.

How about
	double x = (lon + 180.0)/360.0 * pixelswidth;
	double y = (lat +  90.0)/180.0 * pixelsheight;
of course that'll be ugly at far north and south,
but at least it maps the whole space.

> or is there some function in proj4/or anywhere else which
> does this ....??

I don't think this is the sort of thing that proj is for.

Alternatively if you're mapping a small areas of the world
into one page or bitmap, you need to pick a proj4 projection
that works for that small area.

	- Huge

PS I was mistaken in what I wrote before, mercator maps to a
half-cylinder wrapped around the globe.  Your example touches
at long 0.0, so it's still true that your long -104 is "around
the back".


> On 10/6/06, Hugo Tyson <Hugo.Tyson at cpslocates.com> wrote:
> 
> 	proj-bounces at lists.maptools.org <> wrote:
> 	> I want to transform ... this lat/lon 39.739/-104.985...
> 	> to world mercator project ..with the defintition 	>
> 	> World Mercator <54004> +proj=merc +latts=0 +lon0=0
> 	> +k=1.000000 +x0=0 +y0=0 +ellps=WGS84 +datum=WGS84 +units=m no_defs
> 	<> >
> 	> but the method gives an "error"....
> 	>
> 	> this method works fine if i give lat/lon say 42.4182,-74.2643
> 
> 	I've not heard of "world mercator", but looking at the parameters
> 	it's projecting into a plane tangent to the sphere at
> lat 0, lon 0.
> 	So it only works for longitudes from -90 to +90 (though
> these map
> 	to infinities in the projected plane).
> 
> 	Your example of -104 longitude is around the back of the earth,
> 	on the "far side" from the projection plane.   Therefore a
> 	line from the origin through the point never intersects the
> 	projection plane, so no result can be calculated.  Hence your
> 	error.  Lon -74 is in range, so it's OK.

-- 
Hugo Tyson, Cambridge Positioning Systems Ltd
=************************************************************
********************************************************************************
This email and any attachments are confidential. They have been checked for viruses but it is your 
responsibility to conduct your own virus checks.
Cambridge Positioning Systems Limited, 62-64 Hills Road, Cambridge, CB2 1LA, is registered in England 
under company number 2808344.
********************************************************************************




More information about the Proj mailing list