MeteoInfo脚本教程（十五）：地图投影

MeteoInfo

MeteoInfo中地图投影的代码是从Proj4移植来的（Java版用的Proj4J库，代码已经放入了MeteoInfo库中，在软件lib目录中看不到Prj4J库的文件），投影之间的转换需要设置Proj4的投影字符串。比如Lambert Conformal投影表述为：



+proj=lcc   +lat_1=Latitude of first standard parallel
              +lat_2=Latitude of second standard parallel
              +lat_0=Latitude of false origin
              +lon_0=Longitude of false origin
              +x_0=False Origin Easting
              +y_0=False Origin Northing



具体的设置可以看Proj4的帮助文档，一些网页上也有描述，如： http://remotesensing.org/geotiff/proj_list/。

投影时先建立投影字符串，然后创建ProjectionInfo类的对象，再用MapView类的projectLayers(ProjectionInfo projInfo)即可完成投影功能。

下面是两个示例脚本程序：

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2014-11-22
   
4. # Purpose: Read GrADS data and plot
   
5. # Note: Sample
   
6. #-----------------------------------------------------
   
7. #---- Import classes
   
8. print 'Import classes...'
   
9. from org.meteoinfo.layout import MapLayout
   
10. from org.meteoinfo.data.mapdata import MapDataManage
    
11. from org.meteoinfo.data.meteodata import MeteoDataInfo
    
12. from org.meteoinfo.data.meteodata import DrawMeteoData
    
13. from org.meteoinfo.legend import LegendManage
    
14. from org.meteoinfo.legend import LegendType
    
15. from org.meteoinfo.legend import GridLabelPosition
    
16. from org.meteoinfo.layout import LegendStyles
    
17. from org.meteoinfo.shape import ShapeTypes
    
18. from org.meteoinfo.projection import ProjectionInfo
    
19. import os.path
    
20. from java.awt import Color
    
21. from javax.swing import JFrame
    
22. 
    
23. #---- Set directories
    
24. print 'Set directories...'
    
25. baseDir = 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo'
    
26. dataDir = os.path.join(baseDir, 'sample/GrADS')
    
27. mapDir = os.path.join(baseDir, 'map')
    
28. figDir = 'D:/Temp/test'
    
29. 
    
30. #---- Create MapLayout object
    
31. mapLayout = MapLayout()
    
32. mapFrame = mapLayout.getActiveMapFrame()
    
33. 
    
34. #---- Load country layer
    
35. print 'Load country layer...'
    
36. countryLayer = MapDataManage.loadLayer(os.path.join(mapDir, 'country1.shp'))
    
37. lb = countryLayer.getLegendScheme().getLegendBreaks().get(0)
    
38. lb.setDrawFill(False)
    
39. lb.setOutlineColor(Color.blue)
    
40. mapFrame.addLayer(countryLayer)
    
41. 
    
42. #---- Open GrADS data
    
43. print 'Open GrADS data...'
    
44. mdi = MeteoDataInfo()
    
45. mdi.openGrADSData(os.path.join(dataDir, 'model.ctl'))
    
46. 
    
47. #---- Set time index
    
48. mdi.setTimeIndex(2)
    
49. 
    
50. #---- Get pressure grid data
    
51. gdata = mdi.getGridData('PS')
    
52. gdata.extendToGlobal()
    
53. 
    
54. #---- Create legend scheme
    
55. ls = LegendManage.createLegendSchemeFromGridData(gdata, LegendType.GraduatedColor, ShapeTypes.Polygon)
    
56. 
    
57. #---- Create pressure shaded layer
    
58. print 'Create pressure shaded layer...'
    
59. pressLayer = DrawMeteoData.createShadedLayer(gdata, ls, 'Pressure', 'PS', True)
    
60. 
    
61. #---- Add layer
    
62. mapFrame.addLayer(pressLayer)
    
63. 
    
64. #--- Move pressure layer to bottom
    
65. mapFrame.moveLayer(pressLayer, 0)
    
66. 
    
67. #---- Add title
    
68. title = mapLayout.addText('MeteoInfo script demo - Projection', 350, 50, 'Arial', 16)
    
69. 
    
70. #---- Project mapview
    
71. print 'Project mapview...'
    
72. projStr = '+proj=moll +lon_0=105'
    
73. projInfo = ProjectionInfo(projStr)
    
74. mapFrame.getMapView().projectLayers(projInfo)
    
75. 
    
76. #---- Set mapframe
    
77. mapFrame.setDrawNeatLine(False)
    
78. mapFrame.setDrawGridLine(True)
    
79. mapFrame.setDrawGridLabel(False)
    
80. mapFrame.setGridLabelPosition(GridLabelPosition.All)
    
81. mapFrame.setGridXDelt(30)
    
82. mapFrame.setGridYDelt(30)
    
83. 
    
84. #---- Add legend
    
85. legend = mapLayout.addLegend(150, 420)
    
86. legend.setLegendStyle(LegendStyles.Bar_Horizontal)
    
87. legend.setLegendLayer(pressLayer)
    
88. 
    
89. #---- Output figure
    
90. print 'Output figure...'
    
91. mapLayout.paintGraphics()
    
92. #mapLayout.exportToPicture(os.path.join(figDir, 'pressure_test.png'))
    
93. 
    
94. frame = JFrame('MeteoInfo Script Sample', size = (750, 530))
    
95. frame.add(mapLayout)
    
96. frame.visible = True
    
97. print 'Finished!'

#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2014-11-22<br /> # Purpose: Read GrADS data and plot<br /> # Note: Sample<br /> #-----------------------------------------------------<br /> #---- Import classes<br /> print 'Import classes...'<br /> from org.meteoinfo.layout import MapLayout<br /> from org.meteoinfo.data.mapdata import MapDataManage<br /> from org.meteoinfo.data.meteodata import MeteoDataInfo<br /> from org.meteoinfo.data.meteodata import DrawMeteoData<br /> from org.meteoinfo.legend import LegendManage<br /> from org.meteoinfo.legend import LegendType<br /> from org.meteoinfo.legend import GridLabelPosition<br /> from org.meteoinfo.layout import LegendStyles<br /> from org.meteoinfo.shape import ShapeTypes<br /> from org.meteoinfo.projection import ProjectionInfo<br /> import os.path<br /> from java.awt import Color<br /> from javax.swing import JFrame<br /> <br /> #---- Set directories<br /> print 'Set directories...'<br /> baseDir = 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo'<br /> dataDir = os.path.join(baseDir, 'sample/GrADS')<br /> mapDir = os.path.join(baseDir, 'map')<br /> figDir = 'D:/Temp/test'<br /> <br /> #---- Create MapLayout object<br /> mapLayout = MapLayout()<br /> mapFrame = mapLayout.getActiveMapFrame()<br /> <br /> #---- Load country layer<br /> print 'Load country layer...'<br /> countryLayer = MapDataManage.loadLayer(os.path.join(mapDir, 'country1.shp'))<br /> lb = countryLayer.getLegendScheme().getLegendBreaks().get(0)<br /> lb.setDrawFill(False)<br /> lb.setOutlineColor(Color.blue)<br /> mapFrame.addLayer(countryLayer)<br /> <br /> #---- Open GrADS data<br /> print 'Open GrADS data...'<br /> mdi = MeteoDataInfo()<br /> mdi.openGrADSData(os.path.join(dataDir, 'model.ctl'))<br /> <br /> #---- Set time index<br /> mdi.setTimeIndex(2)<br /> <br /> #---- Get pressure grid data<br /> gdata = mdi.getGridData('PS')<br /> gdata.extendToGlobal()<br /> <br /> #---- Create legend scheme<br /> ls = LegendManage.createLegendSchemeFromGridData(gdata, LegendType.GraduatedColor, ShapeTypes.Polygon)<br /> <br /> #---- Create pressure shaded layer<br /> print 'Create pressure shaded layer...'<br /> pressLayer = DrawMeteoData.createShadedLayer(gdata, ls, 'Pressure', 'PS', True)<br /> <br /> #---- Add layer<br /> mapFrame.addLayer(pressLayer)<br /> <br /> #--- Move pressure layer to bottom<br /> mapFrame.moveLayer(pressLayer, 0)<br /> <br /> #---- Add title<br /> title = mapLayout.addText('MeteoInfo script demo - Projection', 350, 50, 'Arial', 16)<br /> <br /> #---- Project mapview<br /> print 'Project mapview...'<br /> projStr = '+proj=moll +lon_0=105'<br /> projInfo = ProjectionInfo(projStr)<br /> mapFrame.getMapView().projectLayers(projInfo)<br /> <br /> #---- Set mapframe<br /> mapFrame.setDrawNeatLine(False)<br /> mapFrame.setDrawGridLine(True)<br /> mapFrame.setDrawGridLabel(False)<br /> mapFrame.setGridLabelPosition(GridLabelPosition.All)<br /> mapFrame.setGridXDelt(30)<br /> mapFrame.setGridYDelt(30)<br /> <br /> #---- Add legend<br /> legend = mapLayout.addLegend(150, 420)<br /> legend.setLegendStyle(LegendStyles.Bar_Horizontal)<br /> legend.setLegendLayer(pressLayer)<br /> <br /> #---- Output figure<br /> print 'Output figure...'<br /> mapLayout.paintGraphics()<br /> #mapLayout.exportToPicture(os.path.join(figDir, 'pressure_test.png'))<br /> <br /> frame = JFrame('MeteoInfo Script Sample', size = (750, 530))<br /> frame.add(mapLayout)<br /> frame.visible = True<br /> print 'Finished!'

1. #-----------------------------------------------------
   
2. # Author: Yaqiang Wang
   
3. # Date: 2014-11-22
   
4. # Purpose: Read GrADS data and plot
   
5. # Note: Sample
   
6. #-----------------------------------------------------
   
7. #---- Import classes
   
8. print 'Import classes...'
   
9. from org.meteoinfo.layout import MapLayout
   
10. from org.meteoinfo.data.mapdata import MapDataManage
    
11. from org.meteoinfo.data.meteodata import MeteoDataInfo, DrawMeteoData
    
12. from org.meteoinfo.legend import LegendManage, LegendType, GridLabelPosition
    
13. from org.meteoinfo.layout import LegendStyles
    
14. from org.meteoinfo.shape import ShapeTypes
    
15. from org.meteoinfo.projection import ProjectionInfo
    
16. from org.meteoinfo.global import Extent
    
17. import os.path
    
18. from java.awt import Color
    
19. from javax.swing import JFrame
    
20. 
    
21. #---- Set directories
    
22. print 'Set directories...'
    
23. baseDir = 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo'
    
24. dataDir = os.path.join(baseDir, 'sample/GrADS')
    
25. mapDir = os.path.join(baseDir, 'map')
    
26. figDir = 'D:/Temp/test'
    
27. 
    
28. #---- Create MapLayout object
    
29. mapLayout = MapLayout()
    
30. mapFrame = mapLayout.getActiveMapFrame()
    
31. mapLayout.getActiveLayoutMap().setWidth(650)
    
32. mapLayout.getActiveLayoutMap().setHeight(380)
    
33. 
    
34. #---- Load country layer
    
35. print 'Load country layer...'
    
36. countryLayer = MapDataManage.loadLayer(os.path.join(mapDir, 'country1.shp'))
    
37. lb = countryLayer.getLegendScheme().getLegendBreaks().get(0)
    
38. lb.setDrawFill(False)
    
39. lb.setOutlineColor(Color.blue)
    
40. mapFrame.addLayer(countryLayer)
    
41. 
    
42. #---- Open GrADS data
    
43. print 'Open GrADS data...'
    
44. mdi = MeteoDataInfo()
    
45. mdi.openGrADSData(os.path.join(dataDir, 'model.ctl'))
    
46. 
    
47. #---- Set time index
    
48. mdi.setTimeIndex(2)
    
49. 
    
50. #---- Get pressure grid data
    
51. gdata = mdi.getGridData('PS')
    
52. gdata.extendToGlobal()
    
53. 
    
54. #---- Create legend scheme
    
55. ls = LegendManage.createLegendSchemeFromGridData(gdata, LegendType.GraduatedColor, ShapeTypes.Polygon)
    
56. 
    
57. #---- Create pressure shaded layer
    
58. print 'Create pressure shaded layer...'
    
59. pressLayer = DrawMeteoData.createShadedLayer(gdata, ls, 'Pressure', 'PS', True)
    
60. 
    
61. #---- Add layer
    
62. mapFrame.addLayer(pressLayer)
    
63. 
    
64. #--- Move pressure layer to bottom
    
65. mapFrame.moveLayer(pressLayer, 0)
    
66. 
    
67. #---- Add title
    
68. title = mapLayout.addText('MeteoInfo script demo - Lambert projection', 350, 30, 'Arial', 16)
    
69. 
    
70. #---- Project mapview
    
71. print 'Project mapview...'
    
72. projStr = '+proj=lcc +lat_1=30 +lat_2=60 +lat_0=0 +lon_0=110'
    
73. projInfo = ProjectionInfo(projStr)
    
74. mapFrame.getMapView().projectLayers(projInfo)
    
75. extent = Extent(78, 130, 15, 53)
    
76. mapFrame.getMapView().zoomToExtentLonLatEx(extent)
    
77. 
    
78. #---- Set mapframe
    
79. mapFrame.setDrawNeatLine(True)
    
80. mapFrame.setDrawGridLine(True)
    
81. mapFrame.setDrawGridLabel(True)
    
82. mapFrame.setGridLabelPosition(GridLabelPosition.LeftBottom)
    
83. mapFrame.setGridXDelt(10)
    
84. mapFrame.setGridYDelt(10)
    
85. 
    
86. #---- Add legend
    
87. legend = mapLayout.addLegend(150, 440)
    
88. legend.setLegendStyle(LegendStyles.Bar_Horizontal)
    
89. legend.setLegendLayer(pressLayer)
    
90. 
    
91. #---- Output figure
    
92. print 'Output figure...'
    
93. mapLayout.paintGraphics()
    
94. #mapLayout.exportToPicture(os.path.join(figDir, 'pressure_test.png'))
    
95. 
    
96. frame = JFrame('MeteoInfo Script Sample', size = (750, 530))
    
97. frame.add(mapLayout)
    
98. frame.visible = True
    
99. print 'Finished!'

#-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2014-11-22<br /> # Purpose: Read GrADS data and plot<br /> # Note: Sample<br /> #-----------------------------------------------------<br /> #---- Import classes<br /> print 'Import classes...'<br /> from org.meteoinfo.layout import MapLayout<br /> from org.meteoinfo.data.mapdata import MapDataManage<br /> from org.meteoinfo.data.meteodata import MeteoDataInfo, DrawMeteoData<br /> from org.meteoinfo.legend import LegendManage, LegendType, GridLabelPosition<br /> from org.meteoinfo.layout import LegendStyles<br /> from org.meteoinfo.shape import ShapeTypes<br /> from org.meteoinfo.projection import ProjectionInfo<br /> from org.meteoinfo.global import Extent<br /> import os.path<br /> from java.awt import Color<br /> from javax.swing import JFrame<br /> <br /> #---- Set directories<br /> print 'Set directories...'<br /> baseDir = 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo'<br /> dataDir = os.path.join(baseDir, 'sample/GrADS')<br /> mapDir = os.path.join(baseDir, 'map')<br /> figDir = 'D:/Temp/test'<br /> <br /> #---- Create MapLayout object<br /> mapLayout = MapLayout()<br /> mapFrame = mapLayout.getActiveMapFrame()<br /> mapLayout.getActiveLayoutMap().setWidth(650)<br /> mapLayout.getActiveLayoutMap().setHeight(380)<br /> <br /> #---- Load country layer<br /> print 'Load country layer...'<br /> countryLayer = MapDataManage.loadLayer(os.path.join(mapDir, 'country1.shp'))<br /> lb = countryLayer.getLegendScheme().getLegendBreaks().get(0)<br /> lb.setDrawFill(False)<br /> lb.setOutlineColor(Color.blue)<br /> mapFrame.addLayer(countryLayer)<br /> <br /> #---- Open GrADS data<br /> print 'Open GrADS data...'<br /> mdi = MeteoDataInfo()<br /> mdi.openGrADSData(os.path.join(dataDir, 'model.ctl'))<br /> <br /> #---- Set time index<br /> mdi.setTimeIndex(2)<br /> <br /> #---- Get pressure grid data<br /> gdata = mdi.getGridData('PS')<br /> gdata.extendToGlobal()<br /> <br /> #---- Create legend scheme<br /> ls = LegendManage.createLegendSchemeFromGridData(gdata, LegendType.GraduatedColor, ShapeTypes.Polygon)<br /> <br /> #---- Create pressure shaded layer<br /> print 'Create pressure shaded layer...'<br /> pressLayer = DrawMeteoData.createShadedLayer(gdata, ls, 'Pressure', 'PS', True)<br /> <br /> #---- Add layer<br /> mapFrame.addLayer(pressLayer)<br /> <br /> #--- Move pressure layer to bottom<br /> mapFrame.moveLayer(pressLayer, 0)<br /> <br /> #---- Add title<br /> title = mapLayout.addText('MeteoInfo script demo - Lambert projection', 350, 30, 'Arial', 16)<br /> <br /> #---- Project mapview<br /> print 'Project mapview...'<br /> projStr = '+proj=lcc +lat_1=30 +lat_2=60 +lat_0=0 +lon_0=110'<br /> projInfo = ProjectionInfo(projStr)<br /> mapFrame.getMapView().projectLayers(projInfo)<br /> extent = Extent(78, 130, 15, 53)<br /> mapFrame.getMapView().zoomToExtentLonLatEx(extent)<br /> <br /> #---- Set mapframe<br /> mapFrame.setDrawNeatLine(True)<br /> mapFrame.setDrawGridLine(True)<br /> mapFrame.setDrawGridLabel(True)<br /> mapFrame.setGridLabelPosition(GridLabelPosition.LeftBottom)<br /> mapFrame.setGridXDelt(10)<br /> mapFrame.setGridYDelt(10)<br /> <br /> #---- Add legend<br /> legend = mapLayout.addLegend(150, 440)<br /> legend.setLegendStyle(LegendStyles.Bar_Horizontal)<br /> legend.setLegendLayer(pressLayer)<br /> <br /> #---- Output figure<br /> print 'Output figure...'<br /> mapLayout.paintGraphics()<br /> #mapLayout.exportToPicture(os.path.join(figDir, 'pressure_test.png'))<br /> <br /> frame = JFrame('MeteoInfo Script Sample', size = (750, 530))<br /> frame.add(mapLayout)<br /> frame.visible = True<br /> print 'Finished!'

小船

正在学习中，谢谢楼主！