NetCDF数据集（.nc）怎么形成矢量图？

tltb09125006

本人是遥感方面的研究生，最近导师要我下载NECP的气象数据，我下到的都是NetCDF数据集（文件后缀名为.nc），之前几天学习了怎么读取该数据，在grads中也可以成图。提取出研究区域的数据后，在ArcGIS中可以通过插值形成矢量图。但是会有好多年的数据需要这样操作。
想求问大家，这个数据有什么其他方式可以，更直接的形成矢量图吗？

MeteoInfo

矢量图指的什么？

tltb09125006

MeteoInfo 发表于 2015-4-17 23:26
矢量图指的什么？

就是在ArcGIS或者ENVI中能读取出来，参与栅格计算或者波段运算的格式

kongfeng0824

以下的代码可把NC输出成ArcGIS支持的tif格式：

1. clc %清屏
   
2. clear %清内存变量
   
3. cd('F:\Data');
   
4. % filename1='PrepCN_1961.nc';
   
5. % ncdisp(filename1)
   
6. % longitude=ncread(filename1,'lon');
   
7. % latitude =ncread(filename1,'lat');
   
8. % precipitation=ncread(filename1,'prep');
   
9. % % precipitation=(1.467e-06).*precipitation+0.048066;
   
10. % precipitation=precipitation.*1000;
    
11. % precipitation(precipitation<50)=0;
    
12. % aa=sum(precipitation(:,:,:),3);
    
13. % pcolor(aa')
    
14. 
    
15. filenames = ls('*.nc');%以字符串数组的形式列出该路径下所有的txt文件
    
16. [m n] = size(filenames);%获取filenames的行列数
    
17. for j = 1:m
    
18. precipitation=ncread(char(filenames(j,:)),'prep');
    
19. precipitation(precipitation(:,:,:)==-999)=0;
    
20. precipitation(precipitation(:,:,:)<25)=0;
    
21. precipitation(precipitation(:,:,:)<50 &precipitation(:,:,:)>=25)=1;
    
22. precipitation(precipitation(:,:,:)>=50)=0;
    
23. aa=sum(precipitation(:,:,:),3);
    
24. clear precipitation;
    
25. Rain=double(aa');
    
26. clear aa;
    
27. 
    
28. ncname=char(filenames(j,1:11));
    
29. R = georasterref('RasterInterpretation','postings');
    
30. R.RasterSize = [69 126];
    
31. R.LatitudeLimits  = [19.75 53.75];
    
32. R.LongitudeLimits = [73.25 135.75];
    
33. R.ColumnsStartFrom='north';
    
34. R.RowsStartFrom='west';
    
35. geotiffwrite(ncname,Rain,R)
    
36. JDT=[ncname,'转换完毕!'];
    
37. disp(JDT)
    
38. end

复制代码

以实现

kongfeng0824

以下的代码可把NC输出成ArcGIS支持的tif格式：

1. clc %清屏
   
2. clear %清内存变量
   
3. cd('F:\Data');
   
4. % filename1='PrepCN_1961.nc';
   
5. % ncdisp(filename1)
   
6. % longitude=ncread(filename1,'lon');
   
7. % latitude =ncread(filename1,'lat');
   
8. % precipitation=ncread(filename1,'prep');
   
9. % % precipitation=(1.467e-06).*precipitation+0.048066;
   
10. % precipitation=precipitation.*1000;
    
11. % precipitation(precipitation<50)=0;
    
12. % aa=sum(precipitation(:,:,:),3);
    
13. % pcolor(aa')
    
14. 
    
15. filenames = ls('*.nc');%以字符串数组的形式列出该路径下所有的txt文件
    
16. [m n] = size(filenames);%获取filenames的行列数
    
17. for j = 1:m
    
18. precipitation=ncread(char(filenames(j,:)),'prep');
    
19. precipitation(precipitation(:,:,:)==-999)=0;
    
20. precipitation(precipitation(:,:,:)<25)=0;
    
21. precipitation(precipitation(:,:,:)<50 &precipitation(:,:,:)>=25)=1;
    
22. precipitation(precipitation(:,:,:)>=50)=0;
    
23. aa=sum(precipitation(:,:,:),3);
    
24. clear precipitation;
    
25. Rain=double(aa');
    
26. clear aa;
    
27. 
    
28. ncname=char(filenames(j,1:11));
    
29. R = georasterref('RasterInterpretation','postings');
    
30. R.RasterSize = [69 126];
    
31. R.LatitudeLimits  = [19.75 53.75];
    
32. R.LongitudeLimits = [73.25 135.75];
    
33. R.ColumnsStartFrom='north';
    
34. R.RowsStartFrom='west';
    
35. geotiffwrite(ncname,Rain,R)
    
36. JDT=[ncname,'转换完毕!'];
    
37. disp(JDT)
    
38. end

复制代码

以实现

kongfeng0824

tltb09125006 发表于 2015-4-18 10:20
就是在ArcGIS或者ENVI中能读取出来，参与栅格计算或者波段运算的格式

以下的代码可把NC输出成ArcGIS支持的tif格式：
clc %清屏
clear %清内存变量
cd('F:\Data');
% filename1='PrepCN_1961.nc';
% ncdisp(filename1)
% longitude=ncread(filename1,'lon');
% latitude =ncread(filename1,'lat');
% precipitation=ncread(filename1,'prep');
% % precipitation=(1.467e-06).*precipitation+0.048066;
% precipitation=precipitation.*1000;
% precipitation(precipitation<50)=0;
% aa=sum(precipitation(:,:,:),3);
% pcolor(aa')

filenames = ls('*.nc');%以字符串数组的形式列出该路径下所有的txt文件
[m n] = size(filenames);%获取filenames的行列数
for j = 1:m
precipitation=ncread(char(filenames(j,:)),'prep');
precipitation(precipitation(:,:,:)==-999)=0;
precipitation(precipitation(:,:,:)<25)=0;
precipitation(precipitation(:,:,:)<50 &precipitation(:,:,:)>=25)=1;
precipitation(precipitation(:,:,:)>=50)=0;
aa=sum(precipitation(:,:,:),3);
clear precipitation;
Rain=double(aa');
clear aa;

ncname=char(filenames(j,1:11));
R = georasterref('RasterInterpretation','postings');
R.RasterSize = [69 126];
R.LatitudeLimits  = [19.75 53.75];
R.LongitudeLimits = [73.25 135.75];
R.ColumnsStartFrom='north';
R.RowsStartFrom='west';
geotiffwrite(ncname,Rain,R)
JDT=[ncname,'转换完毕!'];
disp(JDT)
end

王小燕1991

学习中，正在研究中呢，共同学习

暖气团

kongfeng0824 发表于 2015-4-18 12:50
以下的代码可把NC输出成ArcGIS支持的tif格式：
以实现

如果将非NECP的NetCDF数据转成tif，怎么做呢？
kongfeng0824

周雪

共同学习  共同进步

yinshao

楼主，用那个代码问题解决了么