<html><body><div style="color:#000; background-color:#fff; font-family:arial, helvetica, sans-serif;font-size:12pt"><div><span><div>Mikael and Carl,</div><div><br></div><div>From the <span style="font-weight: bold; font-style: italic;">correct</span> shape file, the point touching the bounding box on the Eastern side is 2345137.5423228294, 1430204.1095800549. Comparing these with numbers from proj with the lat/long inputs, the easting is within 0.1%, and the northing is within 1%, which I believe is close enough (for points picked from on line maps).</div><div><br></div><div>Mikael wrote:</div><div><br></div><div>> you could try this Proj.4 definition: +init=epsg:2992</div><div><br></div><div>The following two return the exactly same numbers when I input -116.463726 45.614818:</div><div><br></div><div>proj +proj=lcc +lon_0=-120.5 +lat_0=41.75 +lat_1=43.0 +lat_2=45.5 +x_0=400000 +y_0=0.0 +units=ft</div><div><br></div><div>proj
+init=epsg:2992</div><div><br></div><div>Now I have to figure out how to merge the change file coordinates into the original file!</div><div><br></div><div>Thanks for solving the Lambert Conic Conformal (2SP) problem!</div><div><br></div><div>Rick</div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "><br></div></span></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "><br></div> <div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "> <div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; "> <div dir="ltr"> <font size="2" face="Arial"> <hr size="1"> <b><span style="font-weight:bold;">From:</span></b> Mikael Rittri <Mikael.Rittri@carmenta.com><br> <b><span style="font-weight: bold;">To:</span></b> PROJ.4 and general Projections Discussions <proj@lists.maptools.org> <br> <b><span style="font-weight: bold;">Sent:</span></b> Thursday, August
16, 2012 6:42 PM<br> <b><span style="font-weight: bold;">Subject:</span></b> Re: [Proj] request help with Lambert Conic Conformal (2SP)<br> </font> </div> <br>
<div id="yiv922366636">
<div><div>Carl,</div><div>In fact, I don't think you are missing anything. </div><div><br></div><div>I made a quick back-of-an-envelope calculation - not using</div><div>proj, just a plain calculator - starting with the projection</div><div>origin, assuming 111120 meters per degree of latitude,</div><div>and <span class="yiv922366636Apple-style-span" style="">111120 * cos(lat) meters per degree of longitude. </span></div><div><br></div><div>For the test point, I got </div><div><br></div><div>Easting = 1408985 ft</div><div>Northing = 2341613 ft. </div><div><br></div><div>That's close enough to your (Carl's) result, with my simple method. </div><div><br></div><div>I may have made a mistake, but Rick: could you double-check your expected coordinates - I am pretty sure they are wrong. </div><div><br></div><div>Regards,</div><div><br></div><div>Mikael
Rittri</div><div>Carmenta</div><div>Sweden</div><div>http://www.carmenta.com</div><div><br>16 aug 2012 kl. 05:30 skrev "Carl Godkin" <<a rel="nofollow" ymailto="mailto:cgodkin@gmail.com" target="_blank" href="mailto:cgodkin@gmail.com">cgodkin@gmail.com</a>>:<br><br></div><blockquote type="cite"><div><br><br><div class="yiv922366636gmail_quote">On Wed, Aug 15, 2012 at 8:00 PM, Rick Graves <span dir="ltr"><<a rel="nofollow" ymailto="mailto:gravesricharde@yahoo.com" target="_blank" href="mailto:gravesricharde@yahoo.com"></a><a rel="nofollow" ymailto="mailto:gravesricharde@yahoo.com" target="_blank" href="mailto:gravesricharde@yahoo.com">gravesricharde@yahoo.com</a>></span> wrote:<br><blockquote class="yiv922366636gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div><div style="font-size: 12pt; font-family: arial, helvetica, sans-serif; "><div><div><font face="arial, helvetica, sans-serif">Hey,</font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">The good news is that Oregon is the last one I need to get working. The bad news is I do not have it working yet.</font></div>
<div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">Here is the metadata from the .prj that I got from the state:</font></div><div><font face="arial, helvetica, sans-serif"><br>
</font></div><div><font face="arial, helvetica, sans-serif">PROJCS["NAD_1983_Oregon_Statewide_Lambert_Feet_Intl",</font></div><div><font face="arial, helvetica, sans-serif">GEOGCS["GCS_North_American_1983",</font></div>
<div><font face="arial, helvetica,
sans-serif">DATUM["D_North_American_1983",</font></div><div><font face="arial, helvetica, sans-serif">SPHEROID["GRS_1980",6378137.0,298.257222101]],</font></div><div><font face="arial, helvetica, sans-serif">PRIMEM["Greenwich",0.0],</font></div>
<div><font face="arial, helvetica, sans-serif">UNIT["Degree",0.0174532925199433]],</font></div><div><font face="arial, helvetica, sans-serif">PROJECTION["Lambert_Conformal_Conic"],</font></div><div><font face="arial, helvetica, sans-serif">PARAMETER["False_Easting",1312335.958005249],</font></div>
<div><font face="arial, helvetica, sans-serif">PARAMETER["False_Northing",0.0],</font></div><div><font face="arial, helvetica, sans-serif">PARAMETER["Central_Meridian",-120.5],</font></div><div><font face="arial, helvetica, sans-serif">PARAMETER["Standard_Parallel_1",43.0],</font></div>
<div><font face="arial, helvetica, sans-serif">PARAMETER["Standard_Parallel_2",45.5],</font></div><div><font face="arial, helvetica,
sans-serif">PARAMETER["Latitude_Of_Origin",41.75],</font></div><div><font face="arial, helvetica, sans-serif">UNIT["Foot",0.3048]],</font></div><div><font face="arial, helvetica, sans-serif">VERTCS["NAD_1983",</font></div>
<div><font face="arial, helvetica, sans-serif">DATUM["D_North_American_1983",</font></div><div><font face="arial, helvetica, sans-serif">SPHEROID["GRS_1980",6378137.0,298.257222101]],</font></div><div>
<font face="arial, helvetica, sans-serif">PARAMETER["Vertical_Shift",0.0],</font></div><div><font face="arial, helvetica, sans-serif">PARAMETER["Direction",1.0],</font></div><div><font face="arial, helvetica, sans-serif"><br>
</font></div><div><font face="arial, helvetica, sans-serif">I initialized proj as follows:</font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">proj +proj=lcc +lon_0=-120.5 +lat_0=41.75 +lat_1=43.0 +lat_2=45.5 +x_0=1312335.958005249 +y_0=0.0
+units=ft</font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">I assume something is not right.</font></div><div><font face="arial, helvetica, sans-serif"><br>
</font></div><div><br></div></div></div></div></blockquote><div><br></div><div>I think there are two things wrong: </div><div>1. Your false easting is here given in feet while I believe proj always wants meters. The number is about 400000m. (See <a rel="nofollow" target="_blank" href="http://trac.osgeo.org/proj/wiki/GenParms">http://trac.osgeo.org/proj/wiki/GenParms</a>).</div>
<div>2. You didn't include an ellipsoid or datum so I'm surprised that proj worked at all since I get an error without +datum=NAD83 or +ellps=GRS80.</div><div><br></div><div>However, correcting both of those, I get X=2344859.57 ft, Y=1434192.31 ft, which isn't what you want either. </div>
<div><br></div><div>Now I want to know what I am missing! Let's hope someone cleverer than me answers.</div><div><br></div><div>carl</div></div>
</div></blockquote><blockquote type="cite"><div><span>_______________________________________________</span><br><span>Proj mailing list</span><br><span><a rel="nofollow" ymailto="mailto:Proj@lists.maptools.org" target="_blank" href="mailto:Proj@lists.maptools.org">Proj@lists.maptools.org</a></span><br><span>http://lists.maptools.org/mailman/listinfo/proj</span></div></blockquote></div></div><br>_______________________________________________<br>Proj mailing list<br><a ymailto="mailto:Proj@lists.maptools.org" href="mailto:Proj@lists.maptools.org">Proj@lists.maptools.org</a><br><a href="http://lists.maptools.org/mailman/listinfo/proj" target="_blank">http://lists.maptools.org/mailman/listinfo/proj</a><br><br> </div> </div> </div></body></html>