资源描述：
《ncl绘图示例（四）兰伯特投影图》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、NCL绘图示例（四）：兰伯特投影图施宁（南京信息工程大学大气科学学院）procedureadd_lc_labels(wks,map,minlat,maxlat,minlon,maxlon,fontheight)locallat_values,nlat,lat1_ndc,lat2_ndc,lon1_ndc,lon2_ndc,slope,txres,lon_values,PI,RAD_TO_DEG,dum_lft,dum_rgt,dum_botbeginPI=3.14159RAD_TO_DEG=180./PI;挑出“
2、较好”的纬度标签lat_values=ispan(toint(minlat),toint(maxlat),10)*1.nlat=dimsizes(lat_values);计算图形左右两边线条的斜率（基于NDC坐标）lat1_ndc=new(1,float)lon1_ndc=new(1,float)lat2_ndc=new(1,float)lon2_ndc=new(1,float)datatondc(map,minlon,lat_values(0),lon1_ndc,lat1_ndc)datatondc(map,mi
3、nlon,lat_values(nlat-1),lon2_ndc,lat2_ndc)slope_lft=(lat2_ndc-lat1_ndc)/(lon2_ndc-lon1_ndc)datatondc(map,maxlon,lat_values(0),lon1_ndc,lat1_ndc)datatondc(map,maxlon,lat_values(nlat-1),lon2_ndc,lat2_ndc)slope_rgt=(lat2_ndc-lat1_ndc)/(lon2_ndc-lon1_ndc)txres=Tru
4、etxres@txFontHeightF=fontheighttxres@txPosXF=0.1dum_lft=new(nlat,graphic);Dummyarraytoholdattachedstrings.dum_rgt=new(nlat,graphic);Dummyarraytoholdattachedstrings.don=0,nlat-1;添加适当的空白lat_label_rgt=""+lat_values(n)+"~S~o~N~";检查是否在北纬、南纬或赤道上if(lat_values(n).lt.0
5、)thenlat_label_lft=lat_values(n)+"~S~o~N~S"lat_label_rgt=lat_label_rgt+"S"endifif(lat_values(n).gt.0)thenlat_label_lft=lat_values(n)+"~S~o~N~N"lat_label_rgt=lat_label_rgt+"N"endifif(lat_values(n).eq.0)thenlat_label_lft=lat_values(n)+"~S~o~N~"endif;添加左边坐标标签txre
6、s@txAngleF=RAD_TO_DEG*atan(slope_lft)-90;适当旋转字体，更加美观dum_lft(n)=gsn_add_text(wks,map,lat_label_lft,minlon,lat_values(n),txres);添加右边坐标标签txres@txAngleF=RAD_TO_DEG*atan(slope_rgt)+90dum_rgt(n)=gsn_add_text(wks,map,lat_label_rgt,maxlon,lat_values(n),txres)enddo;---
7、-------------------------------------------------------------------;添加经度标签delete(txres@txPosXF)txres@txPosYF=-5.0;挑出“较好”的经度标签lon_values=ispan(toint(minlon+10),toint(maxlon-10),10)*1.nlon=dimsizes(lon_values)dum_bot=new(nlon,graphic)don=0,nlon-1;对于每个经度标签，计算其需旋转
8、的角度，以使字体更加美观。datatondc(map,lon_values(n)-0.25,minlat,lon1_ndc,lat1_ndc)datatondc(map,lon_values(n)+0.25,minlat,lon2_ndc,lat2_ndc)slope_bot=(lat1_ndc-lat2_ndc)/(lon1_ndc-lon2_ndc)tx
