2,668 reputation
723
bio website github.com/jlaura
location Flagstaff, AZ
age 32
visits member for 3 years, 6 months
seen Oct 17 at 19:24

I am a geography Ph.D. student at Arizona State University, with interests in FOSS, big data, spatial algorithm implementation in python, visualization, the application of cartographic techniques / spatial analysis to planetary data (Moon, Mars, Titan, etc.), and remote sensing.

My primary research focuses are in the geocomputation domain.


Aug
6
comment ArcPy Add Graph Window
@Luke Thanks for the link, but I am looking for a programmatic solution. This is a good, short term work around.
Aug
6
comment ArcPy Add Graph Window
@PolyGeo 10.2 - poor link in the question. I saw that question, but it appears that that solution places a static graph. I am trying to get a dynamic graph.
Aug
1
comment how to create graphs while running python with ArcGIS 10
instead of import pylab you can use import matplotlib.pyplot as plt. Then the calls will be plt.plot(), etc.
Jul
28
comment ArcPy Multiprocessing 000714 - Working in Debug Mode
@Luke Sounds good. I think what is happening: User installs 64-bit python (via ESRI geoprocessing or otherwise). Windows is setting .py files to open using the more recent installation of python. When background geoprocessing occurs, Arc* is opening the most recent python. So, a repair fixes the issue. So does changing the 'open with' property of the .py file.
Jul
28
comment ArcPy Multiprocessing 000714 - Working in Debug Mode
@Luke Was the reinstallation of the 32-bit Python via a python installer or an Arc installer? I will give a reinstall a try. The only change I made to Arc's python was the installation of SciPy, but perhaps that has caused some issue.
Jul
24
comment ArcGIS 10.2 with MatPlotLib Non-Blocking Plot
@EvilGenius Unfortunate to hear that that is the necessary route. Not surprising, but unfortunate.
Jul
7
comment How to store Large Raster Images of Equipotential Surface Generation Using Monte Carlo Simulation?
While my research focus has shifted significantly, I can now answer this if it of interest to other users?
May
9
comment ArcPy in_memory\ vanishing from ToC
updated - can you update your answer with the datatypes you are passing in? It amazes me that this works in 10.2 and fails in 10.1.
May
9
comment ArcPy in_memory\ vanishing from ToC
Using: res = arcpy.PolygonToRaster_conversion(roi,abundance,rasterized,"MAXIMUM_AREA","",cell‌​size) had no effect. Replacing "MAXIMUM_AREA" with "" also had no effect.
Apr
25
comment create polygon adjacency matrix using python
How about using PySAL? This functionality already exists. You want Rook contiguity. Is the data a lattice? pysal.org/library/weights/weights.html
Jan
28
comment ArcMap Natural Breaks Classification
'You can find the pseudocode here: resources.arcgis.com/content/…; Any insight here? The implementation issues that the post documents exist across variance based classification methods. They are not specific to Fisher-Jenks, nor is the info from Charlie more insightful that going to the original papers that document the algorithm. I appreciate the link, but it offers no more information than either the original work or a few hours spent implementing the algorithm.
Jan
27
comment ArcMap Natural Breaks Classification
Thanks for the info. The link you posted is the one I originally linked to as well. The poster specifically states that this is an ESRI version of Fisher-Jenks. I wonder if any additional information is available - the hand waving in the linked info does not speak to the specific alterations.
Jan
25
comment Unusual results for geoprocessing speed tests
@Aaron how about running the command line scripts as time python script. What is the difference between user + sys and real? If it is quite large, then the time may be I/O related.
Jan
25
comment Unusual results for geoprocessing speed tests
How about the fact that IDLE is also having to get a file pointer? The Arc script might already have the file pointer and/or have the geometries memory mapped.
Jan
4
comment SpatialLite BulkLoad and Merge
Thanks for the heads up and bummed that I missed that post. I did see references to using transactions to bulk insert. This might be the first avenue that I explore as the data is nicely portioned into shapefiles already.
Jan
4
comment SpatialLite BulkLoad and Merge
@HeyOverThere I total agree and would offer the same comment if someone was asking about running queries. I would partition the data as well as 8.2 million rows is getting up there for PostGIS in a single table. This is especially true given all the columns.
Oct
28
comment How do I draw transparent polygons with ArcObjects?
While not an arcobjects solution this GIS.SE might be of interest as it describes trying to set individual polygon alpha values.
Oct
28
comment GDAL command lines on windows
Do these run without wrapping them in a python script? What is the output of gdalinfo <some raster> via the command line?
Oct
27
comment How can I create random concave polygons for test purposes?
@whuber Edited the question to assume a shapefile - so topology is not recorded within the dataset. Do you have a reference to any of the many approaches available with the context of either QGIS or GRASS? Clearly, I can densify my polygons(as topology is not stored) and write custom code to perform alterations in my polygons such that some number of them are concave. This would be in addition to the technique I offered in the question. My goal is to avoid writing such code and leverage some built-in functionality.
Oct
26
comment How can I create random concave polygons for test purposes?
@whuber A sample dataset consisting of just convex polygons does not simulate the convex / concave reality that one sees in applications such as congressional redistricting - this discount multipart polygons and islands (in the sense that the unit has not neighbor, not just in the case of a physical island). Additionally, a generated dataset provides the means for a priori knowledge of the optimal solution. Can you suggest any references to the perturbation method? It sounds as if this will simply introduce slivers at common edges (fixed endpoints with each edge defined by only two points).