Latest Contributions by DarrenWiens2

Hmmm... I'm not sure. I reran your code against my feature class and it works. What coordinate reference system are you using? I was using NAD83 UTM 17, but with a different CRS you may be out of bounds using such a large value for 'max'. Also, did you change the values in your table to meters? The buffer makes buffers of the number of linear units in the dataset, so if your dataset is in meters, the buffers will be meters, if feet, then the buffers will be in feet, etc. Can you also try adding the following line above line 23 to get a count of out_buffs? len(out_buffs) ... View more

Here is a quick way to switch between coordinate reference systems based on a field. Notes: - my zone field is an integer (values are 9 or 10, indicating which NAD83 UTM Zone the point is in). Yours is text, but you can convert it to integer like below (takes zone text value, minus the last character, and converts to integer): cur_wkid = int(wkid_shift + int(row[1][:-1]))‍‍ - the script relies on the wkid number being related to the zone number. In my example, I add 26900 to the zone number to get the wkid number. So, zone 9 becomes wkid = 26909 (the correct wkid for NAD83 UTM Zone 9). In your case, you need to add 31960 to the zone number. So, zone 22 becomes wkid = 31982, and zone 21 becomes wkid = 31981. - change the values for fc, zone_field, x_field, y_field, and wkid_shift to match your data. - you should probably make a back-up copy, just in case things go wrong. Script: >>> fc = 'merge' # your feature class ... fc_sr = arcpy.Describe(fc).spatialReference # get the spatial reference ... zone_field = 'UTM_Zone' # change to your zone field ... x_field = 'Easting' # change to your x field ... y_field = 'Northing' # change to your y field ... wkid_shift = 26900 # change to 31960 ... with arcpy.da.UpdateCursor(fc,['SHAPE@',zone_field,x_field,y_field],spatial_reference=fc_sr) as cursor: # loop through features ...     for row in cursor: ...         cur_wkid = int(wkid_shift + row[1]) # calculate new wkid for the feature ...         proj_point = row[0].projectAs(arcpy.SpatialReference(cur_wkid)).centroid # project a new point ...         row[2] = proj_point.X # get the x value of the new point ...         row[3] = proj_point.Y # get the y value of the new point ...         cursor.updateRow(row) # write the new values to the feature‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

There's no single tool that will do this for you, as far as I know. You generally need additional information (multispectral imagery, terrestrial ecosystem/forestry mapping, etc.) to build a model, incorporating things you know about what you want to model (height, density, smoothness, reflective characteristics, etc.). ... View more

Hmmm... scroll down? ... View more

Rather than hard-coding the name of the OID/FID (sort_field = 'OBJECTID'), you can get it dynamically by describing the data: sort_field = arcpy.Describe(in_points).OIDFieldName‍ ... View more

I don't exactly know what you're trying to do. You're the one who needs to understand the logic of your script. Split() splits a string by a delimiter, which I think is the basic step you're missing, but you need to play with it, understand what's happening, and then change the logic of your script to achieve the final goal. Come up with your new script and if it doesn't do what you think it should, either post and ask for more help, or do some debugging to figure it out. ... View more

I was just showing how you can break apart a given text string into parts based on a delimiter. You can still use GetParameterAsText() to collect the user input (e.g. "x,y,z"), break it apart using split() (['x','y','z']), and make your replacements using those individual pieces. ... View more

On each help page related to arcpy (e.g. arcpy functions, classes, or geoprocessing tools) the Python syntax for the "thing" is near the bottom. In this case, it says: ExportToPNG (map_document, out_png, {data_frame}, {df_export_width}, {df_export_height}, {resolution}, {world_file}, {color_mode}, {background_color}, {transparent_color}, {interlaced}) I don't want to be too cryptic, but you need to do some of the work here. ... View more

There is a resolution parameter: ExportToPNG—Help | ArcGIS Desktop  ... View more

Yes, this goes in Raster Calculator, although I think you want the following expression that will round to the nearest 0.5 (Int truncates, it doesn't round): Int("aspect"*2+0.5)/2.0‍‍ ... View more

Maybe need one more line where Existing CRS = None (but data is already in Orange), then Orange, then Define, then Orange. ... View more