/**
\brief Open existing vector map
Map->name and Map->mapset must be set before.
\param Map poiter to vector map
\param update non-zero for write mode, otherwise read-only
\return 0 success
\return -1 error
*/
int V1_open_old_nat(struct Map_info *Map, int update)
{
char buf[1000];
G_debug(1, "V1_open_old_nat(): name = %s mapset = %s", Map->name,
Map->mapset);
sprintf(buf, "%s/%s", GRASS_VECT_DIRECTORY, Map->name);
dig_file_init(&(Map->dig_fp));
if (update)
Map->dig_fp.file = G_fopen_modify(buf, GRASS_VECT_COOR_ELEMENT);
else
Map->dig_fp.file =
G_fopen_old(buf, GRASS_VECT_COOR_ELEMENT, Map->mapset);
if (Map->dig_fp.file == NULL)
return -1;
if (!(dig__read_head(Map)))
return (-1);
check_coor(Map);
/* set conversion matrices */
dig_init_portable(&(Map->head.port), Map->head.port.byte_order);
/* load to memory */
if (!update)
dig_file_load(&(Map->dig_fp));
return (0);
}
/*!
\brief Save topology file for vector map
\param Map pointer to Map_info structure
\return 1 on success
\return 0 on error
*/
int Vect_save_topo(struct Map_info *Map)
{
struct Plus_head *plus;
char path[GPATH_MAX];
struct gvfile fp;
G_debug(1, "Vect_save_topo()");
/* write out all the accumulated info to the plus file */
plus = &(Map->plus);
dig_file_init(&fp);
Vect__get_path(path, Map);
fp.file = G_fopen_new(path, GV_TOPO_ELEMENT);
if (fp.file == NULL) {
G_warning(_("Unable to create topo file for vector map <%s>"), Map->name);
return 0;
}
/* set portable info */
dig_init_portable(&(plus->port), dig__byte_order_out());
if (0 > dig_write_plus_file(&fp, plus)) {
G_warning(_("Error writing out topo file"));
return 0;
}
fclose(fp.file);
return 1;
}
/*!
\brief Save category index to binary file (cidx)
\param Map pointer to Map_info structure
\return 0 on success
\return 1 on error
*/
int Vect_cidx_save(struct Map_info *Map)
{
struct Plus_head *plus;
char *path;
struct gvfile fp;
G_debug(2, "Vect_cidx_save()");
check_status(Map);
plus = &(Map->plus);
dig_file_init(&fp);
path = Vect__get_path(Map);
fp.file = G_fopen_new(path, GV_CIDX_ELEMENT);
G_free(path);
if (fp.file == NULL) {
G_warning(_("Unable to create category index file for vector map <%s>"),
Vect_get_name(Map));
return 1;
}
/* set portable info */
dig_init_portable(&(plus->cidx_port), dig__byte_order_out());
if (0 > dig_write_cidx(&fp, plus)) {
G_warning(_("Error writing out category index file"));
return 1;
}
fclose(fp.file);
return 0;
}
/**
\brief Open/Create new vector map.
\param Map pointer to vector map
\param name map name
\param with_z 2D or 3D (unused?)
\return 0 success
\return -1 error
*/
int V1_open_new_nat(struct Map_info *Map, const char *name, int with_z)
{
char buf[1000];
struct stat info;
G_debug(1, "V1_open_new_nat(): name = %s", name);
sprintf(buf, "%s/%s", GRASS_VECT_DIRECTORY, name);
/* Set the 'coor' file version */
Map->head.Version_Major = GV_COOR_VER_MAJOR;
Map->head.Version_Minor = GV_COOR_VER_MINOR;
Map->head.Back_Major = GV_COOR_EARLIEST_MAJOR;
Map->head.Back_Minor = GV_COOR_EARLIEST_MINOR;
/* TODO open better */
dig_file_init(&(Map->dig_fp));
Map->dig_fp.file = G_fopen_new(buf, GRASS_VECT_COOR_ELEMENT);
if (Map->dig_fp.file == NULL)
return (-1);
fclose(Map->dig_fp.file);
dig_file_init(&(Map->dig_fp));
Map->dig_fp.file = G_fopen_modify(buf, GRASS_VECT_COOR_ELEMENT);
if (Map->dig_fp.file == NULL)
return (-1);
/* check to see if dig_plus file exists and if so, remove it */
G__file_name(name_buf, buf, GV_TOPO_ELEMENT, G_mapset());
if (stat(name_buf, &info) == 0) /* file exists? */
unlink(name_buf);
G__file_name(name_buf, buf, GRASS_VECT_COOR_ELEMENT, G_mapset());
Map->head.size = 0;
Map->head.head_size = GV_COOR_HEAD_SIZE;
Vect__write_head(Map);
/* set conversion matrices */
dig_init_portable(&(Map->head.port), dig__byte_order_out());
if (!(dig__write_head(Map)))
return (-1);
return 0;
}
/*!
\brief Save spatial index file for vector map
\param Map vector map
\return 1 on success
\return 0 on error
*/
int Vect_save_sidx(struct Map_info *Map)
{
struct Plus_head *plus;
char file_path[GPATH_MAX];
G_debug(1, "Vect_save_spatial_index()");
plus = &(Map->plus);
if (!plus->Spidx_built) {
G_warning(_("Spatial index not available, can not be saved"));
return 0;
}
/* new or update mode ? */
if (plus->Spidx_new == TRUE) {
/* write out rtrees to sidx file */
Vect__get_element_path(file_path, Map, GV_SIDX_ELEMENT);
G_debug(1, "Open sidx: %s", file_path);
dig_file_init(&(plus->spidx_fp));
plus->spidx_fp.file = fopen(file_path, "w+");
if (plus->spidx_fp.file == NULL) {
G_warning(_("Unable to create spatial index file for vector map <%s>"),
Vect_get_name(Map));
return 0;
}
/* set portable info */
dig_init_portable(&(plus->spidx_port), dig__byte_order_out());
if (0 > dig_Wr_spidx(&(plus->spidx_fp), plus)) {
G_warning(_("Error writing out spatial index file"));
return 0;
}
Map->plus.Spidx_new = FALSE;
}
fclose(Map->plus.spidx_fp.file);
Map->plus.Spidx_built = FALSE;
return 1;
}
/*!
\brief Read category index from cidx file if exists
\param Map pointer to Map_info structure
\param head_only read only header of the file
\return 0 on success
\return 1 if file does not exist
\return -1 error, file exists but cannot be read
*/
int Vect_cidx_open(struct Map_info *Map, int head_only)
{
int ret;
char file_path[GPATH_MAX], *path;
struct gvfile fp;
struct Plus_head *Plus;
G_debug(2, "Vect_cidx_open(): name = %s mapset= %s", Map->name,
Map->mapset);
Plus = &(Map->plus);
path = Vect__get_path(Map);
G_file_name(file_path, path, GV_CIDX_ELEMENT, Map->mapset);
if (access(file_path, F_OK) != 0) { /* does not exist */
G_free(path);
return 1;
}
dig_file_init(&fp);
fp.file = G_fopen_old(path, GV_CIDX_ELEMENT, Map->mapset);
G_free(path);
if (fp.file == NULL) { /* category index file is not available */
G_warning(_("Unable to open category index file for vector map <%s>"),
Vect_get_full_name(Map));
return -1;
}
/* load category index to memory */
ret = dig_read_cidx(&fp, Plus, head_only);
fclose(fp.file);
if (ret == 1) {
G_debug(3, "Cannot read cidx");
return -1;
}
return 0;
}
/*!
\brief Save feature index file for vector map
\param Map pointer to Map_info structure
\param offset pointer to Format_info_offset struct
(see Format_info_ogr and Format_info_pg struct for implementation issues)
\return 1 on success
\return 0 on error
*/
int Vect_save_fidx(struct Map_info *Map,
struct Format_info_offset *offset)
{
#ifdef HAVE_OGR
char fname[GPATH_MAX], elem[GPATH_MAX];
char buf[5];
long length;
struct gvfile fp;
struct Port_info port;
if (strcmp(Map->mapset, G_mapset()) != 0 ||
Map->support_updated == FALSE ||
Map->plus.built != GV_BUILD_ALL)
return 1;
length = 9;
sprintf(elem, "%s/%s", GV_DIRECTORY, Map->name);
G_file_name(fname, elem, GV_FIDX_ELEMENT, Map->mapset);
G_debug(4, "Open fidx: %s", fname);
dig_file_init(&fp);
fp.file = fopen(fname, "w");
if (fp.file == NULL) {
G_warning(_("Unable to open fidx file for write <%s>"), fname);
return 0;
}
dig_init_portable(&port, dig__byte_order_out());
dig_set_cur_port(&port);
/* Header */
/* bytes 1 - 5 */
buf[0] = 5;
buf[1] = 0;
buf[2] = 5;
buf[3] = 0;
buf[4] = (char)dig__byte_order_out();
if (0 >= dig__fwrite_port_C(buf, 5, &fp))
return 0;
/* bytes 6 - 9 : header size */
if (0 >= dig__fwrite_port_L(&length, 1, &fp))
return 0;
/* Body */
/* number of records */
if (0 >= dig__fwrite_port_I(&(offset->array_num), 1, &fp))
return 0;
/* offsets */
if (0 >= dig__fwrite_port_I(offset->array,
offset->array_num, &fp))
return 0;
G_debug(3, "Vect_save_fidx(): offset_num = %d", offset->array_num);
fclose(fp.file);
return 1;
#else
G_fatal_error(_("GRASS is not compiled with OGR support"));
return 0;
#endif
}
/* read old 3.0 or 4.0 dig file into array
returns number of elements read into array
or -1 on error */
int read_dig(FILE * Digin, struct Map_info *Mapout,
struct Line **plines, int endian, int att)
{
char buf[100];
struct dig_head In_head;
int lalloc, line = 0, type, portable = 1;
int npoints = 0, nlines = 0, nbounds = 0;
int ndpoints = 0, ndlines = 0, ndbounds = 0, nunknown = 0;
struct Line *lines;
struct line_pnts *nline;
struct line_cats *cat_out;
double dbuf;
int ibuf;
long lbuf;
struct gvfile gvf;
dig_file_init(&gvf);
gvf.file = Digin;
Vect__init_head(Mapout);
/* set conversion matrices */
dig_init_portable(&(In_head.port), endian);
/* Version 3 dig files were not portable and some version 4
* files may be also non portable */
G_message(_("Reading dig file..."));
/* read and copy head */
dig_fseek(&gvf, 0L, SEEK_SET); /* set to beginning */
if (0 >= dig__fread_port_C(buf, DIG4_ORGAN_LEN, &gvf))
return -1;
buf[DIG4_ORGAN_LEN - 1] = '\0';
Vect_set_organization(Mapout, buf);
if (0 >= dig__fread_port_C(buf, DIG4_DATE_LEN, &gvf))
return -1;
buf[DIG4_DATE_LEN - 1] = '\0';
Vect_set_date(Mapout, buf);
if (0 >= dig__fread_port_C(buf, DIG4_YOUR_NAME_LEN, &gvf))
return -1;
buf[DIG4_YOUR_NAME_LEN - 1] = '\0';
Vect_set_person(Mapout, buf);
if (0 >= dig__fread_port_C(buf, DIG4_MAP_NAME_LEN, &gvf))
return -1;
buf[DIG4_MAP_NAME_LEN - 1] = '\0';
Vect_set_map_name(Mapout, buf);
if (0 >= dig__fread_port_C(buf, DIG4_SOURCE_DATE_LEN, &gvf))
return -1;
buf[DIG4_SOURCE_DATE_LEN - 1] = '\0';
Vect_set_map_date(Mapout, buf);
if (0 >= dig__fread_port_C(buf, DIG4_LINE_3_LEN, &gvf))
return -1;
buf[DIG4_LINE_3_LEN - 1] = '\0';
Vect_set_comment(Mapout, buf);
if (0 >= dig__fread_port_C(buf, VERS_4_DATA_SIZE, &gvf))
return -1;
if (buf[0] != '%' || buf[1] != '%') { /* Version3.0 */
In_head.Version_Major = 3;
portable = 0; /* input vector is not portable format */
G_message(_("Input file is version 3."));
}
else {
In_head.Version_Major = 4;
G_message(_("Input file is version 4."));
/* determine if in portable format or not */
if (buf[6] == 1 && (~buf[6] & 0xff) == (buf[7] & 0xff)) { /* portable ? */
portable = 1; /* input vector is portable format */
}
else {
portable = 0; /* input vector is not portable format */
}
}
if (portable == 1) {
G_message(_("Input file is portable."));
}
else {
G_warning(_("Input file is not portable. "
"We will attempt to convert anyway but conversion may fail. "
"Please read manual for detail information."));
}
/* set Cur_Head because it is used by dig__*_convert()
called by dig__fread_port_*() */
dig_set_cur_port(&(In_head.port));
if (0 >= dig__fread_port_L(&lbuf, 1, &gvf))
return -1;
Vect_set_scale(Mapout, (int)lbuf);
if (0 >= dig__fread_port_I(&ibuf, 1, &gvf))
return -1;
Vect_set_zone(Mapout, ibuf);
if (0 >= dig__fread_port_D(&dbuf, 1, &gvf))
return -1; /* W */
if (0 >= dig__fread_port_D(&dbuf, 1, &gvf))
return -1; /* E */
if (0 >= dig__fread_port_D(&dbuf, 1, &gvf))
return -1; /* S */
if (0 >= dig__fread_port_D(&dbuf, 1, &gvf))
return -1; /* N */
if (0 >= dig__fread_port_D(&dbuf, 1, &gvf))
return -1;
Vect_set_thresh(Mapout, dbuf);
/* reading dig file body (elements) */
nline = Vect_new_line_struct();
cat_out = Vect_new_cats_struct();
lalloc = 0;
lines = NULL;
line = 0;
while (1) {
type = read_line(&gvf, nline);
G_debug(3, "read line = %d, type = %d", line, type);
if (type == -2)
break; /* EOF */
switch (type) {
case GV_POINT:
npoints++;
break;
case GV_LINE:
nlines++;
break;
case GV_BOUNDARY:
nbounds++;
break;
case 0: /* dead */
break;
default:
nunknown++;
break;
}
if (!(type & (GV_POINT | GV_LINE | GV_BOUNDARY)))
continue;
if ((type & GV_BOUNDARY) || !att) {
Vect_write_line(Mapout, type, nline, cat_out);
/* reset In_head */
dig_set_cur_port(&(In_head.port));
}
else { /* GV_POINT or GV_LINE */
if (line >= lalloc) {
lalloc += 10000;
lines =
(struct Line *)G_realloc(lines,
lalloc * sizeof(struct Line));
}
lines[line].type = type;
lines[line].n_points = nline->n_points;
lines[line].cat = -1;
lines[line].x =
(double *)G_malloc(nline->n_points * sizeof(double));
lines[line].y =
(double *)G_malloc(nline->n_points * sizeof(double));
memcpy((void *)lines[line].x, (void *)nline->x,
nline->n_points * sizeof(double));
memcpy((void *)lines[line].y, (void *)nline->y,
nline->n_points * sizeof(double));
line++;
}
}
if (att) {
G_message(_("[%d] points read to memory"), npoints);
G_message(_("[%d] lines read to memory"), nlines);
}
else {
G_message(_("[%d] points read and written to output"), npoints);
G_message(_("[%d] lines read and written to output"), nlines);
}
G_message(_("[%d] area boundaries read and written to output"), nbounds);
G_message(_("[%d] dead points skipped"), ndpoints);
G_message(_("[%d] dead lines skipped"), ndlines);
G_message(_("[%d] dead area boundaries skipped"), ndbounds);
G_message(_("[%d] elements of unknown type skipped"), nunknown);
G_message(_("[%d] elements read to memory"), line);
*plines = lines;
return (line);
}
/*!
\brief Open feature index file
\param[in,out] Map pointer to Map_info struct
\param[out] offset pointer to Format_info_offset (OGR or PG)
\return 0 on success
\return -1 on error
*/
int Vect_open_fidx(struct Map_info *Map, struct Format_info_offset *offset)
{
char elem[GPATH_MAX];
char buf[5]; /* used for format version */
long length;
int Version_Major, Version_Minor, Back_Major, Back_Minor, byte_order;
struct gvfile fp;
struct Port_info port;
G_debug(1, "Vect_open_fidx(): name = %s mapset = %s format = %d",
Map->name, Map->mapset, Map->format);
sprintf(elem, "%s/%s", GV_DIRECTORY, Map->name);
dig_file_init(&fp);
fp.file = G_fopen_old(elem, GV_FIDX_ELEMENT, Map->mapset);
if (fp.file == NULL) {
G_debug(1, "unable to open fidx file for vector map <%s>",
Vect_get_full_name(Map));
return -1;
}
/* Header */
if (0 >= dig__fread_port_C(buf, 5, &fp))
return -1;
Version_Major = buf[0];
Version_Minor = buf[1];
Back_Major = buf[2];
Back_Minor = buf[3];
byte_order = buf[4];
/* check version numbers */
if (Version_Major > 5 || Version_Minor > 0) {
if (Back_Major > 5 || Back_Minor > 0) {
G_fatal_error(_("Feature index format version %d.%d is not supported by this release."
" Try to rebuild topology or upgrade GRASS."),
Version_Major, Version_Minor);
return -1;
}
G_warning(_("Your GRASS version does not fully support feature index format %d.%d of the vector."
" Consider to rebuild topology or upgrade GRASS."),
Version_Major, Version_Minor);
}
dig_init_portable(&port, byte_order);
dig_set_cur_port(&port);
/* Body */
/* bytes 6 - 9 : header size */
if (0 >= dig__fread_port_L(&length, 1, &fp))
return -1;
G_debug(4, " header size %ld", length);
G_fseek(fp.file, length, SEEK_SET);
/* number of records */
if (0 >= dig__fread_port_I(&(offset->array_num), 1, &fp))
return -1;
/* alloc space */
offset->array = (int *) G_malloc(offset->array_num * sizeof(int));
offset->array_alloc = offset->array_num;
/* offsets */
if (0 >= dig__fread_port_I(offset->array,
offset->array_num, &fp))
return -1;
fclose(fp.file);
G_debug(3, "%d records read from fidx", offset->array_num);
return 0;
}
