problem with projection mil_os

[geomartino]

Hi

It's the first time I have to deal with a particular projection: Miller Oblated Stereographic (mil_os). According to the documentation, GDAL should be able to convert coordinates using that projection (via OSR).

So far with gdal,OSR in python

from osgeo import ogr, osr mil_os = osr.SpatialReference() mil_os.ImportFromProj4("+proj=mil_os +units=m +ellps=WGS84 +lon_0=20 +lat_0=18 +x_0=0 +y_0=0 +k_0=0.925 +no_defs")

is not working... But with pyproj, I can convert coordinate from mil_os to wgs84

from pyproj import Proj, transform mil_os = Proj(proj='mil_os')

Does anyone have a clue why it's not working with ogr/osr?

Since I want to reproject a shapefile from mil_os to wgs84, I could do it combining ogr/osr with pyproj but I would prefer convert it directly with ogr???

any thoughts?

thanks,

MartinO

[Even Rouault]

This has been improved in trunk (2.0dev) in r24063. Now this should work if you add +wktext at the ned of the your proj4 string

[geomartino]

Replying to rouault:

This has been improved in trunk (2.0dev) in r24063. Now this should work if you add +wktext at the ned of the your proj4 string

I've compile GDAL from trunk and the problem is gone. Thanks a lot for the fix!