/*! \file lib/raster/get_cellhd.c \brief Raster library - Read raster map header (C) 2001-2020 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. \author Original author CERL */ #include #include #include #include #include /*! \brief Read the raster header The raster header for the raster map name in the specified mapset is read into the cellhd structure. If there is an error reading the raster header file, G_fatal_error() is called. Cell header files may contain either grid cell header information or reclass information. If it is a reclass file, it will specify the map and mapset names of the actual grid cell file being reclassed. Rast_get_cellhd(), upon reading reclass information will go read the cell header information for the referenced file. Only one reference is allowed. \param name name of map \param mapset mapset that map belongs to \param[out] cellhd structure to hold cell header info \return void */ void Rast_get_cellhd(const char *name, const char *mapset, struct Cell_head *cellhd) { FILE *fp; int is_reclass; char real_name[GNAME_MAX], real_mapset[GMAPSET_MAX]; const char *detail; /* is_reclass = Rast_is_reclass (name, mapset, real_name, real_mapset); if (is_reclass < 0) { sprintf (buf,"Can't read header file for [%s in %s]\n", name, mapset); tail = buf + strlen(buf); strcpy (tail, "It is a reclass file, but with an invalid format"); G_warning(buf); return -1; } */ is_reclass = (Rast_is_reclass(name, mapset, real_name, real_mapset) > 0); if (is_reclass) { fp = G_fopen_old("cellhd", real_name, real_mapset); if (!fp) { detail = !G_find_raster(real_name, real_mapset) ? _("However, that raster map is missing." " Perhaps, it was deleted by mistake.") : _("However, header file of that raster map can't be" " opened. It seems that it was corrupted after" " creating the reclass raster map."); G_fatal_error(_("Unable to read header file for raster map <%s@%s>. " "It is a reclass of raster map <%s@%s>. %s"), name, mapset, real_name, real_mapset, detail); } } else { fp = G_fopen_old("cellhd", name, mapset); if (!fp) G_fatal_error(_("Unable to open header file for raster map <%s@%s>." " It seems that some previous step failed and" " created an incomplete raster map."), name, mapset); } G__read_Cell_head(fp, cellhd, 1); fclose(fp); }