define __copy__ for headers, and use it in create_from_header

docs/examples.rst

@@ -15,7 +15,7 @@ We use the classification field to filter points, but this can work with the oth
 
     las = pylas.read('pylastests/simple.las')
 
-    new_file = pylas.create()
+    new_file = pylas.create(point_format_id=las.header.point_format_id, file_version=las.header.version)
     new_file.points = las.points[las.classification == 1]
 
     new_file.write('extracted_points.las')

docs/intro.rst

@@ -101,7 +101,7 @@ the same dimensions plus some additional dimensions:
 +----------------------+-----------+--------------+
 | Added dimensions     |   Type    |  Size (bit)  |
 +======================+===========+==============+
-| gps_time             |  Floating |       64     |
+| gps_time             |  Floating |      64      |
 +----------------------+-----------+--------------+
 
 
@@ -122,7 +122,7 @@ the same dimensions plus some additional dimensions:
 +----------------------+-----------+--------------+
 | Added dimensions     |   Type    |  Size (bit)  |
 +======================+===========+==============+
-| gps_time             |  Floating |       64     |
+| gps_time             |  Floating |      64      |
 +----------------------+-----------+--------------+
 | red                  |  unsigned |      16      |
 +----------------------+-----------+--------------+

pylas/headers/rawheader.py

@@ -153,10 +153,21 @@ def point_format_id(self):
     def point_format_id(self, value):
         self._point_data_format_id = value
 
+    @property
+    def point_size(self):
+        return self.point_data_record_length
+
+    @point_size.setter
+    def point_size(self, value):
+        self.point_data_record_length = value
+
     @property
     def are_points_compressed(self):
         return compression.is_point_format_compressed(self._point_data_format_id)
 
+    def __copy__(self):
+        return self.__class__.from_buffer_copy(bytes(self))
+
     def __repr__(self):
         return "<LasHeader({})>".format(self.version)
 

pylas/lib.py

@@ -1,6 +1,7 @@
 """ 'Entry point' of the library, Contains the various functions meant to be
 used directly by a user
 """
+import copy
 
 from . import headers
 from .lasdatas import las12, las14
@@ -81,7 +82,10 @@ def mmap_las(filename):
 
 
 def create_from_header(header):
-    """ Creates a File from an existing header
+    """ Creates a File from an existing header,
+    allocating the array of point according to the provided header.
+    The input header is copied.
+
 
     Parameters
     ----------
@@ -91,14 +95,15 @@ def create_from_header(header):
     -------
     pylas.lasdatas.base.LasBase
     """
+    header = copy.copy(header)
     points = record.PackedPointRecord.zeros(header.point_format_id, header.point_count)
     if header.version >= "1.4":
         return las14.LasData(header=header, points=points)
     return las12.LasData(header=header, points=points)
 
 
 def create_las(*, point_format=0, file_version=None):
-    """ Function to create a new las data object
+    """ Function to create a new empty las data object
 
     .. note::
 

pylastests/test_header.py

@@ -43,3 +43,16 @@ def test_nb_points_return_1_4():
     assert tuple(las.header.number_of_points_by_return) == ((1,) * 14) + (
         len(las.points_data) - 14,
     )
+
+
+def test_header_copy():
+    import copy
+    las = pylas.read(test_common.simple_las)
+    header_copy = copy.copy(las.header)
+
+    assert header_copy.point_format_id == las.header.point_format_id
+    assert header_copy.version == las.header.version
+
+    header_copy.point_format_id = 0
+    assert header_copy.point_format_id != las.header.point_format_id
+    assert header_copy.version == las.header.version

setup.py

@@ -1,19 +1,20 @@
 from setuptools import setup, find_packages
 
-with open('README.rst') as f:
+with open("README.rst") as f:
     readme = f.read()
 
 setup(
-    name='pylas',
-    version='0.1.3',
-    description='Las/Laz reading and writing in python',
+    name="pylas",
+    version="0.1.3",
+    description="Las/Laz reading and writing in python",
     long_description=readme,
-    url='https://github.com/tmontaigu/pylas',
-    author='Thomas Montaigu',
-    author_email='thomas.montaigu@laposte.net',
-    python_requires='>=3',
-    install_requires=['numpy'],
-    license='BSD 3-Clause',
-    packages=find_packages(exclude=('pylastests',)),
-    zip_safe=False
+    url="https://github.com/tmontaigu/pylas",
+    author="Thomas Montaigu",
+    author_email="thomas.montaigu@laposte.net",
+    python_requires=">=3",
+    keywords="las lidar",
+    install_requires=["numpy"],
+    license="BSD 3-Clause",
+    packages=find_packages(exclude=("pylastests",)),
+    zip_safe=False,
 )