data try8 right gs 2200 10000 100 1100 500 500.mat

creation history

cp -i Mouse.aba.R.In.FLAT_CYCLE5.2007-08-28.23460.coord Mouse.aba.R.In.FLAT_CYCLE5.2007-08-28.23460_text.coord
caret_file_convert -text Mouse.aba.R.In.FLAT_CYCLE5.2007-08-28.23460_text.coord
perl -e 'undef $/; $file = <>; $file =~ s/.*EndHeader\n\d+\n//s; print $file' < Mouse.aba.R.In.FLAT_CYCLE5.2007-08-28.23460_text.coord > flat_cycle5.txt

Now, in Octave,

A = load('try8_right_gs_2200_10000_100_1100_500_500.txt');
sampleDist = 100; % microns
flatcycle5 = load('flat_cycle5.txt');

[XI, YI, ZI] = my_griddata(flatcycle5(:,2),flatcycle5(:,3),flatcycle5(:,1),min(flatcycle5(:,2)):sampleDist:max(flatcycle5(:,2)),min(flatcycle5(:,3)):sampleDist:max(flatcycle5(:,3)),'nearest',sampleDist);

% test
D = A(ZI+1,1);
D= reshape(D,size(ZI));
max(max(D - ZI))
%ans = 0


D= zeros([size(ZI) size(A,2)-1]);
 for gene=2:size(A,2)
  D(:,:,gene-1) = A(ZI + 1,gene);
  D(find(0 == ZI) + (gene-2).*prod(size(ZI))) = NaN;
end

% test
im = D(:,:,1); im(isnan(im)) = 0; imagesc(flipud(im))

paintFile = load('paint.txt');
paint = paintFile(ZI + 1, 2);
paint = reshape(paint, size(ZI));
paint(find(0 == ZI)) = NaN;

% test
im = paint; im(isnan(im)) = 0; imagesc(flipud(im))

% test
for i=1:30
  %printf('%d: %d\n', i, length(find(paint == i)))
  j(i) = length(find(paint == i));
end
find(j == 0)

%ans =
%    9   10   11

clear i j

%save 'data_try7_right_gs_2200_10000_300_1100_500_1100.octave_file' D
save -mat-binary 'data_try8_right_gs_2200_10000_100_1100_500_500.mat' D



notNaNIndices = find(~isnan(paint(:)));
paint2 = paint(:);
paint3 = paint2(notNaNIndices);
save('-mat', 'paint.mat', 'paint', 'paint2', 'paint3')

D2 = zeros(size(D,1)*size(D,2),size(D,3));

for i=1:size(D,3)
  dd = D(:,:,i);
  D2(:,i) = dd(:);
  %i
end

notNaNIndices = find(~isnan(D2(:,1)));
D3 = D2(notNaNIndices,:);

D_valid = ~isnan(D(:,:,1));
D_NaNs = isnan(D(:,:,1)).*D(:,:,1); 

D_NaNs_and_ones = D_NaNs + 1; 
D_vert = repmat((1:size(D,1))',1,size(D,2)) .* D_NaNs_and_ones; 
D_horiz = repmat((1:size(D,2)),size(D,1),1) .* D_NaNs_and_ones; 
D_vert2 = D_vert(:); 
D_horiz2 = D_horiz(:); 

% test
max(find(~isnan(D_vert2)) - find(~isnan(D_horiz2))) 
% ans = 0
notNaNIndices = find(~isnan(D_vert2));

D_vert3 = D_vert2(notNaNIndices,:);
D_horiz3 = D_horiz2(notNaNIndices,:);

% test
testColumn = D3(:,1);
im = [];
for i=1:size(testColumn,1)
 im(D_vert3(i),D_horiz3(i),:) = testColumn(i,:);
end
imagesc(flipud(im))


save ('-mat','data_try8_right_gs_2200_10000_100_1100_500_500_start.mat', 'D', 'D2', 'D3', 'D_valid', 'D_NaNs', 'D_NaNs_and_ones', 'D_vert', 'D_horiz', 'D_vert2', 'D_horiz2', 'D_vert3','D_horiz3')

See also [[Mouse.aba.R.In.FLAT_CYCLE5.2007-08-28.23460.coord?]]

CategoryDataSetFile