<html><body><div style="color:#000; background-color:#fff; font-family:arial, helvetica, sans-serif;font-size:12pt"><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "><span>Mikael and Carl,</span></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "><span><br></span></div><div><span>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).<br></span></div><div><span><br></span></div><div><span>Mikael wrote:</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;
"><span style="font-family: 'times new roman', 'new york', times, serif; ">> </span><span style="font-family: 'times new roman', 'new york', times, serif; ">you could try this Proj.4 definition: </span><span style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; ">+init=epsg:2992</span></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "><span style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; "><br></span></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; ">The following two return the exactly same numbers when I input -116.463726 45.614818:</div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "><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<br></div><div><br></div><div>proj +init=epsg:2992<br></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; ">Now I have to figure out how to merge the change file coordinates into the original file!</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; ">Thanks for solving the problem!</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; ">Rick</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; "><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> Rick Graves <gravesricharde@yahoo.com>; PROJ.4 and generalProjections Discussions <proj@lists.maptools.org> <br> <b><span style="font-weight: bold;">Sent:</span></b> Thursday, August 16, 2012 9:52 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="yiv337315236">
<div><div>Rick,</div><div>I don't know about the shapefiles, but you could try this Proj.4 definition</div><div><br></div><div> +init=epsg:2992</div><div><br></div><div>which I think will be equivalent to Carl's definition. <br><br>(I'm not sure how you can verify that it works, but maybe you</div><div>can compare the Oregon data location with other data sources.)</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. 14:44 skrev "Rick Graves" <<a rel="nofollow" ymailto="mailto:gravesricharde@yahoo.com" target="_blank" href="mailto:gravesricharde@yahoo.com">gravesricharde@yahoo.com</a>>:<br><br></div><blockquote type="cite"><div><div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); font-family: arial, helvetica, sans-serif; font-size: 12pt; "><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "><span>Mikael and
Carl,</span></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "><span><br></span></div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; ">Maybe there is a problem with the Oregon shape files, or maybe I made an incorrect assumption about the contents.</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; ">The Oregon redistricting web site is down for me at the moment:</div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "><br></div><div><div><a rel="nofollow" target="_blank" href="http://www.leg.state.or.us/redistricting/">www.leg.state.or.us/redistricting/</a></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; ">But the Oregon shape files are available
here:</div><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "><br></div><div>http://www.srrproject.org/resources/redistricting-shapefiles/<br></div><div><br></div><div>(Click on Oregon to get the link. If you do not instantly know which one is Oregon, you can hover the mouse and a tip will display telling which state is under the pointer. Oregon is on the West coast just North of California.)</div><div><br></div><div>The zip for Oregon includes two sets of files, Congressional_Districts.* and Congressional_Area_Changes.*. </div><div><br></div><div>I assumed that I should go with the latter file.</div><div><br></div><div>I ran diff against the two .prj files, and there is no difference.</div><div><br></div><div>However, I just now discovered that the contents of the shapefiles are radically different.
Here are the bounding boxes:</div><div><br></div><div>Congressional_Districts<br></div><div>(183987.16863517463, 88653.9580052495, 2345137.542322829, 1677526.0958005222)<br></div><div><br></div><div>Congressional_Area_Changes<br></div><div>(427663.998031497, 172986.72867454588, 1005007.5570866163, 1488902.549212602)<br></div></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; ">Sorry for being a beginner, but is there something I should know to do with these two files? Or should I be asking Oregon what is the deal?</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; ">Thanks,</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; ">Rick</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; "><br></div> <div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; "> <div style="font-family: 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 <<a rel="nofollow" ymailto="mailto:Mikael.Rittri@carmenta.com" target="_blank" href="mailto:Mikael.Rittri@carmenta.com">Mikael.Rittri@carmenta.com</a>><br> <b><span style="font-weight:bold;">To:</span></b> PROJ.4 and general Projections Discussions <<a rel="nofollow" ymailto="mailto:proj@lists.maptools.org" target="_blank" href="mailto:proj@lists.maptools.org">proj@lists.maptools.org</a>> <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="yiv337315236">
<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="yiv337315236Apple-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><a rel="nofollow" target="_blank" href="http://www.carmenta.com/">http://www.carmenta.com</a></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"></a><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="yiv337315236gmail_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"></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="yiv337315236gmail_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"></a><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"></a><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><a rel="nofollow" target="_blank" href="http://lists.maptools.org/mailman/listinfo/proj">http://lists.maptools.org/mailman/listinfo/proj</a></span></div></blockquote></div></div><br>_______________________________________________<br>Proj mailing list<br><a rel="nofollow" ymailto="mailto:Proj@lists.maptools.org" target="_blank" href="mailto:Proj@lists.maptools.org"></a><a rel="nofollow" ymailto="mailto:Proj@lists.maptools.org" target="_blank" href="mailto:Proj@lists.maptools.org">Proj@lists.maptools.org</a><br><a rel="nofollow"
target="_blank" href="http://lists.maptools.org/mailman/listinfo/proj"></a><a rel="nofollow" target="_blank" href="http://lists.maptools.org/mailman/listinfo/proj">http://lists.maptools.org/mailman/listinfo/proj</a><br><br> </div> </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><a rel="nofollow" target="_blank" href="http://lists.maptools.org/mailman/listinfo/proj">http://lists.maptools.org/mailman/listinfo/proj</a></span></div></blockquote></div></div><br><br> </div> </div> </div></body></html>