X，Y轴坐标自定义求助！！！

小其其格

• 文章有些图是用MeteoInfo绘制的，最后见刊排版阶段，编辑部要求X坐标的经度数值只在最后
• 
  
• 显示°E，Y坐标只在最后数值显示°N。按照MeteoInfo官网对xticks函数的说明修改的：
  

   xticks(arange(118.5, 122, 1.0), ['118.5','119.5','120.5','121.5E'],fontsize=fontSize, fontname='Times New Roman')

• 但是怎么也出不来结果，求大神指点！！！！
  

    @MeteoInfo

                               
登录/注册后可看大图

                               
登录/注册后可看大图

绘图代码：

1. 
   
2. dataDir = 'E:/wfl/MI/fig2/a~b/'
   
3. fileType = '.bin'
   
4. 
   
5. #字体大小
   
6. fontSize = 36
   
7. 
   
8. import os,sys,datetime
   
9. meteoInfoLabName = os.path.abspath(sys.argv[0])
   
10. print(meteoInfoLabName)
    
11. 
    
12. #读取指定文件夹中的指定文件类型的文件名
    
13. def get_filename(path,fileType):
    
14.     file_name =[]
    
15.     final_name = []
    
16.     for files in os.listdir(path):  #root为目录路径 #dirs为路径下的子目录 #files为路径下的所有非子目录
    
17.         if fileType in files:
    
18.            file_name.append(files.replace(fileType,''))#生成不带‘.bin’后缀的文件名组成的列表
    
19.            final_name = [path +item +fileType for item in file_name]#生成‘.bin’后缀的文件名组成的绝对路径列表
    
20.     return final_name,file_name #输出列表
    
21. 
    
22. 
    
23. 
    
24. def degree_1km(lon, lat):
    
25.     '''
    
26.     Get degree per kilometer at a location.   
    
27.     '''
    
28.     dis_km = distance([lon,lon+1], [lat,lat], islonlat=True) / 1000
    
29.     return 1 / dis_km
    
30. 
    
31. #Location of radar (Longitude/Latitude)
    
32. lon = 121.514
    
33. lat = 30.070
    
34. 
    
35. minlat=28.0
    
36. minlon=118.5
    
37. maxlat=30.5
    
38. maxlon=122.0
    
39. 
    
40. #Get degree per kilometer at the location
    
41. deg_1km = degree_1km(lon, lat)
    
42. 
    
43. 
    
44. 
    
45. 
    
46. #读取多个文件的绝对路径和文件名
    
47. fns,files = get_filename(dataDir,fileType)  #fns为带.bin的文件名的文件列表
    
48. nFiles = len(fns) #文件数量
    
49. print(nFiles)
    
50. print(fns)
    
51. print('\n')
    
52. print(files)
    
53. 
    
54. 
    
55. #读取城市经纬度和名称等信息
    
56. #Read city data of 11 stations from data file
    
57. datafn = dataDir+'CityVars.Zhejiang.txt'
    
58. ncol = numasciicol(datafn)
    
59. nrow = numasciirow(datafn)
    
60. var = asciiread(datafn, shape=(nrow,ncol))
    
61. var_city = var[:,2]
    
62. #Read station name and lon/lat
    
63. stfn = dataDir+'StationNames.Zhejiang.csv'
    
64. table = readtable(stfn, delimiter=',', format='%i%2s%2f',encoding='gb2312')
    
65. stnames = table['Name']
    
66. lat_city = table['LAT']
    
67. lon_city = table['LON']
    
68. 
    
69. 
    
70. #Read Radar datas and plot
    
71. for k in range(nFiles):
    
72.    
    
73.     f = addfile(fns[k])
    
74.     scan = 1
    
75.     rf = f['Reflectivity'][scan,:,:]
    
76.     azi = f['azimuthR'][scan,:]
    
77.     dis = f['distanceR'][:]/1000.0
    
78.     ele = f['elevationR'][scan,:]
    
79. 
    
80.     #Get x/y (kilometers) coordinates of data
    
81.     e = radians(ele)
    
82.     azi = 90 - azi
    
83.     a = radians(azi)
    
84.     nr = rf.shape[0]
    
85.     nd = rf.shape[1]
    
86.     x = zeros((nr + 1, nd))
    
87.     y = zeros((nr + 1, nd))
    
88.     for i in xrange (len(e)):
    
89.         x[i,:] = dis * cos(e[i]) * cos(a[i])
    
90.         y[i,:] = dis * cos(e[i]) * sin(a[i])
    
91.     x[nr,:] = x[0,:]
    
92.     y[nr,:] = y[0,:]
    
93.     rf1 = zeros((nr + 1, nd))
    
94.     rf1[:nr,:nd] = rf
    
95.     rf1[nr,:] = rf[0,:]
    
96. 
    
97.     #km to degree
    
98.     x = x * deg_1km + lon
    
99.     y = y * deg_1km + lat
    
100. 
     
101.     #Plot        
     
102.     #ax = axesm(bgcolor='b')
     
103.     ax = axesm()
     
104.     #Add map layers
     
105.     geoshow('cn_province', edgecolor=None, facecolor=[230,230,230])
     
106.     geoshow('cn_province', edgecolor=[80,80,80])
     
107.    
     
108.     #Plot radar reflectivity
     
109.     levs=[0,5,10,15,20,25,30,35,40,45,50,55,60,65]
     
110.     #color bar of CINRAD
     
111.     cols=[(255,255,255,0),(0,236,238),(0,161,244),(0,0,246),(0,255,0),\
     
112.         (0,200,0),(0,145,0),(255,255,0),(230,193,0),(254,146,0),\
     
113.         (254,0,0),(215,0,0),(192,0,0),(255,0,254),(155,85,199)]
     
114.     layer=pcolorm(x, y, rf1, levs, colors=cols, zorder=1)
     
115.     colorbar(layer, orientation='horizontal',shrink=0.8,fontsize=28,fontname=u'Times New Roman',label='dBZ', labelloc='top')
     
116.     #Plot circles
     
117.     rr = array([50, 100, 150, 200])
     
118.     for r in rr:
     
119.         rd = r * deg_1km
     
120.         ax.add_circle((lon, lat), rd, edgecolor='r')
     
121.     geoshow([lat,lat], [lon-rd,lon+rd], color='r')
     
122.     geoshow([lat+rd,lat-rd], [lon,lon], color='r')
     
123. 
     
124.     #绘制城市名称和描点
     
125.     layer_city = scatterm(lon_city, lat_city, var_city, size=5)
     
126.     layer_city.addfield('Name', 'string', stnames)
     
127.     layer_city.addlabels('Name', fontname=u'楷体', fontsize=fontSize, yoffset=-8)
     
128.    
     
129.     #Set plot
     
130.     if k ==3:
     
131.       minlat=28.0
     
132.       minlon=120.0
     
133.       maxlat=30.5
     
134.       maxlon=123.5
     
135.     xaxis(linewidth = 1.6) #坐标轴粗细
     
136.     yaxis(linewidth = 1.6) #坐标轴粗细
     
137.     #xlim(minlon, maxlon)
     
138.     ylim(minlat, maxlat)
     
139.     yticks(arange(minlat, maxlat,0.5), fontsize=fontSize, fontname='Times New Roman') #arange( start, stop, step )
     
140.     xticks(arange(118.5, 122, 1.0), ['118.5','119.5','120.5','121.5E'],fontsize=fontSize, fontname='Times New Roman')
     
141.    
     
142. 
     
143.     #保存图像
     
144.     #savefig(dataDir+files[k]+'.png', width=1920, height=1080, dpi=600)
     
145.     savefig(dataDir+files[k]+'.eps', width=1920, height=1080, dpi=600)
     
146.     print('Drawing ' +files[k] +'done!')
     
147.     if k < nFiles-1:
     
148.         cla() #清除，不然title要重叠
     
149. 
     
150. 
     
151.    
     
152. 
     

复制代码

出来的图如下，X坐标仍旧是每个数值后面有°E……

                               
登录/注册后可看大图

小其其格

已经解决了，怎么删除帖子呀