我试图从我的C程序中使用某些libgeos(https://github.com/libgeos/geos)函数。下面是我的代码。
#include <stdio.h>
#include <geos_c.h>
int main() {
GEOSContextHandle_t ctx = GEOS_init_r();
GEOSGeometry *inputGeom = GEOSGeomFromWKT("MULTIPOLYGON (((30 20, 45 40, 10 40, 30 20)), ((15 5, 40 10, 10 20, 5 10, 15 5)))");
printf("%s\n", GEOSGeomToWKT(inputGeom));
GEOS_finish_r(ctx);
return 0;
}
但是问题是运行二进制文件时出现分段错误。下面是我用来编译程序的gcc命令。
gcc -g -Wall geos.c -o geos -lgeos_c
请在下面的gdb输出中查看。
(gdb) run
Starting program: /home/srimal/proj/src/geos/geos
warning: the debug information found in "/lib64/ld-2.27.so" does not match "/lib64/ld-linux-x86-64.so.2" (CRC mismatch).
Program received signal SIGSEGV, Segmentation fault.
__strlen_avx2 () at ../sysdeps/x86_64/multiarch/strlen-avx2.S:62
62 ../sysdeps/x86_64/multiarch/strlen-avx2.S: No such file or directory.
如果有人可以帮助我解决此问题,我将非常感激。
答案 0 :(得分:0)
解决了问题。我应该使用initGEOS
和finishGEOS
函数。以下是经过一些改进的正确答案。
#include <stdio.h>
#include <stdlib.h>
#include <stdarg.h>
#include <geos_c.h>
void notice(const char *fmt, ...) {
va_list ap;
fprintf( stdout, "NOTICE: ");
va_start (ap, fmt);
vfprintf( stdout, fmt, ap);
va_end(ap);
fprintf( stdout, "\n" );
}
void log_and_exit(const char *fmt, ...) {
va_list ap;
fprintf( stdout, "ERROR: ");
va_start (ap, fmt);
vfprintf( stdout, fmt, ap);
va_end(ap);
fprintf( stdout, "\n" );
exit(1);
}
int main() {
initGEOS(notice, log_and_exit);
GEOSGeometry *inputGeom = GEOSGeomFromWKT("POINT (30 10)");
if (inputGeom == NULL) {
fprintf(stderr, "NULL GEOM");
return 0;
}
printf("%s\n", GEOSGeomToWKT(inputGeom));
finishGEOS();
return 0;
}