Latest Contributions by JohannesLindner

The error is in this line: var fromPointGeometry = g.paths[0][0];   What this does is take the first point of the first segment of a line feature. It works in the first expression, because you run it on a line feature. It doesn't work in the second expression, because you run it on a point feature and those don't have the path variable. This should work: var g = Geometry($feature) var fsLine = FeatureSetByName($datastore, "ssGravityMain", ["UPELEV"], true) var fromLine = First(Intersects(fsLine , g)) if (fromLine == null) { return } return $feature.RIMELEV - fromLine.UPELEV ... View more

# same goes for python, except you can use and/or/not TrapType = !TrapType! return 'Possum' if TrapType == 'Trapinator' or TrapType == 'Timms' or TrapType == 'Set & forget' else 'Rat' # easier: return 'Possum' if TrapType in ['Trapinator', 'Timms', 'Set & forget'] else 'Rat' ... View more

// test for equal with == // Aracade doesn't have and/or/not, it uses &&/||/! // test for equality has to be repeated var TrapType = $feature.Trap_type; IIf(TrapType == "Trapinator" || TrapType == "Timms" || TrapType == "Set & forget" , 'Possum', 'Rat'); // easier: IIF(Includes(['Trapinator', 'Timms', 'Set & forget'], TrapType), 'Possum', 'Rat'); ... View more

try: from tkinter import filedialog except ImportError: import tkFileDialog as filedialog import csv def interpolate(xy, x): """Linear interpolation xy: list of given data points as [(x, y)] x: x-value of unknown data point Returns the y-value as a linear interpolation between the given data points """ xy = sorted(xy, key=lambda p: p[0]) # sort list by x x0, y0 = [p for p in xy if p[0] <= x][-1] # get data point with next smaller x x1, y1 = [p for p in xy if p[0] >= x][0] # get data point with next larger x y = y0 + (y1 - y0) / (x1 - x0) * (x - x0) # linear interpolation return y def interpolate_empty_values(xy): """Fills empty x or y values by linear interpolation using the known data points. xy: list of data points as [(x, y)], x OR y can be null Returns a list with interpolated values """ known_xy = [(x, y) for x, y in xy if not (x is None or y is None)] known_yx = [(y, x) for x, y in known_xy] filled_xy = [] for x, y in xy: if not (x is None or y is None) or (x is None and y is None): filled_xy.append((x, y)) elif x is None: # Be very careful when interpolating x! # Interpolating x for things like relief plots, time series, # quadratic functions, etc. will give unintended results. x_ = interpolate(known_yx, y) filled_xy.append((x_, y)) elif y is None: y_ = interpolate(known_xy, x) filled_xy.append((x, y_)) return filled_xy # load your csv file csv_file = filedialog.askopenfile(title="Select csv file") # generate [(x, y)] csv_data = [line for line in csv.reader(csv_file)] csv_file.close() header = csv_data[0] csv_data = csv_data[1:] xy = [] for d in csv_data: x = None if d[0] in ['', ' '] else float(d[0]) y = None if d[1] in ['', ' '] else float(d[1]) xy.append((x, y)) # run interpolate_empty_values new_xy = interpolate_empty_values(xy) # write the result into a new csv file new_csv_file = filedialog.asksaveasfile(title="Save csv file") writer = csv.writer(new_csv_file) writer.writerow(header) for xy in new_xy: writer.writerow(xy) new_csv_file.close()   Paste that code into a script and run the script. You don't need arcpy or ArcGIS at all. In fact, expect errors or crashes when using the tkinter methods (filedialog.*) in the ArcGIS python window. ... View more

Putting code in a specially formatted box makes it much more readable. Expand the toolbar in the comment window, then click on "Inserd/Edit code sample":   Arcade has no way to know that you defined the variables in Python.   # If you want to define the variables in the Arcade code, do this: expression = """ var a = 1.51 var b = 1.48 var c = 5.03 when($feature.gridcode == 1, a, $feature.gridcode == 2, b, $feature.gridcode == 3, c, 0) """ arcpy.CalculateField_management(Image_classified, "Dens", expression, "ARCADE") # If you want to define the variable in Python, do this: a = 1.51 b = 1.48 c = 5.03 expression = "when($feature.gridcode == 1, {}, $feature.gridcode == 2, {}, $feature.gridcode == 3, {}, 0)".format(a, b, c) arcpy.CalculateField_management(Image_classified, "Dens", expression, "ARCADE")   ... View more