| New file |
|---|
| 0,0 → 1,58 |
|---|
| <?php |
| /** |
| * Author : Julien Moquet |
| * |
| * Inspired by Proj4php from Mike Adair madairATdmsolutions.ca |
| * and Richard Greenwood rich@greenwoodma$p->com |
| * License: LGPL as per: http://www.gnu.org/copyleft/lesser.html |
| */ |
| /* similar to equi.js FIXME proj4 uses eqc */ |
| class Proj4phpProjEqc { |
| |
| public function init() { |
| |
| if( !$this->x0 ) |
| $this->x0 = 0; |
| if( !$this->y0 ) |
| $this->y0 = 0; |
| if( !$this->lat0 ) |
| $this->lat0 = 0; |
| if( !$this->long0 ) |
| $this->long0 = 0; |
| if( !$this->lat_ts ) |
| $this->lat_ts = 0; |
| if( !$this->title ) |
| $this->title = "Equidistant Cylindrical (Plate Carre)"; |
| |
| $this->rc = cos( $this->lat_ts ); |
| } |
| |
| // forward equations--mapping lat,long to x,y |
| // ----------------------------------------------------------------- |
| public function forward( $p ) { |
| |
| $lon = $p->x; |
| $lat = $p->y; |
| |
| $dlon = Proj4php::$common->adjust_lon( $lon - $this->long0 ); |
| $dlat = Proj4php::$common . adjust_lat( $lat - $this->lat0 ); |
| $p->x = $this->x0 + ($this->a * $dlon * $this->rc); |
| $p->y = $this->y0 + ($this->a * $dlat ); |
| return $p; |
| } |
| |
| // inverse equations--mapping x,y to lat/long |
| // ----------------------------------------------------------------- |
| public function inverse( $p ) { |
| |
| $x = $p->x; |
| $y = $p->y; |
| |
| $p->x = Proj4php::$common->adjust_lon( $this->long0 + (($x - $this->x0) / ($this->a * $this->rc)) ); |
| $p->y = Proj4php::$common->adjust_lat( $this->lat0 + (($y - $this->y0) / ($this->a )) ); |
| return $p; |
| } |
| |
| } |
| |
| Proj4php::$proj['eqc'] = new Proj4phpProjEqc(); |