head.c 4.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130
  1. /*
  2. ****************************************************************************
  3. *
  4. * MODULE: Vector library
  5. *
  6. * AUTHOR(S): Original author CERL, probably Dave Gerdes.
  7. * Update to GRASS 5.7 Radim Blazek.
  8. *
  9. * PURPOSE: Lower level functions for reading/writing/manipulating vectors.
  10. *
  11. * COPYRIGHT: (C) 2001 by the GRASS Development Team
  12. *
  13. * This program is free software under the GNU General Public
  14. * License (>=v2). Read the file COPYING that comes with GRASS
  15. * for details.
  16. *
  17. *****************************************************************************/
  18. #include <string.h>
  19. #include <grass/gis.h>
  20. #include <grass/Vect.h>
  21. int dig__write_head(struct Map_info *Map)
  22. {
  23. unsigned char buf[10];
  24. long length = GV_COOR_HEAD_SIZE;
  25. G_debug(1, "dig__write_head()");
  26. dig_set_cur_port(&(Map->head.port));
  27. dig_fseek(&(Map->dig_fp), 0L, 0);
  28. /* bytes 1 - 5 */
  29. buf[0] = Map->head.Version_Major;
  30. buf[1] = Map->head.Version_Minor;
  31. buf[2] = Map->head.Back_Major;
  32. buf[3] = Map->head.Back_Minor;
  33. buf[4] = Map->head.port.byte_order;
  34. if (0 >= dig__fwrite_port_C(buf, 5, &(Map->dig_fp)))
  35. return (0);
  36. /* bytes 6 - 9 : header size */
  37. if (0 >= dig__fwrite_port_L(&length, 1, &(Map->dig_fp)))
  38. return (0);
  39. /* byte 10 : dimension 2D or 3D */
  40. buf[0] = Map->head.with_z;
  41. if (0 >= dig__fwrite_port_C(buf, 1, &(Map->dig_fp)))
  42. return (0);
  43. /* bytes 11 - 14 : size of coordinate file */
  44. G_debug(1, "write coor size (%ld) to head", Map->head.size);
  45. if (0 >= dig__fwrite_port_L(&(Map->head.size), 1, &(Map->dig_fp)))
  46. return (0);
  47. G_debug(2, "coor body offset %ld", dig_ftell(&(Map->dig_fp)));
  48. return (1);
  49. }
  50. int dig__read_head(struct Map_info *Map)
  51. {
  52. unsigned char buf[10];
  53. struct Port_info port;
  54. dig_fseek(&(Map->dig_fp), 0L, 0);
  55. /* bytes 1 - 5 */
  56. if (0 >= dig__fread_port_C(buf, 5, &(Map->dig_fp)))
  57. return (0);
  58. Map->head.Version_Major = buf[0];
  59. Map->head.Version_Minor = buf[1];
  60. Map->head.Back_Major = buf[2];
  61. Map->head.Back_Minor = buf[3];
  62. Map->head.port.byte_order = buf[4];
  63. G_debug(2,
  64. "Coor header: file version %d.%d , supported from GRASS version %d.%d",
  65. Map->head.Version_Major, Map->head.Version_Minor,
  66. Map->head.Back_Major, Map->head.Back_Minor);
  67. G_debug(2, " byte order %d", Map->head.port.byte_order);
  68. /* check version numbers */
  69. if (Map->head.Version_Major > GV_COOR_VER_MAJOR ||
  70. Map->head.Version_Minor > GV_COOR_VER_MINOR) {
  71. /* The file was created by GRASS library with higher version than this one */
  72. if (Map->head.Back_Major > GV_COOR_VER_MAJOR ||
  73. Map->head.Back_Minor > GV_COOR_VER_MINOR) {
  74. /* This version of GRASS lib is lower than the oldest which can read this format */
  75. G_fatal_error
  76. ("Vector 'coor' format version %d.%d is not supported by this version of GRASS. "
  77. "Update your GRASS.", Map->head.Version_Major,
  78. Map->head.Version_Minor);
  79. return (-1);
  80. }
  81. G_warning
  82. ("Your GRASS version does not fully support vector format %d.%d."
  83. " Consider to upgrade GRASS.", Map->head.Version_Major,
  84. Map->head.Version_Minor);
  85. }
  86. dig_init_portable(&port, Map->head.port.byte_order);
  87. dig_set_cur_port(&port);
  88. /* bytes 6 - 9 : header size */
  89. if (0 >= dig__fread_port_L(&(Map->head.head_size), 1, &(Map->dig_fp)))
  90. return (0);
  91. G_debug(2, " header size %ld", Map->head.head_size);
  92. /* byte 10 : dimension 2D or 3D */
  93. if (0 >= dig__fread_port_C(buf, 1, &(Map->dig_fp)))
  94. return (0);
  95. Map->head.with_z = buf[0];
  96. G_debug(2, " with_z %d", Map->head.with_z);
  97. /* bytes 11 - 14 : size of coordinate file */
  98. if (0 >= dig__fread_port_L(&(Map->head.size), 1, &(Map->dig_fp)))
  99. return (0);
  100. G_debug(2, " coor size %ld", Map->head.size);
  101. /* Go to end of header, file may be written by new version of GRASS with longer header */
  102. dig_fseek(&(Map->dig_fp), Map->head.head_size, SEEK_SET);
  103. return (1);
  104. }