int main()
{
struct ttyrec_header th;
char buf[BUFFER_SIZE];
int s,n,r;
while (1) {
if (read(0, &th, 12)!=12)
return 0;
n=from_little_endian(th.len);
while (n>0) {
s=(n>BUFFER_SIZE)?BUFFER_SIZE:n;
if ((r=read(0, buf, s))<=0) {
fprintf(stderr, "%s\n", r?strerror(errno):"File was truncated");
return 1;
}
if (write(1, buf, r)!=r) {
fprintf(stderr, "Write error\n");
return 1;
}
n-=r;
}
}
return 0;
}
size_t Store::readSize(FILE *file)
{
uint32_t size = 0;
size_t items_read = fread(&size, sizeof(size), 1, file);
from_little_endian(uint32_t, bytes_read);
if(items_read < 1)
{
fseek(file, -1 * items_read, SEEK_CUR);
throw std::runtime_error("Not enough data");
}
return size;
}
bool LASreaderSHP::read_point_default()
{
if (point_count == number_of_points)
{
int i, int_input;
if (fread(&int_input, sizeof(int), 1, file) != 1) { npoints = p_count; return FALSE; }; // record number (BIG)
if (fread(&int_input, sizeof(int), 1, file) != 1) { npoints = p_count; return FALSE; }; // content length (BIG)
if (fread(&int_input, sizeof(int), 1, file) != 1) { npoints = p_count; return FALSE; }; // shape type (LITTLE)
from_little_endian(&int_input);
if (int_input != shape_type)
{
fprintf(stderr, "WARNING: wrong shape type %d != %d in record\n", int_input, shape_type);
}
double double_input;
if (shape_type == 8 || shape_type == 18 || shape_type == 28) // Multipoint
{
if (fread(&double_input, sizeof(double), 1, file) != 1) { npoints = p_count; return FALSE; }; // xmin (LITTLE)
if (fread(&double_input, sizeof(double), 1, file) != 1) { npoints = p_count; return FALSE; }; // ymin (LITTLE)
if (fread(&double_input, sizeof(double), 1, file) != 1) { npoints = p_count; return FALSE; }; // xmax (LITTLE)
if (fread(&double_input, sizeof(double), 1, file) != 1) { npoints = p_count; return FALSE; }; // ymax (LITTLE)
if (fread(&int_input, sizeof(int), 1, file) != 1) { npoints = p_count; return FALSE; }; // number of points (LITTLE)
from_little_endian(&int_input);
number_of_points = int_input;
}
else // or regular point?
{
number_of_points = 1;
}
if (shape_type == 11 || shape_type == 18)
{
if (points_allocated < number_of_points)
{
if (points) delete [] points;
points = new I32[number_of_points*3];
points_allocated = number_of_points;
}
// read points x and y
for (i = 0; i < number_of_points; i++)
{
if (fread(&double_input, sizeof(double), 1, file) != 1) { npoints = p_count; return FALSE; }; // x of point (LITTLE)
from_little_endian(&double_input);
points[3*i+0] = header.get_X(double_input);
if (fread(&double_input, sizeof(double), 1, file) != 1) { npoints = p_count; return FALSE; }; // y of point (LITTLE)
from_little_endian(&double_input);
points[3*i+1] = header.get_Y(double_input);
}
// read points z and write LAS points
if (shape_type == 18) // Multipoint
{
if (fread(&double_input, sizeof(double), 1, file) != 1) { npoints = p_count; return FALSE; }; // zmin (LITTLE)
if (fread(&double_input, sizeof(double), 1, file) != 1) { npoints = p_count; return FALSE; }; // zmin (LITTLE)
}
for (i = 0; i < number_of_points; i++)
{
if (fread(&double_input, sizeof(double), 1, file) != 1) { npoints = p_count; return FALSE; }; // z of point (LITTLE)
from_little_endian(&double_input);
points[3*i+2] = header.get_Z(double_input);
}
}
else
{
if (points_allocated < number_of_points)
{
if (points) delete [] points;
points = new I32[number_of_points*2];
points_allocated = number_of_points;
}
// read points x and y and write LAS points
for (i = 0; i < number_of_points; i++)
{
if (fread(&double_input, sizeof(double), 1, file) != 1) { npoints = p_count; return FALSE; }; // x of point (LITTLE)
from_little_endian(&double_input);
points[2*i+0] = header.get_Z(double_input);
if (fread(&double_input, sizeof(double), 1, file) != 1) { npoints = p_count; return FALSE; }; // y of point (LITTLE)
from_little_endian(&double_input);
points[2*i+1] = header.get_Z(double_input);
}
}
// read points m
if (shape_type == 11 || shape_type == 21 || shape_type == 18 || shape_type == 28)
{
if (shape_type == 18 || shape_type == 28) // Multipoint
{
if (fread(&double_input, sizeof(double), 1, file) != 1) { npoints = p_count; return FALSE; }; // mmin (LITTLE)
if (fread(&double_input, sizeof(double), 1, file) != 1) { npoints = p_count; return FALSE; }; // mmin (LITTLE)
}
for (i = 0; i < number_of_points; i++)
{
if (fread(&double_input, sizeof(double), 1, file) != 1) { npoints = p_count; return FALSE; }; // m of point (LITTLE)
}
}
point_count = 0;
}
if (shape_type == 11 || shape_type == 18)
{
point.set_X(points[3*point_count+0]);
point.set_Y(points[3*point_count+1]);
point.set_Z(points[3*point_count+2]);
}
else
{
point.set_X(points[2*point_count+0]);
point.set_Y(points[2*point_count+1]);
point.set_Z(0);
}
p_count++;
point_count++;
return TRUE;
}
bool LASreaderSHP::open(const char* file_name)
{
if (file_name == 0)
{
fprintf(stderr,"ERROR: fine name pointer is zero\n");
return FALSE;
}
clean();
file = fopen_compressed(file_name, "rb", &piped);
if (file == 0)
{
fprintf(stderr, "ERROR: cannot open file '%s'\n", file_name);
return FALSE;
}
// clean the header
header.clean();
// populate the header as much as it makes sense
sprintf(header.system_identifier, "LAStools (c) by rapidlasso GmbH");
sprintf(header.generating_software, "via LASreaderSHP (%d)", LAS_TOOLS_VERSION);
// maybe set creation date
#ifdef _WIN32
WIN32_FILE_ATTRIBUTE_DATA attr;
SYSTEMTIME creation;
GetFileAttributesEx(file_name, GetFileExInfoStandard, &attr);
FileTimeToSystemTime(&attr.ftCreationTime, &creation);
int startday[13] = {-1, 0, 31, 59, 90, 120, 151, 181, 212, 243, 273, 304, 334};
header.file_creation_day = startday[creation.wMonth] + creation.wDay;
header.file_creation_year = creation.wYear;
// leap year handling
if ((((creation.wYear)%4) == 0) && (creation.wMonth > 2)) header.file_creation_day++;
#else
header.file_creation_day = 111;
header.file_creation_year = 2011;
#endif
header.point_data_format = 0;
header.point_data_record_length = 20;
// initialize point
point.init(&header, header.point_data_format, header.point_data_record_length);
// read SHP header and populate the LAS header with the bounding box
int int_input;
if (fread(&int_input, sizeof(int), 1, file) != 1) return false; // file code (BIG)
from_big_endian(&int_input);
if (int_input != 9994)
{
fprintf(stderr, "ERROR: wrong shapefile code %d != 9994\n", int_input);
return FALSE;
}
if (fread(&int_input, sizeof(int), 1, file) != 1) return false; // unused (BIG)
if (fread(&int_input, sizeof(int), 1, file) != 1) return false; // unused (BIG)
if (fread(&int_input, sizeof(int), 1, file) != 1) return false; // unused (BIG)
if (fread(&int_input, sizeof(int), 1, file) != 1) return false; // unused (BIG)
if (fread(&int_input, sizeof(int), 1, file) != 1) return false; // unused (BIG)
if (fread(&int_input, sizeof(int), 1, file) != 1) return false; // file length (BIG)
from_big_endian(&int_input);
int file_length = int_input;
if (fread(&int_input, sizeof(int), 1, file) != 1) return false; // version (LITTLE)
from_little_endian(&int_input);
if (int_input != 1000)
{
fprintf(stderr, "ERROR: wrong shapefile version %d != 1000\n", int_input);
return FALSE;
}
if (fread(&int_input, sizeof(int), 1, file) != 1) return false; // shape type (LITTLE)
from_little_endian(&int_input);
shape_type = int_input;
if (shape_type != 1 && shape_type != 11 && shape_type != 21 && shape_type != 8 && shape_type != 18 && shape_type != 28)
{
fprintf(stderr, "ERROR: wrong shape type %d != 1,11,21,8,18,28\n", shape_type);
return FALSE;
}
double double_input;
if (fread(&double_input, sizeof(double), 1, file) != 1) return false; // xmin (LITTLE)
from_little_endian(&double_input);
header.min_x = double_input;
if (fread(&double_input, sizeof(double), 1, file) != 1) return false; // ymin (LITTLE)
from_little_endian(&double_input);
header.min_y = double_input;
if (fread(&double_input, sizeof(double), 1, file) != 1) return false; // xmax (LITTLE)
from_little_endian(&double_input);
header.max_x = double_input;
if (fread(&double_input, sizeof(double), 1, file) != 1) return false; // ymax (LITTLE)
from_little_endian(&double_input);
header.max_y = double_input;
if (fread(&double_input, sizeof(double), 1, file) != 1) return false; // zmin (LITTLE)
from_little_endian(&double_input);
header.min_z = double_input;
if (fread(&double_input, sizeof(double), 1, file) != 1) return false; // zmax (LITTLE)
from_little_endian(&double_input);
header.max_z = double_input;
if (fread(&double_input, sizeof(double), 1, file) != 1) return false; // mmin (LITTLE)
from_little_endian(&double_input);
if (fread(&double_input, sizeof(double), 1, file) != 1) return false; // mmax (LITTLE)
from_little_endian(&double_input);
// maybe calculate npoints
if (shape_type == 1)
{
npoints = (file_length-50)/(14);
}
else if (shape_type == 11)
{
npoints = (file_length-50)/(22);
}
else if (shape_type == 21)
{
npoints = (file_length-50)/(18);
}
else if (shape_type == 8)
{
npoints = (file_length-50-20)/(8); // over-estimate (assumes all in one record)
}
else if (shape_type == 18)
{
npoints = (file_length-50-36)/(16); // over-estimate (assumes all in one record)
}
else if (shape_type == 28)
{
npoints = (file_length-50-28)/(12); // over-estimate (assumes all in one record)
}
header.number_of_point_records = (U32)npoints;
header.number_of_points_by_return[0] = (U32)npoints;
// figure out the right scale factor and offset
populate_scale_and_offset();
// populate the quantized bounding box
populate_bounding_box();
p_count = 0;
return TRUE;
}
