Latest Contributions by JohannesBierer

‎05-29-2020

Somehow I'm not able to add all values in a search cursor to a dictionary - one value pair is missing? What is my mistake? def Dict1(self): with arcpy.da.SearchCursor(outIntersect, ["FID_Anbinden1", "FLURSTUECKSTEXT"]) as cursor: for row in cursor: print("{} - {}".format(row[0], row[1])) myKey = row[1] myValue = row[0] myDict[myKey] = myValue # myDict.update({myKey: myValue}) print myDict‍‍‍‍‍‍‍‍‍‍‍ This leads to this output (as well as .update in line 9): 1 - 10125 1 - 10079 1 - 10337/2 2 - 10125 {u'10079': 1, u'10337/2': 1, u'10125': 2} u'10125':1 is missing?

‎05-08-2020

Thank you David, didn't see your answer. But it works with one "/", I only made a big mistake in line 6 🙂

‎05-08-2020

This worked for me: sc = line.replace("[", "") # 082640-000-07332/001 convert to 0826400000733200100 sc2 = sc.replace("]", "") sc21 = sc2.replace("-", "") sc22 = sc21.replace("/", "") sc23 = "{}{}".format(sc22, str("00")) sc24 = sc23.replace("'", "") sc25 = sc24.encode('ascii') sc26 = sc25.replace("u", "") sc3 = "{} = '{}'".format(field, sc26) print sc3‍‍‍‍‍‍‍‍‍‍‍‍‍

‎05-08-2020

Hello, I have an array looking like this 082640-000-07332/001 and would like to convert it to 0826400000733200100 for a search clause? I tried this: sc = line.replace("[", "") # 082640-000-07332/001 convert to 0826400000733200100 sc2 = sc.replace("]", "") sc21 = sc2.replace("-", "") sc22 = sc21.replace("/", "") sc23 = "{}{}".format(sc21, "00") sc3 = "{} = '{}'".format(field, sc23) print sc3‍‍‍‍‍‍‍‍‍ but this code results in this (print sc3): FLURSTUECKSKENNZEICHEN = 'u'08264000007138/000'00' Any help would be great?

‎05-04-2020

Yes like a perimeter/area ratio. I would like to differentiate between compact features I want to preserve and slivers I want to delete after the use of a geoprocessing tool?

‎05-04-2020

I'm searching for a possibility to calculate landscape metric "shape index" in the field calculator? Shape index = 1 when the patch is circular and increases without limit as the patch becomes more irregular.

‎04-08-2020

Hi Norman, Thank you for your answer and the hints concerning the arcpy sa module. Yesterday I used a dem with 5m resolution. Now I get some kind of relatively hight resolution laserdata with 8 points/m². I calculated a dem raster with a resolution of 0.3 m out of this data. The only reason I use fill infront of the resampling is to smoothe the dem slightly in an appropriate manner without changing the data too much. Especially the local hights I don't want to change. With this resolution I get a very rough surface so the fill tool is in my opinion the better solution to smooth the dem than resampling the data at this step? After the fill step I resampled the data to a cellsize of 3m bilinear. Would you use nearest neighbor in this step? After that I go on with flow direction ... sink ... To your questions: 1. I'm not interested in all (small) sinks, I need sinks larger than an area (let's say about 5000m²) which are worth to be considered as rewetting area. 2. Yes I would like to get the total areas of all available sinks that meet the certain > 5000m². 3. The only purpose of changing the resolution is that I get otherwise only millions of very small sinks because of the "hight" resolution of my data. How would you arrange the steps of the calculation? Would you put in a fill step after the resampling to a cellsize of 3m?

‎04-08-2020

You are right, if I resample to strong I get totally different areas (which didn't existed before) after a unknown threshold. The fill tool does a good job to smooth the elevation data. Didn't know the Aggregate Tool, will try it. Thank you for your thoughts.

‎04-08-2020

I'm searching for areas for rewetting, so my intentions is to get some aggregated areas where it's worth to try build barriers to dam the water.

‎04-08-2020

Thank you for the link. I thought if I fill the dem first it would sum the sinks better up. Would you remove the fill calculation? I don't really understand why the small lakes in the area are not all thrown out as sinks? Even in the best result I could reach so far.

‎04-08-2020

Sink results in very different demarkations depending on the resampling of the raster data. Could someone explain why? I don't understand this? In the attached images you could see the differences. I use elevation data with a resolution of 5 meters. Sink06m ... is in my opinion better because nearly all lakes are demarked as sinks. Sink06m is resampled to a cellsize of 6m biliniear. Sink10m is resampled to a cellsize 10 meters, but the results are totally different? Am I doing something wrong? This is how I calculated the areas: import arcpy, os from arcpy.sa import * arcpy.env.overwriteOutput = True arcpy.env.workspace = r"MyPath\DGM.gdb" # Local variables: DGM005_tif = r"MyPath\DGM005_3.tif" Fill_tif1 = "Fill_1" resample = "resample" FlowDir_Fill2 = "FlowDir_1" Sink_FlowDir1 = "Sink_FlowDir1" absenkung = "Absenkung" arcpy.CheckOutExtension("Spatial") outCon2 = Con(Raster(DGM005_tif) > 0, DGM005_tif) # Process: Füllung arcpy.gp.Fill_sa(outCon2, Fill_tif1, "") print("fill") arcpy.Resample_management(Fill_tif1, resample, "8", "BILINEAR") print("resample") # Process: Fließrichtung arcpy.gp.FlowDirection_sa(resample, FlowDir_Fill2, "NORMAL", absenkung) print("flowdir") # Process: Senke arcpy.gp.Sink_sa(FlowDir_Fill2, Sink_FlowDir1) print("sink") arcpy.RasterToPolygon_conversion(Sink_FlowDir1, "{}{}".format("Senke08mfill", os.path.basename(os.path.splitext(DGM005_tif)[0])), "SIMPLIFY", "VALUE") print("Done") arcpy.CheckInExtension("Spatial")‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

‎04-06-2020

Don't know if I'm up to date but I think you have to use make feature layer infront of your select by location. Out of the help files: import arcpy arcpy.env.workspace = "c:/data/mexico.gdb" # Make a layer and select cities which overlap the chihuahua polygon arcpy.MakeFeatureLayer_management('cities', 'cities_lyr') arcpy.SelectLayerByLocation_management('cities_lyr', 'intersect', 'chihuahua')‍‍‍‍‍ you should use: arcpy.env.workspace = r"C:\folder\Data.gdb" and not arcpy.env.workspace = r"C:\folder\Data.gdb\" r"C:\Data.gdb\%Name%_selection" will not work in python code, you will need to create a unique name for your output file.

‎03-24-2020

It's not really clear what you want to do? Feature Class to Excel you can use table to excel tool?

‎03-24-2020

Could this help you? arcgis pro - arcpy.ListFields() displaying numerical field names instead of actual field names - Geographic Information …

‎02-25-2020

Open a map might be possible with os.startfile() ?

Online Status	Offline
Date Last Visited	Sunday

My Ideas

Latest Contributions by JohannesBierer

Missing values in a dictionary?

Re: Replace Stringparts

Re: Replace Stringparts

Replace Stringparts

Re: How to calculate landscape metric "shape index" in the field calculator?

How to calculate landscape metric "shape index" in the field calculator?

Re: Sink results in very different demarkations?

Re: Sink results in very different demarkations?

Re: Sink results in very different demarkations?

Re: Sink results in very different demarkations?

Sink results in very different demarkations?

Re: Arcpy - Select by location - iterate through gdb feature class and export to new fc

Re: arcpy.ListFields

Re: arcpy.ListFields

Re: Is it possible to open a Map or Layout in ArcGIS Pro using Python?

Re: Alternative way to save Feature Set to Shapefi...

Re: Error processing drone imagery over river

Re: How create a line that go through to point and...

Re: What's the future of Concurrent licenses?

Re: What's the future of Concurrent licenses?

Anwendergruppe / UserGroup NRW