Latest Contributions by IanMurray

Have you tried Simplify Line yet? http://desktop.arcgis.com/en/arcmap/10.3/tools/cartography-toolbox/simplify-line.htm http://desktop.arcgis.com/en/arcmap/10.3/tools/cartography-toolbox/how-simplify-line-works.htm#GUID-0321EFE4-CD00-4D24-9814-454D980CAD67 ... View more

Conversely if you do not have a data interoperability license, you can can use LASzip tool free to unzip .laz files. For a small walkthrough: https://gisjames.wordpress.com/2017/04/13/converting-laz-to-las-files/ Link to LASZip: https://laszip.org/ Additional help on StackExchange: https://gis.stackexchange.com/questions/188170/converting-lidar-files-from-laz-to-las-format ... View more

Chris, I don't use ArcEarth, but I have an idea why it might be reading no spatial reference when you bring in the file after using define projection. GeoTIFF is a file format that natively stores its spatial information in a header.  However, the define projection tool doesn't change the header for the the .tif, it just adds a world file with the spatial information (a .tfw).  If ArcEarth is expecting the the spatial information in the header and not a world file and reads no spatial information in the header then it could return no spatial reference error.  The Project Raster tool creates a new raster and when doing so writes the spatial data to the header(At least on a test I did), which is why that one may be coming in correctly. I couldn't find any documentation on how ArcEarth reads spatial information from rasters, but this is what came to mind for me. ... View more

Looking at your screenshot, it looks like the input it failed on was appeared to be a very large numbers of spaces or tabs.  You look to have a condition to screen for a NULL value and for a single blank space, but not for multiple blank spaces or tabs.  Perhaps your first elif should be something like this: elif prop_location.strip() == '':    return None This would check for any amount of possible whitespace that could be in a field and if there are no other characters, return None. ... View more

True but its not even necessarily an accuracy problem, but an outright error problem in some cases.  I've never found anywhere that adequately explains the source of the terrain data in Google Earth, there has been alot of conjecture about its a combination of SRTM with other data in some places, but the fact I often find artifacts from trees would tend to lead me away from that.  There are areas that you can find much higher quality terrain information in Google Earth where they have incorporated LiDAR or other high resolution data, but if you are outside of those areas, how useful the data is becomes highly questionable to me. ... View more

I would consider looking at the Terrain Elevation Layer from ESRI Living Atlas, they regularly incorporate higher resolution data into their Terrain Layer for use.  Obviously you are working in some remote/not well sampled areas, but there should be at least some 24m meter data(as opposed to the 30m SRTM) from the AirBus data they incorporated lately. Coverage Area overview for their various data sources can be viewed here ... View more

Personally I am very leery of the terrain information within Google Earth.  There seems to be some issues with tree affecting the elevations in some areas within Google Earth, as well as some asymptotic like features where there should be none when drawing elevation profiles of a line.  ... View more

https://knowledge.autodesk.com/support/autocad/learn-explore/caas/CloudHelp/cloudhelp/2016/ENU/AutoCAD-Core/files/GUID-9CEED44E-E8D3-43FE-B4CF-D838379C7D78-htm.html https://knowledge.autodesk.com/support/autocad/learn-explore/caas/CloudHelp/cloudhelp/2016/ENU/AutoCAD-Core/files/GUID-14B82899-9C2D-4A34-8A02-49319C6FB38C-htm.html I don't know for sure but the above help seems to indicate it can be done in AutoCAD itself. ... View more

I regularly have co-workers generate CAD drawings with a coordinate system in AutoCAD/Civil 3D, the entire drawing can be brought directly into the correct location in ArcGIS as long as the coordinate system in ArcGIS is the same as the one used in AutoCad.  The data still will not have a projection however, it will however be based on the same coordinate system and come into the proper locaiton.  I usually convert the CAD drawings to a shapefile or GDB by the type (Annotation, Point, Polyline, Polygon, Multipatch) and then define projection for the output files. You will need to check with whomever created the CAD Drawing to see if it had a coordinate system and if so what is was.  ... View more

He has nested cursors 3 deep in from line 127 onward, the one open for that loop, another one on the subsequent with statements on lines 131, 138, 145, 152, and 159, and then within each of those loops call a function opening a third cursor.  That is more than just a small performance hit, see https://gis.stackexchange.com/questions/220107/how-to-improve-performance-of-nested-search-cursors.  ... View more

Yikes that is some nested cursor madness, you should never run one cursor inside of another cursor.  No wonder that gets bogged down.  Please see Richard Fairhurst's blog, it should have some examples that will help out streamlining your processing of data using dictionaries instead of keeping multiple cursors open at a time.  https://community.esri.com/blogs/richard_fairhurst/2014/11/08/turbo-charging-data-manipulation-with-python-cursors-and-dictionaries  Edit: Tagging Richard Fairhurst Further Edit: Didn't see you had troubleshot a similar problem with Richard https://community.esri.com/thread/119381‌ ... View more

Could you give some examples on how you named the files?  I feel like if they are named in a uniform easy to parse manner it would fairly simple to make separate lists for years and months and weeks and then use them as input for Cell Statistics. ... View more

Do your rasters have the same extent and cell size? ... View more

The class break values returns a list, so it would be fairly easy to set the break values to a variable and then access them via a loop or use index values. import arcpy mxd = arcpy.mapping.MapDocument("CURRENT") df = arcpy.mapping.ListDataFrames(mxd, "Lakes")[0] lyr = arcpy.mapping.ListLayers(mxd, "lakene.tif", df)[0] if lyr.symbologyType == "RASTER_CLASSIFIED":   values = lyr.symbology.classBreakValues   for val in values:     print val #checking value     #do something else with the values here.     ‍‍‍‍‍‍‍‍‍‍  You could also just read them to a variable and then do something with them outside of the loop that you had set-up there, the important thing is to reading them into a variable so you have access to them for later in your script.  As the help file says the values for classBreakValues can be read or written, so you can check what they currently are and also have the ability to write new values to the layer symbology in the example you posted.  You want to just read what they currently are. I think the better question is how are you going to set the class break values for every raster prior to doing the processing you want.  A raster file by itself has no symbology at all, only the values stored in it.  A layer file either in a map document or on a hard disk would have to be created for you to actually have break values showing the different symbolization of the raster values.  Your script would have to work with layer files that you have already created, or you can have a single layer file set-up with the symbology and break types and numbers you are looking for and apply them to a raster layer created from each raster and then creating the break values using the reclassify() method seen in the help. Are you planning to use the same symbology method and number of classes each time, or will it vary raster to raster?  How many raster layer break values are you needing to make and read? ... View more

You need to have consistent spacing when indenting for and if statements.  Your for statement is indented 2 spaces, while your if statement is indented 4 spaces from the for statement when it should be 2.  Just make sure they are consistent, and that you aren't mixing spaces and tabs.  https://www.python.org/dev/peps/pep-0008/ for columns in range (xmin, xmax):   pointval= arcpy.GetCellValue_management(sorder,"%s %s" % (columns++,ymin), "1") #indented 2 spaces   if (pointval== 1, pointval==2, pointval==3):     count++ #indented 2 more spaces‍‍‍‍‍‍‍‍ ... View more