int main(int argc, char *argv[])
{
int i;
int dry = FALSE;
int verbose = FALSE;
char* file_name_in = 0;
char* file_name_out = 0;
double xyz_min[3] = {0.0, 0.0, 0.0};
double xyz_max[3] = {0.0, 0.0, 0.0};
double xyz_scale[3] = {0.01,0.01,0.01};
double xyz_offset[3] = {0.0,0.0,0.0};
unsigned int number_of_point_records = 0;
unsigned int number_of_points_by_return[8] = {0,0,0,0,0,0,0,0};
char* parse_string = "xyz";
int file_creation_day = 0;
int file_creation_year = 0;
char* system_identifier = 0;
char* generating_software = 0;
#define MAX_CHARACTERS_PER_LINE 512
char line[MAX_CHARACTERS_PER_LINE];
double xyz[3];
LASPointH point = NULL;
double gps_time;
FILE* file_in = NULL;
char* parse_less = NULL;
LASHeaderH header = NULL;
LASWriterH writer = NULL;
LASError err;
int format = LAS_FORMAT_12;
int xyz_min_quant[3] = {0, 0, 0};
int xyz_max_quant[3] = {0, 0, 0};
double xyz_min_dequant[3] = {0.0, 0.0, 0.0};
double xyz_max_dequant[3] = {0.0, 0.0, 0.0};
for (i = 1; i < argc; i++)
{
if ( strcmp(argv[i],"-h") == 0 ||
strcmp(argv[i],"--help") == 0
)
{
usage();
exit(0);
}
else if ( strcmp(argv[i],"-v") == 0 ||
strcmp(argv[i],"--verbose") == 0
)
{
verbose = TRUE;
}
else if ( strcmp(argv[i],"-d") == 0 ||
strcmp(argv[i],"-dry") == 0 ||
strcmp(argv[i],"--dry") == 0
)
{
dry = TRUE;
}
else if ( strcmp(argv[i],"--parse") == 0 ||
strcmp(argv[i],"-parse") == 0 ||
strcmp(argv[i],"-p") == 0
)
{
i++;
parse_string = argv[i];
}
else if ( strcmp(argv[i],"--scale") == 0 ||
strcmp(argv[i],"-scale") == 0 ||
strcmp(argv[i],"-s") == 0
)
{
i++;
sscanf(argv[i], "%lf", &(xyz_scale[2]));
xyz_scale[0] = xyz_scale[1] = xyz_scale[2];
}
else if ( strcmp(argv[i],"--xyz_scale") == 0 ||
strcmp(argv[i],"-xyz_scale") == 0
)
{
i++;
sscanf(argv[i], "%lf", &(xyz_scale[0]));
i++;
sscanf(argv[i], "%lf", &(xyz_scale[1]));
i++;
sscanf(argv[i], "%lf", &(xyz_scale[2]));
}
else if ( strcmp(argv[i],"--xyz_offset") == 0 ||
strcmp(argv[i],"-xyz_offset") == 0
)
{
i++;
sscanf(argv[i], "%lf", &(xyz_offset[0]));
i++;
sscanf(argv[i], "%lf", &(xyz_offset[1]));
i++;
sscanf(argv[i], "%lf", &(xyz_offset[2]));
}
else if ( strcmp(argv[i],"--input") == 0 ||
strcmp(argv[i],"-input") == 0 ||
strcmp(argv[i],"-i") == 0 ||
strcmp(argv[i],"-in") == 0
)
{
i++;
file_name_in = argv[i];
}
else if ( strcmp(argv[i],"--output") == 0 ||
strcmp(argv[i],"--out") == 0 ||
strcmp(argv[i],"-out") == 0 ||
strcmp(argv[i],"-o") == 0
)
{
i++;
file_name_out = argv[i];
}
else if ( strcmp(argv[i],"--format") == 0 ||
strcmp(argv[i],"-f") == 0 ||
strcmp(argv[i],"-format") == 0
)
{
i++;
if (strcmp(argv[i], "1.0") == 0) {
format = LAS_FORMAT_10;
}
else if (strcmp(argv[i], "1.1") == 0) {
format = LAS_FORMAT_11;
}
else if (strcmp(argv[i], "1.2") == 0) {
format = LAS_FORMAT_12;
}
else {
LASError_Print("Format must be specified as 1.0, 1.1, or 1.2");
}
}
else if ( strcmp(argv[i],"--system_identifier") == 0 ||
strcmp(argv[i],"-system_identifier") == 0 ||
strcmp(argv[i],"-s") == 0 ||
strcmp(argv[i],"-sys_id") == 0)
{
i++;
system_identifier = (char*) malloc(31 * sizeof(char));
strcpy(system_identifier, argv[i]);
}
else if ( strcmp(argv[i],"--generating_software") == 0 ||
strcmp(argv[i],"-generating_software") == 0 ||
strcmp(argv[i],"-g") == 0 ||
strcmp(argv[i],"-gen_soft") == 0)
{
i++;
generating_software = (char*) malloc(31*sizeof(char));
strcpy(generating_software, argv[i]);
}
else if ( strcmp(argv[i],"--file_creation") == 0 ||
strcmp(argv[i],"-file_creation") == 0)
{
i++;
file_creation_day = (unsigned short)atoi(argv[i]);
i++;
file_creation_year = (unsigned short)atoi(argv[i]);
}
else if ( i == argc - 2
&& file_name_in == NULL &&
file_name_out == NULL
)
{
file_name_in = argv[i];
}
else if ( i == argc - 1 &&
file_name_in == NULL &&
file_name_out == NULL
)
{
file_name_in = argv[i];
}
else if ( i == argc - 1 &&
file_name_in &&
file_name_out == NULL
)
{
file_name_out = argv[i];
}
else
{
fprintf(stderr, "ERROR: unknown argument '%s'\n",argv[i]);
usage();
exit(1);
}
}
/* create output file name if none specified and no piped output requested */
if (file_name_out == NULL && file_name_in != NULL)
{
int len = (int)strlen(file_name_in);
file_name_out = strdup(file_name_in);
while (len > 0 && file_name_out[len] != '.')
{
len--;
}
file_name_out[len] = '.';
file_name_out[len+1] = 'l';
file_name_out[len+2] = 'a';
file_name_out[len+3] = 's';
file_name_out[len+4] = '\0';
}
/* make sure that input and output are not *both* piped */
if (file_name_in == NULL && file_name_out == NULL)
{
LASError_Print("both input and output filenames are null!");
usage();
exit(1);
}
file_in = fopen(file_name_in, "r");
if (file_in == NULL)
{
LASError_Print("could not open file to read for first pass");
exit(1);
}
/* create a cheaper parse string that only looks for 'x' 'y' 'z' and 'r' */
parse_less = strdup(parse_string);
for (i = 0; i < (int)strlen(parse_string); i++)
{
if (parse_less[i] != 'x' &&
parse_less[i] != 'y' &&
parse_less[i] != 'z' &&
parse_less[i] != 'r')
{
parse_less[i] = 's';
}
}
do
{
parse_less[i] = '\0';
printf("nuking %d for %c\n", i, parse_less[i]);
i--;
} while (parse_less[i] == 's');
/* first pass to figure out the bounding box and number of returns */
if (verbose) {
fprintf(stderr,
"first pass over file '%s' with parse '%s'\n",
file_name_in,
parse_less);
}
/* read the first line */
while (fgets(line, sizeof(char) * MAX_CHARACTERS_PER_LINE, file_in))
{
point = LASPoint_Create();
if (parse(parse_less, line, xyz, point, &gps_time))
{
/* init the bounding box */
VecCopy3dv(xyz_min, xyz);
VecCopy3dv(xyz_max, xyz);
/* mark that we found the first point */
number_of_point_records = 1;
/* create return histogram */
number_of_points_by_return[LASPoint_GetReturnNumber(point)]++;
/* we can stop this loop */
break;
}
else
{
fprintf(stderr, "WARNING: cannot parse '%s' with '%s'. skipping ...\n",
line,
parse_less);
}
LASPoint_Destroy(point);
point = NULL;
}
/* did we manage to parse a line? */
if (number_of_point_records != 1)
{
fprintf(stderr, "ERROR: could not parse any lines with '%s'\n",
parse_less);
exit(1);
}
/* loop over the remaining lines */
while (fgets(line, sizeof(char) * MAX_CHARACTERS_PER_LINE, file_in))
{
point = LASPoint_Create();
if (parse(parse_less, line, xyz, point, &gps_time))
{
/* update bounding box */
VecUpdateMinMax3dv(xyz_min, xyz_max, xyz);
/* count points */
number_of_point_records++;
/* create return histogram */
number_of_points_by_return[LASPoint_GetReturnNumber(point)]++;
}
else
{
fprintf(stderr, "WARNING: cannot parse '%s' with '%s'. skipping ...\n",
line,
parse_less);
}
LASPoint_Destroy(point);
point = NULL;
}
/* output some stats */
if (verbose)
{
fprintf(stderr,
"npoints %d min %g %g %g max %g %g %g\n",
number_of_point_records,
xyz_min[0],
xyz_min[1],
xyz_min[2],
xyz_max[0],
xyz_max[1],
xyz_max[2]
);
fprintf(stderr,
"return histogram %d %d %d %d %d %d %d %d\n",
number_of_points_by_return[0],
number_of_points_by_return[1],
number_of_points_by_return[2],
number_of_points_by_return[3],
number_of_points_by_return[4],
number_of_points_by_return[5],
number_of_points_by_return[6],
number_of_points_by_return[7]
);
}
/* close the input file */
fclose(file_in);
/* compute bounding box after quantization */
for (i = 0; i < 3; i++)
{
xyz_min_quant[i] = (int)(0.5 + (xyz_min[i] - xyz_offset[i]) / xyz_scale[i]);
xyz_max_quant[i] = (int)(0.5 + (xyz_max[i] - xyz_offset[i]) / xyz_scale[i]);
}
for (i = 0; i < 3; i++)
{
xyz_min_dequant[i] = xyz_offset[i] + (xyz_min_quant[i] * xyz_scale[i]);
xyz_max_dequant[i] = xyz_offset[i] + (xyz_max_quant[i] * xyz_scale[i]);
}
#define log_xor !=0==!
/* make sure there is not sign flip */
for (i = 0; i < 3; i++)
{
if ((xyz_min[i] > 0) log_xor (xyz_min_dequant[i] > 0))
{
fprintf(stderr,
"WARNING: quantization sign flip for %s min coord %g -> %g. use offset or scale up\n",
(i ? (i == 1 ? "y" : "z") : "x"),
xyz_min[i],
xyz_min_dequant[i]
);
}
if ((xyz_max[i] > 0) log_xor (xyz_max_dequant[i] > 0))
{
fprintf(stderr,
"WARNING: quantization sign flip for %s max coord %g -> %g. use offset or scale up\n",
(i ? (i == 1 ? "y" : "z") : "x"),
xyz_max[i],
xyz_max_dequant[i]
);
}
}
#undef log_xor
/* populate the header */
header = LASHeader_Create();
if (system_identifier) LASHeader_SetSystemId(header, system_identifier);
if (generating_software) LASHeader_SetSoftwareId(header, generating_software);
LASHeader_SetCreationDOY(header, file_creation_day);
LASHeader_SetCreationYear(header, file_creation_year);
if (format == LAS_FORMAT_10) {
LASHeader_SetVersionMinor(header, 0);
} else if (format == LAS_FORMAT_11){
LASHeader_SetVersionMinor(header, 1);
} else if (format == LAS_FORMAT_12) {
LASHeader_SetVersionMinor(header, 2);
}
if (strstr(parse_string,"t") && (strstr(parse_string, "R") || strstr(parse_string, "G") ||strstr(parse_string, "B") ) )
{
fprintf(stderr, "Setting point format to 3, overriding version to 1.2 -- RGB + time\n");
LASHeader_SetDataFormatId(header, 3);
LASHeader_SetVersionMinor(header, 2);
}
else if ((strstr(parse_string, "R") || strstr(parse_string, "G") ||strstr(parse_string, "B") ) )
{
fprintf(stderr, "Setting point format to 2, overriding version to 1.2 -- RGB\n");
LASHeader_SetDataFormatId(header, 2);
LASHeader_SetVersionMinor(header, 2);
}
else if (strstr(parse_string,"t"))
{
fprintf(stderr, "Setting point format to 1\n");
LASHeader_SetDataFormatId(header, 1);
}
else
{
LASHeader_SetDataFormatId(header, 0);
}
LASHeader_SetPointRecordsCount(header, number_of_point_records);
LASHeader_SetScale(header, xyz_scale[0], xyz_scale[1], xyz_scale[2]);
LASHeader_SetOffset(header, xyz_offset[0], xyz_offset[1], xyz_offset[2]);
LASHeader_SetMin(header, xyz_min_dequant[0], xyz_min_dequant[1], xyz_min_dequant[2]);
LASHeader_SetMax(header, xyz_max_dequant[0], xyz_max_dequant[1], xyz_max_dequant[2]);
LASHeader_SetPointRecordsByReturnCount(header, 0, number_of_points_by_return[1]);
LASHeader_SetPointRecordsByReturnCount(header, 1, number_of_points_by_return[2]);
LASHeader_SetPointRecordsByReturnCount(header, 2, number_of_points_by_return[3]);
LASHeader_SetPointRecordsByReturnCount(header, 3, number_of_points_by_return[4]);
LASHeader_SetPointRecordsByReturnCount(header, 4, number_of_points_by_return[5]);
/* reopen input file for the second pass */
file_in = fopen(file_name_in, "r");
if (file_in == 0)
{
fprintf(stderr, "ERROR: could not open '%s' for second pass\n",file_name_in);
exit(1);
}
/*
because the output goes to a file we can do everything in a
single pass and compute the header information along the way
*/
/* open output file */
printf("Creating file...\n");
writer = LASWriter_Create(file_name_out, header, LAS_MODE_WRITE);
if (!writer) {
LASError_Print("Could not open file for write mode ");
exit(1);
}
if (verbose) {
fprintf(stderr,
"scanning %s with parse '%s' writing to %s\n",
file_name_in ,
parse_string,
file_name_out
);
}
/* read the first line */
while (fgets(line, sizeof(char) * MAX_CHARACTERS_PER_LINE, file_in))
{
point = LASPoint_Create();
if (parse(parse_string, line, xyz, point, &gps_time))
{
/* init the bounding box */
VecCopy3dv(xyz_min, xyz);
VecCopy3dv(xyz_max, xyz);
/* mark that we found the first point */
number_of_point_records = 1;
/* create return histogram */
number_of_points_by_return[LASPoint_GetReturnNumber(point)]++;
/* compute the quantized x, y, and z values */
LASPoint_SetX(point, xyz[0]);
LASPoint_SetY(point, xyz[1]);
LASPoint_SetZ(point, xyz[2]);
/* write the first point */
err = LASWriter_WritePoint(writer, point);
if (err) {
LASError_Print("could not write point");
exit(1);
}
/* we can stop this loop */
break;
}
else
{
fprintf(stderr, "WARNING: cannot parse '%s' with '%s'. skipping ...\n",
line,
parse_string);
}
LASPoint_Destroy(point);
point = NULL;
}
/* did we manage to parse a line? */
if (number_of_point_records != 1)
{
fprintf(stderr, "ERROR: could not parse any lines with '%s'\n",
parse_less);
exit(1);
}
/* loop over the remaining lines */
while (fgets(line, sizeof(char) * MAX_CHARACTERS_PER_LINE, file_in))
{
point = LASPoint_Create();
if (parse(parse_string, line, xyz, point, &gps_time))
{
/* update bounding box */
VecUpdateMinMax3dv(xyz_min, xyz_max, xyz);
/* count points */
number_of_point_records++;
/* create return histogram */
number_of_points_by_return[LASPoint_GetReturnNumber(point)]++;
/* compute the quantized x, y, and z values */
LASPoint_SetX(point, xyz[0]);
LASPoint_SetY(point, xyz[1]);
LASPoint_SetZ(point, xyz[2]);
/* write the first point */
err = LASWriter_WritePoint(writer, point);
if (err) {
LASError_Print("could not write point");
exit(1);
}
}
else
{
fprintf(stderr, "WARNING: cannot parse '%s' with '%s'. skipping ...\n",
line,
parse_string);
}
LASPoint_Destroy(point);
point = NULL;
}
/* close up stuff */
fclose(file_in);
LASWriter_Destroy(writer);
if (verbose)
{
fprintf(stderr, "done.\n");
}
return 0;
}
int main(int argc, char *argv[])
{
int i = 0;
char* file_name = NULL;
LASReaderH reader = NULL;
LASHeaderH header = NULL;
LASWriterH writer = NULL;
int check_points = FALSE;
int repair_header = FALSE;
int change_header = FALSE;
int repair_bounding_box = FALSE;
int use_stdin = FALSE;
int update_return_counts = FALSE;
int skip_vlr = FALSE;
int wkt = FALSE;
char *system_identifier = NULL;
char *generating_software = NULL;
unsigned char file_creation_day = 0;
unsigned char file_creation_year = 0;
int err = 0;
LASPointSummary* summary = NULL;
for (i = 1; i < argc; i++)
{
if (strcmp(argv[i],"-v") == 0 || strcmp(argv[i],"--version") == 0)
{
char* ver = LAS_GetFullVersion();
fprintf(stderr,"%s", ver);
LASString_Free(ver);
exit(0);
}
else if (strcmp(argv[i],"-h") == 0 || strcmp(argv[i],"--help") == 0)
{
usage();
exit(0);
}
else if (strcmp(argv[i],"--input") == 0
|| strcmp(argv[i],"-input") == 0
|| strcmp(argv[i],"-i") == 0
|| strcmp(argv[i],"-in") == 0)
{
i++;
file_name = argv[i];
}
else if (strcmp(argv[i], "--points") == 0
|| strcmp(argv[i], "--check") == 0
|| strcmp(argv[i], "--check_points") == 0
|| strcmp(argv[i], "-c") == 0
|| strcmp(argv[i], "-points") == 0
|| strcmp(argv[i], "-check") == 0
|| strcmp(argv[i], "-check_points") == 0)
{
check_points = TRUE;
}
else if (strcmp(argv[i], "--nocheck") == 0
|| strcmp(argv[i], "-nocheck") == 0)
{
check_points = FALSE;
}
else if (strcmp(argv[i], "--stdin") == 0
|| strcmp(argv[i], "-ilas") == 0)
{
use_stdin = TRUE;
}
else if (strcmp(argv[i], "--repair") == 0
|| strcmp(argv[i], "-r") == 0
|| strcmp(argv[i], "-repair_header") == 0
|| strcmp(argv[i], "-repair") == 0)
{
repair_header = TRUE;
check_points = TRUE;
}
else if (strcmp(argv[i], "--repair_bb") == 0
|| strcmp(argv[i], "--repair_bounding_box") == 0
|| strcmp(argv[i], "--repair_boundingbox") == 0
|| strcmp(argv[i], "-repair_bb") == 0
|| strcmp(argv[i], "-repair_bounding_box") == 0
|| strcmp(argv[i], "-repair_boundingbox") == 0
|| strcmp(argv[i], "-repair") == 0
|| strcmp(argv[i], "-rb") == 0)
{
repair_bounding_box = TRUE;
check_points = TRUE;
}
else if (strcmp(argv[i],"--system_identifier") == 0
|| strcmp(argv[i],"-system_identifier") == 0
|| strcmp(argv[i],"-s") == 0
|| strcmp(argv[i],"-sys_id") == 0)
{
i++;
system_identifier = (char*) malloc(31 * sizeof(char));
strcpy(system_identifier, argv[i]);
change_header = TRUE;
}
else if (strcmp(argv[i],"--generating_software") == 0
|| strcmp(argv[i],"-generating_software") == 0
|| strcmp(argv[i],"-g") == 0
|| strcmp(argv[i],"-gen_soft") == 0)
{
i++;
generating_software = (char*) malloc(31*sizeof(char));
strcpy(generating_software, argv[i]);
change_header = TRUE;
}
else if (strcmp(argv[i],"--file_creation") == 0
|| strcmp(argv[i],"-file_creation") == 0)
{
/* XXX - mloskot: Consider replacing atoi with strtol,
see http://www.iso-9899.info/wiki/Converting */
i++;
file_creation_day = (unsigned char)atoi(argv[i]);
i++;
file_creation_year = (unsigned char)atoi(argv[i]);
change_header = TRUE;
}
else if (strcmp(argv[i],"--skip_vlr") == 0 || strcmp(argv[i],"--no_vlr") == 0)
{
skip_vlr = TRUE;
}
else if (strcmp(argv[i],"--wkt") == 0)
{
wkt = TRUE;
}
else if (file_name == NULL)
{
file_name = argv[i];
}
else
{
usage();
fprintf(stderr, "ERROR: unknown argument '%s'\n",argv[i]);
exit(1);
}
}
if (use_stdin) {
file_name = "stdin";
}
if (!file_name) {
LASError_Print("No filename was provided to be opened");
usage();
exit(1);
}
reader = LASReader_Create(file_name);
if (!reader) {
LASError_Print("Could not open file ");
exit(1);
}
header = LASReader_GetHeader(reader);
if (!header) {
LASError_Print("Could not get LASHeader ");
exit(1);
}
print_header(stdout, header, file_name, skip_vlr, wkt);
if (change_header) {
if (system_identifier) {
err = LASHeader_SetSystemId (header, system_identifier);
if (err) LASError_Print("Could not set SystemId");
}
if (generating_software) {
err = LASHeader_SetSoftwareId(header, generating_software);
if (err) LASError_Print("Could not set SoftwareId");
}
if ( file_creation_day || file_creation_year) {
err = LASHeader_SetCreationDOY(header, file_creation_day);
if (err) LASError_Print("Could not set file creation day");
err = LASHeader_SetCreationYear(header, file_creation_year);
if (err) LASError_Print("Could not set file creation year");
}
/* We need to wipe out the reader and make a writer. */
if (reader) {
LASReader_Destroy(reader);
reader = NULL;
}
writer = LASWriter_Create(file_name, header, LAS_MODE_APPEND);
if (!writer) {
LASError_Print("Problem creating LASWriterH object");
LASHeader_Destroy(header);
header = NULL;
exit(1);
}
if (writer) LASWriter_Destroy(writer);
writer = NULL;
if (header) LASHeader_Destroy(header);
header = NULL;
}
if (check_points)
{
if (!reader) {
reader = LASReader_Create(file_name);
if (!reader) {
LASError_Print("Could not open file ");
exit(1);
}
}
if (! header) {
header = LASReader_GetHeader(reader);
if (!header) {
LASError_Print("Could not get LASHeader ");
exit(1);
}
}
if (!summary)
summary = SummarizePoints(reader);
print_point_summary(stdout, summary, header);
if (repair_header) {
fprintf(stdout, "\n---------------------------------------------------------\n");
fprintf(stdout, " Repair Summary\n");
fprintf(stdout, "---------------------------------------------------------\n");
if (use_stdin) {
LASError_Print("Cannot update header information on piped input!");
exit(1);
}
if (! header) {
header = LASReader_GetHeader(reader);
if (!header) {
LASError_Print("Could not get LASHeader ");
exit(1);
}
}
if (! repair_bounding_box) {
if ( LASHeader_GetMinX(header) != LASPoint_GetX(summary->pmin) )
repair_bounding_box = TRUE;
if ( LASHeader_GetMinY(header) != LASPoint_GetY(summary->pmin) )
repair_bounding_box = TRUE;
if ( LASHeader_GetMinZ(header) != LASPoint_GetZ(summary->pmin) )
repair_bounding_box = TRUE;
if ( LASHeader_GetMaxX(header) != LASPoint_GetX(summary->pmax) )
repair_bounding_box = TRUE;
if ( LASHeader_GetMaxY(header) != LASPoint_GetY(summary->pmax) )
repair_bounding_box = TRUE;
if ( LASHeader_GetMaxZ(header) != LASPoint_GetZ(summary->pmax) )
repair_bounding_box = TRUE;
}
if (repair_bounding_box) {
fprintf(stdout, " Reparing Bounding Box...\n");
err = LASHeader_SetMin( header,
LASPoint_GetX(summary->pmin),
LASPoint_GetY(summary->pmin),
LASPoint_GetZ(summary->pmin)
);
if (err) {
LASError_Print("Could not set minimum for header ");
exit(1);
}
err = LASHeader_SetMax( header,
LASPoint_GetX(summary->pmax),
LASPoint_GetY(summary->pmax),
LASPoint_GetZ(summary->pmax)
);
if (err) {
LASError_Print("Could not set minimum for header ");
exit(1);
}
}
for (i = 0; i < 5; i++) {
if (LASHeader_GetPointRecordsByReturnCount(header, i) !=
summary->number_of_points_by_return[i])
{
update_return_counts = TRUE;
break;
}
}
if (update_return_counts) {
fprintf(stdout, " Reparing Point Count by Return...\n");
for (i = 0; i < 5; i++) {
LASHeader_SetPointRecordsByReturnCount( header,
i,
summary->number_of_points_by_return[i]);
}
}
if (reader) {
LASReader_Destroy(reader);
reader = NULL;
}
writer = LASWriter_Create(file_name, header, LAS_MODE_APPEND);
if (!writer) {
LASError_Print("Problem creating LASWriterH object for append");
LASHeader_Destroy(header);
header = NULL;
exit(1);
}
LASWriter_Destroy(writer);
writer = NULL;
LASHeader_Destroy(header);
header = NULL;
}
if (summary) {
LASPoint_Destroy(summary->pmin);
LASPoint_Destroy(summary->pmax);
free(summary);
}
}
if (reader) LASReader_Destroy(reader);
if (header) LASHeader_Destroy(header);
#ifdef HAVE_GDAL
/* Various GDAL related cleanups */
#ifdef OSRCleanup
OSRCleanup();
#endif
CPLFinderClean();
CPLFreeConfig();
CPLCleanupTLS();
#endif
return 0;
}