I needed to recursively document a bunch of Shapefiles and geotiffs in a hierarchy of directories (soon file geodatabses). I can’t use Arc*anything so I’m using ogrinfo and gdalinfo in many directories. I started using xargs, but ran into weirdness beyond white characters, and involving bash string manipulation (I’m running on Cygwin in windows…). So I switched approaches and came up with this little ditty.
#!/bin/sh echo "********************** START Shapefiles ***********************" echo "" for x in $(find ./* -type f -name \*.shp); do base=${x##*/} echo "--- Found the Shapefile " $base "and the base layer of " ${base%.*} /bin/GDAL/./ogrinfo.exe -ro -so -al -fields=YES -geom=SUMMARY $x echo "" done echo "" echo "********************** END Shapefiles ***********************" echo "********************** START TIFFs ***********************" echo "" for y in $(find ./* -type f -name \*.tif); do base=${y##*/} echo "--- Found the TIFF file " $base "and the base layer of " ${base%.*} /bin/GDAL/./gdalinfo.exe -approx_stats -mm -noct -proj4 $y echo "" done echo "" echo "********************** END TIFFs ***********************"
which outputs things like this for a TIFF
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 | --- Found the TIFF file melcd_2004_imperviousness.tif and the base layer of melcd_2004_imperviousness Driver: GTiff/GeoTIFF Files: ./Land_Characteristics/impervs/melcd_2004_imperviousness.tif ./Land_Characteristics/impervs/melcd_2004_imperviousness.tfw ./Land_Characteristics/impervs/melcd_2004_imperviousness.aux Size is 65102, 99355 Coordinate System is: PROJCS["UTM", GEOGCS["NAD83", DATUM["North_American_Datum_1983", SPHEROID["GRS 1980",6378137,298.2572221010002, AUTHORITY["EPSG","7019"]], AUTHORITY["EPSG","6269"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4269"]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",-69], PARAMETER["scale_factor",0.9996], PARAMETER["false_easting",500000], PARAMETER["false_northing",0], UNIT["meters",1], AUTHORITY["EPSG","26919"]] PROJ.4 string is: '+proj=utm +zone=19 +datum=NAD83 +units=m +no_defs ' Origin = (336631.500000000000000,5256292.500000000000000) Pixel Size = (5.000000000000000,-5.000000000000000) Metadata: AREA_OR_POINT=Point TIFFTAG_RESOLUTIONUNIT=1 (unitless) TIFFTAG_SOFTWARE=IMAGINE TIFF Support Copyright 1991 - 1999 by ERDAS, Inc. All Rights Reserved @(#)$RCSfile: etif.c $ $Revision: 1.10.1.9.1.9.2.11 $ $Date: 2004/09/15 18:42:01EDT $ TIFFTAG_XRESOLUTION=0.2 TIFFTAG_YRESOLUTION=0.2 Image Structure Metadata: COMPRESSION=CCITTRLE INTERLEAVE=BAND Corner Coordinates: Upper Left ( 336631.500, 5256292.500) ( 71d10' 0.27"W, 47d26'22.55"N) Lower Left ( 336631.500, 4759517.500) ( 71d 0'12.02"W, 42d58'14.82"N) Upper Right ( 662141.500, 5256292.500) ( 66d50'58.28"W, 47d26'23.65"N) Lower Right ( 662141.500, 4759517.500) ( 67d 0'42.12"W, 42d58'15.76"N) Center ( 499386.500, 5007905.000) ( 69d 0'28.13"W, 45d13'28.69"N) Band 1 Block=65102x2 Type=Byte, ColorInterp=Palette Description = Band_1 Min=0.000 Max=1.000 Computed Min/Max=0.000,1.000 Minimum=0.000, Maximum=1.000, Mean=0.992, StdDev=0.091 Metadata: LAYER_TYPE=athematic STATISTICS_MAXIMUM=1 STATISTICS_MEAN=0.99165311862261 STATISTICS_MEDIAN=0 STATISTICS_MINIMUM=0 STATISTICS_MODE=1.9020308565413e-230 STATISTICS_STDDEV=0.090979178661602 Image Structure Metadata: NBITS=1 Color Table (RGB with 2 entries) <GDALRasterAttributeTable /> |
--- Found the TIFF file melcd_2004_imperviousness.tif and the base layer of melcd_2004_imperviousness Driver: GTiff/GeoTIFF Files: ./Land_Characteristics/impervs/melcd_2004_imperviousness.tif ./Land_Characteristics/impervs/melcd_2004_imperviousness.tfw ./Land_Characteristics/impervs/melcd_2004_imperviousness.aux Size is 65102, 99355 Coordinate System is: PROJCS["UTM", GEOGCS["NAD83", DATUM["North_American_Datum_1983", SPHEROID["GRS 1980",6378137,298.2572221010002, AUTHORITY["EPSG","7019"]], AUTHORITY["EPSG","6269"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4269"]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",-69], PARAMETER["scale_factor",0.9996], PARAMETER["false_easting",500000], PARAMETER["false_northing",0], UNIT["meters",1], AUTHORITY["EPSG","26919"]] PROJ.4 string is: '+proj=utm +zone=19 +datum=NAD83 +units=m +no_defs ' Origin = (336631.500000000000000,5256292.500000000000000) Pixel Size = (5.000000000000000,-5.000000000000000) Metadata: AREA_OR_POINT=Point TIFFTAG_RESOLUTIONUNIT=1 (unitless) TIFFTAG_SOFTWARE=IMAGINE TIFF Support Copyright 1991 - 1999 by ERDAS, Inc. All Rights Reserved @(#)$RCSfile: etif.c $ $Revision: 1.10.1.9.1.9.2.11 $ $Date: 2004/09/15 18:42:01EDT $ TIFFTAG_XRESOLUTION=0.2 TIFFTAG_YRESOLUTION=0.2 Image Structure Metadata: COMPRESSION=CCITTRLE INTERLEAVE=BAND Corner Coordinates: Upper Left ( 336631.500, 5256292.500) ( 71d10' 0.27"W, 47d26'22.55"N) Lower Left ( 336631.500, 4759517.500) ( 71d 0'12.02"W, 42d58'14.82"N) Upper Right ( 662141.500, 5256292.500) ( 66d50'58.28"W, 47d26'23.65"N) Lower Right ( 662141.500, 4759517.500) ( 67d 0'42.12"W, 42d58'15.76"N) Center ( 499386.500, 5007905.000) ( 69d 0'28.13"W, 45d13'28.69"N) Band 1 Block=65102x2 Type=Byte, ColorInterp=Palette Description = Band_1 Min=0.000 Max=1.000 Computed Min/Max=0.000,1.000 Minimum=0.000, Maximum=1.000, Mean=0.992, StdDev=0.091 Metadata: LAYER_TYPE=athematic STATISTICS_MAXIMUM=1 STATISTICS_MEAN=0.99165311862261 STATISTICS_MEDIAN=0 STATISTICS_MINIMUM=0 STATISTICS_MODE=1.9020308565413e-230 STATISTICS_STDDEV=0.090979178661602 Image Structure Metadata: NBITS=1 Color Table (RGB with 2 entries) <GDALRasterAttributeTable />
and things like this for a Shapefile
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 | ********************** START Shapefiles *********************** --- Found the Shapefile acfish2.shp and the base layer of acfish2 INFO: Open of <code>./Bio_Eco_Cons/acfish2s/acfish2.shp' using driver </code>ESRI Shapefile' successful. Layer name: acfish2 Geometry: Multi Point Feature Count: 177 Extent: (356761.978694, 4771813.744319) - (653088.978694, 5004932.744319) Layer SRS WKT: PROJCS["NAD_1983_UTM_Zone_19N", GEOGCS["GCS_North_American_1983", DATUM["North_American_Datum_1983", SPHEROID["GRS_1980",6378137.0,298.257222101]], PRIMEM["Greenwich",0.0], UNIT["Degree",0.0174532925199433]], PROJECTION["Transverse_Mercator"], PARAMETER["False_Easting",500000.0], PARAMETER["False_Northing",0.0], PARAMETER["Central_Meridian",-69.0], PARAMETER["Scale_Factor",0.9996], PARAMETER["Latitude_Of_Origin",0.0], UNIT["Meter",1.0]] ACFISH2_ID: Integer (9.0) DMR_ID: String (5.0) ECCM_ID: String (5.0) YEAR: String (4.0) ECCM_SP: String (10.0) USFW_SP: String (10.0) DMR_SP: String (10.0) DMR_NAME: String (30.0) |
********************** START Shapefiles *********************** --- Found the Shapefile acfish2.shp and the base layer of acfish2 INFO: Open of <code>./Bio_Eco_Cons/acfish2s/acfish2.shp' using driver </code>ESRI Shapefile' successful. Layer name: acfish2 Geometry: Multi Point Feature Count: 177 Extent: (356761.978694, 4771813.744319) - (653088.978694, 5004932.744319) Layer SRS WKT: PROJCS["NAD_1983_UTM_Zone_19N", GEOGCS["GCS_North_American_1983", DATUM["North_American_Datum_1983", SPHEROID["GRS_1980",6378137.0,298.257222101]], PRIMEM["Greenwich",0.0], UNIT["Degree",0.0174532925199433]], PROJECTION["Transverse_Mercator"], PARAMETER["False_Easting",500000.0], PARAMETER["False_Northing",0.0], PARAMETER["Central_Meridian",-69.0], PARAMETER["Scale_Factor",0.9996], PARAMETER["Latitude_Of_Origin",0.0], UNIT["Meter",1.0]] ACFISH2_ID: Integer (9.0) DMR_ID: String (5.0) ECCM_ID: String (5.0) YEAR: String (4.0) ECCM_SP: String (10.0) USFW_SP: String (10.0) DMR_SP: String (10.0) DMR_NAME: String (30.0)
Now I just need to see about getting this working with File geodatabases!