title: Slice Large GeoTiff into 1 Degree Tiles
body: <p>I have a bunch of large GeoTiff files (1.4GB, 4 latitude degrees x 8 longitude degrees).</p>
  
  <p>I need to slice each one up into 1 degree lat x 1 degree long tiles (with just a little bit of overage in each).</p>
  
  <p>Here is the gdalinfo output for one of the GeoTiff files:</p>
  
  <pre><code>Driver: GTiff/GeoTIFF
  Files: index.tif
  Size is 38401, 19201
  Coordinate System is:
  GEOGCS["WGS 84",
      DATUM["WGS_1984",
          SPHEROID["WGS 84",6378137,298.257223563,
              AUTHORITY["EPSG","7030"]],
          AUTHORITY["EPSG","6326"]],
      PRIMEM["Greenwich",0],
      UNIT["degree",0.0174532925199433],
      AUTHORITY["EPSG","4326"]]
  Origin = (-112.000104166666674,56.000104166666667)
  Pixel Size = (0.000208333333333,-0.000208333333333)
  Metadata:
    AREA_OR_POINT=Area
  Image Structure Metadata:
    INTERLEAVE=BAND
  Corner Coordinates:
  Upper Left  (-112.0001042,  56.0001042) (112d 0' 0.38"W, 56d 0' 0.37"N)
  Lower Left  (-112.0001042,  51.9998958) (112d 0' 0.38"W, 51d59'59.62"N)
  Upper Right (-103.9998958,  56.0001042) (103d59'59.62"W, 56d 0' 0.37"N)
  Lower Right (-103.9998958,  51.9998958) (103d59'59.62"W, 51d59'59.62"N)
  Center      (-108.0000000,  54.0000000) (108d 0' 0.00"W, 54d 0' 0.00"N)
  Band 1 Block=38401x1 Type=Int16, ColorInterp=Gray
    NoData Value=-32767
  </code></pre>
  
  <p>I want everything to be the same except the size.  So each one needs to be sliced into 32 tiles, each 1 degree x 1 degree.</p>
  
  <p>For instance one would cover:<br>
  -112.0001042, 56.0001042 to -110.9998958, 54.9998958 (essentially -112, 56 to -111, 55).</p>
  
  <p>I see utilities like gdal_retile.py and gdal_grid, but I'm not making progress.  What is the right tool / command line options for this?</p>