[Date Prev][Date Next] [Thread Prev][Thread Next] [Date Index] [Thread Index]

Basemap issue for Debian 11 testing



I have been used Debian as GIS platform for years. Thanks to Debian.

Yesterday I began to test the codes of my book (Python and Open Source GIS) under Debian 11 testing.  https://github.com/bukun/book_python_gis
Now I found an issue.  And a similar link : https://stackoverflow.com/questions/65668542/projerror-x-y-z-and-time-must-be-same-size 

As follows codes:

-------------------------------------------------------------------------------------------------------
>>> from mpl_toolkits.basemap import Basemap
>>> para = {'projection': 'tmerc','lat_0': 0,'lon_0': 3,
...     'llcrnrlon': 1.81,
...     'llcrnrlat': 41.58,
...     'urcrnrlon': 1.84,
...     'urcrnrlat': 41.59}
>>>
>>> mymap = Basemap(**para)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "/usr/lib/python3/dist-packages/mpl_toolkits/basemap/__init__.py", line 1131, in __init__
    self._readboundarydata('gshhs',as_polygons=False)
  File "/usr/lib/python3/dist-packages/mpl_toolkits/basemap/__init__.py", line 1412, in _readboundarydata
    bx, by = self(blons, blats)
  File "/usr/lib/python3/dist-packages/mpl_toolkits/basemap/__init__.py", line 1192, in __call__
    xout,yout = self.projtran(x,y,inverse=inverse)
  File "/usr/lib/python3/dist-packages/mpl_toolkits/basemap/proj.py", line 293, in __call__
    outx,outy = self._proj4(x, y, inverse=inverse)
  File "/usr/lib/python3/dist-packages/pyproj/proj.py", line 182, in __call__
    return self.transform(
  File "/usr/lib/python3/dist-packages/pyproj/transformer.py", line 488, in transform
    self._transformer._transform(
  File "pyproj/_transformer.pyx", line 477, in pyproj._transformer._Transformer._transform
pyproj.exceptions.ProjError: x, y, z, and time must be same size
-------------------------------------------------------------------------------------------------------

Furthermore, i test some other paras. Seems very strange.

Following got same error:

-------------------------------------------------------------------------------------------------------
para = {'projection': 'tmerc','lat_0': 0,'lon_0': 3,
    'llcrnrlon': 1.8,
    'llcrnrlat': 41.5,
    'urcrnrlon': 1.9,
    'urcrnrlat': 41.6}
-------------------------------------------------------------------------------------------------------

Following got no error:

-------------------------------------------------------------------------------------------------------
para = {'projection': 'tmerc','lat_0': 0,'lon_0': 3,
    'llcrnrlon': 1,
    'llcrnrlat': 41,
    'urcrnrlon': 2,
    'urcrnrlat': 42}
para = {'projection': 'tmerc','lat_0': 0,'lon_0': 3,
    'llcrnrlon': 1.8,
    'llcrnrlat': 41,
    'urcrnrlon': 1.9,
    'urcrnrlat': 41}
-------------------------------------------------------------------------------------------------------

Reply to: