Latest Contributions by XanderBakker

There is a nice post by Katie Cullen‌: How to update the Max Record Count of a Feature Service in REST that you might want to look at. Kind regards, Xander ... View more

Maybe you find the necessary information in this website: Python Example - Google Cloud Storage — Google Cloud Platform Kind regards, Xander ... View more

Hi James, Great to hear that you can advance with the previous explanation, but to complete it, I'll continue with... 4) In case you have done a zonal statistics for the semi circle, you will have a count columns which holds the number of cells inside each semi circle. With your "histogram" (frequency distribution), you can sort the data on height (ascending order) and calculate the percentage for each record (remember: area = count * cellsize^2). Next calculate the cumulative area (or percentage) and determine the value that corresponds to the percentage you're looking for (e.g. the 5% you were mentioning). How to do this with some python code, please refer to this thread: https://community.esri.com/message/391944#392135 5) for determining the flood plain you should not take the amount of pixels that are beneath the flood level. Only those pixels that can be reached should be taken into account. An example thread were this is discussed you can find here: Re: Determine Floodplain: Based on known flood level Good luck with your analysis. Kind regards, Xander ... View more

What you can do is to create a Raster Catalog (in Data Management Tools\Raster\Raster Catalog), then use the Workspace to Raster Catalog tool, and after that the raster catalog can be exported (export footprints, in context sensitive menu of the raster catalog in your TOC). This will create a polygon featureclass of the extents. Kind regards, Xander ... View more

ArcGIS Team Python posted a nice example with python to add a rank field (but don't use the category field): Ranking field values | ArcPy Café This post also shows a python solution for your problem: https://community.esri.com/message/128653#128653 Kind regards, Xander ... View more

You could use some Python code to retrieve the JSON from Google Docs and join it using python (arcpy) to a featureclass. And yes, you should probably move this to a different place, since installation support has nothing to do with your question. Kind regards, Xander ... View more

If you are the administrator of the service, you will probably have access to the data. In that case you could use something like the ArcGIS Open Data Help | ArcGIS which has feature to allow users to download the data. In case you are not the owner of the data, you can use the REST API to query the data. A simple way to do this is using some Python code. A small example can be found here: https://community.esri.com/message/436087#436158 Please note the services by default are configured to return no more than 1000 features. If you are not the administrator of the service, you will have to do multiple requests on the REST service to obtain all features. You can provide the REST query with a geometry to restrict the results to that geometry. Kind regards, Xander ... View more

I suppose you would have to create a line from the GPS locations using the Points To Line tool (in Data Management Tools/Features). Next you would probably have to create a route (Linear Referencing Tools/Create Routes) using the time as measurements, followed by Make Route Event Layer (also in Linear Referencing Tools) to convert your samples to locations. Your samples should have the same unit of measurement (time, in seconds or minutes?) to make it work. I myself, would probably use some Python coding to make it work... Kind regards, Xander ... View more

The raster cells are always oriented horizontally and vertically (since they are always rows and columns). The may display rotated if the projection of your data frame is different from the projection of your raster. In ArcMap you can create different formats of data. An Esri grid is normally the result when writing to a folder. A geodatabase raster is the result if you write to a (file)geodatabase. If you use the default geodatabase as output location you will be created rasters inside a file geodatabase. If the focal statistics does not work for you, you may consider creating polygons for the neighborhoods you want to use and do a zonal statistics instead. Kind regards, Xander ... View more

Added an example on working with domains. Pieter Geert van den Beukel‌, next time it might be better to post your question as a thread. This way you can get more valuable answers which are more directed to your specific needs. Kind regards, Xander ... View more

Hi Pierter Geert, I can include some code sample to work with domains in Python, but it would require the use of the data access module. Since you stated "row.getValue()", you may be working with 10.0, the version before the da module was released. Do you have access to the da module? If so, my recomendation would be to switch to the da cursors, which are much faster. I will add some code sample to retrieve the domain description rather than the code. Kind regards, Xander ... View more

While looking for some problems I've encountered I came across this thread. Not sure if anyone is still "listening" to this thread, but thought I would post some code I have to list subtypes and correnponding default field values (and releted domains in the subtype). Maybe it is useful for someone... def main():     # import arceditor     import arcpy     import os     lst_amb = ["DLLO 9.3.1", "TEST 9.3.1", "PROD 9.3.1",                "DLLO 10.1", "TEST 10.1", "PROD 10.1"]     # dictionary with connection files     dct_conn = {"DLLO 9.3.1":r"C:\Users\xbakker\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\Desarrollo 9.3.1 (GEOGENESIS).sde",                 "TEST 9.3.1": r"C:\Users\xbakker\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\Pruebas 9.3.1 (GEOGENESIS).sde",                 "PROD 9.3.1": r"C:\Users\xbakker\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\Produccion 9.3.1 (GEOGENESIS).sde",                 "DLLO 10.1": r"C:\Users\xbakker\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\Desarrollo 10.1.sde",                 "TEST 10.1": r"C:\Users\xbakker\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\Pruebas 10.1.sde",                 "PROD 10.1": r"C:\Users\xbakker\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\Produccion 10.1.sde"}     # set encoding to utf-8     import sys     reload(sys)     print sys.getdefaultencoding()     def_enc = sys.getdefaultencoding()     sys.setdefaultencoding('utf8')     # format of TXT report with subtypes     report = r"D:\Xander\Genesis\Domains\report_subtypes_v4.txt"     with open(report, "w") as f:         f.write("Database\tDataset\tFeatureClass\tSubtypeField\tSubtypeCode\tSubtypeName\tIsDefaultValue\tFieldName\tFieldDefaultValue\tFieldDomain\n")         # loop through databases         for amb in lst_amb:             ws = dct_conn[amb]             arcpy.env.workspace = ws             lst_fds = arcpy.ListDatasets()             # loop through feature datasets             for fds in lst_fds:                 fcs = arcpy.ListFeatureClasses(feature_dataset=fds)                 if len(fcs) != 0:                     # loop through each featureclass en current feature dataset                     for fc in fcs:                         try:                             # list subtypes                             subtypes = arcpy.da.ListSubtypes(os.path.join(ws,fds,fc))                             # colect info on subtype                             for stcode, stdict in subtypes.iteritems():                                 if 'Default' in stdict:                                     IsDefaultValue = stdict['Default']                                 else:                                     IsDefaultValue = ""                                 if 'Name' in stdict:                                     SubtypeName = stdict['Name']                                 else:                                     SubtypeName = ""                                 if 'SubtypeField' in stdict:                                     SubtypeField = stdict['SubtypeField']                                     if SubtypeField != None and SubtypeField != "":                                         # list default field value and field details (like related domain)                                         if 'FieldValues' in stdict:                                             fields = stdict['FieldValues']                                             for fieldname, fieldvals in fields.iteritems():                                                 FieldDefaultValue = fieldvals[0]                                                 if not fieldvals[1] is None:                                                     domainname = fieldvals[1].name                                                 else:                                                     domainname = ""                                                 f.write('{0}\t{1}\t{2}\t{3}\t{4}\t{5}\t{6}\t{7}\t{8}\t{9}\n'.format(                                                          amb, fds, fc, SubtypeField, stcode, SubtypeName, IsDefaultValue,                                                          fieldname, FieldDefaultValue, domainname))                         except:                             f.write('{0}\t{1}\t{2}\t{3}\t{4}\t{5}\t{6}\t{7}\t{8}\t{9}\n'.format(                                     amb, fds, fc, "Error", -1, "Error", "Error",                                     "Error", -1, "Error"))     # restore default settings     sys.setdefaultencoding(def_enc) if __name__ == '__main__':     main() Kind regards, Xander ... View more