|
POST
|
I am trying in my best, i am posting here what i have so far understood, can anyone suggest me how to iterate Friging function to all point data points = arcpy.ListFeatureClasses('*','POINT')
for fc in points:
path = os.path.join(arcpy.env.workspace, fc)
Point = Point_Num[3:]
... View more
11-22-2016
07:50 PM
|
0
|
6
|
3936
|
|
POST
|
Dan , i have made changes in line no 28, as i changed it , for fc in points: but still error i am getting. and error message is "object is not iterable".
... View more
11-22-2016
11:22 AM
|
0
|
1
|
3936
|
|
POST
|
Hi, All Please provide your valuable input to fix the error. Thank you
... View more
11-22-2016
09:45 AM
|
0
|
3
|
3936
|
|
POST
|
Dan , thank you. that link for multiple Z fields in one point feature class, but here i have multiple point feature class and each point feature class have one z field.
... View more
11-22-2016
08:11 AM
|
0
|
0
|
3936
|
|
POST
|
Thank you. yes, i am getting the error at same line which you mentioned, how to fix the error
... View more
11-22-2016
08:00 AM
|
0
|
0
|
3936
|
|
POST
|
I have multiple point feature class (r001_mean, r002_mean.....012_mean) into one folder. I am trying to write a python script to perform Kriging interpolation which will loop all the point feature class and outraster will save as it name of input raster. Each Point feature class hold one Z filed only (GRID_CODE) Below is my script, what i have so far.....its giving error massage of "not defined point_num" import arcpy
from arcpy import env
from arcpy.sa import *
arcpy.env.overwriteOutput = True
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
In_Point = r"D:\NPP_Test\MERRA_TEMP_2012C" #(Point feature name:r001_mean, r002_mean.....r012_mean )
Out_Raster = r"D:\NPP_Test\MERRA_TEMP_2012D"
points = arcpy.ListFeatureClasses()
zFields = "GRID_CODE"
#Kriging Veriable
cellSize = 0.05
lagSize = 0.5780481172534
majorRange = 6
partialSill = 3.304292110
nugget = 0.002701348
kRadius = RadiusFixed(20000, 1)
#Mask region of interest
mask = r"D:\Gujarta Shape file\GUJARATSTATE.shp"
for zField in zFields:
Point = Point_Num[:3]
kModelUniversalObj = KrigingModelUniversal("LINEARDRIFT", lagSize, majorRange, partialSill, nugget)
OutKriging = Kriging(inPointFeatures, zField, kModelUniversalObj, cellSize, kRadius)
#IDWMASk = ExtractByMask(outIDW, mask)
KrigMask = ExtractByMask(OutKriging, mask)
#Save outraster as the same name of input
KrigMask.save("r{}.tif".format(Point_Num))
... View more
11-22-2016
07:23 AM
|
0
|
25
|
10321
|
|
POST
|
xander_bakker. Thank you very much. Now script running perfectly and I got absolutely correct output from above script. Once again thank you for your kind cooperation.
... View more
10-17-2016
05:13 AM
|
0
|
2
|
859
|
|
POST
|
Xander Bakker. Thank you. I am using r001_NPP and r003_NPP, When i mentioned Random_pnts.shp then i am getting error Input : 1 ras1 = r'D:\TEST\Reference_RASTER\NPP_2008\r001_NPP.TIF'
ras2 = r'D:\TEST\Reference_RASTER\NPP_2008\r003_NPP.TIF'
fc_pnt = r'D:\TEST\Reference_RASTER\Random_Points\Random_pnts.shp'
num_points = 1000 Message File Name Line Position
Traceback
<module> <module1> 45
main <module1> 34
RuntimeError: A column was specified that does not exist. Input: 2 And when i trying to put another name (Randompnts) of point feature class using same path of exiting point feature class , then its creating the same name of point feature class as mentioned name (Randompnts) ras1 = r'D:\TEST\Reference_RASTER\NPP_2008\r001_NPP.TIF'
ras2 = r'D:\TEST\Reference_RASTER\NPP_2008\r003_NPP.TIF'
fc_pnt = r'D:\TEST\Reference_RASTER\Random_Points\Randompnts'
num_points = 1000 but after creating point feature class its giving error massage Traceback (most recent call last):
File "D:\Arc-GIS-Python-Script\Scatter_Plot.py", line 45, in <module>
main()
File "D:\Arc-GIS-Python-Script\Scatter_Plot.py", line 30, in main
arcpy.sa.ExtractMultiValuesToPoints(fc_pnt, in_rasters, "NONE")
File "C:\Program Files\ArcGIS\Desktop10.3\ArcPy\arcpy\sa\Functions.py", line 7180, in ExtractMultiValuesToPoints
bilinear_interpolate_values)
File "C:\Program Files\ArcGIS\Desktop10.3\ArcPy\arcpy\sa\Utils.py", line 53, in swapper
result = wrapper(*args, **kwargs)
File "C:\Program Files\ArcGIS\Desktop10.3\ArcPy\arcpy\sa\Functions.py", line 7175, in Wrapper
bilinear_interpolate_values)
File "C:\Program Files\ArcGIS\Desktop10.3\ArcPy\arcpy\geoprocessing\_base.py", line 504, in <lambda>
return lambda *args: val(*gp_fixargs(args, True))
ExecuteError: Failed to execute. Parameters are not valid.
ERROR 000865: Input point features: D:\TEST\Reference_RASTER\Random_Points\Randompnts does not exist.
Failed to execute (ExtractMultiValuesToPoints). Herewith i have attached raster and vector file which i am using for scatter plotting.
... View more
10-16-2016
05:42 AM
|
0
|
4
|
859
|
|
POST
|
Dan Patterson. Thank you. Firstly i tried without dissolve, after read your comment i tried dissolve also, but problem still persist. I just explaining short what i am doing, Firstly i created random points using create random points tool----------> After i copied that Random points path directly from folder and use as "fc_pnt" path --------> Run the script I have noticed, while running the program it is generating same point file in same location, I have followed dissolve also after generating random points but same error i am getting.
... View more
10-12-2016
12:38 AM
|
0
|
2
|
2544
|
|
POST
|
Xander Bakker. Thank you. Excellent idea. Really its very helpful for research community. I have created random points (Number of points:1000) using create random tool under feature class. and I have given that point feature class as fc_pnt path, fc_pnt = r'D:\TEST\Reference_RASTER\Random_Points\Random_pnts.shp' i am getting error, might be i didn't construct fc_pnt data properly Traceback (most recent call last):
File "<module1>", line 45, in <module>
File "<module1>", line 34, in main
RuntimeError: A column was specified that does not exist. Hope i have problem with my point feature data set. may i know how to construct that point feature class. Thank you.
... View more
10-11-2016
10:40 PM
|
0
|
10
|
2544
|
|
POST
|
Dan Patterson Thank you. Yes . This graph require for my research project. I am PhD scholar. I am going to publish a paper How Temperate effects on vegetation growth. So i want a scatter diagram between two variable so i can say to research community by the statistical plot. I have done this task suggestion by you Xander Bakker and @Neil Ayres. Thank you very much.
... View more
10-11-2016
11:10 AM
|
0
|
12
|
2544
|
|
POST
|
xander_bakker, Thank you. For my case i need to compare all the pixels corresponding to same latitude and longitude.
... View more
10-11-2016
07:07 AM
|
0
|
1
|
2544
|
|
POST
|
Dan Patterson, As your last comment we have huge data point. number of columns and rows 2737, 2025 respectively. So whats is the best way to plot those raster.
... View more
10-11-2016
06:01 AM
|
0
|
3
|
7245
|
|
POST
|
Xander Bakker, Dan Patterson , I moved total conversion from there to here. Please cooperate with us. Sorry for the inconvenience cause.
... View more
10-11-2016
05:58 AM
|
0
|
4
|
7245
|
| Title | Kudos | Posted |
|---|---|---|
| 1 | 02-01-2017 09:56 PM | |
| 1 | 01-19-2017 08:03 AM | |
| 1 | 07-31-2017 11:57 AM | |
| 1 | 08-25-2016 09:38 AM | |
| 1 | 08-17-2016 11:13 AM |
| Online Status |
Offline
|
| Date Last Visited |
11-11-2020
02:25 AM
|