Latest Contributions by IngoWolff

Hi David, cheers for tipping me to arcpy.ClearEnvironment("my_workspace")   However, it es not working.   The error message is as follows:   Running script: "K:\6202_Klimawandel_Folgen_Anpassung\620203_Werkbank\Zentrale_Scripte\NFK_mGROWA.py" Traceback (most recent call last): File "K:\6202_Klimawandel_Folgen_Anpassung\620203_Werkbank\Zentrale_Scripte\NFK_mGROWA.py", line 163, in <module> arcpy.ClearEnvironment(LPAWC) File "C:\Program Files\ESRI\ArcGISPro\Resources\ArcPy\arcpy\__init__.py", line 275, in ClearEnvironment return gp.clearEnvironment(environment_name) File "C:\Program Files\ESRI\ArcGISPro\Resources\ArcPy\arcpy\geoprocessing\_base.py", line 126, in clearEnvironment self._gp.ClearEnvironment(*gp_fixargs(args, True))) AttributeError: Object: Cannot clear environment <K:\6202_Klimawandel_Folgen_Anpassung\620203_Werkbank\mGROWA\WHM_Klimaatlas_Neuentwicklungen\LPAWCDate.gdb>     I just added the command for bothe einvironment workspace below my output directory "folder2". When I comment both "ClearEnvironment" commands out, the cript says   Running script: "K:\6202_Klimawandel_Folgen_Anpassung\620203_Werkbank\Zentrale_Scripte\NFK_mGROWA.py" Traceback (most recent call last): File "K:\6202_Klimawandel_Folgen_Anpassung\620203_Werkbank\Zentrale_Scripte\NFK_mGROWA.py", line 172, in <module> output_raster2.save(out_raster2) RuntimeError: ERROR 000875: Output raster: Der Workspace von K:\6202_Klimawandel_Folgen_Anpassung\620203_Werkbank\mGROWA\WHM_Klimaatlas_Neuentwicklungen\LPWCfinal.gdb\L1pawc20230222_prznFK ist ein ungültiger Ausgabe-Workspace.    the output-wokspace is not valid (I have no clue why ths error message comes in German...)   Cheers Ingo ... View more

Hi Dan, ok, the print statement gives the correct output workspace ans anticipated rastername...but only if i comment the .save - command from line 42 (then 43) out.   Cheers! Ingo   ... View more

Apologies! I'm not too often posting scripts... import arcpy from arcpy import env from arcpy.sa import *#Spatial Analyst module import os from datetime import datetime import subprocess import glob arcpy.env.overwriteOutput = True arcpy.env.matchMultidimensionalVariable = False arcpy.CheckOutExtension("Spatial") start_time = datetime.now() # Divide LPAWC*100 by LAPWC # get first set of rasters LPAWC = r"K:\6202_Klimawandel_Folgen_Anpassung\620203_Werkbank\mGROWA\WHM_Klimaatlas_Neuentwicklungen\LPAWCDate.gdb" #'source 1' arcpy.env.workspace = LPAWC rasters_LPAWC = arcpy.ListRasters() rasters_LPAWC = [os.path.join(LPAWC, r) for r in rasters_LPAWC] # necessary for zip, but permits folder2 to get to work... # get second set of rasters LPAWFC = r"K:\6202_Klimawandel_Folgen_Anpassung\620203_Werkbank\mGROWA\WHM_Klimaatlas_Neuentwicklungen\LPAWFC.gdb" #'source 2' arcpy.env.workspace = LPAWFC rasters_LPAWFC = arcpy.ListRasters() rasters_LPAWFC = [os.path.join(LPAWFC, rs) for rs in rasters_LPAWFC] # necessary for zip # folder path2 folder2 = r"K:\6202_Klimawandel_Folgen_Anpassung\620203_Werkbank\mGROWA\WHM_Klimaatlas_Neuentwicklungen\LPWCfinal.gdb"+ "\\" #loop for multiply by 100 and divide by LPAWC to get percent nFK for LPAWC_ras, LPAWFC_ras in zip(rasters_LPAWC, rasters_LPAWFC): r1 = Raster(LPAWC_ras) r2 = Raster(LPAWFC_ras) output_raster2 = 100*r1/r2 out_raster2 = folder2 + arcpy.Describe(r1).baseName + "_prznFK" #'{}{}'.format(r1,"_prznFK") output_raster2.save(out_raster2) print("{} Processed Successfully!".format(out_raster2)) end_time = datetime.now() print('Prozessdauer: {}'.format(end_time - start_time)) ... View more

Hi Joshua, the scripts jumps first into the workspace directory with dir = r"F:\Daten_Wolff\Topographische_Gefaehrdungsanalyse\Test_Source" arcpy.env.workspace = dir‍‍ Then, inside the workspace directory, the list of gdbs is called by gdbList = arcpy.ListWorkspaces('*','FileGDB') for gdb in gdbList: arcpy.env.workspace = gdb #--change working directory to each GDB in list‍‍‍ Now it goes from gdb to gdb. For each gdb, I then apply the classical list of Raster loop, not befor telling the target directory as "folder": folder = r"F:\Daten_Wolff\Topographische_Gefaehrdungsanalyse\BlueSpot_50.gdb"+ "\\" RasterList = arcpy.ListRasters('*_050') for Raster in RasterList: print("processing raster: %s" %os.path.join("env.workspace",Raster)) outRasterdataset = folder + arcpy.Describe(Raster).baseName arcpy.CopyRaster_management (Raster, outRasterdataset, "DEFAULTS", "", "-999", "NONE", "NONE" , "32_BIT_FLOAT", "NONE", "NONE" )‍‍‍‍‍‍‍‍‍‍‍‍ This works in ArcPy 2.7.14 quite nicely so far. I have exatly 5107 gdbs to be searched für fgdbr Raster with the name suffix '*_050' or '*_100' and so on. Not in all of thes gdbs are thes Raster with that exact suffices, but the script shecks them all! So far, from my 5107 gdbs the script copied nicely 4936 fgdbrs with the suffix '*_050' into my targed gdb. I directly added the same loops for more suffices, and its working calm and steadily on my VM. All the best! Ingo ... View more

Ok, it appears I found a solution! Just run everything with ArcPy 2.7.xx instead of ArcPy3.6.xx... for the first time it occured to me, that ArcPy 3.6.xx cannot solve a task ArcPy 2.7.xx can solve. Here the working script, tested on a good dozen of gdbs: import arcpy from arcpy import env import os arcpy.env.overwriteOutput = True arcpy.env.parallelProcessingFactor = "75%" # # #-------------------- #Setting up first directories for source and destination for iterative raster data loop no 1 dir = r"F:\Daten_Wolff\Topographische_Gefaehrdungsanalyse\Test_Source" arcpy.env.workspace = dir gdbList = arcpy.ListWorkspaces('*','FileGDB') for gdb in gdbList: arcpy.env.workspace = gdb #--change working directory to each GDB in list folder = r"F:\Daten_Wolff\Topographische_Gefaehrdungsanalyse\BlueSpot_50.gdb"+ "\\" RasterList = arcpy.ListRasters('*_100') for Raster in RasterList: print("processing raster: %s" %os.path.join("env.workspace",Raster)) outRasterdataset = folder + arcpy.Describe(Raster).baseName arcpy.CopyRaster_management (Raster, outRasterdataset, "DEFAULTS", "", "-999", "NONE", "NONE" , "32_BIT_FLOAT", "NONE", "NONE" ) ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ This might explain, why the error message beforehand told me to go to Esri Support. Possibly a bug in ArcGISPro. However, many thanks for your kind replies! Ingo ... View more

Hi! I got this error-message: arcgisscripting.ExecuteError: Failed to execute. Parameters are not valid. ERROR 000732: Input Raster: Dataset bluespot_T2718199_050 does not exist or is not supported Failed to execute (CopyRaster). I think it is due to the fgdbr Format. But this is rather strange, as it is an ESRI dataformat which should easily be readable in arcpy... ... View more

Exactly, the rasters are file geodatabase raster datasets (fgdbr)... I manage to get all fgdbrs properly displayed with "print", but with the list approach I can't write the fgdbrs into the correct targed directory and withe the walk-approach I can't read the fgdbrs properly ... Would be lovely to get some official help from ‌ as the error code suggests. Cheers Ingo ... View more

What exactly do you mean with another name? Type of raster or just a different part of filename for the selection? As far as I found out, FGDBR Raster are not mentioned within the list of types in Walk—Help | Dokumentation ... Therefore I'm afraid, "Walk" is not really the thing I could use. When I use the List approach it fails to write the rasters into the targed gdb instead of failing of reading the selected rasters...quite an hard nut to crush ... View more

Hi, ok, I tried to use walk instead of the more complicated list-approach. It now seems to go more exactly through the directories, however, I cannot get my rasters within the gdbs selected and copied into a new gdb. import arcpy import os workspace = r"F:\Daten_Wolff\Topographische_Gefaehrdungsanalyse\Test_Source" folder = r"F:\Daten_Wolff\Topographische_Gefaehrdungsanalyse\BlueSpot_50.gdb" + "\\" rasters = [] for dirpath, dirnames, filenames in arcpy.da.Walk(workspace, topdown=True, datatype="RasterDataset"): for filename in filenames: #print("processing Workspace: %s" %os.path.join(dirpath,filename)) if "_050" in filename: print("processing Workspace: %s" %os.path.join(dirpath,filename)) outRasterdataset = folder + str(filename.rsplit(".",0)[0]) arcpy.CopyRaster_management (filename, outRasterdataset, "DEFAULTS", "", "-999", "NONE", "NONE" , "32_BIT_FLOAT", "NONE", "NONE" ) ‍‍‍‍‍‍‍‍‍‍‍‍‍ I only get following error message: Python 3.6.6 |Anaconda, Inc.| (default, Jun 28 2018, 11:27:44) [MSC v.1900 64 bit (AMD64)] on win32 Type "copyright", "credits" or "license()" for more information. >>> RESTART: F:\Daten_Wolff\Scripte\Zentrale_Scripte\Topostar_Copy_raster_5_walk.py processing Workspace: F:\Daten_Wolff\Topographische_Gefaehrdungsanalyse\Test_Source\T2718199.gdb\bluespot_T2718199_050 Traceback (most recent call last): File "F:\Daten_Wolff\Scripte\Zentrale_Scripte\Topostar_Copy_raster_5_walk.py", line 12, in <module> arcpy.CopyRaster_management (filename, outRasterdataset, "DEFAULTS", "", "-999", "NONE", "NONE" , "32_BIT_FLOAT", "NONE", "NONE" ) File "C:\Program Files\ESRI\ArcGISPro\Resources\ArcPy\arcpy\management.py", line 15289, in CopyRaster raise e File "C:\Program Files\ESRI\ArcGISPro\Resources\ArcPy\arcpy\management.py", line 15286, in CopyRaster retval = convertArcObjectToPythonObject(gp.CopyRaster_management(*gp_fixargs((in_raster, out_rasterdataset, config_keyword, background_value, nodata_value, onebit_to_eightbit, colormap_to_RGB, pixel_type, scale_pixel_value, RGB_to_Colormap, format, transform), True))) File "C:\Program Files\ESRI\ArcGISPro\Resources\ArcPy\arcpy\geoprocessing\_base.py", line 498, in <lambda> return lambda *args: val(*gp_fixargs(args, True)) arcgisscripting.ExecuteError: Failed to execute. Parameters are not valid. ERROR 000732: Input Raster: Dataset bluespot_T2718199_050 does not exist or is not supported Failed to execute (CopyRaster). >>> ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ So it seems, that my raster data format FileGeoDatabaseRaster (FGDBR) is not supported? At least it finds the right raster I wanted... What am I missing? Many thanks! Ingo ... View more

Ok, after some iterations on my script I found out, that the gdbList = arcpy.ListWorkspaces('*','FileGDB') for gdb in gdbList: arcpy.env.workspace = gdb #--change working directory to each GDB in list datasetList = arcpy.ListDatasets('*','Raster') #--make a list of all (if any) rastere datasets that exist in current GDB for dataset in datasetList: arcpy.env.workspace = dataset ‍‍‍‍‍‍ part is useless because I just want the FGDBR within each file geodatabase... I shortened my loops accordingly: import arcpy from arcpy import env import os arcpy.env.overwriteOutput = True arcpy.env.parallelProcessingFactor = "75%" # # #-------------------- #Setting up first directories for source and destination for iterative raster data loop no 1 dir = r"F:\Daten_Wolff\Topographische_Gefaehrdungsanalyse\Test_Source" arcpy.env.workspace = dir gdbList = arcpy.ListWorkspaces('*','FileGDB') for gdb in gdbList: print("processing Workspace: %s" %os.path.join("env.workspace",gdb)) arcpy.env.workspace = gdb #--change working directory to each GDB in list folder = r"F:\Daten_Wolff\Topographische_Gefaehrdungsanalyse\BlueSpot_50.gdb"+ "\\" RasterList = arcpy.ListRasters('*_050*') for Raster in RasterList: print("processing raster: %s" %os.path.join("env.workspace",Raster)) outRasterdataset = folder + str(Raster.rsplit("_",0)[0]) arcpy.CopyRaster_management (Raster, outRasterdataset, "DEFAULTS", "", "-999", "NONE", "NONE" , "32_BIT_FLOAT", "NONE", "NONE" ) ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ So technically my script should go through each file geodatabase and then check für FGDBRs ending with *_050 in the filename and copy them into my target file geodatabase named "folder". Sadly I get this error message now: Python 3.6.6 |Anaconda, Inc.| (default, Jun 28 2018, 11:27:44) [MSC v.1900 64 bit (AMD64)] on win32 Type "copyright", "credits" or "license()" for more information. >>> = RESTART: F:/Daten_Wolff/Scripte/Zentrale_Scripte/Topostar_Copy_raster_4.py = processing Workspace: F:\Daten_Wolff\Topographische_Gefaehrdungsanalyse\Test_Source\T2718199.gdb processing raster: env.workspace\bluespot_T2718199_050 Traceback (most recent call last): File "F:/Daten_Wolff/Scripte/Zentrale_Scripte/Topostar_Copy_raster_4.py", line 24, in <module> arcpy.CopyRaster_management (Raster, outRasterdataset, "DEFAULTS", "", "-999", "NONE", "NONE" , "32_BIT_FLOAT", "NONE", "NONE" ) File "C:\Program Files\ESRI\ArcGISPro\Resources\ArcPy\arcpy\management.py", line 15289, in CopyRaster raise e File "C:\Program Files\ESRI\ArcGISPro\Resources\ArcPy\arcpy\management.py", line 15286, in CopyRaster retval = convertArcObjectToPythonObject(gp.CopyRaster_management(*gp_fixargs((in_raster, out_rasterdataset, config_keyword, background_value, nodata_value, onebit_to_eightbit, colormap_to_RGB, pixel_type, scale_pixel_value, RGB_to_Colormap, format, transform), True))) File "C:\Program Files\ESRI\ArcGISPro\Resources\ArcPy\arcpy\geoprocessing\_base.py", line 498, in <lambda> return lambda *args: val(*gp_fixargs(args, True)) arcgisscripting.ExecuteError: ERROR 999999: Something unexpected caused the tool to fail. Contact Esri Technical Support (http://esriurl.com/support) to Report a Bug, and refer to the error help for potential solutions or workarounds. The table was not found. [VAT_bluespot_T2718199_050] No spatial reference exists. The table was not found. [bluespot_T2718199_050] Failed to execute (CopyRaster). >>> ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ So, where is my mistake? Any advice greatly appreciated! Cheers Ingo ... View more

Hi Dan, cheers for your quick answer! Regarding my indentation, I now follow this example arcgis 10.0 - Listing feature classes in multiple geodatabases in folder using ArcPy? - Geographic Information Systems S…  to go throgh my file geodatabases and datasets. My script does look like this: import arcpy from arcpy import env import os arcpy.env.overwriteOutput = True arcpy.env.parallelProcessingFactor = "75%" # # #-------------------- #Setting up first directories for source and destination for iterative raster data loop no 1 dir = r"F:\Daten_Wolff\Topographische_Gefaehrdungsanalyse\Test_Source" arcpy.env.workspace = dir folder = r"F:\Daten_Wolff\Topographische_Gefaehrdungsanalyse\BlueSpot_050.gdb"+ "\\" gdbList = arcpy.ListWorkspaces('*','FileGDB') for gdb in gdbList: arcpy.env.workspace = gdb #--change working directory to each GDB in list datasetList = arcpy.ListDatasets('*','Raster') #--make a list of all (if any) rastere datasets that exist in current GDB for dataset in datasetList: arcpy.env.workspace = dataset #--change working directory to each dataset (if any) in list for Raster in arcpy.ListRasters(): arcpy.env.workspace = dataset for Raster in arcpy.ListRasters('*050*'): print("processing raster: %s" %os.path.join("env.workspace",Raster)) outRasterdataset = folder + str(Raster.rsplit("_",0)[0]) arcpy.CopyRaster_management (Raster, outRasterdataset, "DEFAULTS", "", "-999", "NONE", "NONE" , "32_BIT_FLOAT", "NONE", "NONE" ) ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ There is no error message, but no result either after this script has run through... Thank you for your patience with me newbie! Ingo ... View more

Cheers for that advice, I changed my snipped accordingly ... View more