Reputation
6,062
Top tag
Next privilege 10,000 Rep.
Access moderator tools
Badges
1 14 32
Newest
 Excavator
Impact
~386k people reached

Jul
18
reviewed Approve How to get the geometry type of an empty PostGIS table?
Jul
18
reviewed Approve PostGIS ST_Contains query return empty
Jun
30
reviewed Approve oil-gas tag wiki excerpt
Jun
28
comment Output of Add XY Coordinates?
WGS 84 describes a datum (earth's shape), WGS 84 UTM Zone 30N is a UTM projection with a WGS 84 datum. Your X value is the UTM easting and the Y the UTM northing. You need to use the Define Projection tool to establish its existing coordinate reference system PRIOR to using the Project tool to get it to a GCS WGS 84 datum (aka unprojected data in decimal degrees with a WGS 84 datum)
Jun
28
reviewed Approve Leaflet mouse wheel zoom only after click on map
Jun
28
reviewed Approve OpenLayers Lat Lon Conversion
Jun
26
answered How do I compare two features to test if they are identical?
Jun
26
reviewed Approve Does the USDA SSURGO database have a web-based API that I can query?
Jun
26
comment ArcPy reading a raster in a moving window?
James pardon my ignorance... I guess I am unclear as to your purpose. Ss @Evil_Genius pointed out, once the NumPy array is created from the input raster, there is little need for arcpy or arcmap at all unless you are trying to recreate some spatial analyst function that can't be replicated simply. Your mention of moving windows and that brought to mind neighborhood statistics (focal, block, filter etc) in arcmap. I find arcpy.sa functions are convenience functions that can be replicated in NumPy and NumPy is a strong alternative if one can't afford a Spatial Analyst license.
Jun
25
comment ArcPy reading a raster in a moving window?
Strided_tricks or the SciPy implementation does moving/sliding window operations smoothly and easily. Some of the useful bits are clocked under image analysis, but the outcomes are the same. NumPy, regardless of whether it is on its own or cloaked in SciPy will be fast if you can vectorize your operations as indicated
Jun
22
comment After execute arcgis script, what happen with the variables?
go to your python editor and type ... locals().keys() ... at the prompt. Are variable names still there? If they are, you can reuse them.
Jun
22
comment Compute annual mean of overlapping cells using daily rasters
ps from numpy impot * is spelled incorrectly and highly discouraged... use import numpy as np in that way you always know that np.mean is different than mean in python
Jun
22
comment Compute annual mean of overlapping cells using daily rasters
What is the raster size and shape? Have you managed to get a list of rasters by year using the slicing method? It might be worthwhile creating an numpy (np) array of zeros (see np.zeros) the exact size and shape of the expected input arrays. Convert the arrays one at a time and add them incrementally using the += inplace operator to conserve memory and add them in sequence to the zeros array. Sum/N = avg or check out incremental averaging (assuming you have no nodata values). If your data have nodata values check out np.nansum and np.nanmean as alternatives
Jun
21
comment Why is there a lock when using arcpy.DeleteField_management in arcgis 10.0?
is the word 'this' supposed to be a comment on the line? if so use, '# this' or remove it from the thread
Jun
20
comment Hex grid on 4326 globe (for use in leaflet)
the link and one of the links in it is interesting webpages.sou.edu/~sahrk/dgg/pubs/gdggs03.pdf should you want to stick with some distortion etc. As for covering large areas and finding a projected coordinate system that is great for all things...they don't exist. If you limited your view to smaller areas, then you could go the Google route (earth/map ... one uses a projected system, the other geographic)
Jun
18
answered Hex grid on 4326 globe (for use in leaflet)
Jun
8
awarded  Excavator
Jun
8
revised Exporting feature class into multiple feature classes based on field values using ArcGIS for Desktop?
new version, due to changes in Esri's script locations
Jun
7
comment When building fieldmappings, how do I handle the SHAPE field?
In the arcpy.Describe section, there is a shapeFieldName property
Jun
2
comment Calculate Lat/Long in decimal degrees
Do it once manually from arctoolbox, then go to the Results window, and copy the syntax from there ... assuming it works. In any event dump the @decimaldegrees as suggested KHimba