[Proj] Re: Proj Digest, Vol 51, Issue 2 - axis rotation issue

Neil Robinson halfhaggis at gmail.com
Wed Aug 13 08:22:32 PDT 2008



Gerald I. Evenden wrote:
> On Wednesday 13 August 2008 4:58:50 am you wrote:
> 	...
>   
>> I think it's my fault. I'd had latitudes as S20;  S30; etc. and
>> should've put the negative in.<br>
>>     
>
> My goodness, that make a lot of difference!
>
>   
>> <br>
>> It looks to me that I simply need to multiply the outputs by -1, but
>> I'm loath to do that because I don't understand the underlying
>> mathematics that calculates the projection. <br>
>> Is it really that simple? I could use the "-m" flag and set it to -1,
>> but is that wise?<br>
>> It appears to give me the correct output:<br>
>> <br>
>> ~$ proj +proj=tmerc +ellps=WGS84 +lon_0=27 -m -1 Desktop/proj<br>
>> -55660.46    -0.00<br>
>> -54820.34    1105896.37<br>
>> -52324.06    2212444.34<br>
>> -48243.45    3320218.65<br>
>> -58652.24    2897715.76 <br>
>> <br>
>> Can you advise?<br>
>>     
>
> You have a most strange grid system where the south pole is oriented to the 
> top of the map.  I suppose it makes a difference to someone living in the 
> southern hemisphere but it is a bit odd when compared to most grid 
> orientations.
>   
Don't ask me -- I just live here! I would have been a lot less 
frustrated if the powers that be could've just kept north at the top.
> All you are doing is rotating the x-y (easting-northing) axis by 180 degrees 
> and this can, as you empirically determined, be done my using the '-m -1' 
> option.
>
> As to whether an axis rotation option is libproj4's responsibility or the 
> applications program's (proj4, lproj, etc.) responsibility and thus a "-" 
> option or "+" option respectively is now open for debate.
>   
Indeed. As it stands, qgis, which allows configuration of a custom 
projection using proj4 parameters, doesn't appear to support the -m flag 
or have an equivalent flag or setting I could use.
Since I started this whole investigation to get qgis to work properly, 
my fight is still not over. At least I can measure distances now, which 
is my most important requirement.

I'll see what the qgis people have to say. Thanks for your time

> Any comments from the maptools group?
>
>   
-------------- next part --------------
A non-text attachment was scrubbed...
Name: halfhaggis.vcf
Type: text/x-vcard
Size: 198 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20080813/e3892fef/attachment.vcf>


More information about the Proj mailing list