POST
|
The RMSE does not appear in the Control Point Table like it did in ArcGIS Desktop. Where can I see the RMSE in ArcGIS Pro?
... View more
08-16-2017
03:35 PM
|
1
|
0
|
441
|
POST
|
Cody B., I have used SEE5 in the past, https://www.rulequest.com/see5-win.html, to perform similar classification and regression tree (CART) based classifications. When I use SEE5 it outputs the trees that are being used to distinguish the different classes. I was hoping to teach my students about CART based classification using a small amount of data so that they could dissect the trees. Also, if I were to ever use the Random Tree classifier for my own research I would like to see how the data is being used. It can be informative to know what bands are helping delineate classes in case retraining is necessary. Kyle
... View more
12-15-2016
01:18 PM
|
1
|
0
|
259
|
POST
|
I would like to actually inspect the trees created by the Random Trees Classifier. I have opened the .ecd file created by the training process and some of it is legible, but the bottom part is not. I am guessing those are the trees. I would actually like to see how classes are being separated and what values from each input are being used to split and classify the various classes.
... View more
12-15-2016
09:10 AM
|
0
|
2
|
711
|
POST
|
Thanks to all those who chimed in to help. I have attached the code I finally got to work. At this point in my python education this code makes the most sense in my head. I understand that there is a more efficient way to do this and I am open to learning how if you can point me in the right direction. I explored the .format tool and could not get it to loop the way I wanted to. Thanks again for all the help. # Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *
# Check out any necessary licenses
arcpy.CheckOutExtension("spatial")
# Path names for input and output data
path = "S:\\MODIS\\DATA_v05\\modis016_0250_ndvi\\us\\or\\or"
output_path = "B:\\MODIS_NDVI\\Avg_per_period\\avg"
for period in range (1,24):
#MODIS does not have data for the first three periods of 2000. This will calculate an average for the first three periods ignoring the year 2000. It will also place a 0 in front of the single digit to match the file name.
if period <= 3:
period = ("0"+str(period))
grids01 = Raster(path + "01" + period + "ndvius")
grids02 = Raster(path + "02" + period + "ndvius")
grids03 = Raster(path + "03" + period + "ndvius")
grids04 = Raster(path + "04" + period + "ndvius")
grids05 = Raster(path + "05" + period + "ndvius")
grids06 = Raster(path + "06" + period + "ndvius")
grids07 = Raster(path + "07" + period + "ndvius")
grids08 = Raster(path + "08" + period + "ndvius")
grids09 = Raster(path + "09" + period + "ndvius")
grids10 = Raster(path + "10" + period + "ndvius")
grids11 = Raster(path + "11" + period + "ndvius")
grids12 = Raster(path + "12" + period + "ndvius")
grids = [grids01, grids02, grids03, grids04, grids05, grids06, grids07, grids08, grids09, grids10, grids11, grids12]
avg_temp = (output_path + period + "temp")
avg = (output_path + period)
arcpy.gp.CellStatistics_sa(grids, avg_temp,"MEAN", "DATA")
print ("Cell Stats complete for" + period)
arcpy.gp.Int_sa(avg_temp, avg)
print ("Int complete for" + period)
#This is will place a 0 in front of the single digit to match the file name.
elif period < 10:
period = ("0"+str(period))
grids00 = Raster(path + "00" + period + "ndvius")
grids01 = Raster(path + "01" + period + "ndvius")
grids02 = Raster(path + "02" + period + "ndvius")
grids03 = Raster(path + "03" + period + "ndvius")
grids04 = Raster(path + "04" + period + "ndvius")
grids05 = Raster(path + "05" + period + "ndvius")
grids06 = Raster(path + "06" + period + "ndvius")
grids07 = Raster(path + "07" + period + "ndvius")
grids08 = Raster(path + "08" + period + "ndvius")
grids09 = Raster(path + "09" + period + "ndvius")
grids10 = Raster(path + "10" + period + "ndvius")
grids11 = Raster(path + "11" + period + "ndvius")
grids12 = Raster(path + "12" + period + "ndvius")
grids = [grids00, grids01, grids02, grids03, grids04, grids05, grids06, grids07, grids08, grids09, grids10, grids11, grids12]
avg_temp = (output_path + period + "temp")
avg = (output_path + period)
arcpy.gp.CellStatistics_sa(grids, avg_temp,"MEAN", "DATA")
print ("Cell Stats complete for" + period)
arcpy.gp.Int_sa(avg_temp, avg)
print ("Int complete for" + period)
#These are double digit periods and will match the file name without correction.
else:
period = (str(period))
grids00 = Raster(path + "00" + period + "ndvius")
grids01 = Raster(path + "01" + period + "ndvius")
grids02 = Raster(path + "02" + period + "ndvius")
grids03 = Raster(path + "03" + period + "ndvius")
grids04 = Raster(path + "04" + period + "ndvius")
grids05 = Raster(path + "05" + period + "ndvius")
grids06 = Raster(path + "06" + period + "ndvius")
grids07 = Raster(path + "07" + period + "ndvius")
grids08 = Raster(path + "08" + period + "ndvius")
grids09 = Raster(path + "09" + period + "ndvius")
grids10 = Raster(path + "10" + period + "ndvius")
grids11 = Raster(path + "11" + period + "ndvius")
grids12 = Raster(path + "12" + period + "ndvius")
grids = [grids00, grids01, grids02, grids03, grids04, grids05, grids06, grids07, grids08, grids09, grids10, grids11, grids12]
avg_temp = (output_path + period + "temp")
avg = (output_path + period)
arcpy.gp.CellStatistics_sa(grids, avg_temp,"MEAN", "DATA")
print ("Cell Stats complete for" + period)
arcpy.gp.Int_sa(avg_temp, avg)
print ("Int complete for" + period)
... View more
05-29-2013
07:40 AM
|
0
|
0
|
563
|
POST
|
Do you actually have a "raster" data set named S:\\MODIS\\DATA_v05\\modis016_0250_ndvi\\us\\or0101ndvius ? If so, what format is it, and you might try to remove the "'" from both ends of your grids## = assignments. If I type them in as you have them, I get " ' S:\\MODIS\\DATA_v05\\modis016_0250_ndvi\\us\\or0101ndvius ' " and it might be the extra quotes causing the issue. I don't have the data, so can't really test. If you don't have a raster with that name, that is what it is looking for so you might check to make sure the resultant path actually matches an existing raster. Running under the 64 bit IDE can give this error as well. Normally it is when you are trying to access something through an ArcCatalog connection though (since that only supports 32 bit, the 64 bit one can't find the ArcCatalog/connectionto/ files). R_ Rzufelt, That is the name of my raster. I have attempted various formats of naming and was trying to copy how the command looks when I export it from model builder. I removed the additional quotes and received the same error. The machine I am running is a 64 bit machine, but I am not accessing the data through ArcCatalog. Is there a way to work around the 64 bit problem. I have included the piece of my adjusted script. # Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *
# Set environment settings
env.workspace = "S:\\MODIS\\DATA_v05\\modis016_0250_ndvi\\us\\"
# Check out any necessary licenses
arcpy.CheckOutExtension("spatial")
# Path names for input and output data
path = "S:\\MODIS\\DATA_v05\\modis016_0250_ndvi\\us\\or"
output_path = "B:\\MODIS_NDVI\\Avg_per_period\\avg"
for period in range (1,3):
#MODIS does not have data for the first three periods of 2000. This will calculate an average for the first three periods ignoring the year 2000. It will also place a 0 in front of the single digit to match the file name.
if period <= 3:
period = ("0"+str(period))
grids01 = Raster(path + "01" + period + "ndvius")
grids02 = Raster(path + "02" + period + "ndvius")
grids03 = Raster(path + "03" + period + "ndvius")
grids04 = Raster(path + "04" + period + "ndvius")
grids05 = Raster(path + "05" + period + "ndvius")
grids06 = Raster(path + "06" + period + "ndvius")
grids07 = Raster(path + "07" + period + "ndvius")
grids08 = Raster(path + "08" + period + "ndvius")
grids09 = Raster(path + "09" + period + "ndvius")
grids10 = Raster(path + "10" + period + "ndvius")
grids11 = Raster(path + "11" + period + "ndvius")
grids12 = Raster(path + "12" + period + "ndvius")
#grids = grids01 + ";" + grids02 + ";" + grids03 + ";" + grids04 + ";" + grids05 + ";" + grids06 + ";" + grids07 + ";" + grids08 + ";" + grids09 + ";" + grids10 + ";" + grids11 + ";" + grids12
avg_temp = ('"' + output_path + period + "temp" + '"')
avg = (output_path + period)
outCellStatistics = CellStatistics([grids01, grids02, grids03, grids04, grids05, grids06, grids07, grids08, grids09, grids10, grids11, grids12], "MEAN", "DATA")
outCellStatistics.save(avg_temp)
... View more
05-24-2013
07:05 AM
|
0
|
0
|
563
|
POST
|
I am attempting to loop through a set of MODIS data to compute averages by period. Since 2000 is missing periods 1 - 3 I am just computing the averages of periods 1 -3 for 2001 through 2012. I am then computing the averages for periods 4 - 23 with data from 2000 through 2012. I have written a script to do this, but keep getting an error 000732 about my data not existing or not being supported. I opened one of the grids in Arc and it was fine. I am not sure why I am getting this error and have included the code below. # Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *
# Set environment settings
env.workspace = "S:\\MODIS\\DATA_v05\\modis016_0250_ndvi\\us\\"
# Check out any necessary licenses
arcpy.CheckOutExtension("spatial")
# Path names for input and output data
path = "S:\\MODIS\\DATA_v05\\modis016_0250_ndvi\\us\\or"
output_path = "B:\\MODIS_NDVI\\Avg_per_period\\avg"
for period in range (1,3):
#MODIS does not have data for the first three periods of 2000. This will calculate an average for the first three periods ignoring the year 2000. It will also place a 0 in front of the single digit to match the file name.
if period <= 3:
period = ("0"+str(period))
grids01 = Raster('"' + path + "01" + period + "ndvius" + '"')
grids02 = Raster('"' + path + "02" + period + "ndvius" + '"')
grids03 = Raster('"' + path + "03" + period + "ndvius" + '"')
grids04 = Raster('"' + path + "04" + period + "ndvius" + '"')
grids05 = Raster('"' + path + "05" + period + "ndvius" + '"')
grids06 = Raster('"' + path + "06" + period + "ndvius" + '"')
grids07 = Raster('"' + path + "07" + period + "ndvius" + '"')
grids08 = Raster('"' + path + "08" + period + "ndvius" + '"')
grids09 = Raster('"' + path + "09" + period + "ndvius" + '"')
grids10 = Raster('"' + path + "10" + period + "ndvius" + '"')
grids11 = Raster('"' + path + "11" + period + "ndvius" + '"')
grids12 = Raster('"' + path + "12" + period + "ndvius" + '"')
#grids = grids01 + ";" + grids02 + ";" + grids03 + ";" + grids04 + ";" + grids05 + ";" + grids06 + ";" + grids07 + ";" + grids08 + ";" + grids09 + ";" + grids10 + ";" + grids11 + ";" + grids12
avg_temp = ('"' + output_path + period + "temp" + '"')
avg = (output_path + period)
outCellStatistics = CellStatistics([grids01, grids02, grids03, grids04, grids05, grids06, grids07, grids08, grids09, grids10, grids11, grids12], "MEAN", "DATA")
outCellStatistics.save(avg_temp)
#This is will place a 0 in front of the single digit to match the file name.
elif period < 10:
period = ("0"+str(period))
grids00 = ('"' + path + "00" + period + "ndvius" + '"')
grids01 = ('"' + path + "01" + period + "ndvius" + '"')
grids02 = ('"' + path + "02" + period + "ndvius" + '"')
grids03 = ('"' + path + "03" + period + "ndvius" + '"')
grids04 = ('"' + path + "04" + period + "ndvius" + '"')
grids05 = ('"' + path + "05" + period + "ndvius" + '"')
grids06 = ('"' + path + "06" + period + "ndvius" + '"')
grids07 = ('"' + path + "07" + period + "ndvius" + '"')
grids08 = ('"' + path + "08" + period + "ndvius" + '"')
grids09 = ('"' + path + "09" + period + "ndvius" + '"')
grids10 = ('"' + path + "10" + period + "ndvius" + '"')
grids11 = ('"' + path + "11" + period + "ndvius" + '"')
grids12 = ('"' + path + "12" + period + "ndvius" + '"')
grids = "[" + grids00 + "," + grids01 + "," + grids02 + "," + grids03 + "," + grids04 + "," + grids05 + "," + grids06 + "," + grids07 + "," + grids08 + "," + grids09 + "," + grids10 + "," + grids11 + "," + grids12 + "]"
avg_temp = (output_path + period + "temp")
avg = (output_path + period)
arcpy.gp.CellStatistics_sa(grids,avg_temp,"MEAN", "DATA")
#These are double digit periods and will match the file name without correction.
else:
period = (str(period))
grids00 = ('"' + path + "00" + period + "ndvius" + '"')
grids01 = ('"' + path + "01" + period + "ndvius" + '"')
grids02 = ('"' + path + "02" + period + "ndvius" + '"')
grids03 = ('"' + path + "03" + period + "ndvius" + '"')
grids04 = ('"' + path + "04" + period + "ndvius" + '"')
grids05 = ('"' + path + "05" + period + "ndvius" + '"')
grids06 = ('"' + path + "06" + period + "ndvius" + '"')
grids07 = ('"' + path + "07" + period + "ndvius" + '"')
grids08 = ('"' + path + "08" + period + "ndvius" + '"')
grids09 = ('"' + path + "09" + period + "ndvius" + '"')
grids10 = ('"' + path + "10" + period + "ndvius" + '"')
grids11 = ('"' + path + "11" + period + "ndvius" + '"')
grids12 = ('"' + path + "12" + period + "ndvius" + '"')
grids = "[" + grids00 + "," + grids01 + "," + grids02 + "," + grids03 + "," + grids04 + "," + grids05 + "," + grids06 + "," + grids07 + "," + grids08 + "," + grids09 + "," + grids10 + "," + grids11 + "," + grids12 + "]"
avg_temp = (output_path + period + "temp")
avg = (output_path + period)
arcpy.gp.CellStatistics_sa(grids,avg_temp,"MEAN", "DATA")
... View more
05-23-2013
02:55 PM
|
0
|
6
|
783
|
Title | Kudos | Posted |
---|---|---|
1 | 12-15-2016 01:18 PM | |
1 | 08-16-2017 03:35 PM |
Online Status |
Offline
|
Date Last Visited |
11-11-2020
02:23 AM
|