Yet again, working with that Ohio DEM dataset, I need some reduced resolution versions of the 2.5 foot DEMs. At a grand total of ~600 DEMs for a given county, it was time to batch, this time in Windows Command prompt. The neat thing is the variable substitution in the for-in-do command in the command prompt, where, e.g. for a file variable %f, the file … Continue reading Batch processing resampling of DEMs
I hadn’t seen that PerryGeo’s work has been rolled into GDAL: http://www.gdal.org/gdaldem.html Nice work folks. I especially like the addition of topological position index. I’ll be playing with that… . Continue reading DEM Processing with GDAL Utilities
Alright, so let’s begin again with some LiDAR processing in PostGIS. I really liked this discussion of multipoint conversions of LiDAR and other point datasets, but we won’t address making things more efficient here (yet). Besides, these aren’t just points– we do want to retain their attribution… . Now, any regular database user will recognize my nubeness reading this post. If you are that person, … Continue reading LiDAR processing and analysis in PostGIS (I hope…).
In general, I like the results from the hillshade. One problem that I forgot to take into account– hillshade algorithms typically leave a 1-pixel border along the edge where calculations cannot take place… . For this Ohio dataset, that means that every 5000 feet or so (the x and y dimension of our DEM chunks), we have a 2-pixel wide line, resulting in a subtle … Continue reading Hillshade (cont.)
Using the same LiDAR DEM from which we generated the contours, we can create hillshade tifs using http://www.perrygeo.net/wordpress/?p=7. It compiles easily on a mac, probably even easier on a Linux machine following his directions. Then another simple loop: #!/bin/bash x=0 for f in $( ls *.txt); do x=`expr $x + 1` echo $x $f hillshade $f $f.tif done I’d like this all in a single … Continue reading Hillshade using PerryGeo’s ported GRASS utility
Suppose you had a list of addresses coded like this: 00014 SOUTH ST That you’d prefer coded like this: 14 SOUTH ST If you were addicted to modifying things programmatically in Excel, like I quite guiltily do, you could do this in MS Excel. =(MID(A1,1,FIND(” “,A1)-1))*1 & MID(A1,FIND(” “,A1),LEN(A1)) Everything before the ampersand (&) uses the =mid() function to grab everything up until the first … Continue reading Spreadsheet Fun– Excel
Well, I made some mistakes in the last post, not the least of which is I used the wrong flag for creating an attribute field with elevation. What follows is a little more sophisticated. It takes us from a series of DEM tiles from which I want 2-foot and 5-foot contours (using gdal_contour), and then dumps those shapefiles into PostgreSQL using shp2pgsql. First we prep … Continue reading GDAL Contours (cont.)
Just another quick vignette. From the Ohio Statewide Imagery Program (OSIP) there is a 1-meter DEM for all of Ohio. To get contours from this dataset, one approach is to use GDAL tools, i.e. gdal_contours. As I’m working on a Mac today, I used Kyng Chaos pre-compiled Frameworks: http://www.kyngchaos.com/software:frameworks Then I needed to update my path variable in the BASH shell: export PATH=/Library/Frameworks/GDAL.framework/Programs:$PATH Now we … Continue reading GDAL Contours
At home, I work on a Mac, so when I want to do work here, my scripting language changes. In this case, I prefer BASH as my control script. I downloaded libLAS and did the standard Unix ./configure, make, sudo make install dance in the decompressed libLAS directory ./configure make sudo make install Now I have some tools to manipulate LiDAR las files at home. … Continue reading LiDAR data in BASH using libLAS
Ok, so the INSERT statements were too numerous for inputing the LiDAR point dataset (about a billion points… .) They kept crashing the postgres daemon. So, I used copy from a CSV file: COPY base.cuy_lidar_all FROM ‘c:/path/cuy_lidar_ground_veg.insert’ WITH CSV Keep your fingers crossed… . Continue reading COPY command psql– loading large LiDAR Point Dataset
Smathermather's Weblog 