세계지도-3(일본지도·cartopy판)

6023 단어 UbuntuPython
세계 지도
https://qiita.com/ty21ky/items/c7a93b83cbbfa8a82a17
세계지도-2
https://qiita.com/ty21ky/items/c9e29a5361e35ca7c4f2
계속하다
일본 지도에 현경을 긋다
지도에서basemap과cartopy로
http://www.mwsoft.jp/programming/numpy/plot_map.html
상술한 사이트에는 일본 지도에 현경을 표시하는 방법이 쓰여 있다.상세한 상황은 상술한 사이트를 참조하십시오.
설치
$ git clone https://github.com/OSGeo/proj.4.git
$ cd proj.4
$ git checkout -b 4.9.0 tags/4.9.0
$ ./configure
$ make
$ sudo make install

$ sudo apt-get install libgeos-dev

$ git clone https://github.com/SciTools/cartopy.git
$ cd cartopy
$ git tag

$ git checkout -b v0.14.2 tags/v0.14.2
여기까지 잘 되고 있지만 다음
$ sudo pip install .
이 색상은 색상이 바래집니다.정보를 번역하면 sudo에 -H라고 쓰여 있어서 넣으면 OK.
$ sudo -H pip install .
$ python -c "import cartopy"
여기 또 틀렸어.
libproj.so.0이 없는 것 같습니다.
$ ls -l/usr/local/lib/libp*
-rw-r--r-1root root37918 5월 2 14:21/usr/local/lib/libproj.a
- rwxr-xr-x1root root953 5월 2 14:21/usr/local/lib/libproj.la*
lrwxrwxrwx1rootroot16 5월 2 14:21/usr/local/lib/libproj.so -> libproj.so.0.8.0*
lrwxrwxrwx1rootroot16 5월 2 14:21/usr/local/lib/libproj.so.0 -> libproj.so.0.8.0*
- rwxr-xr-x1root root1775816 5월 2 14:21/usr/local/lib/libproj.so.0.8.0*
여기에 기호 링크가 있습니까?
인터넷에서 검색하면 해외 인터넷 사이트에 빈대 같다고 쓰여 있다.
데비안에서 우분투는 아직 안 된 것 같아요.
인터넷에 쓴 거 많이 만들어도 안돼.
하지만
$ export LD_LIBRARY_PATH=/lib:/usr/lib:/usr/local/lib
이렇게 하면 import가 됩니다.
상술한 사이트를 참고하다
국토 수치 정보 다운로드 서비스
http://nlftp.mlit.go.jp/ksj/
에서 데이터 다운로드
2. 정책 영역 > 행정 영역 > 다운로드할 데이터 선택
에서 설명한 대로 해당 매개변수의 값을 수정합니다.
다운로드 시간이 좋지 않은지 파일 크기가 크지 않지만 시간이 많이 걸린다.(수십 분) 마이크로소프트 이외에 이렇게 긴 시간이 걸린 것은 처음이다.
상기 사이트의 예제 프로그램을 실행하다
test.py
#!/usr/bin/python
# coding: UTF-8

import cartopy.crs as ccrs
import cartopy.io.shapereader as shpreader
import matplotlib.pyplot as plt

# 落としてきた行政区域のshpファイルを指定
fname = '/home/ty/python/map/japan/N03-170101_GML/N03-17_170101.shp'
shapes = list(shpreader.Reader(fname).geometries())

# 東京あたりを描画
ax = plt.axes(projection=ccrs.PlateCarree())
ax.add_geometries(shapes, ccrs.PlateCarree(), edgecolor='black', facecolor='gray', alpha=0.3)
ax.set_extent([139, 141, 35, 36], ccrs.PlateCarree())
plt.show()

처음에 눈치채지 못했는데 이 프로그램은 Basemap과 전혀 상관없어요.

하지만 국토교통성이 만든 지도답게 Basemap의'f'보다 정밀도가 좋다.
현경만 하거나 현경의 색깔을 변경하고 싶지만shp파일을 읽을 수 없기 때문에 현재 조사 중입니다.

좋은 웹페이지 즐겨찾기