Latest Contributions by JohannesBierer

There's as well a identation error at line 14. I thought with line 15 I could do a for loop through each part of the geometrie? ... View more

Maybe it is helpful to create a point feature class of firstPoint == lastPoint places? You could try this script? import arcpy arcpy.env.overwriteOutput = True inFc = r"Yourpath...\Problem Data.gdb\problem_road_isolated" inFcS = r"Yourpath...\Problem Data.gdb\problem_road_isolated_single" outPoint = r"Yourpath...\Problem Data.gdb\\ClosedRing" arcpy.MultipartToSinglepart_management(inFc, inFcS) CoordList = list() with arcpy.da.SearchCursor (inFcS, "SHAPE@") as sCurs: for geom, in sCurs: for i in range(geom.partCount): part = geom.getPart(i) if geom.lastPoint.X == geom.firstPoint.X: CoordList.append("{}; {}".format(geom.lastPoint.X, geom.lastPoint.Y)) pt = arcpy.Point() ptGeoms = [] spatial_reference = arcpy.SpatialReference(3111) for p in CoordList: coord = p.split(";") pt.X = coord[0] print(coord[0]) pt.Y = coord[1] print(coord[1]) ptGeoms.append(arcpy.PointGeometry(pt, spatial_reference)) arcpy.CopyFeatures_management(ptGeoms, outPoint) arcpy.AddXY_management(outPoint)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

Could this be helpful? Value Error when adding adding a layer in ArcGIS Pro arcpy script - Geographic Information Systems Stack Exchange  ... View more

Sometimes a small hint helps, glad you get it to work. ... View more

I don't know how to do this in the model builder - but with python maybe you could define your input data and output with the following code? It's from: datetime - Iterating through a range of dates in Python - Stack Overflow  from datetime import timedelta, date def daterange(start_date, end_date): for n in range(int ((end_date - start_date).days)): yield start_date + timedelta(n) start_date = date(2009, 12, 31) end_date = date(2010, 12, 31) for single_date in daterange(start_date, end_date): d = single_date.strftime("%Y%m%d") d1 = single_date+timedelta(days=1) d2 = d1.strftime("%Y%m%d") input1 = "DSLR_CHIRPS_{}.tif".format(d) input2 = "CHIRPS_{}.tif".format(d2) output = "DSLR_CHIRPS_{}.tif".format(d2) print("{} - {} - {}".format(input1, input2, output))‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

Maybe this could be helpful as well - there are models to download? Building models for automatic landslide-susceptibility analysis, mapping and validation in ArcGIS | SpringerLink  ... View more

Don't forget to backup before ... This code can be used to rename your raster data to MOD60, MOD61. Hopefully. Without terra and aqua, don't know why you want this? import arcpy, os, re workspace = r"yourpath" walk = arcpy.da.Walk(workspace, topdown=True, datatype="RasterDataset") for dirpath, dirnames, filenames in walk: for filename in filenames: rasterE = os.path.join(dirpath, filename) filename1 = re.sub("04_3K.A20190", "", filename) filename2 = filename1[:5] print filename2 rasterRename = os.path.join(dirpath, filename2) arcpy.Rename_management(rasterE, rasterRename)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

After line 13 you could insert: prefix1 = prefix1.replace("O", "Y") Python String replace() Method - Tutorialspoint  ... View more

It can't work because you use different prefix. MOD60_ is not the same as MYD60_. So you could rename MYD60_ to MOD60_ etc. and it should find the rasterdata? Or use replace() in you script? ... View more

Could you post the code you use like this? https://community.esri.com/docs/DOC-8691-posting-code-with-syntax-highlighting-on-geonet  ... View more

Don't know if it could be the reason, you are using spaces in your paths: ...\final 2\... . What if you use ...\final_2\... ? ... View more

I ofter can't find self intersections with the Check Geometry Tool, it is a bit unreliable especially with very small mistakes. I can offer several shapes with self intersections not found by Check Geometry Tool (10.4.1) - unfortunately I don't have pro. ... View more

A working possibility could be to convert to numpy array, then use pandas to convert the array to excel? import arcpy import pandas as pd arcpy.env.overwriteOutput = True inTable = r"yourpath\Export_Output.shp" outTable1 = r"yourpath\export.xls" fieldNames = ["OBJECTID", "PARCEL_NO", "OwnerName", "Address", "City", "State", "ZipCode"] arr = arcpy.da.FeatureClassToNumPyArray(inTable, fieldNames) df = pd.DataFrame(arr) df1 = df.to_excel(outTable1, index=False)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more