利用shp文件编写想要经纬格点数的地形nc文件

无敌大壮壮

1. 
   
2. load "$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/cnmap.ncl"
   
3. load "$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/shapefile_utils.ncl"
   
4. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
   
5. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
   
6. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
   
7. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
   
8. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
   
9. 
   
10. glon1=fspan(90.25,124.75,70);定义一个浮点型数据，从90-112，分为221个格点
    
11. glat1=fspan(23.75,37.25,28);纬度的设置，注意后面的格点+1
    
12. ; print(glon1)
    
13. ; exit
    
14. glon1!0="lon"
    
15. glon1@long_name="lon"
    
16. glon1@units="degrees_east"
    
17. glon1&lon=glon1
    
18. 
    
19. glat1!0="lat"
    
20. glat1@long_name="lat"
    
21. glat1@units="degrees_north"
    
22. glat1@lat=glat1
    
23. 
    
24. 
    
25. nlat=dimsizes(glat1)
    
26. nlon=dimsizes(glon1)
    
27. data=generate_2d_array(10, 10, 0., 100., 0, (/nlat,nlon/))
    
28. 
    
29. ; printVarSummary(data)
    
30. ; exit
    
31. 
    
32. data!0 = "lat"
    
33. data!1 = "lon"
    
34. data&lat = glat1
    
35. data&lon = glon1
    
36. shp_filename ="D:/GIS/graduation/YangtzeRiverArea/YangtzeRiverArea.shp"
    
37. 
    
38. setfileoption("nc","Format","LargeFile")
    
39. fout =addfile("YZ.nc", "c")
    
40. data_mask = shapefile_mask_data(data,shp_filename,True)
    
41. fout->data_mask=data_mask
    
42. printVarSummary(data_mask)
    

复制代码

以上内容是根据csdn这位大佬的内容做的改写   

https://blog.csdn.net/rqjabc/article/details/111051087?ops_request_misc=&request_id=&biz_id=102&utm_term=NCL%20%E7%BC%96%E5%86%99%E4%B8%AD%E5%9B%BD%E5%8C%BA%E5%9F%9Fnc&utm_medium=distribute.pc_search_result.none-task-blog-2~all~sobaiduweb~default-4-111051087.142^v87^control,239^v2^insert_chatgpt&spm=1018.2226.3001.4187

data=generate_2d_array(10, 10, 0., 100., 0, (/nlat,nlon/))   代码中的这个部分可以不用改变  

要改的地方就是这三个地方

游客，如果您要查看本帖隐藏内容请回复

修改完这些就可以得到想要的地形nc文件啦！

lijian

谢谢楼主的分享，感谢！

流浪的心

感谢分享！

随缘

谢谢楼主的分享

3two

感谢大佬分享

yinbolan

厉害厉害厉害厉害厉害厉害厉害

hagdon

非常感谢您的分享，仔细研究一下~

soso

非常有用，在这个基础上可以继续加减

blessing

kanyikan qiaoyiqiao

饼饼

感谢lz感谢lz感谢lz感谢lz