Latest Contributions by JohannesBierer

The following python code should read out of a text file folder paths and project the shapefiles in this folder. Unfortunately this works only with the last (folder) line? Any suggestions why this doesn't work? # Import arcpy module import arcpy, os fname = r"H:\python_scripte\ProjizierenMulti\test.txt" with open(fname) as f:     for line in f:                 print line                 arcpy.env.workspace = line         workspace = arcpy.env.workspace         arcpy.env.overwriteOutput = True         # Use the ListFeatureClasses function to return a list of         #  shapefiles.         featureclasses = arcpy.ListFeatureClasses()         # Copy shapefiles to a file geodatabase         for fc in featureclasses:             print "Projiziere " + fc             outfile = os.path.splitext(fc)[0] + "_E.shp"                             spatialRef = arcpy.Describe(fc).spatialReference             print spatialRef.name             Ftype = arcpy.Describe(fc).datatype             print Ftype                             if Ftype == "ShapeFile":                                 if spatialRef.name == "DHDN_3_Degree_Gauss_Zone_3":                                                             arcpy.Project_management(fc, outfile, "PROJCS['ETRS_1989_UTM_Zone_32N_8stellen',GEOGCS['GCS_ETRS_1989',DATUM['D_ETRS_1989',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',32500000.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',9.0],PARAMETER['Scale_Factor',0.9996],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]", "DHDN_To_ETRS_1989_3", "PROJCS['DHDN_3_Degree_Gauss_Zone_3',GEOGCS['GCS_Deutsches_Hauptdreiecksnetz',DATUM['D_Deutsches_Hauptdreiecksnetz',SPHEROID['Bessel_1841',6377397.155,299.1528128]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Gauss_Kruger'],PARAMETER['False_Easting',3500000.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',9.0],PARAMETER['Scale_Factor',1.0],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]")                     f.close()‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

## Load Data import arcpy import os workspace = r"R:\Daten\geo_daten_Kopie" out_workspace = r"R:\Karto\Bierer2016\OrdnerRDatenGeodatenBereinigen\SIC_Geodaten.gdb\data" import arcpy import os feature_classes = [] for dirpath, dirnames, filenames in arcpy.da.Walk(workspace, datatype="FeatureClass"): ##, type="Polygon"):     for filename in filenames:         feature_classes.append(os.path.join(dirpath, filename))         name = os.path.join(dirpath, filename)         print name         BaseName = os.path.basename(name)         print BaseName         BaseNameCopy = BaseName.split(".shp")         print BaseNameCopy         arcpy.CopyFeatures_management(filename, os.path.join(out_workspace, BaseNameCopy)) The code leads to this error message: Traceback (most recent call last):   File "H:/python_scripte/ListFolderandShapes/ListFeatures_FoldersSubfolders_CopyFeatures_3.py", line 22, in <module>     arcpy.CopyFeatures_management(filename, os.path.join(out_workspace, BaseNameCopy))   File "C:\Python27\ArcGIS10.2\lib\ntpath.py", line 73, in join     elif isabs(b):   File "C:\Python27\ArcGIS10.2\lib\ntpath.py", line 58, in isabs     return s != '' and s[:1] in '/\\' TypeError: 'in <string>' requires string as left operand, not list Anyone has got an idea for a better solution? ... View more

## Load Data import arcpy import os workspace = r"R:\Daten\geo_daten_Kopie" out_workspace = r"R:\Karto\Bierer2016\OrdnerRDatenGeodatenBereinigen\SIC_Geodaten.gdb\data" import arcpy import os feature_classes = [] for dirpath, dirnames, filenames in arcpy.da.Walk(workspace, datatype="FeatureClass"): ##, type="Polygon"):     for filename in filenames:         feature_classes.append(os.path.join(dirpath, filename))         name = os.path.join(dirpath, filename)         print name         BaseName = os.path.basename(name)         print BaseName         BaseNameCopy = BaseName.split(".shp")         print BaseNameCopy         arcpy.CopyFeatures_management(filename, os.path.join(out_workspace, BaseNameCopy)) The code leads to this error message: Traceback (most recent call last):   File "H:/python_scripte/ListFolderandShapes/ListFeatures_FoldersSubfolders_CopyFeatures_3.py", line 22, in <module>     arcpy.CopyFeatures_management(filename, os.path.join(out_workspace, BaseNameCopy))   File "C:\Python27\ArcGIS10.2\lib\ntpath.py", line 73, in join     elif isabs(b):   File "C:\Python27\ArcGIS10.2\lib\ntpath.py", line 58, in isabs     return s != '' and s[:1] in '/\\' TypeError: 'in <string>' requires string as left operand, not list Anyone has got an idea for a better solution? ... View more

Error message in script which worked correctly days before? Any ideas why? line 16, in <module>     quercy = field + "=" + "'" + Loop1 + "'" TypeError: cannot concatenate 'str' and 'int' objects import arcpy import os input = arcpy.GetParameterAsText(0) clipfeature = arcpy.GetParameterAsText(1) output = arcpy.GetParameterAsText(2) clipfeature_Layer = "clipfeature_Layer" arcpy.env.workspace = output arcpy.env.overwriteOutput = True field = "BFN_NR" inputSp = input.split(";") arcpy.MakeFeatureLayer_management(clipfeature, clipfeature_Layer, "", "", "") cursor = arcpy.SearchCursor(clipfeature) for row in cursor:     arcpy.AddMessage(row.getValue(field))     Loop1 = row.getValue(field)     quercy = field + "=" + "'" + Loop1 + "'"     arcpy.SelectLayerByAttribute_management (clipfeature_Layer, "NEW_SELECTION", quercy)     Sel_shape = "Sel_" + Loop1 + ".shp"     arcpy.CopyFeatures_management(clipfeature_Layer, Sel_shape)     for i in inputSp:         arcpy.AddMessage(i)         arcpy.Clip_analysis(i, Sel_shape, "Clip_" + row.getValue(field) + "_" + i) ... View more

I've a table with paths like: R:\Daten\geo_daten\GEOLOGIE\bodenuebersicht\daten\arcview\bodeneigenschaften\shape\bodeneigenschaften.shp and want to copy this Feature somewhere else and delete in the location where it is? What's wrong? import arcpy, numpy arcpy.env.overwriteOutput = True InPath = r"R:\Karto\Geodaten_Archiv.txt" OutPath = r"R:\Daten\geo_daten\geodaten_archiv\SIC_2016_04_07" filename = open(InPath, 'r') arr = arcpy.da.TableToNumPyArray(InPath, "Datenpfad") for line in arr:     print line     arcpy.CopyFeatures_management(line, OutPath) ... View more