Calling all PROJ4 genuises

Chris Hackel cvhackel at YAHOO.COM
Mon Aug 29 12:40:27 PDT 2005


Tyler, thank you *very* much. I'll have to try to
parse through that and learn the parameters better.

My cs2cs is from the proj-4.4.9 release. I'm getting
an error for the 2nd init:
cs2cs +init=epsg:4326 +to +init=102718
Using to definition: init=102718 
Rel. 4.4.9, 29 Oct 2004
<cs2cs>: 
projection initialization failure
cause: no colon in init= string
program abnormally terminated

I'm very interested in getting the proper proj
environment in order to user the simplified version.
Do you know where I could find it?

-Chris

--- Tyler Mitchell <tylermitchell at SHAW.CA> wrote:

> Hi Chris,
> After looking up the EPSG code for a similar
> sounding projection, this gives results closer to
> what you are looking for:
> cs2cs +proj=latlong +ellps=GRS80 +datum=NAD83 +to
> +proj=lcc +lat_1=40.66666666666666
> +lat_2=41.03333333333333 +lat_0=40.16666666666666
> +lon_0=-74 +x_0=300000 +y_0=0 +ellps=GRS80
> +datum=NAD83 +to_meter=0.3048006096012192
> 
> Input: -73.994992 40.717690
> Output: 985638.25       200750.38 0.00
> 
> Where did you get cs2cs and your proj environment
> from?  Depending on this, you may have a shorthand
> way of calling this, if you put a couple lines into
> your epsg code file.  Follow up with us to get that
> done if you are interested.  That would result in a
> much shorter cs2cs line:
> 
> cs2cs +init=epsg:4326 +to +init=102718
> 
> A bit easier to read ;)
> Hope that helps!
> 
> Tyler
> 
> From FWTools\proj_lib\esri:
> # NAD 1983 StatePlane New York Long Island FIPS 3104
> Feet
> <102718> +proj=lcc +lat_1=40.66666666666666
> +lat_2=41.03333333333333 +lat_0=40.16666666666666
> +lon_0=-74 +x_0=300000 +y_0=0 +ellps=GRS80
> +datum=NAD83 +to_meter=0.3048006096012192  no_defs
> <>
> 
> ----- Original Message -----
> From: Chris Hackel <cvhackel at YAHOO.COM>
> Date: Monday, August 29, 2005 12:11 pm
> Subject: [UMN_MAPSERVER-USERS] Calling all PROJ4
> genuises
> 
> > I'm having trouble with coordinate conversion
> (latlong
> > to Lambert/lcc).
> > 
> > I'm using cs2cs in the following way, drawing data
> > from orginfo (below):
> > cs2cs +proj=latlong +ellps=GRS80 +datum=NAD83 +to
> > +proj=lcc +datum=NAD83 +units=ft +lon_0=-74.0
> > +lat_0=40.16666 +lat_1=40.1666666 +lat_2=41.033333
> > +x_0=984250.0 +y_0=0
> > -73.994992 40.717690
> > 3230554.89      200744.22 0.00
> > 
> > These numbers are definitely not what I'm looking
> for.
> > (Extent of data: (910057.005249, 120827.652033) -
> > (1067494.180516, 274629.837850))
> > 
> > Please help! Thank you!
> > -Chris
> > 
> > 
> > 
> >
>
PROJCS["NAD_1983_StatePlane_New_York_Long_Island_FIPS_3104_Feet",
> >    GEOGCS["GCS_North_American_1983",
> >        DATUM["North_American_Datum_1983",
> >           
> > SPHEROID["GRS_1980",6378137.0,298.257222101]],
> >        PRIMEM["Greenwich",0.0],
> >        UNIT["Degree",0.0174532925199433]],
> >    PROJECTION["Lambert_Conformal_Conic_2SP"],
> >    PARAMETER["False_Easting",984250.0],
> >    PARAMETER["False_Northing",0.0],
> >    PARAMETER["Central_Meridian",-74.0],
> >   
> >
> PARAMETER["Standard_Parallel_1",40.66666666666666],
> >   
> >
> PARAMETER["Standard_Parallel_2",41.03333333333333],
> >   
> PARAMETER["Latitude_Of_Origin",40.16666666666666],
> >    UNIT["Foot_US",0.3048006096012192]]
> > 
> 



More information about the MapServer-users mailing list