rotpole projection to ESRI shape file

まず、提供された Python スニペットを実行し、レポートで言及されている pycrs および GDAL の変換経路と併せて Basemap の proj4string 出力を確認します。issue では Basemap のソースファイルもテストも指定されていません。完了とするには、この rotpole 投影を ESRI shapefile 出力で処理するための方法を確認するか、サポートされていない場合は制限事項を文書化する必要があります。

Hi:

Does anyone know how to convert rotpole projection to ESRI format, to write data into a shape file?

I tried to do it using basemap.proj4string, and then pycrs or gdal, but there are some things missing from the proj4string...

I've added +ellps, but +o_proj=longlat +proj=ob_tran do not seem to be recognized...

Here is the proj4 I am trying to use:

+R=6370997.0 +lon_0=-97.0 +o_lat_p=42.5 +units=m +proj=ob_tran +o_lon_p=180.0 +o_proj=longlat +ellps=WGS84

pycrs complains that it does not know the specified name:

Traceback (most recent call last):
  File "/Users/huziy/skynet2/RESCUE/skynet3_rech1/huziy/Netbeans Projects/Python/RPN/src/crcm5/mh_domains/show_domain_with_drainage_area.py", line 213, in <module>
    main()
  File "/Users/huziy/skynet2/RESCUE/skynet3_rech1/huziy/Netbeans Projects/Python/RPN/src/application_properties.py", line 15, in wrapper
    func(*args, **kwargs)
  File "/Users/huziy/skynet2/RESCUE/skynet3_rech1/huziy/Netbeans Projects/Python/RPN/src/crcm5/mh_domains/show_domain_with_drainage_area.py", line 203, in main
    crs = pycrs.parser.from_proj4(bmp.proj4string + " +ellps={}".format("WGS84"))
  File "/Users/huziy/virtualenvs/py3.4/lib/python3.4/site-packages/pycrs/parser.py", line 565, in from_proj4
    raise Exception("The specified projection name could not be found")
Exception: The specified projection name could not be found

Below is my test code (the pj=basemap.proj4string + " +ellps=WGS84"):

    import pycrs
    pj = "+R=6370997.0 +lon_0=-97.0 +o_lat_p=42.5 +units=m +proj=ob_tran +o_lon_p=180.0 +o_proj=longlat  +ellps=WGS84"
    crs = pycrs.parser.from_proj4(pj)
    print(crs.to_esri_wkt())

Cheers