♦️ Add support for reading and saving .xy and .xyz files in the FastPointQuery module.

Author: ZenanH

Project.toml

@@ -5,13 +5,15 @@ version = "0.1.4"
 
 [deps]
 CondaPkg = "992eb4ea-22a4-4c89-a5bb-47a3300528ab"
+DelimitedFiles = "8bb1440f-4735-579b-a4ab-409b98df4dab"
 Downloads = "f43a241f-c20a-4ad4-852c-f6b1247861c6"
 Logging = "56ddb016-857b-54e1-b83d-db4d58db5568"
 PrecompileTools = "aea7be01-6a6a-4083-8856-8a6e6704d82a"
 PythonCall = "6099a3de-0909-46bc-b1f4-468b9a2dfc0d"
 
 [compat]
 CondaPkg = "0.2"
+DelimitedFiles = "1.9"
 PythonCall = "0.9"
 PrecompileTools = "1.2"
 julia = "1.11"

src/FastPointQuery.jl

@@ -1,6 +1,6 @@
 module FastPointQuery
 
-using CondaPkg, Downloads, Logging, PrecompileTools, PythonCall
+using CondaPkg, DelimitedFiles, Downloads, Logging, PrecompileTools, PythonCall
 
 # Python packages
 const np         = PythonCall.pynew()
@@ -39,14 +39,15 @@ function __init__()
     end
 end
 
-include(joinpath(@__DIR__, "fileio/asc.jl"    ))
+include(joinpath(@__DIR__, "fileio/asc.jl"))
 include(joinpath(@__DIR__, "fileio/geojson.jl"))
-include(joinpath(@__DIR__, "fileio/ply.jl"    ))
-include(joinpath(@__DIR__, "fileio/stl.jl"    ))
-include(joinpath(@__DIR__, "fileio/tiff.jl"   ))
+include(joinpath(@__DIR__, "fileio/ply.jl"))
+include(joinpath(@__DIR__, "fileio/stl.jl"))
+include(joinpath(@__DIR__, "fileio/tiff.jl"))
+include(joinpath(@__DIR__, "fileio/xyz.jl"))
 
-include(joinpath(@__DIR__, "utils.jl"     ))
-include(joinpath(@__DIR__, "polygon.jl"   ))
+include(joinpath(@__DIR__, "utils.jl"))
+include(joinpath(@__DIR__, "polygon.jl"))
 include(joinpath(@__DIR__, "polyhedron.jl"))
 
 export get_resource, res_dir

src/fileio/xyz.jl

@@ -0,0 +1,67 @@
+#==========================================================================================+
+| TABLE OF CONTENTS:                                                                       |
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+|  - readxy                                                                                |
+|  - readxyz                                                                               |
+|  - savexy                                                                                |
+|  - savexyz                                                                               |
++==========================================================================================#
+
+export readxy
+export readxyz
+export savexy
+export savexyz
+
+"""
+    readxy(file_dir::P) where P <: String
+
+Description:
+---
+Read the 2D `.xy` file from `file_dir`.
+"""
+function readxy(file_dir::P) where P <: String
+    xy = readdlm(file_dir, ' ')[:, 1:2]
+    return Array{Float64}(xy)
+end
+
+"""
+    readxyz(file_dir::P) where P <: String
+
+Description:
+---
+Read the 3D `.xyz` file from `file_dir`.
+"""
+function readxyz(file_dir::P) where P <: String
+    xyz = readdlm(file_dir, ' ')[:, 1:3]
+    return Array{Float64}(xyz)
+end
+
+"""
+    savexy(file_dir::P, pts::T) where {P <: String, T <: AbstractMatrix}
+
+Description:
+---
+Save the 2D points `pts` to the `.xy` file (`file_dir`).
+"""
+function savexy(file_dir::P, pts::T) where {P <: String, T <: AbstractMatrix}
+    size(pts, 2) == 2 || throw(ArgumentError("The input points should have 2 columns."))
+    open(file_dir, "w") do io
+        writedlm(io, pts, ' ')
+    end
+end
+savexy(pts::T, file_dir::P) where {P <: String, T <: AbstractMatrix} = savexy(file_dir, pts)
+
+"""
+    savexyz(file_dir::P, pts::T) where {P <: String, T <: AbstractMatrix}
+
+Description:
+---
+Save the 3D points `pts` to the `.xyz` file (`file_dir`).
+"""
+function savexyz(file_dir::P, pts::T) where {P <: String, T <: AbstractMatrix}
+    size(pts, 2) == 3 || throw(ArgumentError("The input points should have 3 columns."))
+    open(file_dir, "w") do io
+        writedlm(io, pts, ' ')
+    end
+end
+savexyz(pts::T, file_dir::P) where {P <: String, T <: AbstractMatrix} = savexyz(file_dir, pts)