Latest Contributions by RyanDeBruyn

Note: You can't save .tif format to output GDB workspace.  Use full path when saving or verify your scratchWorkspace is set to file directory if using .tif or other file based extensions. ... View more

@rytz  Thank you for your post and great investigation. With respect to your OP... "When I use Python, the output raster is named with the name of the variable I used to define the name, not the name (defined as a string) itself."  Please understand that this is "by design" in that specifically when working with Spatial Analyst tools and operations in Python you are performing Map Algebra.  Any tool or operator that produces an output raster to the left of the equal sign creates a raster object. For example, in the following expression, out_raster is a raster object. out_raster = Slope("inelevation") When a Raster object is returned from a Map Algebra expression, by default, the object (the variable and associated dataset) is temporary. The temporary dataset associated with a raster object can become permanent by calling the raster object's save method.  If the referenced raster is not made permanent, the variable and the referenced raster dataset will be deleted when the variable goes out of scope, such as when a stand-alone script completes or ArcGIS is closed.  In your case, working within the Interactive Python window inside Pro automatically adds the temporary raster to the display (there is an option to turn this off) and we honour the variable name (raster object).  When you save, the data thus becomes permanent and the output name is updated internally. Hope this helps!  Good luck! Here are some additional resources for further review. Raster Object Working-with-raster-objects  ... View more

You can use the parameter each_interval key to handle whether to calculate a single total insolation value for all locations or multiple values for the specified hour and day interval. NOINTERVAL —A single total radiation value will be calculated for the entire time configuration. This is the default. INTERVAL —Multiple radiation values will be calculated for each time interval over the entire time configuration. The number of outputs will depend on the hour or day interval. For example, for a whole year with monthly intervals, the result will contain 12 output radiation values for each location. Are you creating a new output with unique name for each iteration? ... View more

Make sure your environment workspace and output directories have write capabilities.   You should be able to verify running the tool (or any tool) in the desktop client setting the output and environment workspace/scratchWorkspace respectively. ... View more

Seems like there is some disconnect with the defined spatial reference and how the spatial points are being created.  There are different ways to create input point layers from a set of coordinates if you have them without doing it manually.  Check out the Make XY Event Layer tool  if you have a list of X,Y coordinates.   Can also take advantage or convert Excel To Table (Conversion) tool.   You can then use Project or CopyFeatures (with environment setting)  to set the appropriate coordinate system for you final output.  If you are continuing to have issues I suggest contacting Esri technical support to help you with your data and your process.  ... View more

The error points to a data issue.  You'll need to ensure the point data were created in the same projection as the DEM.  How were the point data created?   ... View more

HI @RabeccaThiessen , Please can you verify if your points and raster have or are in the same spatial reference? Does it work with any other data or if you test a smaller DEM areas?  For example you can zoom into a small area and set environment extent="same as display". ... View more

So you are wanting to change the radiation parameters with each 1/2 hour iteration running the the tool multiple times? You will be setting the time_configuration parameter using the TimeWithinDay({day},{startTime},{endTime}) option which will allow you to set the intervals accordingly.  You could also set the each_interval =  INTERVAL  to create an output for every specified hour_interval.  But as you note this would be only for one radiation setting you specified for that particular run.  To get what you need  in model builder you'd be running the tool multiple times in succession being sure to set the  parameters correctly .  Possibly you could take advantage of creating variable lists or using model value iterators However scripting would be the easiest way.  See python sample #2 in the Point Solar Radiation help topic for how you would set the parameters and run that within a python for loop changing respective values each time.  Also please note  the usage tips regarding transmittivity and diffuse_proportion.   For example,     val = [0.1, 0.2, 0.3] for x in val: transmittivity = x PointsSolarRadiation(inRaster, inPntFC, outFeatures, "", latitude, skySize, timeConfig, dayInterval, hourInterval, "INTERVAL", zFactor, "FROM_DEM", calcDirections, zenithDivisions, azimuthDivisions,"STANDARD_OVERCAST_SKY", diffuseProp, transmittivity, outDirectRad, outDiffuseRad, outDirectDur)     Good luck.  ... View more

Excellent...  Just as a note I ran a test creating Random Points (n=1,000,000)  with FGDB feature class and a single Raster of the continental USA successfully in ArcGIS Pro.  Obviously its not the same data (or system) but the test did run successfully rather quickly (< 1 min). ... View more

@E_haus  I am happy to hear that you have had some successes.  You are still having issues with larger data? Are you able to share your code, the input shapefile and a single raster?  I can have a look on my end?  You can attach up to 9MB size here otherwise direct message me and I can provide you a location to copy data.  Thank you, -Ryan ... View more

Hi @WalkerPowell  Just wondering if you included 3D digital surface model to include the rooftops  for your analysis? For example taken from Lidar? You say the pixels are scatted,  could you provide a screen grab for reference of your data.  I can see variation along any given rooftop for many reasons.   Technically a pole, tree , air conditioner (if burned into the dsm) or adjacent buildings could partially shade an area of the roof part of the time and cause some variation during the analysis (also accounting for edges).  With that I assume you are you trying to identify are the ideal areas for all (possible) rooftops in your study area. Correct? 1. If you have the building footprints you could clip the radiation raster.  You could use CON or RECLASSIFY tools to help identify higher "areas" to start to identify areas of interest.  However you will still have many one cell or scattered pixels. 2. You could calculate the mean insolation  for rooftops using ZonalStatistics and your rooftops as zones. This wouldn't provide you with the best contiguous area to place your panel on any particular rooftop but will give you an better way to classify rooftops. On the scale you are looking at (entire town). Once you identify potential rooftops the location, size and angle of your PV systems starts to become a different question and there are many different ways to deal with optimizing panel placement. Hope this helps a little.  Or let me know if you have any additional questions. ... View more

You could also look into the SA Reclassify/Slice and  Zonal Geometry tools You could do this in a model builder environement or if wanting to script, Python samples are available on the respective tool help documents. Good luck. ... View more

@Deepa  Sounds like you are trying to "cookie cut" a section of the raster using your county feature class boundary. You are correct that you can use Extract By Mask . Additionally you could use Clip Raster and specify the bounding geometry. Note that in this case the output raster is a minimum bounding rectangle of the mask (county features) by default. You can optionally control the analysis extent using the extent environment setting. Every cell location in a raster has a value assigned to it . When information is unavailable for a cell location, the location will be assigned as NoData.   Note that NoData and 0 are not the same—0 is a valid numerical value You can convert to ascii using the Raster to Ascii conversion tool. Good luck ... View more

Thank you.  The inputs you sent ran with recent version of ArcGIS Pro 2.9.2 without issue.   However, I noticed you are using input table (x,y) for the locations.  I don't have a 10.7.1 readily available to but can try with the latest ArcGIS 10.8.2 current release of ArcMap to verify. Perhaps you can convert the table to an input feature class instead.  You should be able to create features using the make-xy-event-layer tool using the spatial reference of your input DEM. Note: The output point feature layer created by this tool is temporary and will not persist after the session ends. You can export this event layer to a feature class on disk using the Copy Features, Feature to Point, or Feature Class to Feature Class tool. let me know if the feature input work for you... ... View more