rotpole projection to ESRI shape file

Bắt đầu bằng cách chạy đoạn mã Python được cung cấp và xem xét đầu ra proj4string của Basemap cùng với các hướng chuyển đổi pycrs và GDAL được đề cập trong báo cáo. Issue không nêu tệp mã nguồn hoặc bài kiểm thử nào của Basemap; để được xem là hoàn tất, cần xác nhận một cách xử lý phép chiếu rotpole này trong đầu ra ESRI shapefile, hoặc ghi lại một hạn chế nếu nó không được hỗ trợ.

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