Latest Contributions by FilipKrál

Hi Justin, I've come across situations when the names of bands were Layer_1, Layer_2, etc., it does not always have to be Band_1, Band_2, etc. Perhaps try to list the individual bands programmatically as shown here: python - Access individual bands and use them in map algebra - Geographic Information Systems Stack Exchange  The main idea is to set your main raster dataset as a workspace and then use arcpy.ListRasters() function. Filip. ... View more

Hi Jose,  I think there is no harm in publishing the two (or more) services from a single map document. As far as I understand, ArcGIS for Server will convert it in to service definition files (.sd) for each service separately. The services will use the .sd files stored in internal ArcGIS server folders and the original .mxd will not be used by the services directly anyway. See PublishMSDToServer—Help | ArcGIS Desktop  Filip. ... View more

Hi Michael and Joseph, The initial piece of code demonstrated how one would go around solving the problem but there are things one needs to add to make it more bulletproof. If I added them in the first place, the code would be too long and hard to understand.  One important assumption in the initial code is that all feature classes in the geodatabase have the same columns. Why else would you want to join them all into one result? I suspect Joseph is seeing the error because the script encountered a feature class that does not fit this assumption. When it comes to empty feature classes, i.e. feature classes with no rows, the code the `for row in rows` loop should not raise an exception as long fc has all the expected columns. The subsequent `del row` statement may raise an exception though. You can to comment out the `del row` statement. One way to deal with these strange cases is to add some try-except statements and handle the failing feature classes in some other way. I'll leave that up to you. Cheers, Filip. ... View more

Hi Jonathan, You are quite right and I owe you an apology. I am glad you pointed this out because it indeed did not make sense to me why one would have to specify the parameter name of field_names when it can be inferred from the position. Anyway, I see the discussion has moved on quite a bit. Filip. ... View more

Hi Eddie, I haven't used the tools you mention much but perhaps the root of the problem you see could be in what was discussed in this thread: https://community.esri.com/thread/159997  Filip. ... View more

I am just guessing here but it does not sound like you're hitting limit on the volume of a response, if there even is any. You could try to adjust the maximum number of records the service should return from a query to something like 10 or less and try your queries again. Have you tried running Check Geometry and Repair Geometry tools on your polygon feature class in ArcGIS for Desktop? Also, exporting the data to GeoJSON in ArcMap may reveal where the problem is. If your polygons are not massive (i.e. if they don't have many vertices), I would try looking at them in GeoJSON also in the Python window like this: import arcpy fc = r'c:\path\to\polygons.shp' with arcpy.da.SearchCursor(fc, "SHAPE@JSON") as sc:     for row in sc:         print(row[0].__geo_interface__)‍‍‍‍‍‍‍‍‍‍‍‍ See also Converting geometries between GeoJSON, esri JSON, and esri Python , especially the last contribution. Filip. ... View more

Hi Michael,  Jonathan Quinn has a very good point about formatting code on GeoNet but he's not quite right about why you are seeing the error. In fact, the variable that refers to the list of fields is called fieldnames consistently throughout the script. I suspect the error is raised by the call to arcpy.da.SearchCursor. It should not be, but anyway, try to replace with arcpy.da.SearchCursor(fc, fieldnames) as sc:‍ with this: with arcpy.da.SearchCursor(fc, field_names=fieldnames) as sc:‍ Filip. ... View more

Hi Graham, Nice to see the good old British National Grid but boy this is a loaded question. Almost all the advanced topics in one. I suggest some changes in the attachment but you still have a long way to go. Some key points are: Result of a geoprocessing tool run is not the output variable but a result object from which you need to pluck out the appropriate variable (hence the .getOutput(0)). In order to return a value from a geoprocessing python script you need to setParameter.  I would stay away from in_memory workspaces until you get everything working and then you can experiment with it. ArcGIS Server should handle the creation of new workspace and scratch workspace for you and then get rid of it after a few weeks or whatever the settings are. The way it is now, result of the script will be a path to the output near table, but I suspect that will not return the content of the table as a result of the REST service. I thing you would need to read the table with a search cursor into something like a RecordSet and return that or dump it to json, then to text, set the text as an output parameter of type Text and then parse this as JSON on the client side. Well... not to discourage you. I did this kind of thing before and it kinda worked, but we never used it in production. It was like scratching your left ear with your right hand behind your head. You may want to consider publishing both the data and the Generate Near Table tool as services and using these in the client fully over the RESTful interface. Best wishes, Filip. ... View more

Hi Mitch,  Thanks for sharing, my favourite is this one https://github.com/NERC-CEH/arcapi I use it all the time. There is also a GeoNet place for that here ArcAPI‌ Filip. ... View more

Hi Michael, You're not supposed to replace 'geometry' or '\n' with anything. The line with 'geometry' lists names of all the columns where column type is not geometry. The '\n' adds a new line character at the end of every line you write to the output csv file. Joseph Head‌ is right. Some print statements will help you figure out what's going on. I recommend adding the following print statements to line 4, 6, 10, and 12, respectively: print("Listing feature classes...") # line 4 print("Found {0} feature classes".format(len(feature_classes))) # line 6 print(field_names) #line 10 print('Writing to {0}'.format(output_file)) # line 12 #Also insert the following just before line 18 and align it with the 'with' keyword:  print(fc) ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ My guess is that the arcpy.ListFeatureClasses() takes a long time but without the print statements you won't know. It is likely that the whole thing is going to take some time. One way to figure out if the process is still doing anything is to check the size of the output csv file every couple of minutes to see if it is still growing. Filip. ... View more

Hi Abdul, Perhaps you are looking for the Extract Values to Points Tool? Also, your question is about arcpy but you filed it into the space for ArcGIS JS API, which is different thing. Arcpy is the Python interface. If you can, move your question to the Python space so more Python developers can see it. Filip. ... View more

Hi Evelyn, If all your layers are included in a single service, you may want to set the scale range visibility of individual layers directly in the map document and republish the service. In ArcGIS Desktop: http://desktop.arcgis.com/en/arcmap/10.3/map/working-with-layers/displaying-layers-at-certain-map-scales.htm In ArcGIS Online? https://doc.arcgis.com/en/arcgis-online/create-maps/set-visibility.htm When you then add that service to your map as a single layer, the layers should turn on and off as you zoom in and out. If you don't have access to the map document, the server, or the map service in ArcGIS Online, you'll need to do it all in the JavaScript client. Define multiple layers, each displaying only a single layer from the service and set the minScale and maxScale properties as appropriate when you create the layers. Also the setScaleRange(minScale, maxScale) may be also useful. Try calling something like myLayer1.setScaleRange(100000, 5000); before you add it to the map. Regardless of scale range visibility, call myLayer1.show(); and myLayer1.hide(); to turn myLayer1 on and off, respectively. Note that if the layer is not within it's visible scale range, nothing will be shown even if the layer is "on". Hope this makes sense. Filip. ... View more

Hi Marcel, I wouldn't rely on OBJECTID to be strictly sequential and increasing. If you created your integer index field just by recalculating the object ID into it, the effect can be the same as using the object id directly. Use the code snippet to assign unique sequential ids that is in field calculator help. If you already did that, then I am not sure where the problem is. Filip. ... View more

To extend Mitch's answer, if your data is in a geodatabase rather than a folder, try the code below. import arcpy file_name_field = 'FNAME' input_gdb_or_folder = r'c:\my.gdb' arcpy.env.workspace = input_gdb_or_folder feature_classes = arcpy.ListFeatureClasses() for fc in feature_classes:     print(fc) # just so you know what the script is processing         # add field to hold the file name if it does not exist     existing_fields = [f.name for f in arcpy.ListFields(fc)]     if file_name_field not in existing_fields:         arcpy.management.AddField(fc, file_name_field, 'TEXT', field_length=200)         # write the file name into each row of the file name filed     with arcpy.da.UpdateCursor(fc, [file_name_field]) as uc:         for row in uc:             uc.updateRow([str(fc)])     del row, uc‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ Filip. ... View more