New file |
0,0 → 1,139 |
<?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 |
*/ |
/******************************************************************************* |
NAME SINUSOIDAL |
|
PURPOSE: Transforms input longitude and latitude to Easting and |
Northing for the Sinusoidal projection. The |
longitude and latitude must be in radians. The Easting |
and Northing values will be returned in meters. |
|
PROGRAMMER DATE |
---------- ---- |
D. Steinwand, EROS May, 1991 |
|
This function was adapted from the Sinusoidal projection code (FORTRAN) in the |
General Cartographic Transformation Package software which is available from |
the U.S. Geological Survey National Mapping Division. |
|
ALGORITHM REFERENCES |
|
1. Snyder, John P., "Map Projections--A Working Manual", U.S. Geological |
Survey Professional Paper 1395 (Supersedes USGS Bulletin 1532), United |
State Government Printing Office, Washington D.C., 1987. |
|
2. "Software Documentation for GCTP General Cartographic Transformation |
Package", U.S. Geological Survey National Mapping Division, May 1982. |
* ***************************************************************************** */ |
|
class Proj4phpProjSinu { |
|
/* Initialize the Sinusoidal projection |
------------------------------------ */ |
public function init() { |
/* Place parameters in static storage for common use |
------------------------------------------------- */ |
#$this->R = 6370997.0; //Radius of earth |
|
if( !$this->sphere ) { |
$this->en = Proj4php::$common->pj_enfn( $this->es ); |
} else { |
$this->n = 1.; |
$this->m = 0.; |
$this->es = 0; |
$this->C_y = sqrt( ($this->m + 1.) / $this->n ); |
$this->C_x = $this->C_y / ($this->m + 1.); |
} |
} |
|
/* Sinusoidal forward equations--mapping lat,long to x,y |
----------------------------------------------------- */ |
public function forward( $p ) { |
|
#$x,y,delta_lon; |
$lon = $p->x; |
$lat = $p->y; |
|
/* Forward equations |
----------------- */ |
$lon = Proj4php::$common->adjust_lon( $lon - $this->long0 ); |
|
if( isset($this->sphere) ) { |
if( !$this->m ) { |
$lat = $this->n != 1. ? asin( $this->n * sin( $lat ) ) : $lat; |
} else { |
$k = $this->n * sin( $lat ); |
for( $i = Proj4php::$common->MAX_ITER; $i; --$i ) { |
$V = ($this->m * $lat + sin( $lat ) - $k) / ($this->m + cos( $lat )); |
$lat -= $V; |
if( abs( $V ) < Proj4php::$common->EPSLN ) |
break; |
} |
} |
$x = $this->a * $this->C_x * $lon * ($this->m + cos( $lat )); |
$y = $this->a * $this->C_y * $lat; |
} else { |
|
$s = sin( $lat ); |
$c = cos( $lat ); |
$y = $this->a * Proj4php::$common->pj_mlfn( $lat, $s, $c, $this->en ); |
$x = $this->a * $lon * $c / sqrt( 1. - $this->es * $s * $s ); |
} |
|
$p->x = $x; |
$p->y = $y; |
|
return $p; |
} |
|
/** |
* |
* @param type $p |
* @return type |
*/ |
public function inverse( $p ) { |
#$lat; |
#$temp; |
#$lon; |
|
/* Inverse equations |
----------------- */ |
$p->x -= $this->x0; |
$p->y -= $this->y0; |
$lat = $p->y / $this->a; |
|
if( isset($this->sphere) ) { |
|
$p->y /= $this->C_y; |
$lat = $this->m ? asin( ($this->m * $p->y + sin( $p->y )) / $this->n ) : ( $this->n != 1. ? asin( sin( $p->y ) / $this->n ) : $p->y ); |
$lon = $p->x / ($this->C_x * ($this->m + cos( $p->y ))); |
} |
else { |
$lat = Proj4php::$common->pj_inv_mlfn( $p->y / $this->a, $this->es, $this->en ); |
$s = abs( $lat ); |
|
if( $s < Proj4php::$common->HALF_PI ) { |
$s = sin( $lat ); |
$temp = $this->long0 + $p->x * sqrt( 1. - $this->es * $s * $s ) / ($this->a * cos( $lat )); |
//temp = $this->long0 + $p->x / ($this->a * cos($lat)); |
$lon = Proj4php::$common->adjust_lon( $temp ); |
} else if( ($s - Proj4php::$common->EPSLN) < Proj4php::$common->HALF_PI ) { |
$lon = $this->long0; |
} |
} |
|
$p->x = $lon; |
$p->y = $lat; |
|
return $p; |
} |
|
} |
|
Proj4php::$proj['sinu'] = new Proj4phpProjSinu(); |