Stereographic 그릴 때 발생하는 문제 해결

[[KS서울날씨청] KS Meteorology]

안녕하세요,

현재 Stereographic 투영법을 통해 지도 위에 contour를 실행하고 있습니다.

win= window(dimension= [1000, 1000])

  map= map('Stereographic', $;Lambert Conformal Conic
    limit= [20, -180, 90, 180] ,CLIP=0, $;
    margin= 0.15, /current, aspect_ratio=0, CENTER_LATITUDE= 90, CENTER_LONGITUDE=150);STANDARD_PAR1=30, STANDARD_PAR2=60
  
    c = CONTOUR( regr[*,*],  longitude[*,*] , latitude[*,*], C_THICK=1.3, C_LABEL_SHOW=1, c_use_label_ORIENTATION= 1, $
    C_LABEL_INTERVAL = 0.63, overplot=map, font_size= 10.1, font_style= 1, $
    C_VALUE= [-0.9:0.9:0.2],GRID_UNITS='degrees', /fill, rgb_table=70)
  
   d = CONTOUR( regr[*,*],  longitude[*,*] , latitude[*,*], C_THICK= 1.3, C_LABEL_SHOW=1, c_use_label_ORIENTATION= 1, $
   C_LABEL_INTERVAL = 0.63, overplot=map, font_size= 10.1, font_style= 1, $
   C_VALUE= [-0.9:0.9:0.2],GRID_UNITS='degrees', color= 'black')
 
  statistic_95_plus = CONTOUR(stat_95[*,*],longitude[*,*] , latitude[*,*], C_THICK= 3,c_linestyle=2, $
   /overplot, C_VALUE=[0.9:1.1:0.1] , C_LABEL_SHOW=0, c_use_label_ORIENTATION= 1, $
   C_LABEL_INTERVAL = 0.6, color= 'yellow',GRID_UNITS='degrees');,, shading=1, SHADE_RANGE= [1400,1000]

 statistic_99_plus = CONTOUR(stat_99[*,*],longitude[*,*] , latitude[*,*], C_THICK= 3, $
   /overplot, C_VALUE=[0.9:1.1:0.1] ,  C_LABEL_SHOW=0, c_use_label_ORIENTATION= 1, $
   C_LABEL_INTERVAL = 0.6, color= 'yellow',GRID_UNITS='degrees');, shading=1, SHADE_RANGE=[1400,1000]
 

  e = CONTOUR( regr[*,*],  longitude[*,*] , latitude[*,*], C_THICK=3.3, C_LABEL_SHOW=1, c_use_label_ORIENTATION= 1, $
    C_LABEL_INTERVAL = 0.63, overplot=map, font_size= 10.1, font_style= 1, $
    C_VALUE= [0:0],color='black',GRID_UNITS='degrees')

  mc= mapcontinents(/countries, thick= 1.4, /hires)

  map.mapgrid.label_position=0
  map.mapgrid.grid_longitude=30
  map.mapgrid.GRID_LATITUDE= 5

  map['latitudes'].label_angle= -1
  map['latitudes'].label_show= 0
  map['20N'].label_show= 1

  map['70N'].label_show= 1
  map['longitudes'].label_angle= 0
  map.mapgrid.thick= 1.7
  map.mapgrid.font_size= 15
  map.mapgrid.transparency= 20
  map.horizon_linestyle=0 ;make a box line
  map.HORIZON_THICK=2
  map.HORIZON_color= 'black
  map.linestyle=1
  map.mapgrid.FONT_STYLE= 'bold'

  map['25N'].hide=1
  map['35N'].hide=1
  map['45N'].linestyle= 6
  map['55N'].hide=1
  map['65N'].hide=1
  map['45N'].label_show= 1

  map['120E'].label_angle= 349
  map['150E'].label_angle= 11
 
  col = COLORBAR(TARGET=c, position= [0.1,0.05,0.9,0.08],BORDER=1, TAPER=1)

위 코드에서 statistic_95_plus (statistic_99_plus)는 통계구간 95 & 99%에 해당되는 곳을 노란색 선으로 표시를 한걸 지시하는겁니다.

그리고 regr을 contour 하는건 그림에서 보이는 바입니다.

제가 드리고 싶은 질문,

1. 이상하게 180도 부분에서 검은색 선이 칠해집니다 (map.horizon_linestyle=0, map.horizon_thick=2) . 그 부분 해결 방법...

2. colorbar 색깔을 역으로 돌리는 방법 (파란색 왼쪽, 빨간색 오른쪽, 가운데 하얀색, 현재는 rgb_table=70 사용)

3. 북위 20, 45, 70도 부분을 좌측 상단 (2사분면) 으로 옮기는 방법

4. 경도 숫자들을 가능한 검은선 밖으로 ㅡ 자 모양으로 배열시키는 법

항상 감사드립니다.

[[KS서울날씨청] KS Meteorology]

1. 이상하게 180도 부분에서 검은색 선이 칠해집니다 (map.horizon_linestyle=0, map.horizon_thick=2) . 그 부분 해결 방법...

--> polar stereographic 으로 해결

2. colorbar 색깔을 역으로 돌리는 방법 (파란색 왼쪽, 빨간색 오른쪽, 가운데 하얀색, 현재는 rgb_table=70 사용)

3. 북위 20, 45, 70도 부분을 좌측 상단 (2사분면) 으로 옮기는 방법

--> label_position = 0.1 으로 해결

4. 경도 숫자들을 가능한 검은선 밖으로 ㅡ 자 모양으로 배열시키는 법

2번 colorbar에서 파란색을 왼쪽, 빨간색으로 오른쪽에 두고, -1~1.1 (0.2) 간격 & -0.1~0.1 (하얀색) 으로 지정하는 법을 알고 싶습니다 ㅠㅠ

2020년 10월 19일 월요일 오후 4시 28분 35초 UTC+9에 KS Meteorology, [KS서울날씨청]님이 작성:

[이상우]

컬러테이블을 커스터마이즈하는데 있어서는 기본적으로 COLORTABLE 함수를 이용하는 방법이 있습니다. 그 외에도 제 IDL 블로그에서 '컬러테이블' 또는 'colortable'로 검색되는 게시물들의 내용들도 한번 참조해보시기 바랍니다.

2020년 10월 20일 화요일 오후 12시 19분 15초 UTC+9에 jungji...@gmail.com님이 작성: