.prj`` contains informations about the coordinates used to encode the shapefiles."]}, {"cell_type": "code", "execution_count": 7, "metadata": {"ExecuteTime": {"end_time": "2016-11-02T11:50:55.776398", "start_time": "2016-11-02T11:50:50.542212"}}, "outputs": [], "source": ["%load_ext pyensae"]}, {"cell_type": "code", "execution_count": 8, "metadata": {"ExecuteTime": {"end_time": "2016-11-02T11:50:59.565358", "start_time": "2016-11-02T11:50:59.548359"}}, "outputs": [{"data": {"text/html": ["\n", "PROJCS[\"NAD_1983_HARN_StatePlane_Washington_North_FIPS_4601_Feet\",GEOGCS[\"GCS_North_American_1983_HARN\",DATUM[\"D_North_American_1983_HARN\",SPHEROID[\"GRS_1980\",6378137.0,298.257222101]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Lambert_Conformal_Conic\"],PARAMETER[\"False_Easting\",1640416.666666667],PARAMETER[\"False_Northing\",0.0],PARAMETER[\"Central_Meridian\",-120.8333333333333],PARAMETER[\"Standard_Parallel_1\",47.5],PARAMETER[\"Standard_Parallel_2\",48.73333333333333],PARAMETER[\"Latitude_Of_Origin\",47.0],UNIT[\"Foot_US\",0.3048006096012192]],VERTCS[\"NGVD_1929\",VDATUM[\"National_Geodetic_Vertical_Datum_1929\"],PARAMETER[\"Vertical_Shift\",0.0],PARAMETER[\"Direction\",1.0],UNIT[\"Foot_US\",0.3048006096012192]]\n", "
"], "text/plain": [""]}, "execution_count": 9, "metadata": {}, "output_type": "execute_result"}], "source": ["%head Centerline.prj"]}, {"cell_type": "markdown", "metadata": {}, "source": ["This means to be converted into another coordinates system: *ESPG:XXXX*. That's what the website [Prj2EPSG](http://prj2epsg.org/search) does. The module [pyproj](https://pypi.python.org/pypi/pyproj/) does the conversion into longitude / latitude."]}, {"cell_type": "code", "execution_count": 9, "metadata": {"ExecuteTime": {"end_time": "2016-11-02T11:54:33.050369", "start_time": "2016-11-02T11:54:33.026243"}}, "outputs": [{"data": {"text/plain": ["(-109.78898903640629, 48.55518509715834)"]}, "execution_count": 10, "metadata": {}, "output_type": "execute_result"}], "source": ["from pyproj import Proj, transform\n", "p1 = Proj(init='epsg:2926') # returned by Prj2ESPG\n", "p2 = Proj(init='epsg:4326') # longitude / latidue\n", "x0, y0 = 1312487.225652799, 231314.96255882084\n", "transform(p1,p2,x0,y0)"]}, {"cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [], "source": []}], "metadata": {"kernelspec": {"display_name": "Python 3", "language": "python", "name": "python3"}, "language_info": {"codemirror_mode": {"name": "ipython", "version": 3}, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.0"}}, "nbformat": 4, "nbformat_minor": 2}