New answers tagged ogr
GDAL supports the use of SQL for all datasources and especially the SQLite dialect http://www.gdal.org/ogr_sql_sqlite.html is quite powerful. This should work for you: ogrinfo -dialect sqlite -sql "select LineNo, StationNo, count(*) from your_layer_name group by LineNo, StationNo order by 1,2" your_tab_file
The code number 102988 is not an EPSG one, but defined by ESRI. If you run gdalsrsinfo on the projection file, it reports: +proj=tmerc +lat_0=31 +lon_0=-111.9166666666667 +k=0.9999 +x_0=213360 +y_0=0 +ellps=GRS80 +units=m +no_defs You have to create a custom CRS with these parameters to work with the dataset. GDAL comes with a similar CRS named ...
You will use GDAL to do this from the GDAL command line or OSGeo4W shell. The syntax for inspecting a raster file is: gdalinfo C:\path\to\raster For more advanced use, if you are using Python/C# bindings to do this then you can deserialize the the output into a dict or object.
Just figured out a way to do this in R, inspired by the link I posted in the comments (which uses outdated functionalities but the right packages); as an example, I'll use the "Southwesternmost" counties in Michigan (shapefile here--Allegan, Berrien, Cass, Kalamazoo, St. Joseph, and Van Buren counties) So, ordering counties alphabetically, the (length>0, ...
As @Jeffrey states, readOGR from the rgdal library imports a CRS if there is one embedded in the shapefile. You can check by (example using a shp I've been playing with): proj4string(india) # from the sp package #  "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0" If this returns NA you can specify a CRS with: proj4string(india) <- ...
Convert polygons to integer raster, same cell size and extent as value raster, snap raster=v.raster. Use rastertonumyarray on both. Create dictionary. Go through both arrays using integer as key of dictionary and updating list of values from main raster. When done go through dictionary, pick associated list, sort, etc. I've done similar thing on 5000×6000 ...
I do not know why you get an error but compare your results with mine. Download first this file https://hub.qgis.org/attachments/8536/Trecks.gdb.zip Unzip (it will create a directory) and run ogrinfo ogrinfo Trecks.gdb Had to open data source read-only. INFO: Open of `trecks.gdb' using driver `OpenFileGDB' successful. 1: Venedigertreck_3D (3D Multi ...
Try opening it in qgis and see if you can access the objects as described. You choose vector>Directory>click on souce type and choose OpenfileGDB or alternatively try to convert it to other formats using ogr. If it fails then upgrade gdal. If qgis can open this then you might be having two versions of gdal
Top 50 recent answers are included