파이톤의 세계지도 -20(지도 이미지에 위치 정보 부여)

mpl_toolkits.베이스맵의 샘플 프로그램을 찾다가 해외 사이트에서 발견됐습니다.tif를 불러오는 프로그램이 자주 있습니다.
프로그램 자체가 실행되고 있지만 tif를 잘 몰라서 샘플을 열고 싶어도 tif를 열 수 없어요??
나는 이전에 tif 파일을 열어본 적이 있다??
(tif가 열리면 tiff가 열리지 않습니다.)
조사하면 GEOTIFF라는 위치 정보를 첨부한 tif 파일인 것 같습니다.
QGIS라는 애플리케이션으로 제작이 가능하기 때문에 설치했습니다.
설치하다.
Ubuntu Studio 18.04시
qgis/bionic 2.18.17 + dfsg-1 am d64 지리정보시스템(GIS)

$ sudo apt install qgis

최신 버전을 설치할 때
Windows·Mac·Linux(사용자 환경에 맞는 QGIS 다운로드)
최신 버전 QGIS 3.2.0'Bonn'22.062018
https://qgis.org/ja/site/forusers/download.html
맵 이미지 만들기
정부기관 홈페이지에서 다운로드하거나 구글맵 등으로 제작한다.
mpl입니다.toolkits.베이스맵으로 제작되었습니다.(이것도 요점을 결정할 때 초점을 맞추면 위도, 경도선이 굵어 요점을 정확하게 결정하기 어렵다.)
한때 인터넷상에서 정밀도가 떨어지는 일본 전역의 지도에서 한 번 해봤지만 실제 위치와 100km 단위로 떨어져 쓸모 있는 것이 되지 못했다.

#!/usr/bin/python3
# coding: UTF-8

import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap

fig = plt.figure(figsize=(10,8))
map = Basemap(projection='cyl',
                 resolution='f',
                 llcrnrlon=123,
                 llcrnrlat=24,
                 urcrnrlon=147,
                 urcrnrlat=46)

map.drawcoastlines(color='lightgray')
map.drawcountries(color='lightgray')
map.fillcontinents(color='white', lake_color='#eeeeee');
map.drawmapboundary(fill_color='#eeeeee')
map.drawparallels(np.arange(24.0, 46, 1.0), labels = [1,0,0,0], fontsize=12)
map.drawmeridians(np.arange(123.0, 147, 2.0), labels = [0,0,0,1], fontsize=12)

plt.show()

지도 이미지에 위치 정보를 부여합니다.
참고 자료
투시 공급선
https://docs.qgis.org/2.18/ja/docs/user_manual/plugins/plugins_georeferencer.html
1. 지도 이미지를 만든다.
구글맵 등 지도를 이용할 때는'이곳에 대해'등을 우클릭해 동경, 북위를 기호로 하는 도로와 도로의 교차점 등을 먼저 조사한다.(4곳)
mpl_toolkits.베이스맵으로 제작할 때 동경, 북위선을 미리 표시한다.
지도상의 도안이 모두 좋은 것은 아니며, 어떤 것들은 GIS에 적합하지 않다고 한다.
Google Map 등은 웹 전자계산기의 도법인 것 같아서 웹 전자계산기의 도법이 없기 때문에 이번에는 전자계산기의 도법 프로젝트='merc로 제작되었습니다.(projection='merc'로 제작하면 오차가 엄청납니다.)

주의

프로젝트='cyl'로 만들어야 합니다.
'시엘'로 만들면 오차가 거의 없어요.
2. QGIS Desktop을 시작합니다.
지도를 보존하다.
프로젝트 (J) 카페를 (A) 카페로 저장하고 지정한 폴더의 이름을 입력하십시오.
****.qgs
3. 위치 정보를 추가하기 위해 투시경을 사용해야 한다.
플러그인(P) - 플러그인 관리 및 설치...카페 [GDAL 시스루 안경]을 선택하고 닫습니다.
4. 래스터 (R) > 투시 디스크 (G) > 투시 디스크 (G) > 투시 디스크 창이 표시됩니다.
파일 - 래스터를 엽니다...카페에서 위치 정보를 부여한 이미지를 엽니다.
공간 참조 시스템 선택 대화 상자를 엽니다.
WGS84를 선택하고 OK 를 클릭하십시오.모니터 지도.
5. 지도에 위치 정보를 추가한다.
'투시경 편집 창' > 포인트 추가 > 을 누르면 지도 좌표 입력 대화상자를 표시하고, 지점의 동경과 북위를 누르면 OK를 누르십시오.카페 > GCP 테이블에 변환 소스와 변환 목적지의 좌표를 표시합니다.
주의)
이 작업이 제대로 진행되지 않으면 지도가 엇갈릴 것이다.
특히 상술한 지도처럼 범위가 넓은 지도는 특히 빗나간다.
정확한 점을 확인하려면 지도를 축소해야 합니다.
이 일을 네 군데 진행하다.오류가 있으면 두 번 클릭하여 수정합니다.


6. 왼쪽 위 모서리의 오른쪽 삼각형에서 원근 시작을 클릭합니다.카페 정보 대화 상자가 표시됩니다.ー＞
OK > 변환에 대한 설정 대화 상자를 표시합니다.ー＞
변환 유형: 다항식 1
재샘플링 방법:가장 가까운 이웃
대상 SRS: 기본 CRS(EPSG: 4326-WGS84)
출력 래스터:폴더 및 파일 이름 지정
OK
왼쪽 위 오른쪽 삼각형의 지구 참조 시작을 다시 클릭하면 GEOTTIFF 파일이 생성됩니다.
계산한 후 지도가 비뚤어지면 그곳에는 붉은 선이 보일 것이다.(포인트 추가 지점)
이 실은 너무 길면 사용할 수 없다.(길수록 오차가 크다)

QGIS Desktop에서 만든 tif 파일 열기
생성된 tif 파일을 열려면 래스터 레이어를 추가합니다.
마우스 커서가 손바닥에 있고 지물 정보를 클릭하면 마우스 커서가 화살표로 바뀐다.
이렇게 하면 위도선과 경도선의 교점에 마우스 커서를 놓으면 창 아래쪽의 좌표가 대체적으로 좋은 값이 된다.
나는 좀 이상하다고 생각한다.
확인
상기 지도상에서 나는 4개국의 음영 파일을 만들어 지도에 그렸지만 포인트를 결정할 때 적당히 잘못한 것 같아서 약간 엇갈렸다.
상기 지도 자체의 범위가 넓기 때문인가, 아니면 Resolution=의 정밀도가 낮고 다각형이 거칠기 때문인가.