宁夏FC计算（2）

% 从表格提取经纬度和数值列，如果不知道行列好的情况下，通过定义分辨率
lat = oidata_B.lat;
lon = oidata_B.lon;
FC = oidata_B.FC;

% 定义输出栅格的分辨率
resolution = 0.01; % 例如，每个像素0.01度
latlim = [min(lat) max(lat)]; % 纬度范围
lonlim = [min(lon) max(lon)]; % 经度范围

% 创建规则的经纬度网格
[lonGrid, latGrid] = meshgrid(lonlim(1):resolution:lonlim(2), ...
                              latlim(1):resolution:latlim(2));

% 使用 griddata 将离散点插值到规则的网格上
FC_grid = griddata(lon, lat, FC, lonGrid, latGrid, 'linear');

% 创建空间参考对象，定义TIFF文件的地理信息
R = georefcells(latlim, lonlim, size(FC_grid));

% 保存为 GeoTIFF 文件
geotiffwrite('output_FC.tif', FC_grid, R);

disp('TIFF文件已成功生成');
 

% 从表格提取经纬度和数值列，这里知道原始数据的行列号，基本就是确定原始数据的分辨率大小。
lat = oidata_B.lat;
lon = oidata_B.lon;
FC = oidata_B.FC;

% 定义输出的行列数
rows = 558;
cols = 380;

% 定义纬度和经度范围
latlim = [min(lat) max(lat)]; % 纬度范围
lonlim = [min(lon) max(lon)]; % 经度范围

% 创建规则的经纬度网格
[lonGrid, latGrid] = meshgrid(linspace(lonlim(1), lonlim(2), cols), ...
                               linspace(latlim(1), latlim(2), rows));

% 使用 griddata 将离散点插值到规则的网格上
FC_grid = griddata(lon, lat, FC, lonGrid, latGrid, 'linear');

% 创建空间参考对象，定义TIFF文件的地理信息
R = georefcells(latlim, lonlim, size(FC_grid));

% 定义输出路径
output_path = 'F:\a\田间持水量\实验\output_FC.tif';

% 保存为 GeoTIFF 文件
geotiffwrite(output_path, FC_grid, R);

disp('TIFF文件已成功生成');