POST
|
Did you ever find any resolution on this? I am running into the same issue.
... View more
06-22-2015
06:45 AM
|
0
|
0
|
466
|
POST
|
Hello, I am conducting a grouping analysis of statewide management units based on a suite of relevant covariates or attributes for each management unit. Further, I have some a priori criteria that already splits the state into 3 groups. The analysis will create subgroups within these 3 parent groups. I have tried running the analysis by splitting the 3 groups into seperate shapefiles and running the analysis for each shapefile. However, the tool does not like it when there is less than 30 features and takes exponentially longer. I think I need to create a spatial weights matrix, but I am not sure how. I tried tricking the tool by creating a date variable and giving features belonging to the same parent group an equivalent date and then creating a spatial weights matrix using the time and space constraint. However, I think this is less than ideal because it does not consider any data that does not have a similar date. Anyway, can anyone offer any advice on how to create a swm that would maintain the 3 parent groups?? Thank you!
... View more
11-21-2014
09:32 AM
|
0
|
0
|
2753
|
POST
|
Hey Christopher, Did you ever resolve this issue? I am running into the same issue now and was wondering if you could add any insight. I suspect it has something to do with the NoData setting, but I am not sure how to correct it. Thanks!
... View more
10-27-2014
12:33 PM
|
0
|
1
|
178
|
POST
|
Hi Gerry and Neil, Thanks for looking into my issue. Either I was unclear in my explanation or I am not understanding your explanation. I am not sure how this will help me make a separate list for each year's winter in the dataset. There are 11 years worth of data in the same folder. The only months I downloaded for each year are 11,12,1,2 and 3 so I don't need to separate out winter months. What I need to separate are the years(e.g., winter of '03-'04). This is a challenge because months 1, 2, and 3 need to be assigned to the previous calendar year. Any thoughts? Thanks again!
... View more
10-24-2014
12:57 PM
|
0
|
2
|
791
|
POST
|
Hello, I am a beginner at python programming, and this is my first post to GeoNet. So here it goes! I have a huge list of rasters downloaded from the SNODAS website for years 2003 - present (intra-annual range of 11/1-3/31). I have a .tif file for everyday during this period which contains the snow depth data for that day. The metadata for each raster is stored in the tif file's name(e.g., us_ssmv11036tS__T0001TTNATS2013111305HP001.dat.reproj.tif ). However, I am only interested in the date (Nov. 13, 2013 in this case). My objective is to split this huge list of rasters into winter-years (i.e., rasters from Jan-Mar of a given year belong to the winter-year of the previous calendar year) so that I can pass each winter-year list of rasters through the greaterThanFrequency tool with a simple loop. Below is a script that works if I download one year of data at a time. I'd appreciate any feedback on how get a list of rasters for each winter-year.
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "E:\DECGIS\snow\NSIDC\projects" \
"\Searchlight\production\data\NSIDC_Data_1021063617170" \
"\Snow_Data_Assimilation_System_(SNODAS)_Data_Products_at_NSIDC"
env.overwriteOutput = True
arcpy.env.snapRaster = "us_ssmv11036tS__T0001TTNATS2013111305HP001.dat.reproj.tif"
arcpy.env.outputCoordinateSystem = "us_ssmv11036tS__T0001TTNATS2013111305HP001.dat.reproj.tif"
arcpy.env.nodata = "-9999"
inRasters = arcpy.ListRasters("*", "TIF")
arcpy.CheckOutExtension("Spatial")
#create raster of constant values for greaterThanFrequency tool
constantValue = "381"
outpath3 = str(env.workspace+"\\threshold4.tif")
cellSize = '8.33333333813986E-03'
outConstRaster = CreateConstantRaster(constantValue, 'INTEGER', cellSize,"us_ssmv11036tS__T0001TTNATS2013111305HP001.dat.reproj.tif")
outConstRaster.save(outpath3)
#run greater than frequency analysis
inValueRaster = "threshold4.tif"
outGTF = GreaterThanFrequency(inValueRaster, inRasters)
# Save the output
outGTF.save(r"E:\DECGIS\snow\output\2013daysGt15.tif")
... View more
10-22-2014
09:10 AM
|
0
|
5
|
2206
|
Online Status |
Offline
|
Date Last Visited |
05-19-2023
02:57 PM
|