las-py is a zero-dependency Python library for parsing .Las file (Geophysical/Canadian well log files).
Currently supports only version 2.0 of LAS Specification. For more information about this format, see the Canadian Well Logging Society web page
-
What's new in 1.1.0
- Export to csv
- Export to csv without rows containing null values
- Bug fixes
-
To Install
$pip insatll las-py
-
Usage
from las_py import Laspy
my_las = Laspy('path_to_las_file.las')
-
Read data
data = my_las.data print(data) #[[2650.0, 177.825, -999.25, -999.25], [2650.5, 182.5, -999.25,-999.25], [2651.0,180.162, -999.25, -999.25], [2651.5, 177.825, -999.25, -999.25], [2652.0, 177.825, -999.25, -999.25] ...]
# get data with rows that has null value stripped data = my_las.data_stripped print(data) #[[2657.5, 212.002, 0.16665, 1951.74597], [2658.0, 201.44, 0.1966, 1788.50696], [2658.5, 204.314, 0.21004, 1723.21204], [2659.0, 212.075, 0.22888, 1638.328], [2659.5, 243.536, 0.22439, 1657.91699]...]
-
Get the log headers
headers = my_las.header print(headers) # ['DEPTH', 'GR', 'NPHI', 'RHOB']
-
Get the log headers descriptions
hds_and_desc = my_las.header_and_descr print(hds_and_desc) # {DEPTH': 'DEPTH', 'GR': 'Gamma Ray', 'NPHI': 'Neutron Porosity','RHOB': 'Bulk density'}
-
Get a particular column, say Gamma Ray log
GR = my_las.column('GR') print(GR) # [-999.25, -999.25, -999.25, -999.25, -999.25, 122.03, 123.14, ...]
# get column with null values stripped GR = my_las.column_stripped('GR') print(GR) # [61.61, 59.99, 54.02, 50.87, 54.68, 64.39, 77.96, ...]
Note this returns the column, after all the data has been stripped off their null values, which means that valid data in a particular column would be stripped off if there is another column that has a null value at that particular row
-
Get the Well Parameters
1. descr - Description/ Full name of the well parameter 2. units - Its unit measurements 3. value - Value
start = my_las.well.STRT.value # 1670.0 stop = my_las.well.STOP.value # 1669.75 null_value = my_las.well.NULL.value # -999.25 # Any other well parameter present in the file, canbe gotten with the same syntax above
-
Get the Curve Parameters
1. descr - Description/ Full name of the log column 2. units - Unit of the log column measurements 3. value - API value of the log column
NPHI = my_las.curve.NPHI.descr # 'Neutron Porosity' RHOB = my_las.curve.RHOB.descr # 'Bulk density' # This is the same for all log column present in the file
-
Get the Parameters of the well
1. descr - Description/ Full name of the log column 2. units - Unit of the log column measurements 3. value - API value of the log column
BHT = my_las.param.BHT.descr # 'BOTTOM HOLE TEMPERATURE' BHT_valaue = my_las.param.BHT.value # 35.5 BHT_units = my_las.param.BHT.units # 'DEGC' # This is the same for all well parameters present in the file
-
Get the number of rows and columns
rows = my_las.row_count # 4 columns = my_las.column_count # 3081
-
Get the version and wrap
version = my_las.version # '2.0' wrap = my_las.wrap # 'YES'
-
Get other information
other = my_las.other print(other) # Note: The logging tools became stuck at 625 metres causing the data # between 625 metres and 615 metres to be invalid.
-
Export to CSV
This writes a csv file to the current working directory, with headers of the well and data section only.
my_las.to_csv('result') # result.csv has been created Successfully!
result.csv
DEPT RHOB GR NPHI 0.5 -999.25 -999.25 -0.08 1.0 -999.25 -999.25 -0.08 1.5 -999.25 -999.25 -0.04 ... ... ... ... 1.3 -999.25 -999.25 -0.08 Or get the version of csv with null values stripped
my_las.to_csv_stripped('clean') # clean.csv has been created Successfully!
clean.csv
DEPT RHOB GR NPHI 80.5 2.771 18.6 -6.08 81.0 2.761 17.4 -6.0 81.5 2.752 16.4 -5.96 ... ... ... ... 80.5 2.762 16.2 -5.06 -
las-py is an MIT-licensed open source project. You can help it grow by becoming a sponsor/supporter. Donate on Patreon