Skip to content
Snippets Groups Projects
Commit 580833ed authored by Layton Aho's avatar Layton Aho
Browse files

Remove kymo language

parent a02fb956
No related branches found
No related tags found
No related merge requests found
...@@ -54,7 +54,7 @@ To create a new synethetic dataset: ...@@ -54,7 +54,7 @@ To create a new synethetic dataset:
- phantomVelocity_mu: the different mean velocity values per class i.e., the number of means determines the number of classes. - phantomVelocity_mu: the different mean velocity values per class i.e., the number of means determines the number of classes.
- phantomVelocity_sigma: the standard deviation for velocity per class. - phantomVelocity_sigma: the standard deviation for velocity per class.
- numOutput: how many phantom images to generate per class. - numOutput: how many phantom images to generate per class.
2. Run phantomVelocitySS\goFigure.m, specifying the saved .mat file as the input. 2. Run phantomVelocitySSF\goFigure.m, specifying the saved .mat file as the input.
## Notes ## Notes
1. The NCD code uses the [FLIF](https://flif.info/) compression library. The thunker file for the FLIF library was generated on Ubuntu 22.04 and **may not work** on other Unix platforms. 1. The NCD code uses the [FLIF](https://flif.info/) compression library. The thunker file for the FLIF library was generated on Ubuntu 22.04 and **may not work** on other Unix platforms.
......
function phantomKymo = createPhantom(phSize,phTracks,phVel_mu,phVel_sigma) function phantomVelcocityOutput = createPhantom(phSize,phTracks,phVel_mu,phVel_sigma)
tmp = zeros(phSize,'uint8'); tmp = zeros(phSize,'uint8');
...@@ -39,6 +39,6 @@ for i = 1:length(trackSrcPts) ...@@ -39,6 +39,6 @@ for i = 1:length(trackSrcPts)
end end
phantomKymo = tmp; phantomVelcocityOutput = tmp;
end end
\ No newline at end of file
...@@ -10,7 +10,7 @@ for i = 1:length(datasetParameters.phantomVelocity_mu) ...@@ -10,7 +10,7 @@ for i = 1:length(datasetParameters.phantomVelocity_mu)
phVel_mu = datasetParameters.phantomVelocity_mu(i); phVel_mu = datasetParameters.phantomVelocity_mu(i);
phVel_sigma = datasetParameters.phantomVelocity_sigma(i); phVel_sigma = datasetParameters.phantomVelocity_sigma(i);
for j = 1:datasetParameters.numKymos for j = 1:datasetParameters.numOutput
tmp = []; tmp = [];
tmp.phantomName = strcat(num2str(i),'_',num2str(j)); tmp.phantomName = strcat(num2str(i),'_',num2str(j));
...@@ -19,7 +19,7 @@ for i = 1:length(datasetParameters.phantomVelocity_mu) ...@@ -19,7 +19,7 @@ for i = 1:length(datasetParameters.phantomVelocity_mu)
tmp.phantomSize = phSize; tmp.phantomSize = phSize;
tmp.numTracks = phTracks; tmp.numTracks = phTracks;
tmp.phantomVelocity = [phVel_mu;phVel_sigma]; tmp.phantomVelocity = [phVel_mu;phVel_sigma];
tmp.phantomKymo = createPhantom(phSize,phTracks,phVel_mu,phVel_sigma); tmp.phantomVelocityOutput = createPhantom(phSize,phTracks,phVel_mu,phVel_sigma);
syntheticDataset = [syntheticDataset;tmp]; syntheticDataset = [syntheticDataset;tmp];
......
function clipLimits = getPhantomClipLimits(syntheticDataset,targetChannel) function clipLimits = getPhantomClipLimits(syntheticDataset,targetChannel)
for ff = 1:length(syntheticDataset) for ff = 1:length(syntheticDataset)
im_ff = syntheticDataset(ff).phantomKymo(:,:,:,targetChannel); im_ff = syntheticDataset(ff).phantomVelocityOutput(:,:,:,targetChannel);
kp = im_ff(find(im_ff)); kp = im_ff(find(im_ff));
......
...@@ -2,9 +2,9 @@ load('sampleDataset_01.mat'); ...@@ -2,9 +2,9 @@ load('sampleDataset_01.mat');
K = 3; K = 3;
tmp_xy(:,:,1) = squeeze(max(syntheticDataset(1).phantomKymo(:,:,:,1),[],2)); tmp_xy(:,:,1) = squeeze(max(syntheticDataset(1).phantomVelocityOutput(:,:,:,1),[],2));
tmp_xy(:,:,2) = squeeze(max(syntheticDataset(11).phantomKymo(:,:,:,1),[],2)); tmp_xy(:,:,2) = squeeze(max(syntheticDataset(11).phantomVelocityOutput(:,:,:,1),[],2));
tmp_xy(:,:,3) = squeeze(max(syntheticDataset(21).phantomKymo(:,:,:,1),[],2)); tmp_xy(:,:,3) = squeeze(max(syntheticDataset(21).phantomVelocityOutput(:,:,:,1),[],2));
figure('Position',[300,200,900,300]); figure('Position',[300,200,900,300]);
layout1 = tiledlayout(1,3,'TileSpacing','compact','Padding','tight'); layout1 = tiledlayout(1,3,'TileSpacing','compact','Padding','tight');
......
...@@ -6,7 +6,7 @@ synthDatasetParameters.phantomSize = [200,200,25]; ...@@ -6,7 +6,7 @@ synthDatasetParameters.phantomSize = [200,200,25];
synthDatasetParameters.phantomTracks = 10; synthDatasetParameters.phantomTracks = 10;
synthDatasetParameters.phantomVelocity_mu = [1,3,5]; synthDatasetParameters.phantomVelocity_mu = [1,3,5];
synthDatasetParameters.phantomVelocity_sigma = [.5,.5,.5]; synthDatasetParameters.phantomVelocity_sigma = [.5,.5,.5];
synthDatasetParameters.numKymos = 50; synthDatasetParameters.numOutput = 50;
syntheticDataset = createSyntheticDataset(synthDatasetParameters); syntheticDataset = createSyntheticDataset(synthDatasetParameters);
......
...@@ -21,14 +21,14 @@ parfor i=1:H ...@@ -21,14 +21,14 @@ parfor i=1:H
continue continue
end end
im_ff = syntheticDataset(ff).phantomKymo(:,:,:,targetChannel); im_ff = syntheticDataset(ff).phantomVelocityOutput(:,:,:,targetChannel);
if ff == gg if ff == gg
d1 = SSF.ncd_ssf_volume(im_ff,im_ff); d1 = SSF.ncd_ssf_volume(im_ff,im_ff);
dxx(i,j) = d1; dxx(i,j) = d1;
else else
im_ff = syntheticDataset(ff).phantomKymo(:,:,:,targetChannel); im_ff = syntheticDataset(ff).phantomVelocityOutput(:,:,:,targetChannel);
im_gg = syntheticDataset(gg).phantomKymo(:,:,:,targetChannel); im_gg = syntheticDataset(gg).phantomVelocityOutput(:,:,:,targetChannel);
d1 = SSF.ncd_ssf_volume(im_ff,im_gg); d1 = SSF.ncd_ssf_volume(im_ff,im_gg);
d2 = SSF.ncd_ssf_volume(im_gg,im_ff); d2 = SSF.ncd_ssf_volume(im_gg,im_ff);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment