Understand shapefiles in different projections (notably OSGB)
[potlatch2.git] / com / gradoservice / proj4as / proj / ProjEqui.as
1 /*******************************************************************************\r
2 NAME                             EQUIRECTANGULAR \r
3 \r
4 PURPOSE:        Transforms input longitude and latitude to Easting and\r
5                 Northing for the Equirectangular projection.  The\r
6                 longitude and latitude must be in radians.  The Easting\r
7                 and Northing values will be returned in meters.\r
8 \r
9 PROGRAMMER              DATE\r
10 ----------              ----\r
11 T. Mittan               Mar, 1993\r
12 \r
13 ALGORITHM REFERENCES\r
14 \r
15 1.  Snyder, John P., "Map Projections--A Working Manual", U.S. Geological\r
16     Survey Professional Paper 1395 (Supersedes USGS Bulletin 1532), United\r
17     State Government Printing Office, Washington D.C., 1987.\r
18 \r
19 2.  Snyder, John P. and Voxland, Philip M., "An Album of Map Projections",\r
20     U.S. Geological Survey Professional Paper 1453 , United State Government\r
21     Printing Office, Washington D.C., 1989.\r
22 *******************************************************************************/\r
23 package com.gradoservice.proj4as.proj\r
24 {\r
25         import com.gradoservice.proj4as.Datum;\r
26         import com.gradoservice.proj4as.ProjConstants;\r
27         import com.gradoservice.proj4as.ProjPoint;\r
28         \r
29         public class ProjEqui extends AbstractProjProjection\r
30         {\r
31                 public function ProjEqui(data:ProjParams)\r
32                 {\r
33                         super(data);\r
34                 }\r
35                 \r
36 \r
37 \r
38   override public function init():void\r
39   {\r
40     if(!this.x0) this.x0=0;\r
41     if(!this.y0) this.y0=0;\r
42     if(!this.lat0) this.lat0=0;\r
43     if(!this.long0) this.long0=0;\r
44     ///this.t2;\r
45   }\r
46 \r
47 \r
48 \r
49 /* Equirectangular forward equations--mapping lat,long to x,y\r
50   ---------------------------------------------------------*/\r
51   override public function forward(p:ProjPoint):ProjPoint \r
52   {\r
53 \r
54     var lon:Number=p.x;                         \r
55     var lat:Number=p.y;                 \r
56 \r
57     var dlon:Number = ProjConstants.adjust_lon(lon - this.long0);\r
58     var x:Number = this.x0 +this. a * dlon *Math.cos(this.lat0);\r
59     var y:Number = this.y0 + this.a * lat;\r
60 \r
61     this.t1=x;\r
62     this.t2=Math.cos(this.lat0);\r
63     p.x=x;\r
64     p.y=y;\r
65     return p;\r
66   }  //equiFwd()\r
67 \r
68 \r
69 \r
70 /* Equirectangular inverse equations--mapping x,y to lat/long\r
71   ---------------------------------------------------------*/\r
72    override public function inverse(p:ProjPoint):ProjPoint\r
73    {\r
74 \r
75     p.x -= this.x0;\r
76     p.y -= this.y0;\r
77     var lat:Number = p.y /this. a;\r
78 \r
79     if ( Math.abs(lat) > ProjConstants.HALF_PI) {\r
80         trace("equi:Inv:DataError");\r
81     }\r
82     var lon:Number = ProjConstants.adjust_lon(this.long0 + p.x / (this.a * Math.cos(this.lat0)));\r
83     p.x=lon;\r
84     p.y=lat;\r
85     \r
86     return p;\r
87   }//equiInv()\r
88                                 \r
89         }\r
90 }