detector response bug fix 2, mnlocs (#40)

Author: dehann

Project.toml

@@ -12,6 +12,7 @@ LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 LoopVectorization = "bdcacae8-1622-11e9-2a5c-532679323890"
 RecursiveArrayTools = "731186ca-8d62-57ce-b412-fbd966d074cd"
 Rotations = "6038ab10-8711-5258-84ad-4b1120ba62dc"
+SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 StatsBase = "2913bbd2-ae8a-5f71-8c99-4fb6c76f3a91"
 
@@ -23,6 +24,7 @@ LoopVectorization = "0.12.173"
 LieGroups = "0.1"
 RecursiveArrayTools = "3.27.0"
 Rotations = "1"
+SparseArrays = "1"
 StaticArrays = "1"
 StatsBase = "0.34"
 julia = "1.10"

src/CameraModels.jl

@@ -3,6 +3,7 @@ module CameraModels
 using LinearAlgebra
 using LieGroups
 using DocStringExtensions
+using SparseArrays
 using StaticArrays
 using FixedPointNumbers
 using StatsBase

src/services/RadianceCorrection.jl

@@ -166,7 +166,7 @@ function solveDetectorResponse(
   #                 logExposure]    zeros(n)]
   #
   neqs = sum(mnlocs.*mnimgs) + 1 + n
-  A = zeros(neqs,n+nlocs)
+  A = spzeros(neqs,n+nlocs)
   b = zeros(neqs)
 
   # Fill in the equations from pixel value observations in A and b,
@@ -184,7 +184,7 @@ function solveDetectorResponse(
         # A2 columns correspond to logExposure of observed pixels from list of locations i
         pixz1 = Z[j][i] + 1
         miloc = offset+n+i
-        a_12 = view(A, k, [pixz1, miloc])
+        a_12 = view(A, k, SA[pixz1, miloc])
         a_12 .= window[pixz1] .* SA[1, -1]
         b[k] = window[pixz1] * logΔTs[j]
         k += 1
@@ -211,7 +211,7 @@ function solveDetectorResponse(
   x = A\b
 
   # Extract gcurve and log exposure values from solution vector
-  gcurve = x[1:n]
+  gcurve = collect(x[1:n])
   logExposure = x[(n+1):end]
   return gcurve, logExposure
 end