Vjylku

Why doesn't a hydraulic lever violate conservation of energy?

What can I do if neighbor is blocking my solar panels intentionally?

One-dimensional Japanese puzzle

how can a perfect fourth interval be considered either consonant or dissonant?

Is an up-to-date browser secure on an out-of-date OS?

Is it ethical to upload a automatically generated paper to a non peer-reviewed site as part of a larger research?

Keeping a retro style to sci-fi spaceships?

60's-70's movie: home appliances revolting against the owners

How many cones with angle theta can I pack into the unit sphere?

Was credit for the black hole image misappropriated?

Solving overdetermined system by QR decomposition

Single author papers against my advisor's will?

Why did Peik Lin say, "I'm not an animal"?

What does Linus Torvalds mean when he says that Git "never ever" tracks a file?

Do warforged have souls?

Can withdrawing asylum be illegal?

Word for: a synonym with a positive connotation?

Does Parliament need to approve the new Brexit delay to 31 October 2019?

Why can't devices on different VLANs, but on the same subnet, communicate?

What's the point in a preamp?

ELI5: Why do they say that Israel would have been the fourth country to land a spacecraft on the Moon and why do they call it low cost?

How to politely respond to generic emails requesting a PhD/job in my lab? Without wasting too much time

Can I visit the Trinity College (Cambridge) library and see some of their rare books

Is there a writing software that you can sort scenes like slides in PowerPoint?

Creating IDW-rasters from each field in large point feature set using ArcPy for ArcGIS Pro? [on hold]

.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty{ margin-bottom:0;
}

0

I am using ArcGIS Pro 2.3.2

I am trying to automate the process of making a raster (IDW tool) in model builder by iterating through fields as z-value. I saw some other question asking the same, and it seems that it couldn't be done that easy in ModelBuilder, and that Python was the way to go.

I have this point feature of 200+ points that all have around 70 fields of Z-values (elements). In the snip below, you can see fields with elements, and also fields with ID's and other that I don't want to make a IDW raster out of.

So for Lithium, it would produce this raster:

Now, how would I go about doing this process for all fields, except the ones without elements?

If it's easier, I wouldn't mind just deleting the non-element rasters afterwards, if I could maybe name the raster set after the field name perhaps.

arcpy arcgis-pro inverse-distance-weighted

share|improve this question

edited 14 mins ago

PolyGeo♦

53.9k1782246

asked 8 hours ago

DOMINUS MIHI ADIUTORDOMINUS MIHI ADIUTOR

70113

put on hold as off-topic by PolyGeo♦ 13 mins ago

This question appears to be off-topic. The users who voted to close gave this specific reason:

• "When seeking help to debug/write/improve code always provide the desired behavior, a specific problem/error and the shortest code (as formatted text, not pictures) needed to reproduce it in the question body. Providing a clear problem statement and a code attempt helps others to help you." – PolyGeo

If this question can be reworded to fit the rules in the help center, please edit the question.

add a comment | 

0

I am using ArcGIS Pro 2.3.2

I am trying to automate the process of making a raster (IDW tool) in model builder by iterating through fields as z-value. I saw some other question asking the same, and it seems that it couldn't be done that easy in ModelBuilder, and that Python was the way to go.

I have this point feature of 200+ points that all have around 70 fields of Z-values (elements). In the snip below, you can see fields with elements, and also fields with ID's and other that I don't want to make a IDW raster out of.

So for Lithium, it would produce this raster:

Now, how would I go about doing this process for all fields, except the ones without elements?

If it's easier, I wouldn't mind just deleting the non-element rasters afterwards, if I could maybe name the raster set after the field name perhaps.

arcpy arcgis-pro inverse-distance-weighted

share|improve this question

edited 14 mins ago

PolyGeo♦

53.9k1782246

asked 8 hours ago

DOMINUS MIHI ADIUTORDOMINUS MIHI ADIUTOR

70113

put on hold as off-topic by PolyGeo♦ 13 mins ago

This question appears to be off-topic. The users who voted to close gave this specific reason:

• "When seeking help to debug/write/improve code always provide the desired behavior, a specific problem/error and the shortest code (as formatted text, not pictures) needed to reproduce it in the question body. Providing a clear problem statement and a code attempt helps others to help you." – PolyGeo

If this question can be reworded to fit the rules in the help center, please edit the question.

add a comment | 

I am using ArcGIS Pro 2.3.2

I am trying to automate the process of making a raster (IDW tool) in model builder by iterating through fields as z-value. I saw some other question asking the same, and it seems that it couldn't be done that easy in ModelBuilder, and that Python was the way to go.

I have this point feature of 200+ points that all have around 70 fields of Z-values (elements). In the snip below, you can see fields with elements, and also fields with ID's and other that I don't want to make a IDW raster out of.

So for Lithium, it would produce this raster:

Now, how would I go about doing this process for all fields, except the ones without elements?

If it's easier, I wouldn't mind just deleting the non-element rasters afterwards, if I could maybe name the raster set after the field name perhaps.

arcpy arcgis-pro inverse-distance-weighted

share|improve this question

edited 14 mins ago

PolyGeo♦

53.9k1782246

asked 8 hours ago

DOMINUS MIHI ADIUTORDOMINUS MIHI ADIUTOR

70113

share|improve this question

edited 14 mins ago

PolyGeo♦

53.9k1782246

asked 8 hours ago

DOMINUS MIHI ADIUTORDOMINUS MIHI ADIUTOR

70113

share|improve this question

edited 14 mins ago

PolyGeo♦

53.9k1782246

asked 8 hours ago

DOMINUS MIHI ADIUTORDOMINUS MIHI ADIUTOR

70113

edited 14 mins ago

PolyGeo♦

53.9k1782246

edited 14 mins ago

PolyGeo♦

53.9k1782246

asked 8 hours ago

DOMINUS MIHI ADIUTORDOMINUS MIHI ADIUTOR

70113

asked 8 hours ago

DOMINUS MIHI ADIUTORDOMINUS MIHI ADIUTOR

70113

1 Answer
1

active

oldest

votes

1

The tricky part is to only list the element fields and no other. You can do this in multiple ways, for example:

• Manually listing the fields which are not element fields, example below. Or manually list all element fields


• Use some wildcard with listfiels, for example if all your element fields start with some letter


• List fields with names shorter than 4 letters

Example (untested, I dont have spatial analyst, you will probably need to look at the syntax for IDW and adjust the parameters):

import arcpy

from os.path import join



fc = r'C:data.gdbfeatureclass' #Change

out_gdb = r'C:outdata.gdb' #Change

not_elementfields = ['OBJECTID','SHAPE','RowID'] #Change



elementfields = [field.name for field in arcpy.ListFields(fc) if f.name not in not_elementfields]

#elementfields = [field.name for field in arcpy.ListFields(fc) if len(f.name)<4]

print(elementfields)



#Run code above first to get the listing correct, then run everything



for field in elementfields:

    print(f'Proccessing IDW for fieldname {field}')

    outIDW = arcpy.sa.Idw(in_point_features=fc, z_field=field, cell_size=5, power=0.5, 

                search_radius= RadiusVariable(12))

    outIDW.save(join(out_gdb,f'IDW_{field}'))

(This is also possible:

import pandas as pd



data = pd.read_csv('https://sciencenotes.org/PDFs/elementlist.csv', names = ['number','symbol','element'])

elements = list(data.symbol)



elementfields = [f.name for f in arcpy.ListFields(fc) if f.name.title() in elements] # title method will camelcase your field names, BE->Be, to match the csv file

)

share|improve this answer

edited 3 hours ago

answered 8 hours ago

BERABERA

17.2k62044

add a comment | 

1

The tricky part is to only list the element fields and no other. You can do this in multiple ways, for example:

• Manually listing the fields which are not element fields, example below. Or manually list all element fields


• Use some wildcard with listfiels, for example if all your element fields start with some letter


• List fields with names shorter than 4 letters

Example (untested, I dont have spatial analyst, you will probably need to look at the syntax for IDW and adjust the parameters):

import arcpy

from os.path import join



fc = r'C:data.gdbfeatureclass' #Change

out_gdb = r'C:outdata.gdb' #Change

not_elementfields = ['OBJECTID','SHAPE','RowID'] #Change



elementfields = [field.name for field in arcpy.ListFields(fc) if f.name not in not_elementfields]

#elementfields = [field.name for field in arcpy.ListFields(fc) if len(f.name)<4]

print(elementfields)



#Run code above first to get the listing correct, then run everything



for field in elementfields:

    print(f'Proccessing IDW for fieldname {field}')

    outIDW = arcpy.sa.Idw(in_point_features=fc, z_field=field, cell_size=5, power=0.5, 

                search_radius= RadiusVariable(12))

    outIDW.save(join(out_gdb,f'IDW_{field}'))

(This is also possible:

import pandas as pd



data = pd.read_csv('https://sciencenotes.org/PDFs/elementlist.csv', names = ['number','symbol','element'])

elements = list(data.symbol)



elementfields = [f.name for f in arcpy.ListFields(fc) if f.name.title() in elements] # title method will camelcase your field names, BE->Be, to match the csv file

)

share|improve this answer

edited 3 hours ago

answered 8 hours ago

BERABERA

17.2k62044

add a comment | 

1

The tricky part is to only list the element fields and no other. You can do this in multiple ways, for example:

• Manually listing the fields which are not element fields, example below. Or manually list all element fields


• Use some wildcard with listfiels, for example if all your element fields start with some letter


• List fields with names shorter than 4 letters

Example (untested, I dont have spatial analyst, you will probably need to look at the syntax for IDW and adjust the parameters):

import arcpy

from os.path import join



fc = r'C:data.gdbfeatureclass' #Change

out_gdb = r'C:outdata.gdb' #Change

not_elementfields = ['OBJECTID','SHAPE','RowID'] #Change



elementfields = [field.name for field in arcpy.ListFields(fc) if f.name not in not_elementfields]

#elementfields = [field.name for field in arcpy.ListFields(fc) if len(f.name)<4]

print(elementfields)



#Run code above first to get the listing correct, then run everything



for field in elementfields:

    print(f'Proccessing IDW for fieldname {field}')

    outIDW = arcpy.sa.Idw(in_point_features=fc, z_field=field, cell_size=5, power=0.5, 

                search_radius= RadiusVariable(12))

    outIDW.save(join(out_gdb,f'IDW_{field}'))

(This is also possible:

import pandas as pd



data = pd.read_csv('https://sciencenotes.org/PDFs/elementlist.csv', names = ['number','symbol','element'])

elements = list(data.symbol)



elementfields = [f.name for f in arcpy.ListFields(fc) if f.name.title() in elements] # title method will camelcase your field names, BE->Be, to match the csv file

)

share|improve this answer

edited 3 hours ago

answered 8 hours ago

BERABERA

17.2k62044

add a comment | 

The tricky part is to only list the element fields and no other. You can do this in multiple ways, for example:

• Manually listing the fields which are not element fields, example below. Or manually list all element fields


• Use some wildcard with listfiels, for example if all your element fields start with some letter


• List fields with names shorter than 4 letters

Example (untested, I dont have spatial analyst, you will probably need to look at the syntax for IDW and adjust the parameters):

import arcpy

from os.path import join



fc = r'C:data.gdbfeatureclass' #Change

out_gdb = r'C:outdata.gdb' #Change

not_elementfields = ['OBJECTID','SHAPE','RowID'] #Change



elementfields = [field.name for field in arcpy.ListFields(fc) if f.name not in not_elementfields]

#elementfields = [field.name for field in arcpy.ListFields(fc) if len(f.name)<4]

print(elementfields)



#Run code above first to get the listing correct, then run everything



for field in elementfields:

    print(f'Proccessing IDW for fieldname {field}')

    outIDW = arcpy.sa.Idw(in_point_features=fc, z_field=field, cell_size=5, power=0.5, 

                search_radius= RadiusVariable(12))

    outIDW.save(join(out_gdb,f'IDW_{field}'))

(This is also possible:

import pandas as pd



data = pd.read_csv('https://sciencenotes.org/PDFs/elementlist.csv', names = ['number','symbol','element'])

elements = list(data.symbol)



elementfields = [f.name for f in arcpy.ListFields(fc) if f.name.title() in elements] # title method will camelcase your field names, BE->Be, to match the csv file

)

share|improve this answer

edited 3 hours ago

answered 8 hours ago

BERABERA

17.2k62044

import arcpy

from os.path import join



fc = r'C:data.gdbfeatureclass' #Change

out_gdb = r'C:outdata.gdb' #Change

not_elementfields = ['OBJECTID','SHAPE','RowID'] #Change



elementfields = [field.name for field in arcpy.ListFields(fc) if f.name not in not_elementfields]

#elementfields = [field.name for field in arcpy.ListFields(fc) if len(f.name)<4]

print(elementfields)



#Run code above first to get the listing correct, then run everything



for field in elementfields:

    print(f'Proccessing IDW for fieldname {field}')

    outIDW = arcpy.sa.Idw(in_point_features=fc, z_field=field, cell_size=5, power=0.5, 

                search_radius= RadiusVariable(12))

    outIDW.save(join(out_gdb,f'IDW_{field}'))

import pandas as pd



data = pd.read_csv('https://sciencenotes.org/PDFs/elementlist.csv', names = ['number','symbol','element'])

elements = list(data.symbol)



elementfields = [f.name for f in arcpy.ListFields(fc) if f.name.title() in elements] # title method will camelcase your field names, BE->Be, to match the csv file

share|improve this answer

edited 3 hours ago

answered 8 hours ago

BERABERA

17.2k62044

answered 8 hours ago

BERABERA

17.2k62044

answered 8 hours ago

BERABERA

17.2k62044