Re: [OpenMap Users] Rotation of map in Lambert projection

From: Rob Ratcliff <rrr6399_at_futuretek.com>
Date: Tue, 13 Feb 2007 10:05:08 -0600

Speaking of the Lambert projection, shouldn't one be able to use this
projection to convert from State-plane coordinates to Lat/Lon and
vice-versa?

I tried it using this approach:

            LatLonPoint center = new LatLonPoint(30.1167,-100.3333);
            LambertConformal lambert = new
LambertConformal(center,1000000.0f, 500,500,
                                        -100.3333333333333,
30.11666666666667, 31.88333333333333,
                                        29.66666666666667,700000,3000000);
            Point2D pt = new Point2D.Double();
            lambert.LLToWorld(30.015643,-98.179725, pt);
            System.out.println(pt.toString());
and get this:
            Point2D.Double[0.039672363330592424, 1.6755878139863907]

I should be getting numbers for easting and northing in the hundreds of
thousands to millions.

Does anybody know how to transform these numbers into Easting and
Northing in meters or feet?

(I need to be able to convert back and forth in order to display some
imagery data from a server that expects the bounding box to be specified
in the state-plane coordinate system.)

Thanks

Rob


p.s. Here is the definition of the SP coordinate system I'm using:

/*
 * 32139 NAD_1983_StatePlane_Texas_Central_FIPS_4203
 *
PROJCS["NAD_1983_StatePlane_Texas_Central_FIPS_4203",GEOGCS["GCS_North_American_1983",DATUM["D_North_Ameri
 *
can_1983",SPHEROID["GRS_1980",6378137,298.257222101]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]
 *
],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["False_Easting",700000],PARAMETER["False_Northing",3000000],
 *
PARAMETER["Central_Meridian",-100.3333333333333],PARAMETER["Standard_Parallel_1",30.11666666666667],PARAMETE
 *
R["Standard_Parallel_2",31.88333333333333],PARAMETER["Latitude_Of_Origin",29.66666666666667],UNIT["Meter",1]]
 *
http://www.cnr.colostate.edu/class_info/nr502/lg3/datums_coordinates/spcs.html
 * lines.append("<ENVELOPE minx='3112576.5229230225'
miny='10112154.75590226' " +
 * "maxx='3112956.9819749305' maxy='10112501.491367795' " +
 * "name='Initial_Extent' />");
*/

--
[To unsubscribe to this list send an email to "majdart_at_bbn.com"
with the following text in the BODY of the message "unsubscribe openmap-users"]
Received on Tue Feb 13 2007 - 12:07:46 EST

This archive was generated by hypermail 2.3.0 : Tue Mar 28 2017 - 23:25:08 EDT