Results 1 to 6 of 6

Thread: cylindrical equal-area projection

Hybrid View

Previous Post Previous Post   Next Post Next Post
  1. #1

    Help cylindrical equal-area projection

    I have a bit of a problem with the cylindrical equal-area projection. I can't figure out how i can project it correctly from longitude , latitude and altitude data to x, y and altitude.

    I need an equal-area projection, because i want to use this for a game map in which i want all the area's to have the same size (that would be fair for the players). But I also want to keep realism, so i need a correct projection of the map.

    Or can i better use an other projection method to achieve this (like for example the mollerweide projection)?

    i hope my question is clear.
    Last edited by haroim; 03-20-2009 at 11:31 AM.

  2. #2
    Guild Artisan su_liam's Avatar
    Join Date
    Aug 2007
    Location
    Port Alberta, Regina(IRL: Eugene, OR)
    Posts
    798

    Post

    Here are some formulas:
    Each parallel has a length of 2*pi*cos(theta), this is the height of the map in y.
    Each meridian has a length of 2*R/cos(theta), this is the width of the map in x.
    The distance of each parallel to the equator is R*sin(phi)/cos(theta).
    Where R is the radius of the generating globe of the chosen area scale, theta is the standard parallel, phi is the latitude of interest.

    For a point on the globe (lat/lon), where lat is the latitude of the point and lon is the longitude of the point(in degrees) and the map is centered at latlon(0,0), west longitudes are negative values as are south latitudes,

    x = 2*pi*cos(theta) * (lon + 180)/360
    y = R*sin(lat)/cos(phi)

    The ratio of width/height of the map is determined from,

    cos(theta) = acos((ratio/pi)^0.5)

    Some standard projections and their standard parallels:
    Lambert - pi/1 ratio, theta=0
    Behrmann - ~2.356:1 ratio, theta=30

    Gall-Peters - pi/2:1 ratio, theta=45

    Hopefully this would be helpful.

    edit: Info from Arthur Robinson, Randall Sale, Joel Morrison, "Elements of Cartography" 4th ed., 1978, John Wiley and Sons: New York.

    The wikipedia describes this well under Gall-Peters projection. Even more generally, but in less detail under Lambert cylindrical equal-area projection.
    Last edited by su_liam; 03-20-2009 at 02:05 PM. Reason: Added wikipedia references, 'cause the wikipedia adds an air of authority

  3. #3

    Default

    thanks for the reaction! i am now looking into it.

  4. #4

    Question

    I have looked into it and i get out the following:

    If i want to make a projection i will need for each (x,y) an altitude. I can find the altitude in a (lon,lat) form. So to get the altitude of (x,y) i will need to know the (lon,lat) of (x,y).

    (x,y)=(lon,lat)

    x = 2*pi*cos(theta) * (lon + 180)/360
    lon = (x*360)/(2*pi*cos(theta))-180

    y = R*sin(lat)/cos(phi)
    lat=invsin(cos(phi)*y/R)

    so (x,y)=(lon,lat)=( (x*360)/(2*pi*cos(theta))-180, invsin(cos(phi)*y/R)) )

    so i should be able to calculate the altitude of (x,y) by taking the altitude of ( (x*360)/(2*pi*cos(theta))-180, invsin(cos(phi)*y/R)) ).

    is this correct? And what is phi? Is it the standard latitude, and what number should i enter?
    Last edited by haroim; 03-21-2009 at 04:02 AM. Reason: invsin(cos(phi*y/R)) must be invsin(cos(phi)*y/R))

  5. #5
    Guild Artisan su_liam's Avatar
    Join Date
    Aug 2007
    Location
    Port Alberta, Regina(IRL: Eugene, OR)
    Posts
    798

    Default

    It looks correct to me... I'll have to check your algebra when I have more time and my brain is working.

    Yes, phi is the standard latitude. It's a constant for a given map. For Lambert cylindrical equal-area it's 0, the equator. For Behrmann, it's 30degrees. For Gall-Peters it's 45 degrees.

  6. #6

    Praise

    Thanks! So the phi is the same as the theta. Should i take the width of my map as R, or should i take the width, calculate it into a sphere, and calculate the radius of the sphere?
    Last edited by haroim; 03-21-2009 at 04:37 AM.

Tags for this Thread

Posting Permissions

  • You may not post new threads
  • You may not post replies
  • You may not post attachments
  • You may not edit your posts
  •