topic Re: ERROR 000875: I cannot get my raster into the dedicated output workspace while using ZIP- loop in Python Questions

ERROR 000875: I cannot get my raster into the dedicated output workspace while using ZIP- loop

IngoWolff — Thu, 23 Mar 2023 12:57:42 GMT

Hello all,

I want to do a raster calculation as done here: https://gis.stackexchange.com/questions/332925/looping-through-two-raster-folders-to-perform-raster-calculation and here: https://gis.stackexchange.com/questions/333082/looping-through-two-raster-folders-to-perform-raster-calculation-gives-error-000

It is more or less running correctly, but the output is written into the first env.workspace or second env.workspace. I simply cannot get the output rasters into my dedicated output directory.

Any clues?

Here my code snipped:

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")

#-----------------------------

# 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, but permits folder2 to get to work...

# folder path2

folder2 = r"K:\6202_Klimawandel_Folgen_Anpassung\620203_Werkbank\mGROWA\WHM_Klimaatlas_Neuentwicklungen\LPWCfinal.gdb"+ "\\" # not working!!!

#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))

Cheers for any help!

Ingo

Re: ERROR 000875: I cannot get my raster into the dedicated output workspace while using ZIP- loop

DanPatterson — Thu, 23 Mar 2023 13:31:38 GMT

indentation and line numbers would help

Code formatting ... the Community Version - Esri Community

Re: ERROR 000875: I cannot get my raster into the dedicated output workspace while using ZIP- loop

IngoWolff — Thu, 23 Mar 2023 13:38:09 GMT

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))

Re: ERROR 000875: I cannot get my raster into the dedicated output workspace while using ZIP- loop

DanPatterson — Thu, 23 Mar 2023 14:12:12 GMT

between line 41 and 42, try throwing in a print statement to print to make sure the preferred outputs are correct

Re: ERROR 000875: I cannot get my raster into the dedicated output workspace while using ZIP- loop

IngoWolff — Thu, 23 Mar 2023 14:22:21 GMT

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

Re: ERROR 000875: I cannot get my raster into the dedicated output workspace while using ZIP- loop

DavidPike — Thu, 23 Mar 2023 17:28:37 GMT

Can you share the output of your prints? I think it may be something to do with the raw strings then how you're concatenating them, it misreads the string and puts it into the workspace, I much prefer to use os.path.join

perhaps try https://pro.arcgis.com/en/pro-app/latest/arcpy/functions/clearenvironment.htm and see if it fails with a helpful message when trying to write to that location.

Re: ERROR 000875: I cannot get my raster into the dedicated output workspace while using ZIP- loop

IngoWolff — Fri, 24 Mar 2023 14:39:39 GMT

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

the output-wokspace is not valid (I have no clue why ths error message comes in German...)

Cheers

Ingo

Re: ERROR 000875: I cannot get my raster into the dedicated output workspace while using ZIP- loop

DavidPike — Fri, 24 Mar 2023 14:45:28 GMT

It's in German as its more efficient.

I think you just want to clear the workspace environment like the sample, try it exactly as below:

arcpy.ClearEnvironment("workspace") # prints None print(arcpy.env.workspace)

Re: ERROR 000875: I cannot get my raster into the dedicated output workspace while using ZIP- loop

DanPatterson — Fri, 24 Mar 2023 15:51:32 GMT

Is there any way you can try the process in a simpler folder without spaces and not beginning with numbers?

Re: ERROR 000875: I cannot get my raster into the dedicated output workspace while using ZIP- loop

IngoWolff — Mon, 27 Mar 2023 07:38:10 GMT

Ok, cheers for the advice with

arcpy.ClearEnvironment("workspace")

It seemed not working though, until I just created a new output path! It sometimes is so odd, why some output .gdbs are easily accepted within ArcCatalog, but not in a stand alone script!

However, this is my script now working, despite long path names which I need for my working routine...

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... arcpy.ClearEnvironment("workspace") # prints None print(arcpy.env.workspace) # 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 arcpy.ClearEnvironment("workspace") # prints None print(arcpy.env.workspace) # folder path folder2 = r"K:\6202_Klimawandel_Folgen_Anpassung\620203_Werkbank\mGROWA\WHM_Klimaatlas_Neuentwicklungen\LPWC_nfk.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" output_raster2.save(out_raster2) print("{} Processed Successfully!".format(out_raster2)) end_time = datetime.now() print('Prozessdauer: {}'.format(end_time - start_time))

Cheers!

Ingo

P.S. German and efficiency is like oil and water, seriously!