That’s great. I’ll be very interested to hear how the search works. You’ve done the difficult part. Here is some sample code:
ff = dir(pattern = "nrrd$")
# it's a good id to set this to some kind of ifentifier for each line
names(ff) <- tools::file_path_sans_ext(basename(ff))
# choice of k may depend somewhat on spatial scale of branches etc in
# your expression patterns
gal4s=nlapply(ff[1:3], dotprops, k=5)
# you can add additional metadata to the neuronlist object like so
gal4s[,'id']=names(gal4s)
gal4s[,]=merge(gal4s[,], metadata_table)
Greg.