Latest Contributions by JohnLay

I'm actually trying to simplify this: arcpy.CreateFeatureDataset_management(Test_gdb, "DoesItWork", "PROJCS['NAD_1983_2011_StatePlane_North_Carolina_FIPS_3200_Ft_US',GEOGCS['GCS_NAD_1983_2011',DATUM['NAD_1983_2011',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Lambert_Conformal_Conic'],PARAMETER['False_Easting',2000000.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',-79.0],PARAMETER['Standard_Parallel_1',34.3333333333333],PARAMETER['Standard_Parallel_2',36.1666666666667],PARAMETER['Scale_Factor',1.0],PARAMETER['Latitude_Of_Origin',33.75],UNIT['Foot_US',0.3048006096012192]],VERTCS['NAVD_1988 Geoid12a',VDATUM['North_American_Vertical_Datum_1988'],PARAMETER['Vertical_Shift',0.0],PARAMETER['Direction',1.0],UNIT['Foot_US',0.3048006096012192]];-121841900 -93659000 3048.00609601219;-100000 10000;-100000 10000;3.28083333333333E-03;0.001;0.001;IsHighPrecision") into something that is a little more manageable and less prone to error. If I replace arcpy.CreateFeatureDataset_management(Test_gdb, "DoesItWork", "PROJCS [ ... ],VERTCS[ ... ]; ... with StatePlane = arcpy.SpatialReference(103122) arcpy.CreateFeatureDataset_management(Test_gdb, "DoesItWork", "PROJCS[StatePlane], ,VERTCS[ ... ]; ... I'm halfway there. I just need to find a way to simplify the vertical spatial reference. I am trying to write a script that will create a feature dataset in order to create terrains. I have to do this 100 times for each of our 100 counties. But, I want the code to be easily manipulated so that it can be absorbed into other projects. ... View more

OK, This almost does what I was looking for (Joshua Bixby‌ and James Crandall‌ I just haven't gotten to your examples yet. James--had some trouble with installing Panda, but am square now) I'm a little lost with it though. Please walk me through the bits I'm missing so that I may apply the info instead of just copy it. magicNumberSet = set([16,17,18]) buildingDict = {}  searchRows = arcpy.da.SearchCursor('L_DAMAGE_RESULTS_WIND', ["BLDG_ID", "HAZARD_ID"])  for searchRow in searchrows:      buildingId, hazardId = searchRow      if buildingId in buildingDict:          buildingDict[buildingId].add(hazardId)      else:          buildingDict[buildingId] = set([hazardId])  yesList = [buildingId for buildingId in buildingDict if magicNumberSet.issubset(buildingDict[buildingId])]  noList = [set(buildingDict.keys()).difference(yesList)] Chris, I get a little lost around line 9. if buildingDict[buildingId].add(hazardId) is appending the value set, I assume is buildingDict[buildingId] = set([hazardId]) creating the first instance. I'm not really sure what "=" is supposed to mean here. My brain automatically goes to one is equal to the other which can't be the case. In line 10, for each key in the dictionary if the value set is a subset of ([16,17,18]) add it to the list. This would mean that sets ([17,18]) and ([17,18,19]) would be excluded from the list, but ([16,17,18,19]) would be added. Using magicNumberSet.issuperset(buildingDict[buildingId]) would mean the reverse is true. By replacing it with magicNumberSet==set(buildingDict[buildingId]) I would essentially be saying that the sets must be equal before being added to the list. Correct? I need to check that there are always 3 buildings and only 3 buildings with a hazard ID of 16, 17, and 18. If there are only 2 buildings or 4 buildings whatever the hazard value, the table fails. Does the order the values appear in the set matter? ... View more

It's not quite that the building must have all three, is is more like there must be 3 buildings with the same ID that each have one of the three hazards. the table would look like this: BLDG_ID     HAZARD_ID        OTHER FIELDS 37013          16                        other info unique to HAZ_ID 16 37013          17                        other info unique to HAZ_ID 17 37013          18                        other info unique to HAZ_ID 18 37014          16                        other info unique to HAZ_ID 16 ... Like I mentioned to Joshua Bixby‌ above, I will play with this some later this afternoon. Thank you both for your suggestions. ... View more

I tried building the dictionary with just the Building ID, but the result was a single Building ID key and a single Hazard ID. I probably built it wrong. But like I said originally, I'm still wrapping my head around dictionaries. e.g.: dict = {1:(u'37013',16), 2: (u'37013', 17), 3: (u'37013', 18), 4: (u'37014', 16), 5: (u'37014', 17), 6: (u'37014', 18),...} became dict = {u'37013':18, u'37014':18,...} so even if the original code worked it would fail b/c there was no Hazard ID 16 or 17. I will give this a look-see a little later this afternoon. I will also give @Chris Snyder' s a look as well. ... View more

Actually, the dictionary would look more like this: dict = {1:(u'37013',16), 2: (u'37013', 17), 3: (u'37013', 18), 4: (u'37014', 16), 5: (u'37014', 17), 6: (u'37014', 18),...} # corresponding with Building ID: Hazard ID. I need to be able to identify that Building 37013 has a row for Hazard ID 16, a row for Hazard 17, and one for Hazard 18. If the table does not meet this criteria, the table fails the check and I need to log the Building ID. ... View more

Thanks James, I haven't had the opportunity to play with numpy yet. I will definitely give this a look. But for the sake of expediency, I'm going to have to go with the Joshua Bixby's solution. ... View more

Yes, but that just makes it too logical and simple . ... View more

Awesome! Thanks. ... View more

da cursors with dictionaries ROCK! From 4 hours to 9 minutes! Thanks Caleb and Richard! The updated bit of code looks like this now: FieldNameList = ["CATTLE", "DAIRY", "POULTRY", "SWINE"] Iteration = 0 Iteration2 = 0 if TABLETYPE == "Animal Operations":     for JoinFile in UPDATELIST_join:         if Iteration > 2:             Iteration = 0             Iteration2 = Iteration2 + 1         layer = CENSUSLayerLIST[Iteration]         ANOPSFIELD = FieldNameList[Iteration2]         readList = ["GEOID10_1", "NumCount"]         valueDict = {r[0]:(r[1:]) for r in arcpy.da.SearchCursor(JoinFile, readList)}         updateList =  ["GEOID10_1", ANOPSFIELD]         with arcpy.da.UpdateCursor(layer, updateList) as cursor:             for row in cursor:                 GEOIDval = row[0]                 if GEOIDval in valueDict:                     row[1] = valueDict[GEOIDval][0]                     cursor.updateRow(row)         Iteration = Iteration + 1 else:     for layer in CENSUSLayerLIST:         JoinFile = UPDATELIST_join[Iteration]         readList = ["GEOID10_1", "NumCount"]         valueDict = {r[0]:(r[1:]) for r in arcpy.da.SearchCursor(JoinFile, readList)}         updateList =  ["GEOID10_1", CENSUSFIELDTRANS]         with arcpy.da.UpdateCursor(layer, updateList) as cursor:             for row in cursor:                 GEOIDval = row[0]                 if GEOIDval in valueDict:                     row[1] = valueDict[GEOIDval][0]                     cursor.updateRow(row)         Iteration = Iteration + 1 Now, just so that I understand what is happening here, valueDict = {r[0]:(r[1:]) for r in arcpy.da.SearchCursor(JoinFile, readList)} is building an array from which with arcpy.da.UpdateCursor(layer, updateList) as cursor:     for row in cursor:         GEOIDval = row[0]         if GEOIDval in valueDict:             row[1] = valueDict[GEOIDval][0]             cursor.updateRow(row) the current cursor row searches for a match. And row[1] = valueDict[GEOIDval][0] means the value of the update cursor row's second field (row[1]) now equals the matched object's next value from the matched part (valueDict[GEOIDval][0]). Correct? ... View more

Here is the full code. I'm going to investigate the da cursors and dictionaries method. I've been aware that calculatefield management is very inefficient for a while; I just couldn't figure out another way to do the same thing. Fieldmappings are new to me, so if there is a more elegant way to write that bit, I'd be interested to learn. SCRIPTPATH = sys.path[0] ROOTFOLDER = os.path.dirname(SCRIPTPATH) CENSUS_FGDB = os.path.join(ROOTFOLDER, "CENSUS.gdb") Census_Blocks = os.path.join(ROOTFOLDER, "CENSUS.gdb\CENSUS_BLOCKS") Census_Group = os.path.join(ROOTFOLDER, "CENSUS.gdb\CENSUS_GROUP") Grid = os.path.join(ROOTFOLDER, "CENSUS.gdb\CENSUS_GRID") CENSUSLIST = [Census_Blocks, Census_Group, Grid] TABLETYPE = arcpy.GetParameterAsText(0) UPDATELAYER = arcpy.GetParameterAsText(1) UPDATEFIELD = arcpy.GetParameterAsText(2)# optional CENSUSFIELD = arcpy.GetParameterAsText(3)# optional usrName = os.getenv('USERNAME') HOME = r'C:\Users\%s' % usrName if CENSUSFIELD == "Buildings":     CENSUSFIELDTRANS = "BLDCNT" elif CENSUSFIELD == "Chemical Sites":     CENSUSFIELDTRANS = "CHEMICAL" elif CENSUSFIELD == "Correctional Facilities":     CENSUSFIELDTRANS = "CORRECT" elif CENSUSFIELD == "Dams":     CENSUSFIELDTRANS = "DAMS" else:     pass if UPDATEFIELD =='#' or not UPDATEFIELD:     UPDATEFIELD = "NOTUSED" DEFAULTGDB = os.path.join(HOME, "Documents", "ArcGIS", "Default.gdb") UPDATELAYERTemp = (os.path.join(DEFAULTGDB, "UPDATELAYERTemp")) desc = arcpy.Describe(UPDATELAYER) type = desc.shapeType arcpy.CopyFeatures_management(UPDATELAYER, UPDATELAYERTemp) # List Fields and delete unnecessary fields = arcpy.ListFields(UPDATELAYERTemp) fieldNameList = [] for field in fields:     if not field.required and not field.name == UPDATEFIELD:         fieldNameList.append(field.name) arcpy.DeleteField_management(UPDATELAYERTemp, fieldNameList) arcpy.AddField_management(UPDATELAYERTemp,"Number","SHORT") arcpy.CalculateField_management(UPDATELAYERTemp,"Number",1,"PYTHON_9.3") # Define Animal Operations ANOPCATTLE = (os.path.join(DEFAULTGDB, "ANOPCATTLE")) ANOPDAIRY = (os.path.join(DEFAULTGDB, "ANOPDAIRY")) ANOPPOULTRY = (os.path.join(DEFAULTGDB, "ANOPPOULTRY")) ANOPSWINE = (os.path.join(DEFAULTGDB, "ANOPSWINE")) exp1 = "\"" + UPDATEFIELD + "\" LIKE '%Dry%' OR \"" + UPDATEFIELD + "\" LIKE '%Beef%'" exp2 = "\"" + UPDATEFIELD + "\" LIKE '%Dairy%' OR \"" + UPDATEFIELD + "\" LIKE '%Milk%'" exp3 = "\"" + UPDATEFIELD + "\" LIKE '%Poultry%'" exp4 = "\"" + UPDATEFIELD + "\" LIKE '%Swine%'" # Define Census Layers CENSUS_Blocks_Layer = (os.path.join(DEFAULTGDB, "CENSUS_Blocks_Layer")) CENSUS_Group_Layer = (os.path.join(DEFAULTGDB, "CENSUS_Group_Layer")) CENSUS_Grid_Layer = (os.path.join(DEFAULTGDB, "CENSUS_Grid_Layer")) CENSUSLayerLIST = [CENSUS_Blocks_Layer, CENSUS_Group_Layer, CENSUS_Grid_Layer] # Make Feature Layers from Census Iteration = 0 for layer in CENSUSLIST:     Newlayer = CENSUSLayerLIST[Iteration]     arcpy.MakeFeatureLayer_management(layer, Newlayer)     Iteration = Iteration + 1 #Define Spatial Join outputs for UPDATELIST_join = [] CENSUS = ["Blocks", "Group", "Grid"] if TABLETYPE == "Animal Operations":     arcpy.Select_analysis(UPDATELAYERTemp, ANOPCATTLE, exp1)     arcpy.Select_analysis(UPDATELAYERTemp, ANOPDAIRY, exp2)     arcpy.Select_analysis(UPDATELAYERTemp, ANOPPOULTRY, exp3)     arcpy.Select_analysis(UPDATELAYERTemp, ANOPSWINE, exp4)     JOINLAYERList = [ANOPCATTLE, ANOPDAIRY, ANOPPOULTRY, ANOPSWINE] else:     JOINLAYERList = [UPDATELAYERTemp] if TABLETYPE == "Lagoons":     field_name = UPDATEFIELD     CENSUSFIELDTRANS = "LAGOONS" else:     field_name = "Number" for file in JOINLAYERList:     Iteration = 0     for layer in CENSUSLayerLIST:         Export_Output = file + "_" + CENSUS[Iteration]         UPDATELIST_join.append(Export_Output)         fieldmappings = arcpy.FieldMappings()         fieldmappings.addTable(layer)         fieldmappings.addTable(file)         FieldIndex = fieldmappings.findFieldMapIndex(field_name)         fieldmap = fieldmappings.getFieldMap(FieldIndex)         field = fieldmap.outputField         field.name = "NumCount"         field.aliasName = "NumCount"         fieldmap.outputField = field         fieldmap.mergeRule = "sum"         fieldmappings.replaceFieldMap(FieldIndex, fieldmap)         for field in fieldmappings.fields:             if field.name not in ["NumCount", "GEOID10_1", "GEOID10_3"]:                 fieldmappings.removeFieldMap(fieldmappings.findFieldMapIndex(field.name))         arcpy.SpatialJoin_analysis(layer, file, Export_Output, "#","#", fieldmappings)         with arcpy.da.UpdateCursor(Export_Output, "NumCount") as cursor:             for row in cursor:                 if row[0] is None:                     cursor.deleteRow()         Iteration = Iteration + 1 FieldNameList = ["CATTLE", "DAIRY", "POULTRY", "SWINE"] Iteration = 0 Iteration2 = 0 if TABLETYPE == "Animal Operations":     for JoinFile in UPDATELIST_join:         if Iteration > 2:             Iteration = 0             Iteration2 = Iteration2 + 1         layer = CENSUSLayerLIST[Iteration]         ANOPSFIELD = FieldNameList[Iteration2]         arcpy.AddJoin_management(layer, "GEOID10_1", JoinFile, "GEOID10_1", "KEEP_COMMON")         arcpy.AddMessage("Updateing " + layer + "' " + ANOPSFIELD)         arcpy.CalculateField_management(layer, ANOPSFIELD, "!NumCount!", "PYTHON_9.3")         arcpy.RemoveJoin_management(layer)         Iteration = Iteration + 1 else:     for layer in CENSUSLayerLIST:         JoinFile = UPDATELIST_join[Iteration]         arcpy.AddJoin_management(layer, "GEOID10_1", JoinFile, "GEOID10_1", "KEEP_COMMON")         arcpy.AddMessage("Updateing " + layer + "' " + CENSUSFIELDTRANS)         arcpy.CalculateField_management(layer, CENSUSFIELDTRANS, "!NumCount!", "PYTHON_9.3")         arcpy.RemoveJoin_management(layer)         Iteration = Iteration + 1 ... View more

That worked. Thanks. ... View more