Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I am attempting to transform data from Lambert Conformal Conic lat/long and I am having some issues with the results.

The shapefile I am attemtping to transform can be found here:

Since I am using proj.4 I used GDAL/srs to convert the WKT to a proj.4 init string. In order to do that I had to change the WKT projection from "Lambert_Conformal_Conic" to "Lambert_Conformal_Conic_SP2".

The WKT I used:


The resulting proj.4 string from the WKT:

"+proj=lcc +lat_1=38.69259533705237 +lat_2=39.02625053604863 +lat_0=40 +lon_0=-77.03638889 +x_0=0 +y_0=0 +ellps=GRS80 +units=m +no_defs"

The proj.4 string I use for WGS84 lat/long:

"+proj=longlat +ellps=WGS84 +no_defs"

When I get a Lambert Conformal Conic point from the shapefile, I transform it to Lat/Long using this line of code:

pj_transform(pj_lcc, pj_latlong, 1, 1, &fX, &fY, NULL );

Where pj_lcc is my Lambert Conformal Conic projection info, pj_latlong is my WGS84 lat/long projection info and fX, fY are the Lambert Conformal Conic point.

An example result from the above transform is:

lcc point: {4842.47, -71287.5} --- long, lat: {-1.34356, 0.686926}

The lat, long should be around 38º51'34"N, 77º02'11"W (center of the shape according to )

Any clues as to what I am doing wrong?

share|improve this question
up vote 3 down vote accepted


ogr2ogr -t_srs EPSG:4326 out.shp in.shp

share|improve this answer
That worked great! Thank you for the quick reply. – afarnham Feb 21 '11 at 20:55

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.