/*! * \brief Write raster range file (floating-point) * * Write the floating point range file <tt>f_range</tt>. This file is * written in binary using XDR format. If there is no defined min/max * in <em>range</em>, an empty <tt>f_range</tt> file is created. * * \param name map name * \param range pointer to FPRange which holds fp range info */ void Rast_write_fp_range(const char *name, const struct FPRange *range) { int fd; char xdr_buf[2][XDR_DOUBLE_NBYTES]; Rast_init(); fd = G_open_new_misc("cell_misc", "f_range", name); if (fd < 0) { G_remove_misc("cell_misc", "f_range", name); G_fatal_error(_("Unable to write range file for <%s>"), name); } /* if range hasn't been updated, write empty file meaning Nulls */ if (range->first_time) { close(fd); return; } G_xdr_put_double(xdr_buf[0], &range->min); G_xdr_put_double(xdr_buf[1], &range->max); if (write(fd, xdr_buf, sizeof(xdr_buf)) != sizeof(xdr_buf)) { G_remove_misc("cell_misc", "f_range", name); G_fatal_error(_("Unable to write range file for <%s>"), name); } close(fd); }
int main(int argc, char *argv[]) { char *name, *mapset; char rname[GNAME_MAX], rmapset[GMAPSET_MAX]; char path[GPATH_MAX]; int row, col, null_fd; unsigned char *null_bits; RASTER_MAP_TYPE map_type; int change_null = 0, create, remove, only_int, only_fp, only_null; int is_reclass; struct GModule *module; struct { struct Option *map; struct Option *setnull; struct Option *null; } parms; struct { struct Flag *f; struct Flag *n; struct Flag *i; struct Flag *c; struct Flag *r; } flags; G_gisinit(argv[0]); module = G_define_module(); module->keywords = _("raster, null data"); module->description = _("Manages NULL-values of given raster map."); parms.map = G_define_standard_option(G_OPT_R_MAP); parms.map->description = _("Name of raster map for which to edit null file"); parms.setnull = G_define_option(); parms.setnull->key = "setnull"; parms.setnull->key_desc = "val[-val]"; parms.setnull->type = TYPE_STRING; parms.setnull->required = NO; parms.setnull->multiple = YES; parms.setnull->description = _("List of cell values to be set to NULL"); parms.setnull->guisection = _("Modify"); parms.null = G_define_option(); parms.null->key = "null"; parms.null->type = TYPE_DOUBLE; parms.null->required = NO; parms.null->multiple = NO; parms.null->description = _("The value to replace the null value by"); parms.null->guisection = _("Modify"); flags.f = G_define_flag(); flags.f->key = 'f'; flags.f->description = _("Only do the work if the map is floating-point"); flags.f->guisection = _("Check"); flags.i = G_define_flag(); flags.i->key = 'i'; flags.i->description = _("Only do the work if the map is integer"); flags.i->guisection = _("Check"); flags.n = G_define_flag(); flags.n->key = 'n'; flags.n->description = _("Only do the work if the map doesn't have a NULL-value bitmap file"); flags.n->guisection = _("Check"); flags.c = G_define_flag(); flags.c->key = 'c'; flags.c->description = _("Create NULL-value bitmap file validating all data cells"); flags.r = G_define_flag(); flags.r->key = 'r'; flags.r->description = _("Remove NULL-value bitmap file"); flags.r->guisection = _("Remove"); if (G_parser(argc, argv)) exit(EXIT_FAILURE); only_int = flags.i->answer; only_fp = flags.f->answer; only_null = flags.n->answer; create = flags.c->answer; remove = flags.r->answer; name = parms.map->answer; mapset = G_find_cell2(name, ""); if (mapset == NULL) G_fatal_error(_("Raster map <%s> not found"), name); is_reclass = (G_is_reclass(name, mapset, rname, rmapset) > 0); if (is_reclass) G_fatal_error(_("Raster map <%s> is a reclass of map <%s@%s>. " "Consider to generate a copy with r.mapcalc. Exiting."), name, rname, rmapset); if (strcmp(mapset, G_mapset()) != 0) G_fatal_error(_("Raster map <%s> is not in your mapset <%s>"), name, G_mapset()); if (parms.null->answer) { if (sscanf(parms.null->answer, "%lf", &new_null) == 1) change_null = 1; else G_fatal_error(_("%s is illegal entry for null"), parms.null->answer); } map_type = G_raster_map_type(name, mapset); if (only_null && G_find_file_misc("cell_misc", "null", name, mapset)) G_fatal_error(_("Raster map <%s> already has a null bitmap file"), name); if (map_type == CELL_TYPE) { if (only_fp) G_fatal_error(_("<%s> is integer raster map (CELL)"), name); if ((double)((int)new_null) != new_null) { G_warning(_("<%s> is integer raster map (CELL). Using null=%d."), name, (int)new_null); new_null = (double)((int)new_null); } } else if (only_int) G_fatal_error(_("<%s> is floating pointing raster map"), name); parse_vallist(parms.setnull->answers, &d_mask); if (G_get_cellhd(name, mapset, &cellhd) < 0) G_fatal_error(_("Unable to read header of raster map <%s>"), G_fully_qualified_name(name, mapset)); if (create) { /* write a file of no-nulls */ null_bits = (unsigned char *)G__allocate_null_bits(cellhd.cols); /* init all cells to 0's */ for (col = 0; col < G__null_bitstream_size(cellhd.cols); col++) null_bits[col] = 0; null_fd = G_open_new_misc("cell_misc", "null", name); G_verbose_message(_("Writing new null file for raster map <%s>..."), name); for (row = 0; row < cellhd.rows; row++) { G_percent(row, cellhd.rows, 1); if (G__write_null_bits(null_fd, null_bits, row, cellhd.cols, 0) < 0) G_fatal_error(_("Error writing null row %d"), row); } G_percent(row, cellhd.rows, 1); close(null_fd); G_done_msg(_("Raster map <%s> modified."), name); exit(EXIT_SUCCESS); } if (remove) { /* write a file of no-nulls */ G_verbose_message(_("Removing null file for raster map <%s>..."), name); null_fd = G_open_new_misc("cell_misc", "null", name); G__file_name_misc(path, "cell_misc", "null", name, mapset); unlink(path); G_done_msg(_("Raster map <%s> modified."), name); exit(EXIT_SUCCESS); } process(name, mapset, change_null, map_type); exit(EXIT_SUCCESS); }