1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586 |
- /****************************************************************************
- *
- * MODULE: r.coin
- *
- * AUTHOR(S): Michael O'Shea - CERL
- * Michael Shapiro - CERL
- *
- * PURPOSE: Calculates the coincidence of two raster map layers.
- *
- * COPYRIGHT: (C) 2006 by the GRASS Development Team
- *
- * This program is free software under the GNU General Public
- * License (>=v2). Read the file COPYING that comes with GRASS
- * for details.
- *
- ***************************************************************************/
- #include <grass/raster.h>
- #include "coin.h"
- int print_coin_hdr(int Conformat)
- {
- char unit_type[20];
- const char *mapset, *location;
- char north[30], south[30], east[30], west[30];
- mapset = G_mapset();
- location = G_location();
- switch (Conformat) {
- case 'a':
- sprintf(unit_type, "acres");
- break;
- case 'h':
- sprintf(unit_type, "hectares");
- break;
- case 'k':
- sprintf(unit_type, "square kilometers");
- break;
- case 'm':
- sprintf(unit_type, "square miles");
- break;
- case 'p':
- sprintf(unit_type, "percent cover");
- break;
- case 'x':
- sprintf(unit_type, "percent of cols");
- break;
- case 'y':
- sprintf(unit_type, "percent of rows");
- break;
- default:
- sprintf(unit_type, "cells");
- break;
- }
- G_format_northing(window.north, north, window.proj);
- G_format_northing(window.south, south, window.proj);
- G_format_easting(window.east, east, window.proj);
- G_format_easting(window.west, west, window.proj);
- fprintf(dumpfile, "\n");
- fprintf(dumpfile, "+%78.78s+\n", midline);
- fprintf(dumpfile, "|%24.24sCOINCIDENCE TABULATION REPORT%25.25s|\n",
- fill, fill);
- fprintf(dumpfile, "|%78.78s|\n", midline);
- fprintf(dumpfile, "| Location: %-16.14sMapset: %-17.15sDate: %-20.20s|\n",
- location, mapset, G_date());
- fprintf(dumpfile, "|%78.78s|\n", fill);
- fprintf(dumpfile, "| Layer 1: %-15.15s-- %-50.49s|\n", map1name, title1);
- fprintf(dumpfile, "| Layer 2: %-15.15s-- %-50.49s|\n", map2name, title2);
- fprintf(dumpfile, "| Mask: %-68.68s|\n", Rast_mask_info());
- fprintf(dumpfile, "|%78.78s|\n", fill);
- fprintf(dumpfile, "| Units: %-68.68s|\n", unit_type);
- fprintf(dumpfile, "|%78.78s|\n", midline);
- fprintf(dumpfile, "| Window:%22.22sNorth: %-10s%31.31s|\n",
- fill, north, fill);
- fprintf(dumpfile, "|%14.14sWest: %-9s%19.19sEast: %-9s%15.15s|\n",
- fill, west, fill, east, fill);
- fprintf(dumpfile, "|%30.30sSouth: %-10s%31.31s|\n", fill, south, fill);
- fprintf(dumpfile, "+%78.78s+\n\n", midline);
- return 0;
- }
|