sensor-vis/vis.R at master · fluo-2/sensor-vis · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
#!/usr/bin/env Rscript
# -*- coding: utf-8 -*-

library(tools)
library(gdata)
library(ggplot2)
library(tikzDevice)
library(reshape2)
library(ggspectra)
source("./handle_data.R")

# read main data source for sheet 1
clean_data <- extract_categorical_data("WL_FWHM")
# create dummy dataframe to initialize wavelength range
wave <- data.frame(wl_center = 100:15000)
# proceed to plot data
tikz("basic_spectral.tex", width=20, height=12, standAlone = TRUE)
g <- ggplot(wave,aes(x=wl_center)) +
  ## wl_guide(alpha=0.8) +
  geom_rect(data=clean_data,aes(xmin=wl_center-(fwhm/2),xmax=wl_center+(fwhm/2),
                ymin=0, ymax=0.49),color="black",fill="red",alpha=0.5,size=1.1) +
  ## geom_vline(xintercept = 685, linetype="dashed", size=1.2) +
  ylab("") +
  xlab("\n Wavelength $\\lambda$ [nm]") +
  ylim(c(0,0.5)) +
  theme_bw() +
  theme(text = element_text(size=30),
        legend.position = "none",
        plot.title = element_text(hjust=0.5),
        axis.ticks.length = unit(0.2, "cm"),
        axis.ticks.y=element_blank(),
        axis.text.y=element_blank(),
        plot.margin = margin(10, 50, 10, 10)) +
  scale_x_continuous(breaks = round(seq(roundUp(min(wave$wl_center),500),
                                        roundUp(max(wave$wl_center),500),
                                        by = 1000),1),
                     expand=expand_scale(mult=c(0.03,0.01))) +
  facet_wrap(Sensor~.,ncol=1)
print(g)
dev.off()
texi2pdf("basic_spectral.tex",clean=TRUE)
file.remove("basic_spectral.tex")
file.rename("basic_spectral.pdf",
            "./img/basic_spectral.pdf")
# additional truncated plot
wave <- data.frame(wl_center = 100:1200)
# proceed to plot data
tikz("basic_spectral_truncated.tex", width=20, height=12, standAlone = TRUE)
g <- ggplot(wave,aes(x=wl_center)) +
  wl_guide(alpha=0.8) +
  geom_rect(data=clean_data,aes(xmin=wl_center-(fwhm/2),xmax=wl_center+(fwhm/2),
                ymin=0, ymax=0.49),color="black",fill="red",alpha=0.5,size=1.1) +
  ## geom_vline(xintercept = 685, linetype="dashed", size=1.2) +
  ylab("") +
  xlab("\n Wavelength $\\lambda$ [nm]") +
  ## xlim(c(min(wave$wl_center),max(wave$wl_center))) +
  ylim(c(0,0.5)) +
  theme_bw() +
  theme(text = element_text(size=30),
        legend.position = "none",
        plot.title = element_text(hjust=0.5),
        axis.ticks.length = unit(0.2, "cm"),
        axis.ticks.y=element_blank(),
        axis.text.y=element_blank(),
        plot.margin = margin(10, 50, 10, 10)) +
  scale_x_continuous(limits=c(min(wave$wl_center),max(wave$wl_center)),
                     breaks = round(seq(roundUp(min(wave$wl_center),100),
                                        roundUp(max(wave$wl_center),100),
                                        by =100),1),
                     expand=expand_scale(mult=c(0,0))) +
  facet_wrap(Sensor~.,ncol=1)
print(g)
dev.off()
texi2pdf("basic_spectral_truncated.tex",clean=TRUE)
file.remove("basic_spectral_truncated.tex")
file.rename("basic_spectral_truncated.pdf",
            "./img/basic_spectral_truncated.pdf")
# without spectrum
tikz("basic_spectral_truncated_nos.tex", width=20, height=12, standAlone = TRUE)
g <- ggplot(wave,aes(x=wl_center)) +
  ## wl_guide(alpha=0.8) +
  geom_rect(data=clean_data,aes(xmin=wl_center-(fwhm/2),xmax=wl_center+(fwhm/2),
                ymin=0, ymax=0.49),color="black",fill="red",alpha=0.5,size=1.1) +
  ## geom_vline(xintercept = 685, linetype="dashed", size=1.2) +
  ylab("") +
  xlab("\n Wavelength $\\lambda$ [nm]") +
  ## xlim(c(min(wave$wl_center),max(wave$wl_center))) +
  ylim(c(0,0.5)) +
  theme_bw() +
  theme(text = element_text(size=30),
        legend.position = "none",
        plot.title = element_text(hjust=0.5),
        axis.ticks.length = unit(0.2, "cm"),
        axis.ticks.y=element_blank(),
        axis.text.y=element_blank(),
        plot.margin = margin(10, 50, 10, 10)) +
  scale_x_continuous(limits=c(min(wave$wl_center),max(wave$wl_center)),
                     breaks = round(seq(roundUp(min(wave$wl_center),100),
                                        roundUp(max(wave$wl_center),100),
                                        by =100),1),
                     expand=expand_scale(mult=c(0.01,0.01))) +
  facet_wrap(Sensor~.,ncol=1)
print(g)
dev.off()
texi2pdf("basic_spectral_truncated_nos.tex",clean=TRUE)
file.remove("basic_spectral_truncated_nos.tex")
file.rename("basic_spectral_truncated_nos.pdf",
            "./img/basic_spectral_truncated_nos.pdf")