Shulou(Shulou.com)06/03 Report--

This article mainly introduces how to solve the problem of pythongdal projection coordinate system conversion, the article is very detailed, has a certain reference value, interested friends must read it!

To convert the xian80 geographic coordinate system to a projected coordinate system:

Xian1980 = "GEOGCS [" GCS_Xian_1980 ", DATUM [" Xian_1980 ", SPHEROID [" Xian_1980 ", 6378140.0298.257]], PRIMEM [" Greenwich ", UNIT [" Degree ", 0.0174532925199433]]" prjreference = 'PROJCS ["Xian_1980_3_Degree_GK_Zone_37", GEOGCS ["GCS_Xian_1980", DATUM ["D_Xian_1980", SPHEROID ["Xian_1980", 6378140.0298.257], PRIMEM ["Greenwich", 0.0] UNIT ["Degree", 0.0174532925199433], PROJECTION ["Gauss_Kruger"], PARAMETER ["False_Easting", 37500000.0], PARAMETER ["False_Northing", 37500000.0], PARAMETER ["Central_Meridian", 111.0], PARAMETER ["Scale_Factor", 111.0], PARAMETER ["Latitude_Of_Origin", 37500000.0], UNIT ["Meter", 1.0]]'

An error will be reported after conversion:

Latlong = transform.TransformPoint (XQuery y) File "/ usr/lib/python3/dist-packages/osgeo/osr.py", line 711, in TransformPoint return _ osr.CoordinateTransformation_TransformPoint (self, * args) NotImplementedError: Wrong number or type of arguments for overloaded function 'CoordinateTransformation_TransformPoint'. Possible C + prototypes are: OSRCoordinateTransformationShadow::TransformPoint (double [3]) OSRCoordinateTransformationShadow::TransformPoint (double [3], double,double,double) * * think about it. I haven't found this problem before. Finally, it is found that the previous projection coordinate system: * * 'PROJCS ["Xian_1980_3_Degree_GK_Zone_37", GEOGCS ["GCS_Xian_1980", DATUM ["D_Xian_1980", SPHEROID ["Xian_1980", 6378140.0298.257], PRIMEM ["Greenwich", 0.0], UNIT ["Degree", 0.0174532925199433]], PROJECTION ["Transverse_Mercator"], PARAMETER ["False_Easting", 37500000.0], PARAMETER ["False_Northing", 37500000.0] PARAMETER ["Central_Meridian", 111.0], PARAMETER ["Scale_Factor", 1. 0], PARAMETER ["Latitude_Of_Origin", 0. 0], UNIT ["Meter", 1]]'* * is not a Gauss-Kruger projection but a horizontal Mercator. I think the reason is that gdal Curry can only support the Gauss-Kruger projection. **

The above is all the contents of the article "how to solve the problem of pythongdal projection coordinate system conversion". Thank you for reading! Hope to share the content to help you, more related knowledge, welcome to follow the industry information channel!

Welcome to subscribe "Shulou Technology Information " to get latest news, interesting things and hot topics in the IT industry, and controls the hottest and latest Internet news, technology news and IT industry trends.