답변 감사합니다.

저도 비슷한 방법을 생각하기는 해보았습니다.
다만, 시도 경계를 서로 다른색으로 칠하고 그림을 저장한 후 다시 불러올 때 채워진 위경도 정보를 가져올 수 있는지 궁금합니다.
저는 시도 구분에 대한 위경도를 얻고 또 다른 격자화된 자료에 구분된 시도 자료에 대한 결과를 얻고 싶다는게 목적입니다.
저도 실력이 부족하여 현재는 시도에 대한 경계 부분의 위경도 정보만 획득할 수 있는 정도입니다..

좋은 의견이 있으셔서 알려주시면 감사하겠습니다.

위경도는 2차원으로 만들어야 하는데
irlat = findgen(401)*0.005+36.5
로 하면 1차원이 되는 문제가 있습니다.

또한 0.005로 하여도 , 가장자리가 다르게 표현됩니다.

limit=[lonmin2,latmin2,lonmax2,latmax2]
이런식으로 되어있습니다. 이정도 인데도 더 좁혀야 하나요??ㅜㅜ
테두리에 원치 않는 결과 때문에 sph_scat함수를 못쓰는거 아닌가 합니다..ㅠㅠ

또한 제가 map_image로 표현을 하는데,
irlat와 irlon을 구하지않으면 어떻게 위경도 값을 입력하죠??