Julian/JuliaMSI
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Fixed scales on plots Topography, Surface and axis on spectrum

Browse Source
This commit is contained in:
Pixelguy14 2025-02-20 15:22:54 -06:00
parent 0456abb70b
commit 3714b404cd
2 changed files with 31 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
app.jl
View File

@ -168,11 +168,13 @@ include("./julia_imzML_visual.jl")
title="Spectrum plot",
xaxis=PlotlyBase.attr(
title="<i>m/z</i>",
showgrid=true
showgrid=true,
tickformat = ".3g"
),
yaxis=PlotlyBase.attr(
title="Intensity",
showgrid=true
showgrid=true,
tickformat = ".3g"
),
margin=attr(l=0,r=0,t=120,b=0,pad=0)
)

29
julia_imzML_visual.jl
View File

@ -253,6 +253,12 @@ function loadContourPlot(interfaceImg::String)
X=repeat(reshape(x, 1, length(x)), length(y), 1)
Y=repeat(reshape(y, length(y), 1), 1, length(x))
# Define tick values and text for colorbars
min_val = minimum(elevation_smoothed)
max_val = maximum(elevation_smoothed)
tickV = range(min_val, stop=max_val, length=8)
tickT = log_tick_formatter(collect(tickV))
layout=PlotlyBase.Layout(
title="2D topographic map of $cleaned_img",
xaxis=PlotlyBase.attr(
@ -271,7 +277,9 @@ function loadContourPlot(interfaceImg::String)
contours_coloring="Viridis",
colorscale="Viridis",
colorbar = attr(
tickformat=".2g"
tickvals = tickV,
ticktext = tickT,
tickmode = "array"
)
)
plotdata=[trace]
@ -291,15 +299,24 @@ function loadSurfacePlot(interfaceImg::String)
img_gray=Gray.(img) # Convert to grayscale
img_array=Array(img_gray)
elevation=Float32.(Array(img_array)) ./ 255.0 # Normalize between 0 and 1
# Smooth the image
sigma=3.0
kernel=Kernel.gaussian(sigma)
elevation_smoothed=imfilter(elevation, kernel)
# Create the X, Y meshgrid coordinates
x=1:size(elevation_smoothed, 2)
y=1:size(elevation_smoothed, 1)
X=repeat(reshape(x, 1, length(x)), length(y), 1)
Y=repeat(reshape(y, length(y), 1), 1, length(x))
# Define tick values and text for colorbars
min_val = minimum(elevation_smoothed)
max_val = maximum(elevation_smoothed)
tickV = range(min_val, stop=max_val, length=8)
tickT = log_tick_formatter(collect(tickV))
# Calculate the number of ticks and aspect ratio for the 3d plot
x_nticks=min(20, length(x))
y_nticks=min(20, length(y))
@ -337,7 +354,9 @@ function loadSurfacePlot(interfaceImg::String)
),
colorscale="Viridis",
colorbar = attr(
tickformat=".2g"
tickvals = tickV,
ticktext = tickT,
nticks=8
)
)
plotdata=[trace3D]
@ -418,11 +437,13 @@ function sumSpectrumPlot(mzmlRoute::String)
title="SUM spectrum plot",
xaxis=PlotlyBase.attr(
title="<i>m/z</i>",
showgrid=true
showgrid=true,
tickformat = ".3g"
),
yaxis=PlotlyBase.attr(
title="Intensity",
showgrid=true
showgrid=true,
tickformat = ".3g"
),
autosize=false,
margin=attr(l=0,r=0,t=120,b=0,pad=0)