Take the 2-minute tour ×
Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It's 100% free, no registration required.

I have a need to create regular points with a coordinate every meter in a county. I've tried several methods including Hawthes Tools in ArcMap9.3, ET GeoWizards, and the standard Toolbox. I have encountered reoccurring instances where all of the processes have froze up.I am assuming that these programs can't handle that data generation. I tried this in QGIS as well. I need help ASAP.

share|improve this question
2  
The average US county would contain three billion (3E9) such points, requiring 48 GB to store as double precision floats. Some counties would require an order of magnitude more. That's an awful lot; that many points wouldn't be distinguishable on any map. Could you maybe back up a little and explain what you need these points for? Perhaps there's a better way. –  whuber Aug 12 '11 at 19:29
2  
@whuber I have the storage capacity (around 200 GB). I'm not really sure what the need is, but my boss said to do it and I really want to keep my job. –  PatrickW Aug 12 '11 at 19:48
1  
OK, everyone can appreciate that. Note that storage is the least of your worries. You already have seen what datasets of this size can do even to industrial-strength GISes like ArcGIS: they might require decades to finish their calculations (I'm not kidding). Perhaps your boss would appreciate a good analysis of tradeoffs between computational resource use (especially computing time and your time) and analytical precision. (It's difficult to imagine any county-level analysis that truly needs one meter precision.) –  whuber Aug 12 '11 at 19:55

2 Answers 2

up vote 7 down vote accepted

Try this:

  1. Convert your county polygons to raster using the ArcGIS Conversion Toolbox > To Raster > Polygon to Raster. Set the cell size to 1m.
  2. Use spatial analyst's Sample Tool (Spatial Analyst toolbox > Extraction tools > sample) to create points for every cell.
  3. Extract only the points that intersect your county polygon.

If you continue to run into problems, try splitting the county polygon into smaller areas and using this same process.

share|improve this answer
2  
+1: if this has to be done, rasters are the way to do it. Moreover, the inherent runlength encoding used by ESRI's raster format will compress even this ~10^10 cell raster into a few megabytes; only the vector output will be large (and tedious to generate). I do worry about what the output is intended for, though... –  whuber Aug 12 '11 at 19:59
1  
PatrickS: Since we have the same name, Patrick the answerer will prepend any comments with "PatrickS:". –  Patrick Aug 12 '11 at 19:59
    
@Patrick Thanks. Didnt catch that. –  PatrickW Aug 12 '11 at 20:16
    
PatrickS: No worries! I'm new too. –  Patrick Aug 12 '11 at 20:21
1  
@PatrickW Welcome to our site! I look forward to seeing lots of great questions and answers with "Patrick"s in them :-). –  whuber Aug 12 '11 at 21:36

Personally, I would get the grid of the MBR of the country: create a 2D array based on the min max coords of the MBR, working how many items would fill in that grid with 1 meter intervals and then loop[ through, creating points on the fly, rounding it off with a delete outside of the county.

So you'll know the Min / Max X, so you'll know how many metres will fill in that grid, so you have your grid of X values. Do the same for the Y and you have that grid.

The following is pseudo code

for y in range(0, rows):
    for x in range(0, cols):
        currentPoint = (y, x) 
        CreatePoint

I reckon that's about 20 lines of python with the erase outside?

share|improve this answer

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.