X
PrevPrev Go to previous topic
NextNext Go to next topic
Last Post 21 Jul 2013 07:08 AM by  anon
map LIMIT in function graphics routine
 2 Replies
Sort:
You are not authorized to post a reply.
Author Messages

anon



New Member


Posts:
New Member


--
21 Jul 2013 07:08 AM
    Hi, I am trying to use IDL's (8.2.3) function graphics for displaying projected imagery and vector data. I am using Albers EA projection with standard parallels between 25 and 0 degrees latitude. The geographic area I am trying to display ranges from -10d lat to 30d lat. However, IDL displays map grids and continent outlines below 0 degree lat only if I lower the second standard parallel. My guess is that IDL uses the standard_par2 as the southern map boundary, which would be an incorrect assumption. Below is an example. Any help is much appreciated. mp = map('albers equal area', semimajor_axis=6378137.d, $ semiminor_axis=6356752.3d, standard_par1=25d, standard_par2=0d, $ center_longitude=110d, center_latitude=-15d, $ false_easting=0d, false_northing=0d, limit=[-10,90,30,140]) m1 = mapcontinents( /countries)

    Deleted User



    New Member


    Posts:
    New Member


    --
    23 Jul 2013 11:19 AM
    This comes from ArcGIS Help 10.1: One method to calculate the standard parallels is by determining the range in latitude in degrees north to south and dividing this range by six. The "one-sixth rule" places the first standard parallel at one-sixth the range above the southern boundary and the second standard parallel minus one-sixth the range below the northern limit. There are other possible approaches. Try : standard_par1=23d, $ standard_par2=-3d, $

    Deleted User



    New Member


    Posts:
    New Member


    --
    23 Jul 2013 04:03 PM
    Hi, thank you for your answer. Surely, the fix to the issue can't be to change the map projection. I am aware of the 1/6's rule for the standard parallels. I didn't think the order mattered, but I checked and it seems that IDL doesn't care. The standard parallels should be placed slightly within the boundary of the area of interest to balance distortions towards the center of the projection and towards the poles. Thus, the map projection is valid beyond the standard parallels but with increasing distortions towards the poles. What seems to happen in IDL (and not other software, e.g. ArcGIS) is that if the second standard parallel is 0 or -0, then the map display seems limited to the northern hemisphere. If the SP is set to a tiny negative value then the map displays as expected. The following code works, so I'm probably going to use this trick for now. But there must be a better solution... mp = map('albers equal area', semimajor_axis=6378137.d, $ semiminor_axis=6356752.3d, standard_par1=25d, standard_par2=-10d^(-6), $ center_longitude=110d, center_latitude=-15d, $ false_easting=0d, false_northing=0d, limit=[-10,90,30,140]) m1 = mapcontinents( /countries)
    You are not authorized to post a reply.