Latest Contributions by Tom_Laue

I add something like this if I want to remind myself that I'm about to edit or overwrite critical data...   answer = input("Do you want to continue? (yes/no): ").lower() if answer == "yes" or answer == "y": print("Continuing...") # Add your code here for further execution else: print("Exiting...") sys.exit() ... View more

Yes, that is correct. You need to use update cursor to edit records whereas a search cursor is only used for read-only.  ... View more

To get unique values in a field, I've used.   import arcpy featureClass = r"C:\Default.gdb\Export_Output_10" fields = ['FEATUREID'] listOfValues=[] with arcpy.da.SearchCursor(featureClass, fields) as cursor: for row in cursor: listOfValues.append(row[0])   Then you can convert the list to a set.   # set() to convert list to set set_of_values = set(listOfValues)     Sets can be subtracted from each other to find the differences. Python Set | difference() - GeeksforGeeks ... View more

I'm only aware of the arcpy command to not export attachments within Pro exports. arcpy.env.maintainAttachments = False Other than writing a script to delete all the attachments after they've been downloading, I'm curious what other's have found as a solution. ... View more

@jzcgis wrote: Honestly, half of my workflow would be cut in half if ESRI would allow for multiple geocoders to be used when geocoding with the option to use your local one first and then anything not matched to use the ESRI geocoder. I have been using this methodology. I first try using one of thecounty Geolocators in our service territory. Then if no results are returned, I default back to the ESRI geolocator.     def ESRIGeolocator(addressString): ## STEPS # street address match check # street name match check # call map grid match geocode_result = arcgis.geocoding.geocode(address=addressString, as_featureset=True) #geocode_result = geocode(address=addressString, as_featureset=True) resultNumber = 0 maxMatchScore = 0 matchScore = 0 bestMatchRecordNumber=-1 for geoResult in geocode_result.features: while resultNumber<len(geocode_result.features): #print("made it here") match_type = geocode_result.features[resultNumber].attributes["Addr_type"] matchAddress = geocode_result.features[resultNumber].attributes['Match_addr'] #print("\t"+matchAddress) #print(match_type) if match_type == "StreetAddress" or match_type == "PointAddress" : #print("made it here") #print(geocode_result.features[resultNumber]) matchScore = geocode_result.features[resultNumber].attributes["Score"] if matchScore > maxMatchScore: maxMatchScore = matchScore bestMatchRecordNumber = resultNumber ## print(maxMachScore) ## print(bestMatchRecordNumber) resultNumber +=1 if bestMatchRecordNumber!=-1 or maxMatchScore>93.1: #print(bestMatchRecordNumber) #print(matchScore) #print(new_geocode_result.features[bestNewMatchRecordNumber]) longitude = geocode_result.features[bestMatchRecordNumber].geometry.x latitude = geocode_result.features[bestMatchRecordNumber].geometry.y bestMatchAddress = geocode_result.features[bestMatchRecordNumber].attributes['Match_addr'] #print(bestMatchAddress) matchScore = geocode_result.features[bestMatchRecordNumber].attributes["Score"] #print("\t"+str(matchScore)+"% geocode match") matchType = geocode_result.features[bestMatchRecordNumber].attributes["Addr_type"] #print("\tMatchtype: "+matchType) #print(geocode_result.features[bestMatchRecordNumber]) else: # no point matches #print("No Matches...Trying with other City Names...") maxMatchScore = 0 citiesInCEG = ['Greenwood','Westfield','Zionsville','Mooresville','Crows Nest','Meridian Hills','Rocky Ripple','Southport','Speedway','Indianapolis','Clermont','Linton','Avon','Brownsburg','Plainfield','Beech Grove','Fishers','Carmel'] for cityName in citiesInCEG: #addressWithoutCity = addressString.split(",")[0]+ ", Indiana, USA" newCityAddress = addressString.split(",")[0]+", "+cityName+ ", Indiana, USA" #print(newCityAddress) new_geocode_result = geocode(address=newCityAddress , as_featureset=True) ## resultNumber = 0 #bestNewMatchRecordNumber=0 for new_geoResult in new_geocode_result.features: matchScore = 0 match_type = new_geocode_result.features[resultNumber].attributes["Addr_type"] matchAddress = new_geocode_result.features[resultNumber].attributes['Match_addr'] matchRegion = new_geocode_result.features[resultNumber].attributes['Region'] #print(matchAddress) if match_type == "PointAddress":# and matchAddress.split(",")[0].strip() == addressWithoutCity.split(",")[0].strip() and matchRegion == "Indiana": #print(new_geocode_result.features[resultNumber]) matchScore = new_geocode_result.features[resultNumber].attributes["Score"] #print(matchScore) if matchScore > maxMatchScore: maxMatchScore = matchScore bestNewMatchRecordNumber = resultNumber if bestMatchRecordNumber!=-1: #print(new_geocode_result.features[bestNewMatchRecordNumber]) longitude = new_geocode_result.features[bestNewMatchRecordNumber].geometry.x latitude = new_geocode_result.features[bestNewMatchRecordNumber].geometry.y bestMatchAddress = new_geocode_result.features[bestNewMatchRecordNumber].attributes['Match_addr'] #print(bestMatchAddress) matchScore = new_geocode_result.features[bestNewMatchRecordNumber].attributes["Score"] #print("\t"+str(matchScore)+"% geocode match") matchType = new_geocode_result.features[bestNewMatchRecordNumber].attributes["Addr_type"] #print("\tMatchtype: "+matchType) resultNumber +=1 ##street match search if bestMatchRecordNumber==-1:# or maxMatchScore<93.1: citiesInCEG = ['Indianapolis','Greenwood','Zionsville','Mooresville','Crows Nest','Meridian Hills','Rocky Ripple','Southport','Speedway','Clermont','Linton','Avon','Brownsburg','Plainfield','Beech Grove','Fishers','Carmel'] resultNumber = 0 maxMatchScore = 0 bestMatchRecordNumber = -1 ########################################### ### street name search ########################################### resultNumber =0 for geoResult in geocode_result.features: ## print("made it here") ## print(resultNumber) match_type = geocode_result.features[resultNumber].attributes["Addr_type"] if match_type == "StreetName": #print(geocode_result.features[resultNumber]) matchScore = geocode_result.features[resultNumber].attributes["Score"] if matchScore > maxMatchScore: maxMatchScore = matchScore bestMatchRecordNumber = resultNumber longitude = geocode_result.features[bestMatchRecordNumber].geometry.x latitude = geocode_result.features[bestMatchRecordNumber].geometry.y matchAddress = geocode_result.features[bestMatchRecordNumber].attributes['Match_addr'] #print(matchAddress) matchScore = geocode_result.features[bestMatchRecordNumber].attributes["Score"] #print(str(matchScore)+"% geocde match") matchType = geocode_result.features[bestMatchRecordNumber].attributes["Addr_type"] #print("\tMatchtype: "+matchType) #print(bestMatchRecordNumber) if bestMatchRecordNumber == -1: for cityName in citiesInCEG: #addressWithoutCity = addressString.split(",")[0]+ ", Indiana, USA" newCityAddress = addressString.split(",")[0]+", "+cityName+ ", Indiana, USA" #print(newCityAddress)#+" K") new_geocode_result = geocode(address=newCityAddress , as_featureset=True) ## newResultNumber = 0 #bestNewMatchRecordNumber=0 for new_geoResult in new_geocode_result.features: newResultNumber = 0 match_type = new_geocode_result.features[newResultNumber].attributes["Addr_type"] matchAddress = new_geocode_result.features[newResultNumber].attributes['Match_addr'] matchRegion = new_geocode_result.features[newResultNumber].attributes['Region'] #print(matchAddress) if match_type == "StreetAddress" or match_type == "StreetName":# and matchAddress.split(",")[0].strip() == addressWithoutCity.split(",")[0].strip() and matchRegion == "Indiana": #print(new_geocode_result.features[resultNumber]) matchScore = new_geocode_result.features[newResultNumber].attributes["Score"] #print(matchScore) if matchScore > maxMatchScore: maxMatchScore = matchScore bestNewMatchRecordNumber = newResultNumber #print(new_geocode_result.features[bestNewMatchRecordNumber]) longitude = new_geocode_result.features[bestNewMatchRecordNumber].geometry.x latitude = new_geocode_result.features[bestNewMatchRecordNumber].geometry.y bestMatchAddress = new_geocode_result.features[bestNewMatchRecordNumber].attributes['Match_addr'] #print(bestMatchAddress) matchScore = new_geocode_result.features[bestNewMatchRecordNumber].attributes["Score"] #print("\t"+str(matchScore)+"% geocde match") matchType = new_geocode_result.features[bestNewMatchRecordNumber].attributes["Addr_type"] #print("\tMatchtype: "+matchType) #print(latitude,longitude) newResultNumber +=1 resultNumber +=1 #print(resultNumber) if maxMatchScore>93: ## try: latlong = str(latitude)+"+"+str(longitude) print(bestMatchAddress) #matchScore = new_geocode_result.features[bestNewMatchRecordNumber].attributes["Score"] print("\t"+str(matchScore)+"% geocde match") #matchType = new_geocode_result.features[bestNewMatchRecordNumber].attributes["Addr_type"] print("\tMatchtype: "+matchType) print(latitude,longitude) return latlong ## except: ## latlong = None else: latlong = None FullAddress = "500 S Capitol Ave" city = "Indianapolois" address = FullAddress address = address.replace(" ","+") url = "https://xmaps.indy.gov/arcgis/rest/services/Locators/IndyStreets/GeocodeServer/findAddressCandidates?Street="+address+"&City="+city+"&ZIP=&Single+Line+Input=&category=&outFields=&maxLocations=&outSR=4326&searchExtent=&location=&distance=&magicKey=&f=json" ctx = ssl.create_default_context() ctx.check_hostname = False ctx.verify_mode = ssl.CERT_NONE html = urllib.request.urlopen(url, context=ctx).read() my_bytes_value = html.decode().replace("'", '"') jsonResult = json.loads(my_bytes_value) #print(jsonResult) my_bytes_value = html.decode().replace("'", '"') jsonResult = json.loads(my_bytes_value) #print(jsonResult) if jsonResult['candidates']: print("result found") #print(jsonResult['candidates'][0]) x = jsonResult['candidates'][0]['location']['x'] y = jsonResult['candidates'][0]['location']['y'] score = jsonResult['candidates'][0]['score'] print(jsonResult['candidates'][0]['address']) #matchScore = new_geocode_result.features[bestNewMatchRecordNumber].attributes["Score"] print("\t"+str(score)+"% geocde match") #matchType = new_geocode_result.features[bestNewMatchRecordNumber].attributes["Addr_type"] #print("\tMatchtype: "+matchType) print(latitude,longitude) latitude = x longitude = y latlong_dict[WO] = str(latitude)+"+"+str(longitude) else: addressString = FullAddress + ", " + city + ", IN" print(addressString) latlong = ESRIGeolocator(addressString)     ... View more

Would Mapbox work? This is what I include in automated emails to plot points on a Mapbox street basemap in the email body via Python.       (replace the TOKEN with your API Mapbox TOKEN) EmailBody += "\n<b>Map of Ticket Location: </b><br>" EmailBody += r'''<img width="600" src="https://api.mapbox.com/styles/v1/mapbox/streets-v11/static/pin-m-marker+285A98('''+str(x)+''','''+str(y)+'''/'''+str(x)+''','''+str(y)+''',13,0/600x300@2x?access_token=TOKEN">'''   ... View more

I have not had a chance to implement this yet. ... View more

I really need this too!!! Exact same issue. ... View more

  I’ve been encountering similar issues.  I have found that ArcGIS Pro (which we license through AGOL) automatically logs itself out of every 90 days. I have many Arcpy scripts running via Task Scheduler 24/7.   Also, periodically (usually once a day at a specific time) scripts will fail saying “RuntimeError: The Product License has not been initialized.”  No idea what is causing this error as this script runs all day every 30 mins and never gets this error the rest of the day.     ESRI Tech Support sent me this suggestion which we have yet to implement…     -------------------------------------------   We’ve done a bit of digging into our documentation and might have a solution for you. The bad news is, it doesn’t look like there’s a good way to handle the Pro license refresh through Python. The good news is, there is a promising option for handling it in your environment directly. Please see below for a configuration step-by-step that should solve your logout issue:   Step 1: Environmental settings Convert your Named license to a Single Use license for this instance of Pro (if you haven’t already). See https://pro.arcgis.com/en/pro-app/latest/get-started/convert-named-user-licenses-to-single-use.htm for more details on how to configure this. This should allow offline access for Pro. Make sure “Sign me in automatically” remains checked when you log into Pro on this machine. Just in case.   Step 2: Caching local user (src: https://community.esri.com/t5/python-questions/arcpy-fails-to-import-when-not-logged-into-host-machine/td-p/553399) Log into the server that will host the script using the username that will be running the script (i.e. the “owner” in task scheduler). Once logged in, open ArcGIS Pro (this will create an esri profile cookie in your AppData folder, that ArcGIS will be able to use to recognize the user). Make sure your license and configurations are set according to step 1 and close Pro.   Step 3: Just in case… (src: https://community.esri.com/t5/arcgis-pro-questions/running-a-scheduled-task-on-a-server-machine/m-p/94134) If the issue persists after running through steps 1 and 2, it means that your Pro configuration requires the application to be open while running scripts against its library, for whatever reason. If this is the case, and only after completing steps 1 and 2, we recommend you include a task before the script in your task scheduler that will open ArcGIS Pro, to remove any weird library locks. This is an edge case but does appear to affect some users, so I wanted to include it in the outline.     ... View more

Hi @Anonymous User  – thanks for the quick reply!   This is for a street cut restoration workflow which has field crews visit a site on two different occasions…   I am using Survey123 for field crews to complete their initial street restoration form including taking a picture of the initial restoration. Another crew will go out later to complete the final restoration and click the existing point (preferably in Field Maps) and fill in the final restoration width, restoration length, restoration notes (optional) and take a final picture. So all points will ultimately have 2 photos (initial and final photo).  Trying to edit existing records with Survey123 is not easy but doable.  I’d rather the second field crew work exclusively in Field Maps since they can see the points they need to visit in the map.  I can hide all the fields except the ones they need to fill out and make those required for them, but I can’t require them to take a picture which we need for our records. I don’t expect any offline editing. ... View more

Christina Miller I was trying to do the exact same thing as you when I found your post.  I’m just shocked ESRI doesn’t provide these kinds of reports out of the box. I spent many hours trying to figure out how to read the webmapID in a given app/dashboard with Python (which I could see in AGOL Assistant JSON view).  Luckily, I eventually found the solution: See line In [13]: https://developers.arcgis.com/python/guide/cloning-content/     Attached is the script I developed.  It searches your ArcGIS Online for any web apps and dashboards and will then list each of those items and its details, the webmap associated to it and its details, and then every layer in that map and it's URL. ... View more