编辑-地球科学江南体育网页版堆栈交换江南电子竞技平台

你不是登录。你的编辑将被放置在一个队列,直到同行评议。

我们欢迎编辑,使读者更容易理解和更有价值。因为社区成员审查编辑,请尽量使文章更好比你怎么发现,例如,通过修正语法或添加额外的资源和超链接。

牧师

标题

身体

我致力于MOD10C2积雪(8天),当我夹MODIS积雪对不同区域的海拔高度对.i.e总趋势的结果。高海拔低显示积雪和低海拔显示更高的积雪。我用matlab inpolygon函数(inpolygon (yq xq,十五,青年志愿))我无法找出犯错误。首先,我以这种方式导入数据hdfvars = {‘Eight_Day_CMG_Snow_Cover’,‘Eight_Day_CMG_Clear_Index’,……“Eight_Day_CMG_Cloud_Obscured”、“Snow_Spatial_QA '};测试管理扫描= ' E: \ \ hdf测试”;dinfo = dir (fullfile(管理扫描' * .hdf '));num_files =长度(dinfo);文件名= fullfile(管理扫描,{dinfo.name});Eight_Day_CMG_Snow_Cover =细胞(num_files, 1); Eight_Day_CMG_Clear_Index = cell(num_files, 1); Eight_Day_CMG_Cloud_Obscured = cell(num_files, 1); Snow_Spatial_QA = cell(num_files, 1); for K = 1 : num_files this_file = filenames{K}; Eight_Day_CMG_Snow_Cover{K} = hdfread(this_file, hdfvars{1}); Eight_Day_CMG_Clear_Index{K} = hdfread(this_file, hdfvars{2}); Eight_Day_CMG_Cloud_Obscured{K} = hdfread(this_file, hdfvars{3}); Snow_Spatial_QA{K} = hdfread(this_file, hdfvars{4}); end Then reshaped data in this way B2 = zeros(3600,7200,24); for i = 1:3600 for j = 1:7200 B2(i,j,1:24) = reshape(Eight_Day_CMG_Snow_Cover{i,j},[1 3 2]); end end Then generate lat lon and subset the area of interest lon = -180:0.05:180; lat = -90:0.05:90; subsetqa = Eight_Day_CMG_Snow_Cover(2001:2817,4801:5741,:); then attempted to extract with elevation polygon, I tried both [monthly] and [8days] fractional snow cover data, elevation polygon are extracted from SRTM DEM. I attempted to use matlab in this way shapefile = 'shapefile.shp' ; S = shaperead(shapefile) ; N = length(S) ; for i = 1:N plot(S(i).X,S(i).Y) hold on end %% lon = load('testlon.mat') ; lon = lon.testlon ; lat = load('testlat.mat') ; lat = lat.testlat ; [X,Y] = meshgrid(lon,lat) ; data = load('testarray.mat') ; data = data.testarray ; [nx,ny,d] = size(data) ; %%Extract data iwant = cell(d,N) ; for i =1:d A = data(:,:,i) ; for j = 1:N idx = inpolygon(X(:),Y(:),S(i).X,S(i).Y) ; iwant{i,j} = A(idx) ; end Then output iwant converted into matrix in this way test = cell2mat(cellfun(@transpose,iwant,'uniform',0)); NOTE: this picture is showing polygon of only one elevation zone and we have six different elevation zones. [![enter image description here][1]][1] [1]: https://i.stack.imgur.com/UOCLB.jpg

标签

编辑总结

取消

1

\ begingroup美元你能添加更多的细节?如何生成海拔区poligon吗?什么你的源头高程数据?你也可以把代码和数据文件的链接你在工作,所以我们可以更好的帮助。首先我认为,一个数据集不正确定位,或高程数据的分辨率太粗了。 \ endgroup美元
- - - - - -卡米洛·Rada
2019年2月2日15:43
\ begingroup美元亲爱的卡米洛·拉扎,问题是编辑和添加了所有要求的事情。寻找你的宝贵的建议和帮助 \ endgroup美元
- - - - - -伊尔凡
2019年2月4日在三
\ begingroup美元看起来不完整的代码。“tp5k6k”是什么?你能添加的图,结果:图;显示亮度图像(数据(:,:1));抓住;情节((1)方式,年代(1).Y, ' r '); \ endgroup美元
- - - - - -卡米洛·Rada
2019年2月4日3
\ begingroup美元更新我的问题。请看看 \ endgroup美元
- - - - - -伊尔凡
2019年2月4日3:31
1

\ begingroup美元这些不是文件的脚本。我的意思是shapefile。轴马力,testlon。垫,testlat。垫子和testarray.mat \ endgroup美元
- - - - - -卡米洛·Rada
2019年2月4日5:10

| 显示6更多的评论

正确的小错误或错误
澄清的意思,而不更改它
添加相关资源或链接
总是尊重作者的意图
不使用编辑与作者回复