無法使用 python rasterio,gdal 打開 jp2（來自哨兵）

Question

我正在嘗試將 python 中的 jp2 光柵產品作為光柵打開，但當我們使用“raterio”和“gdal”包時沒有成功。
我得到這個錯誤：

RasterioIOError: 'b4.jp2' not recognized as a supported file format.
RuntimeError: `b4.jp2' not recognized as a supported file format.

我確實設法使用“glymur” package 將它們作為常規矩陣打開，但隨后我丟失了地理數據。

如果有人知道以光柵形式打開 jp2 的其他 pachage，我將很高興知道！

我試圖遵循這個答案： Python3.7 rasterio library cant open jp2 installed fiona by conda 或： https://gis.stackexchange.com/questions/250833/gdal-missing-jp2-driver-files-how-to-install -them/250837#250837並嘗試安裝： conda create -n mynewenv -c conda-forge gdal rasterio但我仍然遇到相同的錯誤。

b4 = rasterio.open('b4.jp2')
b4 = rasterio.open('b4.jp2', driver= 'JP2OpenJPEG')

Answer 1

這與您嘗試過的非常相似，但是此過程對我有用。 注意：我指定我想要一個 python 3.7 環境。 從 anaconda 提示：

conda create --name testing python=3.7
conda activate testing
conda install -c conda-forge gdal rasterio

驗證正確安裝：

(testing) C:\Users\dan>python
Python 3.7.3 | packaged by conda-forge | (default, Jul  1 2019, 22:01:29) [MSC v.1900 64 bit (AMD64)] :: Anaconda, Inc. on win32
Type "help", "copyright", "credits" or "license" for more information.
>>> import rasterio
>>> import gdal
>>> gdal.VersionInfo()
'3000100'
>>> rasterio.gdal_version()
'3.0.1'

並檢查以確保支持.jp2驅動程序。 從 anaconda 提示：

rio env --formats

它返回一個包含條目的長列表：

JML: OpenJUMP JML
JP2OpenJPEG: JPEG-2000 driver based on OpenJPEG library
JPEG: JPEG JFIF

這是來自conda env export的導出：

name: testing
channels:
  - conda-forge
  - defaults
dependencies:
  - affine=2.3.0=py_0
  - attrs=19.3.0=py_0
  - boost-cpp=1.70.0=h6a4c333_2
  - bzip2=1.0.8=hfa6e2cd_1
  - ca-certificates=2019.9.11=hecc5488_0
  - certifi=2019.9.11=py37_0
  - cfitsio=3.470=hfa6e2cd_2
  - click=7.0=py_0
  - click-plugins=1.1.1=py_0
  - cligj=0.5.0=py_0
  - curl=7.65.3=h4496350_0
  - expat=2.2.5=he025d50_1004
  - freetype=2.10.0=h563cfd7_1
  - freexl=1.0.5=hd288d7e_1002
  - gdal=3.0.1=py37h2fee047_10
  - geos=3.7.2=he025d50_2
  - geotiff=1.5.1=h6e819c3_4
  - gettext=0.19.8.1=hb01d8f6_1002
  - glib=2.58.3=hc0c2ac7_1002
  - hdf4=4.2.13=hf8e6fe8_1002
  - hdf5=1.10.5=nompi_ha405e13_1104
  - icu=64.2=he025d50_1
  - jpeg=9c=hfa6e2cd_1001
  - kealib=1.4.10=hf7dc31f_1005
  - krb5=1.16.3=hdd46e55_1001
  - libblas=3.8.0=14_mkl
  - libcblas=3.8.0=14_mkl
  - libcurl=7.65.3=h4496350_0
  - libffi=3.2.1=h6538335_1006
  - libgdal=3.0.1=h2491c49_10
  - libiconv=1.15=hfa6e2cd_1005
  - libkml=1.3.0=h4ece8bf_1010
  - liblapack=3.8.0=14_mkl
  - libnetcdf=4.7.1=nompi_h8d74e2a_101
  - libpng=1.6.37=h7602738_0
  - libpq=11.5=hb0bdaea_1
  - libspatialite=4.3.0a=hee5a54d_1031
  - libssh2=1.8.2=h642c060_2
  - libtiff=4.0.10=h6512ee2_1003
  - libxml2=2.9.10=h9ce36c8_0
  - lz4-c=1.8.3=he025d50_1001
  - numpy=1.17.3=py37hc71023c_0
  - openjpeg=2.3.1=hb24c2e3_1
  - openssl=1.1.1c=hfa6e2cd_0
  - pcre=8.43=h6538335_0
  - poppler=0.67.0=h92819f6_7
  - poppler-data=0.4.9=1
  - postgresql=11.5=h06f7779_1
  - proj=6.2.0=ha7a8c7b_1
  - pyparsing=2.4.2=py_0
  - python=3.7.3=h510b542_1
  - rasterio=1.1.0=py37h2617b1b_0
  - snuggs=1.4.7=py_0
  - tbb=2018.0.5=he980bc4_0
  - tiledb=1.6.2=h30b6f50_1
  - tk=8.6.9=hfa6e2cd_1003
  - xerces-c=3.2.2=h6538335_1004
  - xz=5.2.4=h2fa13f4_1001
  - zlib=1.2.11=h2fa13f4_1006
  - zstd=1.4.0=hd8a0e53_0
  - intel-openmp=2019.4=245
  - m2w64-expat=2.1.1=2
  - m2w64-gcc-libgfortran=5.3.0=6
  - m2w64-gcc-libs=5.3.0=7
  - m2w64-gcc-libs-core=5.3.0=7
  - m2w64-gettext=0.19.7=2
  - m2w64-gmp=6.1.0=2
  - m2w64-libiconv=1.14=6
  - m2w64-libwinpthread-git=5.0.0.4634.697f757=2
  - m2w64-xz=5.2.2=2
  - mkl=2019.4=245
  - msys2-conda-epoch=20160418=1
  - pip=19.3.1=py37_0
  - setuptools=41.6.0=py37_0
  - sqlite=3.30.1=he774522_0
  - vc=14.1=h0510ff6_4
  - vs2015_runtime=14.16.27012=hf0eaf9b_0
  - wheel=0.33.6=py37_0
  - wincertstore=0.2=py37_0
prefix: C:\Anaconda2\envs\testing

證明它有效：

(testing) C:\Users\RDCRLDDH>python
Python 3.7.5 (default, Oct 31 2019, 15:18:51) [MSC v.1916 64 bit (AMD64)] :: Anaconda, Inc. on win32
Type "help", "copyright", "credits" or "license" for more information.
>>> import rasterio
>>> ds = rasterio.open(r"C:\Users\RDCRLDDH\Desktop\T36QVM_20191023T082009_B01.jp2")
>>> ds.crs
CRS.from_epsg(32636)
>>> ds.read(1)
array([[2372, 2344, 2312, ..., 1423, 1406, 1481],
       [2326, 2292, 2279, ..., 1492, 1462, 1495],
       [2317, 2271, 2261, ..., 1529, 1517, 1511],
       ...,
       [1958, 1946, 1947, ..., 1697, 1734, 1730],
       [1960, 1960, 1967, ..., 1714, 1758, 1721],
       [1959, 1957, 1960, ..., 1769, 1752, 1701]], dtype=uint16)
>>> ds.transform
Affine(60.0, 0.0, 399960.0,
       0.0, -60.0, 2700000.0)
>>>

從 spyder，您可以檢查 spyder 從哪里運行

import sys
print(sys.executable)

您應該會看到如下內容：

C:\Anaconda2\envs\testing\pythonw.exe

安裝目錄C:\Anaconda2\將取決於你安裝 anaconda 的位置，但后面的envs\testing\pythonw.exe應該是一樣的。