git clone 'https://github.com/filonenko-mikhail/cl-ewkb.git'
cl-ewkb is a geospatial library, based on cl-wkb, that implements the OGC Well-Known Binary geographic geometry data model with PostGIS 3d, 4d extensions, and provides WKB and EWKB encoding and decoding functionality. cl-wkb author is J.P. Larocue. Library also provide module cl-wkb, which has CLOS-based API like http://enroutesystems.com/software/cl-wkb/ with PostGIS extensions.
Depends on: ieee-floats flexi-streams
PostGIS extension WKB is described in postgis-1.5/doc/ZMSgeoms.txt All extensions are supported: 3dz, 3dm, 4d, embedded SRID.
point-primitive x y pointz-primitive x y z pointm-primitive x y m pointzm-primitive x y z m point-primtive structs contain coordinates (2d, 3d, 4d)
linear-ring points-primitive linear-ring struct contains array of point-primitive
geometry type srid geometry parent struct for other types and also contains “GEOMETRYCOLLECTION”
point type srid point-primitive line-string type srid points-primitive polygon type srid linear-rings multi-point type srid points multi-line-string type srid line-strings multi-polygon type srid polygons geometry-collection type srid geometries
(decode octets) Decode from wkb sequence (decode-from stream) Decode from wkb stream
Functions for encoding:
(encode object endianness) Encode object to vector
(encode-to object stream endianness) Encode object to stream with endianness
Getting data from postgresql and decode it:
(decode (caar (postmodern:query (:select (:ST_AsEWKB “SRID=4326;LINESTRING(0 0 1 2, 1 1 2 3, 2 2 3 4)”)))))
#(CL-EWKB::LINEAR-RING #(#(CL-EWKB::POINT-PRIMITIVE 0.0d0 0.0d0 CL-EWKB::POINTZ-PRIMITIVE 1.0d0 CL-EWKB::POINTZM-PRIMITIVE 2.0d0) #(CL-EWKB::POINT-PRIMITIVE 1.0d0 1.0d0 CL-EWKB::POINTZ-PRIMITIVE 2.0d0 CL-EWKB::POINTZM-PRIMITIVE 3.0d0) #(CL-EWKB::POINT-PRIMITIVE 2.0d0 2.0d0 CL-EWKB::POINTZ-PRIMITIVE 3.0d0 CL-EWKB::POINTZM-PRIMITIVE 4.0d0))) CL-EWKB::LINESTRING)
Drawing data with opengl (using cl-opengl):
(defun draw-point-primitive (point) (cond ((point-primitive-p point) (gl:vertex (point-primitive-x point) (point-primitive-y point))) ((pointz-primitive-p point) (gl:vertex (point-primitive-x point) (point-primitive-y point) (pointz-primitive-z point))) ((pointm-primitive-p point) (gl:vertex (point-primitive-x point) (point-primitive-y point) 0.0 (pointm-primitive-m point))) ((pointzm-primitive-p point) (gl:vertex (point-primitive-x point) (point-primitive-y point) (pointz-primitive-z point) (pointzm-primitive-m point)))))
(defun draw-point (point) (gl:with-primitives :points (draw-point-primitive (point-primitive point))))
(defun draw-points (line-string) (gl:with-primitives :line-strip (map 'nil (lambda (point) (draw-point-primitive point)) (line-string-points-primitive line-string))))
General drawing function
(defun draw-gisobject (object) (cond ((point-p object) (draw-point object)) ((line-string-p object) (draw-points object))))