topic Re: Variable question in Transportation Questions

Variable question

Anonymous User — Sat, 11 Dec 2021 17:13:54 GMT

Original User: luke.kaim

I am getting errors when I run this bit of code and can not figure out why. Is this because of the variables I am using?

import arcpy, sys, os, re, time
from arcpy import env
from arcpy.sa import *
import arcpy
print "Start time = " + time.asctime()
arcpy.env.overwriteOutput = True
arcpy.env.workspace = "C:\\Users\\Luke Kaim\\Documents\\university_of_maine\\research4\\test7"

# Check out any necessary licenses
arcpy.CheckOutExtension("spatial")




APPROACH_ContourLine_Clip="C:\\Users\\Luke Kaim\\Documents\\university_of_maine\\research4\\test5\\APPROACH_HARBOR_DEPCNT_line.shp"

# Check out any necessary licenses
arcpy.CheckOutExtension("spatial")

APPROACH_ContourLine_Project_shp = "APPROACH_ContourLine_Project.shp"
Generalize_ContourLine5meter_shp = "Generalize_ContourLine5meter.shp"
allinputpoints_shp = "allinputpoints2.shp"
test = "test.shp"





# Process: Simplify Line
arcpy.SimplifyLine_cartography(APPROACH_ContourLine_Clip, Generalize_ContourLine5meter_shp, "POINT_REMOVE", "5 Meters", "RESOLVE_ERRORS", "NO_KEEP", "CHECK")



# Process: Feature Vertices To Points
arcpy.FeatureVerticesToPoints_management(Generalize_ContourLine5meter_shp, test, "ALL")


# Process: Add Field (4)
arcpy.AddField_management(test, "input", "SHORT", "", "", "", "", "NON_NULLABLE", "NON_REQUIRED", "")

# Process: Calculate Field (4)
arcpy.CalculateField_management(test, "input", "4", "PYTHON_9.3", "")


print "start projection"


# Process: Project (5)
arcpy.Project_management(test, APPROACH_ContourLine_Project_shp, "PROJCS['NAD_1983_UTM_Zone_19N',GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',-69.0],PARAMETER['Scale_Factor',0.9996],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]", "NAD_1983_To_WGS_1984_1", "GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]")

print "finsihed projection"

Re: Variable question

Anonymous User — Tue, 06 Dec 2011 19:09:30 GMT

Original User: mzcoyle

Importing arcpy twice can't be good for your health.

Other than that, try posting the errors you got and that may help narrow it down.

Re: Variable question

Anonymous User — Tue, 06 Dec 2011 19:22:45 GMT

Original User: luke.kaim

Thank you for the fast reply. I deleted import arcpy. Feature to vertices is not converting all line segments to points though.

This is the error I am getting. I also am attaching the data as well.

Traceback (most recent call last):
File "C:/Users/Luke Kaim/Documents/university_of_maine/research4/contourlinesTEST.py", line 48, in <module>
arcpy.Project_management(test, APPROACH_ContourLine_Project_shp, "PROJCS['NAD_1983_UTM_Zone_19N',GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',-69.0],PARAMETER['Scale_Factor',0.9996],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]", "NAD_1983_To_WGS_1984_1", "GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]")
File "C:\Program Files\ArcGIS\Desktop10.0\arcpy\arcpy\management.py", line 4927, in Project
raise e
ExecuteError: Failed to execute. Parameters are not valid.
ERROR 000365: Invalid geographic transformation.
Failed to execute (Project).

Thank you
Luke Kaim (SIE)
Lucas.Kaim@maine.edu
(914)263-7866
�??Everything should be made as simple as possible, but no simpler�?� (Albert Einstein).

Re: Variable question

ChrisSnyder — Tue, 06 Dec 2011 19:30:48 GMT

Transforming between NAD83 to WGS84 requires a datum transformation method, and it looks as if you are not supplying that parameter. If you are in the lower 48, "NAD_1983_To_WGS_1984_5" should work.

See: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//003r00000010000000.htm

Re: Variable question

Anonymous User — Tue, 06 Dec 2011 19:33:59 GMT

Original User: mzcoyle

I don't see a geographic transformation defined, this is not optional in the projection you are doing.

Error code 000365

Project tool syntax

You may need this resource to help you choose the correct transformation.

Re: Variable question

Anonymous User — Tue, 06 Dec 2011 20:18:30 GMT

Original User: luke.kaim

That still does not work though.

# Process: Project (5)
arcpy.Project_management(test, APPROACH_ContourLine_Project_shp, "PROJCS['NAD_1983_UTM_Zone_19N',GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',-69.0],PARAMETER['Scale_Factor',0.9996],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]", "NAD_1983_To_WGS_1984_1", "GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]")

Re: Variable question

MathewCoyle — Tue, 06 Dec 2011 20:35:34 GMT

Try executing the projection you want in Arc using the Project tool, and copy the python snippet after it is run to see where there are differences.

Also, avoid using spaces in your file names or paths, that causes nothing but problems.

Re: Variable question

Anonymous User — Tue, 06 Dec 2011 20:36:38 GMT

Original User: luke.kaim

That is what I did.

Re: Variable question

MathewCoyle — Tue, 06 Dec 2011 20:46:20 GMT

Try moving your workspace to something with a shorter path and no spaces.