在不完整的矩形网格点列表上将矩阵设置为一个

Question

我在一个文本文件中有一组数据，其中显示了 81x61 大小的矩形网格的坐标。每行显示一个网格点的经度和纬度。经度从 50.00 变为 80.00（61 个值），纬度从 30.00 变为 70.00（81 个值），顺序如下：

50.00  30.00
51.50  30.00
52.00  30.00
.
.
.
79.50  30.00
80.00  30.00
50.00  31.00
50.50  31.00
51.00  31.00
.
.
.
79.00  70.00
79.50  70.00
80.00  70.00

我还有另一个文本文件，其中包含一些来自上述矩形网格的随机坐标。

我想创建一个大小为 81x61 的矩阵，其元素为 0 和 1，其中 1 对应于第二个文本文件中的坐标。

如何在 Matlab 中编写代码？

---

我需要的小规模示例：

文本文件：

  1  1
  1  2
  1  3
  .
  .
  .
  4  3
  4  4
  4  5 

上述文本文件对应的矩形网格：

1,1  1,2  1,3  1,4  1,5
2,1  2,2  2,3  2,4  2,5
3,1  3,2  3,3  3,4  3,5
4,1  4,2  4,3  4,4  4,5

第二个文本文件：

1  1
1  3
2  4
2  5
3  4
4  1
4  5  

上述文本文件对应的矩阵：

1  0  1  0  0
0  0  0  1  1
0  0  0  1  0
1  0  0  0  1

Answer 1

我自己找到了方法；

假设：经度值最小值为 50.00，纬度值最小值为 30.00

%// First column of the files is Longitude and the second column is Latitude
  fr = fopen('second_file.txt','r');
  lon = textscan(fr,'%f %*[^\n]');
  lon = lon{:};
  fclose(fr);
  fr = fopen('second_file.txt','r');
  lat = textscan(fr,'%*f %f %*[^\n]');
  lat = lat{:};
  fclose(fr);
%// We know that the overall size of the target matrix is 81x61
  overall = zeros(81,61);
%// We assume that the total number of lines in the second file is 1000 (don't know if there is a built-in command to determine that!)
  for k = 1:1000
      i = int16(( lat(k) - 30.00 ) / 0.5 + 1);
      j = int16(( lon(k) - 50.00 ) / 0.5 + 1);
      overall(i,j) = 1;
  end

Answer 2

我想创建一个大小为 81x61 的矩阵，其中元素为 0 和 1 1s 对应于第二个文本文件中的坐标的方式。

这就是问题中的相关信息。答案是 Matlab 函数sub2ind(documentation)。它将 x、y 坐标列表转换为数组索引列表，然后您可以方便地将其设置为一个。

假设您已读取名为 second_file 的 Nx2 矩阵中的第二个文件的内容，以及您在变量 matrix_size (81x61) 中给出的结果矩阵的大小。然后你做：

x = second_file(:, 1);
y = second_file(:, 2);
result = zeros(matrix_size);
index = sub2ind(matrix_size, x, y);
result(index) = 1;

Answer 3

我假设您的第一个文件中各列中的minimum 和maximum 值是纬度和经度的范围。此外，在这两个文件中，第一列都是经度值。

 %// calculating the min and max of latitude and longitude 
 id1  = fopen('first_file.txt');
 A = fscanf(id1,'%f %f',[2 Inf]); 
 long_min = min(A(1,:));
 long_max = max(A(1,:));
 lat_min = min(A(2,:));
 lat_max = max(A(2,:));
 %// calculating output matrix size
 no_row = (lat_max-lat_min)*2+1;   
 no_col = (long_max-long_min)*2+1;
 output = zeros(no_row,no_col);

 id2  = fopen('second_file.txt');
 A = fscanf(id2,'%f %f',[2 Inf]);
 % // converting the values into corresponding indices
 A(1,:) = A(1,:) - long_min;     
 A(2,:) = A(2,:) - lat_min;   
 A = A*2 +1; 
 linear_ind = sub2ind([no_row no_col],A(2,:),A(1,:));
 output(linear_ind) = 1;

第二种方法 我假设在您的第二个文本文件中，第一列条目是纬度，第二列条目是经度。您必须对以下变量进行硬编码：

 long_min : the minimum value among longitude values
 lat_min : the minimum value among latitude values
 long_max : maximum value among longitude values
 lat_max : maximum value among latitude values

这是代码（仅考虑第二个文本文件）

 no_row = (lat_max-lat_min)*2+1;   
 no_col = (long_max-long_min)*2+1;
 output = zeros(no_row,no_col);
 id2  = fopen('second_file.txt');
 A = fscanf(id2,'%f %f',[2 Inf]);
 % // converting the values into corresponding indices
 A(1,:) = A(1,:) - long_min;     
 A(2,:) = A(2,:) - lat_min;   
 A = A*2 +1; 
 linear_ind = sub2ind([no_row no_col],A(2,:),A(1,:));
 output(linear_ind) = 1;