topic Re: Run a python code in ArcGIS in Python Questions

Run a python code in ArcGIS

DuminduJayasekera — Wed, 14 Jan 2015 01:26:55 GMT

Hi,

I need to make it run the following code in arcGIS but this is not written to run in arcgis but in visual studio. I tried but it DOES NOT work. Can someone help me to write thin in python so that I can add as a script in ArcGIS.? I am planning to import this as a script after making this a working code.

Any help is highy appreciated.

Thanks.

------------------------------------------------------------------------

import arcpy

from arcpy import env

from arcpy.sa import *

#Check out the ArcGIS Spatial Analyst extension license

arcpy.CheckOutExtension("Spatial")

arcpy.env.overwriteOutput = True

arcpy.env.scratchWorkspace = "c:\temp\tmp"

env.workspace = "c:\DEMPreProcess"

LULC = Raster("C:\DEMPreProcess\VietnamLULC_Resample.tif")

Simard = Raster("C:\DEMPreProcess\MajoritySTOht.tif")

resultmap = Simard

for i in range(1,14):

if i == 1:

number = 0.4

elif i == 2:

number = 0.4

elif i == 3:

number = 0.4

elif i == 4:

number = 0.4

elif i == 5:

number = 0.4

elif i == 6:

number = 0.3

elif i == 7:

number = 0.3

elif i == 8:

number = 0.3

elif i == 9:

number = 0.3

elif i == 10:

number = 0.3

elif i == 11:

number = 0

elif i == 12:

number = 0.3

elif i == 13:

number = 0

elif i == 14:

number = 0.3

resultmap = Con(LULC == i,Simard*number,resultmap)

resultmap.save("C:\DEMPreProcess\MajoritySTOhtPer.tif")

Re: Run a python code in ArcGIS

DanPatterson_Retired — Wed, 14 Jan 2015 01:38:17 GMT

What is the error message? and what is it supposed to do?

Re: Run a python code in ArcGIS

DuminduJayasekera — Wed, 14 Jan 2015 02:09:09 GMT

Hi Dan,

I am new to python but what I want to write a script in python to take the integer values from i=1 to 14 in VietnamLULC_Resample.tif file and then multiplied by the interger values in the MajoritySTOht.tif file for each corresponding i.

I need to include this process as a pythin script to this model builder.

Thanks and help is appreciated.

Re: Run a python code in ArcGIS

DanPatterson_Retired — Wed, 14 Jan 2015 03:14:33 GMT

It is hard to tell from your code formatting, but I think your last two lines should be indented so that they are within the for loop...currently, they are outside the loop by all appearances. You should use the Advanced Editor, then select >> and choose Syntax Highlighting and format your code in the original thread to Python syntax to make it easier to read and so indentation can be confirmed.

Re: Run a python code in ArcGIS

XanderBakker — Wed, 14 Jan 2015 03:25:42 GMT

Why not reclassify the LULC raster (holding the 14 integer values) using ReclassByASCIIFile with a ASCII file that has the relation between the input and output values like this:

1 : 0.4

2 : 0.4

3 : 0.4

4 : 0.4

5 : 0.4

6 : 0.3

7 : 0.3

8 : 0.3

9 : 0.3

10 : 0.3

11 : 0

12 : 0.3

13 : 0

14 : 0.3

... and then multiply the reclassified raster by the Simard raster.

Re: Run a python code in ArcGIS

DuminduJayasekera — Wed, 14 Jan 2015 18:33:15 GMT

Dan,

The above code is run in Visual Studio but I want to convert it to pythin to include in the model bilder.

Where can I find the Advanced editor and syntax highlighting?

Thanks again for your reply.

Re: Run a python code in ArcGIS

XanderBakker — Wed, 14 Jan 2015 18:41:40 GMT

In the upper right corner of the editor (when editing a post) you will see the link "Use advanced editor". The code you supplied is Python code. It is possible that Visual Studio has been used to run it, but the IDE has no influence on the code itself.

Re: Run a python code in ArcGIS

DuminduJayasekera — Wed, 14 Jan 2015 19:02:53 GMT

OK. But How can I include this code into my model buider environment? and to make it run?

Thanks and highly appreciate your help. I am very new to python and GIS model builder.

Re: Run a python code in ArcGIS

JamesCrandall — Wed, 14 Jan 2015 19:10:22 GMT

The above code is run in Visual Studio

What project type in VS allows you run python code? Or are you executing a .bat file from the application you developed in Visual Studio (there's lots of different types of apps you can develop, VS is just the IDE).

Re: Run a python code in ArcGIS

JamesCrandall — Wed, 14 Jan 2015 19:12:13 GMT

Create a new geoprocessing tool, set the source to your .py file, add the tool to your model.

Re: Run a python code in ArcGIS

XanderBakker — Sun, 12 Dec 2021 07:43:35 GMT

If you can include a small subset of your data I can have a look at the calculation you intent to do, since doing 14 calculations for what can be done in one, doesn't sound correct to me.

You would have to adapt the code to allow for parameters: http://resources.arcgis.com/en/help/main/10.2/index.html#/Understanding_script_tool_parameters/00150000002s000000/

I suppose all things listed below should be change into parameters. So something like this:

env.workspace = arcpy.GetParameterAsText(0)
LULC = Raster(arcpy.GetParameterAsText(1))
Simard = Raster(arcpy.GetParameterAsText(2))

resultmap.save(arcpy.GetParameterAsText(3))

instead of this:

env.workspace = "c:\DEMPreProcess"
LULC = Raster("C:\DEMPreProcess\VietnamLULC_Resample.tif")
Simard = Raster("C:\DEMPreProcess\MajoritySTOht.tif")

resultmap.save("C:\DEMPreProcess\MajoritySTOhtPer.tif")

Re: Run a python code in ArcGIS

DuminduJayasekera — Wed, 14 Jan 2015 19:39:18 GMT

Hi James,

Could you please illustrate how to do that? I am totally new to this.

Thanks again and appreciate it.

Re: Run a python code in ArcGIS

JamesCrandall — Wed, 14 Jan 2015 19:43:20 GMT

http://resources.arcgis.com/en/help/main/10.1/index.html#//001500000006000000http://

Re: Run a python code in ArcGIS

DuminduJayasekera — Wed, 14 Jan 2015 20:03:23 GMT

Thanks Bakker for your willingness to help. I need to somehow include python code process to the blue color "python code process" as a tool as shown in the figure above. Please download the script, and input files (MajoritySTOht.tif, and viet_LULC_30.tif) using the following web links below. The output of the python code process is "MajoritySTOhtPer.tif" (please look in the code) and will be used to link to Focal Statistics process in the model builder.

https://bft.usu.edu/s467k

https://bft.usu.edu/mt9zj

https://bft.usu.edu/wqlyb

Link to download my ArcGIs tool

https://bft.usu.edu/tlxxz

Thanks and highly appreciate your help.

Re: Run a python code in ArcGIS

XanderBakker — Wed, 14 Jan 2015 20:36:10 GMT

A few remarks:

viet_LULC_30.tif is a large raster (2GB) with a resolution of 30 meters, but the detail inside is only using a resolution of 450m (was this raster resampled with a factor 15? If so, why?). I notice that it has values from -1 to 16. In your code you only account for value 1 to 14.Values -1, 0, 15 and 16 will remain untouched. Is that correct?
One of the downloads (WRFV3.6.1.TAR.gz) seems to be C code. I asume this should not be used.
MajoritySTOht.tif is not included in the downloads. Pease attach the version with 450m resolution to keep things small.
The Toolbox Lee_Filtering contains 5 models. I won't touch these, I will simply show how to include the script in the model
The script SRTMCorrectionTreeOffsets.py contains the code you originally po

Re: Run a python code in ArcGIS

XanderBakker — Sun, 12 Dec 2021 07:43:37 GMT

never mind...

Find attached the toolbox with the script tool inside. You can drag and drop this script tool into your model and connect it:

The code of the script is as follows:

import arcpy
from arcpy import env
from arcpy.sa import *

#Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
arcpy.env.overwriteOutput = True
arcpy.env.workspace = "IN_MEMORY"

# input parameters
LULC = arcpy.GetParameterAsText(0) # r"D:\Xander\GeoNet\LULC\data\viet_LULC_450.tif"
Simard = arcpy.GetParameterAsText(1) # r"D:\Xander\GeoNet\LULC\data\MajoritySTOht_450.tif"
ras_LULC = Raster(LULC)
ras_Simard = Raster(Simard)

# output parameter
result = arcpy.GetParameterAsText(2) # r"D:\Xander\GeoNet\LULC\data\MajoritySTOhtPer.tif"

# internal vars
fld_value = "Value"
fld_remap = "remap"

# remap dictionary
dct = {1 : 0.4, 2 : 0.4, 3 : 0.4, 4 : 0.4, 5 : 0.4, 6 : 0.3, 7 : 0.3,
      8 : 0.3, 9 : 0.3, 10 : 0.3, 11 : 0, 12 : 0.3, 13 : 0, 14 : 0.3}

# add field for remap values
arcpy.AddField_management(ras_LULC, fld_remap, "DOUBLE")

# update the remap values base on dictionary
curs = arcpy.UpdateCursor(ras_LULC)
for row in curs:
    val = row.getValue(fld_value)
    if val in dct:
        remap = dct[val]
    else:
        # what should be done with the values in the TIFF that are not in your list?
        remap = 0 # val?
    row.setValue(fld_remap, remap)
    curs.updateRow(row)
del curs, row

# remap the raster
ras_remap = Lookup(ras_LULC, fld_remap)

# create result
resultmap = ras_Simard * ras_remap
resultmap.save(result)

# release the SA license
arcpy.CheckInExtension("Spatial")

Re: Run a python code in ArcGIS

DuminduJayasekera — Wed, 14 Jan 2015 21:36:14 GMT

Hi Xander,

The process used to obtain viet_LULC_30.tif is correct. As you said it was resampled to 30 m. If you open the attribute table, the values are ranging from 0 to 16 and if you see into layer properties min is 0 and max is 16. So, I don't understand where you saw -1. Yes the code was written to use from 1 to 14 because those are the land use types I am interested.

OK. How can I convert that to python code?

Please see the link to download MajoritySTOht.tif file

https://bft.usu.edu/t3vh9

Please open the Vietnam_Test tool that I am using to perform the GIS process and calculations.

Thanks again,

Re: Run a python code in ArcGIS

DuminduJayasekera — Wed, 14 Jan 2015 22:12:18 GMT

Hi Xander,

I was able to import your tool and tried running. See the picture below. I got an error message (See the image)

As far as I see you are using tif files of 450 m resolution. But input rasters are in 30 m resolutions.

I appreciate your help.

Thanks.

Re: Run a python code in ArcGIS

XanderBakker — Thu, 15 Jan 2015 00:25:58 GMT

Yep that is a memory problem (due to the size of your raster). You can change the "IN_MEMORY" on line 8 to an existing file geodatabase or scratch geodatabase, so it doesn't try to do it in memory

Re: Run a python code in ArcGIS

DuminduJayasekera — Thu, 15 Jan 2015 00:34:14 GMT

Thanks again for the reply. If I revise line 8 by like this will it work? or can I give a path? like this c:\temp\tmp

arcpy.env.workspace = "scratch geodatabase"

Please let me know,

Thanks a lot.