static void test_patch_range_compression_lazperf() { int i; int npts = 20; PCPOINTLIST *pl; PCPATCH *pa; PCPATCH *par; char *str; pl = pc_pointlist_make(npts); for ( i = 0; i < npts; i++ ) { PCPOINT *pt = pc_point_make(simpleschema); pc_point_set_double_by_name(pt, "X", i); pc_point_set_double_by_name(pt, "Y", i); pc_point_set_double_by_name(pt, "Z", i * 0.1); pc_point_set_double_by_name(pt, "Intensity", 100 - i); pc_pointlist_add_point(pl, pt); } pa = (PCPATCH*)pc_patch_lazperf_from_pointlist(pl); par = pc_patch_range(pa, 16, 4); str = pc_patch_to_string(par); CU_ASSERT_STRING_EQUAL(str, "{\"pcid\":0,\"pts\":[[15,15,1.5,85],[16,16,1.6,84],[17,17,1.7,83],[18,18,1.8,82]]}"); pcfree(str); pc_patch_free(par); pc_patch_free(pa); pc_pointlist_free(pl); }
static void test_patch_filter_lazperf_zero_point() { PCPOINT *pt; int i; int npts = 5; PCPOINTLIST *pl; PCPATCH_LAZPERF *pal; PCPATCH *pa; // build a list of points pl = pc_pointlist_make(npts); for ( i = 0; i < npts; i++ ) { pt = pc_point_make(simpleschema); pc_point_set_double_by_name(pt, "x", i*2.0); pc_point_set_double_by_name(pt, "y", i*1.9); pc_point_set_double_by_name(pt, "Z", i*0.34); pc_point_set_double_by_name(pt, "intensity", 10); pc_pointlist_add_point(pl, pt); } // build patch lazperf pal = pc_patch_lazperf_from_pointlist(pl); // filter with a resulting patch of 0 point(s) pa = pc_patch_filter((PCPATCH*) pal, 0, PC_BETWEEN, 0.0, 0.0); CU_ASSERT_EQUAL(pa->npoints, 0); pc_patch_free((PCPATCH*) pal); pc_patch_free((PCPATCH*) pa); pc_pointlist_free(pl); }
static void test_patch_pointn_laz_compression() { // 00 endian (big) // 00000000 pcid // 00000000 compression // 00000003 npoints // 0000000800000003000000050006 pt1 (XYZi) // 0000000200000003000000040008 pt2 (XYZi) // 0000000200000003000000040009 pt3 (XYZi) char *hexbuf = "00000000000000000000000003000000080000000300000005000600000002000000030000000400080000000200000003000000040009"; size_t hexsize = strlen(hexbuf); uint8_t *wkb = bytes_from_hexbytes(hexbuf, hexsize); char *str; PCPATCH *pa = pc_patch_from_wkb(simpleschema, wkb, hexsize/2); PCPOINTLIST *li = pc_pointlist_from_patch(pa); PCPATCH_LAZPERF* paz = pc_patch_lazperf_from_pointlist(li); PCPOINT *pt = pc_patch_pointn((PCPATCH*) paz, 2); CU_ASSERT(pt != NULL); str = pc_point_to_string(pt); CU_ASSERT_STRING_EQUAL(str, "{\"pcid\":0,\"pt\":[0.02,0.03,0.04,8]}"); pc_patch_free((PCPATCH *)paz); pc_point_free(pt); pcfree(str); pcfree(wkb); pc_patch_free(pa); pc_pointlist_free(li); }
static void test_patch_lazperf() { PCPOINT *pt; int i; int npts = 400; PCPOINTLIST *pl; PCPATCH_LAZPERF *pal; PCPATCH_UNCOMPRESSED *paul, *pauref; // build a list of points pl = pc_pointlist_make(npts); for ( i = 0; i < npts; i++ ) { pt = pc_point_make(simpleschema); pc_point_set_double_by_name(pt, "x", i*2.0); pc_point_set_double_by_name(pt, "y", i*1.9); pc_point_set_double_by_name(pt, "Z", i*0.34); pc_point_set_double_by_name(pt, "intensity", 10); pc_pointlist_add_point(pl, pt); } // compress the list in a lazperf patch pal = pc_patch_lazperf_from_pointlist( pl ); // get an uncompressed patch from the lazperf patch paul = pc_patch_uncompressed_from_lazperf( pal ); // get an uncompressed patch directly from the pointlist pauref = pc_patch_uncompressed_from_pointlist( pl ); // test the number of points CU_ASSERT_EQUAL(pal->npoints, pauref->npoints); CU_ASSERT_EQUAL(paul->npoints, pauref->npoints); // test bounds CU_ASSERT_DOUBLE_EQUAL(pal->bounds.xmax, pauref->bounds.xmax, 0.0001); CU_ASSERT_DOUBLE_EQUAL(paul->bounds.ymax, pauref->bounds.ymax, 0.000001); // test type CU_ASSERT_EQUAL(pal->type, PC_LAZPERF); CU_ASSERT_EQUAL(paul->type, pauref->type); // test readonly CU_ASSERT_EQUAL(pauref->readonly, paul->readonly); CU_ASSERT_EQUAL(pauref->readonly, pal->readonly); // test datasize CU_ASSERT_EQUAL(paul->datasize, pauref->datasize); // free pc_pointlist_free(pl); pc_patch_free( (PCPATCH*) pal ); pc_patch_free((PCPATCH*) paul); pc_patch_free((PCPATCH*) pauref); }
static void test_wkb_lazperf() { PCPOINT *pt; int i; int npts = 400; PCPOINTLIST *pl; PCPATCH_LAZPERF *pal1, *pal2; PCPATCH_UNCOMPRESSED *pau; uint8_t *wkb1, *wkb2; size_t wkbsize; // build a list of points pl = pc_pointlist_make(npts); for ( i = 0; i < npts; i++ ) { pt = pc_point_make(simpleschema); pc_point_set_double_by_name(pt, "x", i*2.0); pc_point_set_double_by_name(pt, "y", i*1.9); pc_point_set_double_by_name(pt, "Z", i*0.34); pc_point_set_double_by_name(pt, "intensity", 10); pc_pointlist_add_point(pl, pt); } // build patch lazperf pal1 = pc_patch_lazperf_from_pointlist(pl); // get the corresponding wkb wkb1 = pc_patch_lazperf_to_wkb(pal1, &wkbsize); // rebuild a lazperf patch thanks to the wkb pal2 = (PCPATCH_LAZPERF*) pc_patch_lazperf_from_wkb( pal1->schema, wkb1, wkbsize); // get the wkb reference pau = pc_patch_uncompressed_from_pointlist(pl); wkb2 = pc_patch_uncompressed_to_wkb( pau, &wkbsize ); // compare wkb CU_ASSERT_STRING_EQUAL(wkb1, wkb2); // free pc_patch_free((PCPATCH*) pal1); pc_patch_free((PCPATCH*) pal2); pc_patch_free((PCPATCH*) pau); pc_pointlist_free(pl); pcfree(wkb1); pcfree(wkb2); }
static void test_pointlist_lazperf() { PCPOINT *pt; int i; int npts = 400; PCPOINTLIST *pl1, *pl2; PCPATCH_LAZPERF *pch1; PCPATCH_UNCOMPRESSED *pa1, *pa2; char *wkt1, *wkt2; // build a list of points pl1 = pc_pointlist_make(npts); for ( i = 0; i < npts; i++ ) { pt = pc_point_make(simpleschema); pc_point_set_double_by_name(pt, "x", i*2.0); pc_point_set_double_by_name(pt, "y", i*1.9); pc_point_set_double_by_name(pt, "Z", i*0.34); pc_point_set_double_by_name(pt, "intensity", 10); pc_pointlist_add_point(pl1, pt); } // compress the list in a lazperf patch pch1 = pc_patch_lazperf_from_pointlist( pl1 ); // decompress the lazperf patch in a pointlist pl2 = pc_pointlist_from_lazperf(pch1); // test that the string representation of pointlist is equal pa1 = pc_patch_uncompressed_from_pointlist( pl1 ); pa2 = pc_patch_uncompressed_from_lazperf( pch1 ); wkt1 = pc_patch_uncompressed_to_string(pa1); wkt2 = pc_patch_uncompressed_to_string(pa2); CU_ASSERT_STRING_EQUAL(wkt1, wkt2); pc_patch_free((PCPATCH*) pch1 ); pc_patch_free((PCPATCH*) pa1); pc_patch_free((PCPATCH*) pa2); pc_pointlist_free(pl1); pc_pointlist_free(pl2); pcfree(wkt1); pcfree(wkt2); }
static void test_to_string_lazperf() { PCPOINT *pt; int i; int npts = 400; PCPOINTLIST *pl; PCPATCH_LAZPERF *pal; PCPATCH_UNCOMPRESSED *pau; char *str1, *str2; // build a list of points pl = pc_pointlist_make(npts); for ( i = 0; i < npts; i++ ) { pt = pc_point_make(simpleschema); pc_point_set_double_by_name(pt, "x", i*2.0); pc_point_set_double_by_name(pt, "y", i*1.9); pc_point_set_double_by_name(pt, "Z", i*0.34); pc_point_set_double_by_name(pt, "intensity", 10); pc_pointlist_add_point(pl, pt); } // build patch pau = pc_patch_uncompressed_from_pointlist(pl); pal = pc_patch_lazperf_from_pointlist(pl); // get string str1 = pc_patch_uncompressed_to_string(pau); str2 = pc_patch_lazperf_to_string(pal); // compare CU_ASSERT_STRING_EQUAL(str1, str2); // free pc_patch_free((PCPATCH*) pal); pc_patch_free((PCPATCH*) pau); pc_pointlist_free(pl); pcfree(str1); pcfree(str2); }
static void test_patch_set_schema_compression_lazperf() { // init data PCPATCH_LAZPERF *pal; PCPATCH *pat; PCPOINTLIST *pl; PCPOINT *pt; char *str; int i; int npts = 4; // build a patch pl = pc_pointlist_make(npts); for ( i = npts; i >= 0; i-- ) { pt = pc_point_make(simpleschema); pc_point_set_double_by_name(pt, "X", i * 0.1); pc_point_set_double_by_name(pt, "Y", i * 0.2); pc_point_set_double_by_name(pt, "Z", i * 0.3); pc_point_set_double_by_name(pt, "Intensity", 10); pc_pointlist_add_point(pl, pt); } pal = pc_patch_lazperf_from_pointlist(pl); // assign a valid schema to the patch pat = pc_patch_set_schema((PCPATCH*) pal, simpleschema_nointensity, 0.0); CU_ASSERT(pat != NULL); str = pc_patch_to_string(pat); CU_ASSERT_STRING_EQUAL(str, "{\"pcid\":0,\"pts\":[[0.4,0.8,1.2],[0.3,0.6,0.9],[0.2,0.4,0.6],[0.1,0.2,0.3],[0,0,0]]}"); pc_patch_free(pat); pcfree(str); pc_patch_free((PCPATCH*) pal); pc_pointlist_free(pl); }