Latest Contributions by DuminduJayasekera

Hi, I have the following python script to create a TEXT column called "SEC" and then I need to perform a Field calculator operation for the SEC column. I want to read the corresponding string in the "mtrs" column and and fill in SEC column. I have performed this task manually and it worked. I need to slice the "mtrs" string like this "S{}".format(int( !mtrs! [17:19]))" and fill in SEC column.  My script below give me an error in line 14. SyntaxError. invalid syntax.  Can somebody help me to fix this and make this work? Thanks in advance, # Import system modules import arcpy from arcpy import env # Set environment settings env.workspace = "H:/Test/Points.gdb" # Set local variables inFeatures = "plsssection_test_p1" fieldName1 = "SEC" length = 20 expression = ""S{}".format(int( !mtrs! [17:19]))" # Execute AddField twice for two new fields arcpy.AddField_management(inFeatures, fieldName1, "TEXT", length, "", "") arcpy.CalculateField_management(inFeatures, "SEC",expression, "PYTHON_9.3")‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

HI, I need to add "_var86" to the output file in the code below. How can I achieve it in the code below.  Thanks in advance. import arcpy, sys from arcpy import env from arcpy.sa import * # Input data source arcpy.env.workspace = r"C:/Users/Desktop/weekly_avg_project" arcpy.env.overwriteOutput = True # Set output folder OutputFolder = r"C:/Users/Desktop/weekly_avg_project/output" # Loop through a list of files in the workspace NCfiles = arcpy.ListFiles("*.nc") for filename in NCfiles: print("Processing: " + filename) inNCfiles = arcpy.env.workspace + "//" + filename fileroot = filename[0:(len(filename)-3)] TempLayerFile = "SM_amount" outRaster = OutputFolder + "//" + fileroot # Process: Make NetCDF Raster Layer outfile = arcpy.MakeNetCDFRasterLayer_md(filename, "var86", "lon", "lat", TempLayerFile, "", "", "BY_VALUE") proj_ras = arcpy.ProjectRaster_management(outfile, "C:/Users/Desktop/weekly_avg_project/project_ras","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]]", "NEAREST", "0.125 0.125", "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]]") # Process: Copy Raster arcpy.CopyRaster_management(proj_ras, outRaster + ".tif", "", "", "", "NONE", "NONE", "")‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more