北大《空间探测信息处理技术(IDL)》第10章 地图可视化(54P).ppt

1、空间探测信息处理技术,Interactive Data Language (IDL) 第10章：地图可视化 韩培友，IDL可视化分析与应用 ，西北工业大学出版社，2006年. （第12章：地图可视化）,地图制作原理,把不可展的球面投影到可展曲面 可展曲面有圆锥面、圆柱面、平面 设想的地球是透明体，在投影中心（球心）有一点光源S，向四周辐射投影射线，通过球表面（各点A、B、C、D）射到可展面（投影面）上，得到投影点a、b、c，然后再将投影面展开铺平，又作中心投影将其比例尺缩小，从而制成地图。 地图投影基本类型 从投影面类型划分为：圆锥、圆柱、平面（方位投影） 从投影面与地球位置关系划分为：正轴、

2、横轴、斜轴， 或切或割。,正轴切圆锥投影 正轴割圆锥投影,横轴切圆锥投影 横轴割圆锥投影,横轴切圆柱投影 横方位投影,正轴割圆柱投影 斜轴切圆柱投影,斜轴切圆锥投影 正轴切圆柱投影,正方位投影 斜方位投影,地图制作原理,地图制作原理,A flat map may cause distortion in one or more of these conditions. Conformal maintains shape Equal-area maintains area Equidistant maintains distance Directional maintains some direc

3、tions,地图投影,使用MAP_SET建立地图的各种投影 MAP_SET , P0lat, P0lon, Rot, 说明：为地图显示准备各种投影方式 参数： P0lat : 地图投影中心的纬度，-90 90，默认0； P0lon : 地图投影中心的经度，-180180，默认0； Rot : 地图坐标的旋转角度，顺时针为正，默认0； IDLMAP_SET, /CONTINENTS,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Projection Types ,NAME=string : 用名字设置投影方法，或者用下列之一设置： , /AITOFF

4、| , /ALBERS | , /AZIMUTHAL | , /CONIC | , /CYLINDRICAL | , /GNOMIC | , /GOODESHOMOLOSINE | , /HAMMER | , /LAMBERT | , /MERCATOR | , /MILLER_CYLINDRICAL | , /MOLLWEIDE | , /ORTHOGRAPHIC | , /ROBINSON | , /SATELLITE | , /SINUSOIDAL | , /STEREOGRAPHIC | , /TRANSVERSE_MERCATOR ,地图投影,MAP_SET , Lati, Long

5、i, Rotate, Keywords-Projection Types /CYLINDRICAL : cylindrical equidistant projection 最简单的柱面投影，经纬度等间隔。 MAP_SET, /CYLINDRICAL, 0, 0, /GRID, /CONTINENTS, $ TITLE = Simple Cylindrical Equidistant,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Projection Types /MILLER_CYLINDRICAL : Miller Cylindrical pro

6、jection Miller投影的柱面改进，尽可能保持区域形状不变，非等面积、非等距离、非等角度。 MAP_SET, /MILLER_CYLINDRICAL, 0,0,/ISOTROPIC,/GRID, $ LATDEL=15,LONDEL=20, /CONTINENTS, $ TITLE=Miller Cylindrical,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Projection Types /MERCATOR : Mercator projection柱面投影的一种，尽量保持局部形状，不可显示极区80度以上范围。 MAP_SET, /

7、MERCATOR, 0, 0, /ISOTROPIC, $ /GRID, /CONTINENTS, TITLE = Simple Mercator,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Projection Types /TRANSVERSE_MERCATOR : Transverse Mercator projection，又名Gauss-Krueger projection或UTM MERCATOR旋转90投影，即：以纬度线代替子午线柱面投影。 MAP_SET,/TRANSVERSE_MERCATOR,0,-80,/ISOTROPIC,

8、$ ,/GRID, /CONTINENTS, LATDEL=15, LONDEL=15, $ LIMIT=-80, -140, 80, -20, TITLE=UTM Projection,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Projection Types /MERCATOR : Mercator projection如果偏心率不很重要，配合 CENTRAL_AZIMUTH=90, 以及rotate取90度，可以实现Transverse Mercator投影 MAP_SET, /MERCATOR, 0, 110, 90, CENTRAL_A

9、ZIMUTH=90, $ /ISOTROPIC, LIMIT= -5, 70, 70,140, $ /GRID, LATDEL=15, LONDEL=15, /CONTINENTS,$ TITLE = Transverse Mercator,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Projection Types /GOODESHOMOLOSINE : Goodes Homolosine Projection MAP_SET, /GOODESHOMOLOSINE, 0, 0, /ISOTROPIC, /GRID, $ LATDEL=15,LON

10、DEL=20,/HORIZON, E_CONTINENTS=FILL:1, $ TITLE=Goode Homolosine,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Projection Types /ROBINSON : Robinson pseudo-cylindrical projection 投影中心纬度一般设置为0, Arthur Robinson1963提出，是当今最常用的世界地图投影方式。 MAP_SET, /ROBINSON, 0, 150, /ISOTROPIC, /GRID, $ /HORIZON, E_CONTINENTS

11、=FILL:1, TITLE=Robinson,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Projection Types /AITOFF: Aitoff projection 等方位角度投影，距离测量不准确。 MAP_SET, 0, 180, /AITOFF, /GRID, /CONTINENTS, $ /ISOTROPIC, /HORIZON, TITLE= Aitoff Projection,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Projection Types /HAMMER : Ham

12、mer-Aitoff equal area projection MAP_SET,0,180,/HAMMER,/GRID,/CONTINENTS, /ISOTROPIC, $ /HORIZON, TITLE= Hammer-Aitoff Projection,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Projection Types /MOLLWEIDE : Mollweide projection MAP_SET, /MOLLWEIDE, 0, 0, /ISOTROPIC, /GRID, $ LATDEL=20,LONDEL=20,/HORIZ

13、ON,E_CONTINENTS=FILL:1, $ TITLE=Oblique Mollweide,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Projection Types /SINUSOIDAL : sinusoidal projection 投影点纬度（参数1）和旋转角度（参数3）必须为0 MAP_SET, 0, 150, 0, /SINUSOIDAL, /ISOTROPIC, $ /CONTINENTS, LONDEL=20, $ TITLE=Sinusoidal, /HORIZON,地图投影,MAP_SET , Lati, Longi,

14、 Rotate, Keywords-Projection Types /AZIMUTHAL : azimuthal equidistant projection MAP_SET, /AZIMUTHAL, /ISOTROPIC, -90, $ /GRID, LONDEL=20, LATDEL=20, /CONTINENT, $ /HORIZON, TITLE = Polar Azimuthal,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Projection Types /STEREOGRAPHIC : stereographic projectio

15、n 真透视投影，多用于表现极区。 MAP_SET, /STEREOGRAPHIC, 0, 70, /ISOTROPIC, /GRID, $ LATDEL = 20, LONDEL = 20, /HORIZON, /CONTINENT, $ TITLE = Equatorial Stereographic,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Projection Types /CONIC : Lamberts conformal conic projection 可以同时用STANDARD_PARALLELS= 设置12条平行线 MAP_SE

16、T, /CONIC,40,-80,STANDARD_PARALLELS=20,60, $ /ISOTROPIC, LIMIT=0, -260, 90, 100, $ /GRID, LATDEL=15, LONDEL=20, /CONTINENT, $ TITLE= Lamberts Conic,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Projection Types /LAMBERT : Lamberts azimuthal equal area projection MAP_SET, /LAMBERT, 90, 0, -105, /ISOTR

17、OPIC, $ /GRID, LATDEL=20, LONDEL=20, $ /CONTINENTS, E_CONTINENTS=FILL:1, /HORIZON, $ TITLE = Polar Lambert,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Projection Types /GNOMIC : gnomonic projection（球心投影） MAP_SET, /GNOMIC, 35, 110, LIMIT = 10,80,60,140, $ /ISOTROPIC, /GRID, /CONTINENT, $ TITLE = Obl

18、ique Gnomonic,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Projection Types /ORTHOGRAPHIC : orthographic projection 从无限远处透视 MAP_SET, /ORTHOGRAPHIC,20,110,SCALE=50e6,$ /CONTINENTS, /GRID, LONDEL=10, LATDEL=10, $ TITLE = Oblique Orthographic,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Projection Type

19、s /SATELLITE : 卫星投影 Longi, Lati: 星下点位置 SAT_P=P, Omega, Gamma相机设置 P : 观察点的地心距离(单位：Re)； Omega : 相机的倾角； Gamma : 投影面转角。 MAP_SET, /SATELLITE,SAT_P=1.0251,55,150,41.5,-74., $ /ISOTROPIC, /HORIZON, $ LIMIT=39, -74, 33, -80, 40, -77, 41,-74, $ /CONTINENTS, TITLE=Satellite / Tilted Perspective ; Set up the s

20、atellite projection: MAP_GRID,/LABEL,LATLAB=-75,LONLAB=39,LATDEL=1,LONDEL=1 ; Get North vector: p = convert_coord(-74.5, 40.2, 40.5, /TO_NORM) ; Draw North arrow: ARROW, p(0,0), p(1,0), p(0,1), p(1,1), /NORMAL XYOUTS, -74.5, 40.1, North, ALIGNMENT=0.5,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Pro

21、jection Parameters ,CENTRAL_AZIMUTH=degrees_east_of_north ; 设置投影的方位角中心，影响：Cylindrical、Mercator、Miller、Mollweide和Sinusoidal投影方式；默认为0。 , ELLIPSOID=array ; 设置定义椭球的三个元素a, e2, k0，针对Transverse Mercator和Lambert Conic两种投影； a - 迟到面半径（0）; e 椭球的偏心率； k0 子午线中央尺度； 默认为Clarke1866椭球， 即：6378206.4, 0.00676866, 0.9996,

22、地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Projection Parameters ,/ISOTROPIC; X，Y方向具有相同尺度； , LIMIT=vector; 显示地图的区域限制； 1) 四点定义：Latmin, Lonmin, Latmax, Lonmax 2) 八点定义：Lat0, Lon0, Lat1, Lon1, Lat2, Lon2, Lat3, Lon3; 分别对应于下边、左边、上边和右边的点； ,SAT_P=vector ; SATELLITE投影时使用。 ,SCALE=value ; 地图坐标比例，1:scale ,ST

23、ANDARD_PARALLELS=array ; 1至2元素数组，对于圆锥投影，设置1到2个相切（交）的平行线,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Map Characteristics , /ADVANCE ; 在多图绘制中（!P.MULTI）自动移动到下图处理； ,COLOR=index ; 地图边框颜色； , /CONTINENTS , CON_COLOR=index , /HIRES ; 绘制陆地边界，/HIRES表示精确绘制，但不绘制岛屿和湖泊。如果需要，可以使用如下命令详细绘制： MAP_CONTINENTS, /COASTS

24、CON_COLOR=index; 陆地边界颜色； ,MLINESTYLE=0 | 1 | 2 | 3 | 4 | 5 ; 陆地边界线的风格，默认为0 , MLINETHICK=value ; 陆地边界线的粗细，默认为2 , E_CONTINENTS=structure ; 传递更多信息给MAP_CONTINENTS，例如：要填充陆地部分，则： E_CONTINENTS=FILL:1,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Map Characteristics , E_GRID=structure ;传递更多信息给MAP_GRID, 例如：E_

25、GRID=LABEL:2相当于 MAP_GRID, LABEL=2 , /GRID; 绘制坐标（纬度线或子午线） , GLINESTYLE=0|1|2|3|4|5; 坐标的线型 , GLINETHICK=value; 坐标线粗细 , /HORIZON; 同时绘制地平线 , E_HORIZON=structure ; 传递更多信息给/HORIZON，例如：E_HORIZON=FILL:1,COLOR:3,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Map Characteristics ,LATDEL=value ; 纬度坐标线间隔； ,LONDEL

26、=value ; 经度坐标线的间隔； ,LABEL=nlabel every nth gridline ; 坐标线标注 ,LATALIGN=value0.0 to 1.0 ; 纬度线标注的对齐方式，左对齐、居中、右对齐 ,LATLAB=longitude ; 放置纬度标注的经度位置 LONALIGN=value0.0 to 1.0 ; 经度线标注的对齐方式，左对齐、居中、右对齐 ,LONLAB=latitude ; 放置经度标注的纬度位置,地图投影,MAP_SET , Lati, Longi, Rotate, Keywords-Map Characteristics ,/NOBORDER; 无

27、边框 ,/NOERASE; 非擦除，如果设置/ADVANCE将被忽略 ,REVERSE=0 | 1 | 2 | 3 ; =0: 无反转；=1: X反向；=2: Y反向；X,Y均方向。 ,TITLE=string; 地图标题 ,CHARSIZE=value; 字符大小 ,/CLIP; 非裁减 ,/USA ; 绘制美国各州边界 ,XMARGIN=value , YMARGIN=value ; X,Y边框留白 ,POSITION=x1,y1,z2,y2 ,/T3D ,ZVALUE=value 与PLOT同。,获取地图投影信息,MAP_PROJ_INFO 获取当前地图投影或已有地图投影的信息 MAP_

28、PROJ_INFO , iproj ; iproj: 地图投影序号，如果设置/CURRENT，则获取当前投影序号存入iproj ,AZIMUTHAL=variable; 返回1:等方位角投影，0:其它 ,CIRCLE=variable; 返回1: 圆投影；0: 椭圆投影 ,CYLINDRICAL=variable; 返回1:圆柱投影；0:其它 , /CURRENT; 获取当前的投影序号 , LL_LIMITS=variable ; 返回当前地图的矩形边界 , NAME=variable ; 返回投影名称 ,PROJ_NAMES=variable; 返回所有投影名称，以序号排序 ,UV_LIMI

29、TS=variable ; UV平面的矩形框大小 , UV_RANGE=variable ; UV平面的地图范围,图形的地图显示,绘制图形 在使用MAP_SET以后，保存投影参数在!MAP中， 所有的绘图程序 (PLOT, OPLOT, PLOTS, XYOUTS, and CONTOUR),根据!MAP参数进行坐标变换。 X对应于经度，Y对应于纬度。 例example10_,图形的地图显示,例example10_CityM 打开窗口 获得城市坐标数据和名称数据 用MAP_SET建立投影及绘制陆地边界 绘制城市坐标 绘制城市名称,地图绘制,地图绘制

30、 MAP_CONTINENTS , 在MAP_SET中使用以下关键字则自动调用MAP_CONTINENTS /CONTINENTS 或 E_CONTINENTS=structure ,/COASTS ; 绘制海岸线，同时绘制岛屿和湖泊； ,/CONTINENTS ; 仅绘制陆地边界（默认） ,/COUNTRIES ; 绘制政治边界，以1993年为准。不包括海岸线。所以， MAP_CONTINENTS,/COUNTRIES,/COASTS ,/RIVERS ; 绘制河流 MAP_CONTINENTS,/COUNTRIES,/COASTS,/RIVERS ,/HIRES ; 高分辨地图。 ,LIM

31、IT=Latmin, Lonmin, Latmax, Lonmax ; 区域限制,地图绘制,地图绘制 MAP_CONTINENTS , ,FILL_CONTINENTS=1|2 ; 填充大陆区域，1: 填实；2: 线条。 ,SPACING=centimeters ; 线条填充时的间隔； ,ORIENTATION=value ;线条填充时的走向 , MAP_STRUCTURE=structure ; 设置为由MAP_PROJ_INIT返回的参数，结构同!MAP， ,COLOR=index ; ,MLINESTYLE=0 | 1 | 2 | 3 | 4 | 5 ; 边界线型 ,MLINETHICK

32、=value ; 边界线粗细 , /USA; 美国州界,地图绘制,IDLMAP_SET,/GNOMIC,35,110,LIMIT=15,80,55,140, $ IDL/ISOTROPIC, /GRID, TITLE = Oblique Gnomonic IDLMAP_CONTINENTS,地图绘制,IDL MAP_SET,/GNOMIC,35,110,LIMIT=15,80,55,140, $ IDL /ISOTROPIC, /GRID, TITLE = Oblique Gnomonic IDL MAP_CONTINENTS,/countries,地图绘制,IDL MAP_SET,/GNOM

33、IC,35,110,LIMIT=15,80,55,140, $ IDL /ISOTROPIC, /GRID, TITLE = Oblique Gnomonic IDL MAP_CONTINENTS,/rivers,地图绘制,device,decomposed=0 loadct,12 MAP_SET, /GNOMIC, 35, 110, LIMIT = 15,80,55,140, $ /ISOTROPIC, /GRID, TITLE = Oblique Gnomonic MAP_CONTINENTS,/countries, MLINESTYLE=3,MLINETHICK=1.5 MAP_CONT

34、INENTS,/coasts MAP_CONTINENTS,/rivers,color=20,地图经纬度网格,MAP_GRID 在MAP_SET中使用以下关键字则自动调用MAP_GRID /GRID 或 E_GRID=structure， MAP_SET的有关关键字还有： , GLINESTYLE=0|1|2|3|4|5; 坐标的线型 , GLINETHICK=value; 坐标线粗细 MAP_GRID , /BOX_AXES ; 边框坐标，两边标纬度线，上下标经度线； /CLIP_TEXT ; 裁减标注字，如果/BOX_AXES，则被忽略； ,GLINESTYLE=0|1|2|3|4|5;

35、坐标线型； ,GLINETHICK=value; 坐标线粗细； ,/HORIZON; 画出当前的地平线（外框） ,/FILL_HORIZON; 填满地平线 ,INCREMENT=value; 坐表线的绘制精细程度 ,LABEL=nlabel_every_nth_gridline ; 每 n 根坐标线标注一次标签,地图经纬度网格,MAP_GRID MAP_GRID , ,LATDEL=degrees ; 纬度坐标线的间隔，单位：度； ,LATS=vector, LATNAMES=array ; 纬度坐标线及标签名字，用矢量和数组方式定义； ,LONDEL=degrees ; 经度坐标线的间隔，单

36、位：度； ,LONS=vector ,LONNAMES=array ; 经度坐标线及标签名字，用矢量和数组方式定义； ,MAP_STRUCTURE=structure ; 可以不使用!MAP，而直接使用MAP_PROJ_INIT产生的结构； ,/NO_GRID ; 只给出标签而不画出网格线； COLOR; 网格线颜色； CHARSIZE; 标签字大小；,图形的地图显示,获取经纬度坐标 CURSOR, X, Y , Wait | , /CHANGE | , /DOWN | , /NOWAIT | , /UP | , /WAIT , /DATA | , /DEVICE, | , /NORMAL 说

37、明：此时X返回经度，Y返回纬度 例example10_GetC,地图投影计算,MAP_PROJ_INIT 在不影响!MAP的情况下进行地图投影计算 可以使用IDL自己的程序，或 U.S. Geological Surveys General Cartographic Transformation Package (GCTP 2.0) Result = MAP_PROJ_INIT( Projection ,DATUM=value ,/GCTP ,LIMIT=vector ,/RADIANS ,/RELAXED ) 参数 Projection: 投影方法，可以使序号，也可以是名称

38、 Result: 与!MAP相同的结构，供 MAP_PROJ_FORWARD 和 MAP_PROJ_INVERSE 使用。 参见：IDL Reference Guide,地图投影计算,MAP_PROJ_FORWARD 根据地图投影把经纬度坐标（Long, Lati）转换为直角坐标(X,Y). 投影参数可以使用!MAP，也可以由MAP_PROJ_INIT计算的结果。 Result = MAP_PROJ_FORWARD(Longitude, Latitude ; Result : (2,n) 代表 (X,Y) ; Longitude, Latitude: 输入的经纬度 , CONNECTIVITY

39、=vector; 输入连通性参数序号 , MAP_STRUCTURE=value; 地图投影参数结构变量 , POLYGONS=variable; 输出多边形绘制的顶点序号 , POLYLINES=variable; 输出多线绘制的顶点序号 , /FILL; 对多边形在裁减或分裂时进行镶嵌处理 , /RADIANS ); 输入为弧度,地图投影计算,MAP_PROJ_INVERSE 是MAP_PROJ_FORWARD的反运算，直角坐标(X, Y)转换为经纬度坐标（Long, Lati）。 投影参数可以使用!MAP，也可以由MAP_PROJ_INIT计算的结果。 Result = MAP_PROJ

40、_INVERSE (X , Y ; Result : (2, n)数组表达经纬度 , MAP_STRUCTURE=value; 地图投影参数结构变量 , /RADIANS ); 输入为弧度,地图投影计算,例： example10_ex_map_proj_ ; Helper function. Constructs the polyline objects. PRO Ex_Map_AddPolyline, label, gridLon, gridLat, $ sMap, oModel, oContainer, oFont, $ LONGITUDE = longitude .

41、 END ; Main function. Creates a grid over a map ; projection. PRO exam09_Ex_Map_Proj_Forward . END,地图投影计算,MAP_2POINTS函数 计算两个地理坐标点之间的距离、方位角、大圆路径、航线（罗盘线）路径等参数。 Result = MAP_2POINTS(lon0, lat0, lon1, lat1 ; Result : 一般为两元素数组，包含：两点间的距离和方位角（或航线角度）。 ; lon0, lat0, lon1, lat1 : 两点的经纬度 ,DPATH=value ; 两点间的最大角

42、距离 | , /METERS | , /MILES |; 距离单位，默认为角度距离 , NPATH=integer2 or greater ; 使 Result 返回由NPATH定义点数的大圆路径 ; 对于大圆：路径按照等距离间隔； ; 对于航线：路径按照等经度间隔； | , /PARAMETERS ; 使大圆路径返回参数方程sin(c),cos(c),sin(az),cos(az) ; c 为大圆角距离，az为在P0点的方位角 | , RADIANS=value ; 角度单位为弧度，否则为度 , /RADIUS; 计算距离用的地球半径 , /RHUMB ) ; 设置为航线输出,地图投影计算,

43、MAP_2POINTS函数 例：example10_ PRO example10_map2points ; The two points Longitude, latitude in degrees. B = -105.19, 40.02 L = -0.07, 51.30 ; Print the angular distance and azimuth, from B, ; of the great circle connecting the two points: X = MAP_2POINTS(B0, B1, L0, L1) PRINT, X, format=

44、$ (Distance=, F8.4, Degree; Azimuth=, F8.4, Degree) ; The angular distance and course, connecting the two points: X = MAP_2POINTS(B0, B1, L0, L1,/RHUMB) PRINT, X, format= $ (Distance=, F8.4, Degree; course=, F8.4, Degree) DEVICE, DECOMPOSED=1 !P.BACKGROUND=FFFFFFXL !P.COLOR=000000XL MAP_SET, /MOLLWE

45、IDE, 40,-50, /GRID, SCALE=75e6,/CONTINENTS PLOTS, MAP_2POINTS(B0, B1, L0, L1,/RHUMB, DPATH=1) PLOTS, MAP_2POINTS(B0, B1, L0, L1,DPATH=1),linestyle=3 END,地图投影计算,MAP_PATCH函数 把一个图像弯曲到当前地图投影。方法是把图像分解为小三角形。 Result = MAP_PATCH( Image_Orig ; Result : 返回弯曲到地图上的图像或数据集合 ; Image_Orig : 12维数组，行列分别为递增的经度或纬度； , L

46、ons, Lats ; 经度坐标和纬度坐标，可以是向量或2维数组，但需要与Image_Orig的尺寸相配合； , LAT0=value , LAT1=value ; 第一行和最后一行数据的纬度，缺省为-90度和90度 , LON0=value , LON1=value ; 第一列和最后一列数据的纬度，缺省为-180度和180-360/行数 , MAX_VALUE=value; 大于此值的图像数据将设为MISSING , MISSING=value; 定义MISSING数据 , /TRIANGULATE; 三角化 , XSIZE=variable , XSTART=variable ; 返回：X大小和X起点 , YSIZE=variable , YSTART=variable ) ; 返回：Y大小和Y起点,地图投影计算,MAP_PATCH函数 例：example10_ PRO example10_mappatch ; Form a 24 x 24 data